Hostname: page-component-cd9895bd7-dk4vv Total loading time: 0 Render date: 2024-12-26T05:45:57.577Z Has data issue: false hasContentIssue false

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

Published online by Cambridge University Press:  20 January 2023

Shyam H. Menon
Affiliation:
Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia
Christoph Federrath
Affiliation:
Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia ARC Centre of Excellence for Astronomy in Three Dimensions (ASTRO-3D), Canberra, ACT 2611, Australia
Mark R. Krumholz
Affiliation:
Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia ARC Centre of Excellence for Astronomy in Three Dimensions (ASTRO-3D), Canberra, ACT 2611, Australia
Rolf Kuiper
Affiliation:
Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, Albert-Ueberle-Straße 2, 69120 Heidelberg, Germany
Benjamin D. Wibking
Affiliation:
Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia ARC Centre of Excellence for Astronomy in Three Dimensions (ASTRO-3D), Canberra, ACT 2611, Australia
Manuel Jung
Affiliation:
Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany email: shyam.menon@anu.edu.au
Rights & Permissions [Opens in a new window]

Abstract

Core share and HTML view are not available for this content. However, as you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used – but less accurate – closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities.

Type
Contributed Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of International Astronomical Union

References

Buntemeyer, L., Banerjee, R., Peters, T., Klassen, M., & Pudritz, R. E. 2016, Radiation hydrodynamics using characteristics on adaptive decomposed domains for massively parallel star formation simulations. New Astronomy, 43, 4969.CrossRefGoogle Scholar
Castor, J. I. 2004,. Radiation Hydrodynamics. Cambridge University Press.CrossRefGoogle Scholar
Dale, J. E. 2015,. The modelling of feedback in star formation simulations.CrossRefGoogle Scholar
Davis, S. W., Stone, J. M., & Jiang, Y.-F. 2012, A Radiation Transfer Solver for Athena Using Short Characteristics. ApJS, 199(1), 9.CrossRefGoogle Scholar
Gnedin, N. Y. & Ostriker, J. P. 1997, Reionization of the Universe and the Early Production of Metals. ApJ, 486(2), 581598.CrossRefGoogle Scholar
Jiang, Y. F., Stone, J. M., & Davis, S. W. 2012, A godunov method for multidimensional radiation magnetohydrodynamics based on a variable Eddington tensor. Astrophysical Journal, Supplement Series, 199(1).CrossRefGoogle Scholar
Jiang, Y.-F., Stone, J. M., & Davis, S. W. 2013, Saturation of the Magneto-rotational Instability in Strongly Radiation-dominated Accretion Disks. ApJ, 767(2), 148.CrossRefGoogle Scholar
Levermore, C. D. 1984, Relating Eddington factors to flux limiters. J. Quant. Spectrosc. Radiat. Transfer, 31(2), 149160.CrossRefGoogle Scholar
Levermore, C. D. & Pomraning, G. C. 1981, A flux-limited diffusion theory. ApJ, 248, 321334.CrossRefGoogle Scholar
Menon, S. H., Federrath, Ch ., Krumholz, M. R., Kuiper, R., Wibking, B. D., & Jung, M. 2022, VETTAM: a scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method. MNRAS, 512(1), 401.CrossRefGoogle Scholar
Mihalas, D. 1978,. Stellar atmospheres. W.H. Freeman, San Francisco.Google Scholar
Mihalas, D. & Klein, R. I. 1982, On the solution of the time-dependent inertial-frame equation of radiative transfer in moving media to O(v/c). Journal of Computational Physics, 46, 97137.CrossRefGoogle Scholar
Skinner, M. A. & Ostriker, E. C. 2013, A two-moment radiation hydrodynamics module in athena using a time-explicit Godunov method. Astrophysical Journal, Supplement Series, 206(2), 21.CrossRefGoogle Scholar
Stone, J. M., Mihalas, D., & Norman, M. L. 1992, ZEUS-2D: A Radiation Magnetohydrodynamics Code for Astrophysical Flows in Two Space Dimensions. III. The Radiation Hydrodynamic Algorithms and Tests. ApJS, 80, 819.CrossRefGoogle Scholar
Teyssier, R. & Commerçon, B. 2019, Numerical Methods for Simulating Star Formation. Frontiers in Astronomy and Space Sciences, 6, 51.CrossRefGoogle Scholar
Wibking, B. D. & Krumholz, M. R. 2021, Quokka: A code for two-moment AMR radiation hydrodynamics on GPUs. arXiv e-prints, arXiv:2110.01792.Google Scholar