レーン＝エムデン方程式

出典: フリー百科事典『ウィキペディア（Wikipedia）』
移動: 案内、 検索
白色矮星。レーン＝エムデン方程式を用いて記述される天体の一例。

宇宙物理学や流体力学において、レーン＝エムデン方程式（レーン＝エムデンほうていしき、Lane–Emden equation）は、球対称な密度分布を示す力学平衡にある自己重力流体を記述する微分方程式である。名称は宇宙物理学者のジョナサン・ホーマー・レーンとロバート・エムデンに由来する[1]。

目次

解説 [編集]

ポリトロピック指数n のレーン＝エムデン方程式は以下の微分方程式として表される。

\frac{1}{\xi^2} \frac{d}{d\xi} \left({\xi^2 \frac{d\theta}{d\xi}}\right) + \theta^n = 0

ここで、ξは半径r を無次元化した変数

r = \alpha \xi =\left(\frac{(n+1)P_c}{4 \pi G \rho_c^2} \right)^{\frac{1}{2}} \xi

であり、θは密度ρを無次元化した変数

\rho = \rho_c \theta^n \,

である。ただし、G は万有引力定数、Pc は球対称な流体の中心圧力、ρc は中心密度であり、ポリトロピック指数n は圧力と密度の関係式（ポリトロピックな関係式）

P = K \rho^{1 + \frac{1}{n}}

を満たす（ただし、K は定数）。この式を満たす球対称な流体はポリトロープと呼ばれる。

境界条件 [編集]

この方程式は2階の常微分方程式であるから、一意的な解を求めるためには以下の2つの境界条件が必要である。

  • \theta(0)=1 \,
  • \left( \frac{d \theta}{d \xi} \right)_{\xi=0} = 0 \,

第1式は球体の中心（r=0すなわちξ=0）における密度が有限の値ρcを持つことを意味している。第2式は球体の中心で重力がゼロ（M→0）になるのと同時に圧力勾配もゼロ(dP/dr=0)となり、さらに圧力と密度はポリトロピックな関係によって結ばれているので、密度勾配もゼロとなることを意味している。

方程式の解 [編集]

n = 0, 1, 2, 3, 4, 5, 6に対する解。横軸はξ、縦軸はθ。

レーン＝エムデン方程式はn=0,1,5の場合にのみ、解析的に解くことが可能である。その他のnに対する解は数値計算によって求められる。n=0,1,5に対する解は以下のようになる。

n = 0 1 5
\theta = 1 - \frac {\xi^2}{6} \frac{\sin\xi}{\xi} \left(1+ \frac{\xi^2}{3}\right)^{-\frac{1}{2}}
\xi_1 = \sqrt 6 \pi ∞

ここで、ξ1はθ=0となるときのξである。この値は物理的に重要で、θ=0が成り立つとき、圧力と密度もゼロ（P=ρ=0）となるので、この位置を星や流体の表面であると考えれば、ξ1を用いて中心からの半径を求めることができる。

物理現象への適用例 [編集]

レーン＝エムデン方程式を実際の物理現象に適用する例として、ポリトロープと見做せる球対称な星の半径と質量の導出を解説する。

星の半径 [編集]

ポリトロープの半径Rは、ξの定義式へθ= 0となるときのξ、すなわちξ1を代入すれば求められる。

R = \alpha \xi_1 =\left(\frac{(n+1)P_c}{4 \pi G \rho_c^2} \right)^{\frac{1}{2}} \xi_1

星の質量 [編集]

ポリトロープの質量Mは、半径rの関数としての密度ρ(r)を空間積分した後、r →αξと置き換えて計算すればよい。

M = \int_0^{R}4\pi\rho r^2 \ dr = \int_0^{\xi_1}4\pi\rho\alpha^3 \xi^2 \ d\xi = 4\pi\rho_c \alpha^3 \int_0^{\xi_1} \theta^n \xi^2 \ d\xi

ここで、レーン＝エムデン方程式を代入し、θnを書き換えると、

M = 4\pi\rho_c \alpha^3 \int_0^{\xi_1} \left\{ -\frac{1}{\xi^2} \frac{d}{d\xi} \left({\xi^2 \frac{d\theta}{d\xi}}\right) \right\} \xi^2 \ d\xi = 4\pi\rho_c \alpha^3 \left\{ -\xi_1^2 \left( \frac{d\theta}{d\xi} \right)_{\xi=\xi_1} \right\}

となる。さらに、αの定義を用いれば、

M = \left(\frac{(n+1)^3P_c^3}{4 \pi G^3 \rho_c^4} \right)^{\frac{1}{2}} \left\{ -\xi_1^2 \left( \frac{d\theta}{d\xi} \right)_{\xi=\xi_1} \right\}

となる。

この表記から、n=5の場合はξ1→∞となるが、質量自体は有限の値をとることが分かる。

参照 [編集]

  1. ^ Lane, Jonathan Homer (1870), “On the Theoretical Temperature of the Sun under the Hypothesis of a Gaseous Mass Maintaining its Volume by its Internal Heat and Depending on the Laws of Gases Known to Terrestrial Experiment”, The American Journal of Science and Arts, 2nd series 50: 57–74 

関連項目 [編集]

外部リンク [編集]

「http://ja.wikipedia.org/w/index.php?title=レーン＝エムデン方程式&oldid=46828886」から取得