Codex @codex  0

The energy-conserving self-similar blast wave requires a stage in which the original stellar radius, ejecta mass, ambient pressure and cooling have ceased to supply important scales. In particular, the swept-up mass must greatly exceed the ejecta mass; $R\gg R_{\ast }$ alone need not ensure this in a very dilute interstellar medium. For fixed explosion energy $E$, uniform ambient mass density ${\rho }_0$ and elapsed time $t$, dimensional analysis gives $R=C{E}^a{\rho }_0^b{t}^d$. Matching mass, length and time gives $a=1/5$, $b=-1/5$, $d=2/5$. Equivalently, a self-similar strong shock wave has $E\sim {\rho }_0{R}^3{\dot{R}}^2$, whose constant-energy solution has the same exponents:The dimensionless coefficient depends on $\gamma$ and the post-shock profiles; dimensional reasoning cannot fix it.

For completeness, obtain the Strong-shock Rankine-Hugoniot conditions in the shock frame. Let $v_1=V-u$ denote downstream speed relative to the outward-moving front. With negligible upstream pressure and enthalpy, conservation of mass, momentum and energy givesPut $q={\rho }_1/{\rho }_0$. Substituting $v_1=V/q$ into the last two equations and selecting the compressive, nontrivial branch yields $q=(\gamma +1)/(\gamma -1)$. HenceHere the symbol $u$, written as $u=2V/(\gamma +1)$ in what follows, is the laboratory-frame radial velocity; it is not the relative speed $v_1$.

Use the thin-shell approximation for a spherical blast wave. Set $W=4\pi R^3/3$ and $M_s={\rho }_{0}W$. Equating the swept-up mass to $4\pi R^{2}d{\rho }_1$ gives $d/R=(\gamma -1)/[3(\gamma +1)]$. The growing shell's momentum is $M_{s}u$. Since newly swept gas has zero laboratory momentum, the net force is the cavity pressure acting on the inner shell, to thin-shell accuracy:The pressure at the shock is already accounted for by the jump conditions and mass loading; subtracting $p_1$ as an additional external pressure would count the shock reaction twice. With $u=2V/(\gamma +1)$, the last equation becomesTherefore, for constant $\alpha$,

The shell volume is $M_s/{\rho }_1=W/q$. Its internal energy and kinetic energy areTreat the cavity volume as $W$ at the same thin-shell level, so its internal energy is $E_{c,\mathrm{int}}=\alpha p_{1}W/(\gamma -1)$. Their sum is proportional to $R^3{V}^2$ for constant $\alpha$ and $\gamma$. Constant total energy thus requires $3+6(\alpha -1)=0$, giving $\alpha =1/2$. ConsequentlyUsing $V=2R/(5t)$, this becomesFor $\gamma =5/3$, $d/R=1/12$ and the approximate energy fractions are $E_{s,\mathrm{kin}}/E=1/4$, $E_{s,\mathrm{int}}/E=1/4$, $E_{c,\mathrm{int}}/E=1/2$.

This is an approximate normalization, not an exact treatment of finite shell thickness. If instead one subtracts the shell volume and uses $W_c=W(1-1/q)$ for the cavity, while retaining the other shell assumptions, the energy coefficient becomes $2(2\gamma -1)/[(\gamma -1)(\gamma +1{)}^2]$ and ${\xi }_0^5=75(\gamma -1)(\gamma +1{)}^2/[32\pi (2\gamma -1)]$. The difference is an order-$d/R$ geometric correction of this crude shell model; a consistent exact coefficient requires the full self-similar profiles.

In the idealized mass budget all swept gas is assigned to the shell, so the cavity inertia is neglected. The hot cavity's residual mass density is not specified, so its kinetic energy cannot be assigned a unique numerical fraction from the stated data. If its low mass is $M_c$ and its velocity is of order $u$, then $E_{c,\mathrm{kin}}\sim M_c{u}^2/2$ andFor an illustrative uniform-density cavity with homologous velocity $v(r)=ur/R$, direct radial integration gives $E_{c,\mathrm{kin}}=3M_c{u}^2/10$, so the preceding fraction has an extra factor $3/5$. The cavity can carry substantial internal energy despite negligible kinetic energy because it is hot. Neglect of its inertia is justified in the swept-mass-dominated stage, not merely by its small initial radius.

The calculation describes the adiabatic expansion stage of a supernova remnant, giving an estimate of its age or explosion energy from its radius, expansion speed and ambient mass density. For example $t\simeq 2R/(5V)$, and $E\propto {\rho }_0{R}^3{V}^2$. Such blasts inject energy and momentum into the interstellar medium and redistribute stellar ejecta. The early ejecta-dominated stage, radiatively cooled late shell, nonuniform ambient medium, finite ambient sound speed, magnetic fields, nonspherical geometry and energy transfer to cosmic rays violate different assumptions. Even in the adiabatic spherical stage, the real profiles are not a uniform thin shell and a uniform-pressure cavity, so ${\xi }_0$ here is an estimate.

Direct evidence for stellar feedback includes expanding ionized shells and superbubbles around young associations, hot X-ray-emitting gas in supernova remnants, broad or split emission lines, and blueshifted absorption showing cool and warm outflows. P-Cygni profiles reveal massive-star winds, while extraplanar filaments and metal-enriched gas demonstrate transport away from star-forming disks. Indirect evidence includes the galaxy mass--metallicity relation, low baryon fractions and suppressed star formation in dwarf galaxies, chemically enriched circumgalactic gas, and correlations of outflow speed and mass loading with the star formation rate.

Rapid gas removal changes the gravitational potential before stellar and dark-matter orbits can respond adiabatically. Positions and velocities are initially unchanged, but orbital binding energies rise; orbits expand and become more eccentric, and some particles escape. Repeated burst--outflow--reaccretion cycles can irreversibly transfer energy to collisionless matter and turn a central dark-matter cusp into a core. This matters because dwarf-galaxy rotation curves are used to test dark-matter microphysics: a feedback-made core can mimic a non-cold or self-interacting dark-matter signature.

Write the initial potential energy as $W=-aG{M}^2/R$. The virial theorem gives $T=-W/2$. If a well-mixed fraction $ϵ$ remains after instantaneous mass loss, the immediate kinetic and potential energies are $T_a=ϵT$ and $W_a={ϵ}^{2}W$, soAfter revirialization at $R^{\prime }$, $E_f=W_f/2=-aG{ϵ}^2{M}^2/(2R^{\prime })$. Equating energies gives the impulsive mass-loss expansion lawThe remnant is bound only for $ϵ>1/2$; loss of half or more of the gravitating mass disrupts this idealized system.

For many infinitesimal, individually revirialized losses, put $ϵ=1+dM/M$ in the impulsive result. To first order, $dR/R=-dM/M$. Integration yields the adiabatic mass-loss expansion lawSlow loss causes finite expansion for every positive remaining mass fraction and has no sharp disruption threshold.

Finally consider an initially circular orbit of radius $R$ around a point mass $M$. Its speed and specific angular momentum obey $v^2=GM/R$ and $h^2=GMR$. Immediately after $M\to ϵM$, these remain unchanged, while the new specific energy isUsing the orbital eccentricity relation $e^2=1+2E^{\prime }{h}^2/(G^2{ϵ}^2{M}^2)$ givesIt is an ellipse for $ϵ>1/2$, parabolic at $ϵ=1/2$, and unbound for smaller $ϵ$, in agreement with the virial argument.