Numerical simulations of cometary dust

Most observations of comets are done photometrically or spectrophotometrically. The interpretation of the aperture-averaged flux is relatively simple for an isotropic, radially expanding coma of infinite extent - the canonical model. However, the interpretation of the observations is not so clear when the motion of the dust is affected by radiation pressure, or when the emission is time-varying and anisotropic. For example, in a sample of CCD images of 10 comets, Jewitt and Meech (1987, Ap.J. 317, 992) found that the photometric profiles of only three comets were consistent, within the observational errors, with the profiles predicted form the canonical model. Photometric observations with large apertures, however, seem to suggest that the canonical model may be quite adequate (c.f. Osip, Schleicher, and Millis, 1992, Icarus 98, 115). The dust itself is characterized by a size distribution, with size dependencies on the expansion velocity, the scattered and thermal radiation, the response to radiation pressure, and probably the density. How good then are the approximations normally used in determining the production rates of the dust when these effects are present? As part of a program to better understand the dynamics of cometary dust and gas, a computer program has been developed which numerically simulates the emission of both dust and gas from a tilted rotating nucleus.