Physics 556 Stellar Astrophysics Prof. James Buckley
|
|
- Everett Marshall
- 5 years ago
- Views:
Transcription
1 hysics 556 Stellar Astrophysics rof. James Buckley Lecture 8 Convection and the Lane Emden Equations for Stellar Structure
2 Reading/Homework Assignment Read sections 2.5 to 2.9 in Rose over spring break!
3 Stellar Structure
4 Equations of Stellar Structure Hyostatic Equilibrium : d Mass Continuity : dm(r) (r) =4πr 2 ρ(r) = ρ(r) GM(r) r 2 Eqn. of State : gas = nkt(r) = ρ(r) µm u kt(r) rad (r) = 1 at (r)4 3 Radiative transport equation : T 3 κρ L r = 16πac r 2 T 3 Energy dl(r) roduction/conservation= T γ 4πr 2 ρ(r)ɛ(r) Convective Energy Transport: dt d = γ 1
5 Convection ρ i, i,t i ρ o (r + r), o (r + r), T o (r + r) i = o,t i >T o ρ i <ρ o ρ i, i,t i ρo(r + r), o(r + r), To(r + r) ρ i, i,t i ρ o (r), o (r), T o (r) The gas bubble in the star moves fast enough that there is not time to exchange heat at the boundary adiabatic expansion
6 Solar Granulation 30 minute time lapse images of the solar granulation made with the SVST (La alma) (G. Scharmer, et al.). Hot gas rises in center of bright regions, is diverted horizontally (radiates) and sinks back into the Sun in the darker intergranular lanes. The size of the granules range from 250 km (resolution limit) to 2000~km with an average of 1300 km.
7 Adiabatic Expansion Specific heat or heat capacity is defined as the amount of heat that must be added to a substance to raise its temperature by 1 deg K either at constant pressure C or constant temperature C T Enthalpy : H U + V dh = du +(dv + V d ) =(δq dv )+(dv + V d )=δq + Vd δq H C p = dt T Using the ideal gas law V = NkT du = δq dv δq U C v = dt T du = αnk dt C v = αnk d(v)=nkdt du = αnkdt = αd(v)=α( dv + V d ) dh =(α + 1)( dv + V d ) V Nk C p =(α + 1) =(α + 1) T The adiabatic index γ is defined to be: V V and U = α NkT for 2α D.O.F. =(α + 1)Nk γ C /C V =(α + 1)/α
8 Adiabatic Expansion d(v)=nkdt du = αnkdt = αd(v)=α(dv + V d ) For adiabatic expansion, δq =0=dU + dv, substituting for du gives: α( dv + V d )= dv (α + 1) dv = αv d ln b a = α +1 α (α + 1) dv V ln Vb V a = αd b a = Vb V a α+1 α V γ = constant where γ α +1 α = C p C v V γ = const γ mn = const ρ ρ γ = const Substituting = ρ m kt gives Tγ γ = const 1 γ T γ = const
9 Convection r + δr ρ i (r + δr) i (r + δr) T i (r + δr) ρ o (r + δr) o (r + δr) T o (r + δr) Assume no heat flows across the boundary adiabatic expansion Assume that pressure balance is always established by free expansion/contraction of the element o (r) = i (r) Start by considering the adiabatic expansion of the element: o (r + δr)ρ i (r + δr) γ = o (r)ρ i (r) γ (Note i = o ) r ρ i, i,t i d( 0 ρ γ i )=d(const) = 0 d 0 ρ γ i dρ i d 0 = ρ i γ 0 dρ i = dρ i d 0 d 0 = ρ i(r) γ o (r) γ 0 ρ γ 1 i dρ i =0 d o Now, we compare this with the outside density gradient to see if the object continues to be boyant (convective instability) or sinks back down (convective stability). Stability condition : dρ i < dρ 0 dρ i < dρ 0 ρ i(r) d o γ o (r) > dρ o to first order, we can assume that ρ i (r) ρ o (r), substitute, divide both sides by d o / to obtain o dρo Stability condition: < 1 d o γ ρ o
10 Convection For a stellar atmosphere where radiation pressure is negligible, the equation of state is given by the ideal gas law: Differentiating both sides gives = ρ m kt ln = ln ρ + ln T + constant d Starting with the expression for stability = dρ ρ + dt T o ρ o dρo d o < 1 γ or dρ o ρ o 1 γ d o o < do γ o Stability condition : T dt o T o < 0 dt d > γ 1 γ
11 Stability Condition T Stability condition : T dt/ d/ > γ 1 γ Since dt/ = dt/ and d/ = d/ dt γ 1 T < d dt γ 1 T γ < γ Hyostatic Equilibrium : defining g GM r 2 dt d > γ 1 γ dt γ 1 T > γ d (r) V = NkT = ρkt m = ρ(r) GM(r) r 2 d = ρ g d substituting for d/ and dt γ 1 < T (ρ g) γ ρkt/m Condition for stability : dt γ 1 gm < γ k d
12 Convective Transport γ 1 γ Stability condition : T = 1+1/α 1 (α + 1)/α = 1 α +1 dt d > γ 1 γ dt or < dt γ 1 gm < γ k gmn (α + 1)N k = gmn C p = gc p (In chemistry, specific heat c p is defined as heat capacity per unit mass C p /m N) Which matches the stability condition given in some textbooks (e.g., Rose): Stability condition : dt < g c p The δt and velocity (from the boyant force) of the rising elements are determined by the difference between the temperature gradient (outside) and the adiabatic gradient (inside), any difference will result in convective transport that tends to reduce the outside temperature gradient. So, to a reasonable degree of accuracy we can assume that in a convective zone, the temperature gradient takes on almost exactly the adiabatic value, turning our inequality into an equality: Giving us the equation for convective transport T dt d = γ 1 γ
13 Mixed Convection and Radiation Let s assume that the atmosphere starts out in radiative equilibrium, i.e F rad = σt 4 everywhere and 4πr 2 F rad (r) =const = L where L is the total luminosity of the star. Together with the other equations of stellar structure, this determines T (r), (r) and ρ(r) Now say that some distance from the center of the star r 0 the atmosphere becomes unstable to convection - we need to calculate F (r) =F rad + F conv, solve for energy balance, solve for T (r), (r) and ρ(r), check the conditions for convective stability again, and iterate. T i T l m At each step, we can start with the external values T,,ρ and the internal values T i,,ρ i and calculate quantities like: ρ(r) =ρ i (r) ρ(r), T (r) =T i (r) T (r), ρ i ρ Then the flux = mass flow rate excess heat per gram: F conv = ρ v c p T r r 0 To find the velocity of the element v we assume that the convective cell rises through a distance l m (the mixing length) where it looses its identity. The average velocity is found by equating the kinetic energy at r = r 0 + l m with the work done by the boyant force acting through l m on an element with under-density ρ
14 Convective Energy Flux Flux = mass flow rate excess heat per gram: F conv = ρ v c p T Recall that for adiabatic expansion 1 γ T γ i = const T i ρ i ρ T l m d( 1 γ T γ i ) = (1 γ) γ T γ i multiplying by γ T γ i d + γ 1 γ T γ 1 i dt i =0 (1 γ)d + γ dt i =0 d ln T i T i d ln = γ 1 γ r r 0 dt i = γ 1 γ T d The difference in temperature of the element and the surroundings is given by: T dt i δr dt δr = γ 1 γ T d dt δr T δr The excess internal energy per unit mass carried by the gas is u = c p T And the convective energy flux is F c = L c 4πr 2 = c p T δr ρv c u mass flux
15 Convective Energy Flux Now let s find v using the conversion of work to kinetic energy First, lets scale l m by some characteristic scale. The scale height λ of the stellar atmosphere is defined to be the height at which the pressure ops by 1/e Assume that convective elements are constantly being born out of small temperature perturbations with initial velocity v c = 0 at some distance r 0 from the center of the star (the position at which the atmosphere becomes convectively unstable. The elements pick up (or loose) speed as they rise (or sink) until after moving through a distance l m they mix with the surroundings and loose their identity. We can write l m = αλ where α 0.5 to 1.5 The change in kinetic energy of the rising element (per unit volume) is kinetic energy density = 1 2 ρv2 c The boyant force does work per unit volume of dw = ρ g δr in moving through a distance δr. Integrating this and equating this to the change in kinetic energy density gives 1 2 ρv2 c = λ r0 +l m r 0 d/ ρ GM(r) r 2
16 Convective Energy Flux To solve for ρ we use the equation of state: = ρ kt µm u d = dρ ρ + dt T dµ µ (where, to be completely general, we consider the case that even the ionization state or composition of the gas may be different inside and outside of the convective element, i.e. µ = 0) If we have pressure balance between the inside and outside of the gas element, then d =0=dρ ρ + dt T dµ µ dt dρ = ρ T dµ = ρ dt µ T η where η 1 d ln µ d ln T 1 2 ρv2 c = r0 +l m r 0 ρ GM(r) r 2 ρ = η ρ T T 1 2 η ρ dti T dt r0 +l m l m r 0 v = ηg T ( T )1/2 l m g η 2 g ρ T T l2 m
17 Convective Energy Flux substituting v = ηg T ( T )1/2 l m into F conv = L c 4πr 2 = c p T δr ρv c and taking δr l m F conv = ρ c p η g T ( T )3/2 l 2 m Finally, we demand that π(f rad + F conv )=σt 4 and iterate.
18 Convective flux F c = L c 4πr 2 = c p T δr ρv c u mass flux v c = 1 ρ dρ 1 γ 1 d 1/2 GM(r) r 2 1/2 l m 2 and we can show this can be rewritten as v c =( T ) 1/2 GM(r) Tr 2 1/2 l m 2 Combining these equations, and dividing by two (to account for both rising and falling convective blobs???) F c = L c 4πr 2 = c pρ( T ) 3/2 GM(r) Tr 2 1/2 l 2 m 4
19 olytropes In general we have a system of equations to solve for 3 variables (the temperature, pressure and density gradients). Solving this equation depends on the details of the equation of state, opacity and other physics that may change throughout the star. If we have a relationship that instead gave a relationship between pressure and density = (ρ) (with no explicit temperature dependence) we can solve for the structure of the star using the equations of hyostatic equilibrium and mass conservation alone. We call hypothetical stars with a simple relationship between pressure and density polytropes, where = Kρ 1+1/n The equations of hyostatic equilibrium and mass continuity are (i) (ii) d = ρgm(r) r 2 dm(r) = 4πr 2 ρ combining (i) and (ii) we obtain 4πr 2 ρ = d r2 ρg d
20 olytropes d r 2 ρ d r 2 = 4πGρ substituting ρ(r) = ρ c θ n (r) = Kρ (1+1/n) c θ (n+1) d = Kρ(1+1/n) c (n + 1)θ n (r) dθ Define a new dimensionless variable corresponding to the radial distance from the center of the star: r = aξ r d 2 d ρ d a 2 ξ 2 ρ c θ Kρ (1+1/n) r 2 = n c (n + 1)θ n dθ adξ a 3 ξ 2 dξ if we choose a = (n + 1)Kρ (1/n 1) c 4πG = 4πGρ c θ n 1/2 Lane-Emden equation of index n 1 ξ 2 d dθ ξ2 dξ dξ = θn
21 olytropes Analytic solutions of the Lane-Emden equation exist only for a few values of n ξ 2 n = 0, θ = 1 6 n = 1, θ = sin ξ ξ 1 n = 5, θ = ξ2 0.5 For other values of n, we can solve the equations numerically, starting with the boundary conditions ρ = ρ c at r = 0 so θ(0) = 1 d = 0 at r = 0 dθ dξ = 0 ξ=0 Then, rewrite the Lane-Emden equation 1 ξ 2 d dθ ξ2 dξ dξ = θn 2ξ dθ dξ + ξ2 d2 θ dξ + θn ξ 2 = 0 d 2 θ dξ 2 = 2 ξ dθ dξ θn
22 olytropes d 2 θ dξ 2 = 2 ξ dθ dξ θn Now convert this differential equation to a finite difference equation dθ dξ 1 ξ (θ i+1 θ i ) d 2 θ dξ 2 1 ξ 2 (θ i+1 2θ i + θ i 1 ) ξ i ξ θ n θ n i Then, starting at the middle of the star with i = 0, θ = 1, ξ = 0 and dθ/dξ = 0, solve for θ i+1 and repeat until the surface is reached (when θ ops below zero, and hence density ops below zero) The surface of the star R = aξ s can be found by finding the root, θ(ξ s )=0
23 Numerical Solutions θ(ξ) ξ Hint: Start integration at ξ =, i.e., some small positive number not exactly equal to zero!
24 olytropes We can derive some useful physical quantities from the solution of the Lane- Emden equation: Stellar radius : R = aξ s = (n + 1)K 4πG 1/2 ρ (1 n)/2n c ξ s Total mass : M =4πa 3 ρ c ξs 0 θ n ξ 2 dξ using 1 ξ 2 d dθ ξ2 dξ dξ = θn we see that M =4πa 3 ρ c ξs o d ξ 2 dθ dξ and thus M =4π (n + 1)K 4πG 3/2 ρ (3 n)/2n c ξ 2 dθ dξ ξ s
Physics 576 Stellar Astrophysics Prof. James Buckley. Lecture 10
Physics 576 Stellar Astrophysics Prof. James Buckley Lecture 10 Reading/Homework Assignment Makeup class tomorrow morning at 9:00AM Read sections 2.5 to 2.9 in Rose over spring break! Stellar Structure
More informationPhysics 576 Stellar Astrophysics Prof. James Buckley. Lecture 11 Polytropes and the Lane Emden Equation
Physics 576 Stellar Astrophysics Prof. James Buckley Lecture 11 Polytropes and the Lane Emden Equation Reading/Homework Assignment Makeup class tomorrow morning at 9:00AM Read sections 2.5 to 2.9 in Rose
More informationPhysics 556 Stellar Astrophysics Prof. James Buckley. Lecture 9 Energy Production and Scaling Laws
Physics 556 Stellar Astrophysics Prof. James Buckley Lecture 9 Energy Production and Scaling Laws Equations of Stellar Structure Hydrostatic Equilibrium : dp Mass Continuity : dm(r) dr (r) dr =4πr 2 ρ(r)
More informationSummary of stellar equations
Chapter 8 Summary of stellar equations Two equations governing hydrostatic equilibrium, dm dr = 4πr2 ρ(r) Mass conservation dp dr = Gm(r) r 2 ρ Hydrostatic equilibrium, three equations for luminosity and
More information3 Hydrostatic Equilibrium
3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose
More informationSPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours
MSc/MSci Examination Day 28th April 2015 18:30 21:00 SPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours YOU ARE NOT PERMITTED TO READ THE CONTENTS OF THIS QUESTION PAPER UNTIL
More informationASTM109 Stellar Structure and Evolution Duration: 2.5 hours
MSc Examination Day 15th May 2014 14:30 17:00 ASTM109 Stellar Structure and Evolution Duration: 2.5 hours YOU ARE NOT PERMITTED TO READ THE CONTENTS OF THIS QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY
More informationEnergy transport: convection
Outline Introduction: Modern astronomy and the power of quantitative spectroscopy Basic assumptions for classic stellar atmospheres: geometry, hydrostatic equilibrium, conservation of momentum-mass-energy,
More informationStatic Stellar Structure
Most of the Life of A Star is Spent in Equilibrium Evolutionary Changes are generally slow and can usually be handled in a quasistationary manner We generally assume: Hydrostatic Equilibrium Thermodynamic
More informationConvection. If luminosity is transported by radiation, then it must obey
Convection If luminosity is transported by radiation, then it must obey L r = 16πacr 2 T 3 3ρκ R In a steady state, the energy transported per time at radius r must be equal to the energy generation rate
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationCHAPTER 16. Hydrostatic Equilibrium & Stellar Structure
CHAPTER 16 Hydrostatic Equilibrium & Stellar Structure Hydrostatic Equilibrium: A fluid is said to be in hydrostatic equilibrium (HE) when it is at rest. This occurs when external forces such as gravity
More information2. Equations of Stellar Structure
2. Equations of Stellar Structure We already discussed that the structure of stars is basically governed by three simple laws, namely hyostatic equilibrium, energy transport and energy generation. In this
More informationASTR 1040 Recitation: Stellar Structure
ASTR 1040 Recitation: Stellar Structure Ryan Orvedahl Department of Astrophysical and Planetary Sciences University of Colorado at Boulder Boulder, CO 80309 ryan.orvedahl@colorado.edu February 12, 2014
More informationConvection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3
Convection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3 4ac kr L T 3 4pr 2 Large luminosity and / or a large opacity k implies
More informationHandout 18: Convection
Handout 18: Convection The energy generation rate in the center of the sun is ~ 16 erg/gm/s Homework, sum of p-p and CNO production The average energy generation needed to supply the sun s luminosity is
More informationTides in Higher-Dimensional Newtonian Gravity
Tides in Higher-Dimensional Newtonian Gravity Philippe Landry Department of Physics University of Guelph 23 rd Midwest Relativity Meeting October 25, 2013 Tides: A Familiar Example Gravitational interactions
More informationConvection and Other Instabilities. Prialnik 6 Glatzmair and Krumholz 11 Pols 3, 10
Convection and Other Instabilities Prialnik 6 Glatzmair and Krumholz 11 Pols 3, 10 What is convection? A highly efficient way of transporting energy in the stellar interior Highly idealized schematic.
More information9.1 Introduction. 9.2 Static Models STELLAR MODELS
M. Pettini: Structure and Evolution of Stars Lecture 9 STELLAR MODELS 9.1 Introduction Stars are complex physical systems, but not too complex to be modelled numerically and, with some simplifying assumptions,
More informationAstronomy 112: The Physics of Stars. Class 9 Notes: Polytropes
Astronomy 112: The Physics of Stars Class 9 Notes: Polytropes With our discussion of nuclear reaction rates last time, we have mostly completed our survey of the microphysical properties of stellar matter
More informationdp dr = ρgm r 2 (1) dm dt dr = 3
Week 5: Simple equilibrium configurations for stars Follows on closely from chapter 5 of Prialnik. Having understood gas, radiation and nuclear physics at the level necessary to have an understanding of
More informationAy123 Set 1 solutions
Ay13 Set 1 solutions Mia de los Reyes October 18 1. The scale of the Sun a Using the angular radius of the Sun and the radiant flux received at the top of the Earth s atmosphere, calculate the effective
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More information2. Basic assumptions for stellar atmospheres
. Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres
More informationPHYS 633: Introduction to Stellar Astrophysics
PHYS 6: Introduction to Stellar Astrophysics Spring Semester 2006 Rich Townsend (rhdt@bartol.udel.edu Convective nergy Transport From our derivation of the radiative diffusion equation, we find that in
More information= m H 2. The mean mass is 1/2 the mass of the hydrogen atom. Since, the m e << m p, the mean mass of hydrogen is close to 1/2 the mass of a proton.
From Last Time: Mean mass m m = m e + m p 2 = m H 2 The mean mass is 1/2 the mass of the hydrogen atom. Since, the m e
More informationdr + Gm dr Gm2 2πr 5 = Gm2 2πr 5 < 0 (2) Q(r 0) P c, Q(r R) GM 2 /8πR 4. M and R are total mass and radius. Central pressure P c (Milne inequality):
1 Stellar Structure Hydrostatic Equilibrium Spherically symmetric Newtonian equation of hydrostatics: dp/dr = Gmρ/r 2, dm/dr = 4πρr 2. 1) mr) is mass enclosed within radius r. Conditions at stellar centers
More information2. Basic assumptions for stellar atmospheres
. Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres
More informationPolytropic Stars. c 2
PH217: Aug-Dec 23 1 Polytropic Stars Stars are self gravitating globes of gas in kept in hyostatic equilibrium by internal pressure support. The hyostatic equilibrium condition, as mentioned earlier, is:
More informationEquations of Stellar Structure
Equations of Stellar Structure Stellar structure and evolution can be calculated via a series of differential equations involving mass, pressure, temperature, and density. For simplicity, we will assume
More informationThe Equations of Stellar structure
PART IV: STARS I Section 1 The Equations of Stellar structure 1.1 Hyostatic Equilibrium da P P + dp r Consider a volume element, density ρ, thickness and area da, at distance r from a central point. The
More informationOrdinary Differential Equations
Ordinary Differential Equations Part 4 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Ordinary Differential Equations p. 1/23 Two-point Boundary Value Problems NRiC 17. BCs specified at two
More informationM J v2 s. Gm (r) P dm = 3. For an isothermal gas, using the fact that m r, this is GM 2 R = 3N ot M 4πR 3 P o, P o = 3N ot M
1 Stellar Birth Criterion for stability neglect magnetic fields. Critical mass is called the Jean s mass M J, where a cloud s potential energy equals its kinetic energy. If spherical with uniform temperature,
More informationThe Interiors of the Stars
The Interiors of the Stars Hydrostatic Equilibrium Stellar interiors, to a good first approximation, may be understood using basic physics. The fundamental operating assumption here is that the star is
More informationTRANSFER OF RADIATION
TRANSFER OF RADIATION Under LTE Local Thermodynamic Equilibrium) condition radiation has a Planck black body) distribution. Radiation energy density is given as U r,ν = 8πh c 3 ν 3, LTE), tr.1) e hν/kt
More informationModel Atmospheres. Model Atmosphere Assumptions
Model Atmospheres Problem: Construct a numerical model of the atmosphere to estimate (a) Variation of physical variables (T, P) with depth (b) Emergent spectrum in continuum and lines Compare calculated
More informationChapter 16. White Dwarfs
Chapter 16 White Dwarfs The end product of stellar evolution depends on the mass of the initial configuration. Observational data and theoretical calculations indicate that stars with mass M < 4M after
More informationFree-Fall Timescale of Sun
Free-Fall Timescale of Sun Free-fall timescale: The time it would take a star (or cloud) to collapse to a point if there was no outward pressure to counteract gravity. We can calculate the free-fall timescale
More informationρ dr = GM r r 2 Let s assume that µ(r) is known, either as an arbitrary input state or on the basis of pre-computed evolutionary history.
Stellar Interiors Modeling This is where we take everything we ve covered so far and fold it into models of spherically symmetric, time-steady stellar structure. (We ll emphasize insight over rigorous
More informationLinear Theory of Stellar Pulsation
Linear Theory of Stellar Pulsation Mir Emad Aghili Department of Physics and Astronomy, University of Mississippi, University, MS, 38677 Introduction The first pulsating star to be discovered was omicron
More informationThe structure and evolution of stars. Introduction and recap
The structure and evolution of stars Lecture 3: The equations of stellar structure 1 Introduction and recap For our stars which are isolated, static, and spherically symmetric there are four basic equations
More informationFrom Last Time: We can more generally write the number densities of H, He and metals.
From Last Time: We can more generally write the number densities of H, He and metals. n H = Xρ m H,n He = Y ρ 4m H, n A = Z Aρ Am H, How many particles results from the complete ionization of hydrogen?
More informationPhases of Stellar Evolution
Phases of Stellar Evolution Phases of Stellar Evolution Pre-Main Sequence Main Sequence Post-Main Sequence The Primary definition is thus what is the Main Sequence Locus of core H burning Burning Process
More informationStar Formation and Protostars
Stellar Objects: Star Formation and Protostars 1 Star Formation and Protostars 1 Preliminaries Objects on the way to become stars, but extract energy primarily from gravitational contraction are called
More informationAST111 PROBLEM SET 6 SOLUTIONS
AST111 PROBLEM SET 6 SOLUTIONS Homework problems 1. Ideal gases in pressure balance Consider a parcel of molecular hydrogen at a temperature of 100 K in proximity to a parcel of ionized hydrogen at a temperature
More informationP M 2 R 4. (3) To determine the luminosity, we now turn to the radiative diffusion equation,
Astronomy 715 Final Exam Solutions Question 1 (i). The equation of hydrostatic equilibrium is dp dr GM r r 2 ρ. (1) This corresponds to the scaling P M R ρ, (2) R2 where P and rho represent the central
More informationSubstellar Interiors. PHY 688, Lecture 13
Substellar Interiors PHY 688, Lecture 13 Outline Review of previous lecture curve of growth: dependence of absorption line strength on abundance metallicity; subdwarfs Substellar interiors equation of
More informationNuclear Reactions and Solar Neutrinos ASTR 2110 Sarazin. Davis Solar Neutrino Experiment
Nuclear Reactions and Solar Neutrinos ASTR 2110 Sarazin Davis Solar Neutrino Experiment Hydrogen Burning Want 4H = 4p è 4 He = (2p,2n) p è n, only by weak interaction Much slower than pure fusion protons
More information1.1 An introduction to the cgs units
Chapter 1 Introduction 1.1 An introduction to the cgs units The cgs units, instead of the MKS units, are commonly used by astrophysicists. We will therefore follow the convention to adopt the former throughout
More informationOverview of Astronomical Concepts III. Stellar Atmospheres; Spectroscopy. PHY 688, Lecture 5 Stanimir Metchev
Overview of Astronomical Concepts III. Stellar Atmospheres; Spectroscopy PHY 688, Lecture 5 Stanimir Metchev Outline Review of previous lecture Stellar atmospheres spectral lines line profiles; broadening
More informationStars Lab 2: Polytropic Stellar Models
Monash University Stars Lab 2: Polytropic Stellar Models School of Mathematical Sciences ASP3012 2009 1 Introduction In the last lab you learned how to use numerical methods to solve various ordinary differential
More informationStellar Models ASTR 2110 Sarazin
Stellar Models ASTR 2110 Sarazin Jansky Lecture Tuesday, October 24 7 pm Room 101, Nau Hall Bernie Fanaroff Observing the Universe From Africa Trip to Conference Away on conference in the Netherlands
More informationStellar Interiors. Hydrostatic Equilibrium. PHY stellar-structures - J. Hedberg
Stellar Interiors. Hydrostatic Equilibrium 2. Mass continuity 3. Equation of State. The pressure integral 4. Stellar Energy Sources. Where does it come from? 5. Intro to Nuclear Reactions. Fission 2. Fusion
More informationStellar Magnetospheres part deux: Magnetic Hot Stars. Stan Owocki
Stellar Magnetospheres part deux: Magnetic Hot Stars Stan Owocki Key concepts from lec. 1 MagRe# --> inf => ideal => frozen flux breaks down at small scales: reconnection Lorentz force ~ mag. pressure
More informationUNIVERSITY OF SOUTHAMPTON
UNIVERSITY OF SOUTHAMPTON PHYS3010W1 SEMESTER 2 EXAMINATION 2014-2015 STELLAR EVOLUTION: MODEL ANSWERS Duration: 120 MINS (2 hours) This paper contains 8 questions. Answer all questions in Section A and
More informationIntroduction. Stellar Objects: Introduction 1. Why should we care about star astrophysics?
Stellar Objects: Introduction 1 Introduction Why should we care about star astrophysics? stars are a major constituent of the visible universe understanding how stars work is probably the earliest major
More information2. Basic Assumptions for Stellar Atmospheres
2. Basic Assumptions for Stellar Atmospheres 1. geometry, stationarity 2. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres!
More informationAstrophysics Assignment; Kramers Opacity Law
Astrophysics Assignment; Kramers Opacity Law Alenka Bajec August 26, 2005 CONTENTS Contents Transport of Energy 2. Radiative Transport of Energy................................. 2.. Basic Estimates......................................
More informationTHIRD-YEAR ASTROPHYSICS
THIRD-YEAR ASTROPHYSICS Problem Set: Stellar Structure and Evolution (Dr Ph Podsiadlowski, Michaelmas Term 2006) 1 Measuring Stellar Parameters Sirius is a visual binary with a period of 4994 yr Its measured
More informationStellar Interiors ASTR 2110 Sarazin. Interior of Sun
Stellar Interiors ASTR 2110 Sarazin Interior of Sun Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any other materials or any person Bring
More informationStellar Winds: Mechanisms and Dynamics
Astrofysikalisk dynamik, VT 010 Stellar Winds: Mechanisms and Dynamics Lecture Notes Susanne Höfner Department of Physics and Astronomy Uppsala University 1 Most stars have a stellar wind, i.e. and outflow
More informationPAPER 68 ACCRETION DISCS
MATHEMATICAL TRIPOS Part III Tuesday, 2 June, 2009 9:00 am to 11:00 am PAPER 68 ACCRETION DISCS There are THREE questions in total Full marks can be obtained by completing TWO questions The questions carry
More informationAtomic Physics 3 ASTR 2110 Sarazin
Atomic Physics 3 ASTR 2110 Sarazin Homework #5 Due Wednesday, October 4 due to fall break Test #1 Monday, October 9, 11-11:50 am Ruffner G006 (classroom) You may not consult the text, your notes, or any
More informationPHAS3135 The Physics of Stars
PHAS3135 The Physics of Stars Exam 2013 (Zane/Howarth) Answer ALL SIX questions from Section A, and ANY TWO questions from Section B The numbers in square brackets in the right-hand margin indicate the
More informationAy101 Set 1 solutions
Ay11 Set 1 solutions Ge Chen Jan. 1 19 1. The scale of the Sun a 3 points Venus has an orbital period of 5 days. Using Kepler s laws, what is its semi-major axis in units of AU Start with Kepler s third
More informationExercise: A Toy Model for Dust-driven Winds
Astrofysikalisk dynamik, VT 00 Exercise: A Toy Model for Dust-driven Winds Susanne Höfner Department of Physics and Astronomy, Uppsala University Cool luminous giants stars, in particular pulsating AGB
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationASTR-3760: Solar & Space Physics...Spring 2017
ASTR-3760: Solar & Space Physics...Spring 2017 Review material for midterm exam (March 22, 2017) Although I m not recommending full-on memorization of everything in this document, I do think it s important
More informationProblem set: solar irradiance and solar wind
Problem set: solar irradiance and solar wind Karel Schrijver July 3, 203 Stratification of a static atmosphere within a force-free magnetic field Problem: Write down the general MHD force-balance equation
More informationIntroduction to the Sun
Lecture 15 Introduction to the Sun Jiong Qiu, MSU Physics Department Open Q: what physics do we learn about the Sun? 1. Energy - nuclear energy - magnetic energy 2. Radiation - continuum and line emissions;
More informationStellar Spectra ASTR 2120 Sarazin. Solar Spectrum
Stellar Spectra ASTR 2120 Sarazin Solar Spectrum Solar Prominence Sep. 14, 1999 Solar Activity Due to rotation, convection, and magnetic field (Section 7.2 review) Charged Particles in Magnetic Fields
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Quiz 2
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.282 April 16, 2003 Quiz 2 Name SOLUTIONS (please print) Last First 1. Work any 7 of the 10 problems - indicate clearly which 7 you want
More informationwhere v is the velocity of the particles. In this case the density of photons is given by the radiative energy density : U = at 4 (5.
5.3. CONECIE RANSOR where v is the velocity of the particles. In this case the density of photons is given by the riative energy density : U = a 4 (5.11) where a is the riation density constant, and so
More informationMira variables They have very regular light curves, with amplitude V > 2.5 m, and periods Π > 100 days.
5 Pulsations Most AGB stars are variable stars, due to the fact that they pulsate. The periods are long, from 10 to 1000s of days, and they are therefore called Long Period Variables, or LPVs. These long
More informationLecture 2. Relativistic Stars. Jolien Creighton. University of Wisconsin Milwaukee. July 16, 2012
Lecture 2 Relativistic Stars Jolien Creighton University of Wisconsin Milwaukee July 16, 2012 Equation of state of cold degenerate matter Non-relativistic degeneracy Relativistic degeneracy Chandrasekhar
More informationPhysics 576 Stellar Astrophysics Prof. James Buckley. Lecture 13 Thermodynamics of QM particles
Physics 576 Stellar Astrophysics Prof. James Buckley Lecture 13 Thermodynamics of QM particles Reading/Homework Assignment Read chapter 3 in Rose. Midterm Exam, April 5 (take home) Final Project, May 4
More informationAstronomy, Astrophysics, and Cosmology
Astronomy, Astrophysics, and Cosmology Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson VI March 15, 2016 arxiv:0706.1988 L. A. Anchordoqui (CUNY)
More informationThe Virial Theorem for Stars
The Virial Theorem for Stars Stars are excellent examples of systems in virial equilibrium. To see this, let us make two assumptions: 1) Stars are in hydrostatic equilibrium 2) Stars are made up of ideal
More informationdf dz = dp dt Essentially, this is just a statement of the first law in one of the forms we derived earlier (expressed here in W m 3 ) dq p dt dp
A problem with using entropy as a variable is that it is not a particularly intuitive concept. The mechanics of using entropy for evaluating system evolution is well developed, but it sometimes feels a
More informationAy 1 Lecture 8. Stellar Structure and the Sun
Ay 1 Lecture 8 Stellar Structure and the Sun 8.1 Stellar Structure Basics How Stars Work Hydrostatic Equilibrium: gas and radiation pressure balance the gravity Thermal Equilibrium: Energy generated =
More informationThe Chaplygin dark star
The Chaplygin dark star O. Bertolami, J. Páramos Instituto Superior Técnico, Departamento de Física, Av. Rovisco Pais 1, 149-1 Lisboa, Portugal and E-mail addresses: orfeu@cosmos.ist.utl.pt; x jorge@fisica.ist.utl.pt
More informationNote: Each problem is worth 14 points except numbers 5 and 6 which are 15 points. = 3 2
Math Prelim II Solutions Spring Note: Each problem is worth points except numbers 5 and 6 which are 5 points. x. Compute x da where is the region in the second quadrant between the + y circles x + y and
More informationLecture 5-2: Polytropes. Literature: MWW chapter 19
Lecture 5-2: Polytropes Literature: MWW chapter 9!" Preamble The 4 equatios of stellar structure divide ito two groups: Mass ad mometum describig the mechaical structure ad thermal equilibrium ad eergy
More informationPHY-105: Nuclear Reactions in Stars (continued)
PHY-105: Nuclear Reactions in Stars (continued) Recall from last lecture that the nuclear energy generation rate for the PP reactions (that main reaction chains that convert hyogen to helium in stars similar
More information1. The vertical structure of the atmosphere. Temperature profile.
Lecture 4. The structure of the atmosphere. Air in motion. Objectives: 1. The vertical structure of the atmosphere. Temperature profile. 2. Temperature in the lower atmosphere: dry adiabatic lapse rate.
More informationStellar Structure. Observationally, we can determine: Can we explain all these observations?
Stellar Structure Observationally, we can determine: Flux Mass Distance Luminosity Temperature Radius Spectral Type Composition Can we explain all these observations? Stellar Structure Plan: Use our general
More informationPhysics 160: Stellar Astrophysics. Midterm Exam. 27 October 2011 INSTRUCTIONS READ ME!
Physics 160: Stellar Astrophysics 27 October 2011 Name: S O L U T I O N S Student ID #: INSTRUCTIONS READ ME! 1. There are 4 questions on the exam; complete at least 3 of them. 2. You have 80 minutes to
More informationLecture 7.1: Pulsating Stars
Lecture 7.1: Pulsating Stars Literature: KWW chapter 41!" 1 a) Classes of pulsating stars Many stars Intrinsically variable Subset: regular pulsation Type Period (d) RR Lyrae 0.3-0.9 Cepheids 1-50 W Virginis
More information2. Basic assumptions for stellar atmospheres
. Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 20: Stellar evolution: The giant stage 1 Energy transport in stars and the life time on the main sequence How long does the star remain on the main sequence? It will depend on the
More informationdp dr = GM c = κl 4πcr 2
RED GIANTS There is a large variety of stellar models which have a distinct core envelope structure. While any main sequence star, or any white dwarf, may be well approximated with a single polytropic
More informationHomework 1. Astronomy 202a. Fall 2009
Homework 1 Astronomy 0a Fall 009 Solutions Problems: 1. A galaxy has an integrated blue magnitude of 6.5, a rotation velocity (essentially flat) of 107 km s 1, and is ellptical in overall shape with major
More informationStellar Winds. Star. v w
Stellar Winds Star v w Stellar Winds Geoffrey V. Bicknell 1 Characteristics of stellar winds Solar wind Velocity at earth s orbit: Density: Temperature: Speed of sound: v 400 km/s n 10 7 m 3 c s T 10 5
More informationStellar Structure & Evolution
Stellar Structure & Evolution Sommersemester 2007 Ralf Klessen Zentrum für Astronomie der Universität Heidelberg Institut für Theoretische Astrophysik rklessen@ita.uni-heidelberg.de Albert-Überle-Straße
More informationPHYS-333: Problem set #1 Solutions
PHYS-333: Problem set #1 Solutions VERSION of February 10, 2018. 1. Energy flux and magnitude: a. Suppose two objects have energy fluxes of f and f + f, where f f. Derive an approximate expression for
More informationAST Homework V - Solutions
AST 341 - Homework V - Solutions TA: Marina von Steinkirch, steinkirch@gmail.com State University of New York at Stony Brook November, 010 1 (1 point) Derive the homologous form of the luminosity equation
More informationATMO 551a Fall The Carnot Cycle
What is a arnot ycle and Why do we care The arnot ycle arnot was a French engineer who was trying to understand how to extract usable mechanical work from a heat engine, that is an engine where a gas or
More informationAy123 Set 3 solutions
Ay123 Set 3 solutions Mia de los Reyes October 23 218 1. Equation of state and the Chandrasekhar mass (a) Using the Fermi-Dirac distribution for non-relativistic electrons, derive the relationship between
More informationScale analysis of the vertical equation of motion:
Scale analysis of the vertical equation of motion: As we did with the hz eqns, we do for the vertical to estimate the order of magnitude of Dw/ we take the largest of the terms, Dw -- W/T h UW/L W 2 /H
More informationPreliminary Examination: Astronomy
Preliminary Examination: Astronomy Department of Physics and Astronomy University of New Mexico Spring 2017 Instructions: Answer 8 of the 10 questions (10 points each) Total time for the test is three
More information