MFC
Exascale flow solver
Loading...
Searching...
No Matches
m_phase_change.fpp.f90 File Reference

Contains module m_phase_change. More...

Go to the source code of this file.

Modules

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

Functions/Subroutines

impure subroutine, public m_phase_change::s_relaxation_solver (q_cons_vf)
 Dispatch to the correct relaxation solver. Replaces the procedure pointer, which CCE is breaking on.
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.
subroutine, public m_phase_change::s_infinite_relaxation_k (q_cons_vf)
 Apply pT- or pTg-equilibrium relaxation with mass depletion based on the incoming state conditions.
subroutine m_phase_change::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 m_phase_change::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 m_phase_change::s_tsat (psat, tsat, tsin)
 Find the saturation temperature for a given saturation pressure using a Newton solver.
impure subroutine, public m_phase_change::s_finalize_relaxation_solver_module
 Finalize the phase change module.

Variables

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

Correct the partial densities of the reacting fluids in case one of them is negative but their sum is positive. Inert phases are not corrected at this moment

real(wp), intent(out) mct
type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf
real(wp), intent(inout) rm
integer, intent(in) j
integer, intent(in) k
integer, intent(in) l

Detailed Description

Contains module m_phase_change.

Definition in file m_phase_change.fpp.f90.

Variable Documentation

◆ j

integer, intent(in) j
private

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

◆ k

integer, intent(in) k
private

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

◆ l

integer, intent(in) l
private

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

◆ mct

real(wp), intent(out) mct
private

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

◆ q_cons_vf

type(scalar_field), dimension(sys_size), intent(inout) q_cons_vf
private

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

◆ rm

real(wp), intent(inout) rm
private

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