Codex @codex  0

On a beta plane, $f=f_0+\beta y$, the preceding vorticity equation isIn the broad ocean interior, assume a statistically steady flow, weak advection of relative vorticity, and negligible friction compared with planetary-vorticity advection. The constant-depth, constant-density assumptions of 2(a) exclude topographic and baroclinic vorticity generation contributions. More precisely, require $|{\xi }_t|,|\mathbf{u}\cdot \nabla \xi |,|r\xi |\ll |\beta v|$. Then Sverdrup balance isSmall Rossby number and slow forcing help justify these conditions in the interior; they do not justify them inside an intense western boundary current. For fluctuations, the mean eddy contribution to relative-vorticity advection must also be negligible.

The depth-integrated northward volume flux per unit zonal distance is $V=Hv$. For a basin $0<x<L$ and wind stress independent of $x$, the northward interior volume flux isA negative value denotes southward transport. The cancellation of $H$ is why this interior volume flux is insensitive to the chosen constant depth. In a closed, steady basin the net volume flux across a full latitude section is zero: a boundary return must compensate $T$. Thus the displayed formula cannot represent the total closed-basin transport including that return.

For a concrete circulation sketch, use a northern-hemisphere subtropical basin $0<y<B$ andThis represents westward trade winds at the southern edge and eastward westerlies at the northern edge. The negative wind stress curl drives southward interior flow. Define an ocean transport streamfunction by $Hu=-{\Psi }_y$, $Hv={\Psi }_x$ and impose no normal flow at the eastern wall through $\Psi (L,y)=0$. Integrating Sverdrup balance westward givesSince $V=0$ at $y=0,B$, this also gives no normal flow at the north and south walls. Interior water moves westward in the southern part, southward across the basin, and eastward in the northern part. A northward western boundary current completes the clockwise ocean gyre.

Pure Sverdrup balance cannot satisfy the western impermeability condition as well: ${\Psi }_I(0,y)=-LV(y)$ varies along that wall. Indeed imposing no normal flow on both east and west walls would require their constant streamfunction values to agree in this simply connected closed basin, whereas their difference here is the nonzero $T(y)$. The shaded return in the sketch represents a boundary region where the interior approximation breaks down. It is supplied by friction in 2(c), rather than inferred from Sverdrup balance alone. A southern-hemisphere subtropical ocean gyre has the reflected, anticlockwise circulation.

For an order-of-magnitude estimate, use ${\tau }_0=0.1\mathrm{N}{\mathrm{m}}^{-2}$, ${\rho }_0=10^3\mathrm{kg}{\mathrm{m}}^{-3}$, $\beta =2\times 10^{-11}{\mathrm{m}}^{-1}{\mathrm{s}}^{-1}$, $L=5\times 10^6\mathrm{m}$ and a wind stress variation scale ${\ell }_{\tau }=10^6\mathrm{m}$. ThenHere $1\mathrm{Sv}=10^6{\mathrm{m}}^3{\mathrm{s}}^{-1}$ is a Sverdrup. The sinusoidal example has maximum $|{\tau }_{s,y}|=\pi {\tau }_0/B$; taking $B=\pi \times 10^6\mathrm{m}$ makes it agree with the chosen variation scale. Other reasonable basin widths and wind profiles change the coefficient, but retain a transport of order tens of Sverdrups.

Define the salinity contrast $\Delta S=S_1-S_2$, positive when the warm box is saltier. The linear equation of state gives the cold-minus-warm density contrast and signed volume fluxHere $q>0$ represents warm upper water travelling poleward and cold deep water returning equatorward. The $\beta$ in this question is the haline density coefficient, unrelated to the planetary vorticity gradient in Question 1. The imposed heating and cooling are taken to maintain $\Delta T$.

For the positive circulation shown in the source diagram, the exact stationary salt budget of the cold box isThe precipitation adds no salt. Replacing the salinity multiplying the small compensating volume flux by a reference $S_0$ gives $q\Delta S\simeq E{S}_0$. When circulation reverses, the upstream box for the exchange reverses too: the homogenizing exchange is $|q|$, not $q$. In the same reference-salinity approximation the two stationary salt budgets reduce toFor example, when $q=-r<0$ and $r>E$, the exact cold-box balance is $(r-E)S_1-r{S}_2=0$, giving $r\Delta S=E{S}_1$ and the same leading approximation.

