User:Tohline/SSC/Stability Eulerian Perspective
From VisTrailsWiki
Contents 
Stability of Spherically Symmetric Configurations (Eulerian Perspective)
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
A standard technique that is used throughout astrophysics to test the stability of selfgravitating fluids involves perturbing physical variables away from their initial (usually equilibrium) values then linearizing each of the principal governing equations before seeking solutions describing the timedependent behavior of the variables that simultaneously satisfy all of the equations. When the effects of the fluid's self gravity are ignored and this analysis technique is applied to an initially homogeneous medium, the combined set of linearized governing equations generates a wave equation that governs the propagation of sound waves. Here we build on our separate, introductory discussion of sound waves and apply standard perturbation & linearization techniques to spherically symmetric, inhomogeneous and selfgravitating fluids. We will assume that the reader has read this separate introductory discussion and, in particular, understands how the linear wave equation that governs the propagation of sound waves is derived from the set of nonlinear, principal governing equations.
Assembling the Key Relations
Governing Equations and Supplemental Relations
We begin with the set of principal governing equations that provides the foundation for all of our discussions in this H_Book, namely, the
Eulerian Representation
of the Continuity Equation,
Eulerian Representation
of the Euler Equation,
Adiabatic Form of the
First Law of Thermodynamics
.
Poisson Equation

As was done in our separate, introductory discussion of sound waves, we will assume that we are dealing with an ideal gas and supplement this set of equations with a barotropic (polytropic) equation of state,
… with …
which will ensure that the adiabatic form of the first law of thermodynamics is satisfied. When we develop the linearized Euler equation, below, it will be useful to recognize that, assuming is uniform throughout the fluid, we can rewrite this last expression as,






Perturbation then Linearization of Equations
In this Eulerian analysis, we are investigating how conditions vary with time at a fixed point in space, . By analogy with our separate introductory analysis of sound waves, we will write the four primary variables in the form,












where quantities with subscript "0" are initial values — independent of time, but not necessarily spatially uniform, and usually specified via some choice of an initial equilibrium configuration — and quantities with subscript "1" denote variations away from the initial state, which are assumed to be small in amplitude — for example, and . As indicated, we will assume that the fluid configuration is initially stationary and, for simplicity, will not append the subscript "1" to the velocity perturbation. It is to be understood, however, that the velocity, , is small also where, ultimately, this will mean . In what follows, by definition, , , , and are considered to be of first order in smallness, while products of these quantities are of second (or even higher) order in smallness.
Continuity Equation
Substituting the expression for into the lefthand side of the continuity equation and neglecting small quantities of the second order, we have,






where, in the first line, the first term on the righthand side has been set to zero because is independent of time. Hence, we have the desired,
Linearized Continuity Equation



Euler Equation
Next, we turn to the Euler equation and note that the term,

may be altogether neglected because it is of second order in smallness. Substituting the expressions for , , and into the righthand side of the Euler equation and neglecting small quantities of the second order, we have,









where, the binomial theorem has been used to obtain the expression on the righthand side of the second line and, in the last line, the sum of the first pair of terms has been set to zero because the initial configuration is assumed to be in equilibrium. Combining these simplification steps, we have the,
Linearized Euler Equation



First Law of Thermodynamics
In a similar fashion, perturbing the variables in the barotropic equation of state leads to,















Hence, as in our separate introductory discussion of sound waves, we have the,
Linearized Equation of State



Poisson Equation
Finally, plugging the "perturbed" expressions for and into the Poisson equation — which, by its very nature, is a linear equation — we have,






where, in the second line, the sum of the pair of terms on the lefthand side has been set to zero because it is a selfcontained representation of the Poisson equation for the initial unperturbed medium. This gives us the desired,
Linearized Poisson Equation



Summary and Combinations
In summary, the following four linearized equations govern the timedependent physical relationship between the four perturbation amplitudes , , and in selfgravitating fluids:
Linearized Linearized Linearized
Linearized

This set of linearized governing equations is more general than the set of equations that has traditionally been used to describe the propagation of sound waves. The equations governing soundwave propagation can be retrieved by choosing an initially homogeneous medium — in which case, and — and by ignoring the fluid's self gravity, that is, by setting .
As is explicitly demonstrated in our discussion of Bonnor's (1957) research, below, the linearized Euler equation can be combined with the linearized adiabatic form of the first law of thermodynamics to give,



Taking the divergence of this relation, we have,






where, in order to obtain this last expression, we have used the linearized Poisson equation to replace in favor of . Alternatively, taking the time derivative of the linearized continuity equation gives,






Combining these last two expressions, then, gives us a

