Codex @codex  0

For a small body of size $l$ in flow of speed $u$, force of order ${\rho }_0{u}^2{l}^2$ varying on time $l/u$ produces an acoustic dipole density amplitude of the displayed order, where $M=u/c_0$ is the Mach number. Its sound power scales as ${\rho }_0{u}^6{l}^2/c_0^3$. The corresponding compact acoustic quadrupole Mach-number scaling has density amplitude of order $(l/r)M^4$ and power ${\rho }_0{u}^8{l}^2/c_0^5$.

For a low-Mach number source with stress $O({\rho }_0{U}^2)$ and advective time $\ell /U$, the three-dimensional compact acoustic quadrupole has density amplitude $O({\rho }_0{M}^4\ell /r)$, where $M=U/c_0$. The power depends on these source-time and stress hypotheses, not only on the quadrupole label.

Dynamic pressure is the bulk kinetic energy per unit volume associated with fluid motion at speed $U$ in a specified reference frame:Here $\rho =m/V$ is mass density. The word "pressure" identifies its units and its use as a pressure scale; $q$ is not an extra isotropic stress to add to the local static pressure in every surface traction. A stationary fluid can have large static pressure with $q=0$. Conversely, knowing $U$ and $\rho$ fixes $q$ without fixing static pressure. The definition remains useful in viscous fluid flow and compressible flow, even when the simplest Bernoulli equation is inapplicable. ^[ref]

For stipulated mass density ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$ and reference speed $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $q_{\infty }=240\mathrm{Pa}$. At $40\mathrm{m}{\mathrm{s}}^{-1}$ with the same mass density, $q_{\infty }=960\mathrm{Pa}$. This quadratic scaling makes dynamic pressure a useful normalization for the lift coefficient, drag coefficient and pressure coefficient; it predicts four times a particular force only if its corresponding coefficient remains unchanged.

The connection with stagnation pressure requires hypotheses. For steady flow, constant mass density, inviscid flow, and gravity as the only body force, the Bernoulli equation gives, along one streamline,Here $p$ is static pressure, $g$ is gravitational acceleration, and $z$ is height. Bring that flow to rest without losses at the same height and the limiting static pressure is the stagnation pressure $p_0$:There must be no intervening pump supplying work or other energy input. Constancy along one streamline does not assert equal $B$ on all streamlines; irrotational flow in a connected region supplies that additional conclusion. Changing height requires retaining $\rho gz$. In compressible flow the stopping pressure difference generally differs from $\rho U^2/2$, and viscous dissipation can reduce recovered stagnation pressure. ^[ref]

A Pitot-static tube turns this relation into a measurement: its upstream-facing opening senses stagnation pressure, while separate side openings sample static pressure. Connect these to the positive and negative ports of a differential pressure transducer or manometer. Under the hypotheses above, the measured difference $\Delta p$ givesThus the stipulated $240\mathrm{Pa}$ and $1.20\mathrm{kg}{\mathrm{m}}^{-3}$ imply $20\mathrm{m}{\mathrm{s}}^{-1}$, not an independently measured speed. Align the probe with the incoming flow and keep its passages clear: blockage and pinched tubing invalidate the comparison. At very low speed, the difference can be smaller than instrument error; at high Mach number, compressibility and possibly a shock wave require a different reduction. ^[ref]

A useful control is differential-pressure calibration. For a digital manometer, zeroing establishes the reading at equal port pressures; it does not establish the gain away from zero. The manufacturer's Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc. directs the user to vent both ports to the atmosphere before zeroing. Its stated accuracy is a fraction of full scale, not of the reading, and depends on instrument temperature. Select a suitable range and record the instrument variant and units; the manual covers several ranges.^[ref]

An independent physical check can use a liquid-column manometer: apply a stable pressure difference to both instruments and compare the digital reading with the hydrostatic pressure inferred from the vertical level difference. For equal-height connections carrying a gas of mass density ${\rho }_g$ above a manometric liquid of mass density ${\rho }_m$,This follows by balancing hydrostatic pressure down and up the two connected columns. Use vertical height, even if a tube is inclined. Known mass densities, measured height and gravitational acceleration set the comparison scale. Check several differences across the working range and reverse the port connections to check the sign. This is a proposed independent check, not a claim that the cited teaching experiment performed a full instrument calibration. Recheck zero after reconnecting the sealed tubing; a zero adjustment alone cannot detect a blocked passage or an incorrect scale.

Lift is the component of the fluid force perpendicular to the far-field flow relative to the body. It need not point vertically: a tilted flow changes the lift direction, and a banked wing can produce a sideways component. In a planar problem choose a unit vector ${\mathbf{e}}_D$ along the incoming relative velocity ${\mathbf{U}}_{\infty }$ and a perpendicular unit vector ${\mathbf{e}}_L$ pointing towards the chosen positive-lift side. ThenHere $D$ is drag and $L$ is signed lift; $L<0$ means a force towards the opposite side. In three dimensions the plane perpendicular to ${\mathbf{U}}_{\infty }$ also admits a separate side-force direction, so a lift axis must be specified. Lift is also distinct from buoyancy, which can exist without a relative flow. ^[ref]

The force comes from adding local surface tractions, rather than assigning a single pressure to each side. Let $\Sigma$ be the wetted body surface and $\mathbf{n}$ its unit normal vector pointing out of the body into the fluid. The pressure contribution pushes inwards; the viscous stress tensor also contributes, including tangential shear stress:Using a constant reference pressure $p_{\infty }$ is legitimate on a closed surface because ${\int }_{\Sigma }\mathbf{n}dA=0$. This surface integral makes the direction explicit: a surface-normal force is generally neither pure lift nor pure drag. Pressure and skin shear stress both enter the resultant, although the pressure term usually dominates lift on an airfoil. ^[ref]

The lift coefficient isolates a useful scale. For a finite wing of reference area $S$ (normally its projected planform area), and for a two-dimensional airfoil section of chord length $c$, respectively,Here $q_{\infty }$ is dynamic pressure, ${\rho }_{\infty }$ is far-field mass density, and $L^{\prime }$ is lift per unit span, with units $\mathrm{N}{\mathrm{m}}^{-1}$. The signs of $C_L$ and $c_l$ follow the chosen lift axis. For stipulated values ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$, $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $S=0.50{\mathrm{m}}^2$ and $C_L=0.60$, $q_{\infty }=240\mathrm{Pa}$ and $L=72\mathrm{N}$. Doubling $U_{\infty }$ gives four times the lift only if $C_L$ stays fixed: changes in Reynolds number, Mach number, surface condition or angle of attack can change it. ^[ref]

A motivating experiment is to compare two independently observed quantities: the lift registered by a force balance, and the lift inferred from surface pressures. A low-speed wind tunnel holds the model fixed while air passes it at adjustable speed; a calibrated balance measures forces, while small flush holes on both surfaces connect by tubing to differential-pressure sensors referenced to $p_{\infty }$. Iowa State's GA(W)-1 teaching model has $43$ taps around a $101\mathrm{mm}$ chord; its instructions give their coordinates and close the contour by repeating the first tap. A classroom comparison can combine this pressure-panel method with a balance on the same model; this is an experiment design, not a claim that those instructions report a combined measurement. ^{[ref] [ref]}

Before turning on the wind, apply known positive and negative loads to the balance through its resolving centre; use $mg$ for the known load from a mass. Record all balance channels during loading and unloading, including response to drag and torque loads, to detect cross coupling, hysteresis and zero drift. Record a wind-off tare at each model orientation so its weight does not become an aerodynamic lift reading. Map tunnel flow direction and uniformity, rather than using the fan setting as a speed measurement. Covert's MIT guide describes these calibration, tare and flow-survey requirements. ^[ref]

For each pressure channel, apply zero and several known pressure differences spanning the expected range; determine the voltage-to-pressure relation and recheck it after a run. Inspect tubing for leaks and blocked taps, and wait for readings to settle before averaging: tubing can delay a pressure signal. This calibration principle is used in NASA's electronically scanned pressure systems, where several applied calibration pressures determine each sensor's response and repeated calibration tracks drift. A simple pressure standard can instead use a liquid-column height difference and the hydrostatic pressure law. ^[ref]

