よくある質問

葦のような植生は流れによって倒伏するため、抵抗係数 Cd と密生度で表現するより、マニングの粗度係数として考慮するほうが適切とされています。倒伏しない円柱状の樹木として扱う場合は、ソルバーマニュアルの参考文献、水理公式集、河道計画の手引きなどを参照して Cd を決める必要があります。既定値の 0.7 は過去の水理実験に基づく値で、これらの文献のいずれかに由来すると考えられます。植生抵抗の理論的背景は iRIC Academy で解説されています。

Vorticity（渦度）は平面二次元の流速ベクトルから計算される回転の強さを表す量で、時計回りの回転が負、反時計回りが正になります。二次流（断面内の螺旋流）とは別のものです。値が大きい場所ほど平面的な渦運動が強いことを示し、例えばカルマン渦の計算例では渦の発生を可視化するのに使われます。渦度や二次流の理論的背景は iRIC Academy で学ぶことができます。

Nays2DH では、格子の左右岸側の境界からの自由流出は扱えません。側壁を自由流出にしたい場合は Nays2DFlood を使う方法があります。Nays2DFlood では境界条件で側壁を自由流出に設定できます。流れの計算のみであれば両ソルバーの扱いはほぼ同じなので、河床変動を伴わない流況解析なら Nays2DFlood で試してみてください。

iRIC の riv 変換では、国土交通省フォーマット横断データの 1 行目のうち 断面番号のみ を使用し、水際杭高・堤防法尻高・堤内地平均などの値は無視されます。したがって、それらの値が低水路・高水敷の設定に反映されることはありません。低水路と高水敷を区分する値の入力も不要で、1 行目には測点数があれば十分です。急峻な山地河川などで低水路・高水敷の区分がない場合も、特別な設定をせずにそのまま使用できます。

格子をシェープファイルとして直接出力する機能はありませんが、格子と格子属性を CSV 形式でエクスポートして GIS に取り込む方法があります。手順は次のとおりです。

1. オブジェクトブラウザーの「格子」を右クリック
2. 「エクスポート」を選択
3. ファイルの種類で「iRIC 格子 CSV」を選ぶ
4. 保存する

出力された CSV には格子点座標と属性が含まれるので、QGIS などで CSV を点データとして読み込み、必要に応じてシェープファイルに変換してください。

グラフウィンドウの CSV エクスポート機能を使います。手順は次のとおりです。

1. グラフウィンドウを開き、縦断図なら X 軸に I、横断図なら X 軸に J を指定して描画します。
2. 対象とする断面を調整します（縦断図なら J の値、横断図なら I の値を指定）。
3. グラフウィンドウ右下の「CSVエクスポート」をクリックして保存します。

出力された CSV には、指定した時刻・断面の格子点ごとの値が含まれるので、表計算ソフトなどでプロットや平均値の算出ができます。詳細は iRIC GUI ユーザーマニュアルの「グラフウィンドウ」の項に記載されています。

「河川測量データから生成」アルゴリズムでは、中心線の位置に応じて測線間の補間が行われるため、中心線を変えると生成される格子が変わり、当然ながら計算結果も変わります。ただし「正しい中心線」というものは存在せず、横断測量データは測線間の情報が欠けているため、利用できる測量データと現地状況を踏まえて位置を決めていくしかありません。

横断測量データを地形データとして使う場合でも、必ず「河川測量データから生成」を使う必要はなく、「折れ線と格子幅から生成」など他の格子生成アルゴリズムも利用できます。格子作成にはやや試行錯誤が必要なので、複数のアルゴリズムを試して、隣接格子のサイズ差が小さく計算が安定する格子を作ることを推奨します。

Nays2DH の上流端の土砂供給は、基本的に動的平衡状態（上流端断面の平衡流砂量と同じ量の土砂が供給される）として扱われます。そのうえで、計算条件で平衡流砂量に対する供給土砂量の割合を設定でき、この値を調整することで供給量を増減させることができます（ソルバーマニュアルの該当項を参照）。

一方、上流端から供給される土砂の粒径を直接指定する機能は現状ありません。粒径を変えたい場合は、上流端付近の計算格子だけ粒度分布を変えることで近似的に考慮できる可能性があります。

水面の凸凹と遠心力には関係ありません。上流端付近で発生している場合、計算が不安定になる条件、特に上流端条件が厳しい可能性があります。

Nays2DH の上流端条件は「等流」などから選択できますが、等流以外を選んでも概ね等流を想定した扱いになります。そのため、次の点を確認してください。

• 上流端の位置が、等流条件から大きくかけ離れた場所（急な断面変化や湾曲部など）になっていないか
• 上流端付近の河床が逆勾配になっていないか

該当する場合は、上流端の位置をより等流に近い区間へ移す、あるいは地形を見直すことで改善する可能性があります。

このメッセージは、設定した下流端水位が下流端の河床高より低いことを示しています。下流端水位は、下流端断面の河床高よりも高い値にする必要があります。

計算条件の境界条件で下流端水位の設定（一定値または時系列）を見直し、河床高を上回る値に修正してください。時系列で与えている場合は、全時刻で河床高を下回っていないかを確認してください。

Nays2DH は平面二次元モデルで、格子の左右岸境界はスリップ条件として扱われます。計算領域内の左右岸境界セルの水位が、計算領域の外側にもそのまま続いているという想定です。

したがって、水位が横断端点より上がった場合に領域外へ水が広がる計算は行われません。そのような状況が想定される場合は、格子の範囲を十分に広く取る必要があります。境界条件の理論的な背景については iRIC Academy で解説しています。

グラフィックボードのドライバが原因である可能性があります。同様の症状が、ドライバを最新版に更新することで解決した事例があります。

その事例では Intel のオンボードグラフィックスを使用しており、Windows Update ではなく Intel の更新ツールで最新版（ベータ版）を検索・適用したところ解消しました。まずはお使いの GPU メーカーのツールでドライバを最新化して確認してください。

流速の x 成分・y 成分をコンターとして可視化することはできません。値を確認したい場合は、次のいずれかの方法があります。

• CSV に出力して、成分ごとの値を確認する。
• 計算結果の可視化ウィンドウで、オブジェクトブラウザの「スカラー（格子点）」を右クリックし、「属性ブラウザの表示」を選びます。属性ブラウザが表示された状態でコンター上にカーソルを合わせると、その点の計算結果の値（x・y 方向の流速を含む）が数値で表示されます。

Nays2DH では、計算領域の途中から意図的に水が抜けるような条件は設定できません。

途中での流出・流入を扱いたい場合は、Nays2DFlood であればポンプなどの設定で対応できる可能性があります。ただし具体的な計算内容によって適否は変わるため、目的に応じてソルバーの選択を検討してください。

合流点は 1 つのみ作成でき、複数の合流点は設定できません。

合流点を含む格子の作成から計算までの手順については、iRIC の解説動画が公開されているので参考にしてください。

テンプレートマッピングが無効になっているわけではありません。開発側でソースを確認し、実際に点群をインポートして TIN マッピングとテンプレートマッピングを切り替えたところ、結果に違いが生じることが確認されています。

両者で大きな違いが出るのは、点群の分布が不均一な場合、たとえば横断方向には密でも縦断方向にはまばらにしか点がない場合などです。レーザー測量などで取得した、縦横どちらの方向にも十分な密度がある点群では、どちらの手法でもほぼ同じ結果になることがあります。点群の密度が高い場合は、結果が一致していても異常ではありません。

方法は 2 つあります。

• 生成した格子の格子点座標は iRIC 上で直接編集できるので、背景図を見ながら格子点を移動させて局所的に密にすることができます。
• 「折れ線と格子幅から生成」ではなく「汎用格子生成ツール」を使うと、区間ごとに分割数や間隔を細かく制御できるため、うまく切れる可能性があります。このツールは機能が多いので、実際に使いながら機能を確認してみてください。

