4.三次元ソルバー Nays
CUBEの機能強化，及び定常・非定常一次元ソルバー Nays
1Dの開発
（担当：木村一郎）（1）概要
今年度4月にiRICver.
2.
0のリリースを終え，デフォルトでインストールされるソルバーの数も着実に増となった（Nays
2D，Nays2
DFlood，Morpho
2
D，Nays
CUBE，De
lft
3D，River2
D，
えつつあり，現時点で8Fa
STMECK，Storm）．しかし，このうち，流れを三次元的に取り扱うものはNays
CUBEとDe
lft
3Dの2つに過ぎず，鉛直方向の加速度を考慮した純然たる意味での三次元モデルはNays
CUBEの一つだけとなる．
それだけに，現時点でのNays
CUBEの使命は大きく，三次元解析を必要とするユーザーに対して必要な機能を確実かつ安定に供することが強く求められる．そこで，今年度はVer.
2.
0
リリース前後から寄せられたユーザーからの要望をもとに，ソルバーの改良を実施し，機能強化を図った．主な機能強化の内容
は，これまで鉛直二次元ソルバー Nays
CLI
Pで装備されていた浮遊砂解析機能の追加，周期境界条件設置時の流量調整機能の強化の2点である．
一方，河川行政担当者一次元モデルの装備に関する要望が大きいことを考慮し，一次元ソルバーであるNa
ys
1Dを新たに開発した．一次元ソルバーとしては，1.
で述べたMorpho
1
Dも今年度平行して開発が進められているが，Nays
1DとMo
rpho
1
Dの大きな違いは，後者が定常モデル（不等流モデル）であるのに対し，前者は定常モデル，非定常モデルのいずれも装備している点にある．ただし，Morpho1Dは一般断面形状が扱えるに対し，Nays
1
Dは現時点では矩形断面形状のみに対応している．
以下，Nays
CUBE，Nays
1
Dのそれぞれのソルバー開発状況について，節を改めて述べていく．
（2）Nays
CUBEの機能強化（1）・・・浮遊砂解析機能の追加a）背景と基本方針
浮遊砂が河床変動に及ぼす影響は，一般に掃流砂に比べて相対的に小さいとみなされる場合が多いが，水理条件によっては浮遊砂が卓越するケースももちろん存在する．また，河川による物質輸送に着目する場合，浮遊砂の影響が掃流砂に比べて格段に多いケースも報告されている．これらの点を考慮し，これ ま で 鉛 直 二 次 元 ソ ル バ ー Nays
CLI
Pで 装 備 さ れ て い た も の と 同 等 の 浮 遊 砂 輸 送 解 析 モ デ ル を
Naysu
CUBEに組み込むこととした．組み込む浮遊砂解析モデルは河川で最も一般的に用いられるオイラー形の移流拡散方程式をベースに，底面付近の基準面濃度に板倉・岸のモデルを用いたモデルを採用した．
b）浮遊砂の基礎式
浮遊砂輸送方程式は，ソース項の付加されたスカラー移流拡散方程式であり，移動一般曲線座標のもとでは次のようにあらわされる．
（1）
ここに，c：浮遊砂濃度（体積濃度），Dt：浮遊砂拡散係数，wf：浮遊砂沈降速度をそれぞれ表す．右辺は拡散項を表しており，内部セル（底面セル以外のセル）では，次のように分解される．

（2）

7――9（2）式の右辺第３項は浮遊砂ソース項（浮遊項，シンク項は移流項
一方，最下層の底面セルについては，を考慮して，次のように分解される．
中に含まれる）qus（3）

