Differences

This shows you the differences between two versions of the page.

--- un:equations-of-stellar-structure [2025/10/27 00:04] – created asad
+++ un:equations-of-stellar-structure [2025/10/27 00:33] (current) – [Primary and secondary] asad
@@ Line 5: / Line 5: @@
 These equations express conservation of mass, momentum, and energy, together with the laws of heat transfer and the thermodynamic state of stellar matter.
-===== 1. Hydrostatic equilibrium =====
+Of these seven, the **first five** are known as the **primary equations**, because they are **differential relations** that describe how the star’s physical variables vary with radius.
+The **last two** are called the **secondary equations**, because they provide the **thermodynamic and compositional relations** (the *equations of state*) needed to close the system.
+===== - Hydrostatic equilibrium =====
 The inward pull of gravity is balanced by the outward pressure gradient at every layer inside the star:
 $$
-\frac{dP}{dr} = -\rho(r)\, g(r) = -\rho(r)\, \frac{G M(r)}{r^2}.
+\frac{dP}{dr} = -\rho(r)\, g(r) = -\rho(r)\,\frac{G M(r)}{r^2}.
 $$
-Here \(P(r)\) is the local pressure, \(\rho(r)\) the density, \(M(r)\) the enclosed mass, and \(G\) the gravitational constant.
+Here \(P(r)\) is the local pressure, \(\rho(r)\) the density, \(M(r)\) the enclosed mass, \(g(r)=GM(r)/r^2\) the local gravitational acceleration, and \(G\) the gravitational constant.
 This equation ensures that the net force on a small mass element is zero in a stable star.
+It expresses **mechanical equilibrium**, the balance between gravity and pressure.
-===== 2. Conservation of mass =====
+===== - Conservation of mass =====
-The rate at which mass increases with radius equals the density times the area of the spherical shell:
+**Statement.**
 $$
@@ Line 24: / Line 28: @@
 $$
-Integrating from the center to the surface gives the total stellar mass \(M(R)\).
+**Derivation.**
+Consider a thin spherical shell at radius \(r\) with thickness \(dr\).
+The shell’s volume is
+$$
+dV = 4\pi r^2\,dr,
+$$
+so the shell’s mass is
+$$
+dM = \rho(r)\,dV = \rho(r)\,4\pi r^2\,dr.
+$$
+Dividing by \(dr\) gives the radial mass accumulation rate
+$$
+\frac{dM}{dr} = 4\pi r^2 \rho(r),
+$$
+which is the local form of **mass continuity** in spherical symmetry.
+Integrating from the center yields \(M(r)=\int_0^r 4\pi r'^2 \rho(r')\,dr'\), giving the total mass \(M(R)\) at the surface.
-===== 3. Conservation of energy (luminosity distribution) =====
+===== - Conservation of energy (luminosity distribution) =====
-The energy generated inside a shell of thickness \(dr\) adds to the outward luminosity \(L(r)\):
+**Statement.**
 $$
-\frac{dL}{dr} = 4\pi r^2 \rho(r)\, \epsilon(\rho,T),
+\frac{dL}{dr} = 4\pi r^2 \rho(r)\,\epsilon(\rho,T),
 $$
-where \(\epsilon(\rho,T)\) is the **energy generation rate** (in W kg\(^{-1}\)), mainly from nuclear reactions.
+where \(L(r)\) is the luminosity passing outward through radius \(r\), and \(\epsilon(\rho,T)\) is the **energy generation rate per unit mass** (W kg\(^{-1}\)) from nuclear reactions (composition dependence implied: \(\epsilon=\epsilon(\rho,T,\mathrm{C})\)).
-At the stellar surface \(L(R)\) equals the total luminosity of the star.
-===== 4. Temperature gradient: radiative transfer =====
+**Derivation.**
+In steady state, the net increase of luminosity across a thin shell equals the **power generated** within that shell.
+For the same shell \(dV=4\pi r^2 dr\), the mass is \(dM=\rho(r)\,dV\).
+The energy produced per unit time in the shell is
+$$
+d\dot{E} = \epsilon(\rho,T)\,dM = \epsilon(\rho,T)\,\rho(r)\,4\pi r^2 dr.
+$$
+By energy conservation, this equals the outward increase in luminosity,
+$$
+dL = d\dot{E} \;\Rightarrow\; \frac{dL}{dr} = 4\pi r^2 \rho(r)\,\epsilon(\rho,T).
+$$
+This equation describes how luminosity builds up with radius as nuclear energy is released inside the star.
+===== - Temperature gradient: radiative transfer =====
 In regions where energy is transported primarily by radiation, the temperature decreases outward according to
 $$
-\frac{dT}{dr} = -\,\frac{3\kappa\rho L}{16\pi a c T^3 r^2},
+\frac{dT}{dr} = -\,\frac{3\,\kappa(\rho,T)\,\rho\, L}{16\pi a c\, T^3\, r^2},
 $$
