[計算例 2] 単純な湾曲水路の浮遊物追跡

Figure 48 に示すような直線+湾曲+直線水路実験の流量解析およびトレーサーの追跡計算を行う． この実験は水路全幅約1.8mで，左岸側半分が高水敷，右岸側半分が低水路の複断面水路で，右岸側のみが 移動床となっている．平面形状および断面形状の概略を Figure 49 に示す． この実験は千代田水理実験の検討において開発局と寒地土木研究所が実施したものです．

Figure 48 : 実験状況の動画

Figure 49 : 実験水路の形状

以下の計算事例では以下の手順で計算を実施する．

1. Nays2DHにより河床変動計算を実施し，河床変動がほぼ定常に達すた状態の河床形状を得る

2. Nays2d+により準3次元流れ場の計算をする．

3. UTTによりトレーサーの追跡を行う．乱流拡散強度パラーメーターを変更してその影響を比較する．

Nays2DHによる流れと河床変動の計算

ソルバの選択

iRICの起動画面から，[新しいプロジェクト]を選ぶと表示されるソルバの選択画面で， [Nays2dH iRIC3x 1.0 64bit]を選んで[OK]ボタン押すと，

Figure 50 : ソルバーの選択

「無題- iRIC 3.x.xxxx [Nays2DH iRIC3X 1.0 64bit]」と書かれた Windowが現れる．

Figure 51 : 無題

計算格子の作成

[格子]->[格子生成アルゴリズムの選択]で現れるウィンドウで， [2次元単純円弧格子作成ツール(複断面対応版)]を選んで[OK]を押す． (Figure 52 )

Figure 52 : 格子生成アルゴリズムの選択

[格子生成]のグループ[水路形状]，[断面形状]，[追加水路]，[粗度と河床状態]を それぞれ， Figure 53 , Figure 54 , Figure 55 , Figure 56 のようにパラメータを設定して，最後に[格子生成]を押す．

Figure 53 : 格子生成(1)

Figure 54 : 格子生成(2)

Figure 55 : 格子生成(3)

Figure 56 : 格子生成(4)

「マッピングを実行しますか？」と聞かれるので[はい(Y)]を押す． ( Figure 57 )

Figure 57 : 確認

オブジェクトブラウザーの[格子][セルの属性][固定床と移動床]に☑マークを入れると， Figure 58 のように固定床部分が赤，移動床部分が青で示された 格子の図が示される．

Figure 58 : 固定床と移動床

低水路と高水敷の境界の固定床は護岸を想定しているが，本実験での護岸は湾曲部を含むその上下流のみ なので，Figure 59 に示すように[固定床と移動床]をフォーカスして， 直線部の護岸部分(この例では格子番号87より上流の赤色の格子部分)を選んで右クリックして， 属性を[移動床]に変更すし，[OK]を押す．

Figure 59 : 固定床と移動床の属性変更

また下流端は固定床なので，Figure 60 に示すように拡大，回転して 最下流の格子属性を[固定床]に変更する．

Figure 60 : 最下流端の格子を固定床に変更

計算条件の設定

メインメニューから[計算条件]->[設定]で現れる[計算条件]ウィンドウの， [グループ]において， [ソルバー・タイプ]，[境界条件]，[時間]，[河床材料]を それぞれ， Figure 61 , Figure 62 , Figure 63 , Figure 64 のようにパラメータを設定する．

Figure 61 : 計算条件(ソルバー・タイプ)

Figure 62 : 計算条件(境界条件)

Figure 63 : 計算条件(時間)

Figure 64 : 計算条件(河床材料)

なお．Figure 62 の[境界条件]においては， [上流端流量と下流端水位の時間変化]で[Edit]を押して現れる， Figure 65 の[流量時間変化設定ウィンドウ]で時間と流量の関係を与える．

Figure 65 : 計算条件(流量の時間配分の設定)

計算条件の設定が終了したら，[計算条件]ウィンドウで[保存して閉じる]を押す．

計算の実行

計算を実行する前に，メインメニューから[ファイル]->[名前を付けてポロジェクトに保存]を選択して， 新しいフォルダを作成し，そのフォルダを選択することによりプロジェクトを保存しておく． ここでは，[Nays2DH_flow_bed]という名前で保存する．( Figure 66 )

Figure 66 : 計算プロジェクトの保存

