greatitself.
Continue
advanced35 min studyLesson 896 in Physics

Schwarzschild convective stability

Derive local parcel stability and distinguish temperature-gradient instability from composition stabilization and quantitative mixing.

Before you begin
On this page

The question for today

What you’ll understand

  • Explain the sign of buoyancy for a pressure-balanced adiabatic displacement.
  • Apply Schwarzschild and Ledoux criteria in their declared composition domains.
  • Compute a buoyancy frequency without identifying it with a global stellar mode.

Convection begins when a small displacement makes a parcel continue moving away from its original position. The instability is controlled by the parcel’s density relative to its surroundings after it has adjusted its pressure. We assume a local, nearly hydrostatic background, small displacements, fast sound-pressure adjustment, and negligible parcel heat or composition exchange during that displacement. Rotation, magnetic tension, viscosity, and differential thermal or particle diffusion are omitted from this first test. The criterion predicts local onset, rather than the eventual speed, flux, or extent of turbulent mixing.

Follow a parcel and compare it with its environment

Pressure decreases outward, so an upward-moving parcel expands to stay close to the surrounding pressure. Its temperature follows an adiabatic response if little heat is exchanged during the motion. The environment follows the actual background temperature profile instead. With homogeneous composition and ordinary positive thermal expansion, a parcel that stays warmer at equal pressure is less dense and feels upward buoyancy. If the environmental temperature decreases outward too steeply relative to the adiabatic response, the displacement therefore amplifies.

Define the actual gradient as logarithmic temperature change with logarithmic pressure through the background. Define the adiabatic gradient through a material parcel at constant entropy and composition. These are different derivatives even though both are dimensionless. For homogeneous composition the actual profile is locally unstable when its gradient exceeds the adiabatic value. Testing an initially radiative layer replaces the actual gradient by the full-radiative-flux gradient. After convection develops, its actual gradient usually changes, so the initial radiative diagnostic must not be mistaken for the final stratification.

∇=dln⁡Tdln⁡P,∇ad=(∂ln⁡T∂ln⁡P)s,Xi,∇rad>∇ad(homogeneous radiative layer)\nabla=\frac{d\ln T}{d\ln P},\qquad \nabla_{\rm ad}=\left(\frac{\partial\ln T}{\partial\ln P}\right)_{s,X_i},\qquad \nabla_{\rm rad}>\nabla_{\rm ad}\quad\text{(homogeneous radiative layer)}
The inequality is the Schwarzschild onset criterion under the stated adiabatic, composition-homogeneous parcel assumptions. It tests whether radiation alone would require an unstable stratification. Equality is marginal in this ideal local analysis, not a prediction of a finite mixing velocity.↔ Scroll wide equations horizontally.

Composition changes the density comparison

A parcel carries its composition during an adiabatic displacement. If the mean molecular weight increases inward, an upward parcel retains a larger molecular weight than its new surroundings. This can keep it denser even if it is somewhat warmer. The Ledoux criterion includes that stabilizing density contrast through thermodynamic derivatives of the equation of state. For a simple ideal gas both relevant coefficients are one; a gas–radiation mixture or a partially degenerate medium need not share those values.

A region can be unstable under the composition-free Schwarzschild test yet stable under Ledoux when an inward molecular-weight increase supplies sufficient restoring buoyancy. That does not settle all long-term behavior. Thermal diffusion can act faster than particle diffusion, creating slower instabilities and semiconvective or fingering transport under appropriate gradients. Their development needs additional timescales and physics. Likewise, a local stable boundary does not halt every moving convective parcel immediately; inertia and nonlocal flows can produce overshoot beyond the simple onset location.

δ=−(∂ln⁡ρ∂ln⁡T)P,μ,φ=(∂ln⁡ρ∂ln⁡μ)P,T,∇μ=dln⁡μdln⁡P,∇rad>∇ad+φδ∇μ\delta=-\left(\frac{\partial\ln\rho}{\partial\ln T}\right)_{P,\mu},\quad\varphi=\left(\frac{\partial\ln\rho}{\partial\ln\mu}\right)_{P,T},\quad\nabla_\mu=\frac{d\ln\mu}{d\ln P},\quad\nabla_{\rm rad}>\nabla_{\rm ad}+\frac{\varphi}{\delta}\nabla_\mu
The final inequality is Ledoux onset for δ>0 when composition is adequately parameterized by μ. General abundance dependence can require a species-wise expression. Positive ∇μ means μ increases inward because pressure also increases inward. The derivatives must come from the chosen equation of state.↔ Scroll wide equations horizontally.