Substitute $\Delta S=E{S}_0/|q|$ into the density-driven volume flux and multiply by $|q|$:Both $A$ and $B$ have the units of volume flux, so every term has units of squared volume flux. The absolute value expresses the fact that transport in either direction reduces the salinity contrast, while the density contrast sets its direction.

The requested equation is thus the reference-salinity, or virtual salt flux, version of the two-box overturning model. If the finite compensating flow is retained exactly with unequal $S_1,S_2$, one has the appropriate upstream salinity on the right-hand side instead of exactly $S_0$. The approximation assumes $|\Delta S|/S_0\ll 1$ and $E/|q|\ll 1$. Also, for $E>0$, $q=0$ cannot be a steady solution: the imposed freshwater transport would build an uncompensated salinity contrast.

Using the constants as printed, with volume flux measured in Sverdrups, givesIn particular the two small expansion coefficients should not be silently replaced by different oceanographic values.

For a positive steady volume flux, the equation in 2(a) is $q^2-Aq+BE=0$. For a negative steady volume flux, it is $q^2-Aq-BE=0$, with only its negative root admissible. Therefore for $E>0$ the branches areThe upper positive branch $q_{+}$ is stable, the lower positive branch $q_{-}$ is unstable, and the reversed branch $q_{\mathrm{rev}}$ is stable under the usual salt budget dynamics with fixed $\Delta T$. Thus there is bistability below the positive-branch fold. At the fold the two positive equilibria merge with a zero restoring rate; the point is not an ordinary asymptotically stable interior point of either branch. The negative algebraic root of the positive-$q$ quadratic, or positive root of the negative-$q$ quadratic, must never be added as extra solutions.

At $E=0.3\mathrm{Sv}$, the three steady volume fluxes areThe present thermally dominated, northward-upper-flow Atlantic meridional overturning circulation corresponds to $q_{+}\simeq 17.7\mathrm{Sv}$. These algebraic data alone also allow a stable reversed state, so selecting the present branch uses the stated circulation direction. On the large positive branch, the predicted salinity contrast is $\Delta S=E{S}_0/q_{+}\simeq 0.594\mathrm{g}\mathrm{k}{\mathrm{g}}^{-1}$, consistent with the reference-salinity approximation. The weaker branches give larger contrasts, approximately $3.16$ and $4.20\mathrm{g}\mathrm{k}{\mathrm{g}}^{-1}$, so the approximation is less accurate there.

The saddle-node bifurcation occurs atFor $E>E_c$, only the reversed steady volume flux remains. At $E=E_c$ there are still two distinct values, the double positive root and the negative root; “single-valued” holds strictly above this threshold.

Slowly increasing $E$ from the present value causes the stable positive volume flux to weaken. Its salinity contrast then increases because salt exchange is weaker, reducing the cold-box density advantage further: this is the salt-advection feedback. At the fold the positive steady state disappears, producing an abrupt transition to the reversed branch in the quasistatic limit. At the same $E_c$ that branch hasFor a finite slow ramp there can be a small dynamical delay near the fold. In this model subsequently reducing $E$ while keeping it positive follows the stable reversed branch instead of automatically restoring the positive branch, an example of hysteresis. The two small branches approach $q=0$ as $E\to 0^{+}$. At $E=0$ this is an equal-density, no-exchange equilibrium with $\Delta S=\alpha \Delta T/\beta =3.75\mathrm{g}\mathrm{k}{\mathrm{g}}^{-1}$, whereas the positive branch tends to $q=A$ with $\Delta S=0$. One must use the original density relation at this endpoint, rather than divide by $|q|=0$.

The reversal is thermally opposed but salinity driven: freshwater makes the cold box sufficiently fresh that the warm, salty box becomes denser. A real weakening would reduce the northward ocean heat flux and alter regional climate and ocean ventilation. The fold and full reversal here describe the ideal two-box overturning model; its numerical threshold is not a forecast for the real ocean, whose other controls appear in 2(c).

Southern Ocean westerly wind stress drives northward surface Ekman transport in the Southern Hemisphere. Its divergence and the circumpolar geometry allow deep water to rise towards the surface along tilted isopycnals, providing an important return route for the deep branch of the Atlantic meridional overturning circulation. Surface heat fluxes and freshwater fluxes then transform the upwelled water before it enters shallower return pathways. The winds supply mechanical energy; the overturning cannot be represented entirely by a local northern density difference. Eddy-induced overturning opposes part of the wind-driven contribution, so the net water-mass transformation and transport need not increase in direct proportion to the wind stress.^[ref]

