ASTR 474 · Computational Astrophysics · Lecture 07

Lecture 07: Radiative Transfer as a Numerical Problem

Which quantities must be specified to solve the transfer equation along a ray?

14-week course · 3 credits · Week 7

Motivation

Which quantities must be specified to solve the transfer equation along a ray?

Specific intensity changes through emission and absorption. Optical depth is dimensionless and organizes attenuation; a formal solution is available for constant source function, providing a stringent numerical check.

Today’s result must be interpretable as a physical statement with declared data, algorithm, uncertainty, and a check that could fail.

Learning objectives

After this lecture, you can

Core language

Terms for This Lecture

Specific intensity

Definition: Radiative energy transported per projected area, time, solid angle, and frequency.

Why it matters here: It helps answer Which quantities must be specified to solve the transfer equation along a ray? within the numerical model described here.

Relationships: Connect it to Optical depth, the governing expression \(\frac{dI_\nu}{d\tau_\nu}=S_\nu-I_\nu,\quad I(\tau)=I_0e^{-\tau}+S(1-e^{-\tau})\), and the computed diagnostic developed later in the lecture.

Boundary or common confusion: The grey constant-source slab is an approximation. Frequency coupling, scattering source functions, geometry, and non-LTE populations can change the equation being solved.

Physical and mathematical model

State what is being evolved

Specific intensity changes through emission and absorption. Optical depth is dimensionless and organizes attenuation; a formal solution is available for constant source function, providing a stringent numerical check.

\[\frac{dI_\nu}{d\tau_\nu}=S_\nu-I_\nu,\quad I(\tau)=I_0e^{-\tau}+S(1-e^{-\tau})\]

Define the state, coordinates, units, initial/boundary conditions, and domain before interpreting solver output.

Numerical method

Algorithm choice is a modeling decision

State ray direction and boundary intensity; express opacity and emissivity in consistent units; integrate with an adaptive solver; compare with the formal solution at thin and thick optical depth.

  1. Write the model and its assumptions.
  2. Choose a solver and identify its stability, accuracy, and conservation properties.
  3. Record configuration, tolerances or timestep, and output sampling.
  4. Compare against an analytic limit, independent method, or controlled refinement.

Visual reasoning

Formal transfer solution across optical depth

Intensity approaches the source functionA computed formal solution starts at the incident boundary intensity and asymptotes toward a constant source function as optical depth increases. source function S=0.7optical depth tauspecific intensity I boundary I0=0.2
The source-function asymptote and boundary intensity are parameters of the synthetic constant-source slab.

Computed example

One result, with a check

At tau=1.0, the exact intensity is 0.5160603; SciPy DOP853 differs by 9.48e-14.

Inputs and method are generated by src/astr474_computations.py; the machine-readable values are in data/generated/computed_results.json.

The grey constant-source slab is an approximation. Frequency coupling, scattering source functions, geometry, and non-LTE populations can change the equation being solved.

Verification

What would falsify this result?

Do not report a smoother curve as proof of convergence.

Failure modes

Numerical agreement has a domain

The grey constant-source slab is an approximation. Frequency coupling, scattering source functions, geometry, and non-LTE populations can change the equation being solved.

Solver error, discretization error, stochastic error, parameter uncertainty, and model discrepancy are different quantities. Diagnose and report them separately.

Active reasoning

Make a prediction before rerunning

Which quantities must be specified to solve the transfer equation along a ray?

State a measurable outcome, a competing explanation, and the smallest computation that distinguishes them. Include an expected direction or scaling, then test it.

Workflow

From parameters to auditable output

State ray direction and boundary intensity; express opacity and emissivity in consistent units; integrate with an adaptive solver; compare with the formal solution at thin and thick optical depth.

Record the exact command, package versions, data path, random seed if applicable, and output checksum. Keep raw inputs unchanged and generated files separate.

Transfer

Lab and problem-set connection

Lab 05 compares SciPy DOP853 with the analytic slab solution and propagates opacity/temperature perturbations.

Submit plots with labeled axes/units, a concise uncertainty argument, and enough code or intermediate values for a peer to reproduce the claim.

Synthesis

Three durable checks

  1. Does the implemented equation match the intended physical model?
  2. Does the numerical result converge under a controlled method change?
  3. Does the interpretation stay inside the validity and provenance limits?

The grey constant-source slab is an approximation. Frequency coupling, scattering source functions, geometry, and non-LTE populations can change the equation being solved.

Sources

Technical references

Data and provenance details are recorded in reference-log.md. Last updated 2026-09-27.