Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer

Details Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer

2001-03-15

arxiv.org/abs/astro-ph/0103249v1

Astrophys.J.554:615-623,2001

Astronomy and Astrophysics

Astrophysics

24 pages, 8 figures

Scientific paper

10.1086/321336

We describe a general radiative equilibrium and temperature correction procedure for use in Monte Carlo radiation transfer codes with sources of temperature-independent opacity, such as astrophysical dust. The technique utilizes the fact that Monte Carlo simulations track individual photon packets, so we may easily determine where their energy is absorbed. When a packet is absorbed, it heats a particular cell within the envelope, raising its temperature. To enforce radiative equilibrium, the absorbed packet is immediately re-emitted. To correct the cell temperature, the frequency of the re-emitted packet is chosen so that it corrects the temperature of the spectrum previously emitted by the cell. The re-emitted packet then continues being scattered, absorbed, and re-emitted until it finally escapes from the envelope. As the simulation runs, the envelope heats up, and the emergent spectral energy distribution (SED) relaxes to its equilibrium value, without iteration. This implies that the equilibrium temperature calculation requires no more computation time than the SED calculation of an equivalent pure scattering model with fixed temperature. In addition to avoiding iteration, our method conserves energy exactly, because all injected photon packets eventually escape. Furthermore, individual packets transport energy across the entire system because they are never destroyed. This long-range communication, coupled with the lack of iteration, implies that our method does not suffer the convergence problems commonly associated with lambda-iteration. To verify our temperature correction procedure, we compare our results to standard benchmark tests, and finally we present the results of simulations for two-dimensional axisymmetric density structures.

Affiliated with

Also associated with

No associations

Rating

Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Radiative Equilibrium and Temperature Correction in Monte Carlo Radiation Transfer will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-565846