[計算例 2] 単純な湾曲水路の浮遊物追跡
Figure 49 に示すような直線+湾曲+直線水路実験の流量解析およびトレーサーの追跡計算を行う。 この実験は水路全幅約1.8mで、左岸側半分が高水敷、右岸側半分が低水路の複断面水路で、右岸側のみが 移動床となっている。平面形状および断面形状の概略を Figure 50 に示す。 この実験は 寒地土木研究所 の委託により。 (株)建設技術研究所 によって行われたものである。
Figure 49 : 実験状況の動画
Figure 50 : 実験水路の形状
以下の計算事例では以下の手順で計算を実施する。
Nays2DHにより河床変動計算を実施し、河床変動がほぼ定常に達すた状態の河床形状を得る。
Nays2d+により準3次元流れ場の計算をする。
GELATOによりトレーサーの追跡を行う。乱流拡散強度パラーメーターを変更してその影響を比較する。
Nays2DHによる流れと河床変動の計算
ソルバの選択
iRICの起動画面から、[新しいプロジェクト]を選ぶと表示されるソルバの選択画面で、 [Nays2DH]を選んで[OK]ボタン押すと、
Figure 51 : ソルバーの選択
「無題- iRIC 4.x.x.xxxx [Nays2DH iRIC3X 1.0 64bit]」と書かれた Windowが現れる。
Figure 52 : 無題
計算格子の作成
[格子]->[格子生成アルゴリズムの選択]で現れるウィンドウで、 [2次元単純円弧格子作成ツール(複断面対応版)]を選んで[OK]を押す。 (Figure 53 )
Figure 53 : 格子生成アルゴリズムの選択
[格子生成]のグループ[水路形状]、[断面形状]、[追加水路]、[粗度と河床状態]を それぞれ、 Figure 54 、 Figure 55 , Figure 56 , Figure 57 のようにパラメータを設定して、最後に[格子生成]を押す。
Figure 54 : 格子生成(1)
Figure 55 : 格子生成(2)
Figure 56 : 格子生成(3)
Figure 57 : 格子生成(4)
「マッピングを実行しますか？」と聞かれるので[はい(Y)]を押す。 ( Figure 58 )
Figure 58 : 確認
オブジェクトブラウザーの[格子][セルの属性][固定床と移動床]に☑マークを入れると、 Figure 59 のように固定床部分が赤、移動床部分が青で示された 格子の図が示される。
Figure 59 : 固定床と移動床
低水路と高水敷の境界の固定床は護岸を想定しているが、本実験での護岸は湾曲部を含むその上下流のみ なので、Figure 60 に示すように[固定床と移動床]をフォーカスして、 直線部の護岸部分(この例では格子番号101より上流の赤色の格子部分)を選んで右クリックして、 属性を[移動床]に変更すし、[OK]を押す。
Figure 60 : 固定床と移動床の属性変更
また下流端は固定床なので、Figure 61 に示すように拡大、回転して 最下流の格子属性を[固定床]に変更する。
Figure 61 : 最下流端の格子を固定床に変更
計算条件の設定
メインメニューから[計算条件]->[設定]で現れる[計算条件]ウィンドウの、 [グループ]において、 [ソルバー・タイプ]、[境界条件]、[時間]、[河床材料]を それぞれ、 Figure 62 、 Figure 63 , Figure 64 , Figure 65 のようにパラメータを設定する。
Figure 62 : 計算条件(ソルバー・タイプ)
Figure 63 : 計算条件(境界条件)
Figure 64 : 計算条件(時間)
Figure 65 : 計算条件(河床材料)
なお。Figure 63 の[境界条件]においては、 [上流端流量と下流端水位の時間変化]で[Edit]を押して現れる、 Figure 66 の[流量時間変化設定ウィンドウ]で時間と流量の関係を与える。
Figure 66 : 計算条件(流量の時間配分の設定)
計算条件の設定が終了したら、[計算条件]ウィンドウで[OK]を押す。
計算の実行
計算を実行する前に、メインメニューから[ファイル]->[名前を付けてポロジェクトに保存]を選択して、 新しいフォルダを作成し、そのフォルダを選択することによりプロジェクトを保存しておく。 ここでは、[Nays2DH_flow_bed]という名前で保存する。( Figure 67 )
Figure 67 : 計算プロジェクトの保存
メインメニューから[計算]->[実行]を選択すると、今保存したばかりなのにしつこく[保存しますか？] と聞かれるので(Figure 68 ) [はい]を選ぶと計算が開始される。(Figure 69 )
Figure 68 : 「保存しますか？」
Figure 69 : 「計算実行中」
計算が終了したら、メインメニューから[計算結果]->[保存]を選択して計算結果を保存しておく。
計算結果の表示
オブジェクトブラウザーから[iRICZone]->[スカラー]->[ElevationChange(m)]に☑マークを付け、 [ElevationChange(m)]を右クリックして[プロパティ]を選択して、[スカラー設定]を Figure 70 のように設定する。
Figure 70 : 「スカラー設定」
オブジェクトブラウザーから[ベクトル]->[Velocity(ms-1)]に☑マークを付け、 [ベクトル]を右クリックして[プロパティ]を選択して、[ベクトル設定]を Figure 71 のように設定する。
Figure 71 : 「ベクトル設定」
メインメニューに[タイムスケールバー]をゼロに戻し、[アニメーション]->[開始/停止]を選択する ( Figure 72 )
Figure 72 : 「アニメーションの実行」
Figure 73 のようにアニメーションが表示され、河床変動がほぼ定常に達していることが分かる。
Figure 73 : 「河床変動と流速ベクトルのアニメーション」
計算結果のエクスポート
計算で得られた河床形状を次節で行うNays2d+による準3次元流れの計算の境界条件に使用するために 計算結果をテキストファイルにエクスポートする。Figure 74 に示すように、メインメニューから、 [ファイル]->[エクスポート]->[計算結果]を選ぶ.
Figure 74 : 「計算結果のエクスポート(1)」
[計算結果のエクスポート]ウィンドウ( Figure 74 )が表示されるので、 [形式]を[地勢データファイル形式(tpo)]に設定する( Figure 75 )。
Figure 75 : 「計算結果のエクスポート(2)」
出力フォルダは任意の名前で、 [全タイムステップ]の前にあるチェックボックスの☑を外し、[開始][終了]を10,800に設定して、[OK]をクリックすると、計算結果のエクスポートが完了する( Figure 76 )。
Figure 76 : 「計算結果のエクスポート(3)」
エクスポートされた計算結果は、( Figure 77 )に示すように、水深、流速、流砂量、河床高。。。。 など種類別に様々なファイルに保存されているが、このうち次節の計算に使用するのは河床高のみであるので、 [Results_1_Elevation(m).tpo]という名前のファイル以外は不要なので消しても構わない。
Figure 77 : 「計算結果のエクスポート(3)」
Nays2d+による準3次元流れの計算
ソルバの選択
iRICの起動画面から、[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 78 )で[Nays2d+簡単に3次元の流れの計算が出来ます] を選択して[OK]を押す。
Figure 78 : 「ソルバーNays2d+の選択」
計算格子と河床形状のインポートとマッピング
格子のインポート
メインメニューから[インポート]->[格子]を選択し、前記のNays2DHのプロジェクトフォルダー [Nays2DH_flow_bed]の中の[Case1.cgn]を選ぶ。 Figure 79 のような警告が出るが、「余計なお世話です」と心の中で思いながら、 構わず[はい]をクリックすると、格子のインポートが完了する。 ( Figure 80 )
Figure 79 : 「警告」
Figure 80 : 「格子のインポート完了」
河床高のインポート
メインメニューから[インポート]->[地理情報]->[河床高]を選択する( Figure 81 ).
Figure 81 : 「河床高のインポート」
インポートファイルの選択画面が現れる(Figure 82 )ので、前節のNays2dHの計算結果として エクスポートした [Results_1_Elevation(m).tpo]を選択して[開く]。
Figure 82 : 「河床高のインポート(ファイルの選択)」
Figure 83 のようにデータを間引くかどうか来かれるが、特に間引く必要がなければ、そのまま [OK]を押すと[河床高]のインポートが完了する( Figure 84 )
Figure 83 : 「河床高のインポート(間引き設定)」
Figure 84 : 「河床高のインポート完了」
マッピング
インポートした河床高データをインポートした格子上にマッピングを行う。
Figure 85 のように[格子]->[属性のマッピング]->[実行]を選ぶ。
Figure 85 : 「マッピング」
Figure 86 マッピングする[地理情報]を聞かれるので、 [河床高(m)]に☑を入れて[OK]をクリックすると
Figure 86 : 「マッピングする情報の選択」
マッピングが完了する( Figure 87 )ので、[OK]をクリックして 完了。
Figure 87 : 「マッピングの完了」
Nays2d＋の計算条件の設定
メインメニューから[計算条件]->[設定]で現れる[計算条件]ウィンドウの、 [グループ]において、 [流量および下流端水位の設定]、[時間および浸食に関するパラメータパラメータ]、 [境界条件]、[他の計算条件]、[３次元流速分布]を それぞれ、 Figure 88 、 Figure 89 , Figure 90 , Figure 91 , Figure 92 のようにパラメータを設定する。
Figure 88 : 計算条件(流量および下流端水位の設定)
Figure 89 : 計算条件(時間および浸食に関するパラメーター)
Figure 90 : 計算条件(境界条件)
Figure 91 : 計算条件(他の計算条件)
Figure 92 : 計算条件(3次元流速分布)
なお。Figure 88 の[流量および下流端水位の設定]においては、 [流量と下流端水位の時系列]の隣の[Edit]を押して現れる、 Figure 93 の[流量時間変化設定ウィンドウ]で時間と流量の関係を与える。
Figure 93 : 計算条件(流量の時間配分の設定)
計算条件の設定が終了したら、[計算条件]ウィンドウで[OK]を押す。
Nays2d+の計算の実行
計算の実行方法は前節[Nays2DHによる計算の実行]と全く同じで、タダでさえ長いのに これ以上同じことをクドクド書くと嫌われそうなので省略する。 ただ、計算の実行の前に、必ず、プロジェクトを保存しておくことを推奨する。 ここでは、[Nays2d+Flow]とう名前のプロジェクトに保存する。
Figure 94 : プロジェクトの保存(Nays2d+Flow)
計算結果は[Case1.cgn]というCGNSファイルに保存されるが、次のGELATOで使用するのはこの フォルダに保存される[Case1.cgn]を使用する。 計算の実行が終わった時も必ず[計算結果]->[保存]で結果を保存すること。
GELATOによる仮想トレーサーの追跡計算
ソルバの選択
iRICの起動画面から、[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 95 )で[GELATO] を選択して[OK]を押す。
Figure 95 : 「ソルバーGELATOの選択」
格子のインポート
Figure 96 のようにオブジェクトブラウザーの[格子(データなし)]を右クリック して、[インポート]をクリックする。
Figure 96 : 「格子のインポート(1)」
ファイル選択ウィンドウが現れるので、先ほど[Nays2d+]の計算結果を保存したフォルダ [Nays2d+Flow]の中の[Case1.cgn]を選ぶ( Figure 97)
Figure 97 : 「格子のインポート(2)」
お馴染みの Figure 98 の警告が表示されるが、構わず[はい]を押して進むと 格子のインポートが完了する。( Figure 99)
Figure 98 : 「お馴染みの警告」
Figure 99 : 「格子のインポート完了」
GELATOによるトレーサー追跡計算
計算条件の設定
[計算条件]->[設定]で表示される[計算条件]ウィンドウで、
[基本設定]、[プライマリートレーサーの供給条件]、[通常トレーサーの時間設定]、[乱れの影響] をそれぞれ、 Figure 100 、Figure 101 、Figure 102、 Figure 103 のようにパラメータを設定する。ここではまず、[乱れの影響]を考慮しない場合のトレーサー追跡を行う。
Figure 100 : 「GELATOの計算条件(1)」
Figure 101 : 「GELATOの計算条件(2)」
Figure 102 : 「GELATOの計算条件(3)」
Figure 103 : 「GELATOの計算条件(4)」
なお、Figure 100 中の[流れの計算結果を読み込むCGNSファイル]は前節[Nays2d+による流れの計算] の結果を保存した[Nays2d+Flow]プロジェクトフォルダにある[Case1.cgn]を選択する。
計算の実行
メインメニューから[計算]->[実行]を選択すると、「プロジェクトの保存がお勧めされる」ので、 ここは大人しく言うことを聞いて、新たにプロジェクトを保存しておく ( Figure 104).
Figure 104 : 「GELATOプロジェクトの保存(1)」
[ファイルに保存(ipro)]か[プロジェクトとして保存]どちらでも良い。
Figure 105 : 「GELATOプロジェクトの保存(3)」
計算が始まるとお馴染みの Figure 106 この画面が登場し。終了すると、これまたお馴染みの Figure 107 が表示されるので、[OK]を押す。
Figure 106 : 「計算の実行(1)」
Figure 107 : 「計算の実行(2)」
計算結果の表示
メインメニューから[計算結果]->[新しい可視化ウィンドウ(2D)を開く]を選択すると、計算結果が表示される。 (Figure 108 )
Figure 108 : 「可視化ウィンドウ(2D)」
本章冒頭に示した実験の画像 Figure 49 と向きが反対なので、Figure 108 に 矢印で示す(小さく分かりにくいが)90°回転のマークを2回クリックして180°回転させる ( Figure 109).
Figure 109 : 「可視化ウィンドウ(180°回転)」
時刻表示が小さくてメチャ見にくいので、オブジェクトブラウザーの[時刻]->[プロパティ]で [時刻設定]を表示させて( Figure 110) フォントサイズを適当に大きくする。
Figure 110 : 「時刻表示設定」
Figure 111 に示すように[時刻バーを戻し]、メインメニューから [アニメーション]->[実行]でトレーサーの流動アニメーションが開始される ( Figure 112)
Figure 111 : 「アニメーションの開始」
Figure 112 : [トレーサーのアニメーション(乱流拡散強度A=0)」
トレーサーは殆ど拡散せずに、線状に流れて行くのが分かる。
乱流拡散強度の違いの比較
Figure 113 [計算条件]->[設定]の、[グループ][乱れの影響]で、 [ランダムウォークによるセル以下スケールの乱れ考慮]を[する]にして、[Aの値]を[1]にして 再度[計算を実行]、[計算結果のアニメーション表示を]すると、 Figure 114 のようになる。
Figure 113 : 「ランダムウォークパラメータ(A=1)の設定」
Figure 114 : [トレーサーのアニメーション(乱流拡散強度A=1)」
同様に、[A=5]、[A=10]、[A=50]でランダムウォークをやってみると。。。
Figure 115 : [トレーサーのアニメーション(乱流拡散強度A=5)」
Figure 116 : [トレーサーのアニメーション(乱流拡散強度A=10)」
Figure 117 : [トレーサーのアニメーション(乱流拡散強度A=50)」
Figure 49 の実験と比較すると、[A=10]位の感じであることが分かる。
トレーサーのクローン
Figure 118 [計算条件]->[設定]の、[グループ] [トレーサーのクローニング(分割)と再結合]の設定で[クリーニング]を[する]、 [方法の選択]を[トレーサーが1個のセルだけ新規トレーサーを発生させる]、 [最大クローニング世代数]を[20]に設定すし、 [乱れの影響]は[A=10]として再度計算を実行し、結果を表示する( Figure 119 )
Figure 118 : [トレーサークローンの設定」
Figure 119 : [トレーサークローン表示(最大20世代、A=10)」
トレーサーの拡散範囲は Figure 49 の実験動画の緑色染料の拡散範囲に近くなっている。 ここで、オブジェクトブラウザーで、[粒子]->[スカラー]->[Generations]に☑マークを入れると、世代が表示される。 これをアニメーション表示すると、 Figure 120 のようになる。
Figure 120 : [トレーサークローン表示(最大20世代、A=10、世代色別表示)」
基本事項(共通事項) で述べたように、実質の重みは10世代目で \(W=0.00195\)、20世代目だと \(W=0.00000195\) なので、 Figure 120 おける、緑・黄・赤等のトレーサー 濃度は中心部の青色系のトレーサーに比べ対数的に低いことになる。 格子内のトレーサー数に重みを乗じて、実質の濃度を見るには、
1．オブジェクトブラウザーで[スカラー]の☑マークを外す( Figure 121).
Figure 121 : [スカラー」の☑マークを外す
2．オブジェクトブラウザーで[スカラー(セル中心)][Weghted numbers of tracers]に☑マーク を入れる(Figure 122).
Figure 122 : [Weighted numbers of tracers]に☑マークを入れる
3．[Weighted numbers of tracers]を右クリックして[プロパティ]を押す。
Figure 123 : [Weighted numbers of tracers]->[プロパティ]
4．[スカラー設定]ウィンドウで、以下のように設定して[設定]を押す。
Figure 124 : スカラー設定
Figure 125 の[可視化ウィンドウ:2D]が表示されるので、 タイムバーをゼロに戻してメインメニューから[アニメーション]->[開始/停止]を 押すと、 Figure 126 のアニメーションが開始される。
Figure 125 : アニメーションの実行
Figure 126 : 重みを考慮したトレーサー濃度のアニメーション
Figure 49 の実験動画の緑色染料の拡散状況に類似した拡散状況が再現された。
トレーサークローンを利用した流れの可視化
トレーサーのクローンツールを用いた流れの可視化の例を示す。
[Weighted numbers of tracers]の☑マークを外し、メインメニューの[計算条件]->[設定]を開く。 Figure 127 および Figure 128 のように条件を設定して保存する。 ここで、Figure 128 の[トレーサーの無いセルはすべて発生させる] によって、可視化用の多数の トレーサーを発生さる。
Figure 127 : 計算条件の設定(1)
Figure 128 : 計算条件の設定(2)
設定終了後、計算を実行し、 オブジェクトブラウザーの[粒子]と[スカラー]の☑マークを入れて、[Generations]の☑マークを外した後に [アニメーション]->[開始/停止]で Figure 129 のようにアニメーションが表示され、水路全体にトレーサーが 満遍なく配置された可視化となる。
Figure 129 : トレーサーを用いた可視化
魚の遊泳シミュレーション
[計算条件]->[設定]で以下の設定を行う。
Figure 130 : 魚の条件設定(1)
Figure 131 : 魚の条件設定(2)
Figure 132 : 魚の条件設定(3)
Figure 133 : 魚の条件設定(4)
Figure 134 : 魚の条件設定(5)
この条件で、[計算]->[実行]を行った後に、オブジェクトブラウザーで[ポリゴン]->[Fish]->[Type]に☑マークを入れて [アニメーション]->[開始/停止]を選択すると、Figure 134 が再生される。
Figure 135 : アニメーションの設定
Figure 136 : 魚の遊泳アニメーション
NaysDw2による流木の追跡計算
本節では2次元流木追跡ソルバNayswd2により流木の追跡を行う。
ソルバの選択
iRICの起動画面から、[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 137 )で[NaysDw2シンプルな2次元流木追跡ツール] を選択して[OK]を押す。
Figure 137 : 「ソルバーNaysDw2の選択」
格子のインポート
Figure 138 のようにオブジェクトブラウザーの[格子(データなし)]を右クリック して、[インポート]をクリックする。
Figure 138 : 「格子のインポート(3)」
ファイル選択ウィンドウが現れるので、先ほど[Nays2d+]の計算結果を保存したフォルダ [Nays2d+Flow]の中の[Case1.cgn]を選ぶ( Figure 139)
Figure 139 : 「格子のインポート(4)」
Figure 98 の警告が表示されるが、構わず[はい]を押して進むと 格子のインポートが完了する。( Figure 141)
Figure 140 : 「警告」
Figure 141 : 「格子のインポート完了」
計算条件の設定
[計算条件]->[設定]で以下の設定を行う。 Figure 142 [基本設定]の[流況計算結果を読み込むファイル名]を選ぶ
Figure 142 : [基本設定]->[計算結果を読み込むファイル名](1)
Figure 143 前節[Nays2d+]の計算結果の[Case1.cgn]を選ぶ
Figure 143 : [Cgnsファイルの指定]
[基本設定]の他のパラーメタは Figure 144 のようにパラメータを設定する。
Figure 144 : [基本設定]
[流木の供給条件]パラーメタは Figure 145 のようにパラメータを設定する。
Figure 145 : [流木の供給条件]
[流れおよび流木に関する]パラーメタは Figure 146 のようにパラメータを設定する。
Figure 146 : [流れおよび流木に関するパラメーター]
[DEM(個別要素法)パラメーター]は Figure 147 のようにパラメータを設定して、最後に [OK]をクリックする。
Figure 147 : [DEM(個別要素法)パラメーター]
流木追跡計算の実行
Figure 148 メインメニューから[計算]->[実行]を選択。
Figure 148 : [計算]->[実行]
Figure 149 [プロジェクトを保存しますか?]と聞かれるので、[はい]を選んで保存する。
Figure 149 : [プロジェクトを保存しますか？]
[プロジェクトを保存方法の選択]が問われるので、ここでは [プロジェクトとして保存]を選択して[OK]を押し、保存するフォルダ(空のフォルダ)を指定して、[フォルダの選択]をを押す
計算が開始されると、Figure 150 の画面が表示され、終了すると Figure 151 が表示されるので[OK]をクリックする。
Figure 150 : [計算の実行]
Figure 151 : [計算の終了]
流木追跡計算結果の表示
Figure 152 メインメニューから[計算結果]->[新しい可視化ウィンドウ(2D)を開く]を選択
Figure 152 : [可視化ウィンドウの表示]
Figure 153 オブジェクトブラウザーで、[iRICZone][スカラー(格子点)][Res-Velocity]に☑マーク を入れて、[右クリック]->[プロパティ]を選択
Figure 153 : [スカラーの表示]
Figure 154 [スカラー設定]ウィンドウを図のように設定して[OK]を押す。
Figure 154 : [スカラー設定]
Figure 155 タイムバーをゼロに戻し、[アニメーション]->[開始/停止]を押す。
Figure 155 : [アニメーションの開始]
Figure 156 のようなアニメーションが表示される。
Figure 156 : [流木追跡のアニメーション]