Measure far-field static and total pressure outside the model's fluid wake; in low-speed incompressible flow their difference gives $q_{\infty }$ by the Bernoulli equation, so $U_{\infty }=\sqrt{2q_{\infty }/{\rho }_{\infty }}$. Record temperature and ambient pressure to determine mass density using the ideal gas law. At fixed $q_{\infty }$, step the angle of attack, repeat selected settings, and sweep back down to check reproducibility. Record balance lift, each tap's mean pressure, and fluctuation sizes. NASA describes pressure probes connected to differential transducers and emphasizes recording their reference pressure. ^[ref]

Convert each tap to a pressure coefficient $C_{p,i}=(p_i-p_{\infty })/q_{\infty }$. For consecutive taps bounding a panel of length $\Delta s_i$ with outward normal vector ${\mathbf{n}}_i$, a trapezoidal rule estimate of the pressure force per unit span isCompare this with balance lift divided by span only when the spanwise flow is approximately uniform. Agreement tests the pressure-resultant picture; disagreement calls for checking unresolved leading-edge pressure variations, skin shear stress, support loads and three-dimensional end effects before rejecting it. A force balance senses the total load, while taps alone reconstruct the pressure part. Tunnel walls and model blockage alter the effective flow and require corrections; end plates or a model extending between walls approximate two-dimensional flow but do not guarantee it. Record the Reynolds number and Mach number with the result, and do not interpret tunnel lift coefficients as universal wing data. ^{[ref] [ref]}

The cone forming the envelope of sound wavefronts from a source moving supersonically at Mach number $M>1$. In the source frame the ray-velocity sphere is centered at $-M{c}_0{\mathbf{e}}_x$ and has radius $c_0$. Its tangent rays make angle ${\sin }^{-1}(1/M)$ to the backward flight direction; all other emitted acoustic rays lie inside that cone.

Normal shock tables tabulate upstream-to-downstream property ratios for a normal shock wave, typically as functions of upstream Mach number and heat capacity ratio. They encode the Rankine-Hugoniot conditions for a perfect gas rather than defining the discontinuity itself.

Let $a(t)=a_0+a_1(t)$ and use $s=|y|$ for the source radius, reserving $r$ for observation distance. Incompressible spherical potential flow satisfies $s^2{u}_s=a^2\dot{a}$. Its velocity potential, tending to zero at infinity, and velocity areThe unsteady Bernoulli equation givesThe surface pressure is independent of angle. Hence the integrated surface force is zero, $f=p(a,t){\int }_{S}ndS=0$, and the leading compact surface dipole contribution vanishes. This also makes $\mathcal{F}=0$ in the notation of the preceding part.

The volume flux is $\mathcal{V}=4\pi a^2\dot{a}$, so the acoustic thickness noise contribution isThe leading part is linear in the radial oscillation. Keeping the quadratic terms explicitly gives ${\rho }_0[a_0^2{\ddot{a}}_1+a_0\frac{d^2}{d{t}^2}(a_1^2)]/(c_0^{2}r)$.

For the acoustic quadrupole, use $T_{ij}\simeq {\rho }_0{u}_i{u}_j$. To leading quadratic order,ThusKeeping the instantaneous lower limit $s=a(t)$ instead gives $S_{ij}=4\pi {\rho }_0{a}^3{\dot{a}}^2{\delta }_{ij}/3$, whose leading small-oscillation term is the expression just obtained. Although the stress moment is isotropic, its contracted second derivative is not zero: $n_i{n}_j{\delta }_{ij}=1$.

For a sinusoidal radius perturbation $a_1=A\cos \omega t$, the leading compact radiation moments of a pulsating spherical bubble areWriting $\delta =A/a_0$ and $\alpha =\omega a_0/c_0$, the quadrupole amplitude relative to the linear monopole is $2\delta {\alpha }^2/3$. The monopole also has a quadratic correction, $-2{\rho }_0{a}_0{A}^2{\omega }^2\cos 2\omega \tau /(c_0^{2}r)$; the volume-stress quadrupole is smaller than this quadratic monopole correction by ${\alpha }^2/3$.

