偏微分数值解

1、偏微分方程数值解的研究摘要：本文从偏微分方程问题的起源着手，介绍从18世纪起，达朗贝尔首先研究弦振动方 程。然后介绍边界积分方程组的两类常见的数值解方法，投影法和机械求积法。最后介绍两 种非线性偏微分方程的展开解法，Exp-function方法和(G/G)展开法。关键词：投影法：机械求积法；Exp-function方法；(G /G)展开法；0偏微分方程问题的起源微积分创立不久之后，从18世纪初开始，达朗贝尔(Jean Le Rond dAlembert, 1717 一 1783)首先研究了弦振动方程。傅立叶(J. Fourier, 1768 一 1830)第一个 发现并求解了热传导方程，其平衡

2、态则用拉普拉斯方程刻画。这几类方程是最经 典的偏微分方程，由它们而产生大量的求解方法以及通解，但是，不存在求解偏 微分方程的一般方法。虽然19世纪初期柯西(A.L.Cauehy，1789 一 1851)、拉格 朗日(J. L. Lagrange，1736 一 1813)等人解决了一阶偏微分方程的求解问题，但 这种化为一阶常微分方程组求解的基本方法，对于二阶偏微分方程并不适用。从数学研究来讲，数学家们一直相信通解的存在，并且以此为出发点，先求 通解，最后确定常数和函数。1820年柯西首先证明了一阶常微分方程初值问题 解的存在性和唯一性。而且很快从实数推广到复数范围了，对于偏微分方程他只 需要证明

3、一阶偏微分方程组在复数范围内的解的存在性。正是在这里，他创造性 地发明了优函数方法：首先应用幕级数展开式给出形式解，然后与某个己知收敛 的幕级数相比较来证明形式解的收敛性，从而证明了解析解的存在性。30多年之后，魏尔斯托拉斯(K. WeierstraSS，1815 一 1897)的学生俄国数学 家科瓦列夫斯卡娅独立地证明了偏微分方程组柯西问题解析解的存在唯一性，即 柯西一科瓦列夫斯卡娅定理。柯西一科瓦列夫斯卡娅定理是偏微分方程理论中第 一个普遍的存在定理。这样普遍适用的、不依赖于方程类型的性质的一直引人关 注，是偏微分方程理论的一个主要发展方向；但另一方面，18世纪以来对几类 特殊方程的研究结

4、果表明，把微分方程分成不同的类型：双曲型、抛物型和椭圆 型方程。这些是二阶线性偏微分方程的三种基本类型，数学物理中常见的波动方 程、热传导方程及调和方程分别为它们的代表，其定解问题的提法、解的性质以 及求解方法，多数可以推广到创门所代表的三种类型上去。有些方程在区域的一部分可能是双曲型的，而在另外某一部分可能是椭圆型 的，在分界部分或者退化为抛物型的或者是不确定的，这样的方程称为混合型。 1923年F.G特里科米(F. Tricomi)首先研究了这样的方程，并得到了深刻的结果。关于各种各样的定解问题的提法，既有其物理上的依据，在数学上又能按照 一定的准则被证明是合理的，这些准则中最常用的是解的

5、存在性、唯一性以及解 关于定解数据的连续依赖性。如果一个定解问题的解是存在、唯一且连续地依赖 于定解条件，则这个定解问题就称为适定的，此概念是阿达玛(Jacuqes Hdamaard，1565 一 1963)提出来的。双曲型方程柯西问题的现代理论,是由阿达玛对二阶双曲型方程柯西问题的 先驱工作开始的。他通过构造在特征劈锥面上具有奇性的解-基本解来求解柯 西问题，并采用发散积分的有限部分的方法来克服所遇到的奇性困难。他的工作 经过里斯及索伯烈夫等人的发展，对广义函数论的建立是一个重要的推动，而阿 达玛的方法在广义函数论的框架中也得到了更清晰和完善的表达。黎曼(GF.B.Rimenan，1826

6、一 1866)依据所谓的狄里克雷原理断言狄里克雷 问题有解，魏尔斯托拉斯指出了黎曼这个论断的逻辑缺陷，并举出了反例。多年 之后，希尔伯特(D.Hilbert，1862 一 1943)给出了狄里克雷原理的完整无缺的证明。 (狄里克雷问题解的存在性证明，有庞加莱一佩隆扫除法，施瓦兹交替法，差分 法等等。)这方面的研究极大地推动了泛函分析的发展，也使得变分法成为研究 偏微分方程的强有力的工具。拉普拉斯方程的狄里克雷问题解的存在性，波动方程的特征初值问题解的存 在性等具体方程的一般理论问题，在19世纪也同样受到关注并分别推广到更为 一般的一类方程。20世纪30年代彼得罗夫斯基(1.G.Petrovky

7、，1901 一 1973)的工作成为偏微 分方程理论发展的又一个里程碑，他对古典分型理论进行了推广，并同时研究了 每一类方程的特点。1边界积分方程数值解在积分方程的解析解很难求或无法求出时，借助数值方法求解它的近似解就 显得非常重要。从已有的成果来看，方法的应用先于理论分析，也就是说，不少 方法已被提出或应用，但其数学基础的研究则不甚深入。随着应用的愈高要求， 近年人们已把着眼点转向方法的数学原理研究上。我国对奇异积分方程的数值 解法的研究起步较晚，但自82年全国积分方程会议强调开展这方面的研究后， 很快引起了一些学者的兴趣和关注，特别是数值方法的理论分析方面，近年来我 国学者做了不少工作。这

8、里仅考虑方程的常见数值解法。1.1投影法常见的投影法通常有配置法，Galerkin有限元法和最小二乘法。设中.(j = 1,2,.,n)是化的基函数，近似解v ex可表示为这组基函数的线性组合 hh hv =n c 中(x)( 1)j=1投影法都要归结为求解以c, = (j = 1,2,., n)为未知量的线性方程组。配置法：设e D是一组适当选取的配置点(数目为n个)，令Av - f在配置点上恒等于零就得到n个线性方程工Ap (t )c = f (t),i = 1,2,.,n(2)j i j ij=1一般说来，系数矩阵A%(t)是不对称的。配置点的选取要取决于积分算子的性 质和边界的几何状态

9、，它对于离散方程(2)的适定性有决定性的影响。Galerkin方法：在这种方法中，基函数系具有有限支集，系数七是通过下 面的Galerkin有限方程得到的工(Ap ,甲)c = (f，甲),i = 1,2,.,n(3)j i jij=1如果A是自伴的，则上面方程组的系数矩阵是对称的。最小二乘法：这种方法要求| Av - f|取得最小值。如果方程的解存在唯12(. D)一，则最小二乘法将导致一个对称，正定的线性方程组工(Ap.俱 M = (f, Api), i = 1,2,., nj=1以上方法都可以看作Galerkin-Petrov方法的特例。它们采用相同的试探函数，但 选择不同的检验函数。所

10、得到的线性代数方程组的系数矩阵都是满阵。1.2机械求积法选取一个适当的求积公式来近似积分算子Av有A (s) = hLk(k,t )v(t ),t = jh,h = (b一a)/n, j = 1,2,n(4)nj j jj=1然后赋值就可以求出离散矩阵的每一个元素，离散方程为hk(s ,t )v(t ) = f (s ),s = ih,i = 1,2,n(5)i j j i ij=1机械求积法对连续核可以借助Anselone聚紧收敛理论得到它的收敛性证明，但 是对边界积分方程由于核的奇性，则必须另辟蹊径。2两种非线性偏微分方程的展开解法2.1 Exp-function 方法描写浅水波的KdV方

11、程，它的一般形式如下：u + uux +=0KdV方程应用非常广泛，它也是许多领域的孤波现象的模型：如等离子体 声波，弹性杆中纵向色散波，低温下非线性晶格的声子波包的热激发等。我们考虑如下的KdV方程：u + 6uu +uu = 0利用变换u = u（&）,& = kx + wt代入上述方程，则方程转化为常微分方程：wu + 6kuu + k 3u= 0。现假设方程的解具有如下形式:d a exp（诣）f pb exp（用）我们平衡最高阶线性项u”以及最高阶非线性项uu:c exp(d + 7q)& + + c exp(-c-7p)&=q=一 七凉用)8c exp(2d + q)& + c e

12、xp(- 2c- p)& uu = 3寸4乙q b.exp( j& )3于是，我们有：d + 7q = 2d + 6q以及-c-7p = - 2c-6p从而我们得到c,d,p,q的关系为：p = c,q = d。取 p = c = 1, q = d = 1则上式为 u = aiexp(&) + aiexp(&) exp(&) + b + b exp(-&) TOC o 1-5 h z 0-1带回原方程中我们得到exp(-2&) + c exp(-3&) = 03c exp(3&) + c exp(2&) + c exp(&) + c + c exp(-&) + c3210-1-2A = (ex

13、p(&) + b + b exp(-& )4其中0-1c = -wa + 6ka2b + wab k3a b k3a 6ka a301 01 01 001 0c =-6ka2-2wa -8k3a + 2wab - 12kaa+ 8k3ab + 2wab2 -2wa b -4k3a b + 6ka2b2 + 12ka2b1 -11 00 00 01 01 -1c = -18k3a bb + 6wa b b + 18ka2b b + wa b3 + wa b2 6ka2b + k3a b311 0 -11 0 -11 0 -11 00 00 01 0一 k3a b2 + 6ka ab2 18ka

14、a 一 wa b 一 23k3a b 一 5k3a b 一 5wa b +12 0 00 1 00-10-10-1-10-10ka a b 一 12ka a b1 0 -11 -1 0c = -4wa b2 +12ka2b2 + 32k3a b2 - 4k3a b2 -12ka2 - 4wa b-24ka a b + 24ka abb + 4wa b 2 + 4wa b?b + 4k 3a b2b0-1010-101 -11 0 -11 0 -1c = 6ka2b b + 18kaa b2 12ka a b + 5wab b + 12ka ab b wa b3 k3a b3 + wa b b2

15、 6ka a b2-10 011 0101 11 0111101 01 001 00 1 0 TOC o 1-5 h z + wa b2 + 5k3abb2 18ka2b 23k3a b 6wa b b + 18k3a b b kaab 12kaa b 0 110 110011011 011 01110c = 2wa b2b + 4k3a b2b 4k3a b2 b + 2wa b3 8k3a b2 6ka2b221011010 101111016ka2b2 2wa b2 12ka2b + 8k3a b3 + 2wa b2b + 12ka ab1 011111 10 1 0111c = 6ka ab + k 3a b3 k 3a b b 2 + wa b3 wa bb2 6ka 2 b b310 10 11 010 11 0110 1上述式子中假若系数为0,得到如下的代数方程：C3= 0，C2= 0，0 0时:C 2 11 ; 2a12au &) = a c1sinhM 立肉+c2coshn 甲肉 c cosh(1肉+ c sinh1var肉12 c