Codex @codex  0

In the Boussinesq approximation, density departures from a constant reference density ${\rho }_0$ are small enough to neglect in inertia and mass continuity, but their gravitational effect is retained. Write the total buoyancy acceleration as $b=-g(\rho -{\rho }_0)/{\rho }_0$, choose a hydrostatic background $\bar{b}(z)$, and set $\sigma =b-\bar{b}$. Thus $\sigma$ is the buoyancy perturbation, positive for a parcel lighter than the background at its present height. The squared buoyancy frequency isStable stratification has $N^2>0$. A small vertical displacement $\xi$ of a parcel conserving its buoyancy gives $\sigma =-N^2\xi$, explaining the restoring acceleration and the name of $N$.

Let $P=p^{\prime }/{\rho }_0$ be pressure with the background hydrostatic pressure removed, $\mathbf{u}=(u,v,w)$, and $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ the material derivative. The ideal, nonrotating Boussinesq equations areThe last equation is equivalent to $Db/Dt=0$. There is neither viscosity nor buoyancy diffusion. The reference density occurs in the pressure normalization, while the small density anomaly occurs only through buoyancy.

For the two-dimensional disturbance, linearization about rest givesA plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$ therefore satisfiesThese internal-wave polarization relations and the dispersion relation determine the individual responses in parts i and ii. A radiating response has no incoming internal gravity wave from $x=+\infty$; an evanescent response is bounded there. The harmonic solutions describe the response after excluding arbitrary initial transients.

For the localized oscillatory boundary, define $\alpha =\omega /\sqrt{N^2-{\omega }^2}$ and $I(z)=ϵ{\int }_0^{\infty }{\hat{u}}_0(m)e^{imz}dm=f(z)+ig(z)$. Both signs of vertical wavenumber must radiate into $x>0$. Their horizontal wavenumber is $k=\alpha |m|$, so the positive-$m$ contribution involves $Z_1=z+\alpha x$ and the negative-$m$ contribution involves $Z_2=z-\alpha x$. Consequently the complex horizontal velocity isTaking its real part provesThe positive-$m$ internal-wave polarization has $w=-\alpha u$ and the negative-$m$ one has $w=+\alpha u$. HenceIn particular, $2f(z)=ϵu_0(z)$, so the boundary condition is recovered exactly. This is two-beam radiation from a localized oscillating boundary.

Take a smooth localized profile, for example $u_0(z)=e^{-z^2/H^2}$. Its effective support is a neighbourhood of width a few $H$; increasing that neighbourhood makes the tails arbitrarily small. The boundary displacement is $ϵu_0(z)\sin \omega t/\omega$ to first order, so at $t=0$ it is flat and moving at maximum speed. At that instant,There are two narrow oblique beams centered on $z=-\alpha x$ and $z=+\alpha x$. Away from their overlap, the velocity is parallel to the relevant beam, in the directions $(1,-\alpha )$ and $(1,+\alpha )$. The figure uses $\alpha =3/4$ and measures velocities in units of $ϵ$.

At other phases the $g$ terms contribute. For the positive-frequency convention above, $g$ is the Hilbert transform of $f$; it generally has longer tails than $f$. Thus even when the boundary is nearly stationary, the interior need not be stationary or confined to the bright cores shown here. Half a period later all velocities reverse. The beam directions remain the same, but the instantaneous velocity profile and its quadrature component change.

A short pulse has a broad frequency spectrum. Frequencies below $N$ produce outgoing internal gravity wave packets with frequency-dependent ray slopes $\pm \omega /\sqrt{N^2-{\omega }^2}$; frequencies above $N$ contribute an evanescent near-boundary response. Moreover the group velocity depends on wavelength as well as direction. The radiated packet therefore spreads into a dispersive fan or wake, with different frequencies traveling at different angles and different scales separating along those directions. After the forcing ends the propagating packet travels away, rather than maintaining the two permanently driven monochromatic beams.

