35.9 Matter+Radiation Simulations

Here are some simulations that use the RadFLAH (radiation-fluxlimiter-aware hydro) variant of unsplit Hydro. They differ from generic 3T simulations in that there is no phyiscal distinction between the temperatures of electrons and ions; that is, electrons and ions are understood as forming a “matter” component. The “matter” component can have a different temperature from radiation.

35.9.1 Radiation-Inhibited Bondi Accretion

This setup can be used as a 1D spherical version of the radiation-inhibited Bondi accretion problem described by Krumholz et al (2007). The simulation can be set up using the following command:

./setup radflaHD/BondiAccretion -1d -auto +spherical -nxb=16 +mgd mgd_meshgroups=1 \
        species=h1 ManualSpeciesDirectives=True +parallelio +uhd3tr \
        -without-unit=physics/Hydro/HydroMain/unsplit

Note that this setup differs significantly from the one described in Krumholz et al (2007), by using 1D instead of 3D geometry, and by excluding the central region from the domain instead of representing the accreting mass at the origin as a sink (or “star”) particle.

35.9.2 Radiation Blast Wave

This setup can be used to reproduce the Radiation Blast Wave problem described in the Castro II paper by Zhang et al (2011). The simulation can be set up using the following command for using the default MGD solver:

./setup radflaHD/RadBlastWave -1d -auto +spherical -nxb=16 +mgd mgd_meshgroups=1 \
        species=h1 ManualSpeciesDirectives=True +parallelio +uhd3tr

For Case II of Zhang et al (2011), the following variant, using the experimental ExpRelax MGD solver implementation, may give better results:

./setup radflaHD/RadBlastWave -1d -auto +spherical -nxb=16 +mgd mgd_meshgroups=1 \
        species=h1 ManualSpeciesDirectives=True +parallelio +uhd3tr RadTransImpl=ExpRelax

For both cases, the user has to create an appropriate parameter file; several sample *.par files can be found in the Simulation/SimulationMain/radflaHD/RadBlastWave directory.

35.9.3 1D Supernova

A 1D supernova simulation using RadFLAH from which light curves can be generated is included in FLASH4.5. More information will be included here after draft publication of Chatzopoulos and Weide (2017).

Example configuration:

./setup radflaHD/SupernovaRad1D -1d -auto +spherical +mgd mgd_meshgroups=1 -nxb=16 \
 -objdir mesa_sn_2T --maxblocks=1000 +pm4dev +uhd3tr +newMpole \
 species=neut,prot,h1,he3,he4,c12,n14,o16,ne20,mg24,si28,s32,ar36,ca40,ti44,cr48,fe52,fe54,ni56 \
 ManualSpeciesDirectives=True \
 +parallelio -without-unit=physics/Eos/EosMain/multiTemp/Multigamma RadTransImpl=MGD \
 -parfile=flash_working.par