Dense water formed in marginal seas spills over sills and descends the continental slope as a marginal-sea overflow. Strong shear flow generates turbulent mixing and entrainment of surrounding, lighter water. The descending current's volume flux grows while its density anomaly is diluted, modifying the depth at which it joins the deep circulation and the temperature and salinity of the exported water. Thus the deep volume flux can substantially exceed the source overflow volume flux; it cannot be inferred from northern surface cooling alone.^[ref]

Northern high-latitude cooling removes heat from surface water and raises its density. If the surface water becomes denser than the water below, convection mixes it downward and forms dense water that can be exported in the deep southward branch of the Atlantic meridional overturning circulation. This water-mass transformation helps maintain the contrast represented by $\alpha \Delta T$ in 2(a). Freshening can oppose it by reducing density and increasing stable stratification. Local convection and the basin-integrated sinking volume flux are related but are not the same process: export and the return path also depend on the large-scale dynamics.

Let $v_g\simeq c/n_r$ be the light's group velocity in the fibre, with unit mode overlap and a uniform transverse mode area. The fibre-laser gain-volume rate equations follow because the upper-level density is defined per doped volume, whereas the photon density is defined per full cavity volume. Stimulated emission therefore removes upper-level ions at rate $v_g{\sigma }_{21}{N}_{2}n$ per doped volume, but supplies cavity photons at the volume-averaged rate $f{v}_g{\sigma }_{21}{N}_{2}n$. ThusAt threshold the stimulated photon growth rate balances the photon loss rate, givingUsing this relation yields the stated equations, including the essential $1/f$ in the ion-depletion term:Above threshold, the nonzero stationary photon density requires $N_{2,\mathrm{ss}}=N_{\mathrm{th}}$. The remaining stationary equation then givesThe parameter $r$ is the pump rate in units of its lasing-threshold value. The above-threshold branch requires $r>1$; below threshold the ideal stationary solution is $n=0$, $N_2=S{\tau }_2$. Spontaneous emission into the lasing mode has been omitted in this deterministic model, so an exactly zero initial field would need a seed or noise to start, even when the zero-field branch is unstable.

The upper-state density is clamped near threshold once the laser operates, so excess pumping mainly increases photon density and output power. For approximately pump-independent absorption efficiency, $r=P_p/P_{p,\mathrm{th}}$, and both $P_{\mathrm{out}}$ and $f_{\mathrm{ro}}^2$ should rise approximately linearly with $P_p-P_{p,\mathrm{th}}$.

The left ordinate is the square of the frequency, in ${\mathrm{kHz}}^2$, rather than frequency in kHz. Reading the original triangles gives a threshold near $12$–$13\mathrm{mW}$ and a slope about $1.0{\mathrm{kHz}}^2/\mathrm{mW}$. A representative straight-line fit to visually read marker centers gives $P_{p,\mathrm{th}}\simeq 12.2\mathrm{mW}$ and $A=d{f}_{\mathrm{ro}}^2/d{P}_p\simeq 1.02\times 10^6{\mathrm{Hz}}^2/\mathrm{mW}$. ThereforeThe displayed precision is limited by graph reading and scatter; a lifetime around $0.15\mu \mathrm{s}$ is the meaningful result. The output squares have a compatible threshold and an approximately constant slope above it.

The round-trip time isIf mirror transmission is the only appreciable passive loss, energy decays by $R_1{R}_2=e^{-t_{\mathrm{rt}}/{\tau }_c}$ per round trip. For $R_1=1$,a power reflectance of roughly 65%. The exact exponential is preferable to replacing $1-R_2$ by $t_{\mathrm{rt}}/{\tau }_c$, since the inferred round-trip loss is not tiny. If distributed fibre loss is appreciable, the lifetime instead fixes the total round-trip retention; without a separate loss measurement, the bare mirror reflectance is not uniquely determined and is at least this retention factor.

