津波解析(OpenFOAM)

津波解析(OpenFOAM)
２０１５年６月２５日
１２Ｔ０２８
久保栞
1. OpenFOAM
前回まで
•
チュートリアル
•
メッシュの確認
•
圧力分布の表示
•
流線プロット
•
ベクトルプロット
•
ケース形状の変更
•
メッシュ生成(穴あき板の応力解析)
•
ダムの決壊（単純なモデル）
現時点での進捗状況
• チュートリアル
• 計算時間、時間ステップの変更
• さまざまなモデルでの解析（2D，3D）
• グラフの書き方
今後の課題
• グラフの書き方
• さまざまなモデルでの解析(2D，3D)
1.1 メッシュ生成(穴あき板の応力解析)
•
pointsがないと解析できない？(前回)
•
→ solidDisplacementFoamでは不要
圧力分布(σxx)
1.2 ダムの決壊（単純なモデル）
•
右図の状態で解析
•
T=0~3sまで表示（動画）
•
右端に出っ張りを追加するとどうなるか
T=0~3s
deltaT=0.001s
α1=1 ・・・液相
α1=0 ・・・気相
(α1：相比率)
界面において相比率は
0から1の間の任意の値を
とるためくっきりと定義されない
0.7
0.7
0.6
0.6
0.5
0.5
0.4
0.4
0.3
0.3
0.2
0.2
0.1
0.1
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0
0
0.1
0.2
0.3
0.4
エクセルで新たにpointsを作成(今回は右側を採用)
0.5
0.6
0.7
blockMeshDictの変更
•
vertices：ブロックの頂点の座標(9つのブロック)
•
blocks：頂点名とセル番号から定義される
•
hex (0 1 5 4 18 19 23 22) (23 8 1) simpleGrading (1 1 1)
•
形状
•
•
•
頂点名
各方向に必要なセル数
各方向のセル拡大比
boundary：境界の形状とその座標
leftWall { type wall; faces ( (0 18 22 4) (4 22 26 8) (8 26 31 13) ); }
壁の位置
壁の種類
頂点名
液相と気相の分布(T=0s)
wireframe
時間を進めようとすると以下のエラーが発生する
ERROR: In C:¥DBD¥pvs-x64¥paraview¥src¥paraview
¥VTK¥IO¥Geometry¥vtkOpenFOAMReader.cxx, line 6573
vtkOpenFOAMReaderPrivate (00000000107EEA30):
Number of cells/points in mesh and field don't match: mesh = 1862, field = 2268
mapFieldsユーティリティを用いて再度解析し直す？
(mapFields：形状や境界のタイプもしくはその両者が不一致な場を位置づけることができる)
1.3 その他
T型分岐管の流れ
3Dタンクモデル
越流堰
www.youtube.com/watch?v=UiuX7IwJ3tA
https://sites.google.com/site/freshtamanegi/home/openfoam/tutorial
1.4 T型分岐菅の流れ
•
T=0~3s
•
deltaT=0.001

Download Report