有限体积法（7）——迎风格式

薛定谔的alpha狗

已于 2022-06-19 12:36:45 修改

阅读量1.6w

点赞数 10

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

于 2020-07-10 17:51:09 首次发布

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

版权

方程离散

中心差分格式主要的一个不足之处就是它不能辨识流动的方向。在中心差分格式中，左边界的输运量 $\phi$ 的值总是由 $\phi_P$ 和 $\phi_W$ 共同决定。然而，强对流（从左向右流动）的流场中这种处理方式是不合理的，因为左边界处的 $\phi$ 值受到上游节点W的影响要远大于下游节点E的影响。为了克服中心差分格式不具有对流输运特性的缺点，人们提出了迎风格式。迎风格式在确定边界值时会考虑到流动方向的影响，边界处的 $\phi$ 值就直接用上游节点的值来近似，比如在左边界处就是 $\phi_w=\phi_W$ （从左向右流动时）， $\phi_w=\phi_W$ （从右向左流动时），如下所示。
在这里插入图片描述
考虑一维对流扩散问题的离散中推导出的公式 $(9)$ ，
$F_e \phi_e - F_w \phi_w = D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W) \tag{1}$
定义向右流动为正向，向左流动为负向。则当流动为正向时有，
$u_w>0, u_e>0 \Rightarrow F_w>0,F_e>0 \tag{2}$
根据迎风格式，边界处的 $\phi$ 值为
$\phi_w=\phi_W,\phi_e=\phi_P \tag{3}$
带入到方程式 $(1)$ 中，有
$F_e \phi_P - F_w \phi_W = D_e(\phi_E-\phi_P)-D_w(\phi_P-\phi_W) \tag{4}$
重新排列一下，有
$(D_w+D_e+F_e) \phi_P=(D_w+F_w) \phi_W + D_e \phi_E \tag{5}$