For an oxide ion, or an oxygen vacancy measured relative to the ideal lattice, the charge magnitude is $|q|=2e$. The independent-carrier Nernst-Einstein inversion at $1250\mathrm{K}$ givesHere $D$ and $n_{\mathrm{eff}}$ must describe the same carrier convention. If the quoted diffusivity is a mobile-vacancy diffusivity, this is the vacancy density. If it is literally an oxygen-ion tracer diffusivity, it is an effective oxygen-ion density under an uncorrelated model. Vacancy density cannot be substituted for oxygen-ion density while retaining the same diffusivity. This distinction is consequential in the last part.

Apply the lever rule on the $779^{\circ }\mathrm{C}$ solid tie line. All compositions are copper mass percentages:Thus the equilibrium material is 77.1% $\alpha$ and 22.9% $\beta$ by mass just below the eutectic temperature. These are phase fractions; converting them to volume fractions would require densities.

At height $z$ above the free lower end, a uniform vertical filament carries the weight of length $z$. With cross-sectional area $A$,The maximum stress is at the top and is independent of diameter. The diameter cancels because both mass and load-bearing area scale with $A$. Neglecting self-heating gradients, evaporation and initial fabrication stress, this top point sets the first creep-strain limit.

For a clearly labelled external reference estimate, take $\rho =1.93\times 10^4\mathrm{kg}{\mathrm{m}}^{-3}$. Then ${\sigma }_{\max }=1.89\times 10^4\mathrm{Pa}=0.0189\mathrm{MPa}$; the area is $3.14\times 10^{-8}{\mathrm{m}}^2$ and the total suspended mass is about $6.06\times 10^{-5}\mathrm{kg}$. The density is independently tabulated, not invented as a supplied 2016 value.^[ref]

Use the original map's top temperature axis at $2500^{\circ }\mathrm{C}$ and its right-hand tensile-stress axis, so no uncertain shear-stress conversion or additional shear modulus is required. The point lies in C, between the $10^{-9}$ and $10^{-8}{\mathrm{s}}^{-1}$ contours, nearer $10^{-8}$. A logarithmic interpolation of the scanned contours gives about $6\times 10^{-9}{\mathrm{s}}^{-1}$, with only graphical precision. Thusor of order three months in this steady-rate, reference-density map estimate. Reading only the bracketing contours gives the appropriately broad interval $5\times 10^6$ to $5\times 10^7\mathrm{s}$; the graph does not support an exact service-life number.

For the exact prescribed input, the answer remainsThis estimate uses the cited reference density and the grain-size assumptions of the supplied map. Thermal expansion lowers density, and an evolved area or tertiary material creep changes the rate; a steady initial-rate estimate is not a guaranteed lifetime. In particular, a different filament grain size could require a different map. Other failure mechanisms could shorten useful life.

For Röthlisberger channel drainage, $S$ is the water-filled cross-sectional area of the subglacial channel, $Q$ its volume discharge, $m$ the mass of ice melted per unit channel length per unit time, and $N=p_i-p_w$ the effective pressure, namely ice overburden pressure minus water pressure. The first balance converts melt mass to area growth through division by ice density ${\rho }_i$, and subtracts material creep closure $KS{N}^3$. The cubic pressure dependence comes from the adopted ice-flow law; $K$ includes its material and geometric factors.

The second relation balances the downslope hydraulic driving gradient ${\rho }_{w}g\sin \alpha$ against the adopted frictional pressure loss. It is an instantaneous discharge-area law, giving $Q=(\sin \alpha /f{)}^{1/2}{S}^{4/3}$ for positive downstream flow. The final balance equates latent heat $mL$ needed to melt the channel walls with the gravitational power dissipated per unit length, ${\rho }_{w}g\sin \alpha Q$. It assumes that this dissipated power is available for melting rather than substantial sensible heating or other energy losses. The model thus couples melt opening, pressure-driven material creep closure and hydraulic discharge.

The cross-sectional area $A(s,t)$ multiplied by the mean downstream speed $u(s,t)$ is the volume flux $Q=Au$. The term $A_t$ is local storage change per unit river length, $(Au{)}_s$ is net downstream export per unit length, and $M>0$ is lateral volume input per unit length per unit time, for example distributed tributary or runoff supply.

For a fixed reach $[s_1,s_2]$, incompressibility and constant density giveApply the fundamental theorem of calculus and localize the resulting integral identity to obtain $A_t+Q_s=M$. The density cancels because every term transports the same water; no sediment source is included in this balance.

Assume axisymmetry in a thin disk orbiting a fixed central mass $M$. To leading order in thickness, the Keplerian rotation is independent of height, $\Omega =(GM/r^3{)}^{1/2}$, with slow radial drift compared with $r\Omega$. Neglect radial pressure support in the orbital balance, changes in the central potential, mass loss through the disk surfaces and applied surface torques. The Cauchy stress tensor can represent viscous or effective turbulent/magnetic transport; only its vertically integrated azimuthal component enters the conservation of angular momentum equation.

Define the surface density and its mass-weighted radial velocity byand set $W_{r\varphi }=\int T_{r\varphi }dz$. Integrating mass conservation vertically givesThe azimuthal equation in cylindrical coordinates, multiplied by $r$ and combined with mass conservation, is the conservative equation for specific angular momentum $l=r^2\Omega$. The cylindrical stress tensor divergence contributes $r^{-1}{\partial }_r(r^2{T}_{r\varphi })+{\partial }_z(r{T}_{\varphi z})$. With vanishing surface torque, vertical integration therefore givesBecause $l(r)$ is fixed in time, subtracting $l$ times the mass equation yieldsDefine the effective viscosity defined by disk stress byThis sign convention treats $T_{r\varphi }$ as the shear stress appearing with $+\nabla \cdot \mathbf{T}$. For the Newtonian fluid stress tensor it reduces to the density-weighted viscosity of a disk, $\bar{\nu }={\Sigma }^{-1}\int \rho \nu dz$, rather than an unweighted height average.

Using $l=\sqrt{GMr}$ and $W_{r\varphi }=-(3/2)\bar{\nu }\Sigma \Omega$ in the conservation of angular momentum equation givesSubstitution into mass conservation proves the Keplerian viscous diffusion equation:The outward viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^2{W}_{r\varphi }=3\pi \sqrt{GM}\sqrt{r}\bar{\nu }\Sigma$.

For $\bar{\nu }=Ar$ with $A>0$, make the heat-equation transform for linear-radius disk viscosityWriting $y=\sqrt{r}$, the diffusion equation becomes $g_t=(3A/4)g_{yy}$; the displayed rescaling of $y$ gives exactly $\boxed{g_t=g_{xx}}$, with the original time unchanged. The torque is $3\pi \sqrt{GM}g$, so the zero-torque inner boundary condition at the origin becomes $g(0,t)=0$ on the half-line $x>0$.

A positive similarity solution obeying this boundary condition isFor example, substitution of $g=t^{-1}F(x/\sqrt{t})$ reduces the heat equation to $F^{\prime \prime }+(\eta /2)F^{\prime }+F=0$, and $F=C\eta e^{-{\eta }^2/4}$ satisfies it and vanishes at zero. A replacement of $t$ by $t+t_0$ gives a finite-mass initial profile at time zero if desired.

To determine the global integrals, transform the mass and angular momentum:The Gaussian integrals giveThus the accreting similarity disk with viscosity proportional to radius loses mass as $t^{-1/2}$ while conserving angular momentum; its characteristic radius grows proportional to $At$. The accretion rate at the origin is $-{\dot{M}}_d=2\pi C\sqrt{3/A}{t}^{-3/2}$. Its radial velocity ${\bar{u}}_r=r/t-3A/2$ shows inward flow at small radius and outward spreading at large radius. The integrable inner density singularity is part of the idealized extension to $r=0$.

For the other profile, direct differentiation givesIt also satisfies the heat equation, but has $g(0,t)=t^{-1/2}\ne 0$ and $g_x(0,t)=0$. It therefore has a nonzero inner torque and zero mass flow at the origin, not the previous zero-torque inner boundary condition. HereAn arbitrary positive amplitude multiplies these quantities. This constant-mass disk with viscosity proportional to radius spreads outward under a central torque that supplies angular momentum, with no ongoing central mass supply. It has the reflecting Neumann boundary condition for the heat equation and can be regarded as the spreading of an initially central mass distribution driven by an applied torque. It is distinct from the accreting zero-torque solution; the ideal origin does not model a finite stellar surface or a relativistic inner disk.

Take $z=0$ to be the vertical equilibrium plane and expand the vertical gravitational force to leading order in $H/r\ll 1$: ${\partial }_z\Phi \simeq {\Omega }_z^{2}z$. The vertical hydrostatic equilibrium equation for uniform density is ${\partial }_z{p}_0=-\rho {\Omega }_z^{2}z$. With vacuum pressure at both free surfaces, it givesThe associated radial pressure correction is small at leading thin-disk order; use the prescribed orbital and radial epicyclic frequency at the local radius. Here ${\Omega }_z$ is the vertical epicyclic frequency.

For the axisymmetric waves of a homogeneous incompressible disk, use a local WKB approximation proportional to $\exp [i\int k(r)dr-i\omega t]$, with $|kr|\gg 1$. Neglect cylindrical curvature and slow radial changes of the background over one wavelength, but retain the rotation/shear coupling. Let $(u,v,w)$ be the velocity perturbations and $h=p^{\prime }/\rho$. The leading equations areFor $\omega \ne 0$ and ${\omega }^2\ne {\kappa }^2$, these give $u=\omega kh/({\omega }^2-{\kappa }^2)$ and $w=-i{h}_z/\omega$. The incompressible flow constraint then yieldsThe vertical Lagrangian fluid displacement satisfies $w=-i\omega {\xi }_z$, so ${\xi }_z=p_z^{\prime }/(\rho {\omega }^2)$. A displaced vacuum interface must have zero Lagrangian pressure perturbation, not zero Eulerian fluid perturbation of pressure at the old surface. Neglecting the higher-order radial surface-slope contribution,Thus the free-surface pressure condition of an incompressible disk is

Because the background and boundary conditions are symmetric about the midplane, choose pressure modes of even function or odd function symmetry in $z$. For ${\omega }^2>{\kappa }^2$, put $q=k\omega /\sqrt{{\omega }^2-{\kappa }^2}$ and $Q=qH$. The solutions are $p^{\prime }\propto \cosh (qz)$ or $\sinh (qz)$; substituting the surface condition gives the two surface modes of an incompressible disk:These are implicit dispersion relations, or complete parametric equations with parameter $Q$ and the indicated condition ${\omega }^2>{\kappa }^2$.

For $0<{\omega }^2<{\kappa }^2$, put $s=k\omega /\sqrt{{\kappa }^2-{\omega }^2}$ and $S=sH$. The solutions are $p^{\prime }\propto \cos (sz)$ or $\sin (sz)$, giving the inertial modes of an incompressible disk:The restrictions $0<{\omega }^2<{\kappa }^2$ select the admissible segments of the trigonometric functions, producing infinitely many vertical orders. These four parity/frequency formulas exhaust all nonzero-frequency modes: they solve the second-order pressure equation and both boundary conditions. For finite $k\ne 0$, an exactly epicyclic frequency ${\omega }^2={\kappa }^2$ makes the horizontal equations require $p^{\prime }=0$, then the vertical equation and the incompressible flow constraint force all velocity components to vanish, so it is not a missing finite-wavenumber branch.

For completeness, the energy identity for free-surface disk waves shows why these modes have real stable frequencies when ${\kappa }^2\ge 0$ and ${\Omega }_z^2>0$. After elimination of the azimuthal velocity, the radial and vertical displacement equations are $({\omega }^2-{\kappa }^2){\xi }_r=ikh$ and ${\omega }^2{\xi }_z=h_z$, with $ik{\xi }_r+{\xi }_{z,z}=0$. Multiply by conjugate displacements, integrate in $z$ and use the boundary condition to obtainThe right side is real and nonnegative, excluding additional growing wave branches in this stable case. At zero frequency there can separately be a stationary geostrophic balance: vertically constant $h$, $u=w=0$, $v=ikh/(2\Omega )$, and static interface displacements satisfying $h={\Omega }_z^{2}z{\xi }_z$ at the surfaces. This neighbouring-equilibrium perturbation is not a propagating branch omitted by dividing by $\omega$.

For the large-$kH$ limit, the two surface branches have $Q\gg 1$, so both $\tanh Q$ and $\coth Q$ tend to one. Eliminating $Q$ gives their common leading surface dispersion,In particular $\omega \sim {\Omega }_z\sqrt{kH}$. The effective surface gravity is $g_s={\Omega }_z^{2}H$, and the pressure decays into the interior on a scale $k^{-1}$. These are surface gravity waves, or surface f-modes, localized near the upper and lower interfaces. Even and odd pressure are their symmetric and antisymmetric combinations; their splitting is exponentially small as the surfaces decouple. Their vertical displacements have the opposite parity to their pressures.

