よくある質問

上流端の流速分布を「上流端水深から逆算」にしていると、上流端に植生が設定されている場合などに流入の分布が不自然になることがあります。次の対処を試してください。

• 計算条件の上流端の流速分布を「等流計算」に変更する（これで低水路側から流れるようになった事例があります）
• 上流端付近のセルの植生密度を 0 にする、または一度植生なしで計算してみる
• それでも改善しない場合は、高水敷の上流端に障害物セルを設定する

類似の質問が過去のフォーラムにもありますので、あわせて参照してください。

残念ながら、Nays2DH には残り計算時間を表示する機能はありません。計算の進行状況は、ソルバーコンソールに表示される計算時刻や、結果の出力タイムステップ数から見積もってください。

格子のインデックス I（流下方向）が逆になっていることが原因です。Nays2DH を含む多くのソルバーでは I=1 が上流端で、そこから水が流入します。下流端が I=1 になっている格子では、下流端から水が流入する結果になります。

横断測量データ（*.riv など）をインポートする際、測線の並び順が反転していると、「横断測量データから生成」で作られる格子も上下流が反転します。インポート時のダイアログで測線の並び順（上流側から下流側へ / 下流側から上流側へ）を指定できるので、入れ替えて再度格子を作成してください。格子作成時に「開始横断線」が上流側になっていれば正常です。

なお、断面名がすべて数値であれば iRIC は数値の大きさで順序を自動認識しますが、「15+16.0」のように数値以外が混ざった断面名があると自動認識できません。並び順の指定が煩わしい場合は、断面名をすべて数値（例: 15.8）にしてください。

Nays2DH が計算している掃流力（河床せん断力）は、計算結果に出力される ShearStress(Nm-2) です。無次元掃流力は、この値を τ₀ として τ₀ / (s·g·d) で求めます。ここで s は土粒子の水中比重で、計算条件の「+その他」で指定している値を使うのが整合的です。g は重力加速度、d は粒径（m）です。Nays2DH には無次元掃流力を直接出力する機能はないため、結果をエクスポートして表計算ソフトなどで計算するか、オープンソースのソルバーに出力機能を追加する方法があります。なお、流速が 0 の箇所で掃流力が 0 にならない理由については、このスレッドでは明確な回答が得られていません。

規模にもよりますが、砂漣や砂堆といった小規模河床波が計算で再現されている可能性があります。また、計算格子やタイムステップなど河床変動計算に関わる条件設定によっては、河床が不連続なスパイク状になるケースもあります。Nays2DH には河床変動を抑制する負のフィードバック処理は実装されていないため、格子サイズや条件設定側で対応する必要があります。河床波の理論的背景は「現場のための水理学（令和版）」の河床波の章や iRIC Academy で扱っています。

Nays2DH の河床変動計算では、固定床か移動床かの設定が重要です。移動床に設定していると、掃流力の条件に応じて計算中は常に河床変動が考慮されます。「計算条件」→「時間」にある「河床変動の計算開始時刻（秒）」を使えば、計算の初期〜中期は河床変動なしで進め、途中から河床変動を有効にすることは可能です。しかし、最初に河床変動を起こしてから途中で止めるという設定は現状ありません。定常流（一定流量）で河床変動が起こる流量を与え続けると河床は安定しませんが、時間ごとに流量が変わる非定常のハイドログラフを与えれば、流量が小さい時間帯では河床変動は収まります。一度変動した河床が元に戻ることはありません。どうしても河床が安定せず解析が進まない場合は、いったん「固定床」で計算してみるのも一つの方法です。

格子CSVの列構成は、格子生成ツールではなく使用するソルバーに依存し、オブジェクトブラウザーの「格子」以下に表示される属性がそのまま出力されます。ファイルの構造は次のとおりです。

• 1〜2行目: I, J, K 方向の格子点数
• 3行目: ヘッダ行
• 4行目以降: インデックスごとの値

