advanced35 min studyLesson 893 in Physics

Radiative and convective energy transport

Derive optically thick radiative diffusion and identify when an additional convective or conductive flux closure is required.

Before you begin

On this page

The question for today

What you’ll understand

Connect local luminosity, flux, opacity and the radiative temperature gradient.
State the LTE and optical-depth assumptions behind Rosseland diffusion.
Separate an instability criterion from a quantitative convective-flux prescription.

Energy generated near a stellar center must cross surrounding matter before emerging as light. The transport mechanism changes the temperature profile needed to carry that energy, and therefore changes pressure support and reaction rates. Radiation, electron conduction, and bulk fluid motion can all contribute, with different domains of validity. We derive radiative diffusion in an optically thick region with local thermalization and then distinguish its predicted gradient from the actual gradient when convection carries part of the flux. An atmosphere requires separate treatment as photons begin to escape freely.

Diffusion converts a temperature gradient into flux

In a locally thermal radiation field, photon energy density is proportional to the fourth power of temperature. Photons travel short mean free paths between interactions and diffuse down the radiation-energy gradient. Angular averaging supplies a factor of one third in the diffusion coefficient. Combining this with the derivative of radiation energy density yields a flux opposite to the temperature gradient. A negative outward temperature gradient therefore produces positive outward radiative flux, consistently with our luminosity convention.

The opacity here is the Rosseland mean per unit mass, including the appropriate transport effect of absorption and scattering. It is weighted toward relatively transparent frequency windows under the diffusion assumptions. It is not generally the same mean used for optically thin emission or radiative acceleration. Density, composition, and ionization strongly affect opacity; a constant value is a declared idealization, and a simple power-law opacity cannot represent every ionization feature. The photon mean free path must remain short relative to local thermodynamic scale lengths.

F_{\rm rad}=-\frac{4acT^3}{3\kappa_R\rho}\frac{dT}{dr},\qquad \ell_\gamma=\frac1{\kappa_R\rho},\qquad L=4\pi r^2F_{\rm tot}

In SI, κR is m² kg⁻¹, ℓγ is m, and flux is W m⁻². a is the radiation energy-density constant and c the light speed. Diffusion assumes near-isotropy, large optical depth over a scale height, and local thermalization sufficient for the thermal radiation source.↔ Scroll wide equations horizontally.

Combine transport with mechanical balance

If radiation alone must carry the entire specified luminosity, solve the flux relation for the required temperature gradient. Dividing it by hydrostatic balance gives the dimensionless radiative logarithmic gradient. This quantity depends on local pressure, temperature, opacity, enclosed mass, and luminosity. It remains useful as a diagnostic even in a convective region: there it describes the hypothetical gradient needed if radiation carried all the flux, rather than automatically describing the actual stellar stratification.

A larger opacity or luminosity makes radiative transport demand a steeper gradient, all else fixed. This can provoke buoyant instability, depending also on the equation of state and composition gradient. The expression is interpreted by its regular central limit rather than by substituting zero mass and radius independently at the origin. Near the photosphere the diffusion approximation breaks down, so an atmospheric transfer calculation or a specified matching prescription replaces a blind extension of the interior formula to zero optical depth.

\left(\frac{dT}{dr}\right)_{\rm rad}=-\frac{3\kappa_R\rho L}{16\pi ac r^2T^3},\qquad \nabla_{\rm rad}=\frac{3\kappa_RLP}{16\pi acGmT^4}

∇rad means d ln T/d ln P required for radiation to carry all L. The formula combines spherical hydrostatic balance with optically thick radiative diffusion. It is dimensionless; m and L are the local enclosed mass and transported luminosity, not necessarily their surface values.↔ Scroll wide equations horizontally.

A transport law must say how much flux is carried

Convection transports enthalpy through correlated upward and downward motions. Equal mean upward and downward mass flows can carry a net outward heat flux because they have different thermal contents. A local instability criterion identifies whether small displaced parcels grow, but does not determine the resulting flux. A mixing-length prescription adds a model for travel distance, heat exchange, and velocity; its parameters and treatment of boundaries influence the predicted stellar radius, mixing, and lifetime.

