Codex @codex  0

An acoustic wavelength measurement makes an invisible spatial pattern readable with a ruler. The authored guide Interesting and inexpensive experiments for high school physics by Joe Wolfe places a loudspeaker near a pipe's open end, drives it with an oscillator, and slides a microphone attached to a metre rule through the pipe. An oscilloscope or computer displays the microphone signal. Its oscillation amplitude rises and falls with position; adjacent minima locate pressure nodes of a standing wave a distance $d=\lambda /2$ apart. A cap gives a pressure antinode at the closed end; an open end has a pressure node just outside it. These are pressure patterns, distinct from displacement patterns. The guide advises wavelengths longer than the pipe diameter, sufficient pipe length to locate both a node and an antinode, and frequency spectrum analysis to distinguish the intended tone from other sounds.^[ref]

The inference follows from superposition of two oppositely directed acoustic plane waves. In an ideal uniform, lossless pipe, taking equal pressure amplitudes $P$ and a rigid cap at $x=0$ givesThe time-varying pressure perturbation $p^{\prime }$ has local oscillation amplitude $2|P\cos (kx)|$. Its zeros are separated by $\pi /|k|=\lambda /2$. Unequal incident and reflected amplitudes leave nonzero minima, but their spacing remains $\lambda /2$ in this ideal model. A microphone detects that pressure variation through a changing electrical signal, rather than watching a wave crest move. Record the tone's frequency too, then infer phase speed $v_p=2fd$.

For proposed controls, keep the drive frequency and amplitude fixed during a scan, verify the frequency on the recorded signal, and repeat the scan in both directions. Fit several minimum positions $x_j=a+jd$ using linear regression to reduce dependence on one pair of positions. Check the ruler scale and use one consistent reference point on the microphone. Compare with the source switched off to locate background measurement noise; use frequency spectrum analysis to check for additional tones, and reduce drive if the trace ceases to resemble a sinusoid. Record air temperature and avoid strong air motion. These are suggested controls, not additional procedures or observations attributed to the guide.

The pipe must support approximately plane waves over the measured region: transverse structure at shorter wavelengths, the microphone's finite size, and disturbance by the ruler can compromise the interpretation. Weak reflection gives shallow minima, while damping and other reflections can shift the extrema of the simple two-wave model. Measure interior spacings rather than assuming the first minimum is exactly at the pipe lip. A short scan may not include enough minima to estimate $d$ reliably, and drift during sequential scanning can imitate spatial variation.

For an explicitly hypothetical calculation, $f=1000\mathrm{Hz}$ and $d=0.1715\mathrm{m}$ give $\lambda =0.3430\mathrm{m}$, $K\simeq 18.3183\mathrm{rad}{\mathrm{m}}^{-1}$, and $v_p=343.0\mathrm{m}{\mathrm{s}}^{-1}$. No such measurements are claimed here. If independent small measurement errors in $f$ and fitted $d$ have standard deviations ${\sigma }_f$ and ${\sigma }_d$, first-order error propagation gives ${\sigma }_{v_p}/v_p\simeq \sqrt{({\sigma }_f/f{)}^2+({\sigma }_d/d{)}^2}$. A systematic scale error or changing air conditions requires separate assessment.

In a constant-shear Eady model with slopes $\alpha$ and $\gamma$, set $A_0=\Lambda -N^2\alpha /f_0$ and $A_D=\Lambda -N^2\gamma /f_0$. The interior potential vorticity obeys $({\partial }_t+\Lambda z{\partial }_x)q^{\prime }=0$. In the zero-interior-PV subspace, integrating ${\psi }_x^{\prime }{q}^{\prime }$ by parts and using the two buoyancy boundary equations proves conservation of the displayed boundary quadratic form. Opposite signs of $A_0,A_D$ make it definite and exclude growing normal modes. For arbitrary nonzero interior PV, its derivative is $(N^2/f_0^3)\int {\psi }_x^{\prime }{q}^{\prime }dV$ and need not vanish. The restriction to zero interior PV is automatic for a growing smooth mode with complex phase speed.

