Contrib/CompressibleMixingPhaseChangeFoam

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 ( \right )$

$\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$

The mass transfer coefficient is determined at each wall by assuming that the value of T at the wall is zero:

$k_c=-\frac{D_T}{T_b} \frac{\partial T}{\partial y}|_{y=0}$

where $T_b$ is the bulk value of T and y is the direction normal to the wall.

The Sherwood number is then determined by:

$Sh=\frac{k_c d}{D_T}$

where d is the characteristic dimension.

In order to use the solver, you should add $D_T$ , $T_b$ and d to the transportProperties dictionary.

The solver is based on turbFoam, with the transport of a scalar T having a mass diffusion coefficient $D_T$ :

* Phase change model
* Momentum equation
* Energy equation