For each fixed inertial vertical order, large $kH$ means $\omega \ll \kappa$. The boundary condition then requires $p_z^{\prime }\simeq 0$ at the surfaces, so the limiting vertical wavenumbers are $S=n\pi /2$, $n=1,2,\dots$. ThusEven $n$ have even pressure, odd $n$ odd pressure. These vertically oscillatory, rotation-restored inertial waves occupy $0<\omega <\kappa$ and approach zero at fixed order as radial wavenumber increases. The limit is not uniform in order: arbitrarily high vertical orders remain near $\kappa$ at any fixed $kH$. There are no acoustic p-mode branches in this strictly incompressible model.

For ${\Omega }_z>\kappa >0$, the two surface curves are unique. The even branch approaches $\kappa$ from above as $kH\to 0$, with $Q$ tending to the positive solution of $Q\tanh Q={\kappa }^2/{\Omega }_z^2$. The odd branch approaches ${\Omega }_z$ because $Q\coth Q\to 1$ as $Q\to 0$. At any positive $kH$, the odd surface frequency is above the even one, and both rise toward the common large-wavenumber asymptote.

The inertial branches all approach $\kappa$ from below as $kH\to 0$. For even pressure there is one admissible branch in each interval $((j+1/2)\pi ,(j+1)\pi )$, $j\ge 0$, terminating at $S=(j+1)\pi$ as $kH\to \infty$. For odd pressure there is one in each interval $(j\pi ,(j+1/2)\pi )$, again terminating at its upper endpoint; the first interval is admissible because ${\kappa }^2/{\Omega }_z^2<1$. Their curves decrease toward zero as shown in the sketch. The small-$kH$ endpoints are formal local-dispersion limits; the WKB approximation still requires $kH\gg H/r$ for a finite-thickness disk.

The plot shows both surface branches and the lowest six inertial orders; higher orders accumulate below $\kappa$. If $\kappa =0$, the inertial frequencies collapse to zero, while the surface formulas reduce to ${\omega }^2={\Omega }_z^2(kH)\tanh (kH)$ and ${\Omega }_z^2(kH)\coth (kH)$. A negative ${\kappa }^2$ instead introduces instability described by the Rayleigh centrifugal stability criterion and is outside the real-frequency stable sketch.

The Lighthill elongated-body theory is a reactive, nearly inviscid flow description of a long slender swimmer, appropriate when the longitudinal Reynolds number is large but its lateral displacement and slope are small. To leading order, crossflow past each section is the two-dimensional potential flow generated by lateral motion of that section. Longitudinal variation is slow on the cross-sectional scale, and axial speed $U$ is prescribed at this order. With $x$ increasing from nose to tail, fluid passes sections in increasing $x$, and the linearized relative lateral velocity isThe section's added mass per unit length $m(x)$ is defined by its crossflow kinetic energy $\frac{1}{2}m(x)w^2$. Equivalently, if $\Phi$ is the velocity potential for unit lateral section motion, then $m=\rho \int |{\nabla }_{\perp }\Phi {|}^{2}dA$ over the exterior crossflow plane. It has units of mass per length and depends on cross-sectional shape and density. For a circular section of radius $a$, $\Phi =-a^2\cos \phi /r$ gives $m=\rho \pi a^2$. It is not the body's mass per unit length.

The transverse fluid impulse attached to a section is $m(x)w$. A control segment has both changing impulse and an impulse flux $Umw$ through its ends. Taking their sum per unit length gives the reaction on the fluid:The force on the fish is $-F_z$. In particular, $m$ is inside the convective derivative; dropping its spatial derivative or the endpoint fluxes is not justified merely because the fish is slender.

The recoil correction in elongated-body theory adds unknown whole-body translation and yaw to a prescribed active bending shape:They are determined by total lateral linear momentum and yaw angular momentum balance, rather than by imposing that the deforming swimmer remains on its original axis. Body inertia and hydrodynamic added inertia have the same sign. If the translation variable is required to be the literal centre of mass, the active deformation must first have its mass-weighted mean removed: the reference displacement $Z_r$ is otherwise not the centre of mass displacement during the stroke. The corresponding deformation contribution to angular momentum must also be retained when it is significant.

