流路絞り部での蒸気発生(LES)

OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/multiphase/cavitatingFoam/les/throttle

概要

領域 inlet に全圧300×105、領域 outlet に静圧100×105を設定することで、圧力差によって inlet から outlet へ向かって駆動される水の流れを0.0002秒分だけ解析します。壁面 walls はすべり無し条件とし、Z 方向は1メッシュとして2次元問題としています。

モデル モデル

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

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

メッシュ(全体) メッシュ(全体)
メッシュ(絞り部) メッシュ(絞り部)

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

時刻0.00012秒での流速(U) 時刻0.00012秒での流速(U)
時刻0.00012秒での絞り部付近の蒸気の体積率(alpha.vapour) 時刻0.00012秒での絞り部付近の蒸気の体積率(alpha.vapour)
時刻0.00012秒での絞り部付近の流速(U) 時刻0.00012秒での絞り部付近の流速(U)

実行コマンド

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

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

cavitatingFoam

paraFoam

ユーティリティー topoSet、refineMesh を使用してメッシュに3段階の粗密をつけています。計算環境によってはソルバーが途中で停止する場合があります。その場合はコマンド「cavitatingFoam」を再度実行することで計算が再開されます。

計算時間

22分56.57秒 ※シングル、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz