Codex @codex  0

Interstellar shocks arise when relative flow is supersonic: stellar winds, supernova ejecta, colliding clouds and gas moving through strong disturbances can drive them. Neighboring compressive characteristics steepen or a fast flow encounters material it cannot decelerate through upstream sound communication. A thin dissipative front then converts directed kinetic energy into thermal energy and increases entropy.

Integrate mass, normal momentum and total energy conservation through a planar steady shock in its rest frame. With no mass source, negligible external impulse across the thin layer and no radiative energy loss, their normal fluxes are respectivelyFor a perfect gas of constant adiabatic index, $h=\gamma p/[(\gamma -1)\rho ]$. Equality of the fluxes on the two sides yields the three Rankine-Hugoniot conditions in the question. Adiabatic means no energy escapes this control volume; it does not make the shock isentropic.

In the strong-shock limit, let $r={\rho }_2/{\rho }_1=u_1/u_2$ and neglect $p_1$. Momentum conservation gives $p_2={\rho }_1{u}_1^2(1-1/r)$. Substitution into energy conservation givesThe nontrivial compressive root isThe $r=1$ root describes no shock. A monatomic gas has the familiar compression four.

For identical colliding clouds in their center-of-momentum frame, reflection symmetry fixes the contact plane and the uniform compressed material on either side at rest. The unshocked clouds still approach at speeds $\pm V$; their supersonic upstream flow has not received a smooth pressure adjustment ahead of the shock. Let the right shock move to the right at speed $S$. In its own frame the upstream and downstream normal speed magnitudes are $V+S$ and $S$. The same shock compression ratio therefore givesThis is the shock speed in a symmetric cloud collision. For $\gamma =5/3$, the shocks move outward at $V/3$, not at the upstream speed or the relative cloud speed $2V$. As a check, the post-shock pressure is $p_2=(\gamma +1){\rho }_1{V}^2/2$; in the monatomic case $p_2/{\rho }_2=V^2/3$, giving specific internal energy $V^2/2$, exactly the initial bulk kinetic energy per unit cloud mass when brought to rest.

When a shock reaches a cloud's rear boundary, its supply of cold upstream gas ends. If the surrounding pressure is low, the hot compressed gas expands outward and a rarefaction travels back into it, converting thermal energy into motion and cooling the gas. A sustained uniform stationary shocked region is no longer the solution. Interaction with an appreciable exterior medium can instead drive another outward disturbance; the boundary pressure determines that part of the evolution. The stated adiabatic model provides no radiative mechanism that would keep a permanently cold dense sheet at the contact.

Follow a fluid element with position $\mathbf{X}(t)$ satisfying $\dot{\mathbf{X}}=\mathbf{u}(\mathbf{X},t)$. The chain rule gives the material derivativeApply Newton's second law to a material element of volume $dV$ and mass $\rho dV$. The net pressure force is $-\nabla pdV$, and gravity gives $-\rho \nabla \varphi dV$. Dividing by the element mass, neglecting viscosity and using the chain rule componentwise gives the Euler equation

For a static spherical isothermal gas, $p=c_s^2\rho$ and $dp/dr=-GM(r)\rho /r^2$. Trial density $\rho =A/r^2$ gives enclosed mass $M(r)=4\pi Ar$, soThis is the singular isothermal sphere. Its central density diverges, although the enclosed mass tends to zero linearly and its total mass inside a finite radius is finite. It is not the finite-central-density regular isothermal solution.

Truncate at $r_0$ with boundary pressure $p_0$. Then $M=4\pi A{r}_0=2c_s^2{r}_0/G$ and $p_0=c_s^2\rho (r_0)$, givingThe external medium is being treated as a prescribed pressure boundary; its own gravitational field is excluded from the cloud self-energy calculation.

For a monatomic perfect gas, the equilibrium is isothermal but the internal-energy relation still uses its adiabatic index $\gamma =5/3$:The gravitational self-energy follows by assembling shells:Both integrals remain finite at the singular center. ConsequentlyThe pressure-confined singular isothermal gas sphere therefore satisfies the virial theorem including the work of the confining surface pressure. One can see its origin directly by integrating $\mathbf{r}\cdot \nabla p$: integration by parts gives $\oint p\mathbf{r}\cdot d\mathbf{S}-3\int pdV$, whose surface term is $4\pi r_0^3{p}_0$. An isolated vacuum-boundary version would omit that term and could not truncate this same positive-temperature isothermal profile at finite radius. The equality expresses equilibrium force and energy balance, rather than proving dynamical stability.