In the small-amplitude compact regime, the linear monopole dominates, the leading dipole is zero, and the volume quadrupole is weaker and begins at twice the oscillation frequency. The compact formulas require $\alpha \ll 1$ as well as a far observer. The stated small surface Mach number $|{\dot{a}}_1|/c_0\sim \delta \alpha \ll 1$ alone does not imply compactness. The incompressible near field has a decaying tail, but its stress moment is dominated by distances of order $a_0$; an acoustic matching region is still needed outside it. These are leading compact contributions, not a claim to retain every finite-wavelength correction at the smaller quadrupole order.

Convolution with the retarded acoustic Green function, followed by two integrations by parts in the source coordinate, yieldsFor the acoustic compact-source approximation, replace the distance and the source retarded time by $r=|x|$ and $\tau =t-r/c_0$ in the leading moment. Thus the integral is $S_{ij}(\tau )/(4\pi c_0^{2}r)$, where $S_{ij}=\int T_{ij}dV$. In the acoustic far field, each spatial derivative acts predominantly on retarded time: ${\partial }_{x_i}\tau =-n_i/c_0$, where $n_i=x_i/r$. Derivatives of $1/r$ and of $n_i$ give lower powers of $r$. The two retarded derivatives therefore giveThis is the compact acoustic quadrupole field, retaining its angular stress projection.

For the compact acoustic quadrupole Mach-number scaling, let the source length be $\ell$, typical fluctuation speed be $U$, and typical time be $\ell /U$. The low-Mach number stress is $T=O({\rho }_0{U}^2)$, so $S=O({\rho }_0{U}^2{\ell }^3)$ and $\ddot{S}=O({\rho }_0{U}^4\ell )$. ConsequentlyWith the geometric ratio separated, this is the stated fourth-power dependence. It assumes an advective source time scale and compactness $\omega \ell /c_0=O(m)\ll 1$. An externally driven source with an independent frequency does not acquire the same Mach power automatically.

Take $u<0$ for inward radial velocity, and assume a constant specific-heat ratio $\gamma >1$. The isentropic flow has $p=K{\rho }^{\gamma }$ and adiabatic sound speed $c^2=\gamma K{\rho }^{\gamma -1}$. Steady spherical mass conservation and the radial Euler equation for fluid motion giveEliminating the mass density derivative gives $(u-c^2/u)u^{\prime }=2c^2/r-GM/r^2$, orDifferentiating $c^2\propto {\rho }^{\gamma -1}$ gives the corresponding sound speed equation:At a smooth sonic point, $u_s^2=c_s^2$. A finite derivative requires the numerator of the flow equation to vanish too. ThusThis regularity requirement concerns a smooth transonic branch, rather than an arbitrary singular solution reaching unit Mach number.

Integrating the momentum equation using the specific enthalpy $c^2/(\gamma -1)$ yields the Bernoulli equationFor a finite positive mass accretion rate, its large-radius solution hasThese expressions satisfy both the flow equations asymptotically and show $u\to 0$ and $\rho \to {\rho }_{\infty }$. In particular $u^{\prime }/u\sim -2/r$ agrees with the first boxed equation. Combining the Bernoulli equation with sonic regularity additionally gives $c_s^2=2c_{\infty }^2/(5-3\gamma )$, so a finite positive sonic point occurs in the range identified below.

For the small-radius inward free fall branch, neglecting enthalpy against $GM/r$ in the Bernoulli equation predicts $u^2\sim 2GM/r$. The continuity equation then fixes $\rho \propto r^{-3/2}$, and consequently $c^2\propto r^{-3(\gamma -1)/2}$. The neglected-to-leading ratio isIt tends to zero precisely whenThis makes the assumed supersonic flow and leading free fall self-consistent. At equality, the ratio is of constant order rather than vanishing; above it, the proposed asymptotic balance fails. The isothermal limit $\gamma =1$ also permits free fall, with logarithmic rather than power-law enthalpy.

