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$.

Loading noise is the acoustic dipole radiation of the force exerted by a body on its surrounding fluid. For a compact, slowly moving source, the radiating acoustic density perturbation is $n\cdot \dot{\mathcal{F}}(t-R/c_0)/(4\pi c_0^{3}R)$, where $\mathcal{F}$ is the total force. If this total is constant, finite source-delay corrections can still create higher acoustic multipole sources.

Acoustic sources may be organized by the number of spatial derivatives in their effective wave equation forcing. A direct scalar source is an acoustic monopole, a divergence of a force density an acoustic dipole, and a double divergence of a stress tensor an acoustic quadrupole.

Two diametrically opposite rotating blades with equal tangential loads have zero total rotating force. Their leading compact acoustic dipoles cancel, but their retarded positions and moving-surface retarded Jacobians leave an acoustic quadrupole. For line load $2Fr/a^2$ and angular speed $\Omega$, the first radiating density is $2F{\Omega }^{2}a{\sin }^2\theta \sin (2\Omega {\tau }_0)/(3\pi c_0^{4}R)$. The twofold geometry gives frequency $2\Omega$, while projection into the rotor plane gives the amplitude directivity ${\sin }^2\theta$.

Mass and momentum flux defects at a moving interface add surface acoustic monopole and acoustic dipole terms to Lighthill acoustic analogy. For a physical conservative fluid shock wave, the Rankine-Hugoniot conditions set both defects to zero; the shock still affects distributional derivatives of the Lighthill stress tensor. A discontinuity alone is not independent mass or force injection.

The Lighthill acoustic analogy for a moving body has a volume acoustic quadrupole, a surface acoustic monopole from acoustic thickness noise, and a surface acoustic dipole from acoustic loading noise. For an impermeable moving surface, with outward normal $n$, normal surface speed $v_n$, and fluid and body normal velocities equal, the surface mass coefficient is ${\rho }_0{v}_n$ and the surface loading is the force exerted on the fluid. Schematically, the density-source equation isThe density perturbation is understood with the chosen interior extension. Permeable-surface versions have additional mass and momentum flux terms. Source approximations must distinguish local thickness radiation from cancellation of its compact net-volume contribution.

For a fixed impermeable object, the Ffowcs Williams-Hawkings equation has a volume acoustic quadrupole and a surface acoustic dipole, with no mass-injection term. If $F$ is the total force exerted on the fluid and $Q=\int TdV$ the integrated Lighthill stress tensor, the acoustic compact-source approximation in the acoustic far field gives the displayed pressure at retarded time. A force convention directed onto the body reverses the loading sign. The retarded acoustic Green function carries $1/c_0^2$ when the wave operator is ${\partial }_t^2-c_0^2\Delta$.

If the continuity equation has mass source $M$ and the conservative momentum equation has momentum source $S_i$, elimination of momentum gives the displayed wave equation. With no imposed momentum source, its extra term is ${\partial }_{t}M$. A compact source of total mass rate $\mathcal{M}=\int MdV$ radiates the acoustic monopole ${\rho }_M^{\prime }=\dot{\mathcal{M}}(t-r/c_0)/(4\pi c_0^{2}r)$. The momentum source supplies an additional acoustic dipole; specifying mass injection alone does not fix that dipole.

Interpret the surface delta as ${\delta }_s=|\nabla F|\delta (F)$, so integration against it gives a surface integral. If $F$ is signed distance, this equals the notation $\delta (F)$ directly. This convention is needed because an arbitrary rescaling of the level-set function must not change the physical source.

Define the volume flux $\mathcal{V}(t)={\int }_S{V}_{n}dS$ and total pressure force $f(t)={\int }_{S}pndS$. Here $\mathcal{V}$ is the rate of change of body volume, not the volume itself. The surface acoustic monopole contributes, in the acoustic compact-source approximation,The surface acoustic dipole contributesThe second term is a nonradiating near-field term and is smaller in the acoustic far field. Writing $\mathcal{F}(t)=x\cdot f(t)$ with the observation point fixed, the leading augmentation is thereforeThe printed dipole coefficient is missing a factor $1/c_0$. The extra inverse speed follows from differentiating $\tau =t-r/c_0$, and is required dimensionally: the printed expression with $c_0^2$ has the units of density times speed. This corrected acoustic loading noise term also agrees with the far-field acoustic force and stress moments. The monopole coefficient is unchanged.

