iRI
Cユーザ層拡大に向けてのエンジン及び周辺ツールの整備

北海道大学大学院 工学研究院

木村一郎

京都大学 防災研究所

竹林洋史

北海道大学大学院 工学研究院

岩崎理樹

iRI
Cユーザ層拡大に向けてのエンジン及び周辺ツールの整備

Dev
elopmentofsimulationengine
san
dre
late
dtool
stoe
xten
duse
rsofi
RIC.
木村 一郎1・竹林 洋史2・岩崎 理樹3
Ichi
roKI
MURA,Hirosh
iTAKEBAYASHIandToshi
kiIWASAKI
1

北海道大学 工学研究院 准教授
京都大学 防災研究所 准教授
3
北海道大学 工学研究院 博士研究員
2

要旨

国際的な河川数値解析共通プラットフォームとして認知されつつあるiRI
Cについて，更なるユーザー層の拡大を目指し，エンジン群，およびその周辺ツールの開発・整備を実施した．研究は３名の研究者によって共同研究として実施した．
まず，竹林は，Morpho1
DとMorpho2
DforKidsの開発を担当した．Morpho1
Dは，一次元の河床変動解析を行うための解析プログラムであり，新たに非定常流の計算コードを導入し，非定常・定常のどちらでも解析ができるように改良を行った．また，iRI
Cver.
2に対して，Morpho
1
Dの日本語及び英語のマニュアルを開発するとともに，条件設定のためのxml
ファイルの見直しを行った．さらに，小学校高学年以上の児童・生徒が河川の水や土砂の流れについて学ぶた
めのMo
rpho2
DforKids
作成のための問題点の検討を行うとともに，簡素化した
条件設定のためのxml
ファイルの新規開発，簡素化した条件設定に対応したソルバーの開発を行った．
一方，岩崎は，二次元ソルバーであるNays
2Dについて，さらなる機能の充実を実施した．Na
ys
2DはRICNays
時代からのソルバーで，歴史も古く成熟したソルバーと考えられるが，その分利用ユーザーも多く，その要望は多様化しつつある．これに応えるため，次のような３点について機能の強化を実施している．
・混合粒径浮遊砂モデルの導入
・二次流の発達・減衰を考慮した二次流モデルの導入
・プログラムの並列化による計算効率の向上
最後に，木村は，三次元ソルバーであるNays
CUBEについて，機能向上と改良を実施している．Nays
CUBEは，静水圧近似を用いない純然たる３次元モデルという意味ではiRICで唯一のものであり，また，このような三次元ソルバーを有するフリーソフトウエア自体が少ないことから，その使命は大きい．今回の機能強化のポイントは次の４点である．

7――1

・流木モデルの導入
CKスキームの導入
・QUI・非洗掘セル設定機能の追加
・アンダーパス設定機能の追加《キーワード：iRI
C；数値解析；フリーソフトウエア；河川学習ソフトウエア；流砂；河床変動；河川乱流》

8――1

目

次rpho
1
Dの開発とMo
rpho
2
DforKids作成の試み（担当：竹林洋史）
１． 非定常流れの解析を対象としたMo（１） Morpho
1
Dの非定常流れの基礎方程式と解析例a） Mo
rpho
1
Dの非定常流れの基礎方程式b） Mo
rpho
1
Dの非定常流れの解析例（２） Morpho
1
Dの日本語・英語マニュアルの開発
rpho
2
DforKids（３） Mo（４） 本章のまとめ
２． 平面二次元ソルバー Nays
2
Dの機能強化 （担当：岩崎理樹）（１） はじめに（２） 混合粒径浮遊砂モデルの導入a） 土砂連続式b） 掃流砂の計算c） 浮遊砂の計算（３） 二次流の発達・減衰を考慮した二次流モデルの導入a） 計算モデルb） 計算条件と計算ケースc） 計算結果d） 計算結果と線形安定解析の比較による数値計算モデルの評価e） まとめ（４） プログラムの並列化による計算効率の向上
a） 検討条件b） 結果と考察c） まとめ（５） 本章のまとめ
３． 三次元モデルの機能向上（担当：木村一郎）（１） 概要（２） 流木モデルの追加a） 概要b） モデルの基本コンセプトc） 一般曲線座標における球体群の並行移動，回転移動関係式の誘導と流木移動の計算手順d） Nays
CUBEへのモデルの組み込みe） 計算事例
（３） その他の機能向上a） QCUIKスキームの追加b） 河床が掘れないセルの設定c） アンダーパスの設定

