写真:常流をまっすぐ上流へ向かうカワウソ
1次元定常不等流計算を学ぶ際、誰もが最初に習う「鉄則」があります。
- 常流(Subcritical Flow):下流 → 上流 へ計算する
- 射流(Supercritical Flow):上流 ← 下流 へ計算する
実務プログラムでも当然のようにこの方向に固定されていますが、なぜこの計算方向でなければならないのでしょうか?
「エネルギー方程式を解くだけなら、逆から解いても数学的には同じ解になるのでは?」という疑問を持つ方もいるかもしれません。
実際、背水影響を含んだ正しい水深が分かっている状態で、Excel等を用いて逆方向から逐次計算を行っても、両者の計算結果は完全に一致します。
では、なぜ水理学では「常流は下流から、射流は上流から」と計算方向が定められているのか。その理由を、「水深の確定性(出発水位)」と「境界条件(コントロール)」の考え方から分かりやすく解説します。
計算の方向を決定づける最大の理由は、水面波(微小波)の伝播速度と流速の関係、そしてそれに伴う「出発水位(起点となる水深)を確定できるか否か」という点にあります。
不等流計算を正しく進めるためには、計算のスタート地点となる出発水位があらかじめ「既知(確定)」でなければなりません。
常流では流速よりも波の伝播速度の方が速いため、下流側の構造物(ダムや堰)や水位変化の影響(波)が上流へ伝わり、上流側の水位に反映されるのが実現象です。
したがって、下流端の水位(H.W.L.や堰の天端水深など)は物理的な拘束条件として明確に既知とすることができます。しかし、その下流の影響を受けた後の「上流側の水位」をあらかじめ知ることは困難です。
確定できない上流を起点にすることはできないため、「水深が確定している下流から上流へ」向かって計算を進めることになります。
射流では流速が波速を上回るため、下流で何が起きてもその情報は上流へ遡ることができません。上流から流れてくる水の勢いや流入条件がそのまま下流側の水位へ反映されます。
そのため、上流端の条件(水門の開きや急勾配の始点など)は既知として確定できますが、下流側の水位をあらかじめ知ることは困難です。
結果として、「水深が確定している上流から下流へ」向かって計算を進めることになります。
ここで冒頭の疑問に戻ります。
「常流において、下流からの背水影響を受けた『正解の上流水位』が判明しているなら、上流から下流へ計算しても正解になるのではないか?」
結論から言えば、その通りであり、同じ解が得られます。
なぜなら、2断面間の不等流計算は「ベルヌーイの定理(エネルギー保存則)」に沿って計算しているだけだからです。
上流と下流のあいだで成立しているエネルギーの関係式(水頭+損失の釣り合い)を解いているだけなので、すでに正解の出発水位が分かっているなら、上流から下流に向かって計算しても数学的にまったく同じ解に辿り着きます。実際にExcelなどで逆方向から計算しても完全一致するのはこのためです。
では、なぜ世界的な河川解析プログラム(HEC-RASなど)や標準的な手引きにおいて、計算方向が厳格に定められているのでしょうか。
それは「逆方向には計算できない」という意味ではなく、水面形を決定づける『境界条件(コントロール)』の側から計算するのが最も合理的だからです。
波の伝播速度との関係から、常流では下流側の影響(波)が上流へ伝わって水位に反映され、射流では逆に下流側へ反映されるのが実現象である。 したがって、計算の起点となる「出発水位」はあらかじめ既知(確定)である必要があるが、常流において下流の影響を受けた上流側の水位をあらかじめ既知とすることは困難である。射流においても同様に、下流側の水位を既知とすることは困難である。 そのため、常流は水深が確定できる下流から上流へ向けて計算し、射流は上流から下流へ向けて計算するのである。
不等流計算における計算方向のルールは、決して「逆から解くと数学的に計算が破綻する」といった極端な話ではありません。
- 水深の確定性: 流れを支配する境界条件(常流は下流、射流は上流)でしか、最初の出発水位を確定できないため。
- 解の再現性: すでに確定した水面形であれば、逆方向に解を進めても同じエネルギー関係から正しい水深が得られる。
- 計算の合理性: 結果が分からない未知の区間を解く以上、流れを決定づけている原因(支配断面)から順に解いていくのが最も自然で安定したアプローチであるため。
「計算方向の指定は、単なる数式の都合ではなく、物理的な境界条件から水面形を組み立てるための合理的な手順である」と理解しておくと、不等流計算の仕組みがよりすっきりと見えてくるはずです。
