A kinetic model of gas flow in a porous cometary mantle

A numerical study of gas flow through a porous cometary mantle is presented. A kinetic model based on the well-known Test Particle Monte Carlo Method for the solution of rarefied gas dynamics problems is proposed. The physical model consists of two spatial plane regions: the condensed ice phase and a porous dust mantle. The structure of the porous dust layer is described as a bundle of cylindrical inclined channels not crossing each other. A vertical temperature gradient may exist across the dust mantle. The aim is to investigate how the characteristics of molecular flow depend on the capillary length, inclination angle, and temperature gradient. Examples illustrating a significant deviation of some results from equilibrium values are shown. In particular, the gas velocity distribution at both ends of the pore is strongly non-Maxwellian if there is an important temperature contrast across the pore. The emergent gas flow rate is found to vary with the pore length/radius ratio in excellent agreement with Clausing's empirical formula. The degree of collimation of the flow is quantitatively studied as a function of the length/radius ratio, and consequences for the jet force of outgassing through a dust mantle or, indeed, a rough surface are estimated.     

A Monte Carlo investigation of Jovian perturbations on short-period comet orbits

We present statistical distributions of Jovian perturbations on short-period comet orbits resulting from accurate numerical integrations. Our sample of 60, 000 cometary orbits with low inclinations and random orientations is characterized by perihelia between 0 and 7 AU and aphelia between 4 and 13 AU. The perturbations considered are those experienced because of Jupiter's gravitation per orbital revolution by the comets. Regularization and accurate step-length control in the numerical integration gives statistical results appreciably different from those computed by Rickman and Vaghi (1978). Their use of a crude method of integration led to erroneous results for close encounters. Strong asymmetries of the $\text{δ(}\text{1}\text{a}\text{)}$ distributions, in particular for the extreme tails, are observed for perihelion- or aphelion-tangent orbits. These orbits are also shown to experience the strongest energy perturbations on the average. Some results concerning the perturbations of Tisserand parameters are indicated. The perturbation distributions for the angular elements are described and discussed. The role of the minimum distance from Jupiter as an indicator of perturbations is investigated.     

J. Mayo Greenberg (1922–2001) In Memoriam And Dedication

Earth, Moon, and Planets -     

Orbital linkages of Comet 6P/d'Arrest based on its asymmetric light curve

We present an investigation of different models of the nongravitational acceleration on Comet 6P/d'Arrest, as used in orbital linkages spanning 150 years from the discovery of the comet in 1851 until the recent observations made in 2001. Some of our models use the time-shifted g-like function to represent the variation of outgassing rate, but the main thrust is on models using instead a production curve that is fitted to recent light curve observations—mainly those in 1976. We pay special attention to the proper scaling of such a production curve, when applied to other apparitions with a different perihelion distance q, and we find a best fit with a q^−1.6 power-law. Generally, the best fit is found with models, in which the acceleration components are expressed in terms of the angular parameters of the rotating nucleus. We thus find the orientation of the spin axis, and using the orbital evolution we are able to predict a variable time shift of the outgassing curve. The very best results are found when applying this time shift to the light-curve based, angular models. The totality of the 1851-2001 observations can then be linked with a mean residual of less than 4″. This may be brought down to ∼^2″ by solving for individual ‘activity parameters’ of all apparitions, which are multiplicative factors applied to the acceleration amplitudes. These turn out to be within 10% of unity for the best fit. We have also performed a linkage to the observations of Comet 1678 (La Hire) using our models. We find an indication of a secular increase of the amount of asymmetry of the outgassing with respect to perihelion, part of which is due to the variable time shift caused by the orbital evolution.     

Orbits of short period comets captured from the Oort cloud

Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.     

Algorithms for Stellar Perturbation Computations on Oort Cloud Comets