9――1１．非定常流れの解析を対象としたMor
pho1Dの開発とMor
pho2DforKids作成の試み（担当：竹林洋史）（１）Mor
pho1Dの非定常流れの基礎方程式と解析例a）Morpho
1
Dの非定常流れの基礎方程式
4年度）は，定常流れの解析用にMo
rpho1Dを開発した．そこで，本年度（平成25年度）
前年度（平成2rpho1
Dのソースコードを修正した．以下，基礎方程式を示は非定常流れの解析も可能となるように，Moすが，Mo
rpho2
Dと同一の部分は割愛させて頂く．
流れの質量保存則および運動方程式は，以下のようである．

（1）
ここに，tは時間，xは流下方向の座標，Aは本川の流水断面積，Qは本川の流量，Qfは本川から氾濫原への横流出流量，Bfは氾濫原から本川への氾濫流の流入幅，gは重力加速度，Ie は本川のエネルギー勾配であり，以下のようにマニング則によって評価している．

（2）
ここに，n（＝A/S
mはマニングの粗度係数，Rww）は径深である．本解析では，雨期と乾期の流れと流砂に示すを評価するため，河道の横断形状を考慮している．そのため，流水川幅（Bw）と潤辺（Sw）は図1ように評価している．σ xxは乱流応力であり，以下のように評価している．

（3）（4）
は断面平均摩擦速度，hは断面平均水深である．
ここに，ɛ
は渦動粘性係数，κはカルマン定数，u*
急勾配河川でも対応できるように，流砂量は，河床勾配が限界掃流力に与える影響を考慮した芦田・道
上式を用いる2,3,4)．

図1 解析に用いた川幅（Bw ）と潤辺（Sw ）0――2

（5）とrについては後述する．sは土砂の水中比重，tはk
粒径階に対する無次元掃流力である．tはk粒fkkckbb**
径階に対する無次元限界掃流力であり，修正Egiazaroff式 により算定した2)．t*cmは掃流層内の土砂の平均粒径階に対する有効無次元掃流はk粒径に対する無次元限界掃流力であり，岩垣式5)により算定した．t*ek力である2)． uは*emは有効摩擦速度，t*mは掃流砂層内の土砂の平均粒径に対する無次元掃流力である．Kc
限界掃流力の修正関数3,4)であり，一次元場に修正した以下の関係より求める．

（6）
ここに，mcは静止摩擦係数，rは水の比重，qxは流下方向の河床勾配である．本解析は河道横断形状を考慮した一次元河床変動解析であるが，断面流砂量と断面内の流水域の河床材料の粒度は，それぞれ，一つの値としている．fはk粒径階に対する流水域の掃流砂層の土砂の存在率である．掃流砂層内の土kb
砂の粒度は，以下の土砂の質量保存則により求める．

（7）

図2 河床変動による河床横断面形状の変化
ここに，lは河床土砂の空隙率，zは河床位，fは掃流砂層直下の堆積層（遷移層4)）におけるk粒径階
bdakの土砂の存在率，Esdは流砂を全て河床に静止させたときの粘着性土層上の非粘着性土の厚さ，Ebは掃流砂層の厚さである．Edaは掃流砂層直下の堆積層（遷移層）の厚さである．遷移層以下のm番目の堆積層におけるk粒径階の存在率であるfは，河床が浸食され，遷移層に変化した後に，式（7）の関係で変化dmkする．cは鉛直平均の掃流砂の濃度，Ebeは平衡掃流砂層厚さである6)．fは安息角である．rは粘着性土層bb上の非粘着性土の量に関するパラメータである．浮遊砂の沈降速度（wfk）は，Rubeyの式より求めた7)．
nは動粘性係数，cは k粒径階の浮遊砂の基準点濃度である．cは基準点高さにおける k粒径階の平衡
sbksbek
浮遊砂濃度である8)．
k粒径階の断面平均浮遊砂濃度（c）は，以下の関係より求める．
sk（8）
浮遊砂の水深分布は，指数分布を仮定している．
本解析では，一般断面を対象とした解析が可能となっている．そのため，河床は図2に示すように，断1――2

面内の流水域のみ変化するとして以下の河床位方程式より求める．

（9）河床勾配が安息角以上となった場合，河床材料が非粘着性材料であれば，局所河床勾配が安息角以下となるように河床位を補正する（永瀬ら、1996）
表1 解析条件
流れの基礎式

堰越流状態

Case1

定常流

完全越流

Case2

定常流

潜り越流

Case3

非定常流

完全越流

Case4

非定常流

潜り越流

