有限体积法（8）——混合差分格式

最新推荐文章于 2023-11-26 09:41:02 发布

薛定谔的alpha狗

最新推荐文章于 2023-11-26 09:41:02 发布

阅读量2.1k

点赞数 4

分类专栏： CFD 文章标签： CFD 计算流体力学

本文链接：https://blog.csdn.net/weixin_42562856/article/details/107258958

版权

本文深入探讨了混合差分格式，一种结合中心差分和迎风格式优点的数值方法，适用于不同Peclet数下的对流扩散方程离散。通过详细推导，展示了其在计算流体力学中的应用及优势。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

离散推导

Spalding（1972）提出了混合差分格式，该格式结合了中心差分格式和迎风格式的优点。小Pe数情况下（ $P e < 2$ ），使用中心差分格式，它具有二阶计算精度；大Pe数情况下（ $P e > 2$ ），使用迎风格式计算控制体界面对流输运量并忽略扩散作用。虽然迎风格式只有一阶精度，但可较好的反应流动的输运特征。
混合差分格整合了中心差分格式和迎风格式的计算公式，使用分段线性的计算公式来近似通过网格边界面处的通量。通过左边界单位面积通量的混合差分格式计算公式为
$\left. \begin{aligned} q_w &= F_w \left[ \frac{1}{2} \left( 1 + \frac{2}{Pe_w} \right) \phi_W +\frac{1}{2} \left( 1 - \frac{2}{Pe_w} \right) \phi_P \right] ，&-2<Pe_w<2\\ \\ q_w &=F_w \phi_W ，&Pe_w \ge2 \\ \\ q_w &= F_w \phi_P ，&Pe_w \le -2 \end{aligned} \right \}\tag{1}$
公式（1）中, $Pe_w$ 是当地边界面处的Peclet数，例如左边界面处定义为
$Pe_w=\frac{F_w}{D_w}=\frac{(\rho u)_w}{\Gamma_w/\delta x_{WP}} \tag{2}$
$q_w$ 指的是左边界处的总通量，包括对流和扩散，公式（1）第一个式子的推导如下，
$\begin{aligned} &F_e \phi_e - F_w \phi_w = D_e(\phi_E-\phi_P) - D_w(\phi_P -\phi_W) \\ \\ \Rightarrow & [F_e \phi_e - D_e(\phi_E-\phi_P)]-[F_w \phi_w-D_w(\phi_P-\phi_W)] = 0\\ \\ \Rightarrow & q_e - q_w = 0 \end{aligned} \tag{3}$
所以，在左边界处，当 $∣ P e ∣ < 2$ 时，对流项使用中心差分格式离散，即 $\phi_w=(\phi_W+\phi_P)/2$ 。则左边界处通量有，
$\begin{aligned} q_w &=F_w \phi_w -D_w(\phi_P-\phi_W) \\ \\ &=F_w \left( \frac{\phi_W+\phi_P}{2} \right)-D_w(\phi_P-\phi_W) \\ \\ &=\left( \frac{F_w}{2}+D_w \right)\phi_W + \left( \frac{F_w}{2}-D_w \right)\phi_P \\ \\ &=F_w \left[\frac{1}{2} \left( 1+\frac{2}{Pe_w} \right)\phi_W +\frac{1}{2} \left( 1-\frac{2}{Pe_w} \right)\phi_P \right] \end{aligned} \tag{4}$
公式（1）中的后两个式子是省略了扩散项后的通量，并且对流项使用了迎风格式
$\begin{aligned} q_w &= F_w\phi_w - D_w(\phi_P- \phi_W) \\ \\ &=F_w \phi_w\\ \\ &=F_w \phi_W ,当Pe_w \ge2 \\\\ &=F_w \phi_P ,当Pe_w \le-2 \end{aligned} \tag{5}$

由公式（3）可知，对流扩散方程可以写成
$q_e - q_w=0 \tag{6}$
那么把混合离散格式（1）带入到离散方程（6），则
当 $∣ P e ∣ < 2$ 时，
$\begin{aligned} q_e-q_w&=F_e \left[ \frac{1}{2} \left( 1+ \frac{2}{Pe_e} \right)\phi_P +\frac{1}{2} \left( 1- \frac{2}{Pe_e} \right) \phi_E \right] \\ \\ &\qquad -F_w \left[ \frac{1}{2} \left( 1+ \frac{2}{Pe_w} \right)\phi_W +\frac{1}{2} \left( 1- \frac{2}{Pe_w} \right) \phi_P \right] \\ \\ &=0 \end{aligned} \tag{7}$