EDP Sciences
Free access
Volume 417, Number 1, April I 2004
Page(s) 317 - 324
Section Stellar atmospheres
DOI http://dx.doi.org/10.1051/0004-6361:20034473

A&A 417, 317-324 (2004)
DOI: 10.1051/0004-6361:20034473

Improved discretization of the wavelength derivative term in CMF operator splitting numerical radiative transfer

P. H. Hauschildt1, 2 and E. Baron3, 4, 5

1  Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany
    e-mail: yeti@hs.uni-hamburg.de
2  Dept. of Physics and Astronomy & Center for Simulational Physics, University of Georgia, Athens, GA 30602-2451, USA
3  Dept. of Physics and Astronomy, University of Oklahoma, 440 W. Brooks, Rm 131, Norman, OK 73019-0225, USA
4  NERSC, Lawrence Berkeley National Laboratory, MS 50F-1650, 1 Cyclotron Rd, Berkeley, CA 94720-8139, USA
5  Laboratoire de Physique Nucléaire et de Haute Énergies, CNRS-IN2P3, University of Paris VII, 75005 Paris, France

(Received 8 October 2003 / Accepted 17 December 2003 )

We describe two separate wavelength discretization schemes that can be used in the numerical solution of the comoving frame radiative transfer equation. We present an improved second order discretization scheme and show that it leads to significantly less numerical diffusion than the previous scheme. We also show that due to the nature of the second order term in some extreme cases it can become numerically unstable. We stabilize the scheme by introducing a mixed discretization scheme and present the results from several test calculations.

Key words: radiative transfer -- methods: numerical

Offprint request: E. Baron, baron@nhn.ou.edu

© ESO 2004