図3 給水流量の時間変化b）Morpho1Dの非定常流れの解析例
・解析条件
解析の対象とした河川は堰を有する河道である．堰付近については，堰の直上流域（1m上流）にも堰と同じ平面形状で解析格子を設定した．
解析条件を表1に示す．Case1は，流れの基礎方程式として非定常項を無視した定常場の基礎方程式を用いている．下流端水位条件及び上流端供給流量条件は，図3と図4に示している．非定常項を無視しても，各計算時刻において解析区間全体の流量を変化させるため，流量の時間的な変化を考慮した解析も可能である．Case2は，下流端の境界条件以外はCase1と同様である．Case2
では，下流端水位として162.
5mの一定値を与えており，堰地点における河床高さである161.
4mよりも高く設定し，堰地点において潜り越流が形成される条件となっている．Case3は，Case1と同一の条件であり，流れの基礎方程式として非定常場の基礎方程式を用いている．Case4は，Case2と同一の条件であり，流れの基礎方程式とし2――2Ỉ఩㸦m㸧て非定常場の基礎方程式を用いている．

᫬㛫㸦hr㸧

図4 下流端水位の時間変化（Case1,Case3）

図5 最深河床位と水位の縦断分布（Case1）・結果と考察
図5～図9にCa
se1～Case4の1時間，16時間（ピーク流量時），47時間後の水位と最深河床位の縦断分布を示す．図に示すように，全ての計算条件において，安定した流れ及び河床変動の解析結果が得られていることがわかる．また，Ca
se1の定常場の基礎方程式による解析結果とCa
se3の非定常場の基礎方程3――2式による解析結果を比較すると，水面形状の大まかな分布特性は両者で大きな違いは無いが，水面形状
se3の非定常場の基礎方程式による解析結果の方が滑らかになっていることが分かる．これは，がCa
Case3では，流量の空間的な変化によるスムージング効果などが原因と考えられる．また，河床位の計算結果を見ると，両者の違いは非常に小さいことが分かる．これは，河床変動の時間スケールは，流れの変動時間スケールに比べて非常に長いため，洪水流の伝播（流れの非定常性）の影響が河床の変化に現れにくいためである．また，図6及び図8に示すように，完全越流だけでなく潜り越流の場合の流れ及び河床変動の解析も可能である．これらの結果から，非定常の基礎方程式による解析は，流れの状態をより詳しく知りたい場合，計算区間が長く，上流端と下流端で流量が大きく異なる場合，津波遡上や天然ダムの決壊などの非定常性が非常に大きな場での解析での利用が有効と考えられる．

図6 最深河床位と水位の縦断分布（Case3）
図9にCase4の上流端から360mの断面の河床横断形状の時間的な変化を示す．図に示すように，本解析では一般断面形状を対象とした解析を行うとともに，横断方向の格子間隔が不等間隔でも対応可能である．図に示すように，上流端から360mの断面では，時間とともに河道中央部の河床のみが上昇し，河道横断形状が時間的に変化していることが分かる．これは，本解析では，実現象に即して，流水域のみで河床変動を許しているためである．

4――2

図7 最深河床位と水位の縦断分布（Case3）

図8 最深河床位と水位の縦断分布（Case4）5――2

図9 河床横断形状の時間的な変化（Case4,上流から360m断面）（２）Mor
pho1Dの日本語・英語マニュアルの開発
図10に，iRI
Cver
2用に開発したMorpho1
Dの日本語マニュアルを一部抜粋したものを示す．マニュア
ルは，Morpho2
D及びNays
CUBE等と同一フォーマットとなっており，概要，使用されている基礎方程式，
xml
ファイルによる条件設定用ダイアログ，入力データフォーマットの説明の順で記載されている．
Morpho1
Dの特徴としては，横断形状として一般断面形状の利用が可能，混合砂の取り扱い，岩露出域や河川構造物上の流砂量の算定などである．これらについては，基礎方程式で記述するとともに，入力データフォーマットの部分で必要な入力データ形式を詳述している．
（３）Mor
pho2DforKids
防災・減災対策において，インフラ整備と同様に重要となるのが防災教育である．災害は老若男女等しく襲ってくるものであるから，防災教育についても，なるべく低年齢から施していくことが望ましいといえる．子供たちに対する防災教育としては，書物やビデオ教材等が考えられるが，これらはインタラクティブ性を有しない受動的教材である．一方，小学生高学年以上に対しては，これらの受動的教材に加えて，インタラクティブ性のある教材を与え，自ら手を動かしながら学ぶことが望まれる．河川防災教育についてもこの例にもれず，洪水，土石流，豪雨等について，自ら条件を設定し，その結果を能動的に理解する教材の登場が待たれる．i
RI
Cは，これまで専門家の手のものであった河川の数値解析技術の裾野を広げることにある程度成功したが，用語等の制約から，その対象は高校生程度以上と考えるのが適当である．したがって，さらに低年齢の子供たちに対してインタラクティブな教材を提供するには，既往のiRICを低年齢向きに改良した別バージョンで対応せざるを得ないと考えられる．このような目的のための最初の試みとして，Morpho2
Dをベースに，Morpho2
Dforkidsの開発を実施した．
図11に，iRI
Cver2用に開発したMo
rpho2
DforKidsの条件設定画面と通常のMo
rpho2
Dの条件設定画面の比較を示す．従来のMo
rpho2
D条件設定画面と比べて，Morpho2
DforKidsの条件設定画面の大きな違い

