ABSTRACT. The hydrodynamic approach to the modeling of the coUisional regions of the cometary atmospheres is reviewed. A concise account of the existing Monte Carlo simulations is also incorporated. Emphasis is placed on works stimulated by the Halley flyby missions. The basic physical processes are described, based on the analogy with the physics of underexpanded jets. The general structure of the governing equations is presented and their legacy discussed. A detailed intercomparison is made of the application of these methods to the modeling of the comas. The reliability of the conclusions brought forward by existing models is discussed.
. IntroductionIn view of the extreme rarity of past and foreseeable in situ cometary missions, remote sensing of comets has been and will continue to be the essential tool of cometary investigation. Therefore, most properties (and, in particular, the statistical properties) of the solid nucleus of all comets have been and will continue to be inferred from observations of the cometary atmosphere-the coma. With the refinements of the observational techniques, a detailed and reliable understanding of the formation of the coma has emerged as a key objective of cometary physics.The application of gas-dynamic methods to the formation of the coma seems to date only back to the year 1965. Fluid equations were first used in order to account for the acceleration of cometary dust and for the formation of many minor molecules. It is not a simple problem to delineate the precise extent of the collision-dominated part of a coma, i.e., that region where hydrostatic or hydrodynamic methods provide satisfactory results (see Sections 3.3 to 3.5). Primitive considerations, based on evaluations of the ratio of the local collision mean free path (m.f.p.) to the cometocentric distance, suggest that a waterdominated coma has a very small collision zone beyond an heliocentric distance = 3 AU, while that zone covers the whole water survival zone within = 0.5 AU (see, e.g., [13,17]). The essential point, however, is that the distribution of matter in the coma, the dynamical properties of all constituents, and the chemical composition of the gas phase are 937