We investigate different approximate methods of computing the perturbations on the orbits of Oort cloud comets caused by passing stars, by checking them against an accurate numerical integration using Everhart’s RA15 code. The scenario under study is the one relevant for long-term simulations of the cloud’s response to a predefined set of stellar passages. Our sample of stellar encounters simulates those experienced by the Solar System currently, but extrapolated over a time of 10¹⁰ years. We measure the errors of perihelion distance perturbations for high-eccentricity orbits introduced by several estimators – including the classical impulse approximation and Dybczyński’s (1994, Celest. Mech. Dynam. Astron. 58, 1330–1338) method – and we study how they depend on the encounter parameters (approach distance and relative velocity). We introduce a sequential variant of Dybczyński’s approach, cutting the encounter into several steps whereby the heliocentric motion of the comet is taken into account. For the scenario at hand this is found to offer an efficient means to obtain accurate results for practically any domain of the parameter space.     

Observations of Comet 9P/Tempel 1 around the Deep Impact event by the OSIRIS cameras onboard Rosetta

The OSIRIS cameras on the Rosetta spacecraft observed Comet 9P/Tempel 1 from 5 days before to 10 days after it was hit by the Deep Impact projectile. The Narrow Angle Camera (NAC) monitored the cometary dust in 5 different filters. The Wide Angle Camera (WAC) observed through filters sensitive to emissions from OH, CN, Na, and OI together with the associated continuum. Before and after the impact the comet showed regular variations in intensity. The period of the brightness changes is consistent with the rotation period of Tempel 1. The overall brightness of Tempel 1 decreased by about 10% during the OSIRIS observations. The analysis of the impact ejecta shows that no new permanent coma structures were created by the impact. Most of the material moved with . Much of it left the comet in the form of icy grains which sublimated and fragmented within the first hour after the impact. The light curve of the comet after the impact and the amount of material leaving the comet ( of water ice and a presumably larger amount of dust) suggest that the impact ejecta were quickly accelerated by collisions with gas molecules. Therefore, the motion of the bulk of the ejecta cannot be described by ballistic trajectories, and the validity of determinations of the density and tensile strength of the nucleus of Tempel 1 with models using ballistic ejection of particles is uncertain.     

From the Oort cloud to observable short-period comets – I. The initial stage of cometary capture

We investigate the first stage of the dynamical evolution of Oort cloud comets entering the planetary region for the first time. To this purpose, we integrate numerically the motions of a large number of fictitious comets pertaining to two samples, both with perihelion distances up to 5.7 au and random inclinations; the first sample is composed of comets whose orbits have at least one node close to 5.2 au, while the second is not subject to this constraint. We examine the orbits when the comets come to aphelion after their first perihelion passage within the planetary region, and find that there is a clear statistical dependence of the energy perturbations on the Tisserand parameter. There appear to be two main processes, of comparable importance, governing the shortening of semimajor axes to values of less than 1000 au, i.e. planetary close encounters, especially with Jupiter, and indirect perturbations due to the shifting of the motion from barycentric to heliocentric and back; the former process mostly affects comets crossing the ecliptic at about 5.2 au, or on low-inclination orbits, while the latter mostly affects comets of small perihelion distance. This last result may help to understand the relative paucity of Halley-type comets with perihelion distances larger than about 1.5 au.     

A Monte Carlo estimate of the fraction of comets developing into sizeable asteroidal bodies

We present results of Monte Carlo simulations of orbital evolution showing that assuming a steady state there are roughly 50 extinct comets per active one in the Mars-crossing Jupiter family. The large number of extinct comets thus expected compared with the absence of observed apollo or Amor asteroids with aphelion distances greater than 4.2 AU indicates that less than five percent of the extinct comets survive as sizeable asteroidal bodies.Paper presented at the European Workshop on Planetary Sciences, organised by the Laboratorio di Astrofisica Spaziale di Frascati, and held between April 23–27, 1979, at the Accademia Nazionale del Lincei in Rome, Italy.     

The Oort Cloud and long‐period comets

This review starts with a brief historical overview of the subject, after which some recent papers attempting to improve the understanding of comet injection from the Oort Cloud and the origin of new comets are discussed. Special attention is paid to the importance of nongravitational effects in comet orbit determination, the synergy between stellar encounters and the galactic tides for the injection dynamics, and the role of planetary perturbations. The field is thus shown to be advancing rapidly, and brief comments on possible implications for studying the origin of the cloud are made.