[snappyHexMesh] featureAngle > 90の際に生じるエッジの欠けについて

M.TANAKA

unread,

Jul 4, 2022, 6:54:37 AM7/4/22

to OpenFOAM

お世話になります、田中と申します。

snappyHexMeshにおける境界層メッシュの設定について躓いている問題があり、投稿いたします。

以下のチュートリアルを例にとりました。

https://develop.openfoam.com/Development/openfoam/-/tree/master/tutorials/incompressible/pimpleFoam/RAS/rotatingFanInRoom

このケースに対し以下の変更を行い、境界層メッシュの生成を試みました。参考のため、編集したケースフォルダを添付します。

・addLayersControls.layersにパッチを指定

・addLayersControls.meshShrinkerをdisplacementMotionSolverに変更

まず、addLayersControls.featureAngleを60としたところ（画像1枚目）、画像右下の90度エッジや机の270度エッジに見られるように境界層メッシュが消えてしまう箇所が発生しました。

エッジにも境界層メッシュを生成するため、featureAngleを120としたところ（2枚目）、境界層メッシュは期待通りの形に作れたのですが、画像右下に見えるようなセルの欠けがエッジ上で複数生じました。

「エッジの欠けを防ぎつつ、エッジ部にも境界層メッシュを生成する」には、どのような設定が有効でしょうか？

なにかご存知の方がいらっしゃいましたら、教えていただけると幸いです。

参考にしたページ：

https://www.openfoam.com/documentation/guides/latest/doc/guide-meshing-snappyhexmesh-layers.html

http://www.wolfdynamics.com/wiki/meshing_OF_SHM.pdf

↑ "featureAngle 60;"

↑ "featureAngle 120;"

rotatingFanInRoom.tar.gz

E.Mogura

unread,

Jul 9, 2022, 9:53:03 PM7/9/22

to open...@googlegroups.com

E.Mogura です。

質問の答えにはなっていませんが、コメントさせていただきます。

質問の趣旨は、「SHMで境界層レイヤーをきちんと作成する方法はいかに？」ということになろうかと思いますが、よほどシンプルな形状でない限り、経験上難しいというのが回答です。

どうしても（フリーのソフトで）境界層レイヤーをきちんと作成したい、というのであれば、cfMeshで作成することをお薦めします。以下に本例での作成サンプルを提示しておきます。

これは、添付のFreeCADモデルで、拙作のDEXCSランチャー

https://gitlab.com/dexcsof/dexcs-launcher

を使えばごくごく簡単に作成できてしまいますのでお試しあれ。

但し、cfMesh の欠点としては、

・マルチリージョン（本例の場合AMI境界の作成）に対応できない。

・平面上に区分した境界（本例の場合、ドアやoutlet）がいびつになってしまう

といった点があります。いずれも対処方法はありますが、それなりに手間が加わるという点には留意下さい。

2022年7月4日(月) 19:54 M.TANAKA <tanaka19m...@gmail.com>:

--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
このディスカッションをウェブ上で閲覧するには https://groups.google.com/d/msgid/openfoam/4e9c5cf0-f294-46d6-8b6a-999c12290b7an%40googlegroups.com にアクセスしてください。

model.FCStd

M.TANAKA

unread,

Jul 13, 2022, 8:48:58 AM7/13/22

to OpenFOAM

E.Moguraさん

コメントありがとうございます。やはりsHMで商用メッシャーのように表面全体に境界層メッシュを作成するのは困難なのですね。

OSSの範囲内で解析したいと思ってはいますが、形状が複雑な場合、最終的にはPointwiseなどの商用メッシャを使わざるを得ないのかもしれませんね。

ひとまず、教えていただいたcfMeshを使って、２つのメッシュを別々に作ってマージし、マルチリージョン化するという手法にトライしてみます。何かしら形になれば報告したいと思います。

2022年7月10日日曜日 10:53:03 UTC+9 E.Mogura:

M.TANAKA

unread,

Nov 20, 2022, 12:07:10 AM11/20/22

to OpenFOAM

前回の投稿から期間が空いてしまいましたが、進展がありましたので報告します。

私は普段ESI版OpenFOAM(v2206)を使っているのですが、Foundation版(v10)と比較してみたところ、v10ではsnappyHexMeshのfeatureAngleが60 degでもエッジ周りの境界層が正しく生成されることを確認しました。

下記リポジトリにまとめています。

https://gitlab.com/openfoamscripts/investigation/compare-boundary-layer-mesh-quality-between-v10-and-v2206

複雑形状にも境界層メッシュが正しく貼られるのであれば、Foundation版に乗り換えるべきかもしれない、と個人的に考えています。