对于二维不可压流体,写出速度压力形式得连续性方程和运动方程;在如图所示的交错网格下,写出 p(1,j) 满足的差分格式
Latex
重写为
u_{i+1/2,j}^{n+1}=u_{i+1/2,j}^n-\frac{\Delta t}{4\Delta x}\left[u_{i+\frac{3}{2},j}^n\left(u_{i+\frac{3}{2},j}^n+2u_{i+\frac{1}{2},j}^n\right)-u_{i-\frac{1}{2},j}^n\left(u_{i-\frac{1}{2},j}^n+2u_{i+\frac{1}{2},j}^n\right)\right]-\frac{\Delta t}{4\Delta y}\left[\left(u_{i+\frac{1}{2},j}^n+u_{i+\frac{1}{2},j+1}^n\right)\left(v_{i,j+\frac{1}{2}}^n+v_{i+1,j+\frac{1}{2}}^n\right)-\left(u_{i+\frac{1}{2},j-1}^n+u_{i+\frac{1}{2},j}^n\right)\left(v_{i,j-\frac{1}{2}}^n+v_{i+1,j-\frac{1}{2}}^n\right)\right]-\frac{1}{\rho}\frac{\Delta t}{\Delta x}\left(p_{I+1,J}^n-p_{I,J}^n\right)+\frac{\mu\Delt