Homework 16 Problems with Solution - Stellar Structure and Evolution | PHYS 132, Assignments of Physics

Material Type: Assignment; Professor: Martin; Class: STELLAR STRCTR/EVOL; Subject: Physics; University: University of California - Santa Barbara; Term: Unknown 2008;

Typology: Assignments

Pre 2010

Uploaded on 08/30/2009

koofers-user-btd-1
koofers-user-btd-1 🇺🇸

4.5

(2)

10 documents

1 / 2

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Phys 132 Homework 16 Solutions
Professor Crystal Martin
TA: Ellie Hadjiyska
Problem 16 (5.2): (a) We now consider a family of chemically homogeneous stars similar in every
way but for their masses and density profiles which go as
ρ(r) = M/R3Fρ(x)
where x=r/R. We further assume energy is generated by p-p chain, stellar material obeys ideal
gas law, and opacity κρT 3.5. From
dP
dr =ρ(r)GM(r)
r2
where
M(r) = 4πZr
0
dr r2ρ(r) = 4πMZx
0
dx x2Fρ(x)
we get
dP =M2
R4·4πG dx Fρ(x)
x2Zx
0
dx x2Fρ(x)¸
P(r) = ZP(x=r/R)
P(x=1)
dP =M2
R4"4πG Zx=1
x=r/R
dxFρ(x)
x2Zx
0
dx x2Fρ(x)#,
where the expression in the last square bracket is FP(x). From the ideal gas equation of state,
T(r) = µmp
k
P(r)
ρ(r)=M
Rhµmp
kFρ(x)1FP(x)i
where the expression in the square bracket is FT(r). We have the opacity as κ=κ0ρT 3.5and
Hrad =4
3
acT (r)3
ρ(r)
1
κ
dT
dr =4
3
ac
κ0
ρ(r)2T(r)6.5M
R2
dFT
dx
=4
3
ac
κ0
M5.5R2.5Fρ(x)2FT(x)6.5d
dxFT(x)
Then from Lrad(r) = 4πr2Hrad(r) we get
Lrad(r) = M5.5
R0.5·16π
3
ac
κ0
x2µFρ(x)2FT(x)6.5d
dxFT(x)¶¸
where the expression in the square bracket is Frad(x). Now with the nuclear fusion via the PP
chain, we derived in earlier problems the expression ²fusion =ξρT 4where ξis constant. We then
have
Lfusion =M6
R7"4πξ Zx=r/R
x=0
dx x2FT(x)4Fρ(x)2#
pf2

Partial preview of the text

Download Homework 16 Problems with Solution - Stellar Structure and Evolution | PHYS 132 and more Assignments Physics in PDF only on Docsity!

Phys 132 Homework 16 Solutions

Professor Crystal Martin TA: Ellie Hadjiyska

Problem 16 (5.2): (a) We now consider a family of chemically homogeneous stars similar in every way but for their masses and density profiles which go as

ρ(r) = M/R^3 Fρ(x)

where x = r/R. We further assume energy is generated by p-p chain, stellar material obeys ideal gas law, and opacity κ ∝ ρT −^3.^5. From

dP dr =^ −ρ(r)^

GM(r) r^2

where M(r) = 4π

∫ (^) r

0

dr r^2 ρ(r) = 4πM

∫ (^) x

0

dx x^2 Fρ(x)

we get

dP = M

2 R^4

[

− 4 πG dx

Fρ(x) x^2

∫ (^) x

0

dx x^2 Fρ(x)

]

P (r) =

∫ (^) P (x=r/R)

P (x=1)

dP = M

2 R^4

[

4 πG

∫ (^) x=

x=r/R

dx Fρ(x) x^2

∫ (^) x

0

dx x^2 Fρ(x)

]

where the expression in the last square bracket is FP (x). From the ideal gas equation of state,

T (r) = μmp k

P (r) ρ(r)

= M

R

[ (^) μmp k

Fρ(x)−^1 FP (x)

]

where the expression in the square bracket is FT (r). We have the opacity as κ = κ 0 ρT −^3.^5 and

Hrad = −

acT (r)^3 ρ(r)

κ

dT dr =^ −^

ac κ 0 ρ(r)

− (^2) T (r) 6. 5 M R^2

dFT dx = −

ac κ 0 M

  1. (^5) R− 2. (^5) Fρ(x)− (^2) FT (x) 6. 5 d dx FT^ (x)

Then from Lrad(r) = 4πr^2 Hrad(r) we get

Lrad(r) = M

  1. 5 R^0.^5

[

− 16 π 3

ac κ 0

x^2

Fρ(x)−^2 FT (x)^6.^5 d dx

FT (x)

)]

where the expression in the square bracket is Frad(x). Now with the nuclear fusion via the PP chain, we derived in earlier problems the expression ≤fusion = ξρT 4 where ξ is constant. We then have

Lfusion =

M^6

R^7

[

4 πξ

∫ (^) x=r/R

x=

dx x^2 FT (x)^4 Fρ(x)^2

]

where the expression in the square bracket is Ffusion(x).

(b) The power flows are sketched in the plot below. The masses used are in the range 0. 01 − 100 MØ. Now, a star keeps losing energy by converting the gravitational energy into radiation by contraction until nuclear fusion generates additional energy. Therefore, a star will contract until Lrad ∼ Lfusion is reached. This happens when

M^5.^5 R^0.^5

∼ M

6 R^7

⇒ R ∼ M^0.^077

and

L ∼

M^6

R^7 ∼^

M^6

(M^0.^077 )^7 ∼ M

  1. 46

(c) Recall that L ∝ R^2 T (^) eff^4. On the main sequence, the luminosity is generated by nuclear fusion, so L = Lfusion ∝ M^6 /R^7

Then

R^2 T (^) eff^4 ∝

M^6

R^7 ⇒ R ∝^ T^

  1. 058 eff

and hence L ∝ R^2 T (^) eff^4 ∝ T (^) eff^4.^12