http://www.gipsky.com/mediawiki/index.php?title=%E6%81%86%E6%98%9F%E7%B5%90%E6%A7%8B
恆星結構方程式
最簡單和最常用的恆星結構模型是球狀對稱的準靜態模型,假設恆星是處在穩定狀態並且是球型對稱,它包含4個基本的一階偏微分方程:兩個用來描述物質和壓力隨著半徑的變化;另外兩個用來描述溫度與光度隨著半徑的變化。[4]。
在組成的恆星結構方程式 (利用假設的球狀對稱)中,要考慮的各項是物質密度 <math>\rho(r)</math>、溫度 <math>T(r)</math>、總壓力 (氣體的熱壓力加上輻射壓)<math>P(r)</math>、光度<math>l(r)</math>,和在距離中心<math>r</math>,厚度<math>\mbox{d}r</math>的球殼,每單位質量的能量孳生率<math>\epsilon(r)</math>。這顆恆星被假設局部熱力平衡(LTE),所以物質和光子的溫度相等的。雖然局部熱力平衡不是嚴格的被遵守,因為對指定的殼層之下的溫度永遠比上面的高,良好的近似範圍可以估計是光子的平均自由徑 <math>\lambda</math>,因為這個長度遠比考慮的溫度變化為小,也就是說<math>\lambda \ll T/|\nabla T|</math>。
第一種狀態是流體靜力平衡 :恆星內部向外的力量來自壓力梯度與來自萬有引力向內的力量完全平衡。
<math> {\mbox{d} P \over \mbox{d} r} = - { G m \rho \over r^2 } </math>,
此處<math>m(r)</math>是在殼層<math>r</math>內部累積的總質量,G 是萬有引力常數。累積的總質量依據 質量連續方程式隨著距離的增加而增加:
<math> {\mbox {d} m \over \mbox{d} r} = 4 \pi r^2 \rho .</math>
質量連續方程式的積分由恆心中心(<math>r=0</math>) 至恆星半徑為(<math>r=R</math>) 之處時,就可以得到恆星的總質量。
考慮能量從球殼向外傳送的是能量方程式:
<math> {\mbox{d} l \over \mbox{d} r} = 4 \pi r^2 \rho ( \epsilon - \epsilon_\nu )</math>,
此處<math>\epsilon_\nu</math>是單位質量以微中子的形式產生的光度(通常不會與一般的物質產生交互作用而逃離恆星)。在恆星的核心之外,核反應不會發生,也就是沒有能量產生,因此發光度是一個常數。
能量轉移方程式決於能量轉移的模型而有不同的形式。以傳導性的光度轉移(適用於白矮星),能量方程式是:
<math> {\mbox{d} T \over \mbox{d} r} = - {1 \over k} { l \over 4 \pi r^2 },</math>
此處k是 導熱率。
在能量以輻射傳送的情況下,這適用於主序帶上與太陽類似,核心外有殼層包圍著的大質量主序星的恆星內部:
<math> {\mbox{d} T \over \mbox{d} r} = - {3 \kappa \rho l \over 64 \pi r^2 \sigma T^3},</math>
此處<math>\kappa</math>是物質的不透明性,<math>\sigma</math>是 斯特凡-波茲曼常數,並且波茲曼常數被設定為1。
在能量以對流傳送的情況下(適用於所有主序帶上沒有輻射部分的巨星和低質量恆星),已知可用數學公式還沒有一個是嚴謹的。對流的能量傳送通常都使用混合長度理論,混合長度理論對恆星氣體的處理是認為分離的各種不同元素,保留了周圍的溫度、密度和壓力,但是在恆星內能移動的距離就是特徵長度,稱為混合長度[5]。對單原子 理想氣體,混合長度理論認為:
<math> {\mbox{d} T \over \mbox{d} r} = \left(1 - {1 \over \gamma} \right) {T \over P } { \mbox{d} P \over \mbox{d} r},</math>
此處<math>\gamma = c_p / c_v</math>是絕熱指數,在氣體就是比熱的比率(對完全電離的理想氣體,<math>\gamma = 5/3</math>.)。
還需要狀態方程式,關係到其他區域的壓力變化,適用於材料,像是溫度、密度、化學組成等等。相關的狀態也許必須包括理想氣體定律、輻射壓、來自電子簡併的壓力等等。
與一系列的邊界條件組合,這些方程式的一個解可以完整的敘述恆星的行為。典型的邊界條件設由觀測的參數來設定數值,適當的恆星表面(<math>r=R</math>)和中心(<math>r=0</math>):<math>P(R) = 0</math>,意思是恆星表面的壓力為0;<math>m(0) = 0</math>為在恆星的中心沒有質量;如同必須假設密度是有限的,<math>m(R) = M</math>就是恆星的總質量;而<math>T(R) = T_{eff}</math>是恆星表面的有效溫度。
雖然現在的恆星演化模型能描述赫羅圖的主要特徵,重要的改進已經能將與我們再傳送現象的知識及現下仍不確定的因素移除,但最困難的挑戰依然是擾動的數值作業。一些研究小組正在開發的模型可以演算在三度空間的簡易擾動。
|