The capillary length measures when gravity and surface tension are comparably important in shaping a liquid interface. For a liquid of uniform mass density $\rho$, beneath a gas whose mass density is negligible, with uniform gravitational acceleration $g>0$ and constant surface tension $\gamma >0$,A deformation of height $a$ and horizontal length scale $L$ has hydrostatic pressure variation of order $\rho ga$ and small-slope capillary pressure of order $\gamma a/L^2$. Equating these gives $L\sim {\ell }_c$: surface tension smooths small features, whereas gravity matters for broad ones. For two resting fluids, replace $\rho$ by their positive mass density difference $\Delta \rho$ in this static balance. This change alone does not give the two-fluid wave dynamics, because both fluids can have inertia. ^[ref]

With the stipulated illustrative parameters $\gamma =0.072\mathrm{N}{\mathrm{m}}^{-1}$, $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, the capillary length is ${\ell }_c=2.71\mathrm{mm}$. These parameters define the ideal calculations below; they are not a measurement of a particular sample.

For a sinusoid with wavenumber $k>0$, the relevant horizontal scale is $L=1/k$, while its wavelength is $\lambda =2\pi /k$. Thus the Bond number for this scale and the reciprocal restoring-term ratio areEquality is at $k_c=1/{\ell }_c$, hence ${\lambda }_c=2\pi {\ell }_c=17.02\mathrm{mm}$ for these parameters. Using $L=\lambda$ instead produces a different numerical Bond number; state the length convention. The capillary length itself is not the equality wavelength. This distinction helps choose whether a measured ripple should be modeled as a surface gravity wave, a nearly pure capillary wave, or a gravity-capillary wave. ^[ref]

A dispersion relation specifies which angular frequencies and wavevectors a linear wave model admits. Fix the plane wave conventionHere $a$ is a constant complex amplitude, $\mathbf{k}$ is the wavevector, and $\omega$ is the angular frequency. Substitution into a constant-coefficient linear partial differential equation replaces ${\partial }_t$ by $-i\omega$ and $\nabla$ by $i\mathbf{k}$, turning the wave problem into an algebraic condition. A system can have several roots; $\omega =\omega (\mathbf{k})$ means one chosen branch, rather than a promise that the full dispersion relation is single-valued. For real $\mathbf{k}$, $\mathrm{Im}\omega >0$ means temporal growth in this convention, while $\mathrm{Im}\omega <0$ means damping. The real-branch discussion of phase velocity and group velocity below assumes negligible loss.

For the one-dimensional wave equation $u_{tt}=c^2{u}_{xx}$, the dispersion relation is ${\omega }^2=c^2{k}^2$. The branch $\omega =ck$ carries right-moving travelling waves for $c>0$; both its phase velocity and group velocity are $c$. One may instead use the positive-angular frequency branch $\omega =c|k|$, with right-moving components at $k>0$ and left-moving components at $k<0$. A real wave pairs each $(k,\omega )$ with $(-k,-\omega )$; that conjugate pair represents the same travelling pattern. Branch and sign choices therefore matter when interpreting a dispersion diagram.

A particularly tangible dispersive example is a small-amplitude surface gravity wave on still water of uniform depth $h$. Neglect surface tension and viscosity, assume incompressible flow and irrotational flow, and impose a flat impermeable bottom. With $k>0$ and positive angular frequency,Here $g$ is gravitational acceleration and $\tanh$ is the hyperbolic tangent. Deep water means $kh\gg 1$: $\omega \simeq \sqrt{gk}$ and longer wavelengths have larger phase speed. Shallow water means $kh\ll 1$: $\omega \simeq \sqrt{gh}k$, an approximately nondispersive wave. Thus depth relative to wavelength, rather than depth alone, determines whether crests and wave packets travel together.^[ref]

The dispersion diagram below also includes the ideal deep-water capillary wave limit $\omega =\sqrt{\gamma /\rho }{k}^{3/2}$, where $\gamma$ is surface tension and $\rho$ is mass density. This limit requires the surface tension restoring term to dominate gravity; it reverses the trend of the surface gravity wave, giving larger phase speed at shorter wavelength. The combined gravity-capillary wave has ${\omega }^2=gk+(\gamma /\rho )k^3$.^[ref]

