Equation of State Notes
EOS Calls
Initialization
advance_timestep
- Step 1.
Define the average expansion at time
and the new * if ``dpdt_factor``* 0 thenIn
makePfromRhoH
, we compute a thermodynamic for the volume discrepancy term using .
end if
- Step 2.
Construct the provisional edge-based advective velocity,
- Step 3.
React the full state and base state through the first time interval of
In
react_state
burner, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
)
if
evolve_base_state = T
thenIn
make_gamma
, we compute using and .
end if
- Step 4.
- Advect the base state, then the full state, through a time interval ofif
use_thermal_diffusion = T
thenIn advance before
make_explicit_thermal
, we compute the coefficients for using .In
enthalpy_advance
update_scal, we compute above thebase_cutoff_density
using and .In
thermal_conduct
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
else
In
enthalpy_advance
update_scal
, we compute above thebase_cutoff_density
using and .In
enthalpy_advance
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
end if
- Step 5.
React the full state through a second time interval of
In
react_state
burner, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (if use_tfromp = T).
if
evolve_base_state
thenIn
make_gamma
, we compute using and .
end if
- Step 6.
- Define a new average expansion rate at timeif
use_thermal_diffusion
thenIn advance before
make_explicit_thermal
, we compute the coefficients for using .
end if
In
make_S
, we compute thermodynamic variables using .
if
dpdt_factor
0 thenIn
makePfromRhoH
, we compute a thermodynamic for the volume discrepancy term using .
end if
- Step 7.
Construct the final edge-based advective velocity,
- Step 8.
- Advect the base state, then the full state, through a time interval ofif
use_thermal_diffusion = T
thenIn
enthalpy_advance
update_scal
, we compute above thebase_cutoff_density
using and .In
advance
beforethermal_conduct
, we compute the coefficients for using .In
thermal_conduct
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
else
In
enthalpy_advance
update_scal
, we compute above thebase_cutoff_density
using and .In
enthalpy_advance
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
end if
- Step 9.
React the full state and base state through a second time interval of
In
react_state
burner, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
if
evolve_base_state = T
thenIn
make_gamma
, we compute using and .
end if
- Step 10.
- Compute
for the final projection.ifmake_explicit_thermal
thenIn
advance
beforemake_explicit_thermal
, we compute the coefficients for using .
end if
In
make_S
, we compute thermodynamic variables using .
- Step 11.
- Update the velocity.if
dpdt_factor
0 thenIn
makePfromRhoH
, we compute a thermodynamic for the volume discrepancy term using .
end if
- Step 12.
Compute a new
make_plotfile
Temperature Usage
advance_timestep
- Step 1.
Define the average expansion at time
and the new- Step 2.
Construct the provisional edge-based advective velocity,
- Step 3.
React the full state and base state through the first time interval of
In
react_state
burner
, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
- Step 4.
- Advect the base state, then the full state, through a time interval ofif
use_thermal_diffusion = T
thenIn
advance
beforemake_explicit_thermal
, we compute the coefficients for using .In
thermal_conduct
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
else
In
enthalpy_advance
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
end if
- Step 5.
React the full state through a second time interval of
In
react_state
burner
, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
- Step 6.
- Define a new average expansion rate at timeif
use_thermal_diffusion = T
thenIn advance before
make_explicit_thermal
, we compute the coefficients for using .
end if
In
make_S
, we compute thermodynamic variables using .
- Step 7.
Construct the final edge-based advective velocity,
- Step 8.
- Advect the base state, then the full state, through a time interval ofif
use_thermal_diffusion = T
thenIn
advance
beforethermal_conduct
, we compute the coefficients for using .In
thermal_conduct
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
else
In
enthalpy_advance
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
end if
- Step 9.
React the full state and base state through a second time interval of
In
react_state
burner
, we compute and for inputs to VODE using .In
react_state
, we compute using (ifuse_tfromp = F
) or and (ifuse_tfromp = T
).
- Step 10.
- Compute
for the final projection.ifmake_explicit_thermal
thenIn
advance
beforemake_explicit_thermal
, we compute the coefficients for using .
end if
In
make_S
, we compute thermodynamic variables using .
- Step 11.
Update the velocity.
- Step 12.
Compute a new