CoreFlows/Documentation/PhysicalModels/TwoPhase/FiveEqPage.md

   1 The five equation two-fluid model
   2 =================================
   3
   4 The model consists in the phasic mass and momentum balance equations and one mixture total energy balance equation.
   5
   6 The main unknowns are $\alpha$, $P$, $\vec{u}_g$, $\vec{u}_l$ and $T=T_g=T_l$.
   7
   8 The model uses stiffened gas laws $p_g(\rho_g,T)$ and  $p_l(\rho_l,T)$ and linearised internal energy laws $e_g(T)$ and  $e_l(T)$ fitted by either around 1 bar and 373K or around 155 bars and 618K (see \ref pressureEstimate).
   9
  10 $$
  11 \left\{
  12 \begin{array}{ccll}
  13  \frac{\partial m_g}{\partial t}& +& \nabla \cdot \vec{q}_g &= \Gamma_g(h_g,\Phi),\\[1.5ex]
  14 \frac{\partial m_l}{\partial t} &+ &\nabla \cdot \vec{q}_l &= \Gamma_l(h_l,\Phi),\\[1.5ex]
  15 \frac{\partial \vec{q}_g}{\partial t}& +& \nabla \cdot (\vec{q}_g\otimes\frac{\vec{q}_g}{m_g})+ \alpha_g \nabla p&\\[1.5ex]
  16  &+&\Delta p \nabla \alpha_g -\nu_g(\Delta \frac{\vec{q}_g}{m_g}) &= m_g\vec{g}-K\rho_m||\vec{u}_g-\vec{u}_l||(\vec{u}_g-\vec{u}_l)-K_s\delta(x)m_g||\vec{u}_g||\vec{u}_g\\[1.5ex]
  17 \frac{\partial \vec{q}_l}{\partial t}& +& \nabla \cdot (\vec{q}_l\otimes\frac{\vec{q}_l}{m_l})+ \alpha_l \nabla p&\\[1.5ex]
  18 &+&\Delta p \nabla \alpha_l -\nu_l(\Delta \frac{\vec{q}_l}{m_l}) &= m_l\vec{g}-K\rho_m||\vec{u}_l-\vec{u}_g||(\vec{u}_l-\vec{u}_g)-K_s\delta(x)m_l||\vec{u}_l||\vec{u}_l\\[1.5ex]
  19 \partial_t\rho_m E_m&+&\nabla\cdot(\alpha_g\rho_g H_g{}^t\vec{u}_g+\alpha_l\rho_l H_l{}^t\vec{u}_l)&=\Phi+\rho\vec{g}\cdot\vec{u}_m-K\rho_m||\vec{u}_g-\vec{u}_l||^3-K_s\delta(x)(m_g||\vec{u}_g||^3+m_l||\vec{u}||^3)
  20 \end{array}
  21 \right. \nonumber
  22 $$
  23 where the mixture quantities are defined by
  24 $$
  25 \begin{array}{lll}
  26 \rho_m&=&\alpha_g\rho_g+\alpha_l\rho_l\\
  27 \vec{u}_m&=&\frac{\alpha_g\rho_g\vec{u}_g+\alpha_l\rho_l\vec{u}_l}{\alpha_g\rho_g+\alpha_l\rho_l}\\
  28 E_m&=&\alpha_g\rho_g E_g+\alpha_l\rho_l E_l,
  29 \end{array}
  30 $$
  31
  32 whereas the quantities associated to each to phase $k=g,l$ are defined as
  33 - $\alpha_k$ is the phasic volumic presence rate,
  34 - $\rho_k$ is the phasic density,
  35 - $m_k=\alpha_k\rho_k$ is the phasic partial density,
  36 - $\vec u_k$ the phasic velocity,
  37 - $\vec q_k = \rho \vec u$ the phasic momentum,
  38 - $p$ the common phasic pressure,
  39 - $e_k$ the phasic internal energy,
  40 - $E_k = e_k + \frac{||\vec u||^2}{2}$ the phasic total energy,
  41 - $h_k=e_k+\frac{p}{\rho_k}$ the phasic enthalpy
  42 - $H_k=h_k+\frac{1}{2}|\vec{u}_k|^2$ the phasic total enthalpy
  43 - $T$ the common absolute temperature,
  44 - $\nu_k$ the viscosity ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setViscosity),
  45 - $\lambda_k$ the thermal conductivity ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setConductivity),
  46 - $K$ the interphase friction coefficient ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setDragCoeffs),
  47
  48 Geometric and physical source terms are
  49 - $\vec g$ the gravity vector ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setGravity)
  50 - $\Phi(\vec x)$ the heat power received by the fluid ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setHeatPowerField),
  51 - $K_s(\vec x)$ the singular friction function, $\delta_s(\vec x)$ the Dirac delta function with support on the set $s$ ([FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx)::setPressureLossField).
  52
  53 We close the Drift-Model system with a stiffened gas law $p = (\gamma_k -1) \rho_k e_k -\gamma_k p_{0k}$ for each phase and a linearised internal energy law $h_k(T)$ valid around the points $(P=1 bar, T=300K)$ or $(P=155 bars, T=618K)$ depending on the value of the enum \ref pressureEstimate.
  54
  55 For the sake of simplicity, for the moment we consider constant viscosity and conductivity, and neglect the contribution of viscous forces in the energy equation. The constant parameters $\lambda_k, \nu_k,\vec g, K_k$ and the fields $\phi(\vec x),\: \Phi(\vec x),\: K_s(\vec x)$ can be set by the user. The default value for $\phi$ is $\phi=1$.
  56
  57 The phase change is modeled using the formula
  58
  59 $$
  60  \Gamma_g=\left\{\begin{array}{cc}
  61          \frac{\Phi}{\mathcal{L}}&\textrm{ if } T^{sat}\leq T \textrm{ and } 0<\alpha_g<1\\[1.5ex]
  62          0& \textrm{ otherwise }
  63         \end{array}\right..
  64 $$
  65
  66 The parameters $\lambda_k, \nu_k,\vec g, K$ and $\Phi$ can be set by the user.
  67
  68 The class : [FiveEqsTwoFluid](../../../Models/inc/FiveEqsTwoFluid.hxx) implements the equal temperature two fluid model
  69
  70 \subpage Example5EqPage "Here are C and Python example scripts using the five equation two-fluid model "
  71