各行の1〜3列目は I, J, K のインデックス（iRIC上は1始まりですが、CSVでは0始まり）、4〜6列目（X, Y, Z）はその格子点の座標です。7列目以降が属性値で、列名の先頭文字が要素の種類を示します（N_ は格子点、C_ はセル。Nays2DHにはありませんが I_ / J_ はエッジ）。Nays2DHの場合の主な列は次のとおりです。

• N_Elevation: 地形高 (m)
• N_Elevation_zb: 固定床高さ (m)
• C_Obstacle: 障害物セル（0: 通常、1: 障害物）
• C_Fix_movable: 固定床・移動床（0: 移動床、1: 固定床）
• C_vege_density: 植生密生度 (1/m)
• C_vege_height: 植生高さ (m)
• C_roughness_cell: マニングの粗度係数
• C_mix_cell: 河床材料の粒度分布

セルは各方向とも格子点より1つ少ないため、セル属性の I, J, K 最大インデックスの行にはダミー値 0 が入っています。なお、Excel で編集して保存すると表示桁数で値が丸められることがあるため、十分な桁数を表示した状態で作業してください。詳細は iRIC ユーザーマニュアルの「iRIC 格子CSVファイル」および Nays2DH ソルバーマニュアルの格子属性の項を参照してください。

セル (i, j) は、その周囲の4つの格子点 (i, j)、(i+1, j)、(i, j+1)、(i+1, j+1) で構成されます。つまりセルのインデックスは、そのセルの左下（I, J とも小さい側）の格子点のインデックスと同じです。このため、セルの数は I, J 方向とも格子点数より1つ少なくなります。Python などで格子属性（C_ で始まるセル属性）を外部で作成・編集する際は、この対応関係に沿ってセル値を配置してください。

隣り合う格子のサイズが急に変化している箇所があると、計算が発散して停止することがあります。横断測量データから格子を生成する機能は、横断データの河道中心線を格子の中心線として使うため、中心線と左右岸までの距離に差がある場所では左右の格子幅の比が大きくなりがちです。対策としては次の方法があります。

• 横断測量データからの格子生成を使い続ける場合: 横断データの河道中心線を「河道の中心」ではなく「格子範囲の中心」を通るように修正してから格子を生成すると、左右の格子幅のバランスが改善します。
• 他の格子生成ツールを使ってよい場合: 汎用格子生成ツールは格子範囲を自由に変更できるので、堤内地を含めず低水路内だけに格子を作ることができます。「ポアソン方程式を解いて格子を生成」を使うと、比較的簡単に整った格子を作れる場合があります。

隣接格子のサイズ比を小さくする指針は iRIC の「iRICを操作する上での注意点」資料にも記載されています。

現状、障害物セルの範囲だけ計算結果を描画しないという専用の設定はありません。対応方法としては、可視化の設定で表示範囲をしきい値で制限する（たとえば流速 0.1 m/s 以下の領域を描画しない）といった方法で、障害物範囲の表示を目立たなくすることができます。

「Calculation is failure!」は計算が発散したときに表示されます。まず試すべき対策は、計算条件の時間刻み（dt）を小さくすることです。iRIC 公式サイトの FAQ（計算に関する項目）にも同様の対処が記載されています。また、格子間隔が細かすぎる場合には、格子間隔を広く取り直すことで計算が通るようになった例があります。dt の縮小と格子間隔の見直しを組み合わせて調整してみてください。

エラー行そのものではなく、その手前の断面で #x-section の後に宣言した座標点数と、実際に並べた点の数が一致していない可能性が高いです。たとえば実際には147点しか記述していないのに 34 215 のように宣言していると、読み込み側は215点分を読み進めようとし、次の断面のデータを点データとして解釈してしまうため、後続の行で「インデックスの値が不正です」となります。Excel などで *.riv を生成している場合は、各断面の #x-section 行の点数と実際の行数が一致しているかを確認・修正してください。実際にこの点数の修正だけで正常に読み込めるようになった例があります。ファイル形式の詳細は iRIC ユーザーマニュアルの「河川測量データ」の項を参照してください。