m_riemann_solver_hypo_hlld.fpp.f90 File Reference

Contains module m_riemann_solver_hypo_hlld. More...

Go to the source code of this file.

Modules
module	m_riemann_solver_hypo_hlld
	Dual-pass HLLD approximate Riemann solver for hypoelastic flows, with non-conservative interface-velocity coupling.

Functions/Subroutines
integer function	m_riemann_solver_hypo_hlld::f_hlld_wave_zone (s_l, s_lstar, s_m, s_rstar, s_r)
	Classify the position of the interface (xi = 0) in the five-wave HLLD fan: 0 left of S_L, 1..4 for the four inner wedges, 5 right of S_R. Single definition of the wave-fan geometry shared by the flux fold, the NC face-velocity export, and the axisymmetric face-state pick, which must stay consistent. Valid fans only: a fallback face exports the HLL/one-sided traces instead of a fan classification.
subroutine	m_riemann_solver_hypo_hlld::s_hypo_hlld_riemann_solver (ql_prim_rsx_vf, dql_prim_dx_vf, dql_prim_dy_vf, dql_prim_dz_vf, ql_prim_vf, qr_prim_rsx_vf, dqr_prim_dx_vf, dqr_prim_dy_vf, dqr_prim_dz_vf, qr_prim_vf, q_prim_vf, flux_vf, flux_src_vf, flux_gsrc_vf, norm_dir, ix, iy, iz)
	HLLD Riemann solver resolves all 5 waves for the hypoelastic equations: 1 entropy wave, 2 shear stress waves, 2 fast waves.
subroutine	m_riemann_solver_hypo_hlld::s_finalize_nc_iface_vel (nc_iface_vel_vf, norm_dir)
	Copy the hypo interface velocities of a fused-solve pass from Riemann-space buffers to physical-space output arrays, called after the Riemann solver for each sweep direction. The plain routine handles the hat_L-anchored values (nc_iface_vel_rsx; used when hypo_nc_mode_interface is active, and as the hat_L pass of the dual-pass HLLD solve); the _hatR routine handles the hat_R-anchored values (nc_iface_vel_hatR_rsx) of the dual-pass solve. Both routines are emitted from one template so the two passes cannot drift apart.
subroutine	m_riemann_solver_hypo_hlld::s_finalize_nc_iface_vel_hatr (nc_iface_vel_vf, norm_dir)

Variables
real(wp), parameter	m_riemann_solver_hypo_hlld::g_eff_roundoff_factor = 4._wp
	Four-epsilon roundoff band for G_eff: two mixture products, their sum, and the final signed sum. Relative to the operand scale; not a physical stiffness cutoff.

Detailed Description

Contains module m_riemann_solver_hypo_hlld.

Definition in file m_riemann_solver_hypo_hlld.fpp.f90.