一维对流方程的差分格式及稳定性分析

一维对流方程的通用形式为：

$\frac{\partial u}{\partial t}+a\frac{\partial u}{\partial x}=0$

其中，$u(x,t)$为待求解的函数，$a$为常数，表示对流速度。

前差格式

$\frac{u_i^{n+1}-u_i^n}{\Delta t}+a\frac{u_i^n-u_{i-1}^n}{\Delta x}=0$

展开得：

$u_i^{n+1}=u_i^n-a\frac{\Delta t}{\Delta x}(u_i^n-u_{i-1}^n)$

稳定性分析

将$u_i^n=\lambda^n e^{ik\Delta x}$代入上式，得到：

$\lambda=e^{-ika\Delta t}\left(1+\frac{a\Delta t}{\Delta x}(e^{-ik\Delta x}-1)\right)$

当$\left|1+\frac{a\Delta t}{\Delta x}(e^{-ik\Delta x}-1)\right|\leq1$时，格式稳定。

后差格式

$\frac{u_i^{n+1}-u_i^n}{\Delta t}+a\frac{u_i^{n+1}-u_{i-1}^{n+1}}{\Delta x}=0$

展开得：

$u_i^{n+1}=u_i^n-a\frac{\Delta t}{\Delta x}(u_i^{n+1}-u_{i-1}^{n+1})$

稳定性分析

将$u_i^n=\lambda^n e^{ik\Delta x}$代入上式，得到：

$\lambda=\frac{1}{1+a\frac{\Delta t}{\Delta x}(e^{ik\Delta x}-1)}$

当$\left|\frac{1}{1+a\frac{\Delta t}{\Delta x}(e^{ik\Delta x}-1)}\right|\leq1$时，格式稳定。

中心差格式

$\frac{u_i^{n+1}-u_i^n}{\Delta t}+a\frac{u_{i+1/2}^n-u_{i-1/2}^n}{\Delta x}=0$

展开得：

$u_i^{n+1}=u_i^n-a\frac{\Delta t}{\Delta x}(u_{i+1/2}^n-u_{i-1/2}^n)$

其中，$u_{i+1/2}^n$和$u_{i-1/2}^n$需要通过插值得到。一种常用的插值方法是线性插值，即：

$u_{i+1/2}^n=\frac{1}{2}(u_i^n+u_{i+1}^n)$

$u_{i-1/2}^n=\frac{1}{2}(u_{i-1}^n+u_i^n)$

稳定性分析

将$u_i^n=\lambda^n e^{ik\Delta x}$代入上式，得到：

$\lambda=1-ika\frac{\Delta t}{\Delta x}\sin(k\Delta x)$

当$\left|1-ika\frac{\Delta t}{\Delta x}\sin(k\Delta x)\right|\leq1$时，格式稳定。

总结

综上所述，前差、后差和中心差格式的稳定性条件分别为：

• 前差：$\left|1+\frac{a\Delta t}{\Delta x}(e^{-ik\Delta x}-1)\right|\leq1$

• 后差：$\left|\frac{1}{1+a\frac{\Delta t}{\Delta x}(e^{ik\Delta x}-1)}\right|\leq1$

• 中心差：$\left|1-ika\frac{\Delta t}{\Delta x}\sin(k\Delta x)\right|\leq1$

其中，$\Delta t$和$\Delta x$分别表示时间步长和空间步长。稳定性条件的含义是，误差在计算中不会被放大，而是会随时间步长逐渐减小。因此，在选择差分格式时，需要保证其稳定性条件得到满足。