An Alfvén wave is a transverse disturbance restored by magnetic tension. In a uniform conducting fluid with background ${\mathbf{B}}_0=B_0{\mathbf{e}}_z$, a small transverse displacement ${\mathbf{\xi }}_{\perp }(z,t)$ bends the magnetic field lines. Magnetic flux freezing gives ${\mathbf{b}}_{\perp }=B_0{\partial }_z{\mathbf{\xi }}_{\perp }$, and the linear transverse Lorentz force density is $(B_0/{\mu }_0){\partial }_z{\mathbf{b}}_{\perp }$. Thus the Linear Alfvén wave equations reduce toFor this shear polarization, $\delta \rho =\delta p=0$ at linear order. The disturbance transports energy along the background magnetic field; a single traveling wave has equal perturbation magnetic energy and kinetic energy, with ${\mathbf{b}}_{\perp }=\mp \sqrt{{\mu }_0\rho }{\mathbf{u}}_{\perp }$ for propagation along $\pm {\mathbf{B}}_0$. These properties distinguish it from compressive magnetosonic waves.

The ideal magnetohydrodynamics approximation requires strong electrical coupling and wavelengths large compared with microscopic scales. Finite electrical resistivity and viscosity dissipate the wave's energy. In stellar winds and other astrophysical plasmas, Alfvén waves carry disturbances and can contribute to heating when dissipated. The Lundquist liquid-metal Alfvén wave experiment provides a laboratory counterpart: a conducting liquid in a strong applied magnetic field transmits a transverse magnetic disturbance, making the field's mechanical restoring action measurable.

A geometrically thin accretion disk around a central mass $M$ has nearly Keplerian orbits, with ${\Omega }^2\simeq GM/r^3$. Its vertical hydrostatic equilibrium gives $H\sim c_s/\Omega \ll r$. The specific angular momentum $\ell =r^2\Omega \propto r^{1/2}$ decreases inward, so accreting matter must transfer angular momentum to other matter. A viscous stress couples neighboring annuli rotating at different angular velocities.

Writing $\nu$ for effective kinematic viscosity and $\Sigma$ for surface mass density, the outward viscous torque in an accretion disk isSteady mass conservation defines $\dot{M}=-2\pi r\Sigma u_r>0$. If the inner radius $r_{\mathrm{in}}$ has negligible applied torque, angular momentum conservation gives $\mathcal{G}=\dot{M}(\ell -{\ell }_{\mathrm{in}})$, henceWell outside that boundary, $u_r\simeq -3\nu /(2r)$ and the inward drift time is of order $r^2/\nu$, much longer than an orbital period in a thin disk with modest effective stress. Dissipation of orbital energy heats the disk. If that heat escapes locally from two faces, the radiative flux from either face isIntegration over both faces gives $L=GM\dot{M}/(2r_{\mathrm{in}})$ in this Newtonian zero-torque model. Material can lose angular momentum while the disk's outer parts spread outward: outward transport does not imply that most of its mass moves outward. In real ionized disks, magnetorotational instability and magnetic stresses often provide the effective transport; molecular viscosity alone is generally insufficient.

In a radiative shock, a narrow dissipative jump first converts upstream kinetic energy into internal energy. The immediate strong adiabatic compression is $(\gamma +1)/(\gamma -1)$, and the immediate temperature is the result of part 1(c). Downstream gas then radiates, loses specific enthalpy and compresses further. Its shock cooling length is of order $u_2{t}_{\mathrm{cool}}$, with the immediate downstream speed used only as an estimate because the gas slows as it cools.

For the idealization of equal far-upstream and cooled downstream temperatures, unchanged composition and negligible magnetic stress, write $p=c_T^2\rho$. Conserved mass and normal momentum flux across the whole layer giveWith $r={\rho }_f/{\rho }_1=u_1/u_f$, the nontrivial isothermal shock root isIt can greatly exceed the adiabatic strong-shock compression. An upstream adiabatic sound speed instead defines ${\mathcal{M}}_{\mathrm{ad},1}=u_1/(\sqrt{\gamma }{c}_T)$, so the same result is $r=\gamma {\mathcal{M}}_{\mathrm{ad},1}^2$; the two Mach numbers must not be confused. The radiated energy flux is ${\rho }_1{u}_1(u_1^2-u_f^2)/2$ when the two endpoint specific enthalpies agree. The gas's specific entropy rises at the initial jump and falls during cooling, while the escaping radiation carries entropy away. A transverse frozen-in magnetic field, or an imposed temperature floor, can limit the compression. Strongly cooling shocks can therefore collect thin dense shells without violating total conservation of energy.

A stellar wind carries matter away from a star, altering its mass, angular momentum and surroundings. Hot coronal gas can escape through pressure-driven expansion; radiation transfers momentum to gas through continuum or spectral line absorption, or to dust that then drags gas. A magnetic field can transfer rotational energy and angular momentum to an outflow. These mechanisms need not give the same terminal speed or mass loss rate.

For a steady spherical pressure-driven model, conservation of mass gives $\dot{M}=4\pi r^2\rho u$. Combining the radial Euler equation with an isothermal equation of state $p=\rho a^2$ gives the Parker wind equationA smooth transonic branch passes through $u=a$ at $r_s=GM/(2a^2)$, connecting a subsonic inner atmosphere to an escaping supersonic wind. The isothermal Parker wind requires continuing heat input as the gas expands; a polytropic stellar wind instead tracks a changing temperature and must satisfy its own terminal-energy condition. A subsonic isothermal breeze generally tends to a nonzero external pressure, so outward velocity alone does not select the escaping solution.

The solar wind supplies a nearby example of an outward plasma flow carrying a magnetic field. Profiles of spectral lines, especially absorption shifted by the Doppler effect together with emission, diagnose more distant stellar winds. The critical-point condition and boundary heating constrain $\dot{M}$; it is not determined by an escape-speed estimate alone.

An accretion disk forms when inward-moving gas retains enough angular momentum to orbit a compact object. In a geometrically thin Keplerian disk, radial force balance gives $\Omega \simeq \sqrt{GM/r^3}$ and vertical hydrostatic equilibrium gives a thickness $H\sim c_s/\Omega \ll r$. Matter cannot fall far inward while keeping its original specific angular momentum $r^2\Omega \propto r^{1/2}$: stresses must transport angular momentum outward.

The viscous torque in an accretion disk has outward magnitudewhere $\nu$ is the effective kinematic viscosity and $\Sigma$ the surface density. With constant inward accretion rate $\dot{M}>0$ and zero inner torque at $r_{\mathrm{in}}$, conservation of angular momentum givesThe associated standard thin-disk dissipation flux from each face isOrbital gravitational potential energy becomes heat and, in a radiatively efficient disk, escaping radiation. Integrating both faces gives the standard thin-disk luminosity $L=GM\dot{M}/(2r_{\mathrm{in}})$ under this boundary condition. A central stellar surface or an inner torque can release or redistribute additional energy.

The alpha disk prescription $\nu \sim \alpha c_{s}H$ summarizes an effective stress rather than a claim that molecular viscosity alone explains observed rates. In an ionized disk, magnetorotational instability can produce magnetic and fluid stresses. The inward evolution occurs on the viscous timescale $t_{\nu }\sim r^2/\nu$, much longer than an orbit for a thin disk with $\alpha \lesssim 1$. Hot inner annuli and cooler outer annuli naturally give a range of emitted temperatures; the blackbody approximation additionally requires sufficient optical depth.

The Mach number is $\mathcal{M}=u/c_s$, comparing the flow speed relative to a disturbance with the local sound speed. For a body moving supersonically, spherical acoustic disturbances emitted earlier have a tangent envelope with Mach cone half-angle $\sin \theta =1/\mathcal{M}$. Signals cannot propagate arbitrarily far upstream. Compressive characteristics can converge and steepen into a shock wave, within which microscopic dissipation prevents a multivalued fluid state; supersonic motion alone does not require every flow to contain a shock.