Let $R=|\mathbf{x}|$, $\mathbf{n}=\mathbf{x}/R$ and $t_R=t-R/c_0$. Convolving the source terms in the Lighthill acoustic analogy with the retarded acoustic Green function, and moving their spatial derivatives outside the integral, givesThe integrals have compact spatial support. For an acoustically compact source, its size $l$ satisfies $\omega l/c_0\ll 1$, so the source-dependent retardation is negligible at leading order. Also take $R\gg l$. Define $S_{ij}=\int T_{ij}{d}^{3}y$ and $F_i=\int f_i{d}^{3}y$. In the radiation terms the spatial derivatives act on $t_R$, with ${\partial }_{x_i}{t}_R=-n_i/c_0$. Derivatives of $1/R$ produce smaller near-field terms. Thus the far-field acoustic force and stress moments giveUsing $n_i=x_i/R$ gives the equivalent powers of $R$ requested. The force term is an acoustic dipole; the Lighthill stress tensor produces an acoustic quadrupole.

For a small body in low-speed flow, put $M=u/c_0\ll 1$. A typical unsteady force has size $F=O({\rho }_0{u}^2{l}^2)$ and time scale $l/u$, whereas $S=O({\rho }_0{u}^2{l}^3)$. ThereforeThese are the acoustic dipole Mach-number scaling and compact acoustic quadrupole Mach-number scaling. Since $p^{\prime }=c_0^2{\rho }^{\prime }$ and outgoing acoustic intensity is of order $p^{\prime 2}/({\rho }_0{c}_0)$, their sound powers scale asThus an unsteady force on the body can produce the dominant low-Mach number sound. If the integrated force or its time derivative vanishes, that dipole contribution is absent. Force on the body and force on the fluid have opposite signs.

Let the fluid approach a uniform, quiescent state $({\rho }_0,p_0,c_0)$ at infinity. Use the outward unit normal vector $N=\nabla f/|\nabla f|$, pointing from the object into the fluid. Write $H=H(f)$ for the Heaviside step function, $q=\rho -{\rho }_0$, and ${\delta }_S=\delta (f)|\nabla f|$. With the viscous stress tensor ${\tau }_{ij}$, defineThus $T_{ij}$ is the Lighthill stress tensor, and $L_i$ is the surface loading exerted on the fluid. The total force on the object has the opposite sign. Extend $q$ by zero inside the object. Since the boundary is fixed and impermeable, $u\cdot N=0$, and the distributional mass conservation and momentum conservation equations areFor example, the surface contribution to the momentum equation is $(\rho u_i{u}_j+P_{ij})N_j{\delta }_S=L_i{\delta }_S$. Differentiating the first equation in time and subtracting the divergence of the second gives the fixed-body Ffowcs Williams-Hawkings equation:There is no acoustic thickness noise: a fixed impermeable boundary injects no mass. Apply the retarded acoustic Green function $\delta (t-R/c_0)/(4\pi c_0^{2}R)$, where $R=|x-y|$. The supplied surface-delta identity converts the loading term to a surface integral. For an exterior observer, the required integral equation isThis acoustic analogy is an exact rearrangement of the fluid equations; source approximations enter later. In the distant linear acoustic region, $p^{\prime }=c_0^{2}q$.

The acoustic far field requires observer distance $r$ much larger than the source extent $\ell$ and $k_{0}r\gg 1$, so derivatives of the retarded phase dominate derivatives of $1/R$. The acoustic compact-source approximation requires $k_0\ell \ll 1$, so the source's differential propagation delays can be neglected. These conditions concern different lengths and neither implies the other. Put $n=x/r$, $t_r=t-r/c_0$, and define the loading and stress momentsThe leading compact acoustic far field isThe surface term is an acoustic dipole; the volume term is an acoustic quadrupole. The plus sign of the loading term follows from the minus sign of its spatial divergence and the derivative of retarded time.

For the power comparison, assume low Mach number $M=U/c_0$, high Reynolds number, eddy size $\ell$, characteristic time $\tau =\ell /U$, and stress magnitude ${\rho }_0{U}^2$. Keep the prescribed turbulent region and its leading volume stress comparable when inserting the object; neglect changes of source statistics and exceptional cancellations of the leading moments. Compactness follows from $\ell /(c_0\tau )=M\ll 1$. Since $Q\sim {\rho }_0{U}^2{\ell }^3$ and the outgoing acoustic intensity is $⟨p^{\prime 2}⟩/({\rho }_0{c}_0)$, the volume power isIf the object's size is comparable to $\ell$ and its fluctuating resultant loading is $F\sim {\rho }_0{U}^2{\ell }^2$, thenAngular factors and dimensionless loading coefficients have been omitted. The acoustic dipole and acoustic quadrupole interference vanishes in the leading power integrated over a full sphere, by odd angular parity; individual directions can show interference.

