5章_圆柱坐标系下的分离变量法(07.12.28).docVIP

5.圆柱坐标系下的分离变量法 5. 1极坐标系下的拉普拉斯方程 考虑半径为a的一个薄圆盘，已知圆盘内部无热源，边界温度给定，且温度 分布u(x,t)随时间演化已趋于稳定，试求此时的温度分布 u(x, y)。上述定解问题 可表述为 、2u (x, y) = 0 (x, y) D ( u(x, y) x2*2^=g(x, y) (x, y^： ( 其中，v 表示圆盘的边界，即x2・y2二a2，D表示a 围成的内域。 对于二维平面场问题，即物理量的空间分布与z无关，当物体边界为矩形时, 采用直角坐标系比较方便。因为边界方程可方便地用直角坐标表示出来，如 x=0， x=a， y=0， y=b 但当物体边界为圆形时采用极坐标系可大为简化边界方程, 来方便。而在极坐标系下，拉普拉斯方程表示为 从而给问题的求解带 —(r ) .x x 2 +丄旦 汶 r2 2 从而( 2 1 ：： z ：：U、 1 c (r ) 2 t - 0 r :x :x r ： u( r, ) +g() 注意到定解问题(，通常称之为狄里克莱 (Dirichlet )问题，也称第〔边值问题。 若(，则称相应的定解问题为牛曼(Newman)问题, 也称第U边值条件。若(, 则称相应的定解问题为罗宾(Robin)问题。此外，本题研究内域 D中的温度，通 常称为内问题。实际应用中，可能遇到求圆形孔洞外围的温度场或电势场分布问 题，通常称为外问题。 现在回到求解形如(，设 u(r, ) = R(r)：：() 代入( 1 d dR 1 ^d2：：」门 (r ) 2 R 2 =0 r dr dr r :: 两边同除以R； r (u二RG为非零解)得 rd dR 1 d2「 (r ) 2 R dr dr 由于等式左边是关于r的函数，右边是关于「的函数，从而只能有 左边=右边=常数 设这个常数为-’，则得到两个常微分方程 :(p ( H 0 ( 和 r(rR (r)) - ■ R(r) = 0 ( 或者 r2R (r) rR (r)- R(r) = 0 ( 如同在直角坐标系下求解偏微分方程定解问题一样， 我们将首先构造与定解问题 相应的特征值问题，通过求解特征值问题得到平方可积函数空间 L2中的一组完 备正交函数系，再将解按完备正交函数系展开，最终得到级数形式的解表达式。 为此，首先考虑方程(，定解条件需要 2个，但(D r = 0，补充定解条件如下 u(0, M := ( 从而可分离出 R(0) ( 但从r =a处的边界条件中无法分离出关于 R(r)的边界条件。从而无法由 方程(，其定解问题也需要 2个，但(叮C )，但深入考虑本问题的特点后，应 该有u (r, ) = u (r J： 2二) u (r ,)= u (r, 2二) 因为(r,)和(r , 2二)表示同一点。形如(r「，其定解条件应该有4个。除r=a 处的边界条件外，还应该有3个定解条件。有界性条件(，他们或者由问题的物 理性质决定，或者由区域D的几何性质决定。像这样由问题的物理性质决定，或 者由区域D的几何性质决定，而无需在定解问题中明确提出的边界条件， 通常称 为自然边界条件。自然边界条件是隐含在定解问题本身之中的边界条件。 由周期性条件( ：4 ) 一 ：( 2二) ( ：()=：(2二) ( 它们与方程( G(「)= C Dsin、.. ■「 ( 注意到(:•: ( ) 2二，故 、、■ - 0,二 1, -2, 即特征值 ，=n2 (n =0,1,2,…) ( 相应的非平凡解为 ( ) = Cn cosn? Dn sinn ( 由于cosn®和sinn®是线性无关的，所以特征值是 2重简并的(除九=0外)，即 每一个特征值 n(n_1)对应有2个线性无关的特征函数 cosn「(n =1,2, 3…)和 sinn「(n =1,2,3,)。所有正交函数组成函数空间L2(0,2 )完备正交函数系。 )n ( )■ ■ 1,cosn ,sin n ； (n = 1,2,3,) 现在将’n二n2代入(，求解欧拉(Euler)方程 r2R (r) rR (r)—n2R(r)= 0 (0 ： r ：： a) ( (1)当 n =0时，R (r) _ 1 R (r) ~ ~r In R(r) = _ln r c R (r) = Dr 1 R(r) = A B0 In r (2)当n=0时， 令r二d贝U dt 1 R (r )= R (t) 丁 =R (t)- dr r d d 1 1 R (r )= — (R (r)) - (R (t)-)=飞[R (t)-R (t)] dr dr r r 原欧拉方程( R (t)- n2R(t) = 0 从而 R(t) = Anent Bne Bj』 综合式（  Dn = 0