メインメニューから[計算]->[実行]を選択すると，今保存したばかりなのにしつこく[保存しますか？] と聞かれるので(Figure 67 ) [はい]を選ぶと計算が開始される．(Figure 68 )

Figure 67 : 「保存しますか？」

Figure 68 : 「計算実行中」

計算が終了したら，メインメニューから[計算結果]->[保存]を選択して計算結果を保存しておく．

計算結果の表示

オブジェクトブラウザーから[iRICZone]->[スカラー]->[ElevationChange(m)]に☑マークを付け， [ElevationChange(m)]を右クリックして[プロパティ]を選択して，[スカラー設定]を Figure 69 のように設定する．

Figure 69 : 「スカラー設定」

オブジェクトブラウザーから[ベクトル]->[Velocity(ms-1)]に☑マークを付け， [ベクトル]を右クリックして[プロパティ]を選択して，[ベクトル設定]を Figure 70 のように設定する．

Figure 70 : 「ベクトル設定」

メインメニューに[タイムスケールバー]をゼロに戻し，[アニメーション]->[開始/停止]を選択する ( Figure 71 )

Figure 71 : 「アニメーションの実行」

Figure 72 のようにアニメーションが表示され，河床変動がほぼ定常に達していることが分かる．

Figure 72 : 「河床変動と流速ベクトルのアニメーション」

計算結果のエクスポート

計算で得られた河床形状を次節で行うNays2d+による準3次元流れの計算の境界条件に使用するために 計算結果をテキストファイルにエクスポートする．Figure 73 に示すように，メインメニューから， [ファイル]->[エクスポート]->[計算結果]を選ぶ.

Figure 73 : 「計算結果のエクスポート(1)」

[計算結果のエクスポート]ウィンドウ( Figure 73 )が表示されるので， [形式]を[地勢データファイル形式(tpo)]に( Figure 74 ),

Figure 74 : 「計算結果のエクスポート(2)」

出力フォルダは任意の名前で， [全タイムステップ]の前にあるチェックボックスの☑を外し，[開始][終了]を7,200の設定して． [OK]をクリックすると，計算結果のエクスポートが完了する( Figure 75 )．

Figure 75 : 「計算結果のエクスポート(3)」

エクスポートされた計算結果は，( Figure 76 )に示すように，水深，流速，流砂量，河床高．．．． など種類別に様々なファイルに保存されているが，このうち次節の計算に使用するのは河床高のみであるので， [Results_1_Elevation(m).tpo]という名前のファイル以外は不要なので消しても構わない．

Figure 76 : 「計算結果のエクスポート(3)」

Nays2d+による準3次元流れの計算

ソルバの選択

iRICの起動画面から，[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 77 )で[Nays2d+簡単に3次元の流れの計算が出来ます] を選択して[OK]を押す．

Figure 77 : 「ソルバーNays2d+の選択」

計算格子と河床形状のインポートとマッピング

格子のインポート

メインメニューから[インポート]->[格子]を選択し，前記のNays2dHのプロジェクトフォルダー [Nays2DH_flow_bed]の中の[Case1.cgn]を選ぶ． Figure 78 のような警告が出るが，「余計なお世話です」と心の中で思いながら， 構わず[はい]をクリックすると，格子のインポートが完了する． ( Figure 79 )

Figure 78 : 「警告」

Figure 79 : 「格子のインポート完了」

河床高のインポート

メインメニューから[インポート]->[地理情報]->[河床高]を選択する( Figure 80 ).

Figure 80 : 「河床高のインポート」

インポートファイルの選択画面が現れる(Figure 81 )ので，前節のNays2dHの計算結果として エクスポートした [Results_1_Elevation(m).tpo]を選択して[開く]．

Figure 81 : 「河床高のインポート(ファイルの選択)」

Figure 82 のようにデータを間引くかどうか来かれるが，特に間引く必要がなければ，そのまま [OK]を押すと[河床高]のインポートが完了する( Figure 83 )

Figure 82 : 「河床高のインポート(間引き設定)」

Figure 83 : 「河床高のインポート完了」

マッピング

インポートした河床高データをインポートした格子上にマッピングを行う．

Figure 84 のように[格子]->[属性のマッピング]->[実行]を選ぶ．

Figure 84 : 「マッピング」

Figure 85 マッピングする[地理情報]を聞かれるので， [河床高(m)]に☑を入れて[OK]をクリックすると

Figure 85 : 「マッピングする情報の選択」

