莱恩-埃姆登方程

✍ dations ◷ 2025-04-03 12:49:01 #天体物理学,恒星天文学,常微分方程

莱恩-埃姆登方程(Lane–Emden equation)是天文物理中一个表现自引力势能,球对称多方流体的无量纲泊松方程。此方程名字由来于强纳生·荷马·莱恩与罗伯特·埃姆登。此方程的解表示了恒星在半径 r {\displaystyle r} 代表核心的压力与密度。 n {\displaystyle n} 是多方指数;多方指数与代表气体压力及密度的多方方程有关系。

P {\displaystyle P} 是代表压力, ρ {\displaystyle \rho } 则是密度,而 K {\displaystyle K} 则是比例常数。标准的边界条件则是 θ ( 0 ) = 1 {\displaystyle \theta (0)=1} θ ( 0 ) = 0 {\displaystyle \theta '(0)=0} 。因此该方程的解是描述恒星压力和密度与半径的关系,并且给定的多方指数 n {\displaystyle n} 也是多方球的多方指数 n {\displaystyle n} 。流体静力平衡与势能、密度、压力梯度有关;泊松方程与势能、密度有关。

在物理学上,流体静力平衡与势能梯度、密度和压力梯度相关,而泊松方程则可以是势能和密度的关系式。因此如果有一个方程可以进一步指出压力和密度如何互相反映,就可以得到一个解。以上多方气体的特定选项在数学上陈述了这个问题,尤其是该陈述特别简洁并推导出了莱恩-埃姆登方程。这个方程对于恒星等自引力势能气体球是相当有用的近似,但它的假设通常是受到限制。

考虑到自引力势能、流体静力平衡下的球对称流体、质量守恒这些状况,就可使用以下连续性方程:

这里 ρ {\displaystyle \rho } r {\displaystyle r} 的函数。流体静力平衡的公式成为:

m {\displaystyle m} 也是 r {\displaystyle r} 的公式。再一次求导数可得:

这里已经使用一个连续性方程取代质量梯度。再将方程两侧乘上 r 2 {\displaystyle r^{2}} ,并将带有 P {\displaystyle P} 的导数的项置于左侧,方程成为:

方程两侧除以 r 2 {\displaystyle r^{2}} ,在某些意义上这是一维形式所需的方程。此外,如果我们以多变方程 P = K ρ c 1 + 1 n θ n + 1 {\displaystyle P=K\rho _{c}^{1+{\frac {1}{n}}}\theta ^{n+1}} ρ = ρ c θ n {\displaystyle \rho =\rho _{c}\theta ^{n}} 代入,可得到:

将常数聚集并以 r = α ξ {\displaystyle r=\alpha \xi } 取代:

最后得到莱恩-埃姆登方程:

同样地,也可以使用泊松方程进行推导:

我们可以透过以下数学公式以流体静力平衡取代势能梯度:

最后也可以得到莱恩-埃姆登方程。

n {\displaystyle n} 只在3个值时有解析解

如果 n = 0 {\displaystyle n=0} ,方程成为:

重新整理并进行一次积分后的公式成为:

公式两侧都除以 ξ 2 {\displaystyle \xi ^{2}} ,并且再积分一次后得到:

边界条件 θ ( 0 ) = 1 {\displaystyle \theta (0)=1} θ ( 0 ) = 0 {\displaystyle \theta '(0)=0} 暗示积分常数是 C 0 = 1 {\displaystyle C_{0}=1} C 1 = 0 {\displaystyle C_{1}=0}

n = 1 {\displaystyle n=1} ,方程可展开如下:

两端都乘以 ξ 2 {\displaystyle \xi ^{2}} 可得到 k = 1 {\displaystyle k=1} n = 0 {\displaystyle n=0} 的球贝索函数。套用了边界条件以后的解将是:

在经过一连串取代的步骤后,方程可以有进一步的解:

n = 5 {\displaystyle n=5} ,方程的解将是循着径向的无限大值。

一般情形下莱恩-埃姆登方程的解必须以数值积分方式求得。许多数值积分的标准解法要求该问题必须以一阶常微分方程表示,例如:

在这里 ϕ ( ξ ) {\displaystyle \phi (\xi )} 被视为无量纲质量,而质量可使用 m ( r ) = 4 π α 3 ρ c ϕ ( ξ ) {\displaystyle m(r)=4\pi \alpha ^{3}\rho _{c}\phi (\xi )} 表示。相关的边界条件是 ϕ ( 0 ) = 0 {\displaystyle \phi (0)=0} θ ( 0 ) = 1 {\displaystyle \theta (0)=1} 。第一个方程表现了流体静力平衡,而第二个方程则表示质量守恒。

已知如果 θ ( ξ ) {\displaystyle \theta (\xi )} 是莱恩-埃姆登方程的解,那么完整的解方程将是 C 2 / n + 1 θ ( C ξ ) {\displaystyle C^{2/n+1}\theta (C\xi )} 。和这方式相关的解则称为“同调”,而转换的过程是同调性的。如果我们选择不变的变量达到同调性,就可以将莱恩-埃姆登方程降一阶计算。

而这类可选择的变量有多个,一个适当的选择是:

我们可以将相对于 ξ {\displaystyle \xi } 的变量的对数微分,得到:

最后,我们将以上两个方程相除以消去应变量 ξ {\displaystyle \xi } ,留下:

以上即为单一一阶方程。

同调性不变的方程可被视为自主对方程:

这些方程的解的形式可透过以下线性稳定性分析来决定。方程的临界点(当 d V / d log ξ = d U / d log ξ = 0 {\displaystyle dV/d\log \xi =dU/d\log \xi =0} )和雅可比矩阵的特征值、特征矢量如下表所示:

相关

  • 法国电视法国首次播出电视节目是在1931年,这使得法国成为世界上最早播出电视的国家之一。2005年开始,法国开始播出数字电视讯号的电视节目。和其他欧洲国家一样,法国采用DVB-T作为数字
  • 鸽形目small/small鸽形目(学名:Columbiformes)在动物分类学上是鸟纲中的一个目,现在仅包含一科——鸠鸽科。今天分布广阔、种类众多的鸽子与已经灭绝的愚鸠(渡渡鸟)都属于鸠鸽科。沙鸡科(Pteroclidae
  • 圣米格尔-德图库曼坐标:26°49′50″S 65°12′13″W / 26.83056°S 65.20361°W / -26.83056; -65.20361圣米格尔-德-图库曼(西班牙语:San Miguel de Tucumán)是阿根廷城市,为该国最小省份——图
  • 美国联邦储备委员会联邦储备委员会(英语:Federal Reserve Board of Governors,简称联储会)是美国联邦储备系统的主要管理机关。该会负责监管132个联邦储备银行,帮助推行货币政策。联储会理事由总统
  • 选民选民,又称为选举人,是在一个选举中有权参与投票的民众。在普选中,选民就是当地城市的市民或本国公民。在一家上市公司,选民就是股东,包括散户。选民权利除了投票,还有弃权和投弃权
  • 隐藏式阴茎隐匿型阴茎一种病态阴茎。多见于肥胖体形的儿童,其特点为包茎,阴茎周围皮下脂肪很厚,阴茎皮肤发育不良,阴茎海绵体及尿道海绵体发育欠佳,阴茎体不能进入到阴茎皮肤及包皮腔内。因
  • 中华龙鸟属中华龙鸟属(属名:,意为“中国的蜥蜴翅膀”,早期也译为中国蜥翼龙)是目前所发现拥有化石化羽毛痕迹的恐龙中,年代最早而且最原始的,也是辽宁省热河群第一个发现的恐龙化石。中华龙鸟
  • 克罗地亚独立战争克罗地亚战争指的是1991年到1995年之间,克罗地亚从南斯拉夫联邦独立出来时,克罗地亚人和塞尔维亚人之间因民族对立而引发的战争。克罗地亚军的种族灭绝罪行也受到非议,但国际上
  • 玛蒂娜·罗斯布拉特玛蒂娜·罗斯布拉特出生于1954年,是美国的一位律师,作家,和企业家。她1981年毕业于加州大学洛杉矶分校(UCLA),取得了法律与MBA博士学位。毕业后在华盛顿律师事务所从事通讯卫星
  • 极紫外辐射极紫外辐射(英语:Extreme ultraviolet radiation)或高能紫外辐射,简称EUV、XUV,是波长在124nm到10nm之间的电磁辐射,对应光子能量为10eV到124eV。自然界中,日冕会产生EUV。人工EUV