前回のつづきです。
流体解析の流れ
流体解析は以下の流れで行う。
-
Pre-processing
解析に用いる3Dデータを編集したり、空間をセルで分割し、メッシュファイルを作成したりする。 -
Solving
物理モデルや境界条件などの設定をして、流体シミュレーションを実行する。 -
Post-processing
シミュレーション結果を可視化する。
今回は、Pre-procceingについて述べる。Pre-processingは流体解析の中で特に根気と時間を要する段階だ。
DiscoveryでOBJファイルを編集する
対象とする区画を含むOBJファイルを探す
01100_indexmap_op.pdfを確認し、解析したい区画が含まれる標準地域メッシュ を特定する。bldg/lod1から目的の標準地域メッシュのOBJファイル (例：64414278_bldg_6680.obj) を探す。
lod1よりもlod2のほうがより詳細に建物の形状が表現されている。その代わり、データサイズは大きい。Ansys Studentではデータサイズの大きいオブジェクトファイルはメッシュ化の際に不都合なので、lod1を勧める。
DiscoveryにOBJファイルをインポート
Discoveryを起動する。「新規作成」を選択する。OBJファイル (例：64414278_bldg_6680.obj) をDiscoveryにドラッグアンドドロップすることでインポートできる。
y軸とz軸が入れ替わって表示される場合は、「ファセット＞リバースエンジニアリング＞ファセットの向き」で軸の方向を適切に設定する。
OBJファイルをインポートしたら、まずはファセットをチェックする。「ファセット＞クリーンナップ＞ファセットをチェックする」をクリックし、モデルツリーの「ファセット」をクリックする。(左上にある階層構造の図をDiscovery上ではモデルツリーとよぶ)
「メッシュは自己交差しています」などのエラーメッセージが表示された場合は、ファセットを自動修正する。「ファセット＞クリーンナップ＞自動修復」をクリック。このときも先ほどと同様にモデルツリーの「ファセット」をクリックして選択する。再度ファセットをチェックして問題がないことを確認する。エラーが0になっていればOK。（「メッシュには複数のピースがあります」や「メッシュは折り返し部分を含む」などのメッセージは無視しても問題ない。）
次に、ファセットのままでは3Dデータの編集に向かないので、ファセットをソリッドに変換する。モデルツリーのファセットを右クリックし、「ソリッドに変換＞面をマージ」を選択する。データサイズが大きいと、それなりに時間がかかる。
👆変換前
👆変換後
DiscoveryでOBJファイルを編集
建物の編集
解析に不要な建物は削除する。不要な建物をモデルツリーから選択するか、ドラッグやクリックによって直接建物を選択して削除する。
初めてANSYS Fluent (Student版) で流体解析を行う場合は、建物を2-3個だけ残して他は削除することを勧める。なぜなら、建物の数が多いほどメッシュ数が増大し、Student版のメッシュ数の制限を越えてしまうため。まずは、どれくらいの建物の量でメッシュ数がどれだけになるのかという感覚を掴むとよいと思う。
必要な建物だけ残したら、適宜構造を簡略化していく。
以下の図のような小さな突起はシミュレーションの精度に影響を与えるため、「デザイン＞編集＞プル」の機能を使って突起を除去する。地道な作業だが、一つずつ手作業でやっていく。
👆不要な建物を削除し、建物の構造を適宜簡略化した後の様子
計算領域の作成
次に、建物を囲うボックスを作成する。このボックス内が流体解析を行う計算領域となる。
都市の風環境予測のためのCFDガイドブックによると、計算領域の条件として、
- 側面境界、上空境界の大きさは建物の外縁から5H (H：代表的な建物の高さ) 程離れていること
- 閉塞率 (風が吹き込んでくる方向から見て、建物の面積が占める面積の割合) が5%以下であること
が推奨されている。
👆閉塞率（赤線で囲ったところ）
デザイン＞スケッチ の機能を使ってボックスを作成する。
👆風を吹かせる面 (inlet) と風が通り抜ける面 (outlet) は作らない
プリペア＞生成＞ボリュームの抽出 で計算領域を作成する。inletとoutletの枠をクリックし、「シード面」というボタンをクリックしたら、今度は内側の面 (地面など) をクリックする。「内側の面をプレビュー」し、内側の面が赤く塗りつぶされていることを確認する。問題なければチェックボタン
をクリックする。すると、ボックス内を満たす「Volume」というパーツが新たに生成される。
最後に、計算領域に含まれる境界の名前をつける。画面右下あたりの「詳細選択」をクリックする。Volumeのinletの面をクリックした状態で、詳細選択の「＋選択」をクリックすると「Group1」が生成される。これを「inlet」に名前を変更する。同様にして、outlet、top (上空面)、lateral_boundary (側面x2)、ground (地面)、bldg (建物) の境界を作成する。
👈inlet 
👈outlet
👈top 
👈lateral_boundary 
👈ground 
👈bldg シミュレーションに含めるパーツの選択
Volume以外のパーツはシミュレーションに含めないため、左上に表示されているモデルツリーで、アイコン (+/-のようなもの) をクリックして除外する。
👆Volume以外をシミュレーションから除外した様子
ここまで完了したら、Discoveryファイル (拡張子.dsco) として保存する。
メッシュファイルを作成する
次に、Discoveryで作成したファイルからメッシュファイルを作成する。
Fluentを起動する。
Fluentのアイコンをダブルクリックすると、Fluent Launcherが表示される。「Meshing (メッシュモード)」を選択し、Startボタンをクリックし、起動する。
メッシュ化の設定
Workflow TypeはWatertight Geometryを選択する。
表示されているWorkflowに従って進める。
Import Geometry
基本的にデフォルトのままで設定を変更する必要はない。File Nameの横にある「...」をクリックして、Discoveryで作成したファイルを選択する。
Add Local Sizing
以下のように設定する。
[Would you like to add local sizing?] --> [yes]
[Growth Rate] --> [1.2]
[Size Control Type] --> [Face Size] *適宜適切なものを選択
[Target Mesh Size] --> *何通りかメッシュファイルを作成して、適切なサイズを見つける。初めはとりあえずデフォルト値のままにする
解析したい区画周辺を詳細にメッシュで区切りたいので、bldgとgroundを選択して、local meshを作成する。
Generate the Surface Mesh
[Minimum Size] --> *何通りかメッシュファイルを作成して、適切なサイズを見つける。初めはとりあえずデフォルト値のままにする
[Maximum Size] --> *何通りかメッシュファイルを作成して、適切なサイズを見つける。初めはとりあえずデフォルト値のままにする
以降はデフォルト値のまま
Describe Geometry
[Geometry Type] --> [The geometry consists of both fluid and solid regions and/or voids]
以降はデフォルト値のまま
Update Boundaries
以下のように設定する。
| Boundary Name | Boundary Type |
|---|---|
| inlet | velocity-inlet |
| outlet | pressure-outlet |
| top | symmetry |
| lateral_boundary | symmetry |
| ground | wall |
| bldg | wall |
Create Regions
デフォルトのまま
Update Regions
デフォルトのまま
volume-volumeのRegion Typeがfluidになっているか確認する
Add Boundary Layers
デフォルトのまま
Generate the Volume Mesh
[Solver] --> [Fluent]
[Fill With] --> [tetrahedral] *目的のセルの形を選択する。都市の流体解析を行った文献では多くがtetrahydralを採用しているため、ここではそれに従った。
生成されたメッシュファイルを確認
Consoleにメッシュ化の結果が出力される。
まず、セルの数を確認する。
Ansys Studentでは、セル数が1,048,576を超えると、メッシュファイルが保存できないため、もしこの上限値を超えていたらセルサイズを大きくするなどによって、セル数を減らすことを試みる。
次に、minimum Orthogonal Qualityを確認する。この値は0.1以上が望ましい。
minimum Orthogonal Qualityが0.1を下回っているセルはMeshウィンドウで赤く塗りつぶして表示される。
解析したい区画 (ここでは、青線で囲った1ブロック) のシミュレーション精度に影響を与えそうであれば、minimum Orthogonal Qualityが0.1を下回っている部分をDiscoveryで修正する。
セル数やminimum Orthogonal Qualityの条件をクリアしたら、File＞Write＞Mesh でメッシュファイルを保存する。
参考文献
- 都市の風環境予測のためのCFDガイドブック
- Geometry Prep for CFD Simulation Using Ansys Discovery：Discovery上での操作について解説したAnsys公式動画
- Mesh Generation using Ansys Fluent Watertight Workflow：メッシュ作成について解説したAnsys公式動画