Modeling the Cometary Environment using a Fluid Approach

Abstract

Comets are believed to have preserved the building material of the early solar system and to hold clues to the origin of life on Earth. Abundant remote observations of comets by telescopes and the in-situ measurements by a handful of space missions reveal that the cometary environments are complicated by various physical and chemical processes among the neutral gases and dust grains released from comets, cometary ions, and the solar wind in the interplanetary space. Therefore, physics-based numerical models are in demand to interpret the observational data and to deepen our understanding of the cometary environment.

In this thesis, three models using a fluid approach, which include important physical and chemical processes underlying the cometary environment, have been developed to study the plasma, neutral gas, and the dust grains, respectively. Although models based on the fluid approach have limitations in capturing all of the correct physics for certain applications, especially for very low gas density environment, they are computationally much more efficient than alternatives. In the simulations of comet 67P/Churyumov-–Gerasimenko at various heliocentric distances with a wide range of production rates, our multi-fluid cometary neutral gas model and multi-fluid cometary dust model have achieved comparable results to the Direct Simulation Monte Carlo (DSMC) model, which is based on a kinetic approach that is valid in all collisional regimes. Therefore, our model is a powerful alternative to the particle-based model, especially for some computationally intensive simulations.

Capable of accounting for the varying heating efficiency under various physical conditions in a self-consistent way, the multi-fluid cometary neutral gas model is a good tool to study the dynamics of the cometary coma with different production rates and heliocentric distances. The modeled H$_2$O expansion speeds reproduce the general trend and the speed's nonlinear dependencies of production rate and heliocentric distance, which are found in remote observations.

In the multi-fluid dust model, we use a newly developed numerical mesh to resolve the real shaped nucleus in the center and to facilitate prescription of the outer boundary conditions that accommodate the rotating frame. The model studies the effects of the rotating nucleus and the cometary activity in time-dependent simulations for the first time. The result also suggests that the rotation of the nucleus explains why there is no clear dust speed dependence on size in some of the dust observations.

We developed a new multi-species comet MHD model to simulate the plasma environment of comet C/2006 P1 (McNaught) over a wide range of heliocentric distances from 0.17 AU to 1.75 AU, with the constraints provided by remote and in situ observations. Typical subsolar standoff distances of bow shock and contact surface are modeled and presented to characterize the solar wind interaction of the comet at various heliocentric distances. In addition, the model is also the first one to be used to study the composition and dynamics in the distant cometary tail. The results agree well with the measured water group ion abundances from the Ulysses/SWICS 1.7 AU down-tail from the comet and the velocity and temperature measured by Ulysses/SWOOPS.