Work in the shock frame with steady normal velocities $u_1,u_2>0$. Integrating conservation of mass, normal momentum and total energy over a thin pillbox, neglecting gravity across the thin layer, gives the Rankine-Hugoniot conditions for a perfect gas:The energy equation includes pressure work through specific enthalpy. The shock is adiabatic if negligible heat escapes the layer, but its irreversible dissipation means it is not isentropic.

Put $\chi ={\rho }_2/{\rho }_1=u_1/u_2$ and ${\mathcal{M}}_1=u_1/\sqrt{\gamma p_1/{\rho }_1}$. Momentum conservation gives $p_2=p_1+{\rho }_1{u}_1^2(1-1/\chi )$. Substitution in the energy equation produces the unshocked root $\chi =1$ and the compressive rootSince ${\rho }_1{u}_1^2/p_1=\gamma {\mathcal{M}}_1^2$, the strong-shock limit has negligible upstream thermal pressure compared with ram pressure. The Strong-shock Rankine-Hugoniot conditions giveFor unchanged composition, the ideal gas law implies $T_2/T_1=(p_2/p_1)/\chi$, soFor $\gamma =5/3$, the compression tends to four and $T_2\simeq 3\bar{m}{u}_1^2/(16k_B)$, where $\bar{m}$ is mean particle mass.

With constant specific heat $c_v$, entropy production in a perfect-gas shock is measured by the change in specific entropyThe strong-shock expression is thereforeThe exact jump is positive for every compressive shock with ${\mathcal{M}}_1>1$, in accordance with the Second law of thermodynamics; the asymptotic expression is only for ${\mathcal{M}}_1\gg 1$. There is a normalization issue in the printed description: this logarithmic formula with constant $c_v$ describes entropy per unit mass. The true entropy density is $\rho s$, so its difference is ${\rho }_2{s}_2-{\rho }_1{s}_1$, rather than the same logarithmic jump with a constant heat capacity per unit volume. A material parcel's irreversible entropy increase should be assessed using $\Delta s$.

Downstream gas is denser and hotter. For optically thin collisional emission the loss per unit volume is typically $n_e{n}_i\Lambda (T)$, whereas stored thermal energy scales approximately as $nT$. Compression can therefore reduce the radiative gas cooling time; excitation and ionization at the higher temperature can open additional cooling channels. The actual efficiency still depends on the astrophysical cooling function, composition and optical depth. In a radiative shock with separated scales, the temperature first rises across a narrow dissipative layer, then decreases through a much broader cooling region of characteristic length ${\ell }_{\mathrm{cool}}\sim u_2{t}_{\mathrm{cool}}$. Subsonic cooling usually brings further compression with approximately constant pressure. The immediate microscopic shock thickness and the macroscopic cooling length are distinct.

Let state $f$ be the fully cooled gas, with $T_f=T_1$ and the same composition. Set $a^2=p_1/{\rho }_1$, the isothermal sound speed squared. Across the entire shock plus cooling layer, mass and momentum are still conserved, but gas energy alone is not. With ${\chi }_f={\rho }_f/{\rho }_1$,The nontrivial isothermal shock therefore hasIf the upstream Mach number were defined with $a$ instead, this would be ${\rho }_f/{\rho }_1={\mathcal{M}}_{\mathrm{iso}}^2$. Retaining the distinction avoids losing the factor $\gamma$.

The gas's specific entropy rises irreversibly at the hot shock and then decreases as escaping radiation removes heat. Since the final temperature equals the initial temperature,This decrease for the gas is consistent with the Second law of thermodynamics: the emitted radiation and its eventual surroundings carry entropy away. It does not make the microscopic compressive shock reversible. The figure illustrates both the hot entropy jump and the subsequent cooling; its spatial cooling profile is schematic because no cooling law has been specified.