A size threshold also requires a loading model. For a small nonseparating object of size $a\ll \ell$ in an approximately inviscid, smoothly varying incident eddy, a spatially uniform pressure has zero resultant on its closed surface. The first surviving loading is the pressure-gradient or added mass force, $F\sim {\rho }_0{a}^{3}U/\tau ={\rho }_0{U}^2{a}^3/\ell$. Under these explicitly stated assumptions,The two mechanisms become comparable at $a\sim \ell M^{1/3}$; for smaller objects the cases with and without the object have the same leading order of power. This is the small-body loading-noise threshold. It is not a geometry-independent size law: a separated drag force $F\sim {\rho }_0{U}^2{a}^2$ fluctuating on the original eddy time instead gives ${\mathcal{P}}_D/{\mathcal{P}}_Q\sim (a/\ell {)}^4/M^2$ and $a\sim \ell \sqrt{M}$. If both force and time are set by body-scale eddies, $\tau =a/U$, their loading power is ${\rho }_0{U}^6{a}^2/c_0^3$ and comparison with the original volume source gives $a\sim M\ell$. Each estimate states which source time and force it uses.

The physical distinction is momentum conservation. Unforced turbulence in a fluid supplies internal stresses, whose net-force contribution cancels and leaves acoustic quadrupole radiation. The supported object exchanges momentum with the fluid and the external support, permitting a time-dependent resultant force and acoustic dipole radiation. One fewer retarded derivative makes the loading radiation more efficient at low Mach number.

Put ${\rho }^{\prime }=\rho -{\rho }_0$, and take the retarded acoustic Green function for ${\partial }_t^2-c_0^2{\nabla }^2$:This selects the causal, outgoing acoustic density perturbation, with no additional incoming homogeneous wave equation solution. Convolving the acoustic dipole forcing with this Green function and moving its spatial derivative outside the integral givesThe surface delta distribution converts the spatial integral to the moving surface. At fixed surface labels $(p,q)$, setThe radial Mach number enters the moving-surface retarded Jacobian, because $dr/d\tau =-\hat{r}\cdot v$ and henceLet ${\tau }^{\ast }$ be a root of the retarded time equationThe delta change-of-variable rule now givesHere $h_p{h}_q$ is the surface-area factor from the orthogonal surface coordinates. For a subsonic surface, $|v|<c_0$, the retarded equation is monotone in $\tau$ and has one root when the motion is defined for the required past times. For more general motion, sum the displayed contribution over all simple retarded roots. A root with $1-M_r=0$ requires a separate limiting treatment; the simple-root formula does not apply there.

In the Ffowcs Williams-Hawkings equation, the other source types are a surface acoustic monopole associated with acoustic thickness noise, and a volume acoustic quadrupole involving the Lighthill stress tensor. On an impermeable material surface, fluid and surface normal velocities agree. There is no through-surface mass-flux source; the remaining thickness source is $Q={\rho }_0{v}_n$. For a rigid body, the leading acoustic compact-source approximation to that source hasThus there is no leading net-volume acoustic monopole. To neglect thickness radiation beyond that leading cancellation, assume negligible volume displacement, as for ideal thin blades, or that its higher multipoles are small compared with the retained acoustic loading noise. Rigidity alone does not make a moving finite-volume body's local thickness source identically zero.

The volume acoustic quadrupole may be neglected for low Mach number motion when exterior turbulent or nonlinear stresses do not provide a competing strong source. We also assume small linear acoustics perturbations, a uniform reference sound speed, and negligible relevant viscous and entropy sources. These are source-strength approximations, particularly important if a loading contribution itself cancels by symmetry. Under them, the retained acoustic dipole is the force exerted by the object on the fluid, with the sign used in the previous solution.

Let $R=|x|$, $n=x/R$, and ${\tau }_0=t-R/c_0$. In the acoustic far field, $R$ is large compared with the object and $\omega R/c_0\gg 1$. For a source of size $\ell$ with $\omega \ell /c_0\ll 1$ and small surface Mach number, source-dependent delays and the Doppler factor can be neglected to leading order. The surface integral then contains just the total force $\mathcal{F}(\tau )={\int }_{S}FdS$. Differentiating its retarded time, rather than its $R^{-1}$ spreading factor, gives the radiating termThe sign follows from ${\partial }_{x_i}{\tau }_0=-n_i/c_0$. Differentiating $R^{-1}$ or the direction $n$ instead produces the lower-order near field. A constant total force does not radiate at this leading compact order.

