主题
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)

在 Wolfram|Alpha 中探索

参考文献

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.

在 Wolfram|Alpha 中被引用

Lane-Emden 微分方程

请引用本文为

Weisstein, Eric W. "Lane-Emden 微分方程。" 来自 MathWorld--一个 Wolfram Web 资源。 https://mathworld.net.cn/Lane-EmdenDifferentialEquation.html

学科分类