Yıldızların İç Yapısı ve Evrimi

(1)

(2)

Vogt-Russell Teoremi

Yıldız yapısı belirlenirken önce kütlesini ve

kimyasal bileşimini belirleriz. Bundan

sonra 3 yardımcı denklemi ve daha sonra

4 diferansiyel denklemi çözeriz.

Vogt-Russell teoremi eğer yıldızın kimyasal

bileşimi ve kütlesi belli ise bu denklem

sisteminin yalnız bir çözümü vardır der.

Teoremin ispatı yapımamıştır ve bugün

(3)

Boyutsuz Değişkenler

0 0 0 0 0 0 r r R M r r m m R R R M     = → = _ _ = _ _    

Temel yapı denklemlerini çözerek her şeyini bildiğimiz referans yıldızı ₀ indisi ile modelini yapacağımız yıldızı da indissiz gösterelim.

0 0 0 0

dr R dm M

ve

dr = R dm = M

Her iki tarafın diferansiyelini alırsak

Bu ifadeler ile yıldızın içinde yoğunluğun nasıl dağıldığını bulabiliriz.

(4)

Boyutsuz Değişkenler

2 1 4 dr dm = π ρr ₀ ₀ dm M dm = M 0 0 dr R dr = R

Kütle denklemini yazalım ve her iki tarafını son bulduğumuz denklemler ile çarpalım

(5)

Politrop

Bu şekilde her parametre için boyutsuz

değişkenler tanımlayarak ve bazı yaklaşımlar yaparak yıldız yapı denklemlerini çözmek

olasıdır. Ama ne olursa olsun bilgisayarda sayısal olarak çözmek en doğrusu.

Bu 7 yapı denklemini doğru olarak çözmek birçok gökbilimciyi uğraştırmıştır.

Yöntemlerden birini de politrop modeller oluşturmaktır. Daha önce konveksiyon

(6)

Politrop

NkT V P_g = Pg = nkT V N n ≡

İdeal gaz yasasını biliyoruz.

Burada N, parçacık sayısı, n ise parçacık yoğunluğudur. Astrofizikte gaz yasası bu şekilde kullanılmaz, kütle yoğunluğu

kullanılır. Farklı kütlelere sahip parçacıklar varsa o zaman n yerine,

(7)

Politrop

Adyabatik harekette gazın basıncı ile hacmı arasında şöyle bir bağıntı vardı.

Sabit

PV

γ

=

P =_γ Sabit

ρ veya

Burada γ, iki esas özgül ısının oranı olup, adı adyabatik ölçektir. Gaz basıncının bu şekilde ifadesine politropik durum denklemi denir.

V P C C =

γ

H g m T k P µ ρ

= Bu denklem _{nasıl bu hale}

P

_g

=

K

ρ

γ

(8)

Politrop

1 1

n

γ = +

γ _{adyabatik ölçek, n politropik ölçeğe şu} şekilde bağlıdır.

Hidrostatik denge denklemini anımsayalım ve her iki tarafını r2

_ρ

_{(r) ile çarpalım.}

2 ( ) ( ) ( ) dP r G m r r dr r ρ = − 2 ( ) ( ) ( ) r dP r G m r r dr ρ = −

Her iki tarafın r’ye göre diferansiyelini alalım.

(9)

Politrop

2 d r dP dm G dr ρ dr dr   = −     1 r2 d dr r2 ρ dP dr       = −4πGρ

Sağtaraftaki kütle sürekliliği ifadesini yerine koyalım ve r2_{’yi sol tarafa alalım.}

Bu denkleme dikkatle bakarsak basınç gradyenti sadece yoğunluğa bağlıdır. Bu

denklemden hareketle meşhur Lane-Emden denklemini elde edebiliriz.

(10)

Politrop

2 1 2 ( 1) 1 4 n_n n K d r d nG r dr dr ρ _ρ π _ρ −   +   _ _{ = −}   _ _   _ _   1 n n g P Kργ Kρ + = = dP _K n 1 1n d dr n dr ρ ρ + =

Şimdi gaz basıncını politropik ölçek ile

yazalım ve her iki tarafın r’ye göre türevini alalım.

(11)

Politrop

ρ_{(r) için ikinci dereceden bir diferansiyel}

denklem elde edilir. Bu denklemin çözümü iki sınır değer gerektirir.

r=R için ρ=0

r=0 için dρ/dr = 0

Bu son sınır değer merkezde net çekim kuvvetinin sıfır olmasından kaynaklanır. Dolayısıyla dP/dr =0’dan hareketle elde

(12)

(13)

Politrop

Denklemin sağ tarafı boyutsuz olduğu için sol tarafın da boyutsuz olması gerekir yani α _{uzunluk birimindedir. Şimdi boyutsuz bir} değişen daha tanımlayalım

r=αξ deferansiyeli dr = α dξ o zaman denklemimiz, 2 2 1 d d _n d d θ ξ θ ξ ξ ξ   = −    

Bu ifadeye Lane-Emden denklemi denir.

(14)

Politrop

θ_{’nın boyutsuz değişken ξ’ye göre değişimini} verir. Sınır değerlerine bakalım.

(15)

Politrop

Keyfi n değerleri için bu denklemin sayısal çözümü vardır, analitik çözümü yoktur.

P ∝ ρT ve P ∝ ργ 1/( 1) n T γ T ρ _∝ − _∝ n c

ρ ρ θ

=

Adyabatik koşullar altında ideal gaz yasası O zaman,

θ _{değişkenini anımsarsak}

Bu durumda θ sıcaklığın bir göstergesi olacak ve boyutsuz ξ’ye göre değişimi verecektir.

(16)

Politrop

Lane-Emden denkleminin belirli n değerleri için ise analitik çözümü vardır, sadece

n=0,1,5 değerleri için.

n=0 sabit yoğunluğu gösterir. n=0 için

θ

n₌₁

ve

ρ

₌

ρ

_c

θ

n _{den dolayı}

ρ

₌

ρ

c 2 2 6 C ξ θ = − + 2 2 1 1 d d d d θ ξ ξ ξ ξ   = −     2 3 d C d θ ξ ξ = − + ξ 2 1 6 ξ θ = − İntegralini alırsak

ξ

_{=0 da d}

θ

_/d

ξ

_{=0 ikinci sınır koşulu gereği}

C=0. Tekrar integralini alırsak

(17)

Politrop

n=1 için L-E denklemi şu şekli alır. 2 2 2 0 d d d d θ _{θ θ} ξ + ξ ξ + =

Bu denklemi çözmek için yeni değişken

tanımlamalıyız.

χ=θξ

_{o zaman denklem,}

(18)

Politrop

Bu denklem

χ

_=sin

ξ

_{koyarak çözülür. O}

zaman, _sin_ξ

θ

ξ

= Grafik olarak,

ξ=π

_de

θ

_{=0 bu}

ise yıldızın yüzeyi demektir. Ondan sonraki değerlerin fiziksel anlamı

(19)

Politrop

r=αξ değişken dönüşümünü anımsayalım. Burada α, 1/ 2 ₍₁ ₎ 2 ( 1) 4 n n c n K G α ρ π − +   =     1/ 2 2 4 K R G απ π π ∗ = =    

n=1 için yıldızın yarıçapı bulunmuş olur.

(20)

Politrop

L-E denkleminin n=0 çözümünde ρ=ρ_c çıkıyor bu durum ancak dünya benzeri bir gezegenin iç yapısı için kaba bir yaklaşımdır.

Son çözüm n=5 de ise

ξ

_R_{=∞ çıktığına göre}

bu yıldızın yarıçapı sonsuz demektir!. n>5 için tüm çözümlerde yarıçapın sonsuz

çıkacağı gösterilebilir. Bunun anlamı sadece

n<5 için yıldızın yüzeyi vardır. Normal

yıldızlar için sadece iki ilginç çözüm vardır

n=1.5 ve n=3 ve n’nin bu değerleri içinde

(21)

(22)

Politrop

n ξ_R 0 2.45 3.33 × 10-1 1 3.14 1.01 × 10-1 2 4.35 2.92 × 10-2 3 6.90 6.14 × 10-3 4 15.00 5.33 × 10-4 − dθ dξ       ξ =ξ_R

(23)

(24)

Politrop

L-E denklemi sadece bir parametreye yani

politropik ölçeğe yani n’ye bağlıdır. Verilen bir n için yıldızın iç yapısının modeli

oluşturulmuşsa o zaman boyutsuz iki

ölçeklendirilmiş parametre iç yapıyı fiziksel birimlerle göstermeli. Yani K’ya ρ_c’ye M’ye

R’ye bağlı olmalı. Eğer belirli bir n için L-E

(25)

Politrop

1/ 2 ₁ 2 1

(26)

Politrop

Burada N_n boyutsuz sadece n’ye bağlı. O zaman yıldızın ortalama yoğunluğunu da

hesap edebiliriz. 2 3 3 3 3 4 R ort c R M d R d _{ξ ξ} θ ρ ρ ξ π ξ ξ ₌   ≡ = _− _   1 1 4 3 burada 4 n n _n c n c c n n n ort GM P K W W N R ρ ρ π ρ + + _ _ = = _{= } _  

Ayrıca merkezi basıncı da hesaplayabiliriz