有限体积法（5）——对流-扩散方程的离散

薛定谔的alpha狗

已于 2022-06-18 12:11:42 修改

阅读量1w

点赞数 14

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

于 2020-07-09 18:08:48 首次发布

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

版权

方程离散

关于变量 $\phi$ 的输运方程，
$\frac{\partial (\rho \phi)}{\partial t}+ \nabla \cdot (\rho \phi \bold u) = \nabla \cdot (\Gamma \nabla \phi) + S_\phi \tag{1}$

省略时间项就是稳态的对流扩散方程，
$\nabla \cdot (\rho \phi \bold u) = \nabla \cdot (\Gamma \nabla \phi) + S_\phi \tag{2}$

在控制体上积分，并根据奥-高公式，有
$\begin{aligned} \int_{CV} \nabla \cdot (\rho \phi \bold u) dV&= \int_{CV} \nabla \cdot (\Gamma \nabla \phi) dV + \int_{CV} S_\phi dV \\ \\ \Rightarrow \int_{A} \bold n \cdot (\rho \phi \bold u) dA &= \int_{A} \bold n \cdot (\Gamma \nabla \phi) dA + \int_{CV} S_\phi dV \tag{3} \end{aligned}$
考虑一维稳态无源的情况，则对流扩散方程为
$\frac{d}{dx} (\rho u \phi) = \frac{d}{dx} (\Gamma \frac{d \phi}{dx}) \tag{4}$

流动同时也满足连续性方程
$\frac{d}{dx} (\rho u ) = 0 \tag{5}$
对流扩散方程在一维空间上的离散与扩散方程类似，公式 $(1)$ 和 $(4)$ 的对流项和扩散项都有散度，扩散项还有梯度，那么其离散也是包括散度离散和梯度离散。一维空间的网格如下，假设网格是均匀网格，网格间距的表示如图，速度 $u$ 的方向是从左往右的
在这里插入图片描述
散度的离散即把 $\rho \phi \bold u$ 的散度在控制体单元内的体积分转换为 $\rho \phi \bold u$ 在边界上的面积分，
$(\rho u A \phi)_e-(\rho u A \phi)_w=\left( \Gamma A \frac{d \phi}{dx} \right)_e - \left( \Gamma A \frac{d \phi}{dx} \right)_w \tag{6}$

扩散项中梯度的离散用中心差分格式，则对流扩散方程离散为
$(\rho u A \phi)_e-(\rho u A \phi)_w=(\Gamma A)_e \frac{\phi_E-\phi_P}{\delta x_{PE}} - (\Gamma A)_w \frac{\phi_P-\phi_W}{\delta x_{WP}}$
为了表示简单，定义，
$F=\rho u, D=\frac{\Gamma}{\delta x} \tag{7}$
则各边界上，
$\begin{aligned} F_w=(\rho u)_w \qquad F_e=(\rho u)_e \\ \\ D_w=\frac{\Gamma_w}{\delta x_{WP}} \qquad D_e=\frac{\Gamma_e}{\delta x_{PE}} \tag{8} \end{aligned}$
假设各边界面的面积相等，即 $A_w=A_e=A$ ，则离散方程为
$F_e \phi_e - F_w \phi_w = D_e(\phi_E-\phi_P) - D_w(\phi_P-\phi_W) \tag{9}$

最低0.47元/天解锁文章