Cavitation at the narrow flow path (3D, LES)

Update: July 1, 2017
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/multiphase/cavitatingFoam/les/throttle3D

Summary

We set a total pressure of 300×105 in the region "inlet" and a static pressure of 100×105 in the region "outlet". The flow of water will be driven from the "inlet" to the "outlet" by the pressure difference and we calculate the flow for 0.0001 seconds. The region "walls" are set to the no-slip condition.

Model Model

The LES model is used for turbulence calculations.

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

Mesh (entire) Mesh (entire)
Meshes (narrow part) Meshes (narrow part)

The calculation result is as follows.

Flow velocity at 0.0001 sec (U) Flow velocity at 0.0001 sec (U)
Volume ratio of vapor near narrow part at 0.00005 sec (alpha.vapour) Volume ratio of vapor near narrow part at 0.00005 sec (alpha.vapour)
Flow velocity near narrow part at 0.00005 sec (U) Flow velocity near narrow part at 0.00005 sec (U)

Commands

cp -r $FOAM_TUTORIALS/multiphase/cavitatingFoam/les/throttle3D throttle3D
cd throttle3D

cp -r 0.orig 0

blockMesh

cp system/topoSetDict.1 system/topoSetDict
topoSet
refineMesh -dict system/refineMeshDict -overwrite

cp system/topoSetDict.2 system/topoSetDict
topoSet
refineMesh -dict system/refineMeshDict -overwrite

cp system/topoSetDict.3 system/topoSetDict
topoSet
refineMesh -dict system/refineMeshDict -overwrite

decomposePar
mpirun -np 4 cavitatingFoam -parallel
reconstructPar

paraFoam

Three levels of coarseness are applied to the mesh using the utilities topoSet and refineMesh.

Calculation time

17 hours 9 minutes 58.9 seconds *4 parallel, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz