Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mso.anu.edu.au/~geoff/AGD/Sound_Waves.pdf
Дата изменения: Mon Feb 22 14:03:18 2016
Дата индексирования: Sun Apr 10 03:16:11 2016
Кодировка:

Поисковые слова: molecular cloud
Sound
References: L.D. Landau & E.M. Lifshitz: Fluid Mechanics, Chapter VIII F. Shu: The Physics of Astrophysics, Vol. 2, Gas Dynamics, Chapter 8

1 Speed of sound The phenomenon of sound waves is one that can be understood using the fundamental equations of gas dynamics that we have developed so far.
Astrophysical Gas Dynamics: Sound waves 1/52


1.1 Perturbation equations Simplest case: fluid equations without magnetic field and gravity (Euler equations): + vi = 0 t xi v i v i p + vj =­ xj xi t

(1)

Astrophysical Gas Dynamics: Sound waves

2/52


Consider small perturbations to state = 0 , and v i = v 0 i . Without loss of generality we may take v 0 i = 0 . = 0 + v i = v 0 i + v i Expand to first order: v i = 0 + v i = 0 v i + O 2 v i vj = O2 xj
(3) (2)

Astrophysical Gas Dynamics: Sound waves

3/52


The symbol O 2 means second order in the quantities , v i , i.e. terms such as 2 , v i v j , v i etc. The perturbation to the pressure is determined by the perturbations to the density and entropy: p = p s p p p = + s s
(4)

If we take the flow to be adiabatic, then s = 0 . Define the parameter c s by:
2 = p cs 2 p = c s p 0 0 (5)
4/52

s

Astrophysical Gas Dynamics: Sound waves


The perturbation equations become: + 0 v i = 0 t xi 2 = 0 0 vi + cs t xi

(6)

Time derivative of first equation and divergence of second equation: 2 + 0 ------------ v i = 0 t xi t2 2 2 2 0 ------------ v i + c s -------------- = 0 xi t xi xi
Astrophysical Gas Dynamics: Sound waves

2

(7)

5/52


Subtract => 2 2 2 ------------ ­ c s -------------- = 0 xi xi t2
(8)

This the wave equation for a disturbance moving at velocity c s . Speed of sound = c s
(9)

Astrophysical Gas Dynamics: Sound waves

6/52


1.2 Plane-wave solutions Take = A exp i k x ­ t 2 ------------ = ­ 2 t2 2 -------------- = ­ k 2 xi xi 1 2 ­ --- 2 cs 2 2 2 ------------ = ­ k + ------ = 0 2 t2 cs = cs k
Astrophysical Gas Dynamics: Sound waves 7/52

(10)


What does this plane wave represent? z y k n = -k Consider a surface of constant phase = k x ­ t = 0
(11)

x

From the expression for the density, this represents a surface of constant density whose location evolves with time. This surface is planar and for our purposes is best represented in the form: k 0 -- x = --- t + ----k k k
(12)

Astrophysical Gas Dynamics: Sound waves

8/52


From the above equation, we read off the normal to the plane: k n = -k 0 D = --- t + ----k k
(13)

and the perpendicular distance of the plane from the origin is
(14)

Therefore, the speed at which the plane is moving away from the origin is: --- = c s k
(15)

Astrophysical Gas Dynamics: Sound waves

9/52


Thus the solution represents a wave moving at a speed c s in the direction of the wave vector k Solution for velocity The velocity perturbations are derived by considering perturbations of the form: v i = A i exp i k x ­ t
(16)

Note the appearance of a vector A i for the amplitude since v i is a vector.

Astrophysical Gas Dynamics: Sound waves

10/52


The time and spatial derivatives of v i are: v i = ­ i A i exp i k x ­ t t = ik i A exp i k x ­ t xi Substitute into perturbation equation for velocity: 2 = 0 0 vi + cs t xi
(18)

(17)

Astrophysical Gas Dynamics: Sound waves

11/52


2 ­ i 0 A i + c s ik i A = 0

ki 2 A i = c s --

ki k A A 2-- ----- = c s --- --- ----0 k 0 A ----0

(19)

ki = c s --k A ----- Ai = cs 0

The amplitude of the velocity perturbation
(20)

Astrophysical Gas Dynamics: Sound waves

12/52


i.e. Amplitude of Relative amplitude = cs velocity perturbation density perturbation What does small mean? ----- « 1 v i « c s Nature of sound wave Since ki A i = c s --k A ----0
(23)
13/52

(21)

(22)

i.e perturbation velocities much less than the speed of sound.

Astrophysical Gas Dynamics: Sound waves


then the velocity of the oscillating elements of gas are in the direction of the wave-vector, i.e. the wave is longitudinal. 1.3 Numerical value of the speed of sound Take the equation of state p = K s
2 = p cs

