1 / 95

BFC グリッドバージョンの試計算 MEC-model boundary-fitted coordinate grid version

BFC グリッドバージョンの試計算 MEC-model boundary-fitted coordinate grid version. 大阪大学大学院工学研究科 地球総合工学専攻　博士前期課程２年　 栗岡 俊介. 2001.12.07( 金 ) 　第 2 回 MEC ワークショップ. ◆ 発表の手順. １．モデル開発の背景 ２．モデルの概要 ３．支配方程式 ４．数値計算法 ５．モデルの試計算 ６．今後の課題. モデル開発の背景.

vidal
Download Presentation

BFC グリッドバージョンの試計算 MEC-model boundary-fitted coordinate grid version

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. BFCグリッドバージョンの試計算MEC-model boundary-fitted coordinate grid version 大阪大学大学院工学研究科 地球総合工学専攻　博士前期課程２年　 栗岡 俊介 2001.12.07(金)　第2回MECワークショップ

  2. ◆発表の手順 • １．モデル開発の背景 • ２．モデルの概要 • ３．支配方程式 • ４．数値計算法 • ５．モデルの試計算 • ６．今後の課題

  3. モデル開発の背景

  4. これらは海岸線および海底などの境界付近は階段状で近似されており、 3次元解析モデルを境界に埋め込むような場合に静水圧近似モデルの境界付近の計算を精度よく扱う必要がある。 背景 • MECモデルの数値計算プログラム • MEC Variable Mesh バージョン • 直行不等間隔格子を用いたモデル • MEC Nesting Mesh バージョン • 　スケールの異なる直行等間隔格子を用いたモデル

  5. 従来のMECモデルによる格子生成例① 大阪湾 五ヶ所湾 陸地部分（緑色）の格子は計算されない 海岸線が階段状近似 ※第1回MECワークショップCD-ROMより作成

  6. 従来のMECモデルによる格子生成例② 長崎県大村湾 拡大 ※MEC九大グループHPより引用

  7. 目的 本年度はＭＥＣモデルの静水圧近似モデルの部分の計算で使用される計算格子に境界適合格子を用いたバージョンを新たに開発し、試計算を行なうこととした。

  8. モデルの概要

  9. モデルの概要 • MEC Variable Mesh バージョンの静水圧近似　　部分のアルゴリズムを採用した。 • 計算格子に境界適合格子を用いて湾内の海水　部分を構造格子で表現した。 • 潮位の変化に伴って動く移動格子を海水面に　適用した。

  10. モデルの概要 • 計算格子セルを有限体積法によりコントロール　ボリュームとして離散化することで複雑な地形　　を計算格子で表現した。閉境界での境界条件　　の処理が容易かつ高精度に行なえる。 • 計算格子を分割し接合するマルチブロックの方　法を用いることでより複雑な地形表現が可能と　なった。

  11. 支配方程式

  12. ここで、　　　は水平方向流速、　　は鉛直方向流速、　　は圧力、　　は密度、ここで、　　　は水平方向流速、　　は鉛直方向流速、　　は圧力、　　は密度、 　　は水平渦動粘性係数、　　は鉛直動渦粘性係数である。 支配方程式 • Navier-Stokes方程式 • 境界条件 • 連続の式 （海底） （海面）

  13. 計算格子

  14. 計算格子 η y ξ x • 境界適合座標では物理空間で複雑な形状を扱うためにその形状に沿った境界適合格子を用いる。 • 物理空間はデカルト座標系で表され、境界適合格子が正方格子に変換されるような計算空間を考える。 座標変換 (b)計算空間 (a)物理空間 座標変換（２次元の例）

  15. 支配方程式の座標変換

  16. Navier-Stokes方程式 • 水平方向流速u • 水平方向流速v

  17. Navier-Stokes方程式 • Navier-Stokes方程式中の拡散項における応力ついて　は以下のような拡散流束を計算した。

  18. 連続の式 • 鉛直流速w、水位ztについてはMEC Variable Meshバージョン　と同様に以下の連続の式を積分することで求めた。 ※ここで、Ｊはヤコビアン(Jacobian)であり、計算格子セルの体積をあらわす。

  19. 支配方程式の離散化

  20. ●　 格子点(x,y,z) ×　 流場変数(u,v,w,p) 　 　コントロールボリューム k+1 j+1 ζ k η j i-1 ξ i k-1 i+1 j-1 変数配置 格子セル（i,j,k）のコントロールボリューム ’J’

  21. ●　 格子点(x,y,z) ×　 流場変数(u,v,w,p) 　 　コントロールボリューム k+1 j+1 N H k j E E W W S L J-1 k-1 i-1 i i+1 i-1 i i+1 Y-Ｚ 平面 X-Y 平面 コントロールボリューム（移流項）

  22. ●　 格子点(x,y,z) ×　 流場変数(u,v,w,p) 　 　コントロールボリューム k+1 j+1 N H k j E E W W S L J-1 k-1 i-1 i i+1 i-1 i i+1 コントロールボリューム（拡散流束を求める時） Y-Z平面 X-Y 平面

  23. 反変ベクトルの導出 の導出（格子セルのEASTsideの例） ζ η ξ 変数 ( u,v,w) NORTH WEST (i,j,k) EAST SOUTH ξが一定の格子断面の面積ベクトルで面積の大きさを持ち、方向は法線方向である。

  24. 格子セル体積Jの導出 変数 ( u,v,w) NORTH WEST (i,j,k) EAST SOUTH ζ η 格子セル体積Jはセルセンターを通るξ,η,ζ方向のベクトル から求めた。 ξ

  25. 格子セル(i,j,k)で離散化 コード内の支配方程式の離散化（ Navier-Stokes方程式・水平方向流速u） • 実際のコード内でどのような離散化を行なっているかを移流項を例に示す。 • 座標変換後の移流項

  26. 運動量フラックス 変数 ( u,v,w) NORTH (i,j,k) EAST WEST SOUTH 任意のｋにおけるi-j 断面の図 • 格子セル(i,j,k)で離散化

  27. 格子セル(i,j,k)で離散化 鉛直流速ｗの導出 • 連続の式

  28. 流量フラックス 変数 ( u,v,w) 既知 既知 既知 既知 任意のｊにおけるi-k 断面の図 HIGH ② (i,j,2) EAST WEST ① (i,j,1) EAST WEST LOW 既知 →海底からの流入は“0” 連続の式において、水平流速u,vを既知として海底(k=1)から海面まで積分していくことで求まる。

  29. 自由水面の境界条件 • 自由水面の境界条件 移動格子の水面境界を抜ける流れはないことから、自由水面においては次式が成り立つ。 右辺は連続の式から求まる。これはMEC Variable Mesh バージョン と同様、一般的な自由表面の境界条件をあらわしている。

  30. 数値計算法 数値計算法の概要

  31. モデルの試計算

  32. モデルの試計算 • 局部的に水路が狭くなる湾を想定して計算を行なった。(Case01) • きわめて狭い水路をもつ湾を想定し、水路部分を別ブロックに分割して計算を行った。(Case02) • 湾内に構造物があることを想定して計算を行なった。(Case03) • 細長い湾を想定し小さな湾を別ブロックで追加配置した計算を行った。 (Case04)

  33. Case01 計算条件

  34. 計算格子（全体） 流入境界

  35. 計算格子 流入境界

  36. 計算結果（潮位コンター図） 1/4周期後 基準時(5周期目) 3/4周期後 2/4周期後

  37. 計算結果（潮位アニメーション）

  38. 計算結果（水位 zt時系列） 流入境界 湾奥

  39. 計算結果（速度ベクトル・上潮時拡大）

  40. 計算結果（速度ベクトル・下潮時拡大）

  41. 計算結果（速度ベクトル・アニメーション）

  42. 計算結果（速度ベクトル・アニメーション）

  43. Case01の考察 • 局部的に水路が狭くなる湾を想定して計算を行なった。 • 奥の湾は流入境界の振幅0.5mに対して0.2mとなり、約90°の位相遅れが確認できた。これは細い接合部における湾奥の挙動に近いものであり、今後複雑なモデルでの計算を進めていく必要がある。

  44. Case02 計算条件

  45. 計算格子（全体） 流入境界

  46. 計算格子 流入境界 接合部拡大図

  47. 計算結果（潮位コンター図） 基準時(5周期目) 1/4周期後 3/4周期後 2/4周期後

  48. 計算結果（潮位アニメーション）

  49. 計算結果（速度ベクトル・上潮時拡大）

  50. 計算結果（速度ベクトル・下潮時拡大）

More Related