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.

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.

The Rayleigh-Taylor instability extracts gravitational potential energy when a heavy fluid is supported above a lighter one. A small downward excursion of heavy material and upward excursion of light material can grow instead of experiencing restoring buoyancy. The Kelvin-Helmholtz instability instead extracts kinetic energy from a velocity difference across an interface or shear layer. The interface displacement and its induced flow reinforce one another by tapping the background shear; a reversed density ordering is not required.

A useful common model has two inviscid, incompressible, semi-infinite fluids separated at $z=0$, with lower density/velocity ${\rho }_1,U_1$, upper density/velocity ${\rho }_2,U_2$, downward gravity $g$, and surface tension $\sigma$. For a disturbance proportional to $e^{i(kx-\omega t)}$, with $k>0$, the two-fluid interfacial gravity-shear dispersion relation isA negative right side gives exponential growth. With no shear or surface tension and ${\rho }_2>{\rho }_1$, the growth rate of the Rayleigh-Taylor instability is $s=\sqrt{\mathcal{A}gk}$, where $\mathcal{A}=({\rho }_2-{\rho }_1)/({\rho }_2+{\rho }_1)$ is the Atwood number. With equal densities and no restoring force, the growth rate of the Kelvin-Helmholtz instability is $s=k|U_1-U_2|/2$. Heavy-over-light stratification drives Rayleigh-Taylor growth; velocity shear drives Kelvin-Helmholtz growth. A heavy lower fluid supplies gravitational stabilization, and surface tension stabilizes sufficiently short waves in this ideal interface model.

Finite shear-layer thickness, viscosity, compressibility and magnetic tension can change the growth rates and unstable wavelength range. In particular, unlimited growth with increasing $k$ is a limitation of the zero-thickness, zero-surface-tension idealization. Density-interface overturning and jet-boundary mixing are appropriate physical applications of the two mechanisms, but the nonlinear mixed state is not determined by the linear growth rate alone.

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.

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]

Let $z$ increase upwards and let $g>0$ denote the downward gravitational acceleration. The environment obeys hydrostatic equilibrium, $dp/dz=-\rho g$. Displace a small parcel by $\zeta$ while it remains in pressure balance with its surroundings and exchanges negligible heat. For fixed composition and adiabatic index $\gamma$, its density change at constant specific entropy satisfies $d{\rho }_{\mathrm{par}}/dp=\rho /(\gamma p)$. At the displaced position, its density excess over the environment is thereforeThe buoyancy acceleration is $\ddot{\zeta }\simeq -g\Delta \rho /\rho$. HenceStability requires $N^2>0$, equivalentlyAn upward-displaced parcel is then denser than its surroundings and falls back; a downward-displaced parcel is lighter and rises. Equality gives neutral stability in this parcel model. For a perfect gas with constant composition, $ds=c_{v}d\ln (p/{\rho }^{\gamma })$ gives $N^2=(g/c_p)ds/dz$, so a rising specific entropy is stable. This is the parcel form of the Schwarzschild criterion. Rapid pressure equilibration, negligible heat diffusion during displacement and no composition-driven buoyancy are the assumptions behind this test.

Insert the pressure-negligible velocity law into the steady spherical mass flux. With $x=r/R_e$,An adiabatic equation of state $p{\rho }^{-\gamma }=\text{constant}$ and the ideal gas law giveFor a monatomic gas with $\gamma =5/3$,The normalization $T_e$, and hence $p_e=n_e{k}_B{T}_e$, was not specified and remains a boundary parameter. If it is large enough for the pressure gradient to affect acceleration, replace the ballistic velocity by the full adiabatic Bernoulli equation solution. The general mass-flux and thermodynamic relations still hold, but these particular radius powers then need not hold.

For a steady spherical wind, conservation of mass givesThusHere $n_e$ denotes the given number density at the escape radius, rather than electron density; the composition is atomic hydrogen. The reference speed is $v_e\simeq 4.22\times 10^4\mathrm{m}{\mathrm{s}}^{-1}$. This boundary flux remains fixed even if pressure changes the subsequent radial velocity profile.

In the stated gravity-dominated approximation, neglect the pressure force in the radial momentum equation. A steady outward stellar wind then obeysAt $R_e=R_{\mathrm{esc}}$, the escape velocity is $v_e=\sqrt{2G{M}_{\ast }/R_e}$, so $C=0$. The ballistic marginal-escape stellar wind hasThe outward flow decelerates and approaches zero speed at infinity; it is exactly at the escape threshold in this pressure-negligible model. An initial speed above that threshold gives a nonzero asymptotic speed.

