几种简单差分格式

最新推荐文章于 2021-11-25 22:12:45 发布

lk040384

最新推荐文章于 2021-11-25 22:12:45 发布

阅读量1.1w

点赞数 8

分类专栏：程序数学

本文链接：https://blog.csdn.net/lk040384/article/details/105929032

版权

程序同时被 2 个专栏收录

4 篇文章 0 订阅

订阅专栏

数学

3 篇文章 0 订阅

订阅专栏

以对流方程为例：

一阶迎风

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{u_{i+1}^{n}-u_{i}^{n}}{\Delta x}$ （时间向前空间向前）

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{u_{i}^{n}-u_{i-1}^{n}}{\Delta x}$ （时间向前空间向后）

二阶迎风

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{-u_{i+2}^{n}+4u_{i+1}^{n}-3u_{i}^{n}}{2\Delta x}$

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{3u_{i}^{n}-4u_{i-1}^{n}+u_{i-2}^{n}}{2\Delta x}$

三阶迎风

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{-u_{i+2}^{n}+6u_{i+1}^{n}-3u_{i}^{n}-2u_{i-1}^{n}}{6\Delta x}$

$\frac{u_{i}^{n+1}-u_{i}^{n}}{\Delta t}=-a\frac{2u_{i+1}^{n}+3u_{i}^{n}-6u_{i-1}^{n}+u_{i-2}^{n}}{6\Delta x}$

二阶中心

$\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}$ 取成平均值，就得到了Lax-Friedrichs格式

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

在时间和空间上都用二阶中心差分离散，可以得到蛙跳（leap-frog）格式

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

将 $u_{i}^{n+1}$ 在 $u_{i}^{n}$ 进行Taylor级数展开

$u_{i}^{n+1}=u_{i}^{n}+\frac{\partial u}{\partial t}\Delta t+\frac{\partial ^{2}u}{\partial ^{2} t}\frac{(\Delta t)^2}{2!}+O(\Delta t)^3$

根据对流方程有

$\frac{\partial u}{\partial t}=-a\frac{\partial u}{\partial x}$

$\frac{\partial ^{2}u}{\partial ^{2} t}=-a\frac{\partial}{\partial t}(\frac{\partial u}{\partial x})=a^2\frac{\partial ^{2}u}{\partial ^{2} x}$

代入Taylor展开式得

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

用中心差分离散后，可得到Lax-Wendroff格式

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

lk040384

关注

8
点赞
踩
33

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录