有限体积法求解二维方腔流（二）——边界条件处理及SIMPLE和PISO算法

关于动量方程和连续性方程的离散参考这里

1.5. 边界处理

1.5.1. 左边界，$u_w=v_w=0$，不包括边角

• 对流项
  $$\begin{gathered} \sum F_f^n{U}_f^{n+1}=(-u_w^{n}h)U_w^{n+1}+(u_e^{n}h)U_e^{n+1}+(-v_s^{n}h)U_s^{n+1}+(v_n^{n}h)U_n^{n+1} \\ =\frac{h}{4}\left[(u_P^n+u_E^n+v_N^n-v_S^n)U_P^{n+1}+(u_P^n+u_E^n)U_E^{n+1}-(v_P^n+v_S^n)U_S^{n+1}+(v_P^n+v_N^n)U_N^{n+1}\right] \end{gathered}$$

• 扩散项
  $$\begin{gathered} \sum \nu (\nabla U^{n+1}{)}_{f}S=\nu \left[(\nabla U^{n+1}{)}_w{S}_w+(\nabla U^{n+1}{)}_e{S}_e+(\nabla U^{n+1}{)}_s{S}_s+(\nabla U^{n+1}{)}_n{S}_n\right] \\ =\frac{\nu S_f}{d}\left[U_E+U_S+U_N-5U_P\right] \end{gathered}$$

• 压力项
  $$\sum p_{f}S=p_P\left(-h,0\right)+\frac{p_E+p_P}{2}\left(h,0\right)+\frac{p_S+p_P}{2}\left(0,-h\right)+\frac{p_N+p_P}{2}\left(0,h\right)=\frac{h}{2}\left(p_E-p_P,p_N-p_S\right)$$

• 压力泊松方程-左端
  $$\sum (HbyA{)}_f\cdot S=h(HbyA{|}_{e,x}+HbyA{|}_{n,y}-HbyA{|}_{s,y})$$

• 压力泊松方程-右端
  $$C_w=0$$

• 系数
  $$A_P=\frac{V_P}{\Delta t}+\frac{h}{4}(u_P+u_E+v_N-v_S)+\frac{5\nu S_f}{d},A_W=0$$

1.5.2. 右边界，$u_e=v_e=0$，不包括边角

• 对流项
  $$\begin{gathered} \sum F_f^n{U}_f^{n+1}=(-u_w^{n}h)U_w^{n+1}+(u_e^{n}h)U_e^{n+1}+(-v_s^{n}h)U_s^{n+1}+(v_n^{n}h)U_n^{n+1} \\ =\frac{h}{4}\left[(-u_W^n-u_P^n+v_N^n-v_S^n)U_P^{n+1}-(u_P^n+u_W^n)U_W^{n+1}-(v_P^n+v_S^n)U_S^{n+1}+(v_P^n+v_N^n)U_N^{n+1}\right] \end{gathered}$$

• 扩散项
  $$\begin{gathered} \sum \nu (\nabla U^{n+1}{)}_{f}S=\nu \left[(\nabla U^{n+1}{)}_w{S}_w+(\nabla U^{n+1}{)}_e{S}_e+(\nabla U^{n+1}{)}_s{S}_s+(\nabla U^{n+1}{)}_n{S}_n\right] \\ =\frac{\nu S_f}{d}\left[U_W+U_S+U_N-5U_P\right] \end{gathered}$$

• 压力项
  ∑ p f S = p W + p P 2 ( − h , 0 ) + p P ( h , 0 ) + p S + p P 2 ( 0 , − h ) + p N + p P 2 ( 0 , h ) = h 2 ( p P − p W ,   p N − p S ) \sum p_fS=\frac{p_W+p_P}{2}\left( -h,0 \right)+p_P\left( h,0 \right)+\frac{p_S+p_P}{2}\left( 0,-h \right)+\frac{p_N+p_P}{2}\left( 0,h \right) =\frac{h}{2}\left( p_P-p_W,\ p_N-p_S \right) ∑pf​S=2pW​+pP​​(−h,0)+pP​(h,0)+2pS​+pP​​(0,−h)+2p