Application Center - Maplesoft

App Preview:

Reimer's Mass-Loss Rate

You can switch back to the summary page by clicking here.

Learn about Maple
Download Application




Reimer's Mass-Loss Rate **

restart

Problem: Reimers' (1977) mass-loss rate formula is an empirical formula that estimates the rate of mass loss for stars on the asymptotic giant branch (AGB). Use Reimers' formula to estimate the time that it would take a one-solar-mass red-giant star to be reduced to the mass of its degenerate carbon-oxygen core, approximately 0.6 of its original mass.

``

Hints:

 

The mass of the carbon-oxygen core of a one-solar-mass red-giant star is approximate 0.6 of a solar mass.

Use the luminosity formula to calculate the radus of the red giant.

Multiply Reimers' mass-loss formula by M, making the assumption that L, R, and η do not vary significantly with time.

Integrate, and rearrange to solve for M.

Plot M as a function of t (years) for 1 to 600,000.

Calculate the number of years required for the mass to fall to a value of 0.6 solar mass.

 

 

Data:

 

L := 7000

7000

(1)

ls := 3.845*10^26

0.3845000000e27

(2)

radsun := 6.95508*10^8

695508000.0

(3)

T := 3000

3000

(4)

M[0] := 1

1

(5)

eta := 1

1

(6)

sigma := 5.670367*10^(-8)

0.5670367000e-7

(7)

``

 

Useful Equations:

 

NULL

``

``

NULL 

 

 

Solution: Solve the luminosity equation for r to obtain:

 

solve(L*ls = 4*Pi*r^2*sigma*T^4, r)

-0.2159455040e12, 0.2159455040e12

(8)

``

and divide by the solar radius to obtain the radius of the red giant in units of the solar radius.

 

rsu := 2.159455040*10^11/radsun

310.4860103

(9)

 

Multiply Reimers' mass-loss formula by M, making the assumption that L, R, and η do not vary with time. (They do vary with time, but this will give an approximate result.)

 

M*diff(M, t) = -4*10^(-13)*eta*L*rsu

``

int(M, M = M[0] .. M) = int(-4*10^(-13)*eta*L*rsu, t = 0 .. t)

(1/2)*M^2-1/2 = -0.8693608288e-6*t

(10)

M = (M[0]^2-(2*4)*10^(-13)*eta*L*rsu*t)^(1/2)``

M = (1-0.1738721658e-5*t)^(1/2)

(11)

"(->)"

M = (1.-0.17387e-5*t)^(1/2)

(12)

``

Plot the function:

``

plot(sqrt(-0.1738700000e-5*t+1), t = 1 .. 600000, title = "Reimers' Mass-Loss Rate", titlefont = ["ARIAL", 15], labels = ["time in years", "fraction of original mass"], labeldirections = ["horizontal", "vertical"])

``

To calculate the time required for a 1 solar mass star to be reduced to the mass of the degenerate carbon-oxygen core (0.6 solar masses):

 

restart

``

M := .6

``

solve(M = sqrt(-0.17387000000e-5*t+1), t)

368091.1025

(13)

``

It would require about 368,000 years.

 

------------------------------------------------------------------

Reference

 

Reimers, D. (1977). On the absolute scale of mass-loss in red giants. A&A, 61, 217-224.

``

``