m_reactive_burn Module Reference

Condensed-phase reactive burn: a pressure-driven programmed-burn source that converts a "reactant" fluid into a "product" fluid on the multi-fluid model (num_fluids=2, chemistry='F'). The two fluids share the same stiffened-gas EOS (gamma, pi_inf) and differ only in their reference energy qv, so the reactant->product conversion releases (qv_reactant - qv_product) per unit mass through the mixture EOS with no explicit energy source. Because the two fluids are mechanically identical the volume-fraction swap is exact (the product volume fraction is the reaction progress), making this a reactive-Euler/ZND detonation model expressed through the diffuse-interface framework. A shock raises the pressure above rburnpign, the reactant burns, and the energy release sustains the shock – a self-propagating condensed-phase detonation. More...

Functions/Subroutines
subroutine, public	s_compute_reactive_burn (rhs_vf, q_cons_vf, q_prim_vf, bounds)
	Add the programmed-burn reaction source to the continuity and volume-fraction RHS.

Detailed Description

Condensed-phase reactive burn: a pressure-driven programmed-burn source that converts a "reactant" fluid into a "product" fluid on the multi-fluid model (num_fluids=2, chemistry='F'). The two fluids share the same stiffened-gas EOS (gamma, pi_inf) and differ only in their reference energy qv, so the reactant->product conversion releases (qv_reactant - qv_product) per unit mass through the mixture EOS with no explicit energy source. Because the two fluids are mechanically identical the volume-fraction swap is exact (the product volume fraction is the reaction progress), making this a reactive-Euler/ZND detonation model expressed through the diffuse-interface framework. A shock raises the pressure above rburnpign, the reactant burns, and the energy release sustains the shock – a self-propagating condensed-phase detonation.

Function/Subroutine Documentation

s_compute_reactive_burn()

subroutine, public m_reactive_burn::s_compute_reactive_burn	(	type(scalar_field), dimension(sys_size), intent(inout)	rhs_vf,
		type(scalar_field), dimension(sys_size), intent(in)	q_cons_vf,
		type(scalar_field), dimension(sys_size), intent(in)	q_prim_vf,
		type(int_bounds_info), dimension(1:3), intent(in)	bounds )

Add the programmed-burn reaction source to the continuity and volume-fraction RHS.

Parameters

rhs_vf	Right-hand-side accumulator (inout)
q_cons_vf	Conserved variables (partial densities live here)
q_prim_vf	Primitive variables (pressure and volume fractions live here)
bounds	Interior cell bounds

Definition at line 355 of file m_reactive_burn.fpp.f90.