This paper considers the gas phase chemistry in a protoplanetary accretion disk, especially the chemistry initiated in the gas phase by destruction of dust close to the central star. Slow radial particle transport moves gas and dust from the cold outer parts of a protoplanetary accretion disk into its warm central part where chemical reactions in the gas phase are acti- vated. At the same time gases frozen on the surface of dust grains are vaporized and later the dust grains themselves are vaporized or destroyed by chemical surface reactions. These processes initiate a rich chemistry in the protoplanetary accretion disk. The simulation of chemical reactions, as in the case of an accretion disk, mostly leads to a large and stiff system of differ- ential or differential-algebraic equations. For the integration of such systems implicit methods are required. We present an ef- cient BDF-method and give a detailed description of the error and stepsize control and the strategies to minimize the numer- ical effort of the linear algebra problems. Typical applications for chemical processes (chemistry and dust destruction) in an accretion disk are treated with this method. The corresponding code DAESOL turned out to be more robust and much faster than the more conventional code used rst. Some results for the chemistry in a protoplanetary accretion disk are briefly discussed.
[1]
J. Klett,et al.
Microphysics of Clouds and Precipitation
,
1978,
Nature.
[2]
D. Raine,et al.
Accretion power in astrophysics
,
1985
.
[3]
S. Nash,et al.
Numerical methods and software
,
1990
.
[4]
W. M. Haynes.
CRC Handbook of Chemistry and Physics
,
1990
.
[5]
Geochimica et Cosmochimica Acta
,
1992
.
[6]
M. Skrutskie,et al.
in Protostars and Planets III
,
1993
.
[7]
Linda R. Petzold,et al.
Numerical solution of initial-value problems in differential-algebraic equations
,
1996,
Classics in applied mathematics.