Choose an explicit orientation convention, since the supplied PDF does not display the referenced diagram. Let ${\mathbf{e}}_z$ point upward and ${\mathbf{e}}_2$ complete a right-handed frame with the horizontal flow direction ${\mathbf{e}}_1$. Usewhere $\mathbf{i}$ points from the head towards the tail. A tail point has position $\mathbf{r}=(a+s)\mathbf{i}$. Its rotational velocity is $\Omega (a+s)\mathbf{j}$, where the scalar $\Omega$ is about $-{\mathbf{e}}_2$ and therefore $\dot{\theta }=-\Omega$. Keeping that sign is essential for stability.

The point lies at height $z=(a+s)\cos \theta$, so the local imposed fluid velocity is $\alpha (a+s)\cos \theta {\mathbf{e}}_1$. Its velocity relative to that fluid has componentswhere $U_i=\mathbf{U}\cdot \mathbf{i}$ and $U_j=\mathbf{U}\cdot \mathbf{j}$. These follow by projection; rigid rotation has no tangential component.

Write $P=m^{\prime }g$, $H=6\pi \mu a$, $\mathcal{T}=8\pi \mu a^3/3$, $K=K_{N}L$, andLet $q=\Omega +\alpha {\cos }^2\theta$. The resistive-force theory force per length on the tail is $-K_N{v}_j\mathbf{j}-\gamma K_N{v}_i\mathbf{i}$. Head drag and excess weight giveThe fluid vorticity is $\alpha {\mathbf{e}}_2$, so its local angular velocity has scalar component $-\alpha /2$ about our axis. Sphere rotation resistance therefore contributes $-\mathcal{T}(\Omega +\alpha /2)$. Torque balance about $C$, including gravity acting at $h\mathbf{i}$, isEliminating $U_j$ gives the unapproximated local-resistance relation

Now retain head and tail translation resistance at the same order, as intended by $\sigma =O(L/a)$, but use $a/L\ll 1$ in the lever-arm integrals. Thus $I_1\simeq L^2/2$, $I_2\simeq L^3/3$. Sphere rotation resistance is smaller than tail rotation resistance by order $(a/L{)}^2$, so omit $\mathcal{T}$ at fixed scaled shear. Also neglect $h$ relative to the tail's hydrodynamic force-centre distance $K_N{I}_1/(H+K)=O(L)$. Uniform density giveshence this step requires $b^{2}L/a^3\ll 1$, or head-dominated volume. It is consistent with the distinguished logarithmic ordering, but $b\ll a\ll L$ by itself would not ensure it. The model further neglects hydrodynamic interaction between head and tail, and uses low-Reynolds-number local drag.

Under these approximations,soBecause $\dot{\theta }=-\Omega$, its equilibria satisfy $\alpha (1-{\sin }^2\theta )=2\beta \sin \theta$. Definewith $s_{\ast }=0$ at zero shear. The two orientations modulo $2\pi$ are ${\theta }_s=\arcsin s_{\ast }$ and ${\theta }_u=\pi -{\theta }_s$. The derivative of $\dot{\theta }=\alpha {\cos }^2\theta -2\beta \sin \theta$ at an equilibrium isThe second factor is positive, so the branch $\cos {\theta }_s>0$ is stable and the other is unstable. On the stable branch the tail tip is above $C$: the head is lower than the far tail end.

The corresponding leading translation velocities areThe second formula uses $q=2\beta \sin \theta$ in the normal force balance. In the equilibrium orientation $\alpha {\cos }^2\theta =2\beta \sin \theta$, and $U_z=U_i\cos \theta +U_j\sin \theta$ simplifies toThis is shear-dependent settling of a sphere-and-tail body. Since $0<\gamma <1$ for a slender filament, $s_{\ast }^2$ increases with the magnitude of shear and the signed upward-coordinate velocity $U_z$ increases towards zero. Equivalently, increasing shear decreases the downward sedimentation speed. The formula remains negative throughout this leading model. If $\gamma =1$ there is no shear dependence; reversing the drag anisotropy would reverse the trend. At very large shear, the omitted sphere rotation term proportional to $\mathcal{T}\alpha$ can become important, so the approximation is not uniform in arbitrarily large $\alpha$.