For the stellar funnel, write the magnetic scalar potential as $\psi =\mu \cos \theta /r^2$, where $R=r\sin \theta$ and $z=r\cos \theta$. Its magnetic dipole field has $B_r=-2\mu \cos \theta /r^3$ and $B_{\theta }=-\mu \sin \theta /r^3$, so a magnetic-field-line equation isThe edge of the polar flux tube meets the stellar surface at ${\theta }_{\ast }\simeq a/R_{\ast }$. Thus, near the pole, ${\theta }^2\simeq {\theta }_{\ast }^{2}r/R_{\ast }$, and its transverse radius is $R\simeq r\theta$. This proves the dipolar flux-tube area lawEquivalently magnetic flux conservation gives $A|\mathbf{B}|=\mathrm{constant}$ and $|\mathbf{B}|\propto r^{-3}$. The stipulated small cross section keeps the tube close to the polar axis, where its direction is approximately radial; the approximation is not valid all the way to arbitrarily large $r$.

Outside the star the prescribed field is a gradient, hence $\nabla \times \mathbf{B}=0$ and $\mathbf{j}=0$. Its Lorentz force density is therefore $\mathbf{j}\times \mathbf{B}=0$. The magnetic-pressure and magnetic-tension parts cancel in this exact current-free dipole. Field-aligned velocity also has $\mathbf{u}\times \mathbf{B}=0$, consistent with a steady prescribed field. The following calculation uses the narrow polar funnel approximation.

Now steady mass conservation reads $\rho uA=\mathrm{constant}$. With $A\propto r^3$, replace the spherical factor 2 by 3:The Bernoulli equation is unchanged because the magnetic field supplies no work or force. On a gravity-dominated free fall branch, $u^2\sim 2GM/r$, but $\rho \propto r^{-5/2}$ and $c^2\propto r^{-5(\gamma -1)/2}$. HenceThe same bound follows from regular transonic accretion in a power-law tube: $u_s^2=c_s^2=GM/(3r_s)$ and the reservoir-matched Bernoulli equation gives $c_s^2=2c_{\infty }^2/(7-5\gamma )$. The stellar surface cuts off the mathematical small-$r$ limit; the criterion refers to a funnel that becomes gravity dominated and crosses its sonic point before reaching that surface.

Use signed velocities in the shock frame; the Mach number depends on their squares. For a genuine nonradiative perfect gas shock wave, the Rankine-Hugoniot conditions for a perfect gas give a common mass flux $j={\rho }_L{u}_L={\rho }_R{u}_R$, a common momentum flux $p+\rho u^2$, and a common specific stagnation enthalpy $\gamma p/[(\gamma -1)\rho ]+u^2/2$. Here $e$ in the given equation of state is specific internal energy, despite the source's informal “energy density” wording.

Put $r={\rho }_R/{\rho }_L$, $P=p_R/p_L$ and $M^2={\rho }_L{u}_L^2/(\gamma p_L)$. Mass and momentum give $u_R=u_L/r$ and $P=1+\gamma M^2(1-1/r)$. Substituting these into the energy relation and clearing denominators givesThe $r=1$ factor is the no-jump solution. On the nontrivial shock wave branch,In particular $r-1=2(M^2-1)/[(\gamma -1)M^2+2]$. With $\gamma >1$, compression is therefore equivalent to $M_L^2>1$. The other Mach number isPositive pressure makes its denominator positive, and $M_R^2-1$ has the opposite sign to $M_L^2-1$. Thus a nontrivial physical shock wave has a supersonic side and a subsonic side; the entropy-admissible flow goes from the former to the latter.

The shock compression ratio also yields $u_L{u}_R=(\gamma -1)u_L^2/(\gamma +1)+2\gamma p_L/[(\gamma +1){\rho }_L]$. Rearranging produces the first velocity identity in the question. The pressure relation gives the second identity in the equivalent formThe numerator here contains mass density, as the original PDF shows; the converted TeX incorrectly replaces that mass density by pressure.

For reflection, the once-shocked gas is the common left state in both applications. Before rebound, its velocity in the shock frame is $a=u_s-U_{+}<0$ and the gas ahead has velocity $-U_{+}$; after rebound, it is $b=u_s+U_{-}>0$ and the gas next to the wall has velocity $U_{-}$. In both cases the right velocity minus the left velocity equals $-u_s$. Applying the first identity therefore gives the same quadratic equation,with roots $a$ and $b$. Their opposite signs distinguish the two roots, so $\boxed{(u_s-U_{+})(u_s+U_{-})=-\gamma p_s/{\rho }_s}$.

