capillaryRise (blueCFD-core 2024 OF-12)をやってみた

blueCFD-core-2024の計算動作が今までとちがうようなので,ちょっと調べてみました.その時のメモです.


1. tutorialをコピー

 ・場所は "C:\blueCFD-Core-2024\OpenFOAM-12\tutorials\incompressibleVoF\capillaryRise" 

これを作業フォルダにコピーする 


2. Allrun の修正

・動作確認をしやすいように Allrunを調べ,加筆修正する.

 logがあるとAllrunが動作しないので,logを削除する. 

runApplication で計算開始ではなく,foamRun を実行するのは画面に計算の進捗がでるので個人的には好み. 


#!/bin/sh

cd ${0%/*} || exit 1 # Run from this directory 

 rm log* # logがあるとrunApplicaitonを実行しない 

 # Source tutorial run functions

 . $WM_PROJECT_DIR/bin/tools/RunFunctions 

 application=$(getApplication) 

runApplication blockMesh 

runApplication setFields 

#runApplication $application 

foamRun 


3. 試し計算

・計算手順はblockMesh ⏎, setFeilds ⏎, foamRun ⏎ だが

./Allrun で初期化,メッシュ作成,VOF値設定,計算実行が自動化されるので便利.


・ソルバーが incompressibleVOF に変わっているが,実行するのはfoamRun となる.

・OF11からアプリケーション・ソルバー → クラスとして記述されたモジュール・ソルバーを導入したらしい.

(参考)  https://wiki.opencae.or.jp/index.php/OpenFOAM_11_%E3%83%AA%E3%83%AA%E3%83%BC%E3%82%B9%E3%83%8E%E3%83%BC%E3%83%88 


・結果を表示させると frontAndBackがないとワーニングがでる. うるさいので 0\alpha.water.org, p_rgh, U の default 境界を全部 frontAndBack に置き換えるとよい ,


    defaultFaces { type empty; } →   frontAndBack { type empty; }


4. 試し計算2

・Allrun を実行すると runApplication setFields の際に alpha.water.orgをalpha.waterに自動的にコピーするので楽 だが・・・

→しかしparaview で0\alpha.water.org があると時刻0の表示がおかしくなる(VOF値を表示しない)ので,やはり0の外に出す方が個人的にはよいと思う. 

・Allrunの処理:  alpha.water.orgを0フォルダの外に出して,計算最初にコピーするように以下に変更する.


#!/bin/sh 

cd ${0%/*} || exit 1 # Run from this directory 

 # Source tutorial 

clean functions

. $WM_PROJECT_DIR/bin/tools/CleanFunctions 

cleanCase # 計算結果を削除 

rm log* # logがあるとrunApplicaitonを実行しない 

cp alpha.water.orig 0/alpha.water

 # Source tutorial run functions

. $WM_PROJECT_DIR/bin/tools/RunFunctions 

application=$(getApplication) 

runApplication blockMesh 

runApplication setFields # alpha.water.orgをalpha.waterに自動的にコピー 

#runApplication $application 

foamRun 


・待ち遠しいので,controlDict の結果書き出し間隔を細かくして,途中経過をこまめにparaview で確認する. 

writeInterval 0.01;

・0.5 まで計算するのに40分くらいかかった.




 




0コメント

  • 1000 / 1000