主題
Search

亥姆霍茲微分方程——圓柱座標系


圓柱座標系中,尺度因子h_r=1h_theta=rh_z=1,因此拉普拉斯算符由下式給出

 del ^2F=1/rpartial/(partialr)(r(partialF)/(partialr))+1/(r^2)(partial^2F)/(partialtheta^2)+(partial^2F)/(partialz^2).
(1)

嘗試在亥姆霍茲微分方程中進行分離變數法

 del ^2F+k^2F=0
(2)

透過寫作

 F(r,theta,z)=R(r)Theta(theta)Z(z),
(3)

然後結合 (1) 和 (2) 得到

 (d^2R)/(dr^2)ThetaZ+1/r(dR)/(dr)ThetaZ+1/(r^2)(d^2Theta)/(dtheta^2)RZ+(d^2Z)/(dz^2)RTheta+k^2RThetaZ=0.
(4)

現在乘以 r^2/(RThetaZ)

 ((r^2)/R(d^2R)/(dr^2)+r/R(dR)/(dr))+1/Theta(d^2Theta)/(dtheta^2)+(r^2)/Z(d^2Z)/(dz^2)+k^2r^2=0,
(5)

因此方程已被分離。由於根據圓柱座標系的定義,解必須在 theta 中是週期性的,因此 (5) 的第二部分的解必須具有分離常數

 1/Theta(d^2Theta)/(dtheta^2)=-m^2,
(6)

其解為

 Theta(theta)=C_mcos(mtheta)+D_msin(mtheta).
(7)

將 (7) 代回 (5) 得到

 (r^2)/R(d^2R)/(dr^2)+r/R(dR)/(dr)-m^2+(r^2)/Z(d^2Z)/(dz^2)+k^2r^2=0,
(8)

並將兩邊除以 r^2 得到

 1/R(d^2R)/(dr^2)+1/(rR)(dR)/(dr)-(m^2)/(r^2)+1/Z(d^2Z)/(dz^2)+k^2=0.
(9)

(9) 的第二部分的解在 +/-infty 處不能是正弦的,對於物理解而言,因此該微分方程具有分離常數

 1/Z(d^2Z)/(dz^2)=n^2,
(10)

解為

 Z(z)=E_ne^(-nz)+F_ne^(nz).
(11)

將 (11) 代回 (9) 並乘以 R 得到

 (d^2R)/(dr^2)+1/r(dR)/(dr)+(n^2+k^2-(m^2)/(r^2))R=0
(12)

但這只是貝塞爾微分方程的一種修改形式,其解為

 R(r)=A_(mn)J_m(rsqrt(n^2+k^2))+B_(mn)Y_m(rsqrt(n^2+k^2)),
(13)

其中 J_n(x)Y_n(x) 分別是第一類第二類貝塞爾函式。因此,通解為

 F(r,theta,z)=sum_(m=0)^inftysum_(n=0)^infty[A_(mn)J_m(rsqrt(k^2+n^2))+B_(mn)Y_m(rsqrt(k^2+n^2))] 
 ×[C_mcos(mtheta)+D_msin(mtheta)](E_ne^(-nz)+F_ne^(nz)).
(14)

在 Morse 和 Feshbach (1953) 的符號中,分離函式為 f_1(r)=rf_2(theta)=1f_3(z)=1,因此斯塔克爾行列式為 1。

亥姆霍茲微分方程在更一般的 k^2 形式下也是可分離的

 k^2(r,theta,z)=f(r)+(g(theta))/(r^2)+h(z)+k^('2).
(15)

參見

圓柱座標系, 亥姆霍茲微分方程, 亥姆霍茲微分方程——極座標系

使用 探索

參考文獻

Moon, P. 和 Spencer, D. E. Field Theory Handbook, Including Coordinate Systems, Differential Equations, and Their Solutions, 2nd ed. New York: Springer-Verlag, pp. 15-17, 1988.Morse, P. M. 和 Feshbach, H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, pp. 514 和 656-657, 1953.

請引用為

Weisstein, Eric W. "亥姆霍茲微分方程——圓柱座標系。" 來自 —— 資源。 https://mathworld.tw/HelmholtzDifferentialEquationCircularCylindricalCoordinates.html

主題分類