s

­ 1 = p = Ks -----

(24)



kT = -------m
Astrophysical Gas Dynamics: Sound waves 14/52


Therefore: cs = kT -------m
(25)

N.B. The speed of sound in an ideal gas only depends upon temperature. e.g. air at sea level on a warm day = 1.4 T = 300 K = 28.8 c s = 340 m/s T = 10 7 K atmosphere in an elliptical galaxy
7 K , = 0.62 , = 5 T = 10 -(26)

3

c s = 480 km/s
Astrophysical Gas Dynamics: Sound waves

(27)
15/52


2 Subsonic and supersonic flow Disturbances in a gas travel at the speed of sound relative to the fluid. i.e. information in the gas propagates at the sound speed

Sound wave

v = cs

V cs

Astrophysical Gas Dynamics: Sound waves

16/52


A wave travelling at a velocity of v = c s wrt to the fluid cannot send a signal backwards to warn of an impending obstacle. This leads to the formation of shocks (to be considered later). Therefore the speed V = c s is a critical one in gas dynamics. Supersonic Flow Subsonic Flow Mach Number v cs v cs v M = ---cs
(28)

(29)

Nobody ever heard the bullet that killed him -Th. von Karman
Astrophysical Gas Dynamics: Sound waves 17/52


3 Energy & momentum in sound waves 3.1 Expressions for energy density and energy flux Energy density 12 tot = -- v + 2 Expand out the quantities in this equation to second order s = 0 + + s s s 1 2) + s + 1 s 2 + -( -22 s 2 s2
Astrophysical Gas Dynamics: Sound waves

(30)

2

2

2

(31)

18/52


Since we are only considering adiabatic fluctuations, then s = 0 and + 1 2 = s = 0 + -2 2 s
2 (32)

We now evaluate these terms using thermodynamic identities. Since, 1 h kTds = -- d ­ -- d then d = hd + kTds = h s
Astrophysical Gas Dynamics: Sound waves

(33)

(34)

19/52


and 2 = h = ------- 2 s s + p ----------
(35)

1 1 p = ­ ----- + p + -+ ---- 2
2 2 cs cs hh -= ­ -- + -- + ----- = ----

Astrophysical Gas Dynamics: Sound waves

20/52


Hence the total (kinetic + internal) energy is to second order: 1 1 tot = 0 + -- 0 v 2 + h 0 + -- 2 2
2 cs ----- 2 0 (36)

The term h 0 is associated with a change in energy in a given volume due to a change in mass and eventually disappears. Energy flux 1 v 2 + h F E i = v i -2 = 0 h 0 v i + 0 h v i + h 0 v i + O 3
Astrophysical Gas Dynamics: Sound waves 21/52

(37)


Note that the kinetic energy term is of third order. We need to relate h to the variations in other thermodynamic variables. Since 1 kTds = dh ­ -- dp then 1 1 h = -dh = -- dp + kTds p s and
2 cs 1 h = -- p = -----
Astrophysical Gas Dynamics: Sound waves

(38)

(39)

(40)
22/52


Hence, using the previous expression for F E i : F Ei = 0 h 0 v i + h 0 v i = 0 h 0 v i + h 0 v i + c + 0 h vi 2 v i
2 cs ----- 2 0 (41)

s

We combine this with the expression for the energy density: 1 2 + 1 tot = 0 + h 0 + -- 0 v -- 2 2 tot F Ei + =0 t xi
Astrophysical Gas Dynamics: Sound waves

(42)

Because of the energy equation, we have the conservation law:
(43)
23/52


Some of the terms in these expressions can be simplified. First, the leading red term in , 0 is constant and is just the background energy density. Now examine the continuity equation written out to second order: 0 v i + v i = 0 + t xi Multiplying this by h 0 gives: h 0 + h 0 0 v i + h 0 v i = 0 t xi
(45) (44)

Astrophysical Gas Dynamics: Sound waves

24/52


so that when we insert the red and blue terms in the energy conservation equation, they drop out and we are left with t 1 2 + 1 -- -- 0 v 2 2
2 cs 2 + c 2 v = 0 ----- i 0 xi s (46)

Hence, we take for the energy density and its associated energy flux
2 cs sw = 1 v 2 + 1 ----- 2 -- 0 -- 2 2 0 sw 2 F Ei = c s v i = p v i

(47)

Astrophysical Gas Dynamics: Sound waves

25/52


3.2 Energy density and energy flux in a plane wave Express and v i in terms of travelling waves and take the real part: = A exp i k x ­ t A cos k x ­ t ki A - v i = c s --- ----- exp i k x ­ t k 0 ki c s --k A ----- cos k x ­ t 0

(48)

Astrophysical Gas Dynamics: Sound waves

26/52


Hence,
2 2 cs cs sw = 1 v 2 + 1 ----- 2 = ----- A 2 cos 2 k x ­ t -- 0 -- 2 2 0 0 3 cs ki ki sw 2 F Ei = c s v i = ----- A 2 --- os 2 k x ­ t = c s sw ---c o k k

(49)

Average the energy density and flux over a period T = 2 of the wave using: 1T 2 k x ­ t dt = 1 -- cos -T0 2



(50)

Astrophysical Gas Dynamics: Sound waves

27/52


so that the average energy density and flux are:
2 cs sw = 1 ----- A 2 -- 2 0 3 cs ki ki 1 F Ei = -- ----- A 2 --- = c s sw ---k 2 o k

(51)

Notes · The energy flux is in the direction of the wave. · and is equal to the sound speed times the energy density.

Astrophysical Gas Dynamics: Sound waves

28/52


4 Jeans mass - sound waves with self gravity 4.1 Physical motivation There are numerous sites in the interstellar medium of our own galaxy that are the birthplaces of stars. These are cold dense clouds in which stars can form as a result of the process of gravitational collapse. Why is ``cold'' and ``dense'' important? Surprisingly, we can get some idea of this by looking at the propagation of sound waves in a molecular cloud.

Astrophysical Gas Dynamics: Sound waves

29/52


The Horsehead nebula is a dark molecular cloud that is the site of ongoing star formation. The image at the left was obtained at the Kitt Peak National Observatory in Arizona

The image at the right was obtained using the Hubble Space Telescope.
Astrophysical Gas Dynamics: Sound waves 30/52


4.2 Mathematical treatment of gravitational instability Continuity, momentum and Poisson equations: + vi = 0 t xi v i v i 1 p + vj + -+ =0 t xj xi xi 2 = 4 G In order to model self-gravity of the gas we have included Poisson's equation for the gravitational potential.
(52)

Astrophysical Gas Dynamics: Sound waves

31/52


Perturbations: = 0 + vi = vi = 0 + 0 constant
(53)

Astrophysical Gas Dynamics: Sound waves

32/52


Perturbation equations: 0 0 + + 0 v i = + 0 ( v i) + vi t xi t xi xi v i 1 p 0 1 (54) p 0 0 + ----+ ----- p ­ ----+ + = 0 2 t 0 xi 0 xi 0 xi xi xi 2 0 + 2 = 4 G 0 + Zeroth order terms are shown in blue. The term in red is associated with what is sometimes unkindly known as the Jeans Swindle. Jeans neglected this term in his derivation of the Jeans Mass.
Astrophysical Gas Dynamics: Sound waves 33/52


When can we neglect the pressure gradient term? 1 p 0 Compare ----- p with ----. Let us suppose that the back2 0 xi 0 xi ground pressure gradient has a length scale L , and that the perturbations have a length scale l i.e. 1 p ----- p ------0 xi 0 l and p 0 ---- ----2 2 0 xi 0 p0 ----L
(56) (55)

Astrophysical Gas Dynamics: Sound waves

34/52


The pressure gradient term associated with the perturbations is larger than the background pressure gradient term if p p 0 l p 0 2 ---- = ------- » ----- ----- -- « ----- ----- = c 0 -2 0 l 2 L L p 0 c0 0
(57)

So, for perturbation length scales much less than the background lenght scale in the cloud, we can neglect the background pressure gradient.

Astrophysical Gas Dynamics: Sound waves

35/52


Now use the zeroth order equations: 1 p 0 0 ----+ =0 0 xi xi 2 0 = 4 G 0

(58)

Astrophysical Gas Dynamics: Sound waves

36/52


This gives: 0 vi = 0 + 0 ( v i) + t xi xi V i 1 p 0 + ----- p ­ ----+ = 0 2 t 0 xi 0 xi xi 2 = 4 G Simplest case: initial density constant and neglect variation of potential. 0 = constant 0 = constant
Astrophysical Gas Dynamics: Sound waves

(59)

(60) (61)
37/52


The assumption means that the red term above is neglected. The perturbation equations simplify to: + 0 ( v i) = 0 t xi v i 1 + ----- p + = 0 0 xi xi t 2 = 4 G Plane wave solution: = A exp i k x ­ t v i = A i exp i k x ­ t = B exp i k x ­ t
Astrophysical Gas Dynamics: Sound waves

(62)

(63)

38/52


Substitute into perturbation equations: ­ i A + 0 ik i A i = 0
2 cs ­ i A i + ----- ik i A + i k i B = 0 0 (64)

­ k 2 B = 4 GA

Astrophysical Gas Dynamics: Sound waves

39/52


We can arrange these equations into the following set of linear equations: A ­ 0 ki Ai = 0
2 cs ----- k i A ­ A i + k i B = 0 0 (65)

4 GA + k 2 B = 0

Astrophysical Gas Dynamics: Sound waves

40/52


There are a number of ways of dealing with this set of linear algebraic equations for the quantities, A A i and B . One of the best is simply to write them as a set of 5 equations as follows:

Astrophysical Gas Dynamics: Sound waves

41/52



2 cs ----- k 1 0 2 cs ----- k 2 0 2 cs ----- k 3 0

­0 k1 ­0 k2 ­0 k3 0 ­ 0 ­ 0 k1 0 0 k2 A2 = 0 0 A3 0 k3 B k2 A A1

0

0 ­ 0

(66)

0 0

0 0

4G

Astrophysical Gas Dynamics: Sound waves

42/52


This is a homogeneous set of equations that only has a non-trivial solution if the determinant is zero. It is possible to work the determinant out by hand, since it contains a number of zeroes. However, Maple/Mathematica also readily provides the following expression for the determinant:
2 = 2 k 2 ­ 2 + cs k 2 ­ 4 G 0 (67)

Apart from the trivial solutions ( or k = 0 ) to = 0 , the only non-trivial solutions are described by:
2 2 = cs k 2 ­ 4 G 0 (68)

giving
2 = cs k 2 ­ 4 G 0
Astrophysical Gas Dynamics: Sound waves

(69)
43/52


This constitutes an interesting and fundamental difference from normal sound waves because of the additional term relating to self-gravity. Define the Jeans wave number k J by
22 cs kJ = 4 G 0 2 2 2 = cs k 2 ­ kJ (70)

Then if k k J 2 0 i = g say
(71)

Astrophysical Gas Dynamics: Sound waves

44/52


For these imaginary solutions the perturbation in each variable is proportional to exp i t = exp g t .
(72)

One of the imaginary roots corresponds to decaying modes, the other to growing modes, i.e. instability.
2 cs 2 Jeans Length = J = ----- = 2 ----------------kJ 4 G 0 (73)

Astrophysical Gas Dynamics: Sound waves

45/52


Jeans length in a molecular cloud n 10 9 m ­ 3 T 10 K 1
(74)

kT 1 / 2 c s = ------- 0.29 km/s m kT 1 ------------------ 1 / 2 = ------J = mG m 0 kT 1 / 2 2pc ----------- Gn

Astrophysical Gas Dynamics: Sound waves

46/52


Jeans mass The Jeans mass, M J is defined to be the mass in the region defined by the reciprocal of the Jeans wave number:
­ MJ = 0 kJ 3

= 0

2 cs 3 / 2 ----------------- 4 G 0

(75)

3/2 kT = ------------------------------- 1 4 G m 0 / 3 0.4 solar masses
Astrophysical Gas Dynamics: Sound waves 47/52


for the above parameters. Star formation - the modern approach The above discussion is the standard treatment for gravitational instability without the influence of a magnetic field and gives interesting sizes for the initial collapsing region and the mass of the collapsing object. The current attitude to the Jeans mass is that the physics of star formation is more complicated and that magnetic fields and turbulence are involved. However, it is thought that the Jeans Mass is relevant to the masses of molecular cloud cores.

Astrophysical Gas Dynamics: Sound waves

48/52


Physics of the Jean mass

M

R

Consider a gas cloud of mass M and radius R and neglect pressure forces. Then the equation of motion of the cloud is: dv i GM x i ------ = ­ ­ -------- ---dt R2 R
Astrophysical Gas Dynamics: Sound waves

(76)
49/52


Dimensionally, this equation is: R GM R3 1 / 2 --- = -------- Free-fall time = t ff ------- G ­ 1 / 2 (77) GM t2 R2 Now consider the sound crossing time: R t s ---cs
(78)

Astrophysical Gas Dynamics: Sound waves

50/52


If the sound-crossing time is less then the free-fall time, then the collapsing region is able to produce enough pressure to halt the collapse since signals have the time to go from one side of the region to the other. Thus the condition for collapse is R t s t ff ---- G ­ 1 / 2 R ------- cs G Another criterion: Euler's equations with pressure and gravity: dv i 1 p ------ = ­ -­ dt xi xi
(80) 2 1/2 cs

J

(79)

Astrophysical Gas Dynamics: Sound waves

51/52


Collapse will occur if gravitational forces overwhelm the pressure forces. Since p p 2 = = cs xi xi xi then Euler's equations are:
2 2 2 dv i c s c s GM c s GM ------ = ­ ----­ ----- -- ­ -------- ----- ­ -------dt xi xi R R2 R R2 (82) (81)

For collapse the gravitational forces have to win, so that,
2 2 1/2 cs GM c s -------- ----- R ------- J R3 R2 G
Astrophysical Gas Dynamics: Sound waves

(83)
52/52