m_derived_variables.fpp.f90 File Reference

Modules
module	m_derived_variables
	This module features subroutines that allow for the derivation of numerous flow variables from the conservative and primitive ones. Currently, the available derived variables include the unadvected volume fraction, specific heat ratio, liquid stiffness, speed of sound, vorticity and the numerical Schlieren function.

Functions/Subroutines
subroutine, public	m_derived_variables::s_initialize_derived_variables_module
	Computation of parameters, allocation procedures, and/or any other tasks needed to properly setup the module.
subroutine, public	m_derived_variables::s_derive_specific_heat_ratio (q_sf)
	This subroutine receives as input the specific heat ratio function, gamma_sf, and derives from it the specific heat ratio. The latter is stored in the derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_derive_liquid_stiffness (q_sf)
	This subroutine admits as inputs the specific heat ratio function and the liquid stiffness function, gamma_sf and pi_inf_sf, respectively. These are used to calculate the values of the liquid stiffness, which are stored in the derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_derive_sound_speed (q_prim_vf, q_sf)
	This subroutine admits as inputs the primitive variables, the density, the specific heat ratio function and liquid stiffness function. It then computes from those variables the values of the speed of sound, which are stored in the derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_derive_flux_limiter (i, q_prim_vf, q_sf)
	This subroutine derives the flux_limiter at cell boundary i+1/2. This is an approximation because the velocity used to determine the upwind direction is the velocity at the cell center i instead of the contact velocity at the cell boundary from the Riemann solver.
subroutine	m_derived_variables::s_solve_linear_system (a, b, sol, ndim)
	Computes the solution to the linear system Ax=b w/ sol = x.
subroutine, public	m_derived_variables::s_derive_vorticity_component (i, q_prim_vf, q_sf)
	This subroutine receives as inputs the indicator of the component of the vorticity that should be outputted and the primitive variables. From those inputs, it proceeds to calculate values of the desired vorticity component, which are subsequently stored in derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_derive_qm (q_prim_vf, q_sf)
	This subroutine gets as inputs the primitive variables. From those inputs, it proceeds to calculate the value of the Q_M function, which are subsequently stored in the derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_derive_numerical_schlieren_function (q_cons_vf, q_sf)
	This subroutine gets as inputs the conservative variables and density. From those inputs, it proceeds to calculate the values of the numerical Schlieren function, which are subsequently stored in the derived flow quantity storage variable, q_sf.
subroutine, public	m_derived_variables::s_finalize_derived_variables_module
	Deallocation procedures for the module.

Variables
real(wp), dimension(:, :, :), allocatable	m_derived_variables::gm_rho_sf
	Gradient magnitude (gm) of the density for each cell of the computational sub-domain. This variable is employed in the calculation of the numerical Schlieren function.
integer, private	m_derived_variables::flg
	Flagging (flg) variable used to annotate the dimensionality of the dataset that is undergoing the post-process. A flag value of 1 indicates that the dataset is 3D, while a flag value of 0 indicates that it is not. This flg variable is necessary to avoid cycling through the third dimension of the flow variable(s) when the simulation is not 3D and the size of the buffer is non-zero. Note that a similar procedure does not have to be applied to the second dimension since in 1D, the buffer size is always zero.
Finite-difference (fd) coefficients in x-, y- and z-coordinate directions.
Note that because sufficient boundary information is available for all the active coordinate directions, the centered family of the finite-difference schemes is used.
real(wp), dimension(:, :), allocatable, public	m_derived_variables::fd_coeff_x
real(wp), dimension(:, :), allocatable, public	m_derived_variables::fd_coeff_y
real(wp), dimension(:, :), allocatable, public	m_derived_variables::fd_coeff_z