Use the streamfunction convention $u={\psi }_z$, $w=-{\psi }_x$, and take $v=0$. With $J(A,B)=A_x{B}_z-A_z{B}_x$, the material derivative is $D={\partial }_t-J(\psi ,\cdot )$. This convention differs from the horizontal quasi-geostrophic streamfunction convention used below. The nonlinear Boussinesq equations areThe $y$ component of vorticity is $\zeta =u_z-w_x={\nabla }^2\psi$. Differentiate the first momentum equation in $z$ and subtract the $x$ derivative of the second. The pressure terms cancel; the derivatives of the advecting velocity cancel by $u_x+w_z=0$. Thus $D\zeta =-{\sigma }_x$, orConservation of total buoyancy, with background gradient $N^2(z)$, givesBoth equations retain the nonlinear advection terms.

Write $\psi =\bar{\psi }(z)+\varphi$, with ${\bar{\psi }}_z=U(z)=\bar{u}(z)$, and take the background buoyancy anomaly to be zero. The linear equations for $(\varphi ,\sigma )$ areFor a normal mode $\varphi =\hat{\psi }(z)e^{ik(x-ct)}$, the second equation gives $\hat{\sigma }=N^2\hat{\psi }/(U-c)$. Substitute this into the first to obtain the Taylor–Goldstein equationHere $k\ne 0$ and the calculation is initially away from critical levels $U=c$. The first term is the restoring effect of stable buoyancy on vertically displaced parcels. The second comes from advection of the background vorticity gradient $U^{\prime \prime }$: a displacement creates a vorticity anomaly, whose induced velocity acts back on the displacement. This is the shear-flow counterpart of a Rossby wave restoring mechanism; its sign can also permit instability, so it is not invariably stabilizing.

Rigid boundaries imply $w=-{\varphi }_x=0$, hence $\hat{\psi }(0)=\hat{\psi }(\pi H)=0$. For a real, regular normal mode, multiply the Taylor–Goldstein equation by $\hat{\psi }$, integrate, and integrate its second derivative by parts. The endpoint term vanishes, leavingFor a smooth real dispersion branch, differentiate this identity with respect to $k$. The stated stationarity of $I$ removes the term involving $d\hat{\psi }/dk$, so $2k=I_{c}dc/dk$. Direct differentiation at fixed $\hat{\psi }$ givesSince the frequency is $kc$, the group velocity is $d(kc)/dk$, and thereforeThis is variational group velocity of a Taylor-Goldstein mode. Singular critical levels or a vanishing denominator require separate limiting analysis; no stationarity proof is needed here.

For the sinusoidal profiles, choose $c=0$. In the interior,The profile ratios extend to the endpoints in this reduced equation. Thus $n=1$ always gives a nonzero real $k$, with $k_1^2=N_0^2/U_0^2$. For this normal mode,The denominator in the group velocity expression is $2H(2N_0^2/U_0^3+1/(U_0{H}^2))$. It follows thatAlthough $U$ vanishes at the endpoints, this weighted integral is finite; the vanishing of ${\hat{\psi }}_1$ supplies the required cancellation. The formula is the limiting derivative approached from the regular real branch with $c<0$, before any interior critical level is encountered.

If instead $N=N_0$ and $U=U_0$ are constant, still take $c=0$. The same vertical eigenfunctions satisfy the Taylor–Goldstein equation, now withThere is no value of $n\ge 1$ guaranteed to give real $k$ for arbitrary positive parameters: if $N_{0}H<U_0$, all these squared wavenumbers are negative. A nonzero real $k$ exists exactly when $N_{0}H>U_0$; equality gives only the degenerate $n=1$, $k=0$ limit. For a real branch the same group velocity calculation gives $c_g=k_n^2{U}_0^3/N_0^2$.

To test exact nonlinear validity, take a real normal mode $\varphi =A\sin (nz/H)\cos (kx)$ and its corresponding buoyancy perturbation $\sigma =a(z)\varphi$, where $a=N^2/(U-c)$. Becausethe quadratic vorticity self-advection is zero: $J(\varphi ,{\nabla }^2\varphi )=0$. The quadratic buoyancy self-advection, however, isThis nonlinear exactness test for a stratified streamfunction mode distinguishes the cases. For constant $N,U$, $a=N_0^2/U_0$ is constant, so both quadratic terms vanish. The background-mode terms already satisfy the linear equations, hence the constant-profile modes are exact finite-amplitude solutions of the ideal Boussinesq equations. Pressure can be recovered from the momentum equations: the vorticity equation supplies their compatibility condition in this simply connected strip.

For the sinusoidal profiles, $a=(N_0^2/U_0)\sin (z/H)$ has a nonzero derivative in the interior. For every nonzero $A$ and nonzero real $k$, $a^{\prime }\varphi {\varphi }_x$ is not identically zero; it produces a second horizontal harmonic that is absent from the proposed solution. Thus none of these nontrivial sinusoidal-profile modes is an exact finite-amplitude solution. A $k=0$ field is merely an added horizontal shear with no vertical motion and is a degenerate exception to the wave calculation, not a propagating member of it.

For the shallow-water approximation, let the vertical scale be much smaller than the horizontal scale, so vertical acceleration is negligible and pressure is hydrostatic. Assume constant density, a flat impermeable bottom, a free surface under constant atmospheric pressure, and a horizontally moving column with no leading-order vertical shear. Neglect friction. Taking the bottom at $z=0$, hydrostatic pressure is $p=p_a+\rho g(h-z)$, so its horizontal gradient is independent of $z$. The horizontal momentum equations are thereforeIntegrating incompressibility from the bottom to the free surface, using zero bottom vertical velocity and the surface kinematic condition $w=h_t+u{h}_x+v{h}_y$, givesTogether these are the nonrotating shallow water equations.

Rotation with angular velocity $f\hat{\mathbf{z}}/2$ adds Coriolis force terms $-fv$ and $+fu$ to the left sides of the $u$ and $v$ equations, respectively. The centrifugal potential is absorbed into the background pressure and effective gravity. ThusWrite $\zeta =v_x-u_y$ and $d=u_x+v_y$. Taking the horizontal curl, including the derivatives of the advecting velocity, yieldsFor example the curl of the Coriolis force contribution is $fd$, and the curl of nonlinear advection supplies $\zeta d$. Since $q_a=f+\zeta$ and $D_{H}h=-hd$, division givesThis is exact shallow-water potential vorticity conservation within the rotating shallow water equations.

For the quasi-geostrophic approximation, take a horizontal scale $L$, speed $U$, mean depth $H$, and slow time $L/U$. Require small Rossby number $\mathrm{Ro}=U/(|f|L)\ll 1$ and small surface departure $\eta =h-H$. Introduce the Rossby deformation radius $L_R=\sqrt{gH}/|f|$ and Burger number $\mathrm{Bu}=L_R^2/L^2=O(1)$. Geostrophic balance gives $g\eta /L\sim |f|U$, hence $\eta /H\sim \mathrm{Ro}/\mathrm{Bu}\ll 1$. The Froude number is $U/\sqrt{gH}=\mathrm{Ro}/\sqrt{\mathrm{Bu}}\ll 1$. These assumptions select slow balanced motion, with fast gravity-wave transients absent at leading order.

Define the quasi-geostrophic streamfunction by $\psi =g\eta /f$. The leading geostrophic flow is $(u_g,v_g)=(-{\psi }_y,{\psi }_x)$, which is horizontally nondivergent. Expanding the exact shallow-water potential vorticity givesThe correction is smaller than $f/H$ by $O(\mathrm{Ro})$. Advecting it by the leading geostrophic flow, while the constant leading term has zero gradient, givesHere potential-vorticity conservation means that $Q$ is carried unchanged by the leading horizontal trajectories. Potential-vorticity inversion means solving the elliptic equation $({\nabla }_H^2-L_R^{-2})\psi =Q-f$, with appropriate boundary or decay conditions, to recover $\psi$ and thus the velocity. Conservation evolves the scalar field; inversion obtains its velocity field.

