MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_phase_change Module Reference

Phase transition relaxation solvers for liquid-vapor flows with cavitation and boiling. More...

Functions/Subroutines

impure subroutine, public s_relaxation_solver (q_cons_vf)
 Dispatch to the correct relaxation solver. Replaces the procedure pointer, which CCE is breaking on.
impure subroutine, public s_initialize_phasechange_module
 Initialize the phase change module by setting saturation curve coefficients for pT- or pTg-equilibrium.
subroutine, public s_infinite_relaxation_k (q_cons_vf)
 Apply pT- or pTg-equilibrium relaxation with mass depletion based on the incoming state conditions.
subroutine s_infinite_pt_relaxation_k (j, k, l, mfl, ps, p_infpt, q_cons_vf, rhoe, ts)
 Apply pT-equilibrium relaxation for N fluids.
subroutine s_infinite_ptg_relaxation_k (j, k, l, ps, p_infpt, rhoe, q_cons_vf, ts)
 Apply pTg-equilibrium relaxation for N fluids under pT and 2 fluids under pTg-equilibrium. There is a final common p and T during relaxation.
elemental subroutine s_tsat (psat, tsat, tsin)
 Find the saturation temperature for a given saturation pressure using a Newton solver.
impure subroutine, public s_finalize_relaxation_solver_module
 Finalize the phase change module.

Variables

Parameters for the first order transition phase change
integer, parameter max_iter = 1e8_wp
 max # of iterations
real(wp), parameter pcr = 4.94e7_wp
 Critical pressure of water [Pa].
real(wp), parameter tcr = 385.05_wp + 273.15_wp
 Critical temperature of water [K].
real(wp), parameter mixm = 1.0e-8_wp
 Mixture mass fraction threshold for triggering phase change.
integer, parameter lp = 1
 index for the liquid phase of the reacting fluid
integer, parameter vp = 2
 index for the vapor phase of the reacting fluid
Gibbs free energy phase change parameters
real(wp) a
real(wp) b
real(wp) c
real(wp) d

Detailed Description

Phase transition relaxation solvers for liquid-vapor flows with cavitation and boiling.

Function/Subroutine Documentation

◆ s_finalize_relaxation_solver_module()

impure subroutine, public m_phase_change::s_finalize_relaxation_solver_module

Finalize the phase change module.

Definition at line 1179 of file m_phase_change.fpp.f90.

Here is the caller graph for this function:

◆ s_infinite_pt_relaxation_k()

subroutine m_phase_change::s_infinite_pt_relaxation_k ( integer, intent(in) j,
integer, intent(in) k,
integer, intent(in) l,
integer, intent(in) mfl,
real(wp), intent(out) ps,
real(wp), dimension(1:), intent(out) p_infpt,
type(scalar_field), dimension(sys_size), intent(in) q_cons_vf,
real(wp), intent(in) rhoe,
real(wp), intent(out) ts )
private

Apply pT-equilibrium relaxation for N fluids.

Parameters
MFLflag: 0=gas, 1=liquid, 2=mixture

Definition at line 657 of file m_phase_change.fpp.f90.

Here is the caller graph for this function:

◆ s_infinite_ptg_relaxation_k()

subroutine m_phase_change::s_infinite_ptg_relaxation_k ( integer, intent(in) j,
integer, intent(in) k,
integer, intent(in) l,
real(wp), intent(inout) ps,
real(wp), dimension(1:), intent(in) p_infpt,
real(wp), intent(in) rhoe,
type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf,
real(wp), intent(inout) ts )
private

Apply pTg-equilibrium relaxation for N fluids under pT and 2 fluids under pTg-equilibrium. There is a final common p and T during relaxation.

Generic loop iterators

Definition at line 808 of file m_phase_change.fpp.f90.

Here is the caller graph for this function:

◆ s_infinite_relaxation_k()

subroutine, public m_phase_change::s_infinite_relaxation_k ( type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf)

Apply pT- or pTg-equilibrium relaxation with mass depletion based on the incoming state conditions.

Generic loop iterators

Definition at line 403 of file m_phase_change.fpp.f90.

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

◆ s_initialize_phasechange_module()

impure subroutine, public m_phase_change::s_initialize_phasechange_module

Initialize the phase change module by setting saturation curve coefficients for pT- or pTg-equilibrium.

Definition at line 389 of file m_phase_change.fpp.f90.

Here is the caller graph for this function:

◆ s_relaxation_solver()

impure subroutine, public m_phase_change::s_relaxation_solver ( type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf)

Dispatch to the correct relaxation solver. Replaces the procedure pointer, which CCE is breaking on.

Definition at line 373 of file m_phase_change.fpp.f90.

◆ s_tsat()

elemental subroutine m_phase_change::s_tsat ( real(wp), intent(in) psat,
real(wp), intent(out) tsat,
real(wp), intent(in) tsin )
private

Find the saturation temperature for a given saturation pressure using a Newton solver.

Definition at line 1094 of file m_phase_change.fpp.f90.

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

Variable Documentation

◆ a

real(wp) m_phase_change::a
private

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

◆ b

real(wp) m_phase_change::b
private

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

◆ c

real(wp) m_phase_change::c
private

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

◆ d

real(wp) m_phase_change::d
private

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

◆ lp

integer, parameter m_phase_change::lp = 1
private

index for the liquid phase of the reacting fluid

Definition at line 349 of file m_phase_change.fpp.f90.

◆ max_iter

integer, parameter m_phase_change::max_iter = 1e8_wp
private

max # of iterations

Definition at line 345 of file m_phase_change.fpp.f90.

◆ mixm

real(wp), parameter m_phase_change::mixm = 1.0e-8_wp
private

Mixture mass fraction threshold for triggering phase change.

Definition at line 348 of file m_phase_change.fpp.f90.

◆ pcr

real(wp), parameter m_phase_change::pcr = 4.94e7_wp
private

Critical pressure of water [Pa].

Definition at line 346 of file m_phase_change.fpp.f90.

◆ tcr

real(wp), parameter m_phase_change::tcr = 385.05_wp + 273.15_wp
private

Critical temperature of water [K].

Definition at line 347 of file m_phase_change.fpp.f90.

◆ vp

integer, parameter m_phase_change::vp = 2
private

index for the vapor phase of the reacting fluid

Definition at line 350 of file m_phase_change.fpp.f90.