Volume 548, December 2012
|Number of page(s)||9|
|Section||Numerical methods and codes|
|Published online||22 November 2012|
A 3D radiative transfer framework
X. Arbitrary velocity fields in the comoving frame
1 Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany
e-mail: email@example.com; firstname.lastname@example.org
2 Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma, 440 W Brooks, Rm 100, Norman, OK 73019, USA
e-mail: email@example.com; firstname.lastname@example.org
3 Computational Research Division, Lawrence Berkeley National Laboratory, MS 50F-1650, 1 Cyclotron Rd, Berkeley, CA 94720, USA
4 Physics Department, Univ. of California, Berkeley, CA 94720, USA
Received: 4 April 2012
Accepted: 20 October 2012
Aims. General 3D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1D problem with arbitrary flows to those that we have developed for solving the fully 3D relativistic radiative transfer problem for monotonic flows.
Methods. The methods developed for 3D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3D. Simultaneously, the techniques we developed for treating the 1D problem with arbitrary velocity fields are easily adapted to the 3D problem.
Results. The algorithm we present can be used to solve 3D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.
Key words: radiative transfer
© ESO, 2012
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.