マッピングが完了する( Figure 86 )ので，[OK]をクリックして 完了．

Figure 86 : 「マッピングの完了」

Nays2d＋の計算条件の設定

メインメニューから[計算条件]->[設定]で現れる[計算条件]ウィンドウの， [グループ]において， [流量および下流端水位の設定]，[時間および浸食に関するパラメータパラメータ]， [境界条件]，[他の計算条件], [３次元流速分布]を それぞれ， Figure 87 , Figure 88 , Figure 89 , Figure 90 , Figure 91 のようにパラメータを設定する．

Figure 87 : 計算条件(流量および下流端水位の設定)

Figure 88 : 計算条件(時間および浸食に関するパラメーター)

Figure 89 : 計算条件(境界条件)

Figure 90 : 計算条件(他の計算条件)

Figure 91 : 計算条件(3次元流速分布)

なお．Figure 87 の[流量および下流端水位の設定]においては， [流量と下流端水位の時系列]の隣の[Edit]を押して現れる， Figure 92 の[流量時間変化設定ウィンドウ]で時間と流量の関係を与える．

Figure 92 : 計算条件(流量の時間配分の設定)

計算条件の設定が終了したら，[計算条件]ウィンドウで[保存して閉じる]を押す．

Nays2d+の計算の実行

計算の実行方法は前節[Nays2DHによる計算の実行]と全く同じで，タダでさえ長いのに これ以上同じことをクドクド書くと嫌われそうなので省略する． ただ，計算の実行の前に，必ず，プロジェクトを保存しておくことを推奨する． ここでは，[Nays2d+Flow]とう名前のプロジェクトに保存する．

Figure 93 : プロジェクトの保存(Nays2d+Flow)

計算結果は[Case1.cgn]というCGNSファイルに保存されるが，次のUTTで使用するのはこの フォルダに保存される[Case1.cgn]を使用する． 計算の実行が終わった時も必ず[計算結果]->[保存]で結果を保存すること ( Figure 94).

Figure 94 : 計算結果の保存(お忘れなく!)

UTTによる仮想トレーサーの追跡計算

ソルバの選択

iRICの起動画面から，[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 95 )で[UTTとても便利な多機能トレーサー追跡ツール] を選択して[OK]を押す．

Figure 95 : 「ソルバーUTTの選択」

格子のインポート

Figure 96 のようにオブジェクトブラウザーの[格子(データなし)]を右クリック して，[インポート]をクリックする．

Figure 96 : 「格子のインポート(1)」

ファイル選択ウィンドウが現れるので，先ほど[Nays2d+]の計算結果を保存したフォルダ [Nays2d+Flow]の中の[Case1.cgn]を選ぶ( Figure 97)

Figure 97 : 「格子のインポート(2)」

お馴染みの Figure 98 の警告が表示されるが，構わず[はい]を押して進むと 格子のインポートが完了する．( Figure 99)

Figure 98 : 「お馴染みの警告」

Figure 99 : 「格子のインポート完了」

UTTによるトレーサー追跡計算

計算条件の設定

[計算条件]->[設定]で表示される[計算条件]ウィンドウで，

[基本設定]，[通常トレーサーの供給条件], [乱れの影響] をそれぞれ， Figure 100 , Figure 101 Figure 102 のようにパラメータを設定する．ここではまず，[乱れの影響]を考慮しない場合のトレーサー追跡を行う．

Figure 100 : 「UTTの計算条件(1)」

Figure 101 : 「UTTの計算条件(2)」

Figure 102 : 「UTTの計算条件(3)」

なお，Figure 100 中の[流れの計算結果を読み込むCGNSファイル]は前節[Nays2d+による流れの計算] の結果を保存した[Nays2d+Flow]プロジェクトフォルダにある[Case1.cgn]を選択する． ( Figure 103)

Figure 103 : 「流れの計算結果を読み込むCGNSファイルの選択」

計算の実行

メインメニューから[計算]->[実行]を選択すると，「プロジェクトの保存がお勧めされる」ので， ここは大人しく言うことを聞いて，新たにプロジェクトを保存しておく ( Figure 104).

Figure 104 : 「UTTプロジェクトの保存(1)」

Figure 105 [ファイルに保存(ipro)]か[プロジェクトとして保存]どちらでも良いが， ここでは[UTT1.ipro]という名前で保存しておく．

Figure 105 : 「UTTプロジェクトの保存(3)」

