第四章有限差分方法离散方程

第四章有限差分方法离散方程第一页，共四十九页，编辑于2023年，星期五1基本概念在实际问题中，我们所关心的是因变量在空间若干特定位置的数值。

将因变量在给定点的数值直接作为未知数系数，并求解这些数值，作为满足实际需要的解。第二页，共四十九页，编辑于2023年，星期五离散方法比较项目有限差分法有限元法普方法试函数程序难易程度程序灵活性精确性计算效率适宜的方程主要优点主要缺点局部近似很好好差好各类型经济、程序简单较难推广到高阶局部近似好很好好好椭圆型灵活性好不经济总体近似差差很好很好椭圆型精度高不灵活第三页，共四十九页，编辑于2023年，星期五有限差分的概念在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下：（1）区域离散化（2）近似替代（3）逼近求解第四页，共四十九页，编辑于2023年，星期五2差分的基本形式及精度用差商代替微分方程中的导数。以空间导数为例，微分中心为（n,j）向前差分向后差分中心差分第五页，共四十九页，编辑于2023年，星期五

上述几种差分形式可通过Taylor级

数展开的方法，得到前差分和后差分具

有一阶精度;中心差分具有二阶精度。第六页，共四十九页，编辑于2023年，星期五3显式差分与隐式差分显式格式以时间步差——空间中心差对上式离散则因为nΔt时刻的Γ值为已知，可用上式直接计算n+1时刻的Γ值。扩散方程第七页，共四十九页，编辑于2023年，星期五隐式格式

以时间步差——空间在（n+1）层中心差时对扩散方程离散。

未知函数不可能通过上式由已知值直接求解，它必须求解线性方程组才能求出。此为隐式格式。第八页，共四十九页，编辑于2023年，星期五4有限差分格式的

相容性、收敛性及稳定性概念：从偏微分方程建立差分方程时，总是要求τ→0，h→0时差分方程与微分方程充分接近。

作用：研究差分方程与微分方程的关系。

分类：有条件相容和无条件相容有限差分格式的相容性第九页，共四十九页，编辑于2023年，星期五以扩散方程为例：当时间步差——空间中心差得差分方程为若以Taylor级数展开扩散方得：第十页，共四十九页，编辑于2023年，星期五类似于导数的差分形式的截断误差，扩散方程的差分形式的截断误差为o(Δt,Δx2)。如果

Δx，Δt→0时，截断误差o(Δt,Δx2)→0，则称差分方程与原微分方程是相的。当Δx，Δt以任何形式→0时，o(Δt,Δx2)→0，则称无条件相容。当Δx，Δt以某种方式→0时，o(Δt,Δx2)→0，则称有条件相容。第十一页，共四十九页，编辑于2023年，星期五有限差分格式的收敛性概念：指差分方程的解，当Δx，Δt→0时是否逼近原始微分方程的真解。作用：研究差分方程的解是否逼近真解的问题。第十二页，共四十九页，编辑于2023年，星期五有限差分格式的稳定性

概念：指差分方程在求解的过程中，差分方程的解能否保持一致有界。

作用：差分方程的稳定性是其收敛性的充分必要条件，它具有实用价值。

分类：点稳定和步稳定。

第十三页，共四十九页，编辑于2023年，星期五

为了理解稳定性的概念，下面介绍两种类型的不稳定。对流扩散方程

用FTCS离散

在n时刻方程有一个稳定解，由于某种原因存在一个扰动，由该扰动带来解的误差,

假定其为线性叠加即待人上式第十四页，共四十九页，编辑于2023年，星期五则有

对流项扩散项其显式解：将方程变形可得第十五页，共四十九页，编辑于2023年，星期五为了便于讨论，将上述两项的影响分开来讨论

由图可见，1)由于εj+1n>0,εjn<0,εj-1n>0,则4εj+1n>0；由于εjn<0，所以，|εj+1n|<|εjn|扩散项：第十六页，共四十九页，编辑于2023年，星期五即Δεj+1n趋向于校正负的扰动εjn，同理可分析出Δεj+1n+1<0，即Δεj+1n+1正好校正正的扰动εj+1n。εj+1n的幅度小于εjn，扰动会趋向于消失，所以扩散过程有利于计算的稳定。2)εj+1n与Δt有关，与Δt成正比，若Δt很大，随着Δt增大，Δ

|εj+1n|增大；若|εj+1n|<|εjn|，则εj+1n会形成振幅不断增大的振动型过冲，有可能不稳定，属于动态不稳定，可用减小Δt的办法来消除。第十七页，共四十九页，编辑于2023年，星期五对流项：

假定对流速度u＞0，扰动是一个振荡型的。|εjn|，则对j节点有即Δεj+1n与εjn，

|εj+1n|<|εjn|。从而扰动随时间不断的单调增大。结果是不稳定的。称为静态不稳定。第十八页，共四十九页，编辑于2023年，星期五它不能依靠改变参数来消除，只有改变差分差分格式才能避免。

在实际计算中，这种初始误差的产生和分布常常是随机的，若处理不当，会造成计算不稳定。如果方程中对流项与扩散项同时存在时，两者的相互牵制对时间步长的限制条件，取决于对流项与扩散项的相对重要程度。第十九页，共四十九页，编辑于2023年，星期五

