§5.2 Poisson方程的边值问题 示为方程的边值问题均可表Piosson ?? ?íì =ú?ùê?é +?? ?-=D )2()( )1(,)( Mgunu MMhu s ba t 格林公式 入导出积分的工具它们的积分公式需先引 的积分公式,为导出类问题需要先导出它们 为用格林函数法求解这为混合问题 问题为为狄氏函数其中当 -- 11 == ,0, ;0;,0 ba ba Neumen 一、格林公式 1、为何引入格林公式 (1)积分公式:所谓积分公式即解的积分表达式 dxxfyxfdxdu ò=?= )()(:若 òò= =??+??= dxdyyxfyxu yxfdyyudxxudu ),(),( ),( 则 进一步 积分可使问题向前推若不能直接积分,分部 从微分号下解脱出来,挂上积分号将未知函数 xedxdy x cos=如: ò òò += == xdxexe xdexdxexy xx xx sincos coscos)( 则 ò-+= xdxexexe xxx cossincos [ ]xxx exexdxe sincos21cos +=\ò . , 解脱出来的工具 下从微分算符格林公式是将未知函数 方程中均含有我们要求解的三类数理 D D 2 2 2 2 2 2 ),,( zyvxvzyxv ??+??+??=D即 ?)(),( MvMhv 如何得对于 -=D (2) 导出格林公式 于是自然想到由 ò Dt tvdu 上具有连续一阶导数导数在 中具有连续的二阶在和设 - t t),,(),,( zyxvzyxu òò ò ?×?-?×?=D tt ttt vdudvuvdu )(:则 òò ?×?-×?= ts ts vdudvu r 2. 格林第一公式 ò ò ò ×?= ?×?+Ds t t s tt rdvu vduvdu )3(ò ??= s sdnvu )4(: ò òò ??=?×?+D t st stt dnuvvduudv同样 )5()(ò òò ??-??=D-D t st stt dnuvnvuudvvdu 3、格林第二公式 的值或就有可能求得 的边值知值联系起来这意味着若 的边的点与将意义 hvhu n v n uvu n u n vvuvuvu -=D-=D ? ? ? ? ? ? ? ?DD ,,, ,,,,,,)1(: 对称vu,)2( 的解 方程的边值问题就有可能导出求 第二公式使用和对 则的求法后面我们专门会讲易求若 因为为此我们引入点源函数 Poisson zyxu G MMG Greenz)y,G(c,),,( )(G ),( 0 的解。就有可能求得 则由上述格林公式,如已知 中已知一个、。若、含有两个未知函数 因为方程中利用上式显然不足以解 )( ,0, , ),()3( Mhu vv vuvu Mhu -=D =D -=D 即点源产生的场 使它满足中引入在 )6(,)( ),( 0 0 td t ?--=D MMMG MMG 202020 0 )()()( 4 1)4.3.5~3.5(, zzyyxxr rGM ---= =? pt 有则由 二、积分公式---格林函数法 1. 泊松方程的基本积分公式 )()()( ),()()(),( 0 00 MhMMMu MMGMuMuMMG -= D-D× d 的区域积分区间应挖掉 第二公式故为应用 为奇点以注意到 et Green MMMG ,),( 00 :)6()1( uG ×-× òò ò -- - --= ? ?- ? ? ee e tttt tt ttd t )7()()( ])(),([ 0 0 dMhGdMMu dnGMunuMMG即 : .)1().1()(, )( , 现来化简 的解即为必满足则此式 的积分表达的含有化简整理得到 若能由此式导出的和注意此式是由 Mu GMu (1)(6) 0)()( 00 =--\-? MMM dtttt ee 中在Q 0)( 0 =-ò - td ett dMMu于是 的距离到 由 右边 M 4 1)( : 0 0 Mr rGMMG - =?--=D pdQ )( 4 1),( : 0 0 es pe s e e 的距离为上任一点到此时 在 MM MMG Q = spe s e e s s du dnuMMG nò ò = ? ? 4 1 ),(: 0于是 s pe e es dunò= 2 4 0 0 0)( ? =×= e e Mun dsurMu r ò? = sp 200 4 1lim)(由 于是 )(4 1 4 1 0 0 2 Mudu drrudnGu ? == ? ?-= ? ? ò òò e s ss spe sps e ee 有 )7(将以上计算代入 ò ò -= --?? t s t s dMhMMG MudnGunuG )(),( )()(: 0 0于是有 ss t ss t dnGMudnuMMG dMhMMGMu òò ò ? ?- ? ?- = )(),( )(),()(: 0 00于是 积分式表示出 的值可由上中任一点在即 0)( MMhu t-=D 可不换的 但由于对换与而变将上式中随 的值即任一点为了表示得更一般些 GMMGMMG MMzyx M \= ),(),( ,,,, , 00 0 0 )9( )(),( )(),()( 0 00 0 0 00 基本积分公式--- ? ?- ? ?+ = òò ò n GMud n uMMG dMhMMGMu ss t s t 求导是对 取面元、体元分别表示对其中 00 000 .M, Mn dd st 2.泊松方程边值问题的积分公式 )()(1| MfMgu == bs )10()(| )()()2( )1( ?í ì = -=D? ty ü Mfu MhMu s 即 )(|)2( .0)1( Mgu =? = ?í ì sb a 即 中若 (2) (1) GMMG 的 次边界条件此时选择满足第一类齐 0|),( 0 =s )12()( )(),()()9( )11(0|),( ),(),( 0 0 000 0 00 ò ò ? ?- =? ?í ì = -=D s s s s t d dn GMf dMhMMGMu MMG MMMMG 则 即 就可完全确定 的解求出的一旦满足的解 它给出了狄氏问题称为狄式积分公式 (10),)MG(M,)11(, , 0 (10) ?í ì = -=D? ?í ì )(| )()( )2( )1( Mfu MhMu n s ?í ì = =D 0| )M(M,-G)MG(M, 0 0 s d nG 满足此时若选 为狄式格林函数的称满足 )MG(M, 0\ (11) )()(|)2(.0(2) MfMgun ==?= ?í ì sb 即中若 (1) (2) 。,较复杂,在此不讨论即方程和边界条件矛盾 而不会满足稳态方程间而变 随时内温度会不断升高示边界绝热,既然如此 却表点热源产生的影响,而可看成在 在物理上但此式不存在。因为 ),() ( 0M )M(M,-G 0 0 0 MMG n G d t d s -=D =?? =D ò ò + = s t sa t 000 00 )(),(1 ))h(MMG(M,u(M): dMgMMG d可得虽由 (9) ?? ?í ì =ú ? ù ê? é + ? ? -=D )14(0 )13(),(),( ),( 00 0 s ba d GnG MMMMG MMG 满足选 G)2( )1( 的称为定解问题 ?í ì :0)3( 根据前面经验均、中若 1 ?í ì ba(1) (2) )(),(1)( 0 MgMMGnGunGG sa s =??-?? )15()(),(1 )(),()( 000 00 sa t s t dMgMMG dMhMMGMu ò ò + = 得则将 uG ×-× (14)(2) :的解为于是由 ?í ì(1) (2)(9) 3. :积分公式的物理意义 ?í ì = ?-=D )(| )()(: Mfu MMhMu l s二维积分公式 ò ò ? ?- = l dlnGMf dMhMMGMu 0 0 0 000 )( )(),()( s s ?í ì = -=D 0| ),(),( 00 lG MMMMG d其中 ? í ì 和边界面上源产生的场的第二项 体内源产生的场的和第一项 : :(12) (15) 4. 附:证明 ),(),(0 )( 00 MMGMMGG n G Mhu = ?? ?íì =ú?ùê?é +?? -=D 中的 s ba 具有对称性的 为了更普遍我们不妨证 GMfu n u Mhuu ?? ?íì =+?? -=+D )()( )( sba l )()( 0 MhuMhuu -=D?-=+D =l lQ ?? ?íì =ú?ùê?é +?? -=+D )2(0),(),( )1(),(),(),( 11 111 MMGMMnG MMMMGMMG ba dl 设 ?? ?íì =ú?ùê?é +?? -=+D )4(0),(),( )3(),(),(),( 22 222 MMGMMGn MMMMGMMG ba dl 将若不同时为 0,0, 1abaQ [ ]),(),(),(),( ),(),(),(),( 1221 2112 MMGMMGMMGMMG MMGMMGMMGMMG -+ D-D l ),(),(),(),( 1221 MMMMGMMMMG dd -= 并用第二格林函数积分,中对两边在 ttt e d)( - ( ) ( ):,, 12 MMGMMG ×-×(1) (3) s s d MMGnMMG MMGnMMG ò ú ú ú ú ? ù ê ê ê ê ? é ? ?- ? ? ),(),( ),(),( 21 12 )5(),(),( 2112 MMGMMG -= 00 ),(),( ),(),( 21 12 =×+ ú ú ú ú ? ù ê ê ê ê ? é ? ?- ? ? ba s MMGnMMG MMGnMMG :),(),( 212 MMGMMG ×-×又将(2) (4) 0 ),(),( ),(),( 21 12 = ú ú ú ú ? ù ê ê ê ê ? é ? ?- ? ? s MMGnMMG MMGnMMG ),(),( 0 1221 MMGMMG =从而得 的右边等于同样可得 ),(),(0 ),(),( 122 1221 MMGnMMGn MMGMMG ? ?- ? ?×1 = 将若 得代入 b (5) (4)(2) (5) ?? ?í ì =ú ? ù ê? é + ? ? -=D-=+D 具有对称性的 或 满足总之 GGnG MMGMMGG 0 ),()(,( , 00 s ba ddl 三、小结 :)(1 的解为、 Mhu -=D 000 00 ),()(1 )(),()( sa t t dMMGMg dMhMMGMu ò ò + = 的解为其中 ),(),(),( 000 MMMMGMMG d-=D 0 0 0 00 )( )(),()()(| )( 2 s t s ts dnGMf dMhMMGMuMfu Mhu ? ?- =? ?í ì = -=D ò ò的解 、 的解 为其中 ?í ì = -=D 0| ),(),( ),( 00 0 s d G MMMMG MMG 的解为其中 ?? ?í ì =ú ? ù ê? é + ? ? -=D 0 ),(),( ),( 0 0 s ba d GnG MMMMG MMG :)( )( 3 的解为、 ?? ?í ì =ú ? ù ê? é + ? ? -=D Mgunu Mhu s ba 000 00 ),()(1 )(),()( sa t t dMMGMg dMhMMGMu ò ò + = 即可。各类边值问题相应的 只要求出泊松方程的值问题 的边由此可见欲求泊松方程 G ,