CoreFlows/Documentation/PhysicalModels/TwoPhase/DriftModelPage.dox

   1 /*!
   2 \page DriftModelPage The Drift model
   3
   4
   5 [TOC]
   6 The model consists in the steam mass balance equation together with the mixture mass conservation, the mixture momentum balance and mixture energy balance equations. The main unknowns are the steam mass concentration \f$c_v\f$, the pressure \f$P\f$, the mixture velocity \f$\vec{u}_m\f$, and the common temperature \f$T\f$. The model uses stiffened gas laws \f$p_g(\rho_g,T)\f$ and  \f$p_l(\rho_l,T)\f$ as well as  linearised internal energy law \f$e_k(T)\f$ valid around the saturation points \f$(P=1 bar, T=373K)\f$ or \f$(P=155 bars, T=618K)\f$ depending on the value of the enum \ref pressureEstimate.
   7
   8 The drift model is a system of four nonlinear equations taking the following conservative form
   9 \f{eqnarray*}
  10 \left\{\begin{array}{lll}
  11          \partial_t(\phi \rho_m) &+\nabla\cdot(\phi\rho_m\vec{u}_m)&=0\\
  12          \partial_t(\phi m_g)&+\nabla\cdot\phi(m_g\vec{u}_g)&=\phi\Gamma_g(h_m,\Phi)\\
  13          \partial_t(\phi\rho_m\vec{u}_m)&+\nabla\cdot\phi(m_g\vec{u}_g\otimes\vec{u}_g+ m_l\vec{u}_l\otimes\vec{u}_l)+\vec{\nabla}(p\phi)&=p\vec{\nabla}\phi+\phi\rho_m\vec{g}- K_g\phi m_g||\vec{u}_g||\vec{u}_g- K_l\phi m_l||\vec{u}_l||\vec{u}_l- K_s\delta_s(x)\phi\rho_m||\vec{u}_m||\vec{u}_m\\
  14          \partial_t\phi (\rho_m E_m)&+\nabla\cdot\phi(m_g H_g{}^t\vec{u}_g+m_l H_l{}^t\vec{u}_l)&=\Phi+\phi\rho_m\vec{g}\cdot\vec{u}_m- K_g\phi m_g||\vec{u}_g||^3- K_l\phi m_l||\vec{u}_l||^3- K_s\delta_s(x)\phi\rho_m||\vec{u}_m||^3
  15         \end{array}\right.,
  16 \f}
  17 where the mixture quantities are defined by
  18 \f{eqnarray*}
  19 \rho_m&=&\alpha_g\rho_g+\alpha_l\rho_l\\
  20 \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}\\
  21 E_m&=&\alpha_g\rho_g E_g+\alpha_l\rho_l E_l,
  22 \f}
  23
  24 whereas the quantities associated to each to phase \f$k=g,l\f$ are defined as
  25 - \f$\alpha_k\f$ is the phasic volumic presence rate,
  26 - \f$\rho_k\f$ is the phasic density,
  27 - \f$m_k=\alpha_k\rho_k\f$ is the phasic partial density,
  28 - \f$\vec u_k\f$ the phasic velocity,
  29 - \f$\vec q_k = \rho \vec u\f$ the phasic momentum,
  30 - \f$p\f$ the common phasic pressure,
  31 - \f$e_k\f$ the phasic internal energy,
  32 - \f$E_k = e_k + \frac{||\vec u||^2}{2}\f$ the phasic total energy,
  33 - \f$h_k=e_k+\frac{p}{\rho_k}\f$ the phasic enthalpy
  34 - \f$H_k=h_k+\frac{1}{2}|\vec{u}_k|^2\f$ the phasic total enthalpy
  35 - \f$T\f$ the common absolute temperature,
  36 - \f$\nu_k\f$ the viscosity (\ref DriftModel::setViscosity),
  37 - \f$\lambda_k\f$ the thermal conductivity (\ref DriftModel::setConductivity),
  38 - \f$K_k\f$ the phasic regular friction coefficient (\ref DriftModel::setDragCoeffs),
  39
  40 Geometric and physical source terms are
  41 - \f$\vec g\f$ the gravity vector (\ref DriftModel::setGravity)
  42 - \f$\Phi(\vec x)\f$ the heat power received by the fluid (\ref DriftModel::setHeatPowerField),
  43 - \f$\phi(\vec x)\f$ the volumic porosity field (\ref DriftModel::setPorosityField),
  44 - \f$K_s(\vec x)\f$ the singular friction function, \f$\delta_s(\vec x)\f$ the Dirac delta function with support on the set $s$ (\ref DriftModel::setPressureLossField).
  45
  46 We close the Drift-Model system with a stiffened gas law \f$p = (\gamma_k -1) \rho_k e_k -\gamma_k p_{0k}\f$ for each phase and a linearised enthalpy law \f$h_k(T)\f$ valid around the points \f$(P=1 bar, T=300K)\f$ or \f$(P=155 bars, T=618K)\f$ depending on the value of the enum \ref pressureEstimate.
  47
  48 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.
  49
  50 The constant parameters \f$\lambda_k, \nu_k,\vec g, K_k\f$ and the fields \f$\phi(\vec x),\: \Phi(\vec x),\: K_s(\vec x)\f$ can be set by the user. The default value for \f$\phi\f$ is \f$\phi=1\f$.
  51
  52
  53 To close the system we need a drift correlation for the relative velocity:
  54 \f{eqnarray*}
  55 \vec{u}_r=\vec{u}_g-\vec{u}_l=\vec{f}_r(c_g,\vec{u}_m,\rho_m).
  56 \f}
  57 For the moment the only drift correlation available is \f$\vec{u}_g=\vec{u}_l\f$.
  58
  59 The phase change is modeled using the formula
  60 \f{eqnarray*}
  61  \Gamma_g=\left\{\begin{array}{cc}
  62          \frac{\Phi}{\mathcal{L}}&\textrm{ if } T^{sat}\leq T \textrm{ and } 0<\alpha_g<1\\[1.5ex]
  63          0& \textrm{ otherwise }
  64         \end{array}\right..
  65 \f}
  66
  67 For the moment the boiling temperature \f$T^{sat}\f$ is constant and can be changed using the function DriftModell::setSatTemp.
  68
  69 * The class : \ref DriftModel implements the 4 equation drift model
  70
  71 * \subpage ExampleDriftModelPage "Here are C and Python example scripts using the Drift Model  "
  72
  73
  74 */