計算が始まるとお馴染みの Figure 106 この画面が登場し．終了すると，これまたお馴染みの Figure 107 が表示されるので，[OK]を押す．

Figure 106 : 「計算の実行(1)」

Figure 107 : 「計算の実行(2)」

計算結果の表示

メインメニューから[計算結果]->[新しい可視化ウィンドウ(2D)を開く]を選択すると，計算結果が表示される． (Figure 108 )

Figure 108 : 「可視化ウィンドウ(2D)」

本章冒頭に示した実験の画像 Figure 48 と向きが反対なので，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 48 の実験と比較すると，[A=10]位の感じであることが分かる．

トレーサーのクローン

Figure 118 [計算条件]->[設定]の，[グループ] [トレーサーのクローニング(分割)と再結合]の設定で[クリーニング]を[する], [方法の選択]を[トレーサーが1個のセルだけ新規トレーサーを発生させる]， [最大クローニング世代数]を[20]に設定すし， [乱れの影響]は[A=10]として再度計算を実行し，結果を表示する( Figure 119 )

Figure 118 : [トレーサークローンの設定」

Figure 119 : [トレーサークローン表示(最大20世代，A=10)」

トレーサーの拡散範囲は Figure 48 の実験動画の緑色染料の拡散範囲に近くなっている． トレーサーの数が莫大に見えるが，実際は，Figure 120 に示すように， オブジェクトブラウザーで，[粒子]->[スカラー]->[Generations]に☑マークを入れると， 世代が表示される．

Figure 120 : [クローンの世代色別表示」

これをアニメーション表示すると， Figure 121 のようになる．

Figure 121 : [トレーサークローン表示(最大20世代，A=10, 世代色別表示)」

基本事項(共通事項) で述べたように，実質の重みは10世代目で \(W=0.00195\), 20世代目だと \(W=0.00000195\) なので， Figure 121 おける，緑・黄・赤等のトレーサー 濃度は中心部の青色系のトレーサーに比べ対数的に低いことになる． 格子内のトレーサー数に重みを乗じて，実質の濃度を見るには，

1. オブジェクトブラウザーで[スカラー]の☑マークを外す( Figure 122).

Figure 122 : [スカラー」の☑マークを外す

2. オブジェクトブラウザーで[スカラー(セル中心)][Weghted numbers of tracers]に☑マーク を入れる(Figure 123).

Figure 123 : [Weighted numbers of tracers]に☑マークを入れる

3. [Weighted numbers of tracers]を右クリックして[プロパティ]を押す．

Figure 124 : [Weighted numbers of tracers]->[プロパティ]

4. [スカラー設定]ウィンドウで，[自動]の☑マークを外し，[最大値]を[0.1], [最小値]を[1e-08]に設定し，[最小値以下を描画]の☑マークを外し， [カラーマーク]を[手動]にして[設定]を押す．

Figure 125 : スカラー設定

5. [カスタムカラーの設定]が表示されるので，[種類]を[3色]にして． [最大値]を[濃い緑]，[中間値]を[薄目の緑]，[中間値の値]を[1e-06]にして[OK]を押す．

Figure 126 : カスタムカラー

6. [スカラー設定]に戻るので[OK]を押して設定終了．

Figure 127 : スカラー設定（設定終了)

Figure 128 の[可視化ウィンドウ:2D]が表示されるので， タイムバーをゼロに戻してメインメニューから[アニメーション]->[開始/停止]を 押すと， Figure 129 のアニメーションが開始される．

Figure 128 : アニメーションの実行

Figure 129 : 重みを考慮したトレーサー濃度のアニメーション

Figure 48 の実験動画の緑色染料の拡散状況に類似した拡散状況が再現された．

トレーサークローンを利用した流れの可視化

トレーサーのクローンツールを用いた流れの可視化の例を示す．[計算条件]を Figure 130,

[Weighted numbers of tracers]の☑マークを外し，メインメニューの[計算条件]->[設定]を開く． Figure 130 および Figure 131 のように条件を設定して保存する． ここで，Figure 131 の[トレーサーの無いセルはすべて発生させる] によって，可視化用の多数の トレーサーを発生さる．

Figure 130 : 計算条件の設定(1)

Figure 131 : 計算条件の設定(2)

設定終了後，計算を実行し， オブジェクトブラウザーの[粒子]と[スカラー]の☑マークを入れて，[Generations]の☑マークを外した後に [アニメーション]->[開始/停止]で Figure 132 のようにアニメーションが表示され，水路全体にトレーサーが 満遍なく配置された可視化となる．

Figure 132 : トレーサーを用いた可視化

魚の遊泳シミュレーション

[計算条件]->[設定]で以下の設定を行う．

Figure 133 : 魚の条件設定(1)

Figure 134 : 魚の条件設定(2)

Figure 135 : 魚の条件設定(3)

この条件で，[計算]->[実行]を行った後に，オブジェクトブラウザーで[ポリゴン]->[Fish]->[Type]に☑マークを入れて [アニメーション]->[開始/停止]を選択すると, Figure 137 が再生される．

Figure 136 : アニメーションの設定

Figure 137 : 魚の遊泳アニメーション

NaysDw2による流木の追跡計算

本節では2次元流木追跡ソルバNayswd2により流木の追跡を行う．

ソルバの選択

iRICの起動画面から，[新しいプロジェクト]を選ぶと表示されるソルバの選択画面 ( Figure 138 )で[NaysDw2シンプルな2次元流木追跡ツール] を選択して[OK]を押す．

Figure 138 : 「ソルバーNaysDw2の選択」

格子のインポート

Figure 139 のようにオブジェクトブラウザーの[格子(データなし)]を右クリック して，[インポート]をクリックする．

Figure 139 : 「格子のインポート(3)」

ファイル選択ウィンドウが現れるので，先ほど[Nays2d+]の計算結果を保存したフォルダ [Nays2d+Flow]の中の[Case1.cgn]を選ぶ( Figure 140)

Figure 140 : 「格子のインポート(4)」

Figure 98 の警告が表示されるが，構わず[はい]を押して進むと 格子のインポートが完了する．( Figure 142)

Figure 141 : 「警告」

Figure 142 : 「格子のインポート完了」

計算条件の設定

[計算条件]->[設定]で以下の設定を行う． Figure 143 [基本設定]の[流況計算結果を読み込むファイル名]を選ぶ

Figure 143 : [基本設定]->[計算結果を読み込むファイル名](1)

Figure 144 前節[Nays2d+]の計算結果の[Case1.cgn]を選ぶ

Figure 144 : [Cgnsファイルの指定]

[基本設定]の他のパラーメタは Figure 145 のようにパラメータを設定する．

Figure 145 : [基本設定]

[流木の供給条件]パラーメタは Figure 146 のようにパラメータを設定する．

Figure 146 : [流木の供給条件]

[流れおよび流木に関する]パラーメタは Figure 147 のようにパラメータを設定して，最後に [保存して閉じる]をクリックする．

Figure 147 : [流れおよび流木に関するパラメーター]

流木追跡計算の実行

Figure 148 メインメニューから[計算]->[実行]を選択．

Figure 148 : [計算]->[実行]

Figure 149 [プロジェクトを保存しますか?]と聞かれるので，[ハイ]を選んで保存する．

Figure 149 : [プロジェクトを保存しますか？]

Figure 150 [プロジェクトを保存方法の選択]が問われるので，ここでは [プロジェクトとして保存]を選択して[OK]を押す

Figure 150 : [プロジェクトの選択方法の指定]

Figure 151 保存するフォルダ(空のフォルダ)を指定して． [フォルダの選択]をを押す

Figure 151 : [プロジェクトの選択]

計算が開始されると，Figure 152 の画面が表示され，終了すると Figure 153 が表示されるので[OK]をクリックする．

Figure 152 : [計算の実行]

Figure 153 : [計算の終了]

流木追跡計算結果の表示

Figure 154 メインメニューから[計算結果]->[新しい可視化ウィンドウ(2D)を開く]を選択

Figure 154 : [可視化ウィンドウの表示]

Figure 155 オブジェクトブラウザーで，[iRICZone][スカラー(格子点)][Res-Velocity]に☑マーク を入れて，[右クリック]->[プロパティ]を選択

Figure 155 : [スカラーの表示]

Figure 156 [スカラー設定]ウィンドウを図のように設定して[OK]を押す．

Figure 156 : [スカラー設定]

Figure 157 タイムバーをゼロに戻し，[アニメーション]->[開始/停止]を押す．

Figure 157 : [アニメーションの開始]

Figure 158 のようなアニメーションが表示される．

Figure 158 : [流木追跡のアニメーション]