|
MFC
Exascale flow solver
|
|
MFC
Exascale flow solver
|
Phase transition relaxation solvers for liquid-vapor flows with cavitation and boiling. More...
Functions/Subroutines | |
| impure subroutine, public | s_relaxation_solver (q_cons_vf) |
| This subroutine should dispatch to the correct relaxation solver based some parameter. It replaces the procedure pointer, which CCE is breaking on. | |
| impure subroutine, public | s_initialize_phasechange_module |
| The purpose of this subroutine is to initialize the phase change module by setting the parameters needed for phase change and selecting the phase change module that will be used (pT- or pTg-equilibrium). | |
| subroutine, public | s_infinite_relaxation_k (q_cons_vf) |
| This subroutine is created to activate either the pT- (N fluids) or the pTg-equilibrium (2 fluids for g-equilibrium) model, also considering mass depletion, depending on the incoming state conditions. | |
| subroutine | s_infinite_pt_relaxation_k (j, k, l, mfl, ps, p_infpt, q_cons_vf, rhoe, ts) |
| This auxiliary subroutine is created to activate the pT-equilibrium for N fluids. | |
| subroutine | s_infinite_ptg_relaxation_k (j, k, l, ps, p_infpt, rhoe, q_cons_vf, ts) |
| This auxiliary subroutine is created to activate the pTg-equilibrium 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) |
| This auxiliary subroutine finds the Saturation temperature for a given saturation pressure through a newton solver. | |
| impure subroutine, public | s_finalize_relaxation_solver_module |
| This subroutine finalizes 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 water pressure. | |
| real(wp), parameter | tcr = 385.05_wp + 273.15_wp |
| Critical water temperature. | |
| real(wp), parameter | mixm = 1.0e-8_wp |
| threshold for 'mixture cell'. If Y < mixM, phase change does not happen | |
| 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 |
Phase transition relaxation solvers for liquid-vapor flows with cavitation and boiling.
| impure subroutine, public m_phase_change::s_finalize_relaxation_solver_module |
This subroutine finalizes the phase change module.
Definition at line 1321 of file m_phase_change.fpp.f90.
| 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 ) |
This auxiliary subroutine is created to activate the pT-equilibrium for N fluids.
| j | generic loop iterator for x direction |
| k | generic loop iterator for y direction |
| l | generic loop iterator for z direction |
| MFL | flag that tells whether the fluid is gas (0), liquid (1), or a mixture (2) |
| pS | equilibrium pressure at the interface |
| p_infpT | stiffness for the participating fluids under pT-equilibrium |
| q_cons_vf | Cell-average conservative variables |
| rhoe | mixture energy |
| TS | equilibrium temperature at the interface |
Definition at line 697 of file m_phase_change.fpp.f90.
| 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 ) |
This auxiliary subroutine is created to activate the pTg-equilibrium for N fluids under pT and 2 fluids under pTg-equilibrium. There is a final common p and T during relaxation.
| j | generic loop iterator for x direction |
| k | generic loop iterator for y direction |
| l | generic loop iterator for z direction |
| pS | equilibrium pressure at the interface |
| p_infpT | stiffness for the participating fluids under pT-equilibrium |
| rhoe | mixture energy |
| q_cons_vf | Cell-average conservative variables |
| TS | equilibrium temperature at the interface |
Definition at line 869 of file m_phase_change.fpp.f90.
| subroutine, public m_phase_change::s_infinite_relaxation_k | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_cons_vf | ) |
This subroutine is created to activate either the pT- (N fluids) or the pTg-equilibrium (2 fluids for g-equilibrium) model, also considering mass depletion, depending on the incoming state conditions.
| q_cons_vf | Cell-average conservative variables |
Definition at line 400 of file m_phase_change.fpp.f90.
| impure subroutine, public m_phase_change::s_initialize_phasechange_module |
The purpose of this subroutine is to initialize the phase change module by setting the parameters needed for phase change and selecting the phase change module that will be used (pT- or pTg-equilibrium).
Definition at line 380 of file m_phase_change.fpp.f90.
| impure subroutine, public m_phase_change::s_relaxation_solver | ( | type(scalar_field), dimension(sys_size), intent(inout) | q_cons_vf | ) |
This subroutine should dispatch to the correct relaxation solver based some parameter. It replaces the procedure pointer, which CCE is breaking on.
Definition at line 361 of file m_phase_change.fpp.f90.
| elemental subroutine m_phase_change::s_tsat | ( | real(wp), intent(in) | psat, |
| real(wp), intent(out) | tsat, | ||
| real(wp), intent(in) | tsin ) |
This auxiliary subroutine finds the Saturation temperature for a given saturation pressure through a newton solver.
| pSat | Saturation Pressure |
| TSat | Saturation Temperature |
| TSIn | equilibrium Temperature |
Definition at line 1224 of file m_phase_change.fpp.f90.
| real(wp) m_phase_change::a |
Definition at line 341 of file m_phase_change.fpp.f90.
| real(wp) m_phase_change::b |
Definition at line 341 of file m_phase_change.fpp.f90.
| real(wp) m_phase_change::c |
Definition at line 341 of file m_phase_change.fpp.f90.
| real(wp) m_phase_change::d |
Definition at line 341 of file m_phase_change.fpp.f90.
| integer, parameter m_phase_change::lp = 1 |
index for the liquid phase of the reacting fluid
Definition at line 335 of file m_phase_change.fpp.f90.
| integer, parameter m_phase_change::max_iter = 1e8_wp |
max # of iterations
Definition at line 331 of file m_phase_change.fpp.f90.
| real(wp), parameter m_phase_change::mixm = 1.0e-8_wp |
threshold for 'mixture cell'. If Y < mixM, phase change does not happen
Definition at line 334 of file m_phase_change.fpp.f90.
| real(wp), parameter m_phase_change::pcr = 4.94e7_wp |
Critical water pressure.
Definition at line 332 of file m_phase_change.fpp.f90.
| real(wp), parameter m_phase_change::tcr = 385.05_wp + 273.15_wp |
Critical water temperature.
Definition at line 333 of file m_phase_change.fpp.f90.
| integer, parameter m_phase_change::vp = 2 |
index for the vapor phase of the reacting fluid
Definition at line 336 of file m_phase_change.fpp.f90.