Restoring buoyancy has its own frequency scale

Linearize the density contrast for a small displacement and multiply by gravitational acceleration. The resulting restoring acceleration is minus the squared buoyancy frequency times displacement. A positive squared frequency gives oscillatory restoring behavior in this local idealization; a negative value gives exponential growth. The expression uses the actual background gradient and includes the thermal-expansion factor outside the bracket. Leaving out that factor is harmless only for a declared ideal gas with delta equal to one, not for a general equation of state.

The pressure scale height is positive when defined as minus the inverse logarithmic pressure derivative. It connects local pressure stratification with buoyancy and the equation of state. A measured global gravity-mode spectrum depends on the whole buoyancy profile, turning points, and boundary conditions, rather than a single local frequency. Convective growth rates, diffusion, and mode damping require further analysis. The positive-stiffness oscillator instrument can illustrate restoring oscillations, but cannot represent a negative squared buoyancy frequency or simulate nonlinear overturning.

N2=gδHP(∇ad−∇+φδ∇μ),HP=−(dln⁡Pdr)−1,ξ¨=−N2ξN^2=\frac{g\delta}{H_P}\left(\nabla_{\rm ad}-\nabla+\frac{\varphi}{\delta}\nabla_\mu\right),\qquad H_P=-\left(\frac{d\ln P}{dr}\right)^{-1},\qquad \ddot\xi=-N^2\xi
g is the positive inward acceleration in m s⁻², HP is m, and N² is s⁻². ξ is a small outward displacement. This local adiabatic expression omits rotation, magnetic forces and diffusion; it uses the actual ∇, with ∇=∇rad only for a fully radiative background.↔ Scroll wide equations horizontally.

Work it through

A composition gradient changes the verdict

An initially radiative ideal-gas layer has ∇rad=0.45, ∇ad=0.40, ∇μ=0.10, and δ=φ=1. Compare the Schwarzschild and Ledoux tests. With g/HP=10⁻⁶ s⁻², compute N² for that radiative background.

  1. Ignoring composition gives 0.45>0.400.45>0.40, so the Schwarzschild test labels the proposed radiative profile unstable.

  2. The Ledoux threshold is 0.40+0.10=0.50. Since 0.45 is below it, the composition gradient stabilizes adiabatic parcel displacements.

  3. The restoring bracket is 0.40−0.45+0.10=0.050.40-0.45+0.10=0.05, giving N²=5 × 10⁻⁸ s⁻². Its positive sign agrees with the Ledoux verdict.

The result

Schwarzschild alone predicts instability, but Ledoux gives stable adiabatic buoyancy with N²=5 × 10⁻⁸ s⁻²; slower diffusive processes remain outside this test.

The living notebook

Explore the model

Interactive lab

Interactive lab · oscillator

Restoring frequency controls a stable oscillation

Reset and set BOTH Damping ratio ζ to 0 and Natural frequency ω₀ to 1; predict period 2π in the instrument’s time units. Then Reset and set BOTH ζ to 0 and ω₀ to 2; predict period π. Play or Step to compare the displacement and velocity. This fixed positive-stiffness, unit-mass oscillator illustrates a stable buoyancy analogue; it has no negative N², composition gradient, convective flux, or nonlinear overturning.

Restoring frequency controls a stable oscillationUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA positive restoring coefficient gives a frequency scale, whereas an instability requires a different sign and additional dynamics.

Pause & reason

Check your understanding

What can an inward-increasing molecular-weight profile do in the Ledoux criterion?

Your turn

Try it yourself

In a different homogeneous ideal-gas layer use actual ∇=0.30, ∇ad=0.40, δ=1, and g/HP=4 × 10⁻⁶ s⁻². Find N², N, and the period of the corresponding ideal local buoyancy oscillation.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Evaluate the gradient difference first. The oscillation period is two pi divided by the positive angular frequency.

Compare with the solution

The difference is 0.10, so N2=4×10−7 s−2N^2=4\times10^{-7}\ {\rm s^{-2}} and N=2√10 × 10⁻⁴ s⁻¹, about 6.32 × 10⁻⁴ s⁻¹. The ideal period is about 9.93 × 10³ s. This is a local adiabatic parcel period, not the period of an arbitrary observed global mode or a prediction of a turbulent convective turnover time.

Keep the idea

Take it with you

  • Schwarzschild compares radiative and adiabatic gradients for homogeneous composition.
  • Ledoux retains the density effect of a parcel’s conserved composition.
  • Positive N² gives restoring buoyancy, while transport and nonlocal mixing require further physics.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.