理论教育 一阶波动方程差分格式:土建类数值解法

一阶波动方程差分格式:土建类数值解法

时间:2023-10-26 理论教育 版权反馈
【摘要】:一阶波动方程典型的差分格式,除了Lax格式以外,还有偏心格式和LaxWendroff格式等,以下作简要介绍。综合上述分析,差分格式和式是否稳定与微分方程式的系数α的符号有关。故差分格式和式分别称为“右偏心”格式和“左偏心”格式。

一阶波动方程差分格式:土建类数值解法

用两组平行直线x=xi=ilt=tn=分割成矩形网格,其中,lτ分别为空间x和时间t的步长,交点(xitn)称为结点,对t=tn上的全部结点称为差分网格的第n层。

一阶波动方程典型的差分格式,除了Lax格式(1.5.29)以外,还有偏心格式和LaxWendroff格式等,以下作简要介绍。

1.偏心格式

如果对时间步长τ和空间步长x均采用前差商,则对应于方程式(1.1.17)

978-7-111-44528-9-Chapter01-230.jpg

的差分方程写为

uni+1=(1+runi-runi+1 (1.7.1)

其中,978-7-111-44528-9-Chapter01-231.jpg

若对时间采用前差商,而对空间采用后差商,则得到

uni+1=(1-runi+runi-1 (1.7.2)

关于显式差分式(1.7.1)和式(1.7.2)的稳定性,应用稳定性的Fourier分析方法,可得出,当α>0时,显式差分式(1.7.1)不稳定;当α<0,-1≤r≤0时,式(1.7.1)稳定。而对于格式(1.7.2)在α<0时不稳定;当α>0,0≤r≤1时,格式(1.7.2)稳定。根据CFL条件也可得出,当978-7-111-44528-9-Chapter01-232.jpg时,两种格式都是稳定的。

综合上述分析,差分格式(1.7.1)和式(1.7.2)是否稳定与微分方程式(1.1.17)的系数α的符号有关。我们知道,方程式(1.1.17)有一特征线

x-αt=CC为常数)

α>0时,特征线向右倾斜,当α<0时,特征线向左倾斜,因此,所构造的差分格式是否稳定实质上与特征线的走向有关。差分格式(1.7.1)和式(1.7.2)的图分别如图1.7.1和图1.7.2所示。

由稳定性分析可知,若所构造的差分格式与微分方程的特征线的走向一致,则在r满足一定条件下是稳定的;反之,不论r如何都是不稳定的。事实上,当α>0时,过点(in+1)的特征线与t=ti的交点定落在点(in)的左侧;当α<0时,过点(in+1)的特征线与t=ti的交点定落在点(in)的右侧。故差分格式(1.7.1)和式(1.7.2)分别称为“右偏心”格式和“左偏心”格式。

978-7-111-44528-9-Chapter01-233.jpg

图1.7.1 时间前差空间前差(www.daowen.com)

978-7-111-44528-9-Chapter01-234.jpg

图1.7.2 时间前差空间后差

2.LaxWendroff格式(LW格式)

LW格式为

978-7-111-44528-9-Chapter01-235.jpg

为了构造格式(1.7.3),利用对时间的Taylor展开,可得

978-7-111-44528-9-Chapter01-236.jpg

利用微分方程(1.1.17),得到

978-7-111-44528-9-Chapter01-237.jpg

于是,有

978-7-111-44528-9-Chapter01-238.jpg

再利用中心差分公式

978-7-111-44528-9-Chapter01-239.jpg

略去高阶误差项,则得到我们所要构造的差分格式(1.7.3),其截断误差是Oτ2+h2)。

应用稳定性的Fourier分析方法,可得出,当|r|≤1,即978-7-111-44528-9-Chapter01-240.jpg时,LW格式满足CFL条件,其差分问题的解是稳定的。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