For a steady, spherically symmetric radial inflow, conservation of mass makes the mass accretion rate independent of radius:Atomic hydrogen gives $\rho \simeq m_{H}n$. The Newtonian free-fall speed from rest at infinity is $|v_r|=\sqrt{2GM/r}$, obtained by setting the specific mechanical energy $v_r^2/2-GM/r$ to zero. Hence the spherical free-fall accretion rate isUsing the given values gives $\rho \simeq 1.67\times 10^{-21}\mathrm{kg}{\mathrm{m}}^{-3}$, $|v_r|\simeq 2.98\times 10^5\mathrm{m}{\mathrm{s}}^{-1}$, andThe area factor assumes a full spherical inflow with negligible angular momentum. With only a local number density and speed, an arbitrary nonspherical inflow's total rate is not fixed; one instead integrates $-\rho \mathbf{v}\cdot d\mathbf{S}$ over the inflow surface. This is a flux estimate, rather than the pressure-selected Bondi accretion solution.

A uniform conducting fluid supports three propagating magnetohydrodynamic wave branches. Ideal magnetohydrodynamics couples fluid velocity to the magnetic field through the induction equation; the Lorentz force then supplies magnetic tension and magnetic pressure restoring forces. Let $\theta$ be the angle between wave vector $\mathbf{k}$ and background field ${\mathbf{B}}_0$, and define Alfvén speed $v_A=B_0/\sqrt{{\mu }_0{\rho }_0}$ and adiabatic sound speed $c_s^2=\gamma p_0/{\rho }_0$.

The shear Alfvén wave hasIts velocity perturbation is perpendicular to both $\mathbf{k}$ and ${\mathbf{B}}_0$. It bends field lines without first-order compression, so magnetic tension is the restoring force. Its group velocity follows the background field; the displayed phase speed along $\mathbf{k}$ is not the full group speed.

The fast magnetosonic wave and slow magnetosonic wave involve coupled compression and magnetic perturbations, with displacement in the plane of $\mathbf{k},{\mathbf{B}}_0$:The plus sign is the fast branch. Parallel propagation gives characteristic speeds $c_s,v_A$; perpendicular propagation gives fast speed $\sqrt{c_s^2+v_A^2}$ while the other two phase frequencies vanish. Degenerate directions require care in assigning polarizations, rather than claiming three distinct nonzero frequencies everywhere.

These waves transport energy and momentum through magnetized stellar atmospheres and astrophysical plasmas. Their restoring mechanisms distinguish them from ordinary acoustic waves. Resistivity, viscosity, ion-neutral friction and kinetic effects introduce damping or change the model; an ideal undamped dispersion relation alone does not calculate coronal heating.

For a Newtonian circular orbit, $v_{\varphi }^2=GM/r$, so the specific orbital energy is $e=v_{\varphi }^2/2-GM/r=-GM/(2r)$. Half of the released potential energy remains in circular orbital motion. The binding-energy estimate between the specified edges is thereforeUsing $M_{\odot }=1.98847\times 10^{30}\mathrm{kg}$ and one Julian year $=3.15576\times 10^7\mathrm{s}$ gives $\dot{m}=6.30\times 10^{22}\mathrm{kg}{\mathrm{s}}^{-1}$. Thus the radiative luminosity scale is $L\sim 1.9\times 10^{38}\mathrm{W}$, approximately $1.9\times 10^{45}\mathrm{erg}{\mathrm{s}}^{-1}$ or $4.9\times 10^{11}{L}_{\odot }$. The black-hole mass cancels because both edge radii scale with $M$. The specified inner edge is $6R_s$; it must not be silently replaced with a different relativistic orbit radius.

For a finite viscous disc, the precise radiated power also depends on boundary work: torque transport of energy in a steady accretion disk means that orbital binding-energy release need not all be radiated inside the two stated edges. If the intended model is a radiatively efficient steady Keplerian accretion disk with a zero-torque inner boundary condition, integrate the one-face flux $F=3GM\dot{m}[1-\sqrt{r_{\mathrm{in}}/r}]/(8\pi r^3)$ over both faces. It givesThe difference from $L_{\mathrm{bind}}$ is the outward torque power at the outer edge. The question supplies neither boundary torque nor the external material-supply energy budget, so the two calculations state their assumptions explicitly: about $2\times 10^{38}$ W is the binding-energy estimate, while $1.37\times 10^{38}$ W is the finite-disc zero-inner-torque dissipation result. Both neglect radiation advection, winds and relativistic corrections as required for this model.^[ref]