ここに，添え字upperはセル上面の値，lower
はセル下面（底面セルの場合河床面）の値を示す．ソース項（浮遊項）は平衡状態を考慮し，基準面濃度Caと沈降速度wfの積として次のように表わす．
（4）
ここに，基準面濃度は基準面高さ（本モデルでは基準面高さはz=a=Δz|k=1/2：最下層の高さの1/2，としている．）における浮遊砂濃度であり，このモデル化においては板倉・岸モデル（1980）1）を用い，次のように表記される．
（5）
ここに，（6）と表記される．一方，B*，h0，a*，Kはモデル定数であり，これらの値については，アインシュタイン.
143，h0=0.
5，a*=0.
14，K=0.
008と与える．
（1950）にならい，B*=0
砂粒子の沈降速度 wfについては，Rubeyの式により，次のように計算する．
（7）
6 2m/s)を表す．
ここに， S：砂粒子の水中比重（≒1.
65），n：水の動粘性係数 (10②河床の連続式
一般座標系における河床の連続式（エクスナー方程式）を示す．
・デカルト座標系（8）xy
ここに，λ:河床砂の空隙率，qx方向の掃流砂量フラックス成分,qy方向の掃流砂量フラックス成分
B:B:をそれぞれ表わす．qsu:
浮遊砂ソース項，wfcb:
浮遊砂シンク項をそれぞれ表わす．
・一般曲線座標系（9）8――9xh
二次元ヤコビアン,qx方向の掃流砂量フラックスの反変成分,qh方向の掃流砂量フラック
ここに，J:
B:
B:

スの反変成分をそれぞれ表わす.
c）浮遊砂輸送モデルに関するユーザーインターフェースの改良
浮遊砂輸送モデルの導入にあたり，計算条件設定画面のうち，「河床に関する条件の設定」ウインドウ
4.
1ように変更した．
について，図-図4.
1 河床に関する条件設定ダイアログ
（画面）

図4.
1の各設定項目については，表4.
1に示すとおりである．また，計算結果の出力項目については，「浮遊砂濃度」の出力を追加した．

9――9

4.
1 河床に関する条件設定の説明
表#項目

1 平均粒径(m)
（注2）

設定方法

注意事項

河床材料の平均粒径を設定します．単位 単位は(m)に注意．

2 考慮する流砂の 「掃流砂のみ」，「浮遊砂のみ」，「掃流砂と種類浮遊砂両方考慮」のいずれかから選択します．
3 掃流砂量モデル 掃流砂のモデルを，Meyer
Peter
Muller モデルの詳細はⅡ.3を参照．通常は，芦田道の選択（注1） モデル，芦田道上モデル，Kov
acsParker 上モデルをベースに主流方向の河床勾配を考モデルの中から選択します．
慮したKovacsP
arker
モデルを推奨します．
3
34 河床材料の密度 河床材料の密度を[kg/m]
単位で数値で 通常の砂礫の場合は2650 [kg/m]
指定します．

5 河床の空隙率河床の空隙率を設定します．

6 静止摩擦係数

河床材料の静止摩擦係数を設定します． デフォルトは0.
7.

7 動摩擦係数

河床材料の動摩擦係数を設定します．

8 河床水中安息角 河床の水中安息角を設定します．
（度）

デフォルトは0.
4.
通常は変更の必要はありません．
デフォルトは0.
5.
単位は(deg).
デフォルトは35°．この値を大きくとりすぎると計算が不安定になりやすくなります．

9 河 床 変 動 係 数 河 床 変 動 係 数（Morphologicalfactor） 本来河床変動係数は1ですが，流れの変化に比（注4）の値を数値で指定します．河床変動係数 べて河床変動の変化が遅い場合，この値を1よは，1時間ステップ当たり河床変動量に掛 り大きくとることによって計算を効率化できかる倍率です．
ます．
（但し，大きくとりすぎると計算が不安定になることがあるので注意が必要です．）
10 浮遊砂の流入条 上流端からの浮遊砂の流入方法を，「平衡 通常は「平衡状態を仮定したラウス分布」を件（注3）状態を仮定したラウス分布」「浮遊砂濃 選択します．上流端の浮遊砂濃度の計測デー，度を定数で与える」のいずれかから選択 タが得られている場合，あるいは上流端におします．
ける浮遊砂濃度が平衡状態と異なる場合には「浮遊砂濃度を定数で与える」を選択します．
1
1 流入浮遊砂濃度 一つ上の項目で「浮遊砂濃度を定数で与 浮遊砂濃度は体積割合で指定します．たとえ3の河川水中に粒子分（固形分）のえる」を選択した場合，ここでその値を ば，100cm
3存在する場合は，0.
001と指定し体積が0.
1cm指定します．
ます．

注1）掃流砂量モデルの選択
本バージョンで選択可能な掃流砂モデルは，いずれも平衡流砂モデルです．
注2）平均粒径の設定
本バージョンでは均一粒径のモデルを採用しています．混合粒径モデルについては次期バージョン以降で対応予定です．
注3）上流端の浮遊砂濃度分布
この設定は，Ⅲ.
7「境界条件」で主流方向周期境界条件を設定した場合には無視されます．
「浮遊砂濃度を定数で与える」を選択した場合，上流端鉛直方向の濃度分布は一様で与えられます．
注4）河床変動係数の影響
河床が動的な場合（たとえば砂堆の移動を伴う場合など），河床変動は河床の移動速度にも影響を及ぼすことに注意が必要です．すなわち，計算結果における河床の移動速度（砂堆の進行速度など）は物理的な速度に河床変動係数を乗じた値となります．

00――1

図4.
2 周期境界条件流量調整機能追加に伴うインターフェースの変更（3）Nays
CUBEの機能強化（2）・・・周期境界条件の流量調整機能の強化
a）背景と基本方針
主流方向周期境界条件は，流れ場が主流方向に長く一様で，現象が周期性を有する場合に計算を効率化する上で極めて有効な手法である．主流方向周期境界条件を設定すると，下流端からの流出をそのまま上流側の流入とするため，流量を陽に設定することができない．したがって，底面や壁面などの摩擦応力の評価が実問題と異なると，流量が実問題と異なるという問題が生じる．この問題を回避するため，
旧ソルバーでは，重力（すなわち勾配）の影響を加減することで流量を設定流量に動的に一致させるという方法を採用してきた．この手法では，初期の流体体積が保存させるため，平均水深は計算を通じて不変である．しかし，動的に流量が変化する洪水流などの場合には，重力項を恣意的に操作することは物理現象と乖離している．
この問題を回避するため，重力項は変化することなく，水面の上下により動的に流量を設定流量に一致させる新たな流量調整機能を追加することとした．
b）ユーザーインターフェースの変更
上記の機能追加に伴い，「計算条件設定」の「境界条件」設定ウインドウについて，図4.
2に示すように変更を加えた．具体的には，ウインドウの4つ目の設定項目として「周期境界条件設定時流量調整方法」を追加し，新たに追加した「水深で調整」という手法と，旧来の「勾配で調整」という手法をユーザーが選択可能となるよう変更した．
（4）定常・非定常一次元ソルバー Nays
1Dの開発a）開発動機と基本方針
河川実務の現場，特に行政担当者から一次元河川シミュレーション機能の搭載に関する要望がよせられていることを契機に，一次元ソルバーの開発を行うこととした．一次元モデルは解くべき方程式の数は２次元，あるいは３次元ソルバーに比べて少ないが，境界条件や地形形状の扱い等にといて，独自の工夫が必要である．本年度は将来の本格的一次元ソルバー開発の基盤となるシンプルなプロトタイプの01――1

開発を実施した．すなわち，河川横断形状としては矩形形状とし，ユーザーは各KPについて，川幅と平均水深の２つの情報のみを用意することとした．また，河床変動計算については，現時点では均一粒径の掃流砂輸送のみを考慮することとした．一方，流れ場の解析方法としては，定常一次元ソルバーであるMo
rpho
1Dとの差別化を図ることをも考慮して，定常モデル（不等流モデル）に加え，非定常モデル（不定流モデル）も選択可能とした．両者のメリット・デメリットは次の通りである．
【不等流計算】メリット：
・流れに関する時間ステップΔtの制限が無い．したがって計算時間は短い．長期間の計算（数十年以上）が可能
デメリット：
・流れの時間変化の到達が場所的に異なる現象を再現できないため，現象によっては精度が落ちる．
・流れに関するΔt制限は無いが，河床変動が急激に進む場で大きいΔtをとるとは非合理な河床形状が生じ，異常終了する場合がある．
・射流区間では通常の計算方法が使えないので，近似が必要．
【不定流計算】メリット：
・流れの時空間変化を考慮できる．洪水到達の場所的遅れやヒステレシスも再現できる．流れの変化に伴う河床変動の進行についても時間的なラグを生じずに評価できる．
デメリット：
・計算時間間隔Δtに制限があるため，計算所要時間は不等流に比べて大きい．
b）基礎方程式
・定常モデル（不等流モデル）の流れの基礎方程式
定常モデルの流れの基礎方程式は，定常連続式とエネルギー式からなり，次のように記述される．
［連続式］
（10）
［エネルギー式］（11）
ここに，Q:流量，x:流れ方向座標，A:断面積，g:重力加速度，ρ:水の密度，S:潤辺，τb:摩擦応力，z:水位である．摩擦応力を含む左辺第3項については，マニングの粗度係数nを用いて，次のように表現sする．
（12）・非定常モデル（不定流モデル）の流れの基礎方程式
非定常モデルの流れの基礎式は，
［連続式］
（13）

02――1［運動方程式］（14）
ここに，v:主流方向流速，t:時間である．左辺第4項の摩擦応力を含む項は次のようにマニングの粗度係数を用いて次のように表す．
（15）
ここに，Qを陽に掛けた形にしたのは，離散化の際に半陰解法的手法を用いて計算の安定化を図るためである．
c）河床変動モデル
河床変動モデルとしては，掃流砂輸送のみを考慮し，掃流砂モデルとしては，一般的な芦田・道上モ
デル2）と，芦田・道上モデル2）をベースに主流方向の河床勾配の影響を考慮したコバックス・パーカーモデル3）の2つの中から一つをユーザーが選択することとした．これらのモデルの具体的表記を下に示す．
・芦田・道上モデル（1972）：（16）・山口・泉モデル（2003）：（17）
ここに， z河床高，q
流線方向の掃流砂量，K:補正関数，s:流線方向の座標，mc:
動摩擦係数，
τ無次元b:
Bs:掃流力，τ*c:無次元限界掃流力をそれぞれ表している．河床変動については，一般的なエクスナー方程式を用いて計算する．
d）入力データNays
1
Dに必要な入力データは，①河道地形データ②流量データ（ハイドログラフ）③河床平均粒径④平均マニング粗度係数の4つである．このうち，地形データ，流量データ（ハイドログラフ）のデータフォーマットは図4.
3，図4.
4にそれぞれ示す．
上記の①，②のデータについては，図のフォーマットに従ってユーザーがテキスト形式（csv形式）でデータをあらかじめ用意する．③，④については，平均値のみを用意し，計算条件設定画面上から入力する．したがって，外部データを用意する必要はない．なお，流量データに関しては，iRIC上でハイドログラフデータを手入力するオプションも用意した．

03――1

図4.
3 Nays1Dの地形データフォーマット

図4.
4 流量ハイドログラフのフォーマット

図4.
5「河道条件」の設定ウインドウ

04――1

図4.
6「流量と水位の設定」のウインドウ

図4.
7「時間に関する条件」の設定ウインドウ

図4.
8 「その他の条件」の設定ウインドウ

05――1

図4.
9 Nays1D簡易マニュアルの作成

06――1

図4.
9 Nays1D簡易マニュアルの作成（続き）e）計算条件設定ウインドウの構成
計算条件設定ウインドウは①「河道条件」，②「流量と水位の設定」，③「時間に関する条件」，④「その他の設定」の4画面からなる．これらを図4.
5～図4.
8に示す．
なお，異なる計算方法である「不等流モデル」と「不定流モデル」の選択は，図4.
5に示すように「河道条件」の設定ウインドウ内で行うこととした．また，これらのウインドウから，①，②の外部データのリンク先についても指定することとした．
「その他の条件」としては，流砂モデルとして芦田・道上モデル，コバックス・パーカーモデルのいずれかから一つを選択する．
f）利用の手引き資料の作成
Nays
1
Dについては，マニュアル，チュートリアルについては，今年度中の作成には至らなかった．これは，ある程度講習会等でソルバーを使用し，その結果をフィードバックしソルバーを改良することを考慮したためである．しかしながら，ユーザーに対してソルバーのコンセプト，目的，外部データ，計算設定，可視化等の実際を一通り簡明に示した簡易マニュアルを兼ねたプレゼンテーションファイルを作成した．この内容を図4.
9に示す．
（5）まとめと今後の課題
iRI
Cエンジンのうち，三次元ソルバー Na
ys
CUBEに関する2つの機能向上と，一次元ソルバー Na
ys1
Dの新たな開発について，2012
年度の成果を述べてきた．Nays
CUBEについては，今回の浮遊砂機能の追加により，ソルバーの安定性や定常状態での機能についてはほぼ完成したと見て良いが，流量の変化が大きい場合について，今後は次の二点について改良を行っていきたい．
・計算時間間隔の動的設定
・水際移動機能の追加
一方，一次元ソルバー Nays
1
Dについては，本年度は本格的ソルバーのベースとなるプロトタイプの開発に重点を置いた．今後はさらに機能を高めるとともに（浮遊砂の考慮，混合粒径モデルの追加，一般断面の取り扱いを可能にする等）マニュアル，事例集の整備を実施し，Morpho
1
Dとの差別化を明確にしつつもさらに実用的なモデルとしていきたい．

07――1参考文献） It
aku
ra,T.
,andKishi,T.
:Openchannelflowwi
thsuspendedsediments,Proc.Am.Soc.Civ.En
g.
,Vol.
106,
1No.
8,pp.
13251343,1980.
2） 芦田和男,道上正規:移動床流れの抵抗と掃流砂量に関する基礎的研究,土木学会論文集,
第206号,
pp.
5969,1972.
） 山口里実，泉典洋：デューン-平坦床遷移過程にみられる亜臨界分岐現象，土木学会論文集，740号，
3pp.
7594，2003.

08――1