home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #26 / NN_1992_26.iso / spool / sci / astro / 11895 < prev    next >
Encoding:
Internet Message Format  |  1992-11-11  |  4.7 KB
  1. Path: sparky!uunet!charon.amdahl.com!pacbell.com!iggy.GW.Vitalink.COM!cs.widener.edu!eff!sol.ctr.columbia.edu!zaphod.mps.ohio-state.edu!magnus.acs.ohio-state.edu!usenet.ins.cwru.edu!agate!doc.ic.ac.uk!uknet!comlab.ox.ac.uk!oxuniv!wgw
  2. From: wgw@vax.oxford.ac.uk
  3. Newsgroups: sci.astro
  4. Subject: Swift_tuttle's orbit
  5. Message-ID: <1992Nov11.165126.10076@vax.oxford.ac.uk>
  6. Date: 11 Nov 92 16:51:26 GMT
  7. Organization: Oxford University VAX 6620
  8. Lines: 101
  9.  
  10.  
  11.           The following should put all of the recent hoo-hah over 
  12.           comet P/Swift-Tuttle's next apparition into perspective.
  13.  
  14.  
  15. The 2126 apparition of P/Swift-Tuttle
  16. =====================================
  17.  
  18.    As a result of the uncertainty in deriving the nongravitational 
  19. parameters (A1 and A2) and their possible time-dependence, much has been
  20. made of the possibility that P/Swift-Tuttle represents a threat to the
  21. Earth at its next apparition in 2126.   The intersection of the comet's
  22. orbit and that of the Earth at the former's descending node gives rise
  23. to the annual Perseid meteor shower.  In 2126 the Earth will intersect
  24. the shower maximum on August 14.3, when it will have  r = 1.01338 au.
  25. For the comet to be at its descending node ( r = 1.01001 au ) at the
  26. same time requires that it have a perihelion passage time of July 26.43.
  27. Thus, to represent a collison threat to the Earth, the comet's time of
  28. perihelion passage has to be delayed by 20 days with respect to the 
  29. purely gravitational integration (details below).  An integration including
  30. the effects of the nongravitational (outgassing) terms, which is a better 
  31. fit to the previous two apparitions, brings the time of perihelion passage 
  32. forward.  From the A2 = -0.035 integration below we see that the comet will
  33. be at its descending node on 2126 July 18.5, at which time it will be 0.45au
  34. from the Earth.
  35.  
  36.  
  37.  
  38. Details of RADAU integrations of P/Swift-Tuttle's orbit
  39. =======================================================
  40.  
  41. Starting elements taken from IAUC 5636  (epoch 1992 Dec. 4.0).
  42.  
  43. J2000 throughout.   Planetary starting positions and velocities from DE200.
  44. Perturbations by Venus -> Neptune included ; mass of Mercury added to that
  45. of the Sun  (other runs including all nine planets explicitly show that the
  46. error in time of perihelion passage in 2126 does not exceed  0.1 day  when
  47. using this approximation).
  48.  
  49. All runs were done on the VAX 6000/620 of Oxford University Computer Service.
  50.  
  51. For the nongrav run no attempt was made to optimise the values of A1 and A2, 
  52. the values used being chosen to approximate the probable times of perihelion 
  53. passage in 1737 and 1862.
  54.  
  55.  
  56. N-body results  ( no outgassing effects considered )
  57. ==============
  58.  
  59.      epoch          q          e         w         asc         i         nu
  60.  
  61.   2308960.500   0.976371   0.962033   152.8388   139.5157   113.5590   254.7898
  62.   2355360.500   0.979797   0.961454   152.6993   139.4771   113.6880   235.7157
  63.   2401760.500   0.962661   0.962659   152.7805   139.3815   113.5693   129.3062
  64. * 2448960.500   0.958120   0.963590   153.0130   139.4560   113.4300   347.8336
  65.   2497760.500   0.956312   0.963880   153.1270   139.6201   113.4081    12.4523
  66.  
  67.  
  68.  
  69. Nongrav results   ( outgassing effects included : A1 = 0 , A2 = -0.035 )
  70. ===============
  71.  
  72.      epoch          q          e         w         asc         i         nu
  73.  
  74.   2308960.500   0.976695   0.962003   152.8302   139.5090   113.5599   253.4023
  75.   2355360.500   0.980095   0.961423   152.6924   139.4716   113.6843   238.4297
  76.   2401760.500   0.962645   0.962661   152.7804   139.3821   113.5691   129.5740
  77. * 2448960.500   0.958120   0.963590   153.0130   139.4560   113.4300   347.8336
  78.   2497760.500   0.956460   0.963875   153.1236   139.6194   113.4039    22.1289
  79.  
  80.     ( * indicates starting values )
  81.     (    nu is the true anomaly   )
  82.  
  83.  
  84. Thus, for these osculating epochs, we find the following times of successive 
  85. perihelion passages (Gregorian calendar) for the comet,
  86.  
  87.                 epoch             T (n-body)         T (non-grav)
  88.  
  89.            1609 Aug. 14.0      1610 Jan. 20.98     1610 Jan. 27.94    
  90.            1736 Aug. 28.0      1737 July  7.12     1737 June  5.00
  91.            1863 Sept 12.0      1862 Aug. 28.33     1862 Aug. 23.66
  92.         *  1992 Dec.  4.0                1992 Dec. 12.32
  93.            2126 July 15.0      2126 July 6.50      2126 June 29.64
  94.  
  95.  
  96. By way of comparison, Brian Marsden (A.J. vol 78, p.654-662; 1973) gives 
  97. the time of perihelion passage for 1862 as August 23.41.    For 1737 he 
  98. gives a range of solutions with perihelion passage times varying in the 
  99. range June 3.5 to June 16.4.
  100.  
  101.  
  102. NOTE
  103. ====
  104.    It should be stressed that all of the foregoing is based entirely on
  105. the osculating elements given on IAUC 5636 (dated 1992 October 15) which
  106. have presumably already been improved upon.
  107.  
  108.  
  109.  
  110.         Graeme Waddington                        wgw @ vax.ox.ac.uk
  111.