Blowing air into water column with evaporation and dissolution

Update: July 28, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporatingDissolving

Summary

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 air and water and assume that liquid air and gaseous water (vapor) exist.

The properties of each phase of gas and liquid, as well as their constituent substances air and water, 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 both the gas and liquid phases.

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)
Air in the liquid phase at final time (air.liquid) Air in the liquid phase at final time (air.liquid)
Water in the liquid phase at final time (water.liquid) Water in the liquid phase at final time (water.liquid)
Air in the gas phase at final time (air.gas) Air in the gas phase at final time (air.gas)
Water in the gas phase at final time (water.gas) Water in the gas phase at final time (water.gas)

Commands

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

blockMesh

setFields
reactingTwoPhaseEulerFoam

paraFoam

Calculation time

2 hours 2 minutes 27.54 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz