第十四章 格林函数--偏微分方程解的积分表示解偏微分方程主要有两种方法:A: 数理方法中的分离变量法:正交的多项式或无穷级数解,但需要齐次边界条件B: 理论物理中的Green函数方法:既是简单的有理形式解,又允许任意的边界条件!1,Green函数(GF)的意义:物理上:点源产生的场(函数)在时空中的分布特别是它在空间是源函数;在时空是传播函数See below)数学上: 具有点源的偏微分方程在齐次边界条件或者无界区域、初值条件下的解2,GF的分类:边界值GF:即源函数;初始值GF: 即传播函数3,Green函数的性质: 1)对称性:,它与定解问题相关,即与厄米性相关See 4 below) 2)时间传播函数没有对称性:.(因果律引起) 3)存在的必要条件:设方程,若λ是对应齐次方程的本征值,即 和附加齐次边界条件,则不存在这是因为既有点源:矛盾于又无流: 本征值问题存在,但是没有激发,物理上自相矛盾!平面波球面波和柱面波均是Laplace Equation的解,但不是Possion Equation的解球、柱面波分别来自于时(散射问题)渐近行为:4,Green函数的边值条件: 选取边值条件具有人为性,但要求简单并保证算子的厄米性。
1)齐次边值条件:矛盾于上述3,3),详见下面14.1.2末 2)有解性(解收敛):--基本解5,Green函数的用途:偏微分方程的积分解法:1)求, 2)利用迭加原理给出物理问题的积分形式解;Green函数的奇点与元激发的能量和寿命有关 6,Green函数的求法:1) 特殊方法:2)本征函数展开法:相应算子在同一边界下的本征函数作为基矢3)方程齐次化方法:将非齐次项变成边值条件和初值条件4)积分变换法:Laplace Transforms,Fourier Transforms.5) 形式解:算子运算14.1 Green函数与偏微分方程1,定义:Green函数(源函数,影响函数,传播函数,传播子) 数学上,含点源的偏微分方程在一定的边界条件或者初始条件下的解; 物理上,点源在一定物理条件下产生的场这种解(场)在时空中的分布与传播例1, 无界空间Possion Equation: 基本解---无界空间Green函数的叠加(see below for the solution).例2, 无界空间Helmholtz Equation: 例3, 无界空间波动方程:在含时Green函数中,为方便计,我们将它简记为2, Helmholtz Equation and Laplace Equation解的积分形式(在定解问题中求) 设定解问题: 假设已经求出,方法见14.3.作算符运算:得对上式左边利用Gauss公式其中可得如下的Green恒等式:See 童裕孙等,高等数学下册,2nd 版,PP154-156&PP161-162.故有形式解 (5)在已知的前提下,解(5)也不是用表示的最后形式。
这是因为在边界上还未知(最多知道他们在的线性组合—BCs). 幸好的边界条件还可以以多种形式提出;只要选定的边界条件为齐次,则或者其线性组合就可为已知的边界条件,从而最后确定 例如:1) 第一类边界条件: 2) 第二类边界条件: 3)第三类边界条件:得到 这些形式解的前提是要已知或者可求出实际问题中,可能就不存在例如在第二类边界条件中,的问题[Possion Equation, 点源存在,但边界“流”为零-- 物理上不通(既产生又绝缘--矛盾),数学上无解]再例如,构成本征值问题的,不存在相对应的方程(3)和(4). 这是因为即使在边界条件完全相同的情况下,(无源)而方程(3)的源为 在此情况下无解(或适当修改Green函数的意义),其实分离变量法已解决了问题,GF就多此一举了GF要解决更为复杂的物理问题!3. Green 函数的物理意义以为边界的区域,既无论原方程是否齐次(即内有、无源),又无论原边界条件是否齐次(即上有、无源),GF总是定义在内除一点(点源)以外方程的非齐次项处处为零(问题本身总要有非齐次项,如源于边界条件)但要化成齐次边界条件的定解问题的解因此GF是“点源影响函数” 或者 “作用的传播函数”。
对于所讨论的线性方程而言,一旦知道了相应问题的Green函数,只要再做两个积分,把原方程的非齐次项所反应的连续源分布对各点所产生的影响线性迭加起来,便给出原问题的解这是线性迭加原理的最成功应用齐次方程的本征值问题的本征值解可用于表示相应非齐次方程定解问题的GF(求法见最后一节)14.2 Green函数的性质1. Green函数由线性算子和边界条件和初始条件决定:,加上齐次边界条件和初始条件2. Green函数的叠加性1). 非齐次方程特解()+齐次方程通解(=0).例如:2). 3. Green函数的对称性:若算子是厄米的,则由产生的有 特别地,对于实变Green函数,1) Helmholtz方程: 作,则方程右端变为,而左端 利用Guass公式和Green恒等式,上式变为 [这是因为Eqs.(2)和(4)的行列式为零] 故: .2) Green Equations: 作 则方程右端变为,而左端 (因为的厄米性)故 推论,,如果为实变函数总之,物理意义:点源处产生的势传播到处:, 再对整体求和,加上边界条件产生的势,总和为处的势分布4. Green函数的奇异性及其随空间维数降低的减弱性。
1).(有)无界区域Green函数----基本解.a).无界区域: 设解为(详解见下),则基本解为有限形式,且中心对称(与有关)以及关于发散(奇异性与维度有关)b).有界区域:, 和 . 因为基本解有奇异性,所以直接影响分布点源通过边界,i.e.,对处的分布有直接影响2). Helmholtz Equation 的解在附近的奇异性(具体求解见下节):3D: 2D: ,对数发散; 1D: 连续,,导数发散由此可见:Green函数的奇异性随空间维数的降低而减弱14.3 Green函数的求法1. 方程齐次化法:将非齐次项变为边界条件1) Helmholtz Equation:a). 1D: b). 2D: 移动原点且改用极坐标系,其解为是对数分散,或确定 特别(Laplace Eq.): c). 3D: 设散射中心 ,球对称, 请:2)1D波动方程: 冲量定理法: 3)1D 输运方程: 对于第二类边界条件,将上面2)and 3)中的sin换为cos就行了2. 积分变换法:1)Fourier Transform.a). 3D Possion Equation: 作其反演: 再作又因为所以3D Possion Equation在TF以后变为进而积分 b). 3D Helmholtz Equation:如果,则this is the Coulomb screened potential.c). 3D波动方程: (源位于处). Taking FT for t, one has 于上述结果b比较得 故: 为超前/推迟Green函数。
2)Laplace Transform (3+1)D. a). 波动方程:作 LT for t 于是上述方程变为: 解此方程得: 设,则推迟势 令,则上式可表示为(推迟势)这是因为à 这里只有推迟势b). 热传导方程:作LT for t and FT for 再作 Inverse LT and Inverse FT: 3. 本征函数展开法: 2D Possion Equation: 选满足边界条件的正交函数系:正交归一性: 设带入方程得:其中 在此方程两边同乘以,再积分得到: 故 Green函数的解析性研究,即奇点与元激发的能量和寿命有关,见相关专业教材。