-where \(\kappa\) is the **opacity** (m\(^2\) kg\(^{-1}\)), \(a\) is the radiation constant, and \(c\) the speed of light.
+where \(\kappa\) is the **opacity** (m\(^2\) kg\(^{-1}\)), \(a\) the radiation constant, and \(c\) the speed of light.
 A higher luminosity or greater opacity requires a steeper temperature gradient to carry the same energy flux.
-(This is the full radiative transfer equation—its derivation and detailed discussion are given separately in [[un:heat-transfer-in-stars|Heat transfer in stars]].)
+(We **omit the detailed derivation** here; see [[un:heat-transfer-in-stars|Heat transfer in stars]] for a full treatment.)
-===== 5. Temperature gradient: convective transfer =====
+===== - Temperature gradient: convective transfer =====
 When the radiative gradient exceeds the adiabatic one, convection becomes efficient.
-The **adiabatic temperature gradient** is approximately
+The **adiabatic temperature gradient** can be written approximately as
 $$
@@ Line 59: / Line 91: @@
 $$
-where \(\gamma = C_P/C_V\) is the ratio of specific heats, \(\mu\) the mean molecular weight, and \(m_p\) the proton mass.
+where \(\gamma=C_P/C_V\) is the ratio of specific heats, \(\mu\) the mean molecular weight, \(m_p\) the proton mass, \(k\) Boltzmann’s constant, and \(g=GM(r)/r^2\).
-Convection dominates wherever the surrounding radiative gradient is steeper than this adiabatic value.
+Convection dominates wherever \(\left|\frac{dT}{dr}\right|_{\!rad} > \left|\frac{dT}{dr}\right|_{\!ad}\).
-(Full derivations and mixing-length formulations are covered in the article [[un:heat-transfer-in-stars|Heat transfer in stars]].)
+(We **avoid the full derivation** and mixing-length details here; see [[un:heat-transfer-in-stars|Heat transfer in stars]].)
-===== 6. Equation of state =====
+===== - Equation of state =====
-The gas pressure relates density and temperature through the **ideal-gas law**:
+The **equation of state (EOS)** connects macroscopic variables — pressure, density, and temperature — describing the physical state of the stellar gas.
+For a fully ionized ideal gas,
 $$
@@ Line 71: / Line 104: @@
 $$
-where \(m_{av}\) is the mean particle mass and \(\mu\) the mean molecular weight.
+where \(m_{av}\) is the mean particle mass and \(\mu\) the mean molecular weight (composition dependent).
-This equation links the thermodynamic state of the plasma with its composition.
+In real stars, radiation pressure, degeneracy pressure, or Coulomb corrections can be added, but the above relation is the baseline form for most stellar interiors.
-===== 7. Ionization balance (Saha equation) =====
+===== - Ionization balance (Saha equation) =====
 The ionization state of the stellar gas determines both opacity and mean molecular weight.
@@ Line 83: / Line 116: @@
  = \frac{2}{n_e}
    \left(\frac{2\pi m_e kT}{h^2}\right)^{3/2}
-   \frac{G_{i+1}}{G_i}
+   \frac{G_{i+1}}{G_i}\,
    e^{-E_i/(kT)},
 $$
-where \(n_i\) and \(n_{i+1}\) are number densities of successive ions, \(n_e\) is electron density, \(E_i\) the ionization energy,
+where \(n_i\) and \(n_{i+1}\) are number densities of successive ions, \(n_e\) is the electron density, \(E_i\) the ionization energy, and \(G_i\) the partition function of level \(i\).
-and \(G_i\) the partition function of level \(i\).
+This equation determines the degree of ionization and thus sets \(\kappa(\rho,T)\) and \(\mu(\rho,T)\) for the opacity and equation of state.
-This determines the degree of ionization and hence the opacity \(\kappa(\rho,T)\).
+===== Primary and secondary =====
+The **first five equations** — for hydrostatic balance, mass conservation, luminosity, and the two temperature gradients — are **primary** because:
+  * They are **differential equations** that describe the **structural variation** of a star with radius \(r\).
+  * They express the **fundamental conservation laws** of physics:
+    * (1) mechanical equilibrium,
+    * (2) conservation of mass,
+    * (3) conservation of energy, and
+    * (4–5) energy transport by radiation or convection.
+  * Their solutions give the radial profiles of \(P(r)\), \(\rho(r)\), \(T(r)\), \(M(r)\), and \(L(r)\).
+The **last two equations** — the equation of state and Saha ionization relation — are **secondary**, because:
+  * They are **constitutive relations** that specify how matter behaves under given physical conditions.
+  * They provide the **closure** needed to relate pressure, temperature, and density, allowing the primary equations to be solved.
+  * Without them, the system of five differential equations would contain more unknowns than equations and remain indeterminate.
+Together, the primary and secondary equations form a **closed system** of seven equations with seven unknown functions.
-===== Interdependence of the seven equations =====
+===== Interdependence and boundary conditions =====
-These seven equations describe how pressure, mass, luminosity, and temperature vary with radius in a spherically symmetric, non-rotating star.
+These seven equations describe how \(P(r),\,\rho(r),\,T(r),\,L(r)\), and \(M(r)\) vary with radius in a spherically symmetric, non-rotating star.
 They can be summarized as:
