MFC: Simulation
High-fidelity multiphase flow simulation
Loading...
Searching...
No Matches
m_rhs.fpp.f90 File Reference

Functions/Subroutines

program __m_rhs_fpp_f90__
 
subroutine s_initialize_rhs_module ()
 The computation of parameters, the allocation of memory, the association of pointers and/or the execution of any other procedures that are necessary to setup the module.
 
subroutine s_compute_rhs (q_cons_vf, q_prim_vf, rhs_vf, pb, rhs_pb, mv, rhs_mv, t_step)
 
subroutine s_compute_advection_source_term (idir, rhs_vf, q_cons_vf, q_prim_vf, flux_src_n_vf)
 

Relaxed pressure, initial partial pressures, function f(p) and its partial

The purpose of this procedure is to infinitely relax the pressures from the internal-energy equations to a unique pressure, from which the corresponding volume fraction of each phase are recomputed. For conservation purpose, this pressure is finally corrected using the mixture-total-energy equation.

Parameters
q_cons_vfCell-average conservative variables
real(kind(0d0)) pres_relax
 
real(kind(0d0)), dimension(num_fluids) pres_k_init
 
real(kind(0d0)) f_pres
 
real(kind(0d0)) df_pres
 
real(kind(0d0)), dimension(num_fluids) rho_k_s
 
real(kind(0d0)), dimension(num_fluids) alpha_rho
 
real(kind(0d0)), dimension(num_fluids) alpha
 
real(kind(0d0)) sum_alpha
 
real(kind(0d0)) rho
 
real(kind(0d0)) dyn_pres
 
real(kind(0d0)) gamma
 
real(kind(0d0)) pi_inf
 
real(kind(0d0)), dimension(2) re
 
integer i
 
integer j
 
integer k
 
integer l
 
integer q
 
integer iter
 Generic loop iterators.
 
integer relax
 Relaxation procedure determination variable.
 
subroutine s_reconstruct_cell_boundary_values (v_vf, vl_x, vl_y, vl_z, vr_x, vr_y, vr_z, norm_dir)
 The purpose of this subroutine is to WENO-reconstruct the left and the right cell-boundary values, including values at the Gaussian quadrature points, from the cell-averaged variables.
 
subroutine s_finalize_rhs_module ()
 Module deallocation and/or disassociation procedures.
 

Function/Subroutine Documentation

◆ __m_rhs_fpp_f90__()

program __m_rhs_fpp_f90__
Here is the call graph for this function:

◆ s_compute_advection_source_term()

subroutine __m_rhs_fpp_f90__::s_compute_advection_source_term ( integer, intent(in) idir,
type(scalar_field), dimension(sys_size), intent(inout) rhs_vf,
type(vector_field), intent(inout) q_cons_vf,
type(vector_field), intent(inout) q_prim_vf,
type(vector_field), intent(inout) flux_src_n_vf )
private
Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_compute_rhs()

subroutine __m_rhs_fpp_f90__::s_compute_rhs ( type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf,
type(scalar_field), dimension(sys_size), intent(inout) q_prim_vf,
type(scalar_field), dimension(sys_size), intent(inout) rhs_vf,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:, 1:), intent(inout) pb,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:, 1:), intent(inout) rhs_pb,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:, 1:), intent(inout) mv,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:, 1:), intent(inout) rhs_mv,
integer, intent(in) t_step )
private
Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_finalize_rhs_module()

subroutine __m_rhs_fpp_f90__::s_finalize_rhs_module
private

Module deallocation and/or disassociation procedures.

Here is the caller graph for this function:

◆ s_initialize_rhs_module()

subroutine __m_rhs_fpp_f90__::s_initialize_rhs_module
private

The computation of parameters, the allocation of memory, the association of pointers and/or the execution of any other procedures that are necessary to setup the module.

Here is the call graph for this function:
Here is the caller graph for this function:

◆ s_reconstruct_cell_boundary_values()

subroutine __m_rhs_fpp_f90__::s_reconstruct_cell_boundary_values ( type(scalar_field), dimension(iv%beg:iv%end), intent(in) v_vf,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vl_x,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vl_y,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vl_z,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vr_x,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vr_y,
real(kind(0d0)), dimension(startx:, starty:, startz:, 1:), intent(inout) vr_z,
integer, intent(in) norm_dir )
private

The purpose of this subroutine is to WENO-reconstruct the left and the right cell-boundary values, including values at the Gaussian quadrature points, from the cell-averaged variables.

Parameters
v_vfCell-average variables
vL_qpLeft WENO-reconstructed, cell-boundary values including the values at the quadrature points, of the cell-average variables
vR_qpRight WENO-reconstructed, cell-boundary values including the values at the quadrature points, of the cell-average variables
norm_dirSplitting coordinate direction
Here is the call graph for this function:
Here is the caller graph for this function:

Variable Documentation

◆ alpha

real(kind(0d0)), dimension(num_fluids) alpha

◆ alpha_rho

real(kind(0d0)), dimension(num_fluids) alpha_rho

◆ df_pres

real(kind(0d0)) df_pres

◆ dyn_pres

real(kind(0d0)) dyn_pres

◆ f_pres

real(kind(0d0)) f_pres

◆ gamma

real(kind(0d0)) gamma

◆ i

integer i

◆ iter

integer iter

Generic loop iterators.

◆ j

integer j

◆ k

integer k

◆ l

integer l

◆ pi_inf

real(kind(0d0)) pi_inf

◆ pres_k_init

real(kind(0d0)), dimension(num_fluids) pres_k_init

◆ pres_relax

real(kind(0d0)) pres_relax

◆ q

integer q

◆ re

real(kind(0d0)), dimension(2) re

◆ relax

integer relax

Relaxation procedure determination variable.

◆ rho

real(kind(0d0)) rho

◆ rho_k_s

real(kind(0d0)), dimension(num_fluids) rho_k_s

◆ sum_alpha

real(kind(0d0)) sum_alpha