Immediately upstream, free fall from rest at infinity gives $u_1^2=2G{M}_0/r_0$. Treat the surface shock as locally planar and stationary, with negligible upstream pressure. The Strong-shock Rankine-Hugoniot conditions for a monatomic gas, $\gamma =5/3$, give ${\rho }_2/{\rho }_1=4$ and $u_2=u_1/4$. momentum conservation then givesUse the source's gas-constant convention $p=\rho (R^{\ast }/\mu )T$. The immediate postshock temperature isHere adiabatic means that the jump does not lose appreciable heat to the surroundings; the shock is nevertheless irreversible and raises entropy. The gas still has velocity $u_1/4$ immediately downstream, so setting all upstream kinetic energy equal to downstream internal energy would give the wrong coefficient. Later settling, radiative cooling or changing ionization can alter the temperature; they are outside this immediate, fixed-composition perfect-gas jump.

A finite-amplitude compression can steepen because different parts of the wave travel at different characteristic speeds. For a right-moving simple wave in a perfect gas, one may write $dx/dt=c_0+(\gamma +1)u/2$ in terms of the local velocity perturbation $u$. Faster compressive characteristics catch slower ones. When the smooth inviscid description would become multivalued, a physical compression shock replaces it. A rarefaction normally spreads rather than forming an entropy-violating expansion shock.^[ref]

Across a locally planar shock, use velocities relative to its moving front. The Rankine-Hugoniot conditions for a perfect gas conserve mass, momentum and total energy flux:where $h=\gamma p/[(\gamma -1)\rho ]$ is the specific enthalpy. For upstream normal Mach number $M_1>1$,Thus a normal compression shock increases pressure, density and temperature, reduces the normal velocity, and changes a supersonic upstream flow to subsonic downstream flow in the shock frame. In the strong-shock limit, compression approaches $(\gamma +1)/(\gamma -1)$, which is four for a monatomic gas. A shock conserves total energy but generates entropy; an adiabatic shock is not an isentropic compression. For a calorically perfect gas, stagnation temperature is conserved while stagnation pressure falls.^[ref]

Real dissipative processes resolve the narrow front; the jump conditions do not require its detailed internal profile. In a continuum-gas experiment a shock tube can produce a moving shock by releasing a pressure discontinuity. pressure probes separated along the tube measure arrival speed and the pressure jump; distinguishing the shock from the following contact surface is essential. A contact surface can have a density/entropy jump with continuous pressure and normal velocity, whereas a shock carries nonzero mass flux through its front.

A small pressure disturbance propagates because compression raises the pressure, the pressure gradient accelerates nearby fluid, and the induced convergence compresses the next region. For a uniform fluid at rest, linearizing mass and momentum conservation and using the adiabatic pressure-density response givesEliminating ${\rho }^{\prime }$ and ${\mathbf{u}}^{\prime }$ gives ${\partial }_t^2{p}^{\prime }=c_s^2{\nabla }^2{p}^{\prime }$. A plane wave has ${\omega }^2=c_s^2|\mathbf{k}{|}^2$ and longitudinal velocity parallel to $\mathbf{k}$. In a perfect gas,The adiabatic response is appropriate when heat cannot diffuse appreciably out of each compressed region during a period; an isothermal response would give a different derivative. In a moving uniform medium, the acoustic characteristics travel at $U\pm c_s$ along the flow, rather than at $\pm c_s$ in the laboratory.

For a progressive wave, $p^{\prime }={\rho }_0{c}_s{u}_{\parallel }^{\prime }$; pressure, density and longitudinal velocity are in phase. The acoustic energy density is the sum of kinetic and compression energies,Sound transports energy while fluid elements execute small oscillations about their mean positions. The ideal uniform-fluid acoustic mode is nondispersive. Spherical spreading reduces pressure amplitude as $1/r$ when total power is conserved; viscosity and heat conduction attenuate waves, while sufficiently large compressions can steepen into shocks. In a stratified astrophysical gas the uniform result is local: wavelengths short compared with background scale heights justify it, whereas long waves can couple to buoyancy and acoustic-cutoff behavior.