-| Equation | Physical principle | Variable(s) linked |
+^ Equation ^ Type ^ Physical principle ^ Variable(s) linked ^
-|-----------|--------------------|-------------------|
+| (1) \(dP/dr\) | Primary | Hydrostatic balance | \(P, \rho, M\) |
-| (1) \(dP/dr\) | Hydrostatic balance | \(P, \rho, M\) |
+| (2) \(dM/dr\) | Primary | Mass continuity | \(M, \rho\) |
-| (2) \(dM/dr\) | Mass continuity | \(M, \rho\) |
+| (3) \(dL/dr\) | Primary | Energy generation | \(L, \rho, T\) |
-| (3) \(dL/dr\) | Energy generation | \(L, \rho, T\) |
+| (4) \(dT/dr\) | Primary | Radiative transport | \(T, L, \rho, \kappa\) |
-| (4) \(dT/dr\) | Radiative transport | \(T, L, \rho, \kappa\) |
+| (5) \(dT/dr\) | Primary | Convective transport | \(T, g, \gamma, \mu\) |
-| (5) \(dT/dr\) | Convective transport | \(T, g, \gamma, \mu\) |
+| (6) \(P=\rho kT/\mu m_p\) | Secondary | Equation of state | \(P, \rho, T\) |
-| (6) \(P=\rho kT/\mu m_p\) | Equation of state | \(P, \rho, T\) |
+| (7) Saha equation | Secondary | Ionization equilibrium | \(\kappa, \mu, T, \rho\) |
-| (7) Saha equation | Ionization equilibrium | \(\kappa, \mu, T, \rho\) |
-Together, they form a closed system that defines a **stellar model** once the **boundary conditions** are specified:
+A stellar model is obtained by solving these equations with suitable **boundary conditions**, e.g.,
 $$
-M(0) = 0, \quad L(0) = 0, \quad
+M(0)=0,\quad L(0)=0,\quad P(R)=0,\quad T(R)=T_{eff},
-P(R) = 0, \quad T(R) = T_{eff}.
 $$
-By solving them numerically—using observed mass \(M\) and composition \(C\)—one can determine a star’s radius \(R\), luminosity \(L\), internal profiles of \(P(r)\), \(T(r)\), and \(\rho(r)\), and the boundaries between its **radiative** and **convective** zones.
+given the total mass \(M(R)\) and composition of the star.
 ===== Insights =====
-  - The seven equations together express mass, momentum, and energy conservation within a self-gravitating plasma.
+  - The **primary equations** describe the geometry and physical balances within a star; the **secondary equations** describe the state of the material composing it.
-  - Hydrostatic equilibrium ensures mechanical stability; the two energy-transport equations govern thermal equilibrium.
+  - Mass and luminosity equations follow directly from shell geometry: \(dV=4\pi r^2 dr\) and energy generation \(dL=\epsilon\rho dV\).
-  - The EOS and ionization balance connect microscopic particle physics to macroscopic structure.
+  - Hydrostatic equilibrium provides the mechanical backbone; EOS and Saha supply the thermodynamic closure.
-  - Radiative and convective transport determine where a star has radiative cores or convective envelopes.
+  - A steeper required flux (large \(L\)) or greater opacity \(\kappa\) makes the radiative \(dT/dr\) steeper; convection starts when this exceeds the adiabatic gradient.
-  - Solving these equations numerically yields the radius, luminosity, and internal structure for any given stellar mass and composition.
+  - Solving the seven coupled equations yields the full internal profiles of a star and defines its radiative and convective zones.
 ===== Inquiries =====
-  - Write down the physical meaning of each of the seven stellar structure equations.
+  - Which of the seven equations are differential and which are algebraic?
-  - What boundary conditions are required to solve them numerically?
+  - Why are the first five called primary and the last two secondary?
-  - Why is the radiative temperature gradient steeper in high-luminosity stars?
+  - Derive \(\frac{dM}{dr}=4\pi r^2\rho\) and explain its physical meaning.
-  - Under what condition does convection replace radiation as the dominant heat-transfer mechanism?
+  - Explain how the equation of state and Saha relation “close” the system of stellar structure equations.
-  - Explain how the EOS and Saha equation together determine the opacity profile of a star.
+  - Discuss how different opacity or composition affects the balance between radiative and convective energy transport.