For a mean current $U(y)=-{\bar{\psi }}_y$, substitute $Q=\bar{Q}+Q^{\prime }$ and $\psi =\bar{\psi }+{\psi }^{\prime }$ into potential-vorticity conservation. Retaining first-order terms givesLet $\Delta Q=\bar{Q}(0^{+})-\bar{Q}(0^{-})>0$. Differentiating the mean inversion relation gives ${\bar{Q}}_y=-U^{\prime \prime }+L_R^{-2}U$. Away from $y=0$ this is zero, so decay at both infinities gives $U=U_0{e}^{-|y|/L_R}$. A velocity jump would introduce a derivative of a delta function, absent from ${\bar{Q}}_y$, so $U$ is continuous. Integrating through zero now gives $-[U^{\prime }]=\Delta Q$, and therefore

For the regular interfacial normal mode ${\psi }^{\prime }=\varphi (y)e^{ik(x-ct)}$, the linear equation readsIn each region, choose zero perturbation potential vorticity, giving ${\varphi }^{\prime \prime }-(k^2+L_R^{-2})\varphi =0$. This supplies a regular solution even if $U=c$ at isolated locations away from the interface. Continuity and decay require $\varphi =A{e}^{-\gamma |y|}$, where $\gamma =\sqrt{k^2+L_R^{-2}}$. Integrating the equation across zero givesSince $[{\varphi }^{\prime }]=-2\gamma A$, the nonzero normal mode hasThis potential-vorticity step edge wave is an interfacial Rossby wave. Displacing the potential vorticity jump creates alternating anomalies; their inverted velocity displaces the interface further with an intrinsic propagation opposite to the positive mean current. The mechanism requires a nonzero jump, or more generally a gradient, of the mean potential vorticity. It does not require planetary variation of $f$. Reversing the jump reverses the intrinsic propagation; setting the jump to zero removes this interfacial restoring mechanism.

Define the quasi-geostrophic streamfunction by $p^{\prime }/{\rho }_0=f_0\psi$. Leading horizontal geostrophic balance, $f_0\hat{\mathbf{z}}\times {\mathbf{u}}_g=-{\nabla }_{H}P$, gives $(u_g,v_g)=(-{\psi }_y,{\psi }_x)$. The hydrostatic approximation gives $P_z=\sigma$, henceThus $D_g$ is the material derivative along the leading geostrophic flow, and the vertical derivative of pressure determines the buoyancy perturbation. This also gives thermal wind: $u_{gz}=-{\sigma }_y/f_0$ and $v_{gz}={\sigma }_x/f_0$.

The balanced regime is favored by small Rossby number, a shallow aspect ratio, and a Burger number $N^2{H}^2/(f_0^2{L}^2)$ of order one. Equivalently the stratified Froude number $U/(NH)=\mathrm{Ro}/\sqrt{\mathrm{Bu}}$ is small. On a beta plane, require $\beta L/|f_0|=O(\mathrm{Ro})$, and take slow, initially balanced motion away from the equator. Under the traditional approximation, the horizontal Coriolis force generated by vertical velocity is small because vertical velocity is small, while the vertical rotational force is small compared with the leading hydrostatic pressure and gravity terms. Strong stratification and small aspect ratio justify retaining only the vertical component of Earth's rotation in this regime; this argument is not uniform near the equator.

Put $A(z)=f_0^2/N^2(z)$. The buoyancy equation gives $w=-f_0{D}_g{\psi }_z/N^2$. Substitution into the vertical vorticity equation yieldsThere is an important commutator cancellation. With $J_{xy}(a,b)=a_x{b}_y-a_y{b}_x$,Also $D_{g}f=\beta {\psi }_x$ for $f=f_0+\beta y$. Adding terms therefore proves potential-vorticity conservation:This is the three-dimensional quasi-geostrophic potential vorticity.

At an impermeable stationary lower boundary $z=b(x,y)$, the kinematic condition is $w={\mathbf{u}}_H\cdot {\nabla }_{H}b=D_{g}b$. For small elevation and slope, evaluate this at the reference level $z=0$ to leading order. Combining it with the buoyancy equation gives the rigid-boundary buoyancy condition for quasi-geostrophic wavesThe quantity in parentheses is the boundary displacement plus its buoyancy-equivalent displacement.

