iric Full 3D Simulation Engine NaysCUBE & Nays 北海道大学 木村一郎 1
Agenda Part 1: Nays CUBEの基本コンセプト Part 2: Nays CUBEの主な特徴 Part 3: Nays CUBE 計算事例 Part 4: Nays CLIP ( 鉛直二次元モデル ) Part 5: Nays CUBEの基本操作 ( 時間があれば簡単なデモを行います.) + 2
Part I Nays CUBE の基本コンセプト 3
基本コンセプト 水平方向を一般曲線座標鉛直方向をシグマ座標実河川の計算のみならず, 開水路実験の再現も対象とする 茨戸川を三次元格子で分割した例 4
NaysCUBE の名前の由来 物理空間 計算空間 一般曲線座標の物理空間を一辺が 1 の立方体 (CUBE) の集まりからなる, デカルト座標の計算空間に変換して, 計算を実行している. 5
水面変動の計算方法の分類 境界捕捉法 境界適合法 VOF, MAC method, Density Function Method, etc. 柔軟性が高いが計算負荷大河床変動への適用に課題 計算効率が良く, 実河川の計算に適する. 河床変動への適用が容易 6
障害物や植生の取り扱いが容易 水面振動は運動学的条件で計算河床変動は流砂モデルと河床連続式で計算水没 非水没構造物の考慮水没 非水没植生の考慮 7
障害物設定を応用した例 複断面開水路 堰 8
Part II Nays CUBE の主な特徴 9
多彩な境界条件設定 境界条件設定ダイアログ 主流方向 : 通常境界条件, 周期境界条件, 反転周期境界条件横断方向 : 通常境界条件, 周期境界条件, 対称境界条件構造物壁面 : 摩擦考慮, 非考慮 10
多彩な境界条件設定 (1) 1 主流方向の周期境界条件が有効な例 周期境界条件 : 流出部の水理量をすべて流入部に与える周期性を有する現象を一周期分の計算領域で効率よく再現できる. 連続蛇行水路 砂堆 水制群 交互砂州 11
多彩な境界条件設定 (2) 2 主流方向の反転周期境界条件が有効な例 反転周期境界条件 : 流出部の水理量を左右反転して上流に与える半波長で反転する周期性を有する現象を半周期分の計算領域で効率よく再現できる. 連続蛇行流路 交互砂州 魚道など 12
多彩な境界条件設定 (3) 3 横断方向の周期境界条件が有効な例 周期境界条件 : 流出部の水理量をすべて流入部に与える周期性を有する現象を一周期分の計算領域で効率よく再現できる. 横断方向連続構造物 三次元砂堆 13
多彩な境界条件設定 (4) 4 横断方向の対称境界条件が有効な例 対称境界条件 : 左右が対称な流れ場の対称軸を境界とし, 鏡面反射させる境界条件で, 計算を効率化左右岸の一方, あるいは両方に設定できる. FLOW 左岸対称境界条件 FLOW 半分の領域だけで計算 14
多彩な境界条件設定 (5) 植生密生域の平衡区間を反転周期境界条件と対称境界条件の組み合わせで表現した例 15
計算スキーム 乱流モデルなど 移流項スキーム ( 計算精度を大きく左右する ) 一次風上スキーム ( 安定, 低精度 ) 三次 TVD-Muscl スキーム ( 高精度 ) 乱流モデル ( 二次流や非定常現象, 剥離などを再現する上で重要 ) 標準型線形 k-ε モデル ( 安定, 低精度 ) 二次非線形 k-ε モデル ( 高精度, 時として不安定 ) 時間積分 : 二次アダムス バシュフォース法 圧力計算 : 動水圧を考慮した HSMAC 法 16
河床変動の計算方法 現バージョンは掃流砂のみに対応 ( 後述の鉛直二次元ソルバー Nays CLIP では浮遊砂にも対応 ) 掃流砂量モデル ( 通常は, 主流方向の河床勾配の影響を考慮した Kovacs-Parker モデルを推奨 ) MPM (Mayer Peter Muller) モデル 芦田 道上モデル Kovacs-Parker モデル 河床変動に関する計算条件設定ダイアログ 17
Part III Nays CUBE 計算事例 18
事例 1: 某河川の湾曲部の流れ 流れのみの計算第一種二次流の再現 19
事例 2: 連続蛇行部の河床変動 (1) 初期河床 : 平坦周期境界条件 : 無限に連続する蛇行河川における河床変動河岸形状 : サインジェネレイティドカーブ 20
事例 2: 連続蛇行部の河床変動 (2) 21
事例 2: 連続蛇行部の河床変動 (3) 流れの剥離が再現されている 河床変動と水面変動 t=300sec の形状と流線 22
事例 3: 植生密生域の流れと河床変動 FLOW 簡単のため角柱でやってみる. 23
事例 3: 植生密生域の境界条件設定 24
事例 3: 植生密生域の流れと河床変動 ( 動画 ) 25
事例 3: 植生密生域 (t=13sec) 26
事例 4: 第二種二次流 ( 縦筋 ) と河床変動 直線水路の二次流再現 直線人工河川 ( 琵琶湖疏水 ) で観察された第二種二次流 Nezu and Tominaga (1993) 河床に流方向に筋状の凹凸 27
事例 4: 第二種二次流と河床変動 周期境界条件 一様な矩形断面直線開水路を仮定 主流方向の格子洲は,3 セル ( 横断方向 50, 鉛直方向 10) 第二種二次流は乱流の非線形性に起因する現象 乱流モデルは必ず 非線形 k-ε モデル を選ぶ!! 28
事例 4: 第二種二次流と河床変動 周期境界条件 一様な矩形断面直線開水路を仮定 主流方向の格子洲は,3 セル ( 横断方向 50, 鉛直方向 10) 第二種二次流は乱流の非線形性に起因する現象 乱流モデルは必ず 非線形 k-ε モデル を選ぶ!! 29
事例 4: 第二種二次流 ( 乱流モデルの設定 ) 30
事例 4: 第二種二次流 ( 二次流発達の動画 ) 31
事例 4: 第二種二次流 ( 河床変動の動画 ) 32
事例 4: 第二種二次流 (t=400sec) 33
Part IV Nays CLIP ( 鉛直二次元モデル ) 34
NaysCUBE Nays CLIP とは 完全三次元モデル掃流砂のみ計算時間大 三次元空間から鉛直二次元空間を 切り取る というイメージから CLIP と名付けた. Nays Nays CLIP の特徴 鉛直二次元モデル掃流砂 + 浮遊砂計算時間小 Nays CLIP は Nays CUBE から派生したエンジン 鉛直方向の変化が支配的な現象を軽快に再現する. 現 Nays CUBE に含まれていない浮遊砂モデルを装備 ( 但し,Nays CUBE 次期バージョンは浮遊砂モデル装備予定 ) 35
Nays CLIP の概要 適用例 : ダムなどの深い流の流動や土砂移動 Dune のような二次元的河床波 堰を超える流れなど スキーム, 乱流モデルなど NaysCUBE と同一 (Nays CUBE のマニュアルがほぼそのまま適用できる.) 計算格子 等幅で, 横断方向の格子数が 1 の格子を用いる. Nays CLIP で用いる格子の例 36
事例 : ダムの浮遊砂堆砂とゲート排出 循環流の発生 浮遊砂の一部はゲートから排出 堆砂の進行 クレストゲート越流時の堤体底部付近の循環流 この渦と浮遊砂がいかなる相互作用? 37
事例 : ダム堆砂 ( 浮遊砂条件設定 ) 38
事例 : ダムの浮遊砂堆砂とゲート排出 計算格子 浮遊砂モデル 板倉 岸モデル 39
事例 : ダム浮遊砂堆砂ゲート排出 ( 動画 ) 40
Part V Nays CUBE の基本操作 ( 時間があればここで簡単なデモを行います.) 41
計算のフロー Step 1: 計算格子の作成 Step 2: 計算条件の設定 Step 3: 計算の実行 Step 4: 計算結果の可視化 実河川の地形を用いる場合 単純な形状を iric 内で作成することも可能 流量, 計算時間, 河床条件などを設定 計算エンジンが計算を実行します. ユーザは実行ボタンを押すだけです. 計算途中でも, 次の可視化を行うことができます. ベクトル, 流線, 等値面, パーティクル, コンターなどの多彩な描画機能が用意されています. 作業終了 42
おわりに このスライドで示した計算例の計算所要時間は, 事例 1 を除いて,10 分程度以内です.( ノート PC 使用 ) スパコンがなければ三次元計算ができなかったのは完全に過去の話です. ぜひ, この Nays CUBE (& Nays CLIP) を通じて, 三次元計算 ( あるいは鉛直二次元計算 ) にトライしていただけば幸いです. + 43