Radiation, convection and internal structure

Producing energy and moving it are different processes

Energy deposited in stellar material must be transported or stored as part of the energy budget. Radiation transfers energy through the emission, absorption and scattering of light. Convection transports energy through bulk fluid motion, usually with rising and returning material rather than every particle steadily moving outward.

Energy flux F is power passing through unit area, measured in watts per square metre. Luminosity L is power passing through a surrounding surface. For a spherical description with uniform radial flux on that surface, L(r) = 4πr²F(r).

If radiation and convection both contribute, the radiative luminosity is only the radiation-carried part. Replacing it by the total luminosity in a radiative-only relation can give an incorrect temperature gradient. Transport channels and the local energy budget must be accounted for consistently.

A stellar model therefore combines mechanical support with energy generation, losses and transport. The arrangement of transport regions depends on the physical state and evolution; the layout of the present Sun is not a universal diagram for every star.

A diffusion formula has a definite regime

Opacity describes how matter impedes radiation in a specified transport description. In a mass-opacity convention, κ has units m² kg⁻¹. With density ρ, an approximate local interaction length is ℓ = 1/(κρ), provided the chosen opacity is appropriate to that estimate.

Deep in sufficiently opaque material, a local diffusion description can apply when the interaction length is much smaller than the scale over which temperature changes. Local thermodynamic equilibrium means the matter and radiation have the appropriate local thermal relation in that approximation; it does not mean that the temperature is identical throughout the star.

Using the appropriate Rosseland mean opacity κ_R, a radiative diffusion relation is F_rad = −[16σT³/(3κ_Rρ)] dT/dr. Temperature decreasing outward gives positive outward flux. The formula does not describe transparent free streaming or the detailed formation of a photospheric spectrum.

The mean opacity compresses wavelength-dependent interactions into a transport coefficient for this regime. It is not interchangeable with one measured absorption-line depth or an arbitrary telescope-band extinction. Near a surface, a more complete radiative-transfer treatment is needed.

The same opacity change poses two different questions

Construct an original local diffusion case with T = 10⁶ K, ρ = 1000 kg m⁻³, κ_R = 0.02 m² kg⁻¹ and dT/dr = −1 K m⁻¹. Using σ = 5.670374419 × 10⁻⁸ W m⁻² K⁻⁴ at the displayed precision gives F_rad ≈ 1.51210 × 10¹⁰ W m⁻².

Here the interaction-length estimate is 0.05 m, while the temperature scale T/|dT/dr| is 10⁶ m. Their large separation supports the stipulated local optically thick regime. These numbers are not a measured stellar layer or a complete self-consistent stellar structure.

If opacity doubles while T, ρ and the gradient are held fixed, the formula gives half the flux: about 7.56050 × 10⁹ W m⁻². If flux is instead required to remain fixed while the same T and ρ are retained, the gradient must become twice as steep, −2 K m⁻¹.

Those are different conditional comparisons. A real star can adjust its structure or redistribute energy among channels. The algebra alone does not predict that doubling some opacity will halve the entire star’s luminosity.

A rising parcel can cool yet remain warmer than its new surroundings

Consider an original parcel model with the same composition in parcel and environment, ideal-gas behaviour and ratio of heat capacities γ = 5/3. The parcel approximately matches the surrounding pressure as it moves. Adiabatic means it exchanges negligible heat during the displacement, while expansion can still change its temperature.

Let q = P/P₀, where P₀ and T₀ are the common starting pressure and temperature. Rising to lower ambient pressure means q decreases. The adiabatic parcel has T_parcel/T₀ = q^0.4. Compare two stipulated environments: T_stable/T₀ = q^0.25 and T_unstable/T₀ = q^0.60.

At equal pressure and composition, ideal-gas density is inversely proportional to temperature. Relative to the common starting density ρ₀, each density ratio is q divided by its temperature ratio. The complete table supplies all five pressure ratios and the associated temperature and density comparisons.

At q = 0.8, the parcel has cooled to about 0.914610 T₀. It is cooler and about 1.034038 times as dense as the stable environment, tending to sink back. It is warmer and about 0.956352 times as dense as the unstable environment, tending to continue rising. Buoyancy depends on the comparison with the new surroundings, not simply on whether the parcel cooled.

An instability criterion does not determine the carried power

In the compositionally homogeneous limit, compare the environment’s logarithmic gradient ∇ = d ln T/d ln P with the parcel’s adiabatic gradient. The supplied ideal gas has ∇_ad = (γ − 1)/γ = 0.4. A steeper environmental gradient can make an upward-displaced parcel warmer and less dense, giving the buoyant instability illustrated here.

This comparison addresses the onset of an instability under its assumptions. It does not calculate convective velocity, turbulent mixing efficiency or energy flux. A complete model needs a transport prescription and boundary treatment as well as the local stability question.

A composition gradient can change the density comparison and require a broader stability treatment. Rotation, magnetic stresses, heat exchange during motion and departures from the chosen gas model can also matter. The two supplied power-law environments isolate one mechanism rather than covering every stellar situation.

The present standard solar description has a radiative interior below an outer convective envelope. Other stars and evolutionary stages can have different arrangements. Observations and model tests, introduced in the following lessons, help determine where a proposed structure succeeds or needs improvement.

Worked examples

Try the reasoning

1. With κ = 0.02 m² kg⁻¹ and ρ = 1000 kg m⁻³, what local interaction-length estimate follows?

   Show the answer and reasoning

   0.05 m.

   1/(κρ) = 1/(0.02 × 1000) = 1/20 m, using the stated effective-opacity approximation.

2. In the supplied diffusion model, double opacity while keeping temperature, density and gradient fixed. By what factor does flux change?

   Show the answer and reasoning

   It is multiplied by 1/2.

   The radiative flux is inversely proportional to opacity under those held conditions.

3. Does adiabatic upward displacement imply that the parcel’s temperature remains unchanged?

   Show the answer and reasoning

   No.

   Negligible heat exchange does not prevent expansion cooling. In the supplied model its temperature follows q^0.4.

4. At the same pressure and composition in the supplied ideal-gas model, which is less dense: the warmer or cooler parcel?

   Show the answer and reasoning

   The warmer parcel.

   Density is proportional to P/T for fixed composition in this model.

5. Does finding a buoyant instability alone determine how many watts convection transports?

   Show the answer and reasoning

   No.

   A transport prescription and additional physical information are needed beyond the onset criterion.

What to carry forward

• Radiative diffusion requires an appropriate opacity and a sufficiently optically thick local regime.
• Held conditions determine whether an opacity comparison changes the calculated flux or required gradient.
• A parcel’s buoyancy depends on its density relative to the surrounding material, with composition and thermodynamic assumptions stated.