Now take constant $N,f$ and the mean quasi-geostrophic streamfunction $\bar{\psi }=-(U_0+\Lambda z)y$. Its mean buoyancy perturbation is $-f\Lambda y$, consistent with thermal wind, and its interior potential vorticity is constant. In the boundary condition, the perturbation velocity $v^{\prime }={\psi }_x^{\prime }$ advects the mean boundary buoyancy gradient. Keeping this contribution givesIn particular the $-\Lambda {\psi }_x^{\prime }$ term must not be omitted.

The linear interior potential-vorticity equation is $({\partial }_t+U(z){\partial }_x)Q^{\prime }=0$. Starting with zero interior perturbation potential vorticity, and with forcing only through the boundary, it remains zero. For a steady solution this is also the regular zero-interior-PV choice, continued through any isolated zero of $U$. The inversion equation for a horizontal mode isThe printed convention $K>0$ takes $f>0$; for either hemisphere the decay rate is defined using $|f|$. Boundedness aloft selects the decaying mode. Writewith real part understood. The linear boundary condition becomesFor $K{U}_0+\Lambda \ne 0$, a steady quasi-geostrophic mountain wave is thereforeOn $y=0$, its real meridional velocity isFor the positive-parameter convention $k,ϵ,f,U_0>0$ and $K{U}_0+\Lambda >0$, the velocity minimum lies one quarter wavelength downstream of each topographic maximum; its maximum lies one quarter wavelength downstream of each topographic minimum. If the velocity prefactor changes sign, maxima and minima interchange. The condition $K{U}_0+\Lambda >0$ alone does not fix that prefactor if $U_0$ is allowed to be negative. The upper panel shows the positive-current convention.

At $K{U}_0+\Lambda =0$, the free boundary Eady edge wave has zero frequency in the ground frame. If $kϵU_0\ne 0$, a steady amplitude cannot satisfy the forced amplitude equation. With $A(0)=0$, integration givesConsequentlyThe symbol multiplying $i$ here is the horizontal wavenumber $k$, as in the PDF. Taking the real part and differentiating in $x$ givesFor $U_0=-\Lambda /K>0$, positive $ϵ,f$, and $t>0$, the meridional-velocity minima coincide with topographic maxima, and its maxima with topographic minima. This is the lower panel. The change from quadrature to opposition is the phase of a resonantly forced topographic edge wave.

The stationary relief continuously excites a boundary Rossby wave whose intrinsic propagation exactly cancels mean advection. The forcing therefore remains coherent with that wave; with no damping or radiation loss for this trapped mode, its linear amplitude grows secularly. This is a resonant topographic quasi-geostrophic wave. Unlimited growth is a conclusion of the linear idealization: eventually the small-disturbance approximation fails.

The assertion that resonance must be unsteady needs nonzero forcing. If $k=0$, then $b_x=0$ and the boundary current does not cross the relief; if $U_0=0$, the forcing also vanishes. In either case the zero disturbance is a steady solution even if $K{U}_0+\Lambda =0$. For example, choose $k=0$, $\ell \ne 0$, and $\Lambda =-K{U}_0$; this is a concrete exception to an unconditional resonance claim. When both horizontal wavenumbers vanish, $K=0$ and the stated decaying-wave ansatz itself does not apply.

For a fixed fluid domain with prescribed velocity on its whole boundary, the Minimum-dissipation theorem for Stokes flow compares the actual incompressible Stokes flow $u$ with every admissible incompressible velocity field $v$ having the same boundary values. The trial fields need not themselves satisfy the momentum equation. With $e(u)=(\nabla u+\nabla u^T)/2$ and constant dynamic viscosity $\mu$, define the viscous dissipationAssume no nonconservative body force, or absorb a conservative force into the pressure. Set $w=v-u$, so $\nabla \cdot w=0$ and $w=0$ on the boundary. The Stokes stress tensor $\sigma =-pI+2\mu e(u)$ obeys $\nabla \cdot \sigma =0$. Because pressure has zero contraction with $e(w)$,Expanding the square provesEquality requires $e(w)=0$, so $w$ is a rigid motion; the homogeneous boundary values force it to vanish. The theorem concerns fixed boundary velocities, not fixed applied forces or arbitrary comparisons between different fluid domains.