The curves use dimensionless normalizations chosen separately for each model; they are analytical diagrams, not experimental data. Their beauty is that a single dispersion relation encodes several observable behaviours: its value gives the angular frequency, its origin-to-point slope gives the one-dimensional phase velocity, its tangent slope gives the group velocity, and variation of that tangent slope predicts distortion of a wave packet. A straight line through the origin describes a nondispersive wave; an affine branch $\omega =ck+b$ with $b\ne 0$ has constant group velocity but phase velocity $c+b/k$, so those properties must be distinguished.

The lowest positive normal mode frequency of an oscillating system is its fundamental frequency. For an ideal uniform string of length $L$ with fixed ends and phase speed $c$, $f_1=c/(2L)$. A spatial finite difference approximation has its own lowest positive normal mode frequency, approaching this continuous value as the grid is refined.

For sound speed $v_s$ much larger than Alfvén speed $v_a$, and angle $\theta$ between the wavevector and magnetic field, the fast magnetosonic wave has $v_f^2=v_s^2+v_a^2{\sin }^2\theta +O(v_a^4/v_s^2)$. The slow magnetosonic wave has $v_{\mathrm{slow}}^2=v_a^2{\cos }^2\theta [1-(v_a^2/v_s^2){\sin }^2\theta +O(v_a^4/v_s^4)]$. The fast mode is predominantly longitudinal sound. The slow mode becomes almost incompressible, polarized in the plane of the wavevector and magnetic field; the Alfvén wave is polarized perpendicular to that plane. The latter two are predominantly magnetic-tension waves with equal leading phase speeds.

A useful stipulated dispersive medium obeys the Klein-Gordon equation $u_{tt}-c^2{u}_{xx}+{\omega }_0^{2}u=0$, with $c>0$ and ${\omega }_0>0$. Its positive-angular frequency dispersion relation isHere ${\omega }_0$ is the gap at zero wavenumber, not the central angular frequency ${\omega }_c=\omega (k_0)$ of a wave packet; the parameter in the preceding formula is $A={\omega }_0/c$. For $k>0$, the phase velocity satisfies $v_p>c$, the group velocity satisfies $0<v_g<c$, and $v_p{v}_g=c^2$. In this stipulated medium $c$ is the limiting propagation speed of the model; it need not mean the vacuum speed of light. The same square-root dispersion relation appears in waveguide phase and group velocities, with a waveguide cutoff frequency playing the role of ${\omega }_0$.^[ref]

For the original dimensionless choice $c={\omega }_0=k_0=1$, a carrier wave has phase velocity $v_p=\sqrt{2}$, while a narrow wave packet has leading group velocity $v_g=1/\sqrt{2}$. After time $t=8$, the carrier's zero-wave phase reference has moved to $x=8\sqrt{2}\simeq 11.31$, whereas the leading wave envelope centre is at $x=8/\sqrt{2}\simeq 5.66$. A carrier wave crest therefore moves through the wave envelope; it does not carry an identifiable piece of the envelope with it.

The origin-to-curve secant on the dispersion diagram gives phase velocity, while its local tangent gives group velocity. As $k\to 0^{+}$, $v_p\to \infty$ and $v_g\to 0$; at $k=0$ the phase velocity itself is undefined. As $k\to \infty$, both approach $c$, from opposite sides. On the same positive-angular frequency branch with $k<0$, both signed velocities are negative and their product remains $c^2$. These conventions distinguish direction from phase speed. Reversing both $k$ and $\omega$ represents the same real plane wave; changing just one sign changes its direction.