Use the idealized fully developed, laminar cylindrical-flow model implied by the fixed cross-section and zero velocity at its sides. Take the jet length to be $R$, as for a central source reaching the cloud boundary; if the source is off-centre, the actual path length must replace $R$. Let $s$ denote distance from the jet axis and write the axial velocity as $u(s)$. Then $(\mathbf{u}\cdot \nabla )\mathbf{u}=0$, the pressure gradient is $dp/dz=-p_0/R$, and the axial Navier-Stokes equation reduces toRegularity at $s=0$ excludes a logarithmic term. The imposed no-slip boundary condition $u(r_0)=0$ gives the Hagen-Poiseuille flow profileThe kinetic energy crossing each area element per unit time is $(\rho u^2/2)udA$, so the velocity must be cubed before integration. ThusConsequentlyThe mean velocity is $\bar{u}=U_{\max }/2$. Replacing the profile by $\bar{u}$ in $\frac{1}{2}\rho \pi r_0^2{\bar{u}}^3$ underestimates the flux by two; the kinetic energy correction factor is two for this profile. The radius dependence follows from $u\propto r_0^2$ and cross-sectional area $\propto r_0^2$.

The specified $P$ is outgoing bulk kinetic-energy flux, not the pressure-work rate $p_{0}Q$, which scales as $r_0^4$. In a fully developed section the profile is independent of $z$ and pressure work is dissipated as heat; entrance acceleration and energy supply establish the profile. A turbulent or freely entraining astrophysical jet need not have this no-slip parabolic structure. The model assumes a long enough, laminar confined region and neglects entrance losses and gravity.

An accretion disk forms when infalling material has enough angular momentum to orbit before it can reach the central object. In a thin disc dominated by a central mass $M$, the orbital motion is approximately Keplerian rotation:Vertical hydrostatic equilibrium gives the disk scale height $H\sim c_s/{\Omega }_K$; a geometrically thin disc has $H/R\ll 1$. The surface density is $\Sigma =\int \rho dz$, and a steady inward mass flow rate is $\dot{M}=-2\pi R\Sigma v_R>0$.

Moving inward requires losing specific angular momentum. Internal stresses transport angular momentum outward; much of the mass can then accrete while an outer fraction spreads to larger radius. Molecular viscosity alone is generally inadequate for rapid astrophysical accretion. Turbulent and magnetic stresses, including those generated by the magnetorotational instability, provide effective transport. A common parametrization is $\nu =\alpha c_{s}H$ for an effective kinematic viscosity, with $\alpha$ a dimensionless stress parameter.

For a thin Keplerian accretion disk, the evolution equation isThis viscous evolution of an accretion disk occurs on the viscous timescale $t_{\nu }\sim R^2/\nu$, much longer than an orbital period in a thin disc with moderate effective stress.

The orbital specific energy is $-GM/(2R)$, so inward motion releases binding energy. In a steady, optically thick, radiatively efficient disc with zero torque at its inner edge, the emergent flux from each face isThe Stefan–Boltzmann law identifies the local effective temperature. Far from the inner boundary, $T_{\mathrm{eff}}\propto R^{-3/4}$, so the emitted spectrum combines hotter inner and cooler outer annuli. Integrating both faces gives $L_{\mathrm{disc}}=GM\dot{M}/(2R_{\mathrm{in}})$ under these assumptions. This is not a universal spectrum or efficiency: optical depth, inner torque, energy advection and outflows can change the radiation. The essential connection is outward angular-momentum transport, inward accretion, and release of binding energy.

Take positive $\dot{Q}$ to mean net removal of thermal energy per unit mass. The isobaric thermal stability of a heated ideal gas depends on the density change accompanying a temperature disturbance. For an ideal gas of fixed composition, $p=\rho k_{B}T/\bar{m}$, so constant pressure gives $\rho T=\text{constant}$. Therefore the net cooling can be writtenwhere $C>0$ is constant for this pressure. At thermal equilibrium $T=T_0$, $C{T}_0^{\alpha -1}=H$.

