趣味で計算流砂水理

Computational Sediment Hydraulics for Fun

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

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

基礎式です。


\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}

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

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