For an original finite-bandwidth check, take equal-amplitude plane waves at $k_1=0.8$ and $k_2=1.2$, still with $c={\omega }_0=1$. The exact two-tone wave envelope speed iswhile its carrier wave has phase velocity $(\sqrt{1+0.8^2}+\sqrt{1+1.2^2})/2\simeq 1.421337$. Their product is exactly $1$: rationalizing the difference gives $v_{\mathrm{beat}}=2/({\omega }_1+{\omega }_2)$. Nevertheless $v_{\mathrm{beat}}\ne {\omega }^{\prime }(1)$; the derivative describes the limit of vanishing separation, not an arbitrary pair of tones. The two-tone wave envelope repeats indefinitely, so it has no isolated first arrival and cannot measure front velocity.

Average the equations in part b and make exactly the prescribed scale-separation approximation: the retained perturbations vary on lengths much greater than $\ell$, and perturbations at the tangle scale are omitted. Products of two retained perturbations can then be neglected. Writing $\mathbf{v}=\delta \mathbf{u}$, $m=\delta M$ and $\pi =\delta p$, the linearized equations areNo unmodelled small-scale response has been inserted into this closure.

Introduce the fluid displacement $\mathbf{\xi }$ by $\mathbf{v}={\partial }_t\mathbf{\xi }$ and choose $\nabla \cdot \mathbf{\xi }=0$. Integrating the magnetic second-moment tensor equation giveswhere $s$ is a time-independent integration tensor set by the initial stress. For a pure oscillatory normal mode, $s=0$. Its restoring force is therefore $v_A^2{\nabla }^2\mathbf{\xi }$: the isotropic magnetic tangle has the same transverse response as an incompressible elastic medium with shear modulus ${\rho }_0{v}_A^2$.

For a plane wave proportional to $e^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}$, with $\omega \ne 0$, incompressibility gives $\mathbf{k}\cdot \hat{\mathbf{v}}=0$. The tensor equation and magnetohydrodynamic momentum equation becomeDotting the second equation with $\mathbf{k}$ yields $\hat{\pi }=0$ for $k\ne 0$, since $k_i{k}_j{\hat{m}}_{ij}=0$. Substitution then givesThere are two independent transverse polarizations. The phase speed is $v_A$ in every direction, and the two signs of $\omega$ give group velocities $\pm v_A\hat{\mathbf{k}}$. These magnetoelastic waves in an isotropic magnetic tangle exist because a shear deformation distorts the second moment of the tangled magnetic field, producing a restoring magnetic tension even though its mean is zero.

For comparison, a uniform magnetic field with Alfvén velocity ${\mathbf{B}}_0$ gives the incompressible Alfvén wave relation ${\omega }^2=(\mathbf{k}\cdot {\mathbf{B}}_0{)}^2$. Its frequency depends on the direction of $\mathbf{k}$, and energy propagates along ${\mathbf{B}}_0$. Here the averaged restoring response is isotropic: there is no preferred mean-field direction. Both systems resist transverse deformation through magnetic tension, rather than through compressive sound.

The propagating branch is not the entire perturbation space. For example, $\mathbf{v}=0$, $m_{ij}=f(\mathbf{r}){\delta }_{ij}$ and $\pi =f(\mathbf{r})$ solve the linearized equations without propagation, for any small, slowly varying $f$. Such a static stress perturbation in magnetoelastic wave equations balances a pressure gradient. Thus the requested wave conclusion applies to the oscillatory branch; it must not exclude stationary initial stress perturbations.

The denominator of the travelling response vanishes atFor the positive speeds in the paper, a finite travelling response exists whenever $c\ne c^{\ast }$. At $c=c^{\ast }$, the dynamic condition would reduce to $p_1/\rho =0$, which is impossible. The finite-depth dispersion relation of a surface gravity wave with wavenumber $k$ is ${\omega }^2=gk\tanh (kh)$. At $k=1$, its phase speed $\omega /k$ equals $c^{\ast }$: the moving pressure pattern drives a natural surface gravity wave at its own speed, producing resonance.