Adiabaticity alone would instead conserve the full Bernoulli equation $v^2/2+h-G{M}_{\ast }/r=B$. If finite pressure is retained and $v_e$ equals the mechanical escape speed, then $B=h_e$, not zero. With $p=K{\rho }^{\gamma }$ and conservation of mass, the more general implicit solution isIt needs $h_e$ and $\gamma$ in addition to the escape-speed datum. The boxed explicit power law is the leading gravity-dominated result, not an exact solution for an arbitrary warm adiabatic wind.

Let $\mathcal{L}$ be net cooling per unit volume, positive when energy is lost. At thermal equilibrium $\mathcal{L}=0$. For a slow perturbation that maintains pressure balance, the linearized thermal equation isThe heat capacity factor is positive. A positive derivative restores equilibrium: hotter gas cools more strongly, while cooler gas is heated. This is isobaric thermal stability of a heated ideal gas.

For fully ionized hydrogen, both protons and electrons contribute to the ideal gas law, so $P=2n{k}_{B}T$. Treat the stated constant heat input $\Gamma$ as a fixed heating rate per unit volume. Substituting $n=P/(2k_{B}T)$ into the cooling givesAt equilibrium $T_0$, its derivative isWith a temperature-independent recombination coefficient, the derivative is $-\Gamma /T_0$: the equilibrium is thermally unstable to isobaric perturbations. More generally, any nonincreasing ${\alpha }_B(T)$ also gives instability; stability requires its logarithmic slope to exceed one. Heating at fixed volume and pressure-balanced perturbations are distinct assumptions: the density falls as a parcel is heated, reducing the quadratic-density cooling.

No temperature law for ${\alpha }_B$ is supplied, so the displayed slope criterion states the general result without inventing one. A constant incident ultraviolet flux need not imply constant local volumetric heating in a self-consistent ionization-balance calculation. If the heating instead varies with density, temperature or neutral fraction, its derivative must be included; the conclusion here uses the explicit constant-$\Gamma$ model.

For a uniform sphere, $M=4\pi \rho R^3/3$ and its gravitational potential energy is $\Omega =-3G{M}^2/(5R)$. The ideal-gas kinetic energy is $E_{\mathrm{kin}}=3M{k}_{B}T/(2\bar{m})$. Equating these magnitudes using the virial theorem gives $GM/R=5k_{B}T/\bar{m}$. Eliminating $R$ yields the uniform-sphere pressure-work collapse criterionIn particular, for fixed composition,The same Jeans mass scaling follows from $c_s^2{k}^2\sim 4\pi G\rho$ in the Jeans instability dispersion relation and a characteristic size $k^{-1}$. Pressure communicates support in a sound-crossing time; sufficiently strong self-gravity wins when the gravitational time is shorter.

The numerical coefficient depends on the equilibrium, boundary and wavelength-to-cloud conventions. A sphere with both exactly uniform density and uniform temperature is not an exact hydrostatic self-gravitating ideal-gas equilibrium without additional stresses. The uniform-sphere argument establishes the characteristic mass scaling; it does not prove an exact stability boundary for that nonexistent equilibrium.

Hydrostatic equilibrium and dimensional scaling give characteristic stellar pressure $p\sim G{M}^2/R^4$ and mass density $\rho \sim M/R^3$. For a fixed-index stellar polytrope with $p=K{\rho }^{1+1/n}$, eliminating density gives the polytropic mass-radius relationThe dimensionless coefficient is determined by the corresponding Lane-Emden equation and boundary conditions. Both the index and $K$ must be held fixed to infer this mass dependence.

Nonrelativistic electron degeneracy has $n=3/2$, hence $R\propto M^{-1/3}$ for cold fixed-composition white dwarfs: adding mass makes the object smaller. Relativistic electron degeneracy approaches $n=3$; the mass scale becomes proportional to $(K/G{)}^{3/2}$ independently of the radius, giving the ideal Chandrasekhar limit. The limiting model signals loss of the ordinary stable cold sequence; it does not determine a unique radius at that mass.

For a gas-pressure-supported star, the virial theorem instead gives characteristic temperature $k_{B}T\sim G\bar{m}M/R$. Specifying a common characteristic temperature gives $R\propto M$, but stellar temperatures are not generally identical. A main-sequence star also requires nuclear energy generation and radiative or convective transport. Consequently there is no universal mass-radius exponent applicable to every main-sequence mass and composition. An isentropic convective model and a cold degeneracy-supported model can share index $3/2$ while having quite different physical meanings and values of $K$.

For steady inviscid motion of a barotropic fluid in a conservative gravitational potential $\Phi$, the Euler equations for an inviscid fluid imply the Bernoulli equationalong each streamline. The barotropic enthalpy function $h$ has units of energy per unit mass. For an isentropic ideal gas, $p=K{\rho }^{\gamma }$ gives $h=\gamma p/[(\gamma -1)\rho ]$ up to an additive constant. For constant density, the corresponding term is $p/\rho$. Pressure work converts between kinetic, enthalpy and gravitational energy.