Efficient deep convection often keeps the actual gradient close to its adiabatic value, while inefficient surface convection can require substantial superadiabaticity. Electron conduction can dominate energy transport in dense degenerate regions. Diffusive radiative and conductive coefficients add when they multiply the same temperature gradient, whereas convective flux requires its own closure. A static spherical background is an average description of this multidimensional motion. It should not be interpreted as proof that every parcel is motionless or that one universal coefficient predicts all stellar heat transport.

F_{\rm tot}=-(K_{\rm rad}+K_{\rm cond})\frac{dT}{dr}+F_{\rm conv},\qquad K_{\rm rad}=\frac{4acT^3}{3\kappa_R\rho}

Krad and Kcond have SI units W m⁻¹ K⁻¹, while Fconv is W m⁻². This scalar local form assumes isotropic diffusion. It separates a constitutive diffusion coefficient from a convective flux requiring additional dynamical or phenomenological information.↔ Scroll wide equations horizontally.

Work it through

A local radiative-flux toy

At one optically thick layer take an already evaluated Krad=2 × 10¹⁰ W m⁻¹ K⁻¹, negligible conduction and convection, radius r=10⁸ m, and dT/dr=−10⁻² K m⁻¹. Find outward flux and luminosity. At fixed local Krad, predict the gradient required if luminosity doubles.

The flux is minus coefficient times gradient: $F=2\times10^8\ {\rm W\,m^{-2}}$ . Its positive sign indicates outward flow.
Multiplying by spherical area gives $L=8\pi\times10^{24}\ {\rm W}$ , approximately 2.51 × 10²⁵ W.
At the same radius and coefficient, doubling luminosity doubles flux. The required gradient becomes −2 × 10⁻² K m⁻¹, provided the diffusion assumptions and stated fixed-coefficient approximation remain valid.

The result

The layer carries 2 × 10⁸ W m⁻² and about 2.51 × 10²⁵ W; twice that luminosity requires twice the gradient magnitude in this local toy.

The living notebook

Explore the model

Interactive lab

Interactive lab · heat

Diffusion damps shorter thermal structure faster

Reset and set BOTH Spatial mode n to 1 and Diffusivity κ to 0.5; after one Step (t=0.5), predict amplitude exp(−0.25). Then Reset and set BOTH n to 2 and κ to 0.5; after one Step predict exp(−1). The amplitude and L² norm readouts test the n² diffusion rate. This is one analytic periodic one-dimensional mode, without stellar geometry, varying opacity, radiation transport, nuclear heating, or convection.

Current temperature modeInitial sineCurrent amplitude

Current amplitude: 1
Spatial L² norm: 1.772
Decay rate · κn²: 3.6 /s

Spatial mode n3Diffusivity κ0.4

t = 0 s

The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. 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 forDiffusion has a scale-dependent relaxation rate, while a stellar profile also needs sources, geometry, and closure.

Pause & reason

Check your understanding

Interactive checks need JavaScript. You can still read the explanation.

Read the checkpoint explanation

∇rad is a hypothetical full-radiative-flux gradient. The actual gradient is determined by the combined transport, whereas ∇ad and ∇μ describe different quantities.

Your turn

Try it yourself

In a separate local diffusion toy, Krad=3 × 10⁹ and Kcond=1 × 10⁹ W m⁻¹ K⁻¹. The total outward flux is 10⁸ W m⁻², with no convection. Find dT/dr and the fraction carried by conduction.

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

Give me a hint

Add the two coefficients before solving for the common gradient, then compute the conductive contribution alone.

Compare with the solution

The total coefficient is 4 × 10⁹ W m⁻¹ K⁻¹, giving $dT/dr=-0.025\ {\rm K\,m^{-1}}$ . Conduction carries 2.5 × 10⁷ W m⁻², one quarter of the total; radiation carries the remaining three quarters. Adding reciprocal opacities can express the same operation if an effective conductive opacity is defined consistently, but adding the raw opacities would give the wrong transport coefficient.

Keep the idea

Take it with you

Rosseland diffusion needs optical thickness and a locally thermal radiation field.
The radiative gradient diagnoses what radiation alone would require at the specified luminosity.
Convection needs a flux closure, and electron conduction supplies a separate diffusive channel.

Read further

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

Pols, 5.2–5.3 and 5.5.1–5.5.4: diffusion, opacity, conduction, convective transport and mixing

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

Lesson not marked complete.