The restriction concerns a bounded travelling response. To distinguish this from existence of a solution of an initial value problem, set $T=\tanh h$, $P=p_1/\rho$, and ${\omega }_{\ast }^2=gT$. Ifthe linearized free-surface boundary conditions reduce to the forced harmonic oscillatorFor example, if the forcing is switched on at $t=0$ from a flat surface with zero initial velocity, then, for $c\ne {\omega }_{\ast }$,This satisfies $H(0)=\dot{H}(0)=0$ and contains both free oscillations and the forced component. At $c={\omega }_{\ast }$ a solution of the initial value problem still exists, but its particular solution containswhich grows linearly with time. Thus the ideal inviscid flow model predicts resonant growth rather than a finite constant amplitude. Additional free oscillations can set the initial conditions but cannot remove that secular growth.

For a plane wave with amplitude $\alpha$, wavenumber $k>0$ and angular frequency $\omega >0$, the dispersion relation $\omega =\omega (k)$ specifies which frequencies the governing equations permit. A point of fixed wave phase obeys $kx-\omega t=\text{constant}$, so the phase speed isIt is the speed of individual crests. For a narrow wave packet, two nearby plane waves with differences $\Delta k$ and $\Delta \omega$ have a slowly varying wave envelope depending on $\Delta kx-\Delta \omega t$. Taking the limit gives the group velocityThis is the envelope's propagation speed and, in the present linear lossless wave model, the energy transport speed for a narrow packet.

A dispersive wave has phase speed dependent on wavenumber, so different components can change their relative phases and alter the packet's shape. A nondispersive wave on the chosen propagation branch has constant phase speed over the relevant band:Indeed $d{c}_p/dk=(c_g-c_p)/k$, so equality throughout a band is equivalent to constant phase speed there. Equality at an isolated wavenumber does not make the entire wave family nondispersive.

Define the dimensionless parameterThe stipulated hierarchy $kh\ll \delta \ll 1$ implies both $kh\ll 1$ andThus the finite-depth correction uses $\tanh (kh)\sim kh$, while the restoring effect of bending stiffness is small compared with gravity. The dispersion relation becomesHence, to leading order, these surface gravity waves are nondispersive waves withFor clarity, the first corrections to the phase speed areThey exhibit the weak residual wave dispersion at finite wavenumber. The claimed nondispersive behavior is the asymptotic leading-order result, not exact equality for the complete relation.

In the deep-water limit $kh\gg 1$, $\tanh (kh)\sim 1$. Put $\beta =B/\rho >0$. The squared phase speed of the flexural-gravity wave is thenIt diverges as $k\to 0^{+}$ and as $k\to \infty$ in this limiting model. Its unique stationary point satisfiesThe second derivative $2g/k^3+6\beta k$ is positive. Using $g/k_{\ast }=3\beta k_{\ast }^3$ gives the minimum phase speed of a flexural-gravity wave:This is the deep-water minimum: applying it to a finite depth requires $k_{\ast }h\gg 1$. Otherwise the $\tanh (kh)$ factor must be retained in the minimization. Within the deep-water approximation, phase speed and group velocity coincide at $k_{\ast }$, although the wave family remains dispersive.

A fixed wave phase satisfies $kx-\omega t=\text{constant}$, so the signed phase velocity is $v_p=\omega /k$. The group velocity $v_g=d\omega /dk$ gives the leading translation of the wave envelope of a narrow-band wave packet on a differentiable dispersion relation branch.

Here $v_p=k/m$. For the conventional positive branch $m>0$, $\omega >0$ and $k>0$,On the opposite propagation branch $k<0$ it is $-\sqrt{\omega /m}$, so $\sqrt{\omega /m}$ is the phase speed for either sign of $k$. More generally, for real $m\ne 0$ use $v_p=k/m$ with $k^2=m\omega$, allowing the sign of $\omega$ to follow that of $m$. The $k=0$ constant solution has no distinguished phase velocity.

The hydrostatic buoyancy anomaly is $b^{\prime }=f_0{\psi }_z$, and its leading adiabatic evolution is $D_g{b}^{\prime }/Dt+N_0^{2}w=0$. The small vertical velocity advects the background stratification; vertical advection of the perturbation buoyancy is higher order. ConsequentlyThe impermeability condition at a stationary lower boundary $z=\alpha y$ is $w=\alpha v$. Using $v_g={\psi }_x$ gives the topographic buoyancy boundary condition for quasi-geostrophic flowFor application at the reference plane, the actual small-height condition is $|\alpha |L/D\ll 1$. Slopes scaled by $D/L$ are of Rossby-number order, so this topographic contribution is retained alongside the small ageostrophic vertical velocity. The dimensional slope should not be compared to the Rossby number without also specifying the aspect ratio.

