偏微分方程的数值解方法及源程序.pdf

《偏微分方程的数值解方法及源程序.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第二十章偏微分方程的数值解自然科学与工程技术中种种运动发展过程与平衡现象各自遵守一定的规律。这些规律的定量表述一般地呈现为关于含有未知函数及其导数的方程。我们将只含有未知多元函数及其偏导数的方程，称之为偏微分方程。方程中出现的未知函数偏导数的最高阶数称为偏微分方程的阶。如果方程中对于未知函数和它的所有偏导数都是线性的，这样的方程称为线性偏微分方程，否则称它为非线性偏微分方程。初始条件和边界条件称为定解条件，未附加定解条件的偏微分方程称为泛定方程。对于一个具体的问题，定解条件与泛定方程总是同时提出。定解条件与泛定方

2、程作为一个整体，称为定解问题。§1偏微分方程的定解问题各种物理性质的定常（即不随时间变化）过程，都可用椭圆型方程来描述。其最典型、最简单的形式是泊松(Poisson)方程22∂u∂uΔu=+=f(x,y)（1）22∂x∂y特别地，当f(x,y)≡0时，即为拉普拉斯(Laplace)方程，又称为调和方程22∂u∂uΔu=+=0（2）22∂x∂y带有稳定热源或内部无热源的稳定温度场的温度分布，不可压缩流体的稳定无旋流动及静电场的电势等均满足这类方程。Poisson方程的第一边值问题为22⎧∂u∂u⎪+=f(x,y)(

3、x,y)∈Ω22⎨∂x∂y（3）⎪u(x,y)

4、=ϕ(x,y)Γ=∂Ω⎩(x,y)∈Γ其中Ω为以Γ为边界的有界区域，Γ为分段光滑曲线，ΩUΓ称为定解区域，f(x,y),ϕ(x,y)分别为Ω,Γ上的已知连续函数。第二类和第三类边界条件可统一表示成⎛∂u⎞⎜+αu⎟(x,y)∈Γ=ϕ(x,y)（4）⎝∂n⎠其中n为边界Γ的外法线方向。当α=0时为第二类边界条件，α≠0时为第三类边界条件。在研究热传导过程，气体扩散现象及电磁场的传播等随时间变化的非定常物理问题时，常常会遇到抛物型方程。其最简单的形式为一维热传导方程2∂

5、u∂u−a=0(a>0)（5）2∂t∂x方程（5）可以有两种不同类型的定解问题：初值问题（也称为Cauchy问题）-240-2⎧∂u∂u⎪−a=0t>0,−∞

6、件u(0,t)=g(t),u(l,t)=g(t)称为第一类边界条件。第二类和12第三类边界条件为⎡∂u⎤−λ(t)u=g(t),0≤t≤T⎢∂x1⎥1⎣⎦x=0（8）⎡∂u⎤+λ(t)u=g(t),0≤t≤T⎢∂x2⎥2⎣⎦x=l其中λ(t)≥0,λ(t)≥0。当λ(t)=λ(t)≡0时，为第二类边界条件，否则称为第三1212类边界条件。双曲型方程的最简单形式为一阶双曲型方程∂u∂u+a=0（9）∂t∂x物理中常见的一维振动与波动问题可用二阶波动方程22∂u2∂u=a（10）22∂t∂x描述，它是双曲型方程的典型

7、形式。方程（10）的初值问题为⎧∂2u∂2u2⎪=at>0,−∞0,0

8、定解问题是适定的。可以证明，上面所举各种定解问题都是适定的。§2偏微分方程的差分解法差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则

9、差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：（i）选取网格；（ii）对微分方程及定解条件选择差分近似，列出差分格式；（iii）求解差分格式；（iv）讨论差分格式解对于微分方程解的收敛性及误差估计。下面我们只对偏微分方程的差分解法作一简要的介绍。2.1椭圆型方程第一边值问题的差分解法以Poisson方程（1）为基本模