Blowing air into water column with evaporation

Update: July 28, 2017
OpenFOAM 4.x

Case directory



Model geometry Model geometry

We calculate the mixing of air into the water column for 100 seconds. Air flows in from the region "inlet" at the bottom of the water column at (0, 0.1, 0) m/s and out from the region "outlet". The calculation is performed as a two-dimensional problem with one mesh in the z-direction.

We consider the phase change of water and assume that gaseous water (vapor) exists.

The physical properties of the gaseous and liquid phases and their constituent substances air and H2O are set in the files "constant/phaseProperties", "constant/thermophysicalProperties.gas" and "constant/thermophysicalProperties.liquid". In the file "constant/phaseProperties", the "multiComponentPhaseModel" type is specified for the gas phase and the "purePhaseModel" type is specified for the liquid phase.

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

Mesh Meshes

The calculation result is as follows.

Volume ratio of liquid (alpha.liquid)
Volume ratio of liquid at final time (alpha.liquid) Volume ratio of liquid at final time (alpha.liquid)
Vapor at final time (H2O.gas) Vapor at final time (H2O.gas)


cp -r $FOAM_TUTORIALS/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating bubbleColumnEvaporating
cd bubbleColumnEvaporating


mpirun -np 6 reactingTwoPhaseEulerFoam -parallel


Calculation time

1 hour 7 minutes 59.42 seconds *6 parallel, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz