P--kpT +- at4 (3) ,AH 3. 10l'-7(r) =- IK?dP, 2. M(r)/M' P ~O' L(r) the amount of radiant energy crossing a spherical surface of radius r,

Similar documents
Convection When the radial flux of energy is carried by radiation, we derived an expression for the temperature gradient: dt dr = - 3

CHAPTER 8 ENTROPY GENERATION AND TRANSPORT

Fundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres

Entropy generation and transport

Equations of Stellar Structure

Stellar structure Conservation of mass Conservation of energy Equation of hydrostatic equilibrium Equation of energy transport Equation of state

Handout 18: Convection

9.1 Introduction. 9.2 Static Models STELLAR MODELS

GS388 Handout: Radial density distribution via the Adams-Williamson equation 1

EVOLUTION OF STARS: A DETAILED PICTURE

Convection. If luminosity is transported by radiation, then it must obey

Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow

Free-Fall Timescale of Sun

3 Hydrostatic Equilibrium

Astronomy 404 October 9, 2013

Physics 556 Stellar Astrophysics Prof. James Buckley. Lecture 9 Energy Production and Scaling Laws

SPA7023P/SPA7023U/ASTM109 Stellar Structure and Evolution Duration: 2.5 hours

ASTM109 Stellar Structure and Evolution Duration: 2.5 hours

Stellar Pulsations and Variability

partium superarum plerumque compressis, intergeniculis ramorum primariorum,

PHYS 633: Introduction to Stellar Astrophysics

Lecture 7: Stellar evolution I: Low-mass stars

Model Atmospheres. Model Atmosphere Assumptions

Ay 1 Lecture 8. Stellar Structure and the Sun

The Euler Equation of Gas-Dynamics

UNIVERSITY OF SOUTHAMPTON

2. Equations of Stellar Structure

Exam II: Solutions. UNIVERSITY OF ALABAMA Department of Physics and Astronomy. PH 125 / LeClair Spring 2009

PHYS 432 Physics of Fluids: Instabilities

dp dr = GM c = κl 4πcr 2

The Interiors of the Stars

Astro Instructors: Jim Cordes & Shami Chatterjee.

Introduction. Stellar Objects: Introduction 1. Why should we care about star astrophysics?

High Energy Astrophysics

BERNOULLI EQUATION. The motion of a fluid is usually extremely complex.

Physics 576 Stellar Astrophysics Prof. James Buckley. Lecture 10

P M 2 R 4. (3) To determine the luminosity, we now turn to the radiative diffusion equation,

Star Formation and Protostars

[2 marks] Show that derivative of the angular velocity. What is the specific angular momentum j as a function of M and R in this Keplerian case?

Lecture Notes for PHY 405 Classical Mechanics

NOTES ON BINARY STARS. H1. ON THE LIGHT OFECLIPSING BINARIES OUT OF ECLIPSE

while the Planck mean opacity is defined by

Evolution Beyond the Red Giants

FLOWS IN GAS PIPELINES

Where does Bernoulli's Equation come from?

PHYSICS 218 EXAM 2 Tuesday, October 26, 2010

read 9.4-end 9.8(HW#6), 9.9(HW#7), 9.11(HW#8) We are proceding to Chap 10 stellar old age

Stellar structure and evolution. Pierre Hily-Blant April 25, IPAG

Topics ASTR 3730: Fall 2003

First Law of Thermodynamics

Linear Theory of Stellar Pulsation

The structure and evolution of stars. Introduction and recap

Stellar Interiors. Hydrostatic Equilibrium. PHY stellar-structures - J. Hedberg

2. Stellar atmospheres: Structure

Adiabatic expansion Isothermal compression Adiabatic compression

PROCEEDINGS. The present investigation will show that these two line elements, together

6.2 Governing Equations for Natural Convection

Outline. Today we will learn what is thermal radiation

The atmosphere in motion: forces and wind. AT350 Ahrens Chapter 9

Overview spherical accretion

THIRD-YEAR ASTROPHYSICS

Energy transport: convection

Introduction to the Sun

(2) secondary radiations induced by them in the earth's atmosphere or

Classical Mechanics. FIG. 1. Figure for (a), (b) and (c). FIG. 2. Figure for (d) and (e).

Stars and HR Diagrams

F q. Gas at radius R (cylindrical) and height z above the disk midplane. F z. central mass M

White Dwarfs - Degenerate Stellar Configurations

1 Introduction to Governing Equations 2 1a Methodology... 2

The Tropical Atmosphere: Hurricane Incubator

Mechanical Engineering. Postal Correspondence Course HEAT TRANSFER. GATE, IES & PSUs

Astrophysical Fluid Dynamics (2)

Chapter 2. Energy and the First Law of Thermodynamics

Please choose the letter corresponding to the best answer to each problem (5 pts each).

DIVIDED SYLLABUS ( ) - CLASS XI PHYSICS (CODE 042) COURSE STRUCTURE APRIL

ON ISING'S MODEL OF FERROMAGNETISM

Conductors and Dielectrics

Exercise 1: Vertical structure of the lower troposphere

Physics 556 Stellar Astrophysics Prof. James Buckley

3. Stellar radial pulsation and stability

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Chapter 7 Neutron Stars

High-Energy Astrophysics Lecture 6: Black holes in galaxies and the fundamentals of accretion. Overview

COX & GIULI'S PRINCIPLES OF STELLAR STRUCTURE

11/19/08. Gravitational equilibrium: The outward push of pressure balances the inward pull of gravity. Weight of upper layers compresses lower layers

1.1 Motivation. 1.2 The H-R diagram

Nuclear Reactions and Solar Neutrinos ASTR 2110 Sarazin. Davis Solar Neutrino Experiment

n v molecules will pass per unit time through the area from left to

Introduction to the School

Lecture 14: The Sun and energy transport in stars. Astronomy 111

10/17/2012. Stellar Evolution. Lecture 14. NGC 7635: The Bubble Nebula (APOD) Prelim Results. Mean = 75.7 Stdev = 14.7

AP Physics C Mechanics Objectives

Phases of Stellar Evolution

Summary of stellar equations

Physics 576 Stellar Astrophysics Prof. James Buckley. Lecture 11 Polytropes and the Lane Emden Equation

ATMO 551a Fall The Carnot Cycle

EART164: PLANETARY ATMOSPHERES

t KH = GM2 RL Pressure Supported Core for a Massive Star Consider a dense core supported by pressure. This core must satisfy the equation:

The Origin of Gravity and the Impetus for Celestial Motions

AST1100 Lecture Notes

Transcription:

572 ASTRONOMY: S. CHANDRASEKHAR PROC. N. A. S. Species type: tetrasporic, Herb. Univ. Calif. No. 564574; St. James, Cape Peninsula, South Africa; Coll. Ecol. Surv. F.-123. Jania digitata is the largest species of the genus yet known. The general branching is dichotomous-cymoid, but oftentimes three or more branches arise digitately, in one plane, from the apices of the intergenicula. THE STABILITY OF THE RADIATIVE GRADIENTS IN THE INTERIOR OF A STAR By S. CHANDRASE}HAR YERKES OBSERVATORY OF THIE UNIVERSITY OF CHICAGO Communicated October 11, 1937 If a spherically symmetrical distribution of matter is in hydrostatic equilibrium and if, further, radiative equilibrium obtains, then the radiative temperature gradient is determined by [see Note 1 at the end of the paper] dt 1 K-c dp (1) T 4 K(r) P where K is the opacity coefficient, 17 L (r) I * ( ) = 1 XP (2) 10l'-7(r) =- IK?dP, 2 M(r)/M' P ~O' L(r) the amount of radiant energy crossing a spherical surface of radius r, and M(r) is the mass interior to r. Finally P--kpT +- at4 (3),AH 3 where k, a, 1u and H are respectively the Boltzmann constant, the Stefan- Boltzmann (radiation) constant, the molecular weight and the mass of the proton. To examine the stability of the radiative gradient we suppose that an element of mass bm originally at temperature T, density p and pressure P suffers a sudden increase of temperature of amount ST > 0. Then this element would exert a pressure of a definite amount SP > 0 on its surroundings thus expanding and becoming less dense than its immediate neighborhood. The element Sm would then experience a force tending to displace it into regions of lower density. During such a movement the element would continue to expand, the temperature altering meantime.

VOL. 23, 1937 ASTRONOMY: S. CHANDRASSEKHAR We shall now make the following assumptions: (1) at each instant of time the element Am expands (or contracts) to such an extent that the pressure exerted on the element by the surrounding material is the same as the element exerts on the surrounding material; (2) that the process of expansion (or contraction) takes place adiabatically; (3) that the viscous forces restricting the movement of the element Am can be neglected. We shall first examine the consequences of these assumptions and later subject them to a closer scrutiny. By our second assumption, since the expansion of Am takes place adiabatically, we should have [see Note 2 at the end of the paper]: dt r_ -1 dp (4) T r2 P where r72 r2 =1 + (32 ~ + 3(y (4-3#))Qy-l) - 1)(1 - )(4 + ) 5 k 11-4 3P=-^ pt; (l-3)p=- at4. (6) 1.H 3 (4) governs the rate of change of the temperature of the element 3m as it moves outward into regions of lower density, expanding meantime. Comparing (1) and (4) we see that the temperature of 3m as it moves outward alters at a rate different from that of its immediate surroundings, because according to our first assumption P alters in the same way for both 3m and the surroundings. Let us now suppose that r2-1 1 K17 r2 4 "(r) Then it follows that the element 3m after moving outward for a certain distance will find itself at the same temperature, pressure and density as its surroundings at that point and consequently the original disturbance eventually dies out. In the same way, if the element bm originally suffers a decrease of temperature of amount 3T then it would become denser than its immediate surroundings and would consequently sink to regions of higher density. If we make our three assumptions as before, then the adiabatic compression it experiences as it sinks to regions of higher density increases the temperature of bm at a rate greater than does its surroundings. Again it will soon find itself at a point where 3m and its neighborhood at that point have the same pressure, density and temperature. Thus in either case-i.e., either a positive or a negative increment OT of the element 3m-the dis- 573

574 ASTRONOMY: S. CHANDRASEKHAR PROC. N. A. S. turbance dies out if (7) is true. In this sense, the radiative equilibrium is stable if the radiative gradient is less than the corresponding adiabatic gradient-a result first proved by Schwarzschild, but the assumptions on which it can be proved are generally not made explicit. On the other hand if r2-1 1 K? (8) P2 4 K-(r) then if the element bm suffers an increase (decrease) of temperature, then as before it would move outward (downward) to regions of lower (higher) density but now the temperature of bm decreases (increases) less than does that of its surroundings, and hence is always at a temperature higher (lower) than its surroundings. Hence we have proved on the basis of our assumptions that the radiative equilibrium is stable or unstable according as 4 r2 1 > or < _ 07 (9) r2 "q~~k(r) The following table gives the values of 4(r2 - i)/r2 for different values of (1-1 -, 4(r2-1)/r2 -, 4 (r2 -i)/r2 I - d3 4 (r2- )/r? 0 1.6 0.4 1.065 0.8 1.004 0.1 1.304 0.5 1.039 0.9 1.000 0.2 1.177 0.6 1.022 1.0 1.0000 0.3 1.108 0.7 1.010...... Suppose now that we have initially a situation in which the radiative gradient obtains but one which exceeds the adiabatic gradient. By our above discussion the radiative gradient is unstable in the sense that a slight alteration of the local temperature will give rise to a system of ascending and descending currents which would have the effect of reducing the existing radiative temperature gradient. It is clear that eventually a steady state would be set up though it is not a priori clear what the nature of that steady state would be. We shall continue our discussion on the basis of the three fundamental assumptions we have already made. It is generally assumed that eventually the temperature gradient that would be set up would very approximately be the adiabatic one. Actually, of course, the adiabatic gradient would not be rigorously set up. For there will always be a certain amount of transport of heat energy by radiation and the presumption is that eventually there would be a "slight" excess of the currents in the direction of increasing temperature such that the net transport of heat energy due to the currents is positive in the direction of increasing temperature and is further of an amount exactly equal to the transport of heat energy by the radiation in the opposite direction. Now the

VOL. 23, 1937 ASTRONOMY: S. CHANDRASEKHAR 575 transport of heat energy by the radiation will be governed by the equation der = cdpr = - 4rr() pdr. (10) If this were so, it follows that the law of the conservation of momentunm would require the excess currents in the direction of increasing temperature to transfer momentum of an amount exactly that transferred by radiation in the opposite direction. Hence on these premises the equations of equilibrium are clearlv k d GM(r) IAH d p)= - r2 (11) dp r2 dt P r2-it ~~~~~~~(12) It should be remarked that in the literature3 equation (12) for the temperature gradient has been mistakenly combined with d (H PT + 1 at4) = GM(r) P. (13) Equations (12) and (13) together form an inconsistent set of equations, for the use of (13) as the equation of hydrostatic equilibrium already implies the existence of the radiative temperature gradient. The simplest way of seeing this is perhaps the following: Equation (13) assumes that dpr(v)dv represents the rate of transfer of momentum across an elementary infinitesimal cylindrical slab at r due to radiation in the frequency interval (v, v + dv). On the other hand it is a general theorem that the mechanical force exerted by radiation is given by - K,pFdr/c. These two statements together imply precisely the existence of the radiative temperature gradient which contradicts the use of (12) in combination with (13). We have now to examine the fundamental assumptions on which our discussion has so far been carried. The first important point to realize is that if the radiation pressure is important-in other words if f8 and (1-I ) are of comparable orders of magnitude then the system of "convection" currents set up must be quite "violent." For, as we have seen, the excess currents in the direction of increasing temperature transfer momentum of amount exactly equal to that transferred by radiation in the opposite direction, and hence the resulting pressure due to mass motion must be of the same order as the gas kinetic pressure, and this implies a state of quite active stirring. It is then a priori probable that the inertia of motions is still quite large, in which case the elements of ascending and descending masses would experience considerable viscous friction (contrary to one of our hypotheses that the viscous friction could be neglected) which in turn

576 ASTRONOMY: S. CHANDRASEKHAR PROC. N. A. S. would communicate to the elements appreciable amounts of heat contradicting our assumption (2) that the expansion (or contraction) of the elements take place adiabatically. Finally our hypothesis (1) that at each instant of time the element bm expands (or contracts) to such an extent that the pressure exerted on the element Am by the surrounding material is the same as the element exerts on the surrounding material, assumes that during the movements of the element of mass the equalization of density in the element due to the diffusion of the surrounding material is negligible, which need not in general be the case. On the whole, it must be concluded that we cannot at present make any definite predictions about the nature of the "convection currents" that would be set up if the radiation pressure is important and if further the radiative gradient exceeds the adiabatic. It, however, seems probable that in the limit of vanishing radiation pressure the equations (11) and (12) may approximate very closely to the actual equations of equilibrium. In conclusion, I wish here to record my appreciation of the valuable discussions I have had on the subject of this note with Professors B. Stromgren and G. P. Kuiper. NOTE 1. The equations of equilibrium for material in radiative equilibrium are dp GM(r) dpr KL(r) Po' (ii). dr 47rcr2 k 1 P=Pr + P = pt+ at4. (iii) Dividing (ii) by (i) and remembering our definition of q we have or integrating which can also be written as By (iv) and (vi) we have or L dpr = 4cGM P, (iv) L P Pr = 47cGM JoXdP, (v) L Pr= K(r). P. (vi) dp Kq7 dp Pr K?7(r) P' (vii)

VOL. 23, 1937 PHYSIOLOG Y: H. B. TUKE Y 577 which is the equation quoted. NOTE 2. For an adiabatic change we should have dt 1 K7 dp T 4 K?(r) P' (Vii) I-uV dv + (- dt + PdV =0, (i) where U is the intemal energy of the system. For our case we have U = at4v + cvt; P = (cp - cv) T + -at4, (ii) V 3 where cp and cv are the two specific heats of the gas. From (i) and (ii) we derive that (2.+ (12p7 1 )dt + -1 Pe ) T + (4pr d + P) - = 0. (iii) On the other hand r2 is defined by the relation or according to (ii) dp + r2 dt P 1-i2 0 (iv) [4Pr + P +1r2 (Pr +P) ]T -P dv =0. (v) Comparing (iii) and (v) we derive the expression for r2 quoted. A condition equivalent to (9) has been given by H. Siedentopf (A. N., 244, 278 (1932)). NOTE 3. Bierman, L., Zeits. f. Ap., 5, 128 (1932); A. N., 257, 270 (1935). PLANT PARTS OF DECIDUOUS FRUITS WHICH GIVE EVIDENCE OF BEING FAVORABLE FOR TISSUE CULTURE* By H. B. TUKEY NEW YORK STATE AGRICULTURAL EXPERIMENT STATION, GENEvA, NEW YORK Communicated October 11, 1937 White,6'7 in an early review of plant tissue cultures, called attention to the fact that, as in animal tissue cultures, certain plant tissues might be found to respond to treatment more readily than others. Among these he suggested young embryos, root tips and seed primordia. Root tips have since proved well suited to culture and manipulation in vitro, as called to attention by Robbins,2 and as developed.by White,8 Robbins and Bartley' and others. More recently Bonner' has shown that the