The basic quasi-geostrophic streamfunction is $\overline{\psi }=-\Lambda zy$. It gives $U=\Lambda z$, basic buoyancy gradient ${\overline{b}}_y=-f_0\Lambda$, and spatially constant interior quasi-geostrophic potential vorticity. Linearizing its conservation law therefore givesThe printed equation omits the shear-advection term. This cannot be omitted for general disturbances: an interior anomaly proportional to $\cos [k(x-\Lambda zt)]$ is materially conserved but not independent of time at a fixed location. The intended zero-interior-PV Eady model modes do have $q^{\prime }=0$ and satisfy both expressions; the ensuing mode calculation is valid in that subspace.

At either boundary, linearizing the geostrophic material derivative includes $v^{\prime }{\overline{\psi }}_{zy}=-\Lambda {\psi }_x^{\prime }$. LetThe two boundary equations areThese are exactly the respective bottom and top slope corrections to the rigid-boundary buoyancy condition for quasi-geostrophic waves.

Assume periodicity horizontally and first consider $q^{\prime }=0$. Multiply its defining equation by ${\psi }_x^{\prime }$ and integrate over the layer. The two horizontal Laplacian terms vanish by integration by parts; the remaining volume term after one vertical integration is an $x$-derivative. ThereforeFor nonzero $A_0,A_D$, the boundary equations giveThe top advection term integrates to zero. Subtracting proves conservation of the boundary pseudomomentum of a sloping Eady layer:The original PDF uses squared vertical derivatives and the upper slope $\gamma$ in the first denominator; these details are corrupted in the converted TeX. For arbitrary nonzero interior PV the same calculation instead gives $d\mathcal{P}/dt=(f_{0}s{)}^{-1}\int {\psi }_x^{\prime }{q}^{\prime }dV$, which need not vanish. For example, at an instant take ${\psi }^{\prime }=\cos kx+(z/D{)}^2\sin kx$, independent of $y$; the zonal average of ${\psi }_x^{\prime }{q}^{\prime }$ is $-ks/D^2$. Thus the boundary-only conservation statement also needs the zero-interior-PV qualification.

For an exponentially growing smooth normal mode, this qualification is automatic: the interior equation is $ik(\Lambda z-c)\hat{q}=0$, and a phase speed with nonzero imaginary part cannot equal the real basic velocity. Therefore $\hat{q}=0$. If $A_0{A}_D<0$, the conserved boundary quadratic form is definite, so a nonzero exponentially growing boundary amplitude is impossible. Instability consequently requiresfor $\Lambda \ne 0$. When one $A$ vanishes, the corresponding boundary derivative is simply advected and must vanish for a growing mode. The other boundary then yields a real single-edge-wave speed; if both vanish, the zero-interior-PV Neumann problem has no nonzero growing mode. These degenerate cases are neutral and should not be handled by dividing by zero.

For equal slopes, put $\zeta =z/D$, $C=c/(\Lambda D)$, $M=N_0|k|D/|f_0|$, and $a=N_0^2\alpha /(f_0\Lambda )-1$. The normal-mode interior equation is ${\Phi }_{\zeta \zeta }-M^2\Phi =0$ for the zero-interior-PV modes. ThusThe coefficient equations areSetting their determinant to zero, without dividing by $C$ or $a$, givesso the Eady model with parallel sloping boundaries hasFor positive $k,f_0$, $M$ is the printed $\mu$, and $a$ is the printed $\tilde{\alpha }$. The expression is even in the sign of that vertical scale. For $a\ge 0$ the radicand is at least $1/4$, so both speeds are real and these modes have no exponential growth. For equal slopes the earlier necessary condition is merely $a^2>0$, which is satisfied in many stable cases and is not sufficient. At $a=0$ the speeds are $0$ and $\Lambda D$; the invariant with nonzero denominators cannot be used, but the determinant still correctly gives neutral modes. No claim excluding transient amplification is needed.