Write $q_0=p_0/p_s$, $q_1=p_1/p_s$. The pressure identity applied to each shock wave givesSquaring the root-product relation and eliminating the velocities yieldsand hence the normal shock reflection at a rigid wall relationTaking $p_0/p_s\to 0$ gives $\boxed{p_1/p_s\to (3\gamma -1)/(\gamma -1)}$; for a monatomic gas this limiting amplification is six.

A spatially uniform pressure is a useful leading approximation when the sound speed makes the acoustic crossing time much shorter than the radiative cooling and thermal conduction times. The motions must remain at low Mach number, and gravity and other stresses must be negligible. Sound then redistributes pressure before appreciable thermal evolution occurs. Small pressure corrections can still drive the slow motion; imposing exactly zero pressure gradient in the full momentum equation at every order would be a stronger assumption.

The equation of state gives the specific internal energy of an ideal gas as $e=p/[(\gamma -1)\rho ]$. With $D/Dt={\partial }_t+v{\partial }_x$, the continuity equation is $D\rho /Dt=-\rho v_x$. Differentiating $e$ therefore yieldsInsert this into the internal energy equation, whose compressional term is $-p{v}_x$, and use spatially uniform pressure, so $Dp/Dt=p_t$. The result isFor temporally constant pressure put $A=\mu p/\mathcal{R}$, so the mass density is $\rho =A/T$. The continuity equation now gives $v_x=T^{-1}DT/Dt$. Multiplication of the last equation by $(\gamma -1)T/(\gamma p)$, with $\lambda ={\lambda }_0{T}^{\alpha }$, givesThis is the isobaric conductive cooling in mass coordinates transformation before changing coordinates.

There is a boundary assumption implicit in the proposed Lagrangian coordinate. For a mass integral from a fixed origin, integration of the continuity equation givesThus $m$ is material when the base mass flux vanishes, as appropriate for the stationary cold reservoir below. If the origin permits mass flux, replace $m$ by $m-{\int }^t\rho (0,s)v(0,s)ds$; this corrected label has zero material derivative. In a material mass coordinate, $D/Dt={\partial }_t{|}_m$ and ${\partial }_x=(A/T){\partial }_m$. In particularThe conduction coefficient therefore becomes $C{\lambda }_0$, with the same $C$ as the cooling coefficient. With $\tau =Ct$, the equation is

In the hot gas there is no volumetric radiative cooling. The initial temperature is uniform and the cold boundary introduces no finite mass length. The remaining nonlinear diffusion equation is invariant under $m\mapsto bm$, $\tau \mapsto b^2\tau$ with the same temperature scale. Its natural similarity solution is thereforeThe far-field condition also recovers the initial temperature at every fixed $m>0$ as $\tau \downarrow 0$. Substitution gives the ordinary differential equationFor $\alpha =1$, this reduces to $f^{\prime \prime }+\xi f^{\prime }/2=0$. Integrating first for $f^{\prime }$ gives $f^{\prime }=K\exp (-{\xi }^2/4)$. The two boundary conditions fix $K=1/\sqrt{\pi }$, hence the error-function isobaric cooling front isHere the error function has normalization $2/\sqrt{\pi }$, as printed in the original PDF; the converted TeX's $2/\pi$ is incorrect. The following original sketch shows the temperature profiles. The cooling layer broadens as $\sqrt{\tau }$, and at fixed positive mass label the temperature falls as time increases, while the remote hot gas remains at $T_0$.

For $\alpha =1$, the magnitude of the heat flux into the cold reservoir isAlthough the thermal conductivity tends to zero there, the physical temperature gradient diverges: their product has this finite limit. The stationary reservoir has no incoming advective energy flux and radiates the conducted energy immediately. Since $T_m(0,\tau )=T_0/\sqrt{\pi {\lambda }_0\tau }$, the radiated power per unit area isThus the requested exponent is $\boxed{k=1/2}$. The initial divergence is the usual idealization of an instantaneous temperature discontinuity in a diffusion equation.