現状の仕様が原因である可能性が高いです。

• 点群から作成される TIN は、点群データセットごと（読み込んだファイルごと）に個別に作成されます。
• マッピングはオブジェクトブラウザの下側にあるデータから順に上書きされ、最も上側のデータが優先されます。

そのため、分割した点群同士が上空から見て重なっていると、意図したマッピングになりません。点群データを 1 つのファイルにまとめてから再度読み込むと、想定どおりにマッピングできるはずです。なお、この仕様は開発側でも問題と認識されており、将来のバージョンで改良される予定です。

並列計算で長時間の計算を行うと、エラーメッセージが出ないまま CPU 稼働が途中で 0 になって止まることがあります。その場合は、計算条件の並列化に関する設定で CPU 数を 1（並列計算しない）にすると解決する可能性があります。

並列化を外しても改善しない場合は、タイムステップや格子サイズなど計算の安定性に関わる条件も合わせて見直してください。

上流端では等流計算が行われ、設定した流量が流れるような範囲に自動的に流れが与えられます。そのため、基本的には計算上妥当な範囲に流量が配分されていると考えてください。幅を直接指定する設定はありません。

強制的に狭い範囲だけに流れを与えたい場合は、上流端の流したくない部分のセルに構造物（障害物）を設定する方法があります。構造物を設定したセルには流れが発生しなくなるため、実質的に流入幅を制限できます。

流下時間を直接計測する機能は Nays2DH にはありません。代わりに、可視化ウィンドウでパーティクル（粒子）を投入し、上流から下流へ到達するまでの時間を追う方法があります。パーティクルの投入位置は格子点上であれば任意に指定できるため、任意の2点間の流下時間の目安として利用できます。ただし、計算結果の出力間隔が大きすぎるとパーティクルの追跡がうまくいかない場合があるため、出力間隔を小さめに設定してから試すことをおすすめします。

実河川のように川幅や勾配が場所ごとに変わる場合、基本的には不等流計算が適しています。ただし Nays2DH の不等流計算は、断面形が複雑な場合や常流・射流が混在する場合に耐えられるほど頑健には作られていません。そのため、不等流計算で初期水面形をうまく作れない場合や、計算開始直後に発散してしまう場合には、いったん等流計算で初期条件を与え、その影響が消えるまで慣らし計算を行うという方法が用意されています。等流計算は初期条件のためのものなので、全断面で計算開始時と同じ流量を用いて等流水深を求めています。

Nays2DH は OpenMP で並列化されており、複数コアを使って計算できます。設定手順は次のとおりです。

1. 計算条件の設定でソルバータイプを「アドバンス（Advanced）」にする
2. 「+その他（Others）」のグループを開く
3. 「並列計算のスレッド数」に使用するコア数を入力する

なお、PC の最大コア数を指定しても最速になるとは限りません。計算規模によって適切なスレッド数は異なるため、いくつか試して比較することをおすすめします。

エラー表示なしにソルバーが止まったように見えて実際には進まない現象は、複数の CPU コアを使った並列計算で発生することが報告されています。この状態になるとプロジェクトファイルが壊れてしまい、途中までの結果を読み込めないことがあります。可能であれば、並列計算のスレッド数を 1 にして単一コアで実行してみてください。また、ウイルス対策ソフトの干渉やディスク容量不足など、実行環境に起因する場合もあることが分かっており、すべてのケースが解決できているわけではありません。長期間の計算では、こまめにプロジェクトを保存しておくことも有効です。

できます。格子生成アルゴリズムとして「横断測量データから生成」を使う場合は、次の手順で横断面ウィンドウ上に格子点の生成位置が黄色い四角で表示されます。

1. 河川測量データをインポートする
2. 横断面ウィンドウを開く
3. オブジェクトブラウザで「格子生成条件」を選び、右クリックメニューから「格子生成アルゴリズムの選択」を選ぶ
4. 「横断測量データから生成」を選択して OK を押す

