Mixing
莱恩－埃姆登方程
n = 0, 1, 2, 3, 4, 5, 6 時的莱恩－埃姆登方程解。
莱恩－埃姆登方程（Lane–Emden equation）是天文物理中一個表現自重力位能，球對稱多方流體的無因次泊松方程。此方程式名字由來於強納生·荷馬·萊恩與羅伯特·埃姆登。此方程式的解表示了恆星在半徑 r{displaystyle r}
- 1ξ2ddξ(ξ2dθdξ)+θn=0{displaystyle {frac {1}{xi ^{2}}}{frac {d}{dxi }}left({xi ^{2}{frac {dtheta }{dxi }}}right)+theta ^{n}=0}
當
- ξ=r(4πGρc2(n+1)Pc)12{displaystyle xi =rleft({frac {4pi Grho _{c}^{2}}{(n+1)P_{c}}}right)^{frac {1}{2}}}
以及
- ρ=ρcθn{displaystyle rho =rho _{c}theta ^{n},}
下標 c 代表核心的壓力與密度。n{displaystyle n}
- P=Kρ1+1n{displaystyle P=Krho ^{1+{frac {1}{n}}},}
P{displaystyle P}
目录
1 應用
2 推導
2.1 以流體靜力平衡推導
2.2 以泊松方程推導
3 方程式解
3.1 解析解
3.1.1 n=0{displaystyle n=0}
3.1.2 n=1{displaystyle n=1}
3.1.3 n=5{displaystyle n=5}
3.2 數值解
4 同調變數
4.1 同調不變方程
4.2 拓撲結構不變的同調方程
5 相關條目
6 延伸閱讀
7 參考資料
8 外部連結
應用
在物理學上，流體靜力平衡與位能梯度、密度和壓力梯度相關，而泊松方程則可以是位能和密度的關係式。因此如果有一個方程式可以進一步指出壓力和密度如何互相反映，就可以得到一個解。以上多方氣體的特定選項在數學上陳述了這個問題，尤其是該陳述特別簡潔並推導出了莱恩－埃姆登方程。這個方程式對於恆星等自重力位能氣體球是相當有用的近似，但它的假設通常是受到限制。
推導
以流體靜力平衡推導
考慮到自重力位能、流體靜力平衡下的球對稱流體、質量守恆這些狀況，就可使用以下連續性方程式：
- dmdr=4πr2ρ{displaystyle {frac {dm}{dr}}=4pi r^{2}rho }
這裡 ρ{displaystyle rho }
- 1ρdPdr=−Gmr2{displaystyle {frac {1}{rho }}{frac {dP}{dr}}=-{frac {Gm}{r^{2}}}}
m{displaystyle m}
- ddr(1ρdPdr)=2Gmr3−Gr2dmdr=−2ρrdPdr−4πGρ{displaystyle {begin{aligned}{frac {d}{dr}}left({frac {1}{rho }}{frac {dP}{dr}}right)&={frac {2Gm}{r^{3}}}-{frac {G}{r^{2}}}{frac {dm}{dr}}\&=-{frac {2}{rho r}}{frac {dP}{dr}}-4pi Grho end{aligned}}}
這裡已經使用一個連續性方程式取代質量梯度。再將方程式兩側乘上 r2{displaystyle r^{2}}
- r2ddr(1ρdPdr)+2rρdPdr=ddr(r2ρdPdr)=−4πGr2ρ{displaystyle r^{2}{frac {d}{dr}}left({frac {1}{rho }}{frac {dP}{dr}}right)+{frac {2r}{rho }}{frac {dP}{dr}}={frac {d}{dr}}left({frac {r^{2}}{rho }}{frac {dP}{dr}}right)=-4pi Gr^{2}rho }
方程式兩側除以 r2{displaystyle r^{2}}
- 1r2ddr(r2Kρc1n(n+1)dθdr)=−4πGρcθn{displaystyle {frac {1}{r^{2}}}{frac {d}{dr}}left(r^{2}Krho _{c}^{frac {1}{n}}(n+1){frac {dtheta }{dr}}right)=-4pi Grho _{c}theta ^{n}}
將常數聚集並以 r=αξ{displaystyle r=alpha xi }
α2=(n+1)Kρc1n−1/4πG{displaystyle alpha ^{2}=(n+1)Krho _{c}^{{frac {1}{n}}-1}/4pi G},
最後得到莱恩－埃姆登方程：
- 1ξ2ddξ(ξ2dθdξ)+θn=0{displaystyle {frac {1}{xi ^{2}}}{frac {d}{dxi }}left({xi ^{2}{frac {dtheta }{dxi }}}right)+theta ^{n}=0}
以泊松方程推導
同樣地，也可以使用泊松方程進行推導：
- ∇2Φ=1r2ddr(r2dΦdr)=−4πGρ{displaystyle nabla ^{2}Phi ={frac {1}{r^{2}}}{frac {d}{dr}}left({r^{2}{frac {dPhi }{dr}}}right)=-4pi Grho }
我們可以透過以下數學公式以流體靜力平衡取代位能梯度：
- dΦdr=1ρdPdr{displaystyle {frac {dPhi }{dr}}={frac {1}{rho }}{frac {dP}{dr}}}
最後也可以得到莱恩－埃姆登方程。
方程式解
解析解
n{displaystyle n}
n=0{displaystyle n=0}
如果 n=0{displaystyle n=0}
- 1ξ2ddξ(ξ2dθdξ)+1=0{displaystyle {frac {1}{xi ^{2}}}{frac {d}{dxi }}left({xi ^{2}{frac {dtheta }{dxi }}}right)+1=0}
重新整理並進行一次積分後的公式成為：
- ξ2dθdξ=C1−13ξ3{displaystyle xi ^{2}{frac {dtheta }{dxi }}=C_{1}-{frac {1}{3}}xi ^{3}}
公式兩側都除以 ξ2{displaystyle xi ^{2}}
- θ(ξ)=C0−C1ξ−16ξ2{displaystyle theta (xi )=C_{0}-{frac {C_{1}}{xi }}-{frac {1}{6}}xi ^{2}}
邊界條件 θ(0)=1{displaystyle theta (0)=1}
n=1{displaystyle n=1}
當 n=1{displaystyle n=1}
- d2θdξ2+2ξdθdξ+θ=0{displaystyle {frac {d^{2}theta }{dxi ^{2}}}+{frac {2}{xi }}{frac {dtheta }{dxi }}+theta =0}
兩端都乘以 ξ2{displaystyle xi ^{2}}
- θ(ξ)=sinξξ{displaystyle theta (xi )={frac {sin xi }{xi }}}
n=5{displaystyle n=5}
在經過一連串取代的步驟後，方程式可以有進一步的解：
- θ(ξ)=11+ξ2/3{displaystyle theta (xi )={frac {1}{sqrt {1+xi ^{2}/3}}}}
當 n=5{displaystyle n=5}
數值解
一般情形下莱恩－埃姆登方程的解必須以數值積分方式求得。許多數值積分的標準解法要求該問題必須以一階常微分方程表示，例如：
- dθdξ=−ϕξ2{displaystyle {frac {dtheta }{dxi }}=-{frac {phi }{xi ^{2}}}}
- dϕdξ=θnξ2{displaystyle {frac {dphi }{dxi }}=theta ^{n}xi ^{2}}
在這裡 ϕ(ξ){displaystyle phi (xi )}
同調變數
同調不變方程
已知如果 θ(ξ){displaystyle theta (xi )}
而這類可選擇的變數有多個，一個適當的選擇是：
- U=dlogmdlogr=ξ3θnϕ{displaystyle U={frac {dlog m}{dlog r}}={frac {xi ^{3}theta ^{n}}{phi }}}
和
- V=dlogPdlogr=(n+1)ϕξθ{displaystyle V={frac {dlog P}{dlog r}}=(n+1){frac {phi }{xi theta }}}
我們可以將相對於 ξ{displaystyle xi }
- 1UdUdξ=1ξ(3−n(n+1)−1V−U){displaystyle {frac {1}{U}}{frac {dU}{dxi }}={frac {1}{xi }}(3-n(n+1)^{-1}V-U)}
和
1VdVdξ=1ξ(−1+U+(n+1)−1V){displaystyle {frac {1}{V}}{frac {dV}{dxi }}={frac {1}{xi }}(-1+U+(n+1)^{-1}V)}.
最後，我們將以上兩個方程式相除以消去應變量 ξ{displaystyle xi }
- dVdU=−VU(U+(n+1)−1V−1U+n(n+1)−1V−3){displaystyle {frac {dV}{dU}}=-{frac {V}{U}}left({frac {U+(n+1)^{-1}V-1}{U+n(n+1)^{-1}V-3}}right)}
以上即為單一一階方程式。
拓撲結構不變的同調方程
同調性不變的方程式可被視為自主對方程式：
- dUdlogξ=−U(U+n(n+1)−1V−3){displaystyle {frac {dU}{dlog xi }}=-U(U+n(n+1)^{-1}V-3)}
和
- dVdlogξ=V(U+(n+1)−1V−1){displaystyle {frac {dV}{dlog xi }}=V(U+(n+1)^{-1}V-1)}
這些方程式的解的形式可透過以下線性穩定性分析來決定。方程式的臨界點（當 dV/dlogξ=dU/dlogξ=0{displaystyle dV/dlog xi =dU/dlog xi =0}
| 臨界點 | 特徵值 | 特徵向量 | ||
|---|---|---|---|---|
(0,0){displaystyle (0,0)} | 3{displaystyle 3} | −1{displaystyle -1} | (1,0){displaystyle (1,0)} | (0,1){displaystyle (0,1)}  |
(3,0){displaystyle (3,0)} | −3{displaystyle -3} | 2{displaystyle 2} | (1,0){displaystyle (1,0)} | (−3n,5+5n){displaystyle (-3n,5+5n)}  |
(0,n+1){displaystyle (0,n+1)} | 1{displaystyle 1} | 3−n{displaystyle 3-n} | (0,1){displaystyle (0,1)} | (2−n,1+n){displaystyle (2-n,1+n)}  |
(n−3n−1,2n+1n−1){displaystyle left({frac {n-3}{n-1}},2{frac {n+1}{n-1}}right)} | n−5±Δn2−2n{displaystyle {frac {n-5pm Delta _{n}}{2-2n}}} | (1−n∓Δn,4+4n){displaystyle left(1-nmp Delta _{n},4+4nright)}  | ||
相關條目
- 恆星結構
- 埃姆登－錢德拉塞卡方程
- 錢德拉塞卡白矮星方程
延伸閱讀
Lane, Jonathan Homer, 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, 1870, 50: 57–74 .
參考資料
^ Chandrasekhar, Subrahmanyan. An introduction to the study of stellar structure. Chicago, Ill.: University of Chicago Press. 1939.
^ Horedt, Georg P. Topology of the Lane-Emden equation. A&A. 1987, 117 (1-2): 117–130 [2012-06-27]. 引文格式1维护：日期与年 (link)
外部連結
- 埃里克·韦斯坦因. Lane-Emden Differential Equation. MathWorld.
- Horedt, George Paul ( 1986 ) 'Seven-digit tables of Lane-Emden functions' PDF ( 5.9MB ), Astrophysics and Space Science vol. 126, no. 2, Oct. 1986, p. 357–408. ( ISSN 0004-640X ). Collected at the Smithsonian/NASA Astrophysical Data System.
Comments
Post a Comment