Global moedeling of comets: nucleus, neutral and ionized coma of comets 67P/Churyumov-Gerasimenko and 46P/Wirtanen. Preparations for the ROSETTA radio science investigations

Models of the thermal behaviour of a cometary nucleus, the evolution of the neutral gas coma, the ionized cometary coma and of the interaction of the cometary plasma with the solar wind are studied in this work. The general aim is to develop a global model of the comet and its environment in order to characterize the physical conditions around comets 67P/Churyumov-Gerasimenko and 46P/Wirtanen with respect to the heliocentric distance. The model of the heat diffusion within the cometary nucleus is one-dimensional. A grid of one-dimensional models is distributed over the nucleus in order to determine the temperature distribution and the sublimation characteristics of the comet on the whole surface of the comet. A heat balance equation is applied as boundary condition on the surface. Many parameters that have to be accounted for in a heat diffusion model are not precisely known to date. The variation of these parameters within reasonable limits yields a wide range of possible results. The heat diffusion within the cometary nucleus is derived from an energy conservation equation that includes heat conduction through the porous cometary material and heat convection due to the transport of latent heat by the gas phase within the nucleus. Model results are evaluated by a comparison of modeled and observed global gas production rates. Exemplary maps of the local temperature distribution and local sublimation rates at particular heliocentric distances are also provided. The neutral gas coma of the comet is modeled with a hydrodynamic approximation. The acceleration of the spacecraft due to the gas mass flux is evaluated with the model results. The ionized coma of a comet can also have an effect on the carrier signal. A one-dimensional model of the plasma density at the comet-sun axis is developed. The assumption of photochemical equilibrium is not necessarily justified within the coma of weak outgassing comets. The continuity equation of the plasma density has to be solved without this simplifying assumption. A model of the electron temperature profile is also generated. The transition from a regime where electrons are effectively cooled to a region with temperatures of the electron fluid similar to solar wind levels is assumed to set in at the position of the thermal electron collisionopause. The interaction of the cometary plasma with the solar wind is also studied. (orig.)