概念：法国数学家傅里叶发现，任何周期函数都可以用正弦函数和余弦函数构成的无穷级数来表示（选择正弦函数与余弦函数作为基函数是因为它们是正交的），后世称为傅里叶级数，是一种特殊的三角级数。公式如下：其中kj:波数λ=2π/kj相角φ=kjΔjωj：周期

傅里叶级数第二十页，共四十九页，编辑于2023年，星期五

设，上式可整理为扩散方程为：用FTCS离散扩散方程用VonNeumann法对扩散方程（FTCS）格式进行稳定分析第二十一页，共四十九页，编辑于2023年，星期五方程的解用傅里叶分量可写成其中：kx为波数；λ为波长，λ=2π/kjx，当λ→

∞，kx→0.所以kx=0代表直线；定义相角θ=kxΔx,是波数为kx的分量在时刻n的幅度函数。将解的傅里叶分量带人差分方程第二十二页，共四十九页，编辑于2023年，星期五消去得，利用Euler公式：第二十三页，共四十九页，编辑于2023年，星期五定义，G为幅度因子由上式可得

可见，G=G（θ），由于不同的θ值代表不同的分量，所以幅度因子对于不同的傅里叶分量有不同的值。根据VonNeumann法的定义，要使方程的解保持有界，对于所有的θ值都应该有︳G︳≦1

相当于第二十四页，共四十九页，编辑于2023年，星期五

因为定义域在（-1,1）要使定义域在上式成立，只有

d≦1/2即为保持差分方程计算稳定的条件。因为，较稳定条件也可写成

不存在δ﹤0的情况。d≦1/2，-1≦1-4d第二十五页，共四十九页，编辑于2023年，星期五用VonNeumann法对对流扩散方程的（FTCS）格式进行稳定分析对流扩散方程：用FTCS离散第二十六页，共四十九页，编辑于2023年，星期五用VonNeumann进行稳定分析消去得第二十七页，共四十九页，编辑于2023年，星期五幅度因子实部虚部可见G是一个复变量。当c→0时，即对流扩散方程化成纯扩散方程，这个幅度因子即转化成实幅度因子，与前述讨论的结果一致。稳定条件可以从模︳G︳来讨论

第二十八页，共四十九页，编辑于2023年，星期五是的函数，当（-1,1）内变化时，的变化状态，进一步判别差分的稳定性和稳定条件为了研究G（θ）取得权值的条件，将对求二阶导数，得由数学分析可知，当时，取得极小值。第二十九页，共四十九页，编辑于2023年，星期五在边界上如果在边界上满足稳定条件，则整个计算过程都是稳定的。dc→第三十页，共四十九页，编辑于2023年，星期五第三十一页，共四十九页，编辑于2023年，星期五差分方程稳定的必要充分条件：这种方法对一维、二维的问题都适用。是一种普遍方法。第三十二页，共四十九页，编辑于2023年，星期五5守恒性物理概念：如果对一个差分方程在定义域的任一有限空间内作求和运算（相当于连续问题中对微分方程积分），所得的表达式满足该区域上物理量守恒的关系时，称差分格式具有守恒性。第三十三页，共四十九页，编辑于2023年，星期五

在Δt时间间隔内流入与流出某一区域中的通量之差等于该时间间隔中该区域内的Φ增量。第三十四页，共四十九页，编辑于2023年，星期五如何控制守恒性从控制容积积分法建立差分方程的角度看，满足下列两个条件时差分格式具有守恒性：

1）控制方程守恒型的。

2）在同一界面上各物理量（Φ及有关物性）及Φ的一阶导数是连续的。所谓连续，这里指的是从界面两侧的两个控制容积来写出的该界面的值是相等的。第三十五页，共四十九页，编辑于2023年，星期五6迁移性物理概念：如果对流项的某种差分格式仅能使扰动沿着流动方向传递，则这种格式称为具有迁移性。第三十六页，共四十九页，编辑于2023年，星期五7误差7.1误差的类型：1）截断误差是指以差分代替导数时略去的部分，与函数性质及变量有关。严格来说，是差分方程的误差而不是解的误差。2）舍入误差是指在差分方程的求解的过程中，特别是迭代求解时，由于计算机计算长度的限制而引起的。舍入误差的计算较困难，必须用概率分析的方法来定量估计。第三十七页，共四十九页，编辑于2023年，星期五7.2误差的物理概念从物理概念上理解，误差可分成阻尼误差、相位误差、混淆误差、传输误差四种。1）阻尼误差是由于离散引入的隐含的人工粘性所造成的误差。2）相位误差是由于离散的作用引起的。在有限差分计算中，不同的傅里叶分量有不同的对流流速，波长越大的分量对流流速越接近于正确的速度u，而波长较短的分量则以低于或高于u的相速传送这样各个傅里叶分量间的相位关系产生畸变，造成解的相位误差。第三十八页，共四十九页，编辑于2023年，星期五3）混淆误差是由于傅里叶分量的相互作用，能量逆转重新补充到长波分量中，从而使长波分量产生畸变甚至引起计算的不稳定。4）传输误差是由于不正确的对出流边界条件的额外要求引出的另一种类型的误差。第三十九页，共四十九页，编辑于2023年，星期五8举例将对流方程用（FTCS）显示差分离散，并分析其性能。用（FTCS）离散用VonNeumann法分析稳定性第四十页，共四十九页，编辑于202