L2补充完整推导

最新推荐文章于 2024-09-21 16:35:48 发布

乜乜学习日记

最新推荐文章于 2024-09-21 16:35:48 发布

阅读量272

点赞数

分类专栏： hydraulic

本文链接：https://blog.csdn.net/Wzzero3202/article/details/124175958

版权

算法

hydraulic 专栏收录该内容

8 篇文章 8 订阅

订阅专栏

j点为h的计算点，j-1、j+1为q的计算点

连续方程：

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

$\frac{\partial A}{\partial t} = B_{s}\frac{\partial h}{\partial t}$

$\frac{\partial h}{\partial t}\approx \frac{1}{\Delta t}\left ( h_{j}^{n+1}-h_{j}^{n} \right )$

$\frac{\partial Q}{\partial x}\approx \frac{1}{x_{j+1}-x_{j-1}}[\frac{1}{2}(Q_{j+1}^{n}+Q_{j+1}^{n+1}-\frac{1}{2}(Q_{j-1}^{n}+Q_{j-1}^{n+1})]$

$\frac{B_{s}}{\Delta t}(h_{j}^{n+1}-h_{j}^{n})+\frac{1}{x_{j+1}-x_{j-1}}\left [ \frac{1}{2}(Q_{j+1}^{n}+Q_{j+1}^{n+1})-\frac{1}{2}(Q_{j-1}^{n}+Q_{j-1}^{n+1}) \right ]=0$

$A_{j}^{*}Q_{j-1}^{n+1}+B_{j}^{*}h_{j1}^{n+1}+C_{j}^{*}Q_{j+1}^{n+1}=D_{j}^{*}$

$\left\{\begin{matrix} A_{j}^{*} =\frac{1}{4\Delta x}\\ B_{j}^{*} =\frac{B}{\Delta t}\\ C_{j}^{*} =-\frac{1}{4\Delta x}\\ D_{j}^{*} =-\frac{1}{4\Delta x}Q_{j+1}^{n}+\frac{B}{\Delta t}h_{j}^{n}+\frac{1}{4\Delta x}Q_{j-1}^{n} \end{matrix}\right.$

运动方程：

$\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+g\frac{\partial h}{\partial x}+g(S_{f}-S_{0})=0$

$\frac{\partial Q}{\partial t}+\frac{\partial }{\partial x}\left ( \frac{\partial aQ^{2}}{A} \right )+gA\frac{\partial h}{\partial x}+\frac{g\left | Q \right |Q}{C^{2}AR}=0$ (假设S0=0)

$a=\frac{A}{Q_{2}}\int_{A}^{}u^{2}dA$

$\frac{\partial }{\partial x}(\frac{\partial aQ^{2}}{A}) =\frac{2aQ}{A}\frac{\partial Q}{\partial x}-\frac{aQ^{2}}{A^{2}}(B+h\frac{\partial B}{\partial h})\frac{\partial h}{\partial x}+\frac{Q^{2}}{A}\frac{\partial a}{\partial h}\frac{\partial h}{\partial x}$

$\frac{\partial Q}{\partial x}=-B_{s}\frac{\partial h}{\partial t}$

$\frac{\partial }{\partial x}(\frac{\partial aQ^{2}}{A}) =-\frac{2aQB_{s}}{A}\frac{\partial h}{\partial t}-\frac{aQ^{2}}{A^{2}}(B+h\frac{\partial B}{\partial h})\frac{\partial h}{\partial x}+\frac{Q^{2}}{A}\frac{\partial a}{\partial h}\frac{\partial h}{\partial x}$

$\frac{\partial Q}{\partial t}-\frac{2aQB_{s}}{A}\frac{\partial h}{\partial t}-[\frac{aQ^{2}}{A^{2}}(B+h\frac{\partial B}{\partial h})-\frac{Q^{2}}{A}\frac{\partial a}{\partial h}-gA]\frac{\partial h}{\partial x}+\frac{g\left | Q \right |Q}{C^{2}AR}=0$

$\left\{\begin{matrix} D2=\frac{2aQB_s}{A}\\ D3 = \frac{aQ^{2}}{A^{2}}(B+h\frac{dB}{dh})-\frac{Q^2}{A}\frac{da}{dh}-gA\\ D4 = \frac{g}{C^{2}AR} =\frac{n^{2}g}{R^{4/3}A}\end{matrix}\right.$

$\frac{1}{\Delta t}(Q_{j-1}^{n+1}-Q_{j-1}^{n})-\frac{D2}{\Delta t}[\frac{x_{j}-x_{j-1}}{x_{j}-x_{j-2}}(h_{j-2}^{n+1}-h_{j-2}^{n})+\frac{x_{j-1}-x_{j-2}}{x_{j}-x_{j-2}}(h_{j}^{n+1}-h_{j}^{n})]-\frac{D3}{2(x_{j}-x_{j-2})}[(h_{j}^{n+1}+h_{j}^{n})-(h_{j-2}^{n+1}+h_{j-2}^{n})]+D4\left | Q_{j-1}^{n} \right |Q_{j-1}^{n+1}=0$

$A_{j-1}h_{j}^{n+1}+B_{j-1}Q_{j-1}^{n+1}+C_{j-1}h_{j-2}^{n+1}=D_{j-1}$

$\left\{\begin{matrix} A_{j-1}=-\frac{D2}{\Delta t}\frac{x_{j-1}-x_{j-2}}{x_{j}-x_{j-2}}-\frac{D3}{2(x_{j}-x_{j-2})}=-\frac{D2}{2\Delta t}-\frac{D3}{4\Delta x}\\ B_{j-1}=\frac{1}{\Delta t}+D4\left | Q_{j-1}^{n} \right |\\ C_{j-1}=-\frac{D2}{\Delta t}\frac{x_{j-1}-x_{j-2}}{x_{j}-x_{j-2}}+\frac{D3}{2(x_{j}-x_{j-2})}=-\frac{D2}{2\Delta t}+\frac{D3}{4\Delta x}\\ D_{j-1}=(-\frac{D2}{\Delta t}\frac{x_{j}-x_{j-1}}{x_{j}-x_{j-2}}+\frac{D3}{2(x_{j}-x_{j-2})})h_{j}^{n}+\frac{1}{\Delta t}Q_{j-1}^{n}+(-\frac{D2}{\Delta t}\frac{x_{j}-x_{j-1}}{x_{j}-x_{j-2}}-\frac{D3}{2(x_{j}-x_{j-2})})h_{j-2}^{n}\\ =(-\frac{D2}{2\Delta t}+\frac{D3}{4\Delta x})h_{j}^{n}+\frac{1}{\Delta t}Q_{j-1}^{n}+(-\frac{D2}{2\Delta t}-\frac{D3}{4\Delta x})h_{j-2}^{n} \end{matrix}\right.$