The rotation introduces angular frequency $\Omega$ and, after combining the two blades, harmonics such as $2\Omega$. The acoustic compact-source approximation requires the propagation time $a/c_0$ to be small compared with ${\Omega }^{-1}$. ThusIt also bounds every blade element's Mach number by $\mu$. Consequently $|1-M_r|=1+O(\mu )$, and its leading value is one. The separate acoustic far field condition is $\Omega R/c_0\gg 1$.

Choose the positive rotation sense so that a first blade at phase $\alpha$ has radial and tangential unit vectorsIntegrating the given line force from $r=0$ to $a$ yields ${\mathcal{F}}_1=F{e}_{\varphi }(\Omega \tau )+D{e}_x$. Its axial component is constant, while ${\dot{\mathcal{F}}}_1=-F\Omega e_r$. Since $n=\cos \theta e_x+\sin \theta e_y$, the compact acoustic dipole sound from this blade isThe other blade has phase $\alpha +\pi$ and contributes the opposite rotating force. At a common compact retarded time, their total force is $2D{e}_x$, soThis cancellation calls for the next source-delay correction; it does not mean that the complete moving-source field vanishes.

Apply the moving-interface conservation jump identity first to mass, then to each momentum component. Write $j_i=\rho u_i$ and ${\Pi }_{ij}=\rho u_i{u}_j+p{\delta }_{ij}-{\sigma }_{ij}$, with all quantities understood piecewise on the two sides. Define the jumps of flux relative to the moving shock wave byThe global distributional conservation equations areDifferentiate the first in time and subtract the divergence of the second. With ${\rho }^{\prime }=\rho -{\rho }_0$ and the piecewise Lighthill stress tensor $T_{ij}=\rho u_i{u}_j+(p^{\prime }-c_0^2{\rho }^{\prime }){\delta }_{ij}-{\sigma }_{ij}$, this givesDerivatives act on the complete distributions, including their moving support. The acoustic quadrupole term represents momentum-stress fluctuations throughout the volume. The time derivative of the surface mass-flux defect is an acoustic monopole, representing injection or removal of mass/volume. The divergence of the surface momentum-flux defect is an acoustic dipole, representing a force sheet. This is the distributional acoustic analogy across a moving interface.

For an actual freely propagating fluid shock wave with no singular mass or momentum supply, the Rankine-Hugoniot conditions give $Q=0$ and $L_i=0$. Such a shock does not acquire independent monopole and force-sheet sources merely because it is discontinuous. Its effects remain in the distributional derivatives of $T_{ij}$, including singular derivatives of its jump. Nonzero surface sources are appropriate for an interface with exchange/forcing or for a formulation that omits one side of the fluid.

A shock does not, by itself, justify retaining the $m^4$ scaling. That estimate required a low-Mach number stress $O({\rho }_0{U}^2)$ varying on the slow time $\ell /U$. Fast shock motion, short time scales, or thermodynamic deviations can invalidate that estimate. If these same compact, slow-source assumptions remain valid for the integrated Lighthill stress tensor, its quadrupole estimate still follows, even distributionally. There is no universal replacement power deducible from the mere presence of a shock; nor should vanished physical flux defects be treated as additional independent radiation sources.

For a nonseparating small body in a smooth nearly inviscid eddy of speed $U$ and length $\ell$, the fluctuating resultant is an added mass or pressure-gradient force $F\sim {\rho }_0{a}^3{U}^2/\ell$. A uniform pressure integrates to zero over a closed surface. With time $\ell /U$, its acoustic dipole power divided by the eddy's acoustic quadrupole power is $(a/\ell {)}^6/M^2$, where $M=U/c_0$ is the Mach number. Equal powers therefore give the displayed threshold. An area-scaled separated drag force with the same eddy time instead gives $a/\ell \sim \sqrt{M}$; a body-scale time gives a different criterion. The threshold depends on the specified loading and time models.

For a small rigid sphere translating with velocity $\mathbf{V}(t)$, the integrated thickness source ${\rho }_0\int \mathbf{V}\cdot {\mathbf{n}}_{s}dS$ vanishes. The first spatial moment is instead ${\rho }_0\int y_j\mathbf{V}\cdot {\mathbf{n}}_{s}dS={\rho }_0\mathcal{V}{V}_j$, where $\mathcal{V}$ is the sphere volume. Expanding the retarded time across the compact source gives the displayed first nonzero thickness field. It is an acoustic dipole, at the same order as acoustic loading noise for a translating sphere. The added mass of a sphere supplies half this volume-mass coefficient through the loading term; their sum is $3{\rho }_0\mathcal{V}\mathbf{n}\cdot \ddot{\mathbf{V}}/(8\pi c_0^{3}r)$.