The first law of thermodynamics at constant pressure includes the expansion work, so the relevant heat capacity is the specific heat capacity at constant pressure, $c_p$, rather than the constant-volume value. Linearizing its temperature equation givesHenceFor optically thin Bremsstrahlung, $\alpha =1/2$, soA hotter parcel expands and becomes less dense. Its cooling per unit mass, proportional to $T^{-1/2}$ at fixed pressure, decreases while its cosmic ray heating stays fixed, making it hotter still. A colder, denser parcel cools faster and runs away in the opposite direction. Here $H$ has units $\mathrm{W}\mathrm{k}{\mathrm{g}}^{-1}$ and $c_p{T}_0$ units $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$, so $s$ is an inverse time. The conclusion assumes the disturbance remains close to pressure equilibrium and neglects heat conduction and changes of composition. The relevant $\alpha$ is the temperature exponent in the supplied cooling law.

Write the initial mass density as $\rho =mn$. The ideal gas law gives pressure $p=n{k}_{B}T$, and the sphere has $M=4\pi \rho R^3/3$. Here $k_B$ is the Boltzmann constant and $G$ is the gravitational constant. First obtain the gravitational potential energy by adding concentric shells. The mass inside radius $r$ is $M(r)=M{r}^3/R^3$, soIn a small contraction $dR<0$ at fixed total mass, the gravitational energy released is $-dW=3G{M}^2|dR|/(5R^2)$. The positive work required against the gas pressure is $-pdV=4\pi R^{2}p|dR|$. Equating these givesUse $R=(3M/(4\pi mn){)}^{1/3}$ to eliminate the radius. The uniform-sphere pressure-work collapse criterion isThe units are kilograms: $(k_{B}T/(Gm){)}^{3/2}$ has units $\mathrm{k}{\mathrm{g}}^{3/2}{\mathrm{m}}^{-3/2}$, and ${\rho }^{-1/2}$ supplies $\mathrm{k}{\mathrm{g}}^{-1/2}{\mathrm{m}}^{3/2}$. At fixed initial $n$ and $T$, $M>M_{\mathrm{crit}}$ makes gravity stronger than the pressure-work requirement in this model, allowing contraction.

This is a Jeans mass estimate for homologous contraction of an initially uniform sphere. It is not an exact stability limit for a hydrostatic gas sphere: uniform mass density and uniform temperature would give uniform pressure, rather than the gradient required by hydrostatic equilibrium. The numerical coefficient therefore differs from conventions based on a Jeans instability wavelength or a pressure-confined isothermal sphere. Subsequent contraction also depends on heat exchange; isothermal pressure support during gravitational collapse weakens relative to gravity as the radius decreases.

The local conservation laws. The first equation is the radial Euler equation: the material acceleration $Du/Dt={\partial }_{t}u+u{\partial }_{r}u$ is caused by the pressure gradient force per unit mass, $-{\partial }_{r}p/\rho$. It neglects gravity and viscosity. The second equation is mass conservation in spherical coordinates,The term $2\rho u/r$ appears because the area of a spherical surface grows as $r^2$.

To obtain the third equation, let $e=p/[(\gamma -1)\rho ]$ be the specific internal energy of an ideal gas with constant heat capacity ratio $\gamma$. The first law of thermodynamics for a smooth adiabatically moving element isDifferentiating $e$ and collecting terms givesand henceThus the entropy parameter $p/{\rho }^{\gamma }$ is conserved along each smooth pathline; for fixed composition, specific entropy differs from $c_v\ln (p/{\rho }^{\gamma })$ only by a constant. This is local isentropic flow behind the shock. A physical shock wave itself generates entropy and must be treated with jump conditions, rather than by applying this smooth equation across the discontinuity.

The front-radius scaling. In the self-similar blast wave approximation there is no scale set by ambient pressure, cooling, gravity, or the initial source radius. Write $R\propto E^a{\rho }_0^b{t}^d$. Using $[E]=\mathsf{M}{\mathsf{L}}^2{\mathsf{T}}^{-2}$ and $[{\rho }_0]=\mathsf{M}{\mathsf{L}}^{-3}$, dimensional analysis givesThus $a=1/5$, $b=-1/5$, $d=2/5$ andThe dimensionless coefficient ${\xi }_0$ depends on $\gamma$ and the resolved similarity profiles; dimensions alone do not fix it. Define the length scale $\ell (t)=(E{t}^2/{\rho }_0{)}^{1/5}$, so $\xi =r/\ell$ and the front is at the fixed coordinate $\xi ={\xi }_0$.