趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning

数値計算とか河川工学とかプログラミングのことを書いています

2017-09-12

続：改めて不等間隔格子の河川流一次元計算

冷静に考えたら、座標変換ですね。

基礎式です。

$\frac{\partial A}{\partial t} + \frac{\partial \xi}{\partial x}\frac{\partial Q}{\partial \xi} = 0$

$\frac{\partial Q}{\partial t} + \frac{\partial \xi}{\partial x}\frac{\partial }{\partial \xi} \left( \frac{Q^2}{A} \right) = -gA\frac{\partial \xi}{\partial x}\frac{\partial H}{\partial \xi} -gAI_e$

$\frac{\partial \xi}{\partial x} = \frac{1}{\Delta x_{i+1} + \Delta x_{i}}\left(\frac{ \Delta x_{i}}{ \Delta x_{i+1} }(\xi_{i+1} - \xi_{i} ) + \frac{ \Delta x_{i+1}}{ \Delta x_{i} }(\xi_{i} - \xi_{i-1} )\right)$

$\Delta x_{i+1} = x_{i+1} - x_{i}$

$\Delta x_{i} = x_{i} - x_{i-1}$

実河川の問題では必須のような気がします。

今度テスト計算をしてみます。