Boiling at high temperature wall

Update: August 10, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling

Summary

Model geometry Model geometry

We calculate the boiling of water when the channel wall is heated for 2 seconds. The calculation is performed as a two-dimensional problem with one mesh in the Z direction.

The water flows in from the region "inlet" and out from the region "outlet". The wall1 is adiabatic and the wall2 is heated using the "fixedMultiPhaseHeatFlux condition" in the file "constant/T.liquid" as follows.

wall2
{
    type            fixedMultiPhaseHeatFlux;
    relax           0.5;
    q               uniform 100000;
    phase           "liquid";
    value           uniform 370;
}		

The meshes are as follows, and the number of mesh is 1875.

Mesh Meshes
Mesh Meshes (zoomed)

The calculation result is as follows.

Volume ratio of the gas at the final time (alpha.gas) Volume ratio of the gas at the final time (alpha.gas)
Temperature of the gas at the final time (T.gas) Temperature of the gas at the final time (T.gas)

We can see that boiling is occurring on the high temperature side.

Commands

cp -r $FOAM_TUTORIALS/multiphase/reactingTwoPhaseEulerFoam/RAS/wallBoiling wallBoiling
cd wallBoiling

blockMesh

reactingTwoPhaseEulerFoam

paraFoam

Calculation time

8 minutes 57.66 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz