一维对流方程的蛙跳格式的截断误差是什么,怎么算的
时间: 2023-05-25 16:05:44 浏览: 166
二维对流扩散方程恒稳的蛙跳积分格式 (2002年)
一维对流方程的蛙跳格式的截断误差为:
$\tau_i^n = \frac{\partial u}{\partial t} - \frac{u_i^{n+1} - u_i^n}{\Delta t} + a\frac{u_{i+1}^n - u_{i-1}^n}{2\Delta x}$
其中,$u_i^n$是时间$t_n$和空间$x_i$处的数值解,$\Delta t$和$\Delta x$分别是时间和空间的离散步长,$a$是对流速度。
蛙跳格式是一种二阶精度的时间离散方法,所以它的截断误差也是二阶的。对于一维对流方程:
$\frac{\partial u}{\partial t} + a\frac{\partial u}{\partial x} = 0$
应用蛙跳格式,我们有:
$\frac{u_i^{n+1} - u_i^n}{\Delta t} + a\frac{u_{i+1}^n - u_{i-1}^n}{2\Delta x} = O(\Delta t^2) + O(\Delta x^2)$
将这个式子代入截断误差的定义式中,我们可以得到:
$\tau_i^n = O(\Delta t^2) + O(\Delta x^2) - a\frac{\partial u}{\partial x} + O(\Delta t \Delta x)$
由于蛙跳格式是一个二阶格式,所以$\Delta t^2$和$\Delta x^2$项的系数是$O(1)$,即忽略不计的。由此,我们可以得到,一维对流方程的蛙跳格式的截断误差为:
$\tau_i^n = O(\Delta t \Delta x) - a\frac{\partial u}{\partial x}$
其中,$O(\Delta t \Delta x)$是截断误差的主要部分,而$a\frac{\partial u}{\partial x}$是精确解中的一项,因此,截断误差的主要来源是时间和空间离散化的不精确。
阅读全文