有限体积法（9）——高阶差分格式：QUICK格式

最新推荐文章于 2023-10-12 18:32:21 发布

薛定谔的alpha狗

最新推荐文章于 2023-10-12 18:32:21 发布

阅读量1w

点赞数 9

分类专栏： CFD 文章标签： cfd

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

版权

迎风格式和混合格式只有一阶计算精度，虽然迎风格式使用起来非常稳定并且满足输运性要求，但一阶精度容易导致数值扩散的误差，可以通过使用高阶离散格式来降低这些误差。高阶格式一般需要使用更多的节点值，通过考虑更大范围内的影响来降低误差。中心差分格式具有二阶精度，但它的稳定性很差，也不满足输运特性。考虑到向中心差分格式这种不能包含流动方向信息的格式是不稳定的，我们就需要寻找那些包含流动方向信息的高阶差分格式。下面就介绍这类差分格式：QUICK格式。

方程离散

QUICK(Quadratic upstream interpolation for convective kinetics)格式是由英国学者Leonard（1979）提出的用于计算控制体界面值的二次插值计算格式。QUICK格式利用控制体界面两侧的三个点进行插值，包括两个位于界面两侧的相邻节点和一个位于上游侧的远邻点，如下图。
在这里插入图片描述
当 $u_w>0,u_e>0$ 时，通过节点WW，W和P的拟合曲线用于计算控制体左侧界面的参数 $\phi_w$ ，而通过节点W、P和E的拟合曲线用于计算控制体右侧界面的参数 $\phi_e$ 。当 $u_w<0,u_e<0$ 时，则用节点W、P和E来计算左侧界面的 $\phi_w$ ，而用节点P、E和EE来计算右侧界面的 $\phi_e$ 。若用下标 $i - 1$ 和 $i$ 分别代表界面左、右两侧的相邻节点，用下标 $i - 2$ 代表上游的远邻点，那么QUICK格式的公式为：
$\phi_{face}=\frac{6}{8}\phi_{i-1}+\frac{3}{8}\phi_i-\frac{1}{8}\phi_{i-2} \tag{1}$
当 $u_w>0$ 时， $w$ 界面左右两侧相邻节点为W和P，上游远邻点为WW，如上图。该界面处参数插值的近似公式为
$\phi_w=\frac{6}{8}\phi_W+\frac{3}{8}\phi_P-\frac{1}{8}\phi_{WW} \tag{2}$
当 $u_e>0$ 时， $e$ 界面左右两侧相邻节点为P和E，上游远邻点为W，则该界面处参数插值的近似公式为
$\phi_e=\frac{6}{8}\phi_P+\frac{3}{8}\phi_E-\frac{1}{8}\phi_{W} \tag{3}$
理论上，扩散项的梯度也可以使用二次拟合曲线的导数来近似，但是对于均匀网格的情况，过两点的抛物线在其中点处的导数和过这两个点的直线的斜率相等，所以这里扩散项依旧使用中心差分格式来离散。那么对于前面的一维对流扩散方程，如果 $F_w>0,F_e>0$ ，我们使用上述的QUICK格式来离散对流项，用中心差分格式来离散扩散项，得到其离散方程
$\begin{aligned} \left[ F_e\left(\frac{6}{8}\phi_P+\frac{3}{8}\phi_E-\frac{1}{8}\phi_{W} \right) -F_w\left(\frac{6}{8}\phi_W+\frac{3}{8}\phi_P-\frac{1}{8}\phi_{WW} \right) \right]\\\\ =D_e(\phi_E-\phi_P) -D_w(\phi_P-\phi_W) \end{aligned} \tag{4}$
整理之，
$\begin{aligned} \left[ D_w-\frac{3}{8}F_w +D_e+\frac{6}{8}F_e \right]\phi_P=\left[ D_w+\frac{6}{8}F_w+\frac{1}{8}F_e \right]\phi_W \\\\+ \left[ D_e-\frac{3}{8}F_e \right] \phi_E - \frac{1}{8}F_w\phi_{WW} \end{aligned} \tag{5}$
写出标准形式，
$a_P\phi_P = a_W\phi_W +a_E\phi_E+a_{WW}\phi_{WW} \tag{6}$
系数：

最低0.47元/天解锁文章