Hydrostatic equilibrium explained

In fluid mechanics, hydrostatic equilibrium (hydrostatic balance, hydrostasy) is the condition of a fluid or plastic solid at rest, which occurs when external forces, such as gravity, are balanced by a pressure-gradient force.^[1] In the planetary physics of Earth, the pressure-gradient force prevents gravity from collapsing the planetary atmosphere into a thin, dense shell, whereas gravity prevents the pressure-gradient force from diffusing the atmosphere into outer space.^{[2] [3]} In general, it is what causes objects in space to be spherical.

Hydrostatic equilibrium is the distinguishing criterion between dwarf planets and small solar system bodies, and features in astrophysics and planetary geology. Said qualification of equilibrium indicates that the shape of the object is symmetrically rounded, mostly due to rotation, into an ellipsoid, where any irregular surface features are consequent to a relatively thin solid crust. In addition to the Sun, there are a dozen or so equilibrium objects confirmed to exist in the Solar System.

Mathematical consideration

For a hydrostatic fluid on Earth:$$dP\; =\; -\; \backslash rho(P)\; \backslash ,\; g(h)\; \backslash ,\; dh$$

Derivation from force summation

Newton's laws of motion state that a volume of a fluid that is not in motion or that is in a state of constant velocity must have zero net force on it. This means the sum of the forces in a given direction must be opposed by an equal sum of forces in the opposite direction. This force balance is called a hydrostatic equilibrium.

The fluid can be split into a large number of cuboid volume elements; by considering a single element, the action of the fluid can be derived.

There are three forces: the force downwards onto the top of the cuboid from the pressure, P, of the fluid above it is, from the definition of pressure,$$F\_\backslash text\; =\; -\; P\_\backslash text\; A$$Similarly, the force on the volume element from the pressure of the fluid below pushing upwards is$$F\_\backslash text\; =\; P\_\backslash text\; A$$

Finally, the weight of the volume element causes a force downwards. If the density is ρ, the volume is V and g the standard gravity, then:$$F\_\backslash text\; =\; -\backslash rho\; g\; V$$The volume of this cuboid is equal to the area of the top or bottom, times the height – the formula for finding the volume of a cube.$$F\_\backslash text\; =\; -\backslash rho\; g\; A\; h$$

By balancing these forces, the total force on the fluid is$$\backslash sum\; F\; =\; F\_\backslash text\; +\; F\_\backslash text\; +\; F\_\backslash text\; =\; P\_\backslash text\; A\; -\; P\_\backslash text\; A\; -\; \backslash rho\; g\; A\; h$$This sum equals zero if the fluid's velocity is constant. Dividing by A,$$0\; =\; P\_\backslash text\; -\; P\_\backslash text\; -\; \backslash rho\; g\; h$$Or,$$P\_\backslash text\; -\; P\_\backslash text\; =\; -\; \backslash rho\; g\; h$$P_top − P_bottom is a change in pressure, and h is the height of the volume element—a change in the distance above the ground. By saying these changes are infinitesimally small, the equation can be written in differential form.$$dP\; =\; -\; \backslash rho\; g\; \backslash ,\; dh$$Density changes with pressure, and gravity changes with height, so the equation would be:$$dP\; =\; -\; \backslash rho(P)\; \backslash ,\; g(h)\; \backslash ,\; dh$$

Derivation from Navier–Stokes equations

Note finally that this last equation can be derived by solving the three-dimensional Navier–Stokes equations for the equilibrium situation where$$u\; =\; v\; =\; \backslash frac\; =\; \backslash frac\; =\; 0$$Then the only non-trivial equation is the

z -equation, which now reads$$\backslash frac\; +\; \backslash rho\; g\; =\; 0$$Thus, hydrostatic balance can be regarded as a particularly simple equilibrium solution of the Navier–Stokes equations.

Derivation from general relativity

By plugging the energy–momentum tensor for a perfect fluid$$T^\; =\; \backslash left(\backslash rho\; c^\; +\; P\backslash right)\; u^\backslash mu\; u^\backslash nu\; +\; P\; g^$$into the Einstein field equations$$R\_\; =\; \backslash frac\; \backslash left(T\_\; -\; \backslash frac\; g\_\; T\backslash right)$$and using the conservation condition$$\backslash nabla\_\backslash mu\; T^\; =\; 0$$one can derive the Tolman–Oppenheimer–Volkoff equation for the structure of a static, spherically symmetric relativistic star in isotropic coordinates:$$\backslash frac\; =\; -\backslash frac\; \backslash left(1+\backslash frac\backslash right)\; \backslash left(1+\backslash frac\backslash right)\; \backslash left(1\; -\; \backslash frac\backslash right)^$$In practice, Ρ and ρ are related by an equation of state of the form f(Ρ,ρ) = 0, with f specific to makeup of the star. M(r) is a foliation of spheres weighted by the mass density ρ(r), with the largest sphere having radius r:$$M(r)\; =\; 4\backslash pi\; \backslash int\_0^r\; dr\text{'}\; \backslash ,\; r\text{'}^2\; \backslash rho(r\text{'}).$$Per standard procedure in taking the nonrelativistic limit, we let, so that the factor$$\backslash left(1+\backslash frac\backslash right)\; \backslash left(1+\backslash frac\backslash right)\; \backslash left(1-\backslash frac\; \backslash right)^\; \backslash rightarrow\; 1$$Therefore, in the nonrelativistic limit the Tolman–Oppenheimer–Volkoff equation reduces to Newton's hydrostatic equilibrium:$$\backslash frac\; =\; -\backslash frac\; =\; -g(r)\backslash ,\backslash rho(r)\backslash longrightarrow\; dP\; =\; -\; \backslash rho(h)\backslash ,g(h)\backslash ,\; dh$$(we have made the trivial notation change h = r and have used f(Ρ,ρ) = 0 to express ρ in terms of P).^[4] A similar equation can be computed for rotating, axially symmetric stars, which in its gauge independent form reads:$$\backslash frac\; -\; \backslash partial\_i\; \backslash ln\; u^t\; +\; u\_t\; u^\backslash varphi\backslash partial\_i\backslash frac=0$$ Unlike the TOV equilibrium equation, these are two equations (for instance, if as usual when treating stars, one chooses spherical coordinates as basis coordinates

(t,r,\theta,\varphi) , the index i runs for the coordinates r and \theta ).

Applications

Fluids

The hydrostatic equilibrium pertains to hydrostatics and the principles of equilibrium of fluids. A hydrostatic balance is a particular balance for weighing substances in water. Hydrostatic balance allows the discovery of their specific gravities. This equilibrium is strictly applicable when an ideal fluid is in steady horizontal laminar flow, and when any fluid is at rest or in vertical motion at constant speed. It can also be a satisfactory approximation when flow speeds are low enough that acceleration is negligible.

Astrophysics and planetary science

From the time of Isaac Newton much work has been done on the subject of the equilibrium attained when a fluid rotates in space. This has application to both stars and objects like planets, which may have been fluid in the past or in which the solid material deforms like a fluid when subjected to very high stresses.In any given layer of a star there is a hydrostatic equilibrium between the outward-pushing pressure gradient and the weight of the material above pressing inward. One can also study planets under the assumption of hydrostatic equilibrium. A rotating star or planet in hydrostatic equilibrium is usually an oblate spheroid, that is, an ellipsoid in which two of the principal axes are equal and longer than the third.An example of this phenomenon is the star Vega, which has a rotation period of 12.5 hours. Consequently, Vega is about 20% larger at the equator than from pole to pole.

In his 1687 Philosophiæ Naturalis Principia Mathematica Newton correctly stated that a rotating fluid of uniform density under the influence of gravity would take the form of a spheroid and that the gravity (including the effect of centrifugal force) would be weaker at the equator than at the poles by an amount equal (at least asymptotically) to five fourths the centrifugal force at the equator.^[5] In 1742, Colin Maclaurin published his treatise on fluxions, in which he showed that the spheroid was an exact solution. If we designate the equatorial radius by

r _{e,} the polar radius by r _{p,} and the eccentricity by \epsilon, with\epsilon = \sqrt{ 1- r ^{2 },} _{e} he found that the gravity at the poles is^[6] \begin{align} g _{p} & =4 \pi (r _{p})/(r _{e}) (\epsilon r _{e-r} _{p\arctan(\epsilon} r _{e/r} _{p)})/(\epsilon ^{3}) G\rho \\ & =3 (\epsilon r _{e-r} _{p\arctan(\epsilon} r _{e/r} _{p)})/(\epsilon ^{3} r ^{3} _{e}) GM \\ \end{align} where G is the gravitational constant, \rho is the (uniform) density, and M is the total mass. The ratio of this to g _{0,} the gravity if the fluid is not rotating, is asymptotic tog _{p/g} _{0\sim} 1+ (1 { 15 }\epsilon)/(1+ ^{2\sim}) (2)/(15) f as \epsilon goes to zero, where f is the flattening:f = (r _{e-r} _{p})/(r _{e}) .

The gravitational attraction on the equator (not including centrifugal force) is

\begin{align} g _{e} & = (32 \left( (1 {r)/(_{e}))/(r) ^{p} - (\epsilon r _{e-r} _{p\arctan(\epsilon} r _{e/r} _{p)})/(\epsilon ^{3} r ^{2} _{e} r _{p})} \right)GM\\ & = (32 (r _{e\arctan(\epsilon} r _{e/r} _{p)-\epsilon} r _{p})/(\epsilon ^{3} r ^{3} _{e}))/(GM) \\ \end{align}

Asymptotically we have:

g _{e/g} _{0\sim} 1- (1 { 30 }\epsilon)/(1- ^{2\sim}) (1)/(15) f

Maclaurin showed (still in the case of uniform density) that the component of gravity toward the axis of rotation depended only on the distance from the axis and was proportional to that distance, and the component in the direction toward the plane of the equator depended only on the distance from that plane and was proportional to that distance. Newton had already pointed out that the gravity felt on the equator (including the lightening due to centrifugal force) has to be

(r _{p})/(r _{e}) g _{p} in order to have the same pressure at the bottom of channels from the pole or from the equator to the centre, so the centrifugal force at the equator must beg _{e - (r _{p})/(r _{e})} g _{p\sim} (25 \epsilon)/(^{2g}) ^{e\sim (45 fg)/(_{e.})}

Defining the latitude to be the angle between a tangent to the meridian and axis of rotation, the total gravity felt at latitude

\phi (including the effect of centrifugal force) isg(\phi) = (g _{p(1-f)})/(\sqrt{ 1- ( 2 f - f ^{2)\sin} ^{2\phi}) }.

This spheroid solution is stable up to a certain (critical) angular momentum (normalized by

M\sqrt{G\rho r _{e}} ), but in 1834 Carl Jacobi showed that it becomes unstable once the eccentricity reaches 0.81267 (or f reaches 0.3302).Above the critical value the solution becomes a Jacobi, or scalene, ellipsoid (one with all three axes different). Henri Poincaré in 1885 found that at still higher angular momentum it will no longer be ellipsoidal but piriform or oviform. The symmetry drops from the 8-fold D point group to the 4-fold C, with its axis perpendicular to the axis of rotation.^[7] Other shapes satisfy the equations beyond that, but are not stable, at least not near the point of bifurcation.^{[7] [8]} Poincaré was unsure what would happen at higher angular momentum, but concluded that eventually the blob would split in two.

The assumption of uniform density may apply more or less to a molten planet or a rocky planet, but does not apply to a star or to a planet like the earth which has a dense metallic core. In 1737 Alexis Clairaut studied the case of density varying with depth.^[9] Clairaut's theorem states that the variation of the gravity (including centrifugal force) is proportional to the square of the sine of the latitude, with the proportionality depending linearly on the flattening (

f ) and the ratio at the equator of centrifugal force to gravitational attraction. (Compare with the exact relation above for the case of uniform density.) Clairaut's theorem is a special case, for an oblate spheroid, of a connexion found later by Pierre-Simon Laplace between the shape and the variation of gravity.^[10]

If the star has a massive nearby companion object then tidal forces come into play as well, distorting the star into a scalene shape when rotation alone would make it a spheroid. An example of this is Beta Lyrae.

Hydrostatic equilibrium is also important for the intracluster medium, where it restricts the amount of fluid that can be present in the core of a cluster of galaxies.

We can also use the principle of hydrostatic equilibrium to estimate the velocity dispersion of dark matter in clusters of galaxies. Only baryonic matter (or, rather, the collisions thereof) emits X-ray radiation. The absolute X-ray luminosity per unit volume takes the form

l{L} _{X=Λ(T} _{B)\rho} ^{2} ^{B} where T _{B} and \rho _{B} are the temperature and density of the baryonic matter, and Λ (T) is some function of temperature and fundamental constants. The baryonic density satisfies the above equation $$p\_B(r+dr)-p\_B(r)=-dr\backslash frac\backslash int\_0^r\; 4\backslash pi\; r^2\backslash ,\backslash rho\_M(r)\backslash ,\; dr.$$The integral is a measure of the total mass of the cluster, with r being the proper distance to the center of the cluster. Using the ideal gas law p _{B=kT} _{B\rho} _{B/m} _{B} (k is the Boltzmann constant and m _{B} is a characteristic mass of the baryonic gas particles) and rearranging, we arrive at$$\backslash frac\backslash left(\backslash frac\backslash right)=-\backslash frac\backslash int\_0^r\; 4\backslash pi\; r^2\backslash ,\backslash rho\_M(r)\backslash ,\; dr.$$Multiplying by r ^{2 /\rho} _{B(r)} and differentiating with respect to r yields$$\backslash frac\backslash left[\backslash frac\{r^2\}\{\backslash rho\_B(r)\}\backslash frac\{d\}\{dr\}\backslash left(\backslash frac\{kT\_B(r)\backslash rho\_B(r)\}\{m\_B\}\backslash right)\backslash right]=-4\backslash pi\; Gr^2\backslash rho\_M(r).$$If we make the assumption that cold dark matter particles have an isotropic velocity distribution, then the same derivation applies to these particles, and their density \rho _{D=\rho} _{M-\rho} _{B} satisfies the non-linear differential equation$$\backslash frac\backslash left[\backslash frac\{r^2\}\{\backslash rho\_D(r)\}\backslash frac\{d\}\{dr\}\backslash left(\backslash frac\{kT\_D(r)\backslash rho\_D(r)\}\{m\_D\}\backslash right)\backslash right]=-4\backslash pi\; Gr^2\backslash rho\_M(r).$$With perfect X-ray and distance data, we could calculate the baryon density at each point in the cluster and thus the dark matter density. We could then calculate the velocity dispersion \sigma ^{2} _{D} of the dark matter, which is given by $$\backslash sigma^2\_D=\backslash frac.$$ The central density ratio \rho _{B(0)/\rho} _{M(0)} is dependent on the redshift z of the cluster and is given by$$\backslash rho\_B(0)/\backslash rho\_M(0)\backslash propto\; (1+z)^2\backslash left(\backslash frac\backslash right)^$$ where \theta is the angular width of the cluster and s the proper distance to the cluster. Values for the ratio range from 0.11 to 0.14 for various surveys.^[11]

Planetary geology

See also: List of gravitationally rounded objects of the Solar System.

The concept of hydrostatic equilibrium has also become important in determining whether an astronomical object is a planet, dwarf planet, or small Solar System body. According to the definition of planet adopted by the International Astronomical Union in 2006, one defining characteristic of planets and dwarf planets is that they are objects that have sufficient gravity to overcome their own rigidity and assume hydrostatic equilibrium. Such a body will often have the differentiated interior and geology of a world (a planemo), though near-hydrostatic or formerly hydrostatic bodies such as the proto-planet 4 Vesta may also be differentiated and some hydrostatic bodies (notably Callisto) have not thoroughly differentiated since their formation. Often the equilibrium shape is an oblate spheroid, as is the case with Earth. However, in the cases of moons in synchronous orbit, nearly unidirectional tidal forces create a scalene ellipsoid. Also, the purported dwarf planet is scalene due to its rapid rotation, though it may not currently be in equilibrium.

Icy objects were previously believed to need less mass to attain hydrostatic equilibrium than rocky objects. The smallest object that appears to have an equilibrium shape is the icy moon Mimas at 396 km, whereas the largest icy object known to have an obviously non-equilibrium shape is the icy moon Proteus at 420 km, and the largest rocky bodies in an obviously non-equilibrium shape are the asteroids Pallas and Vesta at about 520 km. However, Mimas is not actually in hydrostatic equilibrium for its current rotation. The smallest body confirmed to be in hydrostatic equilibrium is the dwarf planet Ceres, which is icy, at 945 km, whereas the largest known body to have a noticeable deviation from hydrostatic equilibrium is Iapetus being made of mostly permeable ice and almost no rock.^[12] At 1,469 km Iapetus is neither spherical nor ellipsoid. Instead, it is rather in a strange walnut-like shape due to its unique equatorial ridge.^[13] Some icy bodies may be in equilibrium at least partly due to a subsurface ocean, which is not the definition of equilibrium used by the IAU (gravity overcoming internal rigid-body forces). Even larger bodies deviate from hydrostatic equilibrium, although they are ellipsoidal: examples are Earth's Moon at 3,474 km (mostly rock),^[14] and the planet Mercury at 4,880 km (mostly metal).^[15]

In 2024, Kiss et al. found that has an ellipsoidal shape incompatible with hydrostatic equilibrium for its current spin. They hypothesised that Quaoar originally had a rapid rotation and was in hydrostatic equilibrium, but that its shape became "frozen in" and did not change as it spun down due to tidal forces from its moon Weywot.^[16] If so, this would resemble the situation of Iapetus, which is too oblate for its current spin.^{[17] [18]} Iapetus is generally still considered a planetary-mass moon nonetheless,^[19] though not always.^[20]

Solid bodies have irregular surfaces, but local irregularities may be consistent with global equilibrium. For example, the massive base of the tallest mountain on Earth, Mauna Kea, has deformed and depressed the level of the surrounding crust, so that the overall distribution of mass approaches equilibrium.

Atmospheric modeling

In the atmosphere, the pressure of the air decreases with increasing altitude. This pressure difference causes an upward force called the pressure-gradient force. The force of gravity balances this out, keeping the atmosphere bound to Earth and maintaining pressure differences with altitude.

See also

• List of gravitationally rounded objects of the Solar System; a list of objects that have a rounded, ellipsoidal shape due to their own gravity (but are not necessarily in hydrostatic equilibrium)
• Statics
• Two-balloon experiment

External links

• Strobel, Nick. (May, 2001). Nick Strobel's Astronomy Notes.
• by Richard Pogge, Ohio State University, Department of Astronomy

Notes and References

1. Book: White, Frank M. . Fluid Mechanics . 2008 . McGraw-Hill . New York . 978-0-07-128645-9 . Pressure Distribution in a Fluid . 63, 66.
2. Book: Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-scale Circulation. 9781139459969. Vallis. Geoffrey K.. 6 November 2006. Cambridge University Press .
3. Book: Ocean Circulation in Three Dimensions. 9780521768436. Klinger. Barry A.. Haine. Thomas W. N.. 14 March 2019. Cambridge University Press.
4. Book: Zee. A.. Einstein gravity in a nutshell. 2013. Princeton University Press . Princeton . 9780691145587. 451–454.
5. https://en.wikisource.org/wiki/The_Mathematical_Principles_of_Natural_Philosophy_(1846)/BookIII-Prop2 Propositions X-XXIV (Motions of celestial bodies and the sea)
6. Book: Colin Maclaurin . A Treatise on Fluxions . 1742 . 125 . Maclaurin does not use modern notation but rather gives his results in geometric terms. The gravity results are in article 646. At one point he makes an erroneous statement equivalent to d(\tan\theta - \theta)/d\tan\theta = \tan ^{2\theta} but his subsequent statements are correct.
7. Henri Poincaré . Les formes d'équilibre d'une masse fluide en rotation . Revue Général des Sciences Pures et Appliquées . 1892 .
8. Web site: Gallery : The shape of Planet Earth . Josleys.com . 2014-06-15.
9. Clairaut. Alexis. Colson. John. 1737. An Inquiry concerning the Figure of Such Planets as Revolve about an Axis, Supposing the Density Continually to Vary, from the Centre towards the Surface. 103921. Philosophical Transactions.
10. See Sir George Stokes . On Attractions, and on Clairaut's Theorem . The Cambridge and Dublin Mathematical Journal . 1849 . 194–219 . Sir George Stokes .
11. Book: Weinberg. Steven. Cosmology. 2008. Oxford University Press. New York. 978-0-19-852682-7. 70–71.
12. Thomas . P.C. . July 2010 . Sizes, shapes, and derived properties of the saturnian satellites after the Cassini nominal mission . Icarus . 208 . 1 . 395–401 . 10.1016/j.icarus.2010.01.025 . 2010Icar..208..395T . 23 December 2018 . https://web.archive.org/web/20181223003125/http://www.ciclops.org/media/sp/2011/6794_16344_0.pdf . dead .
13. Castillo-Rogez. J. C.. Matson, D. L. . Sotin, C. . Johnson, T. V. . Lunine, Jonathan I. . Thomas, P. C. . Iapetus' geophysics: Rotation rate, shape, and equatorial ridge. Icarus. 2007. 190. 1. 179–202. 10.1016/j.icarus.2007.02.018. 2007Icar..190..179C.
14. Garrick-Bethell . I. . Wisdom . J . Zuber . MT . Evidence for a Past High-Eccentricity Lunar Orbit . Science . 4 August 2006 . 313 . 5787 . 652–655 . 10.1126/science.1128237 . 16888135 . 2006Sci...313..652G . 317360 .
15. Sean Solomon, Larry Nittler & Brian Anderson, eds. (2018) Mercury: The View after MESSENGER. Cambridge Planetary Science series no. 21, Cambridge University Press, pp. 72–73.
16. etal . C. . Kiss . T. G. . Müller . G. . Marton . R. . Szakáts . A. . Pál . L. . Molnár . The visible and thermal light curve of the large Kuiper belt object (50000) Quaoar . Astronomy & Astrophysics . March 2024 . 684 . A50 . 10.1051/0004-6361/202348054 . 2401.12679 . 2024A&A...684A..50K.
17. Cowen, R. (2007). Idiosyncratic Iapetus, Science News vol. 172, pp. 104–106. references
18. 10.1016/j.icarus.2010.01.025. Thomas. P. C.. July 2010. Sizes, shapes, and derived properties of the saturnian satellites after the Cassini nominal mission. Icarus. 208. 1. 395–401. 2010Icar..208..395T. 2015-09-25. 2018-12-23. https://web.archive.org/web/20181223003125/http://www.ciclops.org/media/sp/2011/6794_16344_0.pdf. dead.
19. Emily Lakdawalla et al., What Is A Planet? The Planetary Society, 21 April 2020
20. Chen . Jingjing . Kipping . David . 2016 . Probabilistic Forecasting of the Masses and Radii of Other Worlds . The Astrophysical Journal . 834 . 1 . 17 . 10.3847/1538-4357/834/1/17. 1603.08614 . 119114880 . free .