格子点の位置が表示されるのは、プリプロセッサで「格子生成条件」が選択されている間だけで、それ以外のときは標高点が赤い点で表示されます。

他のアルゴリズムで生成した格子の横断形状を確認したい場合は、格子生成後にオブジェクトブラウザで「格子」→「格子形状」を選び、確認したい位置の格子点をクリックして選択し、右クリックメニューの「横断面ウィンドウを開く」を選ぶと、その点を含む横断線の形状が表示されます。

全く堆積していない固定床では本来「平均粒径」という概念がありませんが、Nays2DH ではその値をゼロにせず、設定した粗度に対応する粗度高さのような値を出力しています。そのため固定床部分でも 0 以外の平均粒径が表示されます。また、ごくわずかでも堆積が生じている場合は粒径の計算が行われるため、時間的に値が変化しているように見える可能性があります。値が時間的に変化しているかどうかを確認し、微小な堆積の有無と合わせて判断してください。

「等流計算」を選ぶと、上流端で与えた各時刻の流量と下流端断面の標高データから等流計算を行い、河床高と等流水深の和を下流端水位として与えます。「自由流出」は、下流端のひとつ上流側の断面の水位を下流端にコピーし、それを境界条件として計算する方式です。

使い分けの目安として、等流水深は上流端の流量に対応して決まるため、河川延長が長い場合には上流の流量変化が下流に伝わる前に下流端水位だけが先に変化してしまうことがあります。そのような場合は自由流出も選択肢になります。なお、河川延長が比較的短いケースでは両者でほとんど同じ水位になるという報告もあります。

現状の iRIC には、流速や水深が指定範囲にあるセルの面積を集計する機能はありません。計算結果を CSV などにエクスポートし、別途プログラムを組んで集計する必要があります。参考として、Nays2DH の Python 版の解説とソースコードが iRIC のサイト（i-ric.org/yasu/nbook2/）で公開されているので、コード内で面積計算を行いたい場合はこれを適宜変更して利用する方法もあります。

浮遊砂量そのものを直接出力する項目はありませんが、浮遊砂濃度・流速・水深は結果として出力されています。これらとメッシュ（格子）の幅を掛け合わせれば浮遊砂量を求めることができるので、計算結果をエクスポートして別途計算する方法を検討してください。

iRIC 3.0.19 リビジョン 6747 以前で、Open Street Map と国土地理院地図が表示されなくなる不具合が発生していました。原因は iRIC が地図表示に使用している外部 Web サービス側の変更で、iRIC 3.0.19 リビジョン 6750 で修正されています。同じ症状が出る場合は、まず iRIC を 3.0.19 リビジョン 6750 以降に更新してください。なお、地図タイルは外部サービスに依存しているため、その後も表示が変わる可能性があります。

iRIC の GUI から Nays2DH の上流端供給土砂量を時系列などで任意に与える機能は用意されていません。対応方法としては、Nays2DH の Python 版の解説とソースコード（i-ric.org/yasu/nbook2/ で公開）を利用し、上流端の流砂量の境界条件部分をコードで自由に調整する方法が案内されています。また、iRIC-UC 会員向けの有料サポート（uc.i-ric.org）で個別に相談することもできます。

現在のところ、Nays2DH には粒径別の流砂量を出力する機能はありません。出力されるのは掃流砂・浮遊砂といった区分の流砂量までです。粒径別の値が必要な場合は、ソースコードを改変して独自に出力を追加する必要があります。

地形データとして x, y, z の点群データがあれば計算は可能です。DTM から作った標高点や、岸線・中心線に沿った標高点も、x, y, z 形式のテキストファイルに変換すれば iRIC にインポートできます。手順は、オブジェクトブラウザの「地理情報」→「地形高」を右クリックし、「インポート」→「点群データ」を選ぶ方法です。横断線に沿っていない線状のデータも、点群として扱えば読み込めます。インポート後に格子を作成して計算に進んでください。QGIS で作成したデータについては、シェープファイルなどであればインポートできます。