that describes the timevariation at any point in space of the fractional density fluctuation,
in a selfgravitating, barotropic fluid. For purposes of comparison with the wave equation that governs the behavior of sound waves, it is worth rewriting this expression in the form,



where, as before,
Comparison with Classic Research Publications
James Jeans (1902 & 1928)
James H. Jeans (1902, Philosophical Transactions of the royal Society of London. Series A, 199, 1) used precisely this type of perturbation and linearization analysis when he first derived what is now commonly referred to as the Jeans Instability. For example, if our discussion is restricted only to fluctuations in the radial coordinate direction of a spherically symmetric configuration — in which case and — our expression for the linearized Euler equation exactly matches equation (12) from Jeans (1902), which, for purposes of illustration, is displayed in the following framed image.
Linearized Euler Equation as Derived and Presented by Jeans (1902) 


The correspondence between the righthandsides of equation (12) from Jeans (1902) and our derived expression for the linearized Euler equation is clear after accepting the following variable mappings: The lefthand side of equation (12) from Jeans (1902) also matches the lefthand side of our linearized Euler equation, although this may not be immediately apparent. In the paper by Jeans, is not a component of the velocity vector but is, rather, the radial displacement of a fluid element. Hence, 
A broader discussion of gravitational instability in the context of the formation of "great nebulae" (i.e., galaxies) and stars appears in Chapter XIII — specifically, pp. 337342 — of the book by Jeans (1928) titled, "Astronomy and Cosmogony." The governing wave equation for selfgravitating fluids that we have derived, above, appears as equation (314.6) in this published discussion by Jeans (1928), although the term on the lefthand side involving does not appear, presumably because Jeans assumed that the initial, unperturbed medium was homogeneous. In an effort to facilitate comparison with our derived expression, equation (314.6) from Jeans (1928) has been reprinted here as a framed image.
Wave Equation for SelfGravitating Fluids as Derived and Presented by Jeans (1928) 

Note that Jeans (1928) uses , without the subscript "0", to represent the initial density, and , rather than "G", for the Newtonian gravitational constant. 
W. B. Bonnor (1957)
Above, in our opening layout of the governing equations and supplemental relations, we pointed out that,



If we make this substitution in our linearized Euler equation, and also use the linearized first law of thermodynamics to replace in favor of , we obtain,









This is the version of the linearized Euler equation that was derived by Bonnor (1957, MNRAS, 117, 104) in the section of his paper that addresses the growth of Newtonian, selfgravitating fluctuations on a static (cosmological) background. The two equation images reproduced in the following outlined box document Bonnor's (1957) initial expression (his equation 2.1) for the nonlinear Euler equation and his derived expression (equation 2.7) for the linearized Euler equation. After allowing for the identified variable mappings, Bonnor's two expressions precisely match, respectively, the form of the nonlinear Euler equation that is included among our set of principal governing equations, and this last derived form of our linearized Euler equation.
Bonnor's (1957, MNRAS, 117, 104) Derivation 


Original nonlinear Euler Equation 

Linearized Euler Equation 

The correspondence between these two equations from Bonnor (1957) and our derived expressions is clear after accepting the following variable mappings: 
Bonnor then proceeded to combine the full set of linearized governing equations, in the manner we have detailed above, into a wave equation that is appropriately modified to handle selfgravitating fluids. In an effort to facilitate comparison with our derived expression, Bonnor's (1957) modified wave equation (2.10) has been reprinted here as a framed image.
Wave Equation for SelfGravitating Fluids as Derived and Presented by Bonnor (1957) 

See Also
 Part II of Spherically Symmetric Configurations: Stability
 Wave Equation
 Sound Waves and Gravitational Instability — class notes provided online by David H. Weinberg (The Ohio State University)
Scratch Work
Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
 Go Home 
Part 1















Linearized

Replace in favor of in Euler equation.
Take the divergence of this entire expression, then use linearized Poisson equation to replace in favor of .
Rearrange terms in the linearized equation of continuity, then take the partial timederivative of the entire expression.
Subtract the step #2 expression from the step #3 expression.
If homentropic as well, then,
where:
Wave Equation for SelfGravitating Fluids as Derived and Presented by Jeans (1928) 

Part 2
① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨ ⑩
Linearized

 Term involving :





















Part 3



Wave Equation for SelfGravitating Fluids as Derived and Presented by P. Ledoux & Th. Walraven (1958) 

Wave Equation for SelfGravitating Fluids as Derived and Presented by S. Rosseland (1969) 

Part 4
↵
⏎
For any scalar variable, , the relationship between a Lagrangian (total) and Eulerian (partial) timederivative is,









© 2014  2020 by Joel E. Tohline 