Modeling the dust size distribution in comets with dust fragmentation
A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.