Math 2.0

Astronomy and astrophysics Stellar structure Proposition

Hydrostatic equilibrium

Builds on Euler equationsPoisson's equationNewton's law of universal gravitationLane–Emden equationVirial theorem

Statement

Let a fluid with density ρ\rho, pressure PP and velocity vv fill a region of ℝ3\R^3, in a gravitational potential Φ\Phi. Suppose that it obeys the Euler equations ρ(∂tv+(v⋅∇)v)=−∇P−ρ∇Φ\rho (\partial_t v + (v \cdot \nabla) v) = -\nabla P - \rho \nabla \Phi. Suppose that the fluid is at rest, that is, v=0v = 0 everywhere and at all times.

  1. Then ∇P=−ρ∇Φ\nabla P = -\rho \nabla \Phi.
  2. Suppose moreover that Φ\Phi solves Poisson's equation ΔΦ=4πGρ\Delta \Phi = 4 \pi G \rho, where GG is the gravitational constant. Suppose also that the region is a ball centred at 00, and that ρ\rho, PP and Φ\Phi are smooth and depend only on r=∣x∣r = \abs{x}. Let m(r)=∫0r4πs2ρ(s) dsm(r) = \int_0^r 4 \pi s^2 \rho(s) \, ds be the mass inside radius rr. Then for r>0r > 0, dPdr=−Gm(r)ρ(r)r2.\frac{dP}{dr} = -\frac{G m(r) \rho(r)}{r^2}.

Remarks

Part 1 is the Euler equation with v=0v = 0. For part 2, integrate Poisson's equation over the ball of radius rr and apply the divergence theorem. This gives 4πr2Φ′(r)=4πGm(r)4 \pi r^2 \Phi'(r) = 4 \pi G m(r), so Φ′(r)=Gm(r)/r2\Phi'(r) = G m(r)/r^2, which is Newton's law of gravitation for the mass inside rr. Example: in a ball of uniform density ρ0\rho_0 and radius RR with P(R)=0P(R) = 0, integrating gives the central pressure P(0)=23πGρ02R2P(0) = \frac{2}{3} \pi G \rho_0^2 R^2. Example: in a thin layer near the ground with constant g=∣∇Φ∣g = \abs{\nabla \Phi}, part 1 gives dP/dz=−ρgdP/dz = -\rho g, the barometric law of a column of air or water. The hypothesis v=0v = 0 is needed: in a rotating star the centrifugal term adds to ∇Φ\nabla \Phi and the surfaces of equal pressure are no longer spheres.

Going deeper

For a student

In the equations of stellar structure, hydrostatic equilibrium is the equation for dP/drdP/dr. Combined with a relation between pressure and density, such as a polytropic law P=Kρ1+1/nP = K \rho^{1 + 1/n}, it reduces to a single second-order equation, the Lane–Emden equation. Multiplying by 4πr34 \pi r^3 and integrating over the star gives the virial theorem for a self-gravitating gas: 3∫P dV3 \int P \, dV equals minus the gravitational energy. In general relativity the equation is replaced by the Tolman–Oppenheimer–Volkoff equation, which leads to an upper bound on the mass of neutron stars.

Reading path

To understand this page from the ground up, read these 85 pages in order. The foundations come first.

  1. Derivative concept
  2. Differential equation concept
  3. Euclidean space object
  4. Euclidean vector concept
  5. Euler equations concept
  6. Function concept
  7. Interval concept
  8. Lane–Emden equation concept
  9. Limit point of a set definition
  10. Local maximum definition
  11. Ordinary differential equation concept
  12. Partition of an interval definition
  13. Poisson's equation concept
  14. Sequence definition
  15. Velocity concept
  16. Fréchet derivative definition
  17. Cross product definition
  18. Ordinary differential equation definition
  19. Function of several variables concept
  20. Vector-valued function concept
  21. Lipschitz condition definition
  22. Classical mechanics concept
  23. Limit of a function definition
  24. Local minimum definition
  25. Mesh of a partition definition
  26. Riemann sum definition
  27. Real sequence definition
  28. Subsequence definition
  29. Limit of a function of several variables concept
  30. Momentum concept
  31. Partial derivative concept
  32. Solution of an ordinary differential equation definition
  33. Inertial frame of reference concept
  34. Continuity at a point definition
  35. Derivative definition
  36. Riemann integrable function definition
  37. Bounded sequence definition
  38. Limit of a sequence definition
  39. Monotone sequence definition
  40. Continuity in R^n definition
  41. Differentiable function of several variables definition
  42. Partial derivative definition
  43. Initial value problem definition
  44. Newton's laws of motion statement
  45. Continuous function definition
  46. Antiderivative definition
  47. Riemann integral definition
  48. Fermat's theorem on stationary points statement
  49. Convergent sequence definition
  50. Monotone convergence theorem for sequences statement
  51. Every sequence has a monotone subsequence statement
  52. Differentiable on an open set definition
  53. Gradient definition
  54. Higher partial derivatives definition
  55. Picard–Lindelöf hypotheses definition
  56. Intermediate value theorem statement
  57. Definite integral definition
  58. Bolzano–Weierstrass theorem statement
  59. Iterated derivative definition
  60. Picard operator definition
  61. Additivity of the integral over intervals statement
  62. Boundedness theorem statement
  63. Monotonicity of the integral statement
  64. C^k function characterised by continuous derivatives statement
  65. Integral equation form of an initial value problem statement
  66. An iterate of the Picard operator is a contraction statement
  67. Extreme value theorem statement
  68. Vector field definition
  69. Picard–Lindelöf theorem statement
  70. Mean value theorem for integrals statement
  71. Rolle's theorem statement
  72. Conservative vector field definition
  73. Mean value theorem statement
  74. First fundamental theorem of calculus statement
  75. Classical mechanics definition
  76. Potential energy definition
  77. Momentum definition
  78. Second fundamental theorem of calculus statement
  79. Kinetic energy definition
  80. Angular momentum definition
  81. Fundamental theorem of calculus statement
  82. Central force definition
  83. Virial theorem statement
  84. Newton's law of universal gravitation statement
  85. Hydrostatic equilibrium statement

Where this leads

22 pages build on this one.

Further reading

  1. Web https://en.wikipedia.org/wiki/Hydrostatic_equilibrium
About this page