チュートリアルケース(motorBike LES のメッシュ作成について)
223 views
Skip to first unread message
N.K
Nov 15, 2022, 4:23:28 AM11/15/22
to OpenFOAM
状況説明
LESのmotorBikeのチュートリアルケースを流用して、固定翼を解析することを検討しています。できる限りチュートリアルケースと同じ状態でエラー箇所を絞るために
チュートリアルケースで与えられるmotorbike obj gz に合わせる形で、CADから固定翼のobjモデルを生成しgz形式に圧縮
モデル名もmotorbike obj gz として
３Ｄモデル以外は全てチュートリアルケースと同一の条件にしたのですが
mpirun -np (設定コア数) renumberMesh -overwrite -parallel において
次のエラーが出てしまいました。
参考書等を見る限り、BlockMesh SnappyHexMeshでの直方体の空間メッシュやlocation In Meshは問題がないように思われ且つ
チュートリアルケースの３Dデータでは問題なくrenumberMeshされたのですが
３Dモデルデータ以外をチュートリアルケースから変えていないにも関わらず、検索エラーが出るということは
SnappyHexMeshDictの設定の寄与等でPatchField を正確に出力できていないということでしょうか(featureAngle等)
初心者であるため頓珍漢なことを書いているかもしれませんが
アドバイスがあれば、よろしくお願いします。
以下を参考にコマンドは実行しています
チュートリアルケースの
N.K
Nov 15, 2022, 4:26:14 AM11/15/22
to OpenFOAM
追伸 ュートリアルケースの３Dデータでは問題なくrenumberMeshされたのですが＞ cp $FOAM_TUTORIALS/resources/geometry/motorBike.obj.gz constant/triSurface/ から得られるバイク型の３Dモデルを指しています
2022年11月15日火曜日 18:23:28 UTC+9 N.K:
kominami
Nov 15, 2022, 6:18:48 AM11/15/22
to OpenFOAM
NKさん
kominamiです。
1)
motorbike obj gzは、obj形式のの3D-CADデーターの圧縮書庫です。もともとのチュートリアルのモーターバークの形状ファイルの容量が大きいため圧縮書庫になっています。NKさんの翼形状は単純な形のため、ファイルの中身が圧縮形式でないのにも関わらずファイル名だけを変えてしまっていませんか？
2)
＞以下を参考にコマンドは実行しています
Allrunスクリプトでなくで、OpenFOAM端末に一行ずつコマンド入力をしているという理解で良いですか？
単純に、renumberMesh コマンドならどうなりますか？（パラレル実行よりもエラーメッセージで表示される情報が判りやすいと思います。不足しているpatch名が表示されるかもしれません。
renumberMeshの前にparaview(paraFoam)でメッシュは表示されますか？ この段階で正常に表示されないのならば、メッシュが正常に作成されていない可能性が高いと思います。
3)
エラーメッセージの前後、特に前の情報も必要です。
＝＝＝＝＝＝＝＝＝＝＝＝
・OpenFOAM Google Groupに投稿する際に抑えるポイント
2022年11月15日火曜日 18:26:14 UTC+9 N.K:
Message has been deleted
N.K
Nov 15, 2022, 7:54:34 AM11/15/22
to OpenFOAM
To: OpenFOAM
kominamiさん 回答ありがとうございます。説明不足が多く恐れ入ります。
1) ＞翼形状は、motorBike と名付けたobjファイルをCADで出力後にgzipコマンドでgzに変更しています。
圧縮した目的はチュートリアルケースとの違いを極力無くす事です。
(元々はstlファイルを読み込ませていたのですが、stlファイルで生じたエラーがobj gzでも再現するのかを確認するためです。）
2) OpenFOAM端末に一行ずつコマンド入力をしているという理解で良いですか？
＞はい、一行ずつコマンド入力をしてエラーが生じた処理を特定したいと考えています。
renumberMeshの前にparaview(paraFoam)でメッシュは表示されますか？
＞renumberMeshの直前にfoamコマンドで Case Type decomposed Case で読み込んだ場合に、先ほどの投稿でのメッシュが出力されています
単純に、renumberMesh コマンドならどうなりますか？
＞実行結果を添付させていただきます。
3) ＞ RenumberMesh 8並列のコードについて実行したログを txtデータに記入したものを添付させていただきます。
また 実行環境はOSがWindows10 (WSL)でUbuntu20.04.LTS OpenFoamは v2012となっています
2022年11月15日火曜日 20:18:48 UTC+9 kominami:
N.K
Nov 15, 2022, 8:19:46 AM11/15/22
to OpenFOAM
こちらがケースディレクトリ全体をアップロードしたGoogleドライブになります。
2022年11月15日火曜日 21:54:34 UTC+9 N.K:
haruka....@gmail.com
Nov 16, 2022, 10:11:54 AM11/16/22
to OpenFOAM
アップロードされたケースディレクトリの「motorBike/motorBike/constant/triSurface/motoroBike.obj.gz」をダウンロードして、gzip展開するとファイル名が「wing.obj」になるのでこのあたりが問題になってそうな感じがします。
OpenFOAMがmotorBike.obj.gz内のmotorBike.objを参照しようとしてwing.objしかないのでエラーが出てる可能性があるように思います。
$mv wing.obj motorBike.obj && gzip motorBike.obj
とかしてmotoroBike.obj.gzを作り直して試すといいのではないでしょうか？
2022年11月15日火曜日 22:19:46 UTC+9 N.K:
N.K
Nov 19, 2022, 11:12:15 AM11/19/22
to OpenFOAM
H-Tsubotaさん ありがとうございます
Windows上でobjファイル名を変更したのか、とファイル名が「wing.obj」になっていたみたいです。
しかしながらgz展開後motorBike.objになる状態で再試行したのですが
https://www.xsim.info/articles/OpenFOAM/tutorials/incompressible-pisoFoam-les-motorBike.html に記載されている
mpirun -np 8 renumberMesh -overwrite -parallel のコマンドにおいてエラーが再現されました。
decomposeParDictは領域分割数等は使用しているCPUに一致しておりSnapyHexMeshも実行できていることからdecomposeParDictには問題がないと思いますが、
チュートリアルケースでは
Reading volScalarField: k nuTilda nut p
Reading volVectorField: U
Reading volVectorField: U
と読み込んでいるため３Dモデルの出力に問題がありそうです。
手探りで試行錯誤を進めみます。
2022年11月17日木曜日 0:11:54 UTC+9 H-Ts...@xsim.info:
N.K
Nov 20, 2022, 8:41:58 PM11/20/22
to OpenFOAM
デフォルトのチュートリアルケースについてpatch名を書き直したところ
potentialFoamまでは実行できたため
何故最初のチュートリアルケースがうまくいったかは不明のままですが 一応の解決はしました。
皆さまコメントありがとうございました。
2022年11月20日日曜日 1:12:15 UTC+9 N.K:
Reply all
Reply to author
Forward
0 new messages