Go to the source code of this file.

Modules
module	mod_rhd_phys
	Radiation-Hydrodynamics physics module Module aims at solving the Hydrodynamic equations toghether with the zeroth moment of the radiative transfer equation. A closure is provided by the flux limited diffusion (FLD)-approximation in the mod_fld.t module. See [1]Moens, N., Sundqvist, J. O., El Mellah, I., Poniatowski, L., Teunissen, J., and Keppens, R., “Radiation-hydrodynamics with MPI-AMRVAC . Flux-limited diffusion”, Astronomy and Astrophysics, vol. 657, 2022. doi:10.1051/0004-6361/202141023. For more information. Another possible closure in the works is the anisotropic flux limited diffusion approximation (AFLD) described in mod_afld.t.

Functions/Subroutines
subroutine, public	mod_rhd_phys::rhd_phys_init ()
	Initialize the module.

subroutine, public	mod_rhd_phys::rhd_check_params

subroutine, public	mod_rhd_phys::rhd_set_mg_bounds
	Set the boundaries for the diffusion of E.

subroutine, public	mod_rhd_phys::rhd_check_w (primitive, ixil, ixol, w, flag)
	Returns logical argument flag where values are ok.

subroutine, public	mod_rhd_phys::rhd_to_conserved (ixil, ixol, w, x)
	Transform primitive variables into conservative ones.

subroutine, public	mod_rhd_phys::rhd_to_primitive (ixil, ixol, w, x)
	Transform conservative variables into primitive ones.

subroutine, public	mod_rhd_phys::rhd_get_csound2 (w, x, ixil, ixol, csound2)
	Calculate the square of the thermal sound speed csound2 within ixO^L. csound2=gamma*p/rho.

subroutine, public	mod_rhd_phys::rhd_get_pthermal (w, x, ixil, ixol, pth)
	Calculate thermal pressure=(gamma-1)(e-0.5m**2/rho) within ixO^L.

subroutine, public	mod_rhd_phys::rhd_get_pradiation (w, x, ixil, ixol, prad)
	Calculate radiation pressure within ixO^L.

subroutine, public	mod_rhd_phys::rhd_get_ptot (w, x, ixil, ixol, ptot)
	calculates the sum of the gas pressure and max Prad tensor element

subroutine, public	mod_rhd_phys::rhd_get_tgas (w, x, ixil, ixol, tgas)
	Calculates gas temperature.

subroutine, public	mod_rhd_phys::rhd_get_trad (w, x, ixil, ixol, trad)
	Calculates radiation temperature.

double precision function, dimension(ixo^s), public	mod_rhd_phys::rhd_kin_en (w, ixil, ixol, inv_rho)

Variables
logical, public, protected	mod_rhd_phys::rhd_energy = .true.
	Whether an energy equation is used.

logical, public, protected	mod_rhd_phys::rhd_thermal_conduction = .false.
	Whether thermal conduction is added.

type(tc_fluid), allocatable, public	mod_rhd_phys::tc_fl

type(te_fluid), allocatable, public	mod_rhd_phys::te_fl_rhd

logical, public, protected	mod_rhd_phys::rhd_radiative_cooling = .false.
	Whether radiative cooling is added.

type(rc_fluid), allocatable, public	mod_rhd_phys::rc_fl

logical, public, protected	mod_rhd_phys::rhd_dust = .false.
	Whether dust is added.

logical, public, protected	mod_rhd_phys::rhd_viscosity = .false.
	Whether viscosity is added.

logical, public, protected	mod_rhd_phys::rhd_gravity = .false.
	Whether gravity is added.

logical, public, protected	mod_rhd_phys::rhd_particles = .false.
	Whether particles module is added.

logical, public, protected	mod_rhd_phys::rhd_rotating_frame = .false.
	Whether rotating frame is activated.

integer, public, protected	mod_rhd_phys::rhd_n_tracer = 0
	Number of tracer species.

integer, public, protected	mod_rhd_phys::rho_
	Index of the density (in the w array)

integer, dimension(:), allocatable, public, protected	mod_rhd_phys::mom
	Indices of the momentum density.

integer, dimension(:), allocatable, public, protected	mod_rhd_phys::tracer
	Indices of the tracers.

integer, public, protected	mod_rhd_phys::e_
	Index of the energy density (-1 if not present)

integer, public, protected	mod_rhd_phys::p_
	Index of the gas pressure (-1 if not present) should equal e_.

integer, public, protected	mod_rhd_phys::r_e
	Index of the radiation energy.

integer, public, protected	mod_rhd_phys::te_
	Indices of temperature.

integer, public, protected	mod_rhd_phys::tcoff_
	Index of the cutoff temperature for the TRAC method.

double precision, public	mod_rhd_phys::rhd_gamma = 5.d0/3.0d0
	The adiabatic index.

double precision, public	mod_rhd_phys::rhd_adiab = 1.0d0
	The adiabatic constant.

double precision, public, protected	mod_rhd_phys::small_r_e = 0.d0
	The smallest allowed radiation energy.

logical, public, protected	mod_rhd_phys::rhd_trac = .false.
	Whether TRAC method is used.

integer, public, protected	mod_rhd_phys::rhd_trac_type = 1

double precision, public, protected	mod_rhd_phys::he_abundance =0.1d0
	Helium abundance over Hydrogen.

character(len=8), public	mod_rhd_phys::rhd_radiation_formalism = 'fld'
	Formalism to treat radiation.

character(len=8), public	mod_rhd_phys::rhd_pressure = 'Trad'
	In the case of no rhd_energy, how to compute pressure.

logical, public, protected	mod_rhd_phys::rhd_radiation_force = .true.
	Treat radiation fld_Rad_force.

logical, public, protected	mod_rhd_phys::rhd_energy_interact = .true.
	Treat radiation-gas energy interaction.

logical, public, protected	mod_rhd_phys::rhd_radiation_diffusion = .true.
	Treat radiation energy diffusion.

logical, public, protected	mod_rhd_phys::rhd_radiation_advection = .true.
	Treat radiation advection.

logical, public, protected	mod_rhd_phys::rhd_partial_ionization = .false.
	Whether plasma is partially ionized.

double precision, public	mod_rhd_phys::kbmpmua4
	kb/(m_p mu)* 1/a_rad**4,

double precision, public, protected	mod_rhd_phys::h_ion_fr =1d0
	Ionization fraction of H H_ion_fr = H+/(H+ + H)

double precision, public, protected	mod_rhd_phys::he_ion_fr =1d0
	Ionization fraction of He He_ion_fr = (He2+ + He+)/(He2+ + He+ + He)

double precision, public, protected	mod_rhd_phys::he_ion_fr2 =1d0
	Ratio of number He2+ / number He+ + He2+ He_ion_fr2 = He2+/(He2+ + He+)

double precision, public, protected	mod_rhd_phys::rr =1d0

logical, public, protected	mod_rhd_phys::eq_state_units = .true.

Modules

Functions/Subroutines

Variables