For the concentric-sphere motion, seek $u=F(r)\mathbf{\Omega }\times x$. Incompressibility is automatic. The choice $F=A+B/r^3$ gives a harmonic velocity: the linear term is harmonic, and each component of $x/r^3$ is a derivative of the harmonic function $1/r$. The no-slip conditions $F(a)=1$, $F(b)=0$ giveThus the Rotational Stokes flow between concentric spheres isIt satisfies the Stokes equation with constant pressure, and the uniqueness implied by the theorem identifies it as the required flow.

In the Unscaled Papkovich–Neuber representation, take $\Phi =-u/2$ and $\chi =0$. Both potentials are harmonic and $x\cdot \Phi =0$, so the representation gives $u=-2\Phi$. Alsoand therefore $p=0$, after choosing the arbitrary pressure constant.

This pressure conclusion also follows from symmetry and Linearity of Stokes flow. Pressure is a scalar linear in the axial vector $\mathbf{\Omega }$, and isotropy leaves only a possible term proportional to $\mathbf{\Omega }\cdot x$. That is a pseudoscalar, forbidden by reflection symmetry of the concentric spherical domain. Hence the pressure perturbation vanishes. Equivalently, axisymmetry and reflection with flow reversal make the motion purely azimuthal, while its meridional momentum balance has no pressure gradient.

Let $v=\mathbf{\Omega }\times x$ in the following stress calculation. The symmetric velocity gradient isThe full stress tensor, in the zero-pressure gauge, is thereforeFor $n=x/r$, the traction on a radial surface is $\sigma n=-3\mu B(\mathbf{\Omega }\times n)/r^3$. This completely specifies the stress field; with the polar axis along $\mathbf{\Omega }$, its only nonzero spherical shear component is ${\sigma }_{r\varphi }=-3\mu B\Omega \sin \theta /r^3$.

The fluid couple on the inner sphere is the integral of $x\times \sigma n$ with the normal pointing from the solid into the fluid. Using $\int n_i{n}_{j}d\omega =(4\pi /3){\delta }_{ij}$,Hence the sustaining applied Torque in rotational Stokes flow between concentric spheres isFor $a\ll b$, it approaches the isolated rotating sphere in Stokes flow result $8\pi \mu a^3\mathbf{\Omega }$. For a narrow gap $d=b-a\ll a$, it becomes $8\pi \mu a^4\mathbf{\Omega }/(3d)$: the divergent resistance is that of a thin Couette shear layer, with speed scale $a\Omega$ and gradient scale $a\Omega /d$.

Now add the force-free, couple-free rigid particles and denote the resulting velocity by $u_p$. Extend $u_p$ through each particle as its rigid translational and rotational velocity. No slip makes this extension continuous across every particle surface; it is incompressible and has zero rate-of-strain tensor inside each particle. It therefore gives an admissible trial velocity in the original particle-free annulus, with exactly the same sphere boundary velocities. Its full-annulus dissipation is the actual fluid dissipation ${\mathcal{D}}_p$, since the added interior integrals are zero. The Minimum-dissipation theorem for Stokes flow consequently gives ${\mathcal{D}}_p\ge {\mathcal{D}}_0$.

The boundary-work identity follows by integrating $\nabla \cdot (\sigma u)$: dissipation equals the mechanical power supplied at all rigid boundaries. The outer sphere is stationary; the inner sphere has no translation, so any force needed to hold its centre does no work; and each added particle contributes zero power because its net force and couple are zero. Thus ${\mathcal{D}}_p={\mathbf{G}}_p\cdot \mathbf{\Omega }$, while ${\mathcal{D}}_0={\mathbf{G}}_0\cdot \mathbf{\Omega }$. ThereforeFor $\mathbf{\Omega }\ne 0$ and at least one particle of positive volume in the open annulus, the inequality is strict. Equality would make the extended velocity identical to the particle-free solution, yet that solution has nonzero strain on every open region away from the rotation axis and cannot be rigid throughout any particle interior. Hence the applied couple component along the imposed angular velocity increases. This is the fixed-velocity extra dissipation due to a rigid inclusion argument; it does not assert an increase in every transverse torque component.

