中科大计算热物理讲义05对流扩散方程的数值方法.docxVIP

中科大计算热物理讲义05对流扩散方程的数值方法

对流扩散方程是描述流体中物质浓度随时间和空间变化的偏微分方程。在数值求解对流扩散方程时，主要方法包括有限差分法、有限元法和有限体积法等。以下是对流扩散方程数值方法的详细内容：

一、有限差分法

1.显式格式

对于一维对流扩散方程：

\[\frac{\partialu}{\partialt}+a\frac{\partialu}{\partialx}=D\frac{\partial^2u}{\partialx^2}\]

可以采用显式格式进行离散：

\[\frac{u_i^{n+1}u_i^n}{\Deltat}+a\frac{u_{i+1}^nu_{i1}^n}{2\Deltax}=D\frac{u_{i+1}^n2u_i^n+u_{i1}^n}{(\Deltax)^2}\]

其中，\(u_i^n\)表示第i个网格点在n时刻的浓度值。

2.隐式格式

对于一维对流扩散方程，隐式格式可以表示为：

\[\frac{u_i^{n+1}u_i^n}{\Deltat}+a\frac{u_{i+1}^{n+1}u_{i1}^{n+1}}{2\Deltax}=D\frac{u_{i+1}^{n+1}2u_i^{n+1}+u_{i1}^{n+1}}{(\Deltax)^2}\]

该格式需要求解一个三对角线性方程组。

二、有限元法

1.加权余量法

对于一维对流扩散方程，加权余量法的基本思想是寻找一组函数，使得方程的加权余量在所有可能的函数空间中为零。具体求解过程如下：

构造试探函数\(\phi_i(x)\)和权函数\(w_i(x)\)；

对方程两边乘以权函数\(w_i(x)\)并在定义域上积分；

利用分部积分和边界条件，得到加权余量方程；

采用Galerkin方法，选择试探函数与权函数相同，得到线性方程组。

2.变分法

变分法的基本思想是将偏微分方程转化为一个泛函极值问题。对于一维对流扩散方程，变分法求解过程如下：

构造泛函\(J(u)\)，使其包含原方程和边界条件；

对泛函\(J(u)\)进行变分，得到EulerLagrange方程；

求解EulerLagrange方程，得到原方程的近似解。

三、有限体积法

有限体积法是一种基于守恒原理的数值方法。对于一维对流扩散方程，有限体积法求解过程如下：

1.将求解区域划分为有限个小体积单元；

2.在每个小体积单元上，对原方程进行积分；

3.利用积分守恒原理，得到离散方程；

4.求解离散方程，得到原方程的近似解。

在上述三种数值方法中，有限差分法计算简单，但精度较低；有限元法精度较高，但计算复杂；有限体积法在守恒性方面具有优势，适用于复杂流动问题。实际应用中，应根据具体问题和计算需求选择合适的数值方法。