Let the equilibrium ideal gas be at rest, with uniform ${\rho }_0,p_0$ and ${\mathbf{B}}_0=B_0{\mathbf{e}}_z$. This local wave calculation neglects gravity, viscosity and magnetic diffusion. Define the adiabatic sound speed $c\equiv v_s=(\gamma p_0/{\rho }_0{)}^{1/2}$, the Alfvén speed $v_a=B_0/({\mu }_0{\rho }_0{)}^{1/2}$, and a wavevector $\mathbf{k}=k(\sin \theta ,0,\cos \theta )$, where $k=|\mathbf{k}|$ and $\theta$ is its angle to the equilibrium magnetic field. Perturbations have dependence $e^{i\mathbf{k}\cdot \mathbf{x}-i\omega t}$ and phase speed $v=\omega /k$.

The linearized ideal magnetohydrodynamic equations for mass conservation, pressure, induction and momentum giveThe relation between pressure and mass density follows from eliminating $\mathbf{k}\cdot {\mathbf{u}}^{\prime }$ between the two material evolution equations; it applies to the nonzero-frequency wave sector. The induction result automatically has $\mathbf{k}\cdot \mathbf{b}=0$.

The $y$ MHD wave polarization is perpendicular to both the field and wavevector. It has $\delta p=\delta \rho =0$ and obeysThus the transverse Alfvén wave hasIn the $x$-$z$ plane the two remaining velocity components instead obey the symmetric matrix eigenproblemTaking its determinant gives $v^4-(c^2+v_a^2)v^2+c^2{v}_a^2{\cos }^2\theta =0$. Solving this quadratic equation in $v^2$ derives both magnetoacoustic wave phase speeds,The plus sign is the fast magnetosonic wave, the minus sign the slow magnetosonic wave. The radicand is nonnegative, since it equals $[(c^2-v_a^2{)}^2+4c^2{v}_a^2{\sin }^2\theta ]/4$. Both squared speeds are nonnegative. Each nonzero mode has forward and backward propagation; an advected entropy disturbance is a separate zero-frequency mode in this static equilibrium.

For the strong-collision ion-neutral Alfvén modes, put $a=k_{\parallel }{v}_A$ and assume $0<|a|/\mu \ll 1$. Rewrite the cubic asThe two slow roots obey $\omega =O(a)$, so the first approximation is ${\omega }_0^2=a^2/2$. Set $\omega ={\omega }_0+\delta \omega$ and retain the next order, of size $a^3$:ConsequentlyThese are weakly damped Alfvén waves. Their phase speed is $v_A/\sqrt{2}=B_0/\sqrt{4\pi ({\rho }_i+{\rho }_n)}$, since collisions make the neutral mass participate in the oscillation. Their ambipolar damping rate is ${\Gamma }_A=a^2/(8\mu )$, with ${\Gamma }_A/|\mathrm{Re}{\omega }_{\pm }|\ll 1$.

For the remaining root, at $a=0$ the cubic is ${\omega }^2(\omega +2i\mu )$, so start from ${\omega }_{d0}=-2i\mu$. Denote the cubic by $F(\omega ,a)$. Its value and derivative at this root areThe first correction is therefore $\delta {\omega }_d=i{a}^2/(4\mu )$, givingThis branch is purely damped. Indeed, the real growth-rate polynomial for $s=-i\omega$ is $s^3+2\mu s^2+a^{2}s+\mu a^2$, and its simple real root near $-2\mu$ stays real for sufficiently small real $a^2$. It represents rapid decay of ion-neutral relative motion. In the $a\to 0$ limit, the other two frequencies tend to zero; all three branches are accounted for.

