Accretion requires an angular momentum budget
Derive a steady disk surface density from torque transport and an explicit inner boundary.
On this page
The question for today
What you’ll understand
- Distinguish inward mass transport from outward angular momentum transport.
- Derive the steady zero-torque boundary factor.
- Explain what an effective viscosity and an alpha prescription assume.
Gravity supplies an energy reservoir, but gravity alone does not let rotating gas spiral steadily inward. A parcel on a smaller circular orbit needs less specific angular momentum. A disk must therefore redistribute angular momentum internally or lose it through its surfaces. We will derive the simplest steady internal-transport model, while keeping its boundary condition separate from the physical mechanism that produces the stress.
An inward flow can carry angular momentum outward
Use cylindrical radius r around a Newtonian point mass M, neglecting disk self-gravity. A geometrically thin disk has sound speed much smaller than orbital speed, so radial pressure support is a small correction and the orbital angular velocity is approximately Keplerian. The specific angular momentum l=√(GMr) increases outward even though the angular velocity decreases. Thus a parcel moving inward must shed angular momentum. Dissipating orbital energy without a corresponding torque does not solve this constraint.
Let Σ be mass per disk-plane area, integrated vertically, and let vᵣ be negative for inward motion. Define the positive inward mass rate as minus circumference times surface density times radial velocity. An effective kinematic viscosity ν represents a vertically averaged internal stress; it need not be molecular viscosity. In a Keplerian shear, the resulting positive outward torque is proportional to νΣl. Torque has units of angular momentum per time, equivalently energy, but torque transport and local heat generation remain different physical quantities.
The inner boundary fixes the integration constant
In a steady disk without mass loss, the mass rate is independent of radius. Conservation makes the net inward angular momentum flux, advective transport minus the outward torque, constant. If the internal torque vanishes at an inner radius Rᵢₙ, that constant is the mass rate times the angular momentum carried through the inner edge. Substituting the Keplerian torque gives the familiar boundary factor. Far outside the edge, the factor approaches one; close to the edge it cannot be dropped.
The zero-torque condition is a model assumption, rather than a consequence of circular motion. A magnetic connection to a star, appreciable stress through a black-hole inner region, or a disk wind changes the angular momentum balance. Even in the chosen model, extrapolating the thin steady radial-velocity expression all the way to its inner endpoint is unsafe: the assumed slow, nearly circular flow eventually requires an inner transition. A formal divergence from vanishing surface density is not a measured infinite speed.
A closure for stress is not its microscopic explanation
An often useful convention is ν=αcₛH with vertical scale height H=cₛ/Ω and dimensionless α. We use precisely that convention; alternative stress-based definitions can carry factors such as two thirds. It gives a redistribution estimate tᵥᵢₛc≈r²/ν. Molecular collisions are generally too slow to explain rapidly evolving astrophysical disks. In sufficiently ionized disks, magnetic tension and differential rotation can generate magnetorotational instability and correlated Maxwell and Reynolds stresses. Its operation depends on field strength, conductivity, and available wavelengths, so specifying α does not prove a particular turbulence mechanism.
A simple reservoir comparison isolates the mass budget. In the queue lab, Reset and set both rates to arrival three and service three, then Step twice: the fluid backlog stays four. Reset and set arrival three and service two, then Step twice: it becomes five. Interpret this only as inflow minus outflow into an annulus. The instrument has no angular momentum or stress variable, and its separately labeled stochastic residence-time formula is not a disk accretion time.
Work it through
Find a disk surface density and inward velocity
For a hypothetical steady disk with zero inner torque, take Rᵢₙ=10⁸ m, r=4Rᵢₙ, ν=10¹⁰ m²/s and Ṁ=6π×10¹⁴ kg/s. Within the stated Newtonian thin-disk model, find the boundary factor, Σ and vᵣ.
The boundary factor is 1−√(Rᵢₙ/r)=1−1/2=1/2. This retains the inner-edge angular momentum flux.
The product νΣ is (6π×10¹⁴)/(3π) times 1/2=10¹⁴ kg/s. Dividing by ν gives Σ=10⁴ kg/m².
Mass conservation gives vᵣ=−Ṁ/(2πrΣ)=−(6π×10¹⁴)/(2π×4×10⁸×10⁴)=−75 m/s. Its sign agrees with the inward-rate convention.
The factor is 1/2, Σ=10⁴ kg/m², and vᵣ=−75 m/s. These are conditional model values, not observed disk measurements.
The living notebook
Explore the model
Interactive lab · queue
Separate a mass reservoir from its torque mechanism
Reset and set BOTH rates to λ=3,μ=3; Step twice to t=1: predict fluid backlog 4. Reset and set λ=3,μ=2; Step twice: predict 5. This inflow–outflow analogy contains no disk stress or angular momentum, and the M/M/1 residence-time readout is not a viscous timescale.
- Fluid backlog
- 4
- Net rate · λ − μ
- −1 /s
- M/M/1 mean residence
- 1 s
Fluid backlog starts at 4 and follows max(0, 4 + (λ − μ)t). The separate M/M/1 residence-time formula assumes Poisson arrivals, exponential service and λ < μ. 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 steady mass reservoir still needs an independent angular momentum transport mechanism.
Pause & reason
Check your understanding
Your turn
Try it yourself
Adopt ν=αcₛH and H=cₛ/Ω with α=0.1 and H/r=0.05. Estimate tᵥᵢₛcΩ and the number of orbital periods in tᵥᵢₛc. Use tᵥᵢₛc≈r²/ν.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Substitute ν=αH²Ω. One orbital period is 2π/Ω, not 1/Ω.
Compare with the solution
The dimensionless estimate is tᵥᵢₛcΩ≈(r/H)²/α=20²/0.1=4000. Dividing by 2π gives about 637 orbits. The estimate expresses a separation of transport and orbital times in this alpha convention; it is not an exact draining time for every disk boundary condition.
Keep the idea
Take it with you
- Accreting circular gas must lose specific angular momentum.
- The inner torque fixes the steady surface-density boundary factor.
- An alpha viscosity closes averaged equations without specifying all plasma microphysics.
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.