6――2

は，専門用語をわかりやすい言葉で表現している点と条件設定項目の量が減っている点である．今回は，小学校高学年以上の児童・生徒を対象としているため，簡単な漢字は用いている．小学校低学年も対象とした場合は，全てひらがなでの表記となるが，その場合は小学校高学年・中学生・高校生にとっては読みづらくなると思われ，小学校低学年用のトランスレーションファイルは別途用意する方が良いと思われる．
条件設定の簡素化及び基礎的な現象の勉強・理解のため，河床材料は一様粒径のみを対象とした．また，植生や障害物の設定条件は残したが，固定床の設定ウィンドウは削除し，移動床のみの解析を対象とした．さらに，計算途中の状況をコンソールに表示する時間間隔などは，タイムステップの1000倍に固定したり，流れのタイムステップに対する河床変動のタイムステップの比を1で固定したりするなど，容易に条件設定ができ，シンプルな条件で解析ができることを優先した．

7――2

図10 i
RI
Cv
er
2用に開発したMor
pho1Dの日本語マニュアル8――2

図10 i
RI
Cv
er
2用に開発したMor
pho1Dの日本語マニュアル（続き）

9――2

図1
1 Mor
pho2DforKidsの条件設定画面と通常のMor
pho2Dの条件設定画面の比較

0――3

図1
1 Mor
pho2DforKidsの条件設定画面と通常のMor
pho2Dの条件設定画面の比較（４）本章のまとめ
Morpho1
Dに非定常流れの計算コードを導入し，非定常・定常のどちらでも解析ができるように改良を行った．また，iRI
Cver2に対して，Morpho
1
Dの日本語及び英語のマニュアルを開発するとともに，条件設定のためのxml
ファイルの見直しを行った．さらに，小学校高学年以上の児童・生徒が河川の水や土砂の流れについて学ぶためのMo
rpho2
DforKids
作成のための問題点の検討を行うとともに，簡素化した
条件設定のためのxml
ファイルの新規開発，簡素化した条件設定に対応したソルバーの開発を行った．
参考文献（第１章分）
1） 竹林洋史，藤田正治：粘着性・非粘着性土混在場における一般断面による一次元河床変動解析，砂防学会誌，第64巻 第2号，pp.3
14，2011．
2） 芦田和男，道上正規:移動床流れの抵抗と掃流砂量に関する基礎的研究,土木学会論文報告集，第206号，pp.
5969，1972.
3） Kovacs,A.andPa
rker,G.
：Anewv
ector
ialbedloadformulationanditsa
pplicationtothetimeev
olutionofstraigh
triv
erch
annels,J.Fl
uidMech.Vol.2
67,pp.1
53183，1994.
4） B.Y.Liu：StudyonSedimentTranspor
tandBedEvolu
tioni
nCo
mpoun
dChannels.
,京都大学学位論文，
1991．

1――3） 岩垣雄一：限界掃流力の流体力学的研究，土木学会論文集，第41号，pp.
121，1956.
5） Egashira,S.an
dAshida,K.Uni
fiedviewofthemechanicsofdebrisflowan
dbedload,Adv
ancesin
6Micromechanicso
fGranularMaterials,(EditedbyH.
H.
Shenetal.
)Elsev
ier,pp.3
91400,1992.
7） Rubey,W.W.
：Settlin
gvelocitieso
fg
ravel,sandandsiltparticles,Amer
icanJ.ofSc
ience,Vol.25,pp.325338，1933．
） Lane,E.W.andKalinske,A.A.
：Engineerin
gcalcu
lati
o
nofsu
s
pendedsedi
ment,” Trans.A.
G.
U.
,Vol.22，81941．
9） 永瀬恭一，道上正規，檜谷治：狭窄部を持つ山地河川の河床変動計算，水工学論文集，第40巻，pp.
887892，1996.

2――3