Multi-fluid model of comet 1P/Halley

A new three-dimensional magnetohydrodynamic model of the coma of a comet has been developed and applied to simulations of a Halley-class coma using the solar-wind conditions of the Giotto flyby of Halley in 1986. The code developed for high-performance parallel processing computers, combines the high spatial resolution of smaller than 1 km grid spacing near the nucleus, with a large computational domain that enables structures nearly 10 million km down the comet tail to be modeled. Ions, neutrals, and electrons are considered as separate interacting fluids. Significant physical processes treated by the model include both photo and electron impact ionization of neutrals, recombination of ions, charge exchange between solar-wind ions and cometary neutrals, and frictional interactions between the three fluids considered in the model. A variety of plasma structures and physical parameters that are the output of this model are compared with relevant Giotto data from the 1986 Halley flyby.