固定床高さを地形高（移動床の高さ）とは別に与えるには、次の 2 段階の設定が必要です。

1. オブジェクトブラウザの「地理情報」で「固定床高さ」の地形データを読み込み（またはポリゴンなどで作成し）、格子にマッピングする
2. 計算条件の設定の「ソルバータイプ」の一番下にある固定床高さの与え方で、「固定床高さのデータを有効にする」を選択する

これにより任意の固定床高さを指定できます。多機能格子ツールで水路を作ると地形高には自動で勾配がつきますが、固定床高さは別途データとして与える必要があります。

点群データをインポートしたあとに「ポアソン方程式を解いて格子を作成」を選ぶと「横断測量データが見つかりません。中心線を自分で定義してください」と表示されますが、そのまま OK を押して中心線を引き、「左岸線・右岸線の生成」で岸線を作成し、「格子生成」で格子を作る手順で問題ありません。Nays2DH は構造格子を使うモデルなので、どのアルゴリズムで作成しても四角形の格子であれば計算できます。ただし、極端に変形した格子は計算精度が低下するため注意してください。格子の作り方については iRIC サイトのダウンロードページにある資料（iRIC notes）も参考になります。

計算条件の設定で、下流端水位を決める等流計算に使う勾配を「自動で計算する」にしている可能性があります。この設定では与えられた初期地形から勾配を算出するため、上流側の断面を変更しただけでも計算される勾配がわずかに変わり、その結果として下流端水位が変化します。現況と計画で下流端水位を揃えたい場合は、勾配を自動計算ではなく数値で直接与える設定にすると解決すると考えられます。

固定床高さの扱いには次のルールがあります。

• 地理情報のマッピングでは、ポリゴンを設定していない空白部分には各項目の初期値（固定床高さなら 0 m、固定床と移動床なら移動床、マニング粗度なら 0.03 など）が入ります。
• 一部にでも「固定床高さ」が設定されていると、「固定床と移動床」の設定に関わらず全体が固定床高さを持つ扱いになり、明示的に設定していない部分には固定床高さの最大値が入ります。
• 「地形高」が「固定床高さ」より低い場所では、その部分の固定床高さは地形高の値に更新されます。ただしこの更新はプリプロセッサの「格子点の属性」→「固定床高さ」の表示には反映されません。

そのため、上流端だけに固定床高さ 0 m を与えたケースでは、ほぼ平らな 0 m の固定床（窪み部分だけ地形高まで低い）の上に勾配のついた移動床が載った状態として計算されます。一部だけを固定床にしたい場合は、この仕様を踏まえて固定床高さを全域で意図した値に設定してください。

計算結果のエクスポートメニュー（VTK、CSV、tpo、シェープファイルなど）には CGNS の項目がありませんが、プロジェクトをフォルダ形式で保存すると CGNS ファイルが得られます。手順は次のとおりです。

1. 左上の「ファイル」メニューを開く
2. 「名前を付けてプロジェクトを保存」を選ぶ
3. 保存先として空のフォルダを選択する

保存したフォルダ内に Case1.cgns が出力されるので、これを EvaTRiP の入力ファイルとして読み込んでください。

現在の仕様では、粒度分布を与える領域は最大 10 個までで、それ以上に分割して設定することはできません。10 個を超える区分が必要な場合は、粒度分布の似た領域をまとめるなどして 10 個以内に収める必要があります。

現地の河床材料の粒度分布は、平面的にも深さ方向にも異なります。混合粒径で異常な洗掘が起きる場合、粗粒化しているべき場所に粗い材料が正しく与えられていない、あるいは実際には存在しないはずの細粒分が河床の深さ方向に与えられていて侵食が止まらない、といった可能性が考えられます。土質試験の結果をそのまま全域・全層に適用するのではなく、粒度分布の空間分布や深さ方向の層構成を見直し、条件を変えて試算してみてください。