Write $\mathbf{u}={\partial }_t\mathbf{\xi }$. In the uniform equilibrium, the linearized ideal magnetohydrodynamic equations reduce toIntegrating from the undisplaced reference state, and using the uniform background magnetic field, givesThese are the perturbations induced by the fluid displacement; independent time-independent changes to the reference state are excluded. The linear Lorentz force density is $(\nabla \times \delta \mathbf{B})\times \mathbf{B}/{\mu }_0$. Its magnetic pressure and magnetic tension parts giveFor a Fourier mode, put $s=\mathbf{k}\cdot \mathbf{\xi }$, $q=\mathbf{B}\cdot \mathbf{\xi }$ and $K=\mathbf{k}\cdot \mathbf{B}$. Then $\delta p=-i\gamma ps$ and $\delta \mathbf{B}=i(K\mathbf{\xi }-\mathbf{B}s)$. Substitution givesIf $s=q=0$, the fluid displacement is perpendicular to both the wave vector and the magnetic field. Only magnetic tension restores it, andThis is the Alfvén wave. The Alfvén velocity is the field-directed propagation vector: $\omega =\pm \mathbf{k}\cdot {\mathbf{v}}_a$. The signed phase velocity normal to the wavefront is $\pm v_a\cos \theta \hat{\mathbf{k}}$; the group velocity is $\pm {\mathbf{v}}_a$. The distinction matters for oblique propagation. For generic directions the Alfvén wave has one transverse polarization, with $\delta \rho =\delta p=0$.

To obtain the other magnetohydrodynamic waves, take the scalar products with $\mathbf{k}$ and $\mathbf{B}$:The determinant condition for a nonzero pair $(s,q)$ isThus the fast magnetosonic wave and slow magnetosonic wave haveTheir fluid displacements lie in the plane of $\mathbf{k}$ and $\mathbf{B}$ and are generally compressive. In the fast magnetosonic wave, gas and magnetic pressure provide the stronger restoring combination; in the slow magnetosonic wave, their perturbations oppose one another. The phase speeds depend on direction. They satisfy $v_{\mathrm{sl}}\le \min (c_s,v_a)$ and $v_{\mathrm{f}}\ge \max (c_s,v_a)$. For parallel propagation the two speeds are $c_s$ and $v_a$, with a degeneracy between a transverse branch and the Alfvén wave. For perpendicular propagation $v_{\mathrm{f}}^2=c_s^2+v_a^2$ and $v_{\mathrm{sl}}=0$; the latter is a nonpropagating limiting disturbance. These special directions require interpreting the polarizations by continuity rather than assuming three distinct nonzero frequencies.

Let $b$ denote buoyancy in a stably stratified Boussinesq approximation. Geostrophic balance and hydrostatic balance imply the thermal wind relation $f_0{U}_z=-{\bar{b}}_y$. Consequently the specified vertical shear requiresFor northern-hemisphere $f_0>0$ and $\Lambda >0$, the mean potential temperature decreases northward and increases upward: ${\bar{\Theta }}_y=-({\Theta }_0/g)f_0\Lambda$, ${\bar{\Theta }}_z=({\Theta }_0/g)N^2$. The corresponding mean isentropes slope upward toward the north, with $dz/dy=f_0\Lambda /N^2$. Stable potential-temperature stratification does not require the actual temperature to increase upward; the dry adiabatic lapse also enters that comparison.

In the rigid-lid constant-shear baroclinic model, perturbation quasi-geostrophic potential vorticity isThe mean interior potential-vorticity gradient vanishes on this constant-$f_0$, constant-$N$, linear-shear basic state. Linear advection of that perturbation therefore becomes $({\partial }_t+\Lambda z{\partial }_x)q^{\prime }=0$, the supplied interior equation. It retains slow balanced motion and omits friction, planetary-vorticity gradients and nonlinear disturbance advection.

Substitution of ${\psi }^{\prime }=\Psi (z)\cos (ly)e^{ik(x-ct)}$ givesA possible branch has zero interior perturbation potential vorticity, soBoth supplied applications have $f_0>0$, so this equals the source's expression. For an unstable mode the nonzero imaginary part of the phase speed prevents a real critical level $\Lambda z=c$; for real speeds the displayed solution remains a valid zero-interior-PV branch away from, and smoothly through, any critical level. Taking the real part gives the physical disturbance.