A converging subsonic nozzle can accelerate a fluid at the expense of its enthalpy; in the incompressible horizontal approximation, higher speed accompanies lower static pressure. Together with conservation of mass, this explains a Venturi effect measurement. Compressible nozzle flow additionally needs its equation of state and sonic conditions. In a stellar wind, the potential term is $-GM/r$, so enthalpy or mechanical energy can be spent escaping gravity. A pressure-negligible marginal-escape wind has $v^2/2-GM/r=0$.

The constant need not be the same on different streamlines. Irrotational barotropic flow in a connected region gives a common value; a general rotational flow does not. A pump, viscosity, a radiating shock, or heat addition changes the relevant energy balance. In particular, adiabaticity alone does not allow one to drop the enthalpy term.

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.

Work in an inertial centre of mass frame and write $I={\sum }_i{m}_i|{\mathbf{r}}_i{|}^2$. Differentiation and Newton's second law giveThe pairwise Newtonian gravitational potential energy isFor each pair, combining the two force contributions gives $-G{m}_i{m}_j/|{\mathbf{r}}_i-{\mathbf{r}}_j|$. Therefore ${\sum }_i{\mathbf{r}}_i\cdot {\mathbf{F}}_i=\Omega$, andFor a stationary configuration, or after a long-time average of a bounded system whose endpoint term vanishes, the virial theorem givesThe brackets may be suppressed when $T,\Omega$ denote equilibrium averages. A collapsing cloud or an eccentric Kepler orbit need not obey $2T+\Omega =0$ instantaneously. Confining surface pressure, magnetic stress or other external forces add their own virial terms; the displayed result concerns an isolated gravitational particle system.

Let all velocity amplitudes refer to the same positive normal direction, from medium 1 into medium 2. For a right-going linear sound wave, its pressure amplitude is $p^{\prime }=Z{u}^{\prime }$ with acoustic impedance $Z=\rho c_s$; for a left-going wave it is $p^{\prime }=-Z{u}^{\prime }$. Continuity of normal velocity and pressure at the interface givesEliminating $u_t$ yields the normal-incidence acoustic velocity reflection coefficientThe common adiabatic index and mean particle mass $\bar{m}$ give $c_{s,j}=\sqrt{\gamma k_B{T}_j/\bar{m}}$, soA negative coefficient reverses the reflected velocity perturbation's phase in this fixed coordinate system. The reflected pressure coefficient has the opposite sign. Pressure balance also gives ${\rho }_1{T}_1={\rho }_2{T}_2$, so the impedance ratio can equivalently be written $Z_2/Z_1=\sqrt{{\rho }_2/{\rho }_1}$. These are normal-incidence, small-amplitude results without a surface mass or tension.

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 local thermal imbalance can amplify a gas perturbation even without an initial displacement of an interface. Let $\mathcal{L}(\rho ,T)$ denote net cooling per unit mass, including heating with a negative sign, and suppose $\mathcal{L}({\rho }_0,T_0)=0$. When the sound-crossing time is short compared with the thermal evolution time, pressure equilibrates and an ideal gas perturbation is approximately isobaric. The linearized thermal equation isThus a negative constant-pressure derivative produces isobaric thermal instability: a cooler parcel experiences a net loss that cools it further, while a warmer parcel heats further. For example, isobaric thermal stability of a heated ideal gas with $\mathcal{L}=A\rho T^{\alpha }-H$ gives $\sigma =(1-\alpha )H/(c_p{T}_0)$ and instability for $\alpha <1$. Different heating laws and fixed-density perturbations give different criteria; thermal conduction damps sufficiently short-scale temperature contrasts. Such feedback helps a gas separate into denser cool and rarer warm regions when stable thermal branches exist.

The Rayleigh-Taylor instability instead releases gravitational potential energy when a heavier fluid lies above a lighter fluid in a downward gravitational field. For two inviscid incompressible half-spaces with no surface tension or magnetic field, a horizontal perturbation of wavenumber $k$ hasThe Atwood number measures the density contrast. A descending heavy-fluid finger and an ascending light-fluid bubble both lower the total gravitational potential energy, increasing the original displacement. Interchanging the layers gives an oscillatory stable interface. An accelerating or decelerating astrophysical shell has an analogous effective acceleration, allowing mixing at suitable ejecta interfaces in a supernova remnant. Finite interface width, viscosity, surface tension and magnetic tension modify the short-wavelength behavior; the ideal sharp-interface formula does not select a finite fastest-growing wavelength.

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.

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.