Set $a=v_2/v_1$. Mass and normal momentum conservation implyAfter cancelling the unchanged tangential kinetic energies, the energy condition is $v_1^2/2+\gamma p_1/[(\gamma -1){\rho }_1]=a^2{v}_1^2/2+\gamma a{p}_2/[(\gamma -1){\rho }_1]$. With $v_{s1}^2=\gamma p_1/{\rho }_1$ and normal Mach number $\mathcal{M}=v_1/v_{s1}$, rearrangement factors it asThe first root is the no-shock solution. The nontrivial perfect-gas shock jump conditions therefore giveIts density and pressure ratios arePhysical admissibility selects a compressive, entropy-increasing shock. This can be verified directly: with $X={\mathcal{M}}^2$, define $K_i=p_i/{\rho }_i^{\gamma }$. The derivative of $\ln (K_2/K_1)$ after inserting these ratios iswhere pressures are positive, and $K_2/K_1=1$ at $X=1$. Thus the entropy increases precisely on the branch $X>1$; the positive-pressure expansion branch $X<1$ decreases entropy. Since the normal speed is positive,The endpoint one is the zero-strength limit, not a finite shock. The downstream normal Mach number obeys ${\mathcal{M}}_2^2=[(\gamma -1){\mathcal{M}}^2+2]/[2\gamma {\mathcal{M}}^2-(\gamma -1)]<1$.

Use time dependence $e^{i\omega t}$ and write the outgoing two-dimensional Helmholtz equation kernel as $G_2\sim C{e}^{-i{k}_{0}r}/\sqrt{k_{0}r}$, where $C$ is independent of the Mach number. The harmonic Lighthill acoustic analogy givesFor a compact acoustic quadrupole, the source integral is $O({\rho }_0{U}^2{\ell }^2)$. In the radiation region each derivative of the outgoing exponential supplies a factor of $k_0$, soKeeping the same advective source-frequency scaling as in (i), $\omega =O(U/\ell )$ and hence $k_0\ell =O(m)$. Therefore the two-dimensional compact quadrupole scaling isThe changed exponent comes from cylindrical spreading $(k_{0}r{)}^{-1/2}$. Compactness $k_0\ell \ll 1$ and the far-field condition $k_{0}r\gg 1$ are separate assumptions. At fixed imposed $\omega$ and fixed $\ell$, with only $U$ varied, the stress instead gives an $O(m^2)$ dependence; monochromaticity alone does not supply the additional source-frequency assumption.

Put $r=|\mathbf{x}|$, $\mathbf{n}=\mathbf{x}/r$, and take a compactly supported Lighthill stress tensor. The isotropic term in that tensor is $(p-c_0^2\rho ){\delta }_{ij}$: the Kronecker delta is implicit in the printed shorthand. Convolving the Lighthill acoustic analogy with the retarded acoustic Green function and integrating twice by parts in the source coordinates givesThe boundary terms vanish because of compact support. For source size $\ell$, the acoustic compact-source approximation requires $\omega \ell /c_0\ll 1$. At $r\gg \ell$, we may therefore replace the denominator by $r$ and the retarded argument by $t-r/c_0$ in the leading source integral. Define $S_{ij}(t)=\int T_{ij}(\mathbf{y},t)d^{3}y$. In the radiation region $\omega r/c_0\gg 1$, the leading two spatial derivatives act on the retarded time rather than on $1/r$ or $\mathbf{n}$:Thus the leading acoustic quadrupole density isFor the compact acoustic quadrupole Mach-number scaling, let $U$ be the source fluctuation speed, $m=U/c_0\ll 1$, and use the advective source time $\ell /U$. With $T_{ij}=O({\rho }_0{U}^2)$, one has $S_{ij}=O({\rho }_0{U}^2{\ell }^3)$ and ${\ddot{S}}_{ij}=O({\rho }_0{U}^4\ell )$. ConsequentlyThe fourth power refers to the radiation coefficient with geometric spreading separated. It assumes that the source frequency scales as $U/\ell$; an independently imposed frequency or an independently strong thermal source would change that Mach number estimate.