主題
Search

Lane-Emden 微分方程


LaneEmden

Lane-Emden 微分方程是一個二階常微分方程,產生於恆星內部結構的研究,也稱為多方氣體微分方程。其表示式為:

 1/(xi^2)d/(dxi)(xi^2(dtheta)/(dxi))+theta^n=0
(1)
 1/(xi^2)(2xi(dtheta)/(dxi)+xi^2(d^2theta)/(dxi^2))+theta^n=(d^2theta)/(dxi^2)+2/xi(dtheta)/(dxi)+theta^n=0
(2)

(Zwillinger 1997, pp. 124 and 126)。它具有以下邊界條件

theta(0)=1
(3)
[(dtheta)/(dxi)]_(xi=0)=0.
(4)

上方展示了theta(xi) 對於 n=0、1、2、3 和 4 的解。 n=0、1 和 5 的情況可以解析求解(Chandrasekhar 1967, p. 91);其他情況必須透過數值方法獲得解。

對於 n=0 (gamma=infty),Lane-Emden 微分方程為:

 1/(xi^2)d/(dxi)(xi^2(dtheta)/(dxi))+1=0
(5)

(Chandrasekhar 1967, pp. 91-92)。直接求解得到:

 d/(dxi)(xi^2(dtheta)/(dxi))=-xi^2
(6)
 intd(xi^2(dtheta)/(dxi))=-intxi^2dxi
(7)
 xi^2(dtheta)/(dxi)=c_1-1/3xi^3
(8)
 (dtheta)/(dxi)=(c_1-1/3xi^3)/(xi^2)
(9)
 theta(xi)=intdtheta=int(c_1-1/3xi^3)/(xi^2)dxi
(10)
 theta(xi)=theta_0-c_1xi^(-1)-1/6xi^2.
(11)

邊界條件 theta(0)=1 給出 theta_0=1c_1=0,因此:

 theta_1(xi)=1-1/6xi^2,
(12)

並且 theta_1(xi)拋物線

對於 n=1 (gamma=2),微分方程變為:

 1/(xi^2)d/(dxi)(xi^2(dtheta)/(dxi))+theta=0
(13)
 d/(dxi)(xi^2(dtheta)/(dxi))+thetaxi^2=0,
(14)

這是球貝塞爾微分方程

 d/(dr)(r^2(dR)/(dr))+[k^2r^2-n(n+1)]R=0
(15)

其中 k=1n=0,因此解為:

 theta(xi)=Aj_0(xi)+Bn_0(xi).
(16)

應用邊界條件 theta(0)=1 得到:

 theta_2(xi)=j_0(xi)=(sinxi)/xi,
(17)

其中 j_0(x)第一類球貝塞爾函式(Chandrasekhar 1967, p. 92)。

對於 n=5,進行 Emden 變換:

theta=Ax^omegaz
(18)
omega=2/(n-1),
(19)

這會將 Lane-Emden 方程簡化為:

 (d^2z)/(dt^2)+(2omega-1)(dz)/(dt)+omega(omega-1)z+A^(n-1)z^n=0
(20)

(Chandrasekhar 1967, p. 90)。經過進一步的推導(此處未重現),方程變為:

 (d^2z)/(dt^2)=1/4z(1-z^4)
(21)

然後,最終得到:

 theta_5(xi)=(1+1/3xi^2)^(-1/2).
(22)

在 中探索

參考文獻

Chandrasekhar, S. 恆星結構研究導論。 New York: Dover, pp. 84-182, 1967.Iyanaga, S. and Kawada, Y. (Eds.). 數學百科詞典。 Cambridge, MA: MIT Press, p. 908, 1980.Seshadri, R. and Na, T. Y. 工程邊值問題中的群不變性。 New York: Springer-Verlag, p. 193, 1985.Zwillinger, D. 微分方程手冊,第 3 版。 Boston, MA: Academic Press, pp. 124 and 126, 1997.

在 中被引用

Lane-Emden 微分方程

請引用本文為

Weisstein, Eric W. "Lane-Emden 微分方程。" 來自 --一個 資源。 https://mathworld.tw/Lane-EmdenDifferentialEquation.html

學科分類