|
MFC
Exascale flow solver
|
|
MFC
Exascale flow solver
|
Contains module m_variables_conversion. More...
Go to the source code of this file.
Modules | |
| module | m_variables_conversion |
| Conservative-to-primitive variable conversion, mixture property evaluation, and pressure computation. | |
Functions/Subroutines | |
| subroutine, public | m_variables_conversion::s_convert_to_mixture_variables (q_vf, i, j, k, rho, gamma, pi_inf, qv, re_k, g_k, g) |
| Dispatch to the s_convert_mixture_to_mixture_variables and s_convert_species_to_mixture_variables subroutines. Replaces a procedure pointer. | |
| subroutine, public | m_variables_conversion::s_compute_pressure (energy, alf, dyn_p, pi_inf, gamma, rho, qv, rhoyks, pres, t, e_e_in, pres_mag) |
| Compute the pressure from the appropriate equation of state. | |
| subroutine, public | m_variables_conversion::s_convert_mixture_to_mixture_variables (q_vf, i, j, k, rho, gamma, pi_inf, qv) |
| Convert mixture variables to density, gamma, pi_inf, and qv for the gamma/pi_inf model. Given conservative or primitive variables, transfers the density, specific heat ratio function and the liquid stiffness function from q_vf to rho, gamma and pi_inf. | |
| subroutine, public | m_variables_conversion::s_convert_species_to_mixture_variables (q_vf, k, l, r, rho, gamma, pi_inf, qv, re_k, g_k, g) |
| Convert species volume fractions and partial densities to mixture density, gamma, pi_inf, and qv. Given conservative or primitive variables, computes the density, the specific heat ratio function and the liquid stiffness function from q_vf and stores the results into rho, gamma and pi_inf. | |
| subroutine, public | m_variables_conversion::s_convert_species_to_mixture_variables_kernel (rho_k, gamma_k, pi_inf_k, qv_k, alpha_k, alpha_rho_k, re_k, g_k, g) |
| Host- and device-callable conversion kernel for species and mixture variables. | |
| impure subroutine, public | m_variables_conversion::s_initialize_variables_conversion_module (store_mixture_fields, enforce_density_floor, preserve_qbmm_number, lagrange_beta_index) |
| Initialize the variables conversion module. | |
| subroutine, public | m_variables_conversion::s_initialize_mv (qk_cons_vf, mv) |
| Initialize bubble mass-vapor values at quadrature nodes from the conserved moment statistics. | |
| subroutine, public | m_variables_conversion::s_initialize_pb (qk_cons_vf, mv, pb) |
| Initialize bubble internal pressures at quadrature nodes using isothermal relations from the Preston model. | |
| subroutine, public | m_variables_conversion::s_convert_conservative_to_primitive_variables (qk_cons_vf, q_t_sf, qk_prim_vf, ibounds) |
| Convert conserved variables (rho*alpha, rho*u, E, alpha) to primitives (rho, u, p, alpha). Conversion depends on model_eqns: each model has different variable sets and EOS. | |
| impure subroutine, public | m_variables_conversion::s_convert_primitive_to_conservative_variables (q_prim_vf, q_cons_vf) |
| Convert primitives (rho, u, p, alpha) to conserved variables (rho*alpha, rho*u, E, alpha). | |
| subroutine, public | m_variables_conversion::s_convert_primitive_to_flux_variables (qk_prim_vf, fk_vf, fk_src_vf, is1, is2, is3, s2b, s3b, dir_idx_in, dir_flg_in, hll_u_interface_in) |
| Convert primitive variables to Eulerian flux variables. | |
| subroutine, public | m_variables_conversion::s_compute_species_fraction (q_vf, k, l, r, alpha_rho_k, alpha_k) |
| Compute partial densities and volume fractions. | |
| impure subroutine, public | m_variables_conversion::s_finalize_variables_conversion_module () |
| Deallocate fluid property arrays and post-processing fields allocated during module initialization. | |
| subroutine, public | m_variables_conversion::s_compute_mixture_coefficients (alpha_rho_k, alpha_k, rho_k, gamma_k, pi_inf_k, qv_k) |
| Mixture coefficients of one state. Under bubbles_euler with num_fluids == 1 the sole advection slot aliases the void fraction (eqn_idxalf == eqn_idxadvend), so alpha is not a composition there and the coefficients are the liquid's. Clipping stays with callers; it differs between solvers and cannot coincide with that case, as mpp_lim requires num_fluids > 1. | |
| subroutine, public | m_variables_conversion::s_compute_mixture_coefficients_dt (dalpha_rho_dt, dadv_dt, alpha_rho, adv, drho_dt, dgamma_dt, dpi_inf_dt, dqv_dt) |
| Time derivative of the mixture coefficients, mirroring s_compute_mixture_coefficients. | |
| subroutine, public | m_variables_conversion::s_compute_energy (pres, alpha_rho_k, alpha_k, vel_sum, e) |
| Total energy per unit volume, thermodynamic terms only. Callers add magnetic and elastic energy, which are not equation-of-state terms. The chemistry and relativistic branches use a different relation and stay open-coded. | |
| subroutine | m_variables_conversion::s_reference_curve (rho, i, p_ref, e_ref, dp_drho, de_drho, g0, dg0) |
| The reference curve of a state-dependent EOS at rho: p_ref, e_ref, their d/drho, and Gamma_G with its d/drho. A new family adds one case here and nothing else. | |
| logical function, public | m_variables_conversion::f_is_state_dependent (i) |
| Whether the EOS of fluid i is a family whose coefficients vary with density. | |
| logical function | m_variables_conversion::f_has_isentropic_reference (i) |
| True when fluid i's reference curve is itself an isentrope (de_ref = -p_ref d(1/rho), which holds for JWL and Vinet but not for the Mie-Gruneisen Hugoniot) and its Gruneisen coefficient is constant. Those two together make the isentrope through any state closed-form, so it never has to be integrated. | |
| impure real(wp) function | m_variables_conversion::f_hugoniot_compression_limit (c0, s, s2, s3) |
| The largest compression a cubic Hugoniot fit can represent. mu(u_p) = u_p/(u_s - u_p) rises, peaks where c0 = s2 u_p^2 + 2 s3 u_p^3, and falls after; only the rising branch is a physical shock. Returns a huge value for the linear fit, which never turns over. Host-side: called once per fluid at initialization. | |
| subroutine, public | m_variables_conversion::s_eos_coefficients (rho, i, gamma, pi_inf, dpi, dgamma) |
| Gamma, Pi, dPi/drho and dGamma/drho of fluid i at density rho, the coefficients of rho e = Gamma p + Pi(rho). Stiffened and ideal gas keep the constants resolved at init, bit for bit. | |
| real(wp) function, public | m_variables_conversion::f_isentrope_exponent (gamma) |
| Exponent of the stiffened-gas isentrope p + B = const rho**n. Precomputed per fluid as isentrope_n. | |
| real(wp) function, public | m_variables_conversion::f_isentrope_pressure (pi_inf, gamma) |
| Reference pressure of that isentrope. Precomputed per fluid as isentrope_B. | |
| real(wp) function, public | m_variables_conversion::f_sg_thermal (pres, rho_or_t, n, b, cv) |
| Stiffened-gas thermal law p + B = (n - 1)*cv*rho*T. Pass rho to get T, or T to get rho. | |
| real(wp) function, public | m_variables_conversion::f_mixture_temperature (alpha_rho_k, pres, gamma_k, pi_inf_k) |
| Thermal-equilibrium mixture temperature for stiffened gas, from primitives. Algebraically identical to the conservative form in m_phase_change's s_infinite_pt_relaxation_k, T = (rho*e + p - sum(alpha_rho_i*qv_i)) / sum(alpha_rho_i*cv_i*n_i), because rho*e = gamma_mix*p + pi_inf_mix + sum(alpha_rho_i*qv_i) in MFC's stored variables. | |
| subroutine, public | m_variables_conversion::s_phase_coefficients (alpha_rho, alpha, i, rho, gamma, pi_inf, dpi, dgamma) |
| Coefficients of phase i at its own density alpha_rho/alpha: the per-cell dispatch when some fluid's EOS is state dependent, the constants resolved at init otherwise (bit for bit). | |
| real(wp) function | m_variables_conversion::f_c2_from_coefficients (rho, pres, gamma, pi_inf, dpi, dgamma) |
| c^2 = [((Gamma + 1) p + Pi)/rho - dPi/drho - p dGamma/drho]/Gamma, the frozen speed of one phase. | |
| subroutine | m_variables_conversion::s_phase_c2 (rho, pres, i, c2) |
| Frozen sound speed squared of one phase at (rho, p) from its own coefficients. These helpers are subroutines, not functions: a device function that calls a device subroutine is a pattern no other backend-tested code in MFC uses. | |
| subroutine | m_variables_conversion::s_ode_slope (kind, i, x, y, dydx) |
| Slope of the ODE kind for fluid i: dp/drho = c^2 along an isentrope (x = rho, y = p), or the reference temperature dT/dV = (de_ref/dV + p_ref)/c_v - Gamma_G T/V (x = V, y = T), the Maxwell relation applied to e = e_ref + c_v (T - T_ref). | |
| subroutine | m_variables_conversion::s_rk4 (kind, i, x0, y0, x1, y) |
| Fixed-step classical RK4 for the ODE kind from (x0, y0) to x1. | |
| subroutine, public | m_variables_conversion::s_phase_pressure_on_isentrope (pres, rho, xi, i, p_isen) |
| Pressure of phase i after the isentropic density change rho -> xi rho: closed form for the constant-coefficient families, integrated for a state-dependent EOS (the star states it serves are close to rho). | |
| subroutine, public | m_variables_conversion::s_phase_temperature (rho, pres, i, t) |
| Temperature of phase i at (rho, p): the stiffened-gas relation, or T_ref(rho) + (e - e_ref)/c_v. | |
| subroutine, public | m_variables_conversion::s_phase_density_on_isentrope (i, rho_from, p_from, p_to, rho_to, c2_to) |
| Density of phase i on the isentrope through (rho_from, p_from) at p_to, and c^2 there: Newton on the pressure integrator, whose slope is c^2. The relaxation's own Newton wraps this, so a few steps suffice. | |
| subroutine, public | m_variables_conversion::s_phase_internal_energy (pres, alpha, alpha_rho, i, e_phase) |
| Internal energy per unit volume of phase i at pressure pres: alpha (Gamma p + Pi) + alpha_rho qv, with the coefficients at the phase's own density. | |
| subroutine, public | m_variables_conversion::s_phase_bulk_modulus (pres, alpha, alpha_rho, i, blkmod) |
| Bulk modulus rho c^2 of phase i at pressure pres: f_bulk_modulus for a constant-coefficient fluid, bit for bit, minus the reference-curve terms rho (dPi/drho + p dGamma/drho)/Gamma otherwise. | |
| real(wp) function, public | m_variables_conversion::f_elastic_energy (tau, g, is_shear) |
| Elastic strain energy of one stress component, doubled for a shear component: the tensor stores it once, the energy counts both off-diagonal entries. Zero without a shear modulus. | |
| real(wp) function, public | m_variables_conversion::f_hypoelastic_energy (q_cons_vf, j, k, l, rho, g) |
| Hypoelastic strain energy at one cell, summed over the stress components. | |
| real(wp) function, public | m_variables_conversion::f_pressure (e_int, gamma, pi_inf, qv) |
| Pressure of a stiffened gas from its internal energy density - the inverse of s_compute_energy. Callers subtract the kinetic, magnetic and elastic energy first; none of those are equation-of-state terms. | |
| real(wp) function, public | m_variables_conversion::f_bulk_modulus (pres, gamma, pi_inf) |
| Isentropic bulk modulus. Takes coefficients rather than a fluid index, so a mixture - whose effective gamma and pi_inf come from s_compute_mixture_coefficients - is the same call as a single fluid. Elastic callers add their own shear term. | |
| real(wp) function, public | m_variables_conversion::f_relativistic_enthalpy (pres, rho, gamma) |
| Relativistic specific enthalpy, h = 1 + (Gamma + 1)p/rho. Ideal gas only: the stiffness does not appear, so a fluid with a nonzero pi_inf is not represented here (the validator refuses that combination). | |
| subroutine, public | m_variables_conversion::s_compute_speed_of_sound (pres, rho, gamma, pi_inf, adv, c, alpha_rho) |
| Speed of sound of a thermodynamic state. Enthalpy is not an argument: for a real state H, |u|^2 and qv all cancel out of c^2 = ((Gamma + 1)p + Pi)/(Gamma rho). Averaged states, whose enthalpy is a free input, use the _avg variant. | |
| subroutine, public | m_variables_conversion::s_compute_speed_of_sound_avg (pres, rho, gamma, pi_inf, qv, vel_sum, h, c_c, adv, c, alpha_rho) |
| Speed of sound of an interface-averaged state. An average of two states is not a state - its enthalpy is not the one its pressure and density imply - so the caller supplies H, |u|^2 and qv. Only the enthalpy-reading branches differ from s_compute_speed_of_sound; keep the condition below in step with the branch list there. | |
| subroutine, public | m_variables_conversion::s_compute_fast_magnetosonic_speed (rho, c, b, norm, c_fast, h) |
| Compute the fast magnetosonic wave speed from the sound speed, density, and magnetic field components. | |
Variables | |
| real(wp), dimension(:), allocatable | m_variables_conversion::gs_vc |
| integer, dimension(:), allocatable | m_variables_conversion::bubrs_vc |
| real(wp), dimension(:,:), allocatable | m_variables_conversion::res_vc |
| integer | m_variables_conversion::is1b |
| integer | m_variables_conversion::is2b |
| integer | m_variables_conversion::is3b |
| integer | m_variables_conversion::is1e |
| integer | m_variables_conversion::is2e |
| integer | m_variables_conversion::is3e |
| logical | m_variables_conversion::enforce_density_floor_vc = .false. |
| logical | m_variables_conversion::preserve_qbmm_number_vc = .false. |
| integer | m_variables_conversion::lagrange_beta_index_vc = 0 |
| real(wp), dimension(:,:,:), allocatable, public | m_variables_conversion::rho_sf |
| Scalar density function. | |
| real(wp), dimension(:,:,:), allocatable, public | m_variables_conversion::gamma_sf |
| Scalar sp. heat ratio function. | |
| real(wp), dimension(:,:,:), allocatable, public | m_variables_conversion::pi_inf_sf |
| Scalar liquid stiffness function. | |
Contains module m_variables_conversion.
Definition in file m_variables_conversion.fpp.f90.