Difference between revisions of "Contrib/CompressibleMixingPhaseChangeFoam"

Revision as of 19:58, 28 December 2012

Solver for two fluids with phase change (for example - water <---> steam), pressure and temperature density dependence

Model Equations

=== Equation of state ===
=== Liquid volume transport ===

Let us consider transport of liquid (heavy phase) volume fraction $\alpha_l$ : $\frac{\partial \alpha_l \rho_l}{\partial t} + \nabla \cdot \left ( \alpha_l \rho_l \mathbf{U} \right ) = \dot m_l$

$\frac{\partial{T}}{\partial t} + \nabla \cdot \left(\mathbf{U} T\right) - \nabla \cdot \left( (D_T + \nu_{turb}/Sc_{turb}) \nabla T \right) = 0$

* Phase change model
* Momentum equation
* Energy equation

Solver sources and tutorials located here

@@ Line 2: / Line 2: @@
 == Model Equations ==
-  * Equation of state
+  === Equation of state ===
-  * Liquid volume transport
+  === Liquid volume transport ===
 Let us consider transport of liquid (heavy phase) volume fraction <math>\alpha_l</math>:
 <math>
 \frac{\partial \alpha_l \rho_l}{\partial t} + \nabla \cdot
@@ Line 19: / Line 18: @@
 \frac{\partial{T}}{\partial t} + \nabla \cdot \left(\mathbf{U} T\right) - \nabla \cdot \left( (D_T + \nu_{turb}/Sc_{turb}) \nabla T \right) = 0
 </math>
-The mass transfer coefficient is determined at each wall by assuming that the value of ''T'' at the wall is zero:
-<math>k_c=-\frac{D_T}{T_b} \frac{\partial T}{\partial y}|_{y=0}</math>
-where <math>T_b</math> is the bulk value of ''T'' and ''y'' is the direction normal to the wall.
-The Sherwood number is then determined by:
-<math>Sh=\frac{k_c d}{D_T}</math>
-where ''d'' is the characteristic dimension.
-In order to use the solver, you should add <math>D_T</math>, <math>T_b</math> and ''d'' to the transportProperties dictionary.
-The solver is based on turbFoam, with the transport of a scalar ''T'' having a mass diffusion coefficient <math>D_T</math>: