Radiation implicit step.
More...
Main routines used for the implicit integration of the radiation-matter interaction terms. Four different implicit methods are here implemented:
- RADIATION_FIXEDPOINT_RAD: Solves the implicit step following Takahashi & Oshuga (2013), by iterating the radiation fields.
- RADIATION_FIXEDPOINT_GAS: fixed-point algorithm based on iterations of the MHD fields.
- RADIATION_NEWTON_GAS: Solves the implicit step by means of the Newton's method, performing iterations of the matter fields (p_g,u^i), i.e., the matter pressure and spatial components of its 4-velocity. The components of the Jacobian are computed numerically by carrying small variations of these fields.
RADIATION_NEWTON_RAD: Same as RADIATION_NEWTON_RAD, but performing iterations of the radiation fields.
In all methods the convergence criterion considers both radiation and MHD fields unless RADIATION_FULL_CONVERGENCE == NO, in which case only the variation of the iterated fields is controlled.
- Authors
- J. D. Melon Fuksman (fuksm.nosp@m.an@m.nosp@m.pia.d.nosp@m.e)
- Date
- Dec 03, 2020
References
- Melon Fuksman, J. D., and Mignone, A. 2019, "A Radiative Transfer Module for Relativistic Magnetohydrodynamics in the PLUTO Code" ApJS, 242, 20.
- Takahashi, H. R., & Ohsuga, K. 2013, "A Numerical Treatment of Anisotropic Radiation Fields Coupled With Relativistic Resistive Magnetofluids", ApJ, 772, 127.