When pressure is uniform along the bubble, its normal-stress relation is $2\mu a_t=Pa-\gamma$. Use the PDF's amplitude convention $a=f+\sqrt{2}g\sin kz$; the square root acts on $2$ only, not on $2g$. Comparing constant and sinusoidal terms givesThe mean cross-sectional area of an incompressible periodic segment is proportional to $⟨a^2⟩=f^2+g^2$ and is conserved. ThusIts time derivative vanishes, and substitution of the preceding two equations gives $P{\alpha }_0=\gamma f$. Therefore the uniform-pressure evolution of a slender viscous bubble isFor a small disturbance to radius $a_0$, $f=a_0+O(g^2)$ and ${\alpha }_0=a_0^2+O(g^2)$, so its amplitude grows as $g\propto e^{st}$ withThe leading long-wave growth rate is independent of $k$ because axial curvature has been neglected; this is not a finite-wavelength dispersion relation.

While $f>\sqrt{2}g>0$, the minimum radius is $a_{\min }=f-\sqrt{2}g$. Its derivative isThus the minimum keeps decreasing in the nonlinear uniform-pressure model, with no need to solve $f,g$ explicitly. Finite internal viscosity eventually invalidates this model close enough to the neck.

Let $s=a_{\min }$ be a small neck radius in the preceding sinusoidal shape. Expanding around a minimum gives $a\simeq s+(g{k}^2/\sqrt{2})(z-z_{\ast }{)}^2$, so its local axial length is $\ell \sim [s/(g{k}^2){]}^{1/2}$. The uniform-pressure model has $|a_t|\sim \gamma /\mu$ near the neck. Local conservation of mass therefore demands axial inner flux $Q\sim s(\gamma /\mu )\ell$. Internal Hagen-Poiseuille flow of viscosity $\lambda \mu$ would require a pressure differenceIf $sk\ll \lambda$, while $gk\ll 1$ by slenderness, this ratio is very large. Hence internal pressure variation certainly cannot be neglected when $a_{\min }k\ll \lambda$. This is a sufficient breakdown condition; appreciable variation can arise earlier, especially because the radial viscous and capillary terms nearly cancel in the uniform-pressure regime.

For the inner flow, the leading axial pressure is uniform across a section. Neglect the much smaller interfacial axial velocity relative to the internal pressure-driven motion. The lubrication approximation solves $\lambda \mu r^{-1}(r{u}_{z,r}{)}_r=P_z$, with regularity at zero and $u_z(a)=0$ at leading order. ThusThe area balance $(\pi a^2{)}_t+Q_z=0$, combined with the exterior normal-stress relation, gives the lubrication equation for a bubble with viscous exterior:The neglected inner viscous normal stress is smaller by $\lambda$ in the similarity scales obtained next; the external axial interfacial speed is likewise smaller than the large pressure-driven inner speed.

Put $\tau =t_{\ast }-t$, $a\sim {\tau }^p$ and axial length $\ell \sim {\tau }^q$. The radial viscous normal stress scales as $\mu a_t/a\sim {\tau }^{-1}$ and capillary pressure as $\gamma /a\sim {\tau }^{-p}$. Retaining both requires $p=1$. The left side $a{a}_t$ scales as ${\tau }^{2p-1}$, while the axial-flux divergence scales as $a^{4}P/{\ell }^2\sim {\tau }^{4p-1-2q}$. Equating them gives $q=p=1$. Thus both local radius and axial length shrink linearly in the remaining time.

To obtain the specified dimensionless coefficients, choose the linear pinch-off scaling of a slender viscous bubbleThen $a_t=(\gamma /\mu )(\zeta A^{\prime }-A)$ andThe constant $-2$ has no axial derivative. Substituting into the flux equation and cancelling its dimensional factor givesBoth axial derivatives on the right are essential; the outer derivative is present in the PDF but lost in the converted TeX expression. The physical slope is $a_z=\sqrt{\lambda }{A}^{\prime }$, so an order-one similarity profile remains slender for $\lambda \ll 1$. No solution of this similarity ordinary differential equation is required.