流路絞り部での蒸気発生（3次元、LES）

OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/multiphase/cavitatingFoam/les/throttle3D

概要

領域 inlet に全圧300×105、領域 outlet に静圧100×105を設定することで、圧力差によって inlet から outlet へ向かって駆動される水の流れを0.0001秒分だけ解析します。壁面 walls はすべり無し条件とします。

モデル モデル

乱流計算には LES モデルを使用します。

メッシュは以下の通りで、メッシュ数は1659700です。

メッシュ（全体） メッシュ（全体）
メッシュ（絞り部） メッシュ（絞り部）

計算結果は以下の通りです。

時刻0.0001秒での流速（U） 時刻0.0001秒での流速（U）
時刻0.00005秒での絞り部付近の蒸気の体積率（alpha.vapour） 時刻0.00005秒での絞り部付近の蒸気の体積率（alpha.vapour）
時刻0.00005秒での絞り部付近の流速（U） 時刻0.00005秒での絞り部付近の流速（U）

実行コマンド

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

ユーティリティー topoSet、refineMesh を使用してメッシュに3段階の粗密をつけています。

計算時間

17時間09分58.9秒　※4並列、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz