Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: daixieit

(c) P. Artymowicz, solely for use by students of ASTB23 course at UTSc.

Posted: 19 Sep, due: 28 Sep.

ASTB23 (STARS AND PLANETS) PROBLEM SET #1 .

Points in the square brackets give the relative weight with which the problems count toward the inal score for the assign- ment. If you need any physical constants you may ind them in the textbook or on internet.  Always state and properly check the physical units in your formulae *before* plugging in constants and input data.  cf.  methodology of solutions to problems in tutorial notes.

1    [25 p.] wattage of the sun, and a thought experiment

(a) Derive the mass of the sun Mc based on Earth,s orbital motion around the sun.  Assume trajectory to be circular.  Set up an equation for the linear speed Earth,s vk  (a.k.a.  keplerian speed), by equating the

centripetal acceleration vk(2)/r, where r = 1 AU = 149.5 mln km = 1 AU (as we derived in the tutorial),

with the acceleration due to gravity of the sun. To stablishvk , use the period of motion (take it to be P = 1 year) and the known orbital radius.  Gravitational constant equals G = 6.6743e-11 m3/(kg s2).  compare the result with solar mass found in wikipedia.

(b) Based on the lux of radiation falling on Earth (1360..1370 W/m2), in the tutorial we derived the luminosity (wattage) of the sun. Assume the number is Lc = 3.85 . 1026 W.

How many watts does 1 kg of the sun,s interior produce, on average, from nuclear reactions? (Use Lc and Mc.) compare the sun with a car battery, which supplies power 2kW and weighs 10 kg. Which body is more eficient and by what factor (to within an order of magnitude,i.e. a round power of 10)?

(c) Estimate to within order of magnitude the working time of each object.  Assume 2*current age of the sun, discharge time of battery follows from its capacity of 20 kWh (kilowatt-hours).  compare (in J/kg) the energy per unit mass emitted over the usable lifetime of the two objects. compare the orders of magnitude.

conduct a thought experiment, where you endow each kg of the sun with eficiency of energy release qual to the car battery (200 W). What would be the main-sequence lifetime of the sun?  Would life on Earth as we know it exist then? Would it be possible anywhere in the Solar System?

Hint: Assume the total energy emitted over lifetime remains the same, no matter how fast it is released; it only depends on the initial amount of ,fuel,.

2    [15 p.] what is the surface temperature of the sun?

Base your calculation on sun,s known luminosity Lc, and assume that the surface is at temperature T and emits as blackbody according to the law:  F = σT4, where F is the lux from the surface measured in W/m2, and σ = 5.67037e — 8 W/m2/k4  is the Stefan-Boltzmann constant.

Hint: We know the mean distance to the sun (d = 1 AU = 149.5 mln km) and its mean angular diameter on the sky: 0.5332 degrees. Derive sun,s physical radius R inkm or m, do not quote it from internet.

3    [30 p.] Radius of a famous star

The star Beta pictoris, a.k.a. β pic is, as the second letter β of Greek alphabet reveals, the second-brightest star in the southern constellation of painter,s Easel (pictor in Latin). It is famous for its dusty disk, which is an analog to the young solar system.

Find out how many times larger is the physical radius of that star (R*) than the solar radius (Rc), knowing its surface temperature Te   = 825o K (estimated from the peak wavelength of spectrum), dis- tance d* = 19.3 pc (obtained from parallax), and its apparent magnitude equal to m = 3.86 mag (making it visible in good conditions but inconspicuous to the naked eye).

Adopt the apparent magnitude of the sun equal to mc = —26.73 mag, and sun,s effective temperature Tc = 578o K.

Reminder:  m = —2.5logI/Io  is the magnitude of a star (unit is mag or magnitudo, essentially, like radian, it is 1). I is the observed brightness, which is simply the lux of radiation received on Earth. If m andI don,t carry any subscripts (perhaps the name of the ilter cutting out a part of the spectrum only), we mean that they are bolometric, which means describe a total over all wavelengths.  For a zero-magnitude star like the nearby, iducial star vega, I = Io  = 2.188 . 1o 8  w/m2, which, if ever needed, deines the constant Io.

4    [30 p.] Linear density star

consider a star with mass M and radius R. Hydrostatic equilibrium equation reads

dp/dr = — Gm(r)p(r)/r2

where m(r) is mass inside radius rand p(r) spatial density of gas at radius r.

warmup consideration:

we are going to consider a more realistic density proile next. But to see how the methods work, let,s consider irst a constant density ”star”.  Divide the star mentally into shells of thickness dr and density p(r), and assume p(r) = pc = const .

Mass function m(r) can be obtained by integration of a known density distribution p(r) over radius:

m(r) = lordm = 4πlorr2p(r)dr

(under the integral sign a math purist would replace r with some other similarly named dummy variable. I do not.) In our constant density case p(r) = p = const ., integration gives

m(r) = (4π/3)r3p .

That,s obvious to anyone who knows the volume of a sphere! (This is how you compute it using calculus). we have M m(R) = (4π/3)R3pc, from which a constant in the density law can be expressed through

global stellar parameters M and R:

pc =

M

(4π/3)R3

It,s just the mass of a sphere divided by its volume, i.e.  the mean density < p 持= Total Mass/volume. This expression applies to all stars with given their M and R.  Their central density may be much higher than mean density for realistic density proiles.

pressure P(r) can be obtained by integration of the hydrostatic equation with its minus sign in front, from a starting point to layer r.  This starting point is actually r = R, the surface, because it,s there that pressure starts at value zero and keeps building up as we go deeper.  This is an approximation for sure, because in real stars neither pressure nor density fall to exactly zero at what we consider stellar surface (photosphere). But they,re so small compared with their central and mean values that approximating them by zero at the surface should not worry us.

The central pressure, in a general case, is Pc  P(0):

P(0) = Gl0Rm(r)p(r)r-2 dr.

In the special case of constant density, the integral is very simple, and the result reduces to

P(0) =  .

we,ve gone through the dimensions of central pressure in other estimates done in lectures, so we know that this formula has right units of force divided by area.  plugging in the M and R of the sun, we obtain 1.3e+14 pa (pa is s.I. unit of pressure, 1 pa = 1 N/m2).

* * *

Your task is to take a somewhat more realistic (linear) density proile

p(r) = pc(1 -x)

where x = r/R, and R is star,s radius, and redo all the calculations.

A. start from m(r). Find central density such that m(R) = M.  show that a spherical body with linear density falloff has a central density pc  that is 4 times larger than the mean density < p 持.

B. perform calculation of P(r) function. check the units. Find Pc P(0) for this linear density model. compute the value assuming solar M and R.

c. knowing the equation of state of ideal gas, combine Pc  and pc  into an estimate of Tc.  check the units again. You may assume that mean molecular weight of solar-composition gas is μ = 1.2. Evaluate the temperature assuming solar values for mass and radius. compare the result with an estimate obtained in the textbook from the virial theorem.  By what factor is it different from the virial estimate?  Are your Tc and the texbook T exactly the same physical quantity?