J/A+A/559/A80 Radiation hydrodynamics integrated in the PLUTO code (Kolb+, 2013)
Radiation hydrodynamics integrated in the PLUTO code.
Kolb S.M., Stute M., Kley W., Mignone A.
<Astron. Astrophys., 559, A80 (2013)>
=2013A&A...559A..80K 2013A&A...559A..80K (SIMBAD/NED BibCode)
ADC_Keywords: Models
Keywords: radiative transfer - hydrodynamics - accretion, accretion disks
Abstract:
The transport of energy through radiation is very important in many
astrophysical phenomena. In dynamical problems the time-dependent
equations of radiation hydrodynamics have to be solved. We present a
newly developed radiation-hydrodynamics module specifically designed
for the versatile magnetohydrodynamic (MHD) code PLUTO. The solver is
based on the flux-limited diffusion approximation in the
two-temperature approach. All equations are solved in the co-moving
frame in the frequency-independent (gray) approximation. The
hydrodynamics is solved by the different Godunov schemes implemented
in PLUTO, and for the radiation transport we use a fully implicit
scheme. The resulting system of linear equations is solved either
using the successive over-relaxation (SOR) method (for testing
purposes) or using matrix solvers that are available in the PETSc
library. We state in detail the methodology and describe several test
cases to verify the correctness of our implementation. The solver
works in standard coordinate systems, such as Cartesian, cylindrical,
and spherical, and also for non-equidistant grids. We present a new
radiation-hydrodynamics solver coupled to the MHD-code PLUTO that is a
modern, versatile, and efficient new module for treating complex
radiation hydrodynamical problems in astrophysics. As test cases,
either purely radiative situations, or full radiation-hydrodynamical
setups (including radiative shocks and convection in accretion disks)
were successfully studied. The new module scales very well on parallel
computers using MPI. For problems in star or planet formation, we
added the possibility of irradiation by a central source.
Description:
This is the modified version of the radiation.c module implemented in
Kolb et al. (2013A&A...559A..80K 2013A&A...559A..80K).
The module requires a set of opacity tables (Plank and Rosseland
opacities, and radiative losses). These tables has to be read by a
file named opacity.c. It is also mandatory to include a table for the
partial derivative of the radiative losses with respect to the
temperature.
To download the full code please go to
'http://cerere.astropa.unipa.it/progetti_ricerca/HPC/resources.htm'.
For further informations please contact the authors.
File Summary:
--------------------------------------------------------------------------------
FileName Lrecl Records Explanations
--------------------------------------------------------------------------------
ReadMe 80 . This file
radiation.zip 512 38 Modified version of the radiation.c module
--------------------------------------------------------------------------------
See also:
http://cerere.astropa.unipa.it/progetti_ricerca/HPC/resources.htm : Code
Acknowledgements:
Salvatore Colombo, salvatore.colombo(at)inaf.it
(End) Patricia Vannier [CDS] 22-Jul-2019