Astronomy and Astrophysics – Astrophysics
Scientific paper
2008-12-09
2009MNRAS.396.1383P
Astronomy and Astrophysics
Astrophysics
21 pages, 17 figures, published in MNRAS
Scientific paper
10.1111/j.1365-2966.2009.14843.x
We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code {\small GADGET}. It is based on a fast, robust and photon-conserving integration scheme where the radiation transport problem is approximated in terms of moments of the transfer equation and by using a variable Eddington tensor as a closure relation, following the `OTVET'-suggestion of Gnedin & Abel. We derive a suitable anisotropic diffusion operator for use in the SPH discretization of the local photon transport, and we combine this with an implicit solver that guarantees robustness and photon conservation. This entails a matrix inversion problem of a huge, sparsely populated matrix that is distributed in memory in our parallel code. We solve this task iteratively with a conjugate gradient scheme. Finally, to model photon sink processes we consider ionisation and recombination processes of hydrogen, which is represented with a chemical network that is evolved with an implicit time integration scheme. We present several tests of our implementation, including single and multiple sources in static uniform density fields with and without temperature evolution, shadowing by a dense clump, and multiple sources in a static cosmological density field. All tests agree quite well with analytical computations or with predictions from other radiative transfer codes, except for shadowing. However, unlike most other radiative transfer codes presently in use for studying reionisation, our new method can be used on-the-fly during dynamical cosmological simulation, allowing simultaneous treatments of galaxy formation and the reionisation process of the Universe.
Petkova Margarita
Springel Volker
No associations
LandOfFree
An implementation of radiative transfer in the cosmological simulation code GADGET 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 An implementation of radiative transfer in the cosmological simulation code GADGET, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and An implementation of radiative transfer in the cosmological simulation code GADGET will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-526748