Codex @codex  0

The force scale is $F\sim \rho g{a}^3$. At the distant free surface, the Stokeslet's pressure and viscous normal-stress scales are $F/d^2$. Gravity restores a displaced surface by a pressure change $\rho gh$, soThus $h/d=O((a/d{)}^3)\ll 1$, consistent with treating the boundary as approximately flat and impermeable at leading order. This uses the stipulated gravity-dominated restoring balance.

The liquid-air interface has negligible tangential traction. The stress-free image of a normal Stokeslet places an opposite point force at the reflection point, a distance $2d$ above the bubble centre. Across the flat surface the combined normal velocity is odd and the tangential velocity is even, so the normal velocity and tangential stress vanish there. A same-radius inviscid drop of density $2\rho$ has downward excess weight $4\pi \rho g{a}^3/3$, exactly opposite the bubble's upward buoyancy. Its exterior Stokeslet therefore supplies this image, giving the asserted full-space comparison to leading order in $a/d$.

The image's incident axial velocity at the bubble isTo first reflection, the bubble translates with this locally uniform incident velocity plus its isolated rise speed $F/(4\pi \mu a)$. Hence the bubble rise near a distant stress-free free surface isOnly the first relative correction is retained; the image's spatial variation and additional reflections affect higher terms.

The customary capillary-controlled small-deformation condition is a small capillary number,equivalently a small Bond number $\mathrm{Bo}=\rho g{a}^2/\gamma$. This is a sufficient condition, not a necessary condition for the special isolated spherical solution in part (b). More precisely, the image's strain is $O(F/(\mu d^2))$, so its deformation-producing normal stress compared with $\gamma /a$ is $O(\mathrm{Bo}(a/d{)}^2)$. When exploiting the exact isolated solution, that smaller combination controls the leading boundary-induced distortion. If capillarity at the distant free surface is also retained, the gravity-dominated deflection scaling requires $\gamma /(\rho g{d}^2)\ll 1$; $a\ll \sqrt{\gamma /(\rho g)}\ll d$ allows both the usual spherical-bubble criterion and gravity-dominated surface restoration.

Let the undisturbed water surface be $z=0$, the oil's upper surface $z=\eta (x,t)$ and its lower surface $z=b(x,t)$, so $h=\eta -b$. Leading hydrostatic pressure in the nearly inviscid water is $p_w=-{\rho }_{w}gz$ relative to atmospheric pressure; in the oil it is $p_H=\rho g(\eta -z)$. Equality at the interface gives $\rho gh=-{\rho }_{w}gb$, henceThus the horizontal hydrostatic driving gradient isThis reduced gravity accounts for most of the oil thickness being below the undisturbed water level.

Because the oil is much more viscous than the water and air, tangential tractions on both faces are negligible. In a thin layer, the leading horizontal velocity is consequently independent of depth: the floating extensional viscous gravity current is a plug, rather than the parabolic profile of a no-slip gravity current. Integrating incompressibility $u_x+w_z=0$ between the moving material boundaries, and using their kinematic boundary condition, gives conservation of oil volume,To derive the factor $4$ in the extensional balance, incompressibility gives $w_z=-u_x$. The Newtonian fluid stress tensor has ${\sigma }_{zz}=-p-2\mu u_x$. Normal-stress balance relative to the hydrostatic part therefore requires $p=p_H-2\mu u_x$. The horizontal stress then satisfiesThus the longitudinal tensile force in an oil slice, relative to its hydrostatic pressure, is $4\mu h{u}_x$ per unit width. Balancing its difference across the slice with the horizontal driving force $-h{p}_{H,x}dx$ yieldsThe coefficient is the planar extensional Trouton ratio, arising from both the direct horizontal strain and the pressure correction required by the transverse normal-stress condition.

Using the pressure gradient above, integration in $x$ givesThe stated nose stress condition sets $C(t)=0$, not an arbitrary imposed inlet tension. ConsequentlyFor a material slice injected at time $t_0$, let $\tau =t-t_0$ be its age. The thickness equation along its Lagrangian trajectory is $dh/dt=-h{u}_x=-k{h}^2$, with $h=h_0$ at injection. DefineThen $h=h_0/(1+\tau /T)$. The slice injected during $d{t}_0$ has volume $h_0{u}_{0}d{t}_0$ per unit width, so at fixed observation time $t$,Since a newly injected slice has $x(t,t)=0$, integration and subsequent differentiation along a material slice give the constant-flux floating extensional gravity current:The oldest slice forms the nose. Measuring time from the start of injection, $t_0=0$ there, soIts speed ${\dot{x}}_N=u_0(1+t/T)$ is exactly the material velocity at the nose, as required.

For a fixed position $x$ behind the nose, the age needed to reach it is determined by $x/u_0=\tau +{\tau }^2/(2T)$, independently of the observation time. Hence both thickness and velocity become stationary there once the nose has passed. Eliminating $\tau$ givesTheir product is $h_0{u}_0$, so the steady profile carries precisely the inlet volume flux; integrating it to $x_N(t)$ gives the total injected volume $h_0{u}_{0}t$.

The lengthis the horizontal distance over which reduced-gravity driving causes an order-one change in inlet speed and thickness by viscous extension. The associated age is $T=8\ell /u_0$; the exact profile contains $x/(4\ell )$. It is an extensional adjustment length, not a shear-diffusion scale. The result applies while the stated thin-layer and negligible-inertia approximations remain valid.

Take positive $x$ in the snail's direction of travel, and let $y$ measure height above the rigid plane. The wave moves at laboratory velocity $U-V$. In the wave frame the plane's horizontal velocity is $V-U$ and the foot's is $V$. The foot is stationary as a shape, so its material velocity is tangent to $y=h(x)$: to leading lubrication theory order its components are $(V,V{h}_x)$. Apply the no-slip boundary condition to both surfaces.

The leading Stokes flow equations are $p_y=0$ and $\mu u_{yy}=p_x$. Their solution is the local Couette-Poiseuille flowThe normal-velocity condition on each steady surface and incompressibility make the wave-frame flux independent of $x$. Integrating the velocity givesPut $C=V-U/2$. Periodicity of pressure gives $0={\int }_0^L{p}_{x}dx=12\mu (C{I}_2-q{I}_3)$, and henceThe shear traction exerted by the fluid on the plane, positive in the $x$ direction, is

The snail is force-free horizontally: its gait consists of internal motion and it has no externally applied horizontal drive. Uniform locomotion therefore requires zero net horizontal fluid force on the foot per wavelength. The stresses on opposite periodic fluid faces cancel, so horizontal momentum balance also makes the net force on the plane zero. ThusSubstitution of $q$ proves the lubrication model of a crawling snail:The Cauchy-Schwarz inequality applied to $h^{-1/2}$ and $h^{-3/2}$ gives $0<\alpha \le 1$, with equality exactly for a uniform gap. Thus the denominator is positive, and a nonuniform waveform gives forward motion; the flat gap gives $U=0$.

For the sinusoidal shape, the given integrals yieldTherefore $\alpha =1/(1+A^2/2)$ andFor small amplitude, $U\sim 3A^{2}V$: the leading locomotion is quadratic, and there is no propulsion from a flat foot. As $A\to 1^{-}$, $U\to V$, so the wave becomes nearly stationary in the laboratory frame while the snail advances relative to it. However, the minimum gap tends to zero; finite speed in this limit does not imply finite dissipation, and $A=1$ itself is excluded by the noncontact model.

For the general waveform, use the mechanical work identity for viscous dissipation. Pressure does no boundary work because the boundary motion is tangential in the wave frame. The plane's velocity is constant, and its total force is zero, so it contributes zero net work. The upper boundary supplies, at leading lubrication order,The zero plane-force relation is $\mu U{I}_1=\frac{1}{2}\int h{p}_{x}dx$. ConsequentlyAs a direct check, integrating the squared shear rate gives $\Phi ={\int }_0^L[\mu U^2/h+h^3{p}_x^2/(12\mu )]dx=\mu I_1[U^2+12C^2(1-\alpha )]$, which reduces to the same result using the force-free speed relation. Zero horizontal force on the foot does not mean zero gait work: the tangential material velocity along its sloping surface includes vertical motion.

For prescribed speed $U>0$ and sinusoidal shape, eliminate $V=U(1+2A^2)/(3A^2)$ to obtainWriting $y=A^2$, the logarithmic derivative of the amplitude-dependent factor isThe only root in $0<y<1$ is $y=1/2$. The dissipation diverges at both endpoints, so this is the global minimum and gives the optimal sinusoidal waveform for viscous snail locomotion:The magically driven flat glider has simple Couette flow $u=Uy/h_0$ in the laboratory frame, and dissipation ${\Phi }_{\mathrm{glide}}=\mu L{U}^2/h_0$ per wavelength. HenceThe comparison is at the same $U$, mean thickness $h_0$ and wavelength $L$; the flat glider is externally propelled, whereas the snail must power its internal waveform.

With the streamfunction convention $(u,v)=({\psi }_y,-{\psi }_x)$, disturbance vorticity is $\omega =-\Delta \psi$. The basic Couette flow has constant vorticity $-1$. Linearizing the two-dimensional Euler equations for an inviscid fluid in vorticity form givesThe second term vanishes, soThis derivation uses incompressibility and an inviscid homogeneous fluid, without introducing a viscous decay term.

For the initial streamwise Fourier mode, write $\psi =e^{i\alpha x}a(y,t)$. Its transported Laplacian has amplitudeThe decaying Green function on the line isThis is the negative of the one-dimensional modified Helmholtz Green function for $-{\partial }_y^2+{\alpha }^2$: it is continuous and its first derivative with respect to $y$ has jump one at $y=s$. Consequently $({\partial }_y^2-{\alpha }^2)G=\delta (y-s)$, and inversion yieldsUnder the usual decay of the initial profile and its derivatives, the initial integral returns $\varphi$. There is no nonzero homogeneous solution of $a^{\prime \prime }-{\alpha }^{2}a=0$ decaying at both infinities.

For the large-time calculation, put $q=\alpha t$ and $H(s)=G(y,s)F(s)$, keeping the observation point $y$ fixed. Assume sufficient localization to justify the following integrations by parts; for example, $F$ may be a Schwartz function. The function $H$ is continuous, but its derivative in $s$ has jumpsince $G_s(y,y+)-G_s(y,y-)=1$. The Fourier decay from a Green-function derivative jump givesHere the regular second derivative is computed separately on either side of $s=y$. If it is integrable, its Fourier transform tends to zero by the Riemann-Lebesgue lemma. HenceThe cusp in the Green function, rather than a stationary point of the linear phase, supplies this leading contribution.

To obtain the streamwise velocity, differentiate the integral itself. This avoids an unjustified differentiation of a pointwise remainder. The kernel $G_y$ has a jump of $-1$ as a function of $s$ at $s=y$, so one integration by parts givesThe transverse velocity follows exactly from $v=-i\alpha \psi$, while the Laplacian equation gives vorticity without an asymptotic approximation:Thus the generic orders are $u=O(t^{-1})$, $v=O(t^{-2})$, and $\omega =O(1)$, at each fixed $y$. If $F(y)=0$, the leading coefficients vanish and faster decay can occur. For real physical disturbances, take the real part of these complex Fourier modes. The vorticity generally keeps oscillating rather than approaching a pointwise limit, even while the velocity decays. This is linear inviscid damping in unbounded Couette flow caused by phase mixing.

For arbitrary initial vorticity, the characteristics satisfy $\dot{x}=y$, $\dot{y}=0$, and $\dot{\omega }=0$. Therefore the exact solution isPutting ${\omega }_0=-F(y)e^{i\alpha x}$ recovers the preceding nondecaying mode. In general the shear is area preserving, so its transport preserves every finite Lp norm of vorticity. Local pointwise behavior for arbitrary ${\omega }_0$ depends on its spatial profile: an initially localized profile can move away from a fixed point. There is no universal pointwise nonzero limit or $t^{-2}$ decay of vorticity for arbitrary initial data.

The stated smoothness and decay of $\varphi$ alone do not guarantee the claimed algebraic rates. The intended integration by parts calculation requires control of derivatives and tails; convergence of the displayed integral alone is insufficient. Here is a counterexample to the literal unrestricted rate, with $\alpha =1$. Choose a nonnegative smooth function $r$ supported in $(-1/4,1/4)$ with $\int r=1$, and a sufficiently rapidly increasing sequence $t_n$. SetChoose $\sum e^n/t_n<\infty$. Then $F$ is smooth and integrable, $\varphi$ is smooth and tends to zero at both infinities, and every evolution integral converges absolutely. At the observation point $(x,y)=(0,0)$, all supports have $s>0$, soAt $t=t_n$, the $n$th summand is $-1/(2t_n)$. The Fourier transform of $r$ decays faster than any inverse power. Inductively choose $t_n$ so large that all earlier summands at $t_n$ together have modulus less than $1/(8t_n)$, and choose subsequent times so that ${\sum }_{m>n}1/t_m<1/(4t_n)$. The later summands then also have modulus less than $1/(8t_n)$. Thus $|\psi (0,0,t_n)|\ge 1/(4t_n)$, contradicting $O(t^{-2})$. This construction meets the printed assumptions; what it lacks is uniform derivative localization. The preceding boxed decay formulas are the correct intended result under the explicitly stated sufficient localization hypothesis.

Let $\chi =H(F)$ use the Heaviside step function, and ${\delta }_s=|\nabla F|\delta (F)$, so that ${\mathbf{n}}_s=\nabla F/|\nabla F|$ points from the body into the fluid. The kinematic boundary condition gives $F_t+\mathbf{u}\cdot \nabla F=0$ on the surface. In the exterior, mass conservation and the Euler equations for an inviscid fluid areProducts with $\chi$ are understood distributionally, using ${\partial }_t\chi =F_t\delta (F)$ and ${\partial }_i\chi =F_i\delta (F)$. With ${\rho }^{\prime }=\rho -{\rho }_0$, the extended continuity equation isThe corresponding extended momentum equation isbecause its advective surface term is $\rho u_i(F_t+\mathbf{u}\cdot \nabla F)\delta (F)=0$. Differentiate continuity in time and substitute the divergence of momentum. Subtracting $c_0^2{\nabla }^2({\rho }^{\prime }\chi )$ gives the Ffowcs Williams-Hawkings equationwhere the inviscid Lighthill stress tensor isIf $p$ denotes total pressure, this tensor includes the constant ambient pressure $p_0$. Its contribution $p_0{\nabla }^2\chi$ cancels the corresponding surface term exactly. For radiation integrals it is convenient to subtract $p_0$ from both terms, making $T_{ij}=\rho u_i{u}_j+(p^{\prime }-c_0^2{\rho }^{\prime }){\delta }_{ij}$ with $p^{\prime }=p-p_0$. The displayed source identity is unchanged by this consistent pressure-gauge subtraction.

Assume no incoming sound and no free initial acoustic field. Put $R_y=|\mathbf{x}-\mathbf{y}|$, ${\tau }_y=t-R_y/c_0$, $Q={\rho }_0{u}_n{\delta }_s$, and $P_i=p^{\prime }{n}_{s,i}{\delta }_s$. Convolution with the retarded acoustic Green function gives the three source integrals, for an observer outside the body:These are the surface thickness, surface loading and volume quadrupole contributions. The retarded surface Dirac delta distributions must remain inside the integrals: integrating them over a moving retarded surface produces the moving-surface retarded Jacobian, and is not generally the same as simply substituting one retarded time into an instantaneous surface integral.

For an acoustic compact-source approximation, suppose the dominant sources have size $\ell$, source frequency scale $\omega$, and $\omega \ell /c_0\ll 1$. Also require $r=|\mathbf{x}|\gg \ell$ and the acoustic far field $\omega r/c_0\gg 1$. Then $R_y\simeq r-\mathbf{n}\cdot \mathbf{y}$, $1/R_y\simeq 1/r$, and the leading source integrals use the common time $\tau =t-r/c_0$. Slow source motion relative to $c_0$ is implicit in this leading retardation approximation. DefineThus $m$ is net displaced mass per unit time, $\mathbf{f}$ is the fluid pressure force on the body, and $S_{ij}$ is the integrated excess Lighthill stress tensor. The minus sign in $\mathbf{f}$ fixes the requested loading convention; the force on the fluid is opposite. On the radiating $1/r$ terms, ${\partial }_{x_i}\simeq -(n_i/c_0){\partial }_t$, givingSpatial derivatives of $1/r$ produce near-field terms of orders $r^{-2}$ and $r^{-3}$, omitted here. Neglecting source delays is a separate approximation from discarding those observer near-field terms. If a leading integrated source vanishes by symmetry, its first nonzero delay correction can be as important as another retained multipole.

For the translating sphere, write its center as $\mathbf{X}(t)=a_1\cos (\omega t)\mathbf{e}$, with $\mathbf{e}$ a unit vector, and its velocity as $\mathbf{V}=\dot{\mathbf{X}}$. Let $\mathcal{V}=4\pi a^3/3$, $\eta =a\omega /c_0\ll 1$ and $d=a_1/a\ll 1$. Impermeability gives $u_n=\mathbf{V}\cdot {\mathbf{n}}_s$. Since its volume is constant,in the leading compact formula. The individual local thickness sources have scale ${\rho }_0(a_1/r){\eta }^2$, but their net monopole cancels.

At first order in amplitude and leading order in $\eta$, the exterior incompressible velocity potential is $\varphi =-a^3\mathbf{V}\cdot {\mathbf{r}}_s/(2r_s^3)$, where ${\mathbf{r}}_s$ is measured from the mean center. It has ${\partial }_{r_s}\varphi =\mathbf{V}\cdot {\mathbf{n}}_s$ at $r_s=a$. Linear Bernoulli equation gives $p^{\prime }=-{\rho }_0{\varphi }_t=({\rho }_{0}a/2)\dot{\mathbf{V}}\cdot {\mathbf{n}}_s$ there. Since $\int n_{s,i}{n}_{s,j}dS=(4\pi a^2/3){\delta }_{ij}$,This is the added mass of a sphere. The explicit leading loading contribution isIt scales as ${\rho }_0(a_1/r){\eta }^3$. In homentropic small-amplitude flow the volume stress is quadratic to leading order, $T_{ij}=O({\rho }_0{V}^2)$ over volume $O(a^3)$, so the volume acoustic quadrupole has scaleIt is smaller than the linear dipole by $O(d\eta )=O(V/c_0)$ and can contain a second harmonic.

For completeness, the exact thickness-source integral has a nonzero first delay correction even though $m=0$. Its first spatial moment isExpand its integrand at ${\tau }_y=\tau +\mathbf{n}\cdot \mathbf{y}/c_0$. The leading nonzero translating-sphere acoustic thickness dipole is thereforeIt too scales as ${\rho }_0(a_1/r){\eta }^3$. Thus zero net-volume monopole does not mean zero thickness radiation. The full leading linear sound, summing thickness and loading, isThe leading compact equation gives the requested zero $\dot{m}$ and added-mass $\dot{f}$ terms; this extra delay calculation is necessary to obtain the complete linear field when its nominal monopole has cancelled.

Suppress the common time factor $e^{i\omega t}$ and invert the stated Fourier transform with $(2\pi {)}^{-1}\int \Phi (k,y)e^{-ikx}dk$. Use limiting absorption principle $k_0\mapsto k_0-i0$, and choose $\gamma =\sqrt{k^2-k_0^2}$ with positive real part for the evanescent spectrum. For $|k|<k_0$ its outgoing continuation is $\gamma =+i\sqrt{k_0^2-k^2}$. The pole $-k_0$ lies above the inversion contour and $+k_0$ below it. These prescriptions fix both radiation and pole signs.

The even interior scattered transform solves ${\Phi }_{yy}-{\gamma }^2\Phi =0$, so it is $A(k)\cosh (\gamma y)$. The exterior transform is $B(k)e^{-\gamma (|y|-b)}$. Let $D^{-}(k)$ be the transform of their common upward normal derivative on $y=b$. On $x>0$ the rigid plates set both derivatives to zero; on $x<0$ continuity requires them to agree. Hence this derivative is supported on the negative half-line and is analytic below the transform contour. It follows thatThe difference between exterior and interior scattered traces on $y=b$ is $-D^{-}/L$, where $L=\gamma \sinh (\gamma b)e^{-\gamma b}$. At the open negative half-line the total traces must match, so this difference equals the incident wave $e^{i{k}_{0}x}$. Its negative-half transform is $-i/(k+k_0)$. The positive-half trace difference has an unknown upper-half-plane transform $C^{+}$, giving the Wiener-Hopf equationFactor $L=L^{+}{L}^{-}$, divide by $L^{-}$, and subtract the forcing pole using $C_0=L^{+}(-k_0)$:The left side is analytic below and the right side above; the pole in the last quotient is removable. Analytic continuation makes their common value entire. The finite-energy edge condition and algebraic factor growth make this entire remainder vanish: the allowed normal-velocity edge singularity is at most inverse-square-root, so its half-line transform is $O(|k{|}^{-1/2})$, while $L^{\pm }=O(|k{|}^{1/2})$. The same edge condition excludes a growing polynomial remainder. This pole subtraction in a Wiener-Hopf equation yieldsThe corresponding exterior result for the Wiener-Hopf solution for an open parallel-plate waveguide isBoth forms are useful: the first makes the common normal derivative transparent, and the second removes the apparent factor singularity at a grazing saddle.

The final unheaded request has a sign error in its printed exponent. To see this from the derived solution rather than an outgoing-decay assertion alone, simplify the interior transform toNear $k=-k_0$, $C_0/L^{+}(k)=1+O(k+k_0)$ and $e^{-\gamma b}\cosh (\gamma y)=1-b\gamma +O({\gamma }^2)$ for fixed $|y|<b$. ThusFor $x<0$, inversion of the first term gives $-e^{i{k}_{0}x}$, cancelling the continued incident potential. The remaining branch-point term uses $\gamma \sim i\sqrt{2k_0}(k+k_0{)}^{1/2}$ on the specified contour, so it is $b\sqrt{2k_0}(k+k_0{)}^{-1/2}$. Its inverse transform at $x=-r$ is proportional to $e^{-i{k}_{0}r+i\pi /4}/\sqrt{\pi r}$. Therefore the upstream decay of an open parallel-plate waveguide isThe amplitude decays as $(-x{)}^{-1/2}$, not $(-x{)}^{1/2}$. The positive exponent is present in the original PDF, so it cannot be proved under the stated outgoing problem. The incident pole cancellation and the explicit branch expansion give the corrected conclusion.

Fix time dependence $e^{-i\omega t}$ and use pressure amplitude for $\psi$. Write the incident acoustic plane wave as ${\psi }_i=A{e}^{i\alpha x-i\beta z}$, where $\alpha =k\cos \theta$, $\beta =k\sin \theta >0$. “Perfectly reflecting” specifies zero absorbed energy but does not uniquely select a boundary condition: a rigid wall has a Neumann boundary condition, whereas a pressure-release surface has a Dirichlet boundary condition. We first give the rigid-wall result, and then the alternative pressure-release result.

For the rigid surface the normal-velocity condition, without an immaterial normalizing factor, isExpand $\psi ={\psi }^{[0]}+{\psi }^{[1]}+{\psi }^{[2]}+\cdots$ in the height amplitude, scaling a fixed differentiable profile. The flat total field is ${\psi }^{[0]}=A{e}^{i\alpha x}(e^{-i\beta z}+e^{i\beta z})$. Taylor expansion givesLet $\hat{h}(q)=\int h(x)e^{-iqx}dx$ and $\beta (q)=\sqrt{k^2-q^2}$ on the outgoing sheet: it is positive real for propagating modes and positive imaginary for evanescent modes. Fourier transforming this boundary derivative gives $2A(k^2-\alpha q)\hat{h}(q-\alpha )$. The outgoing angular spectrum therefore yields the small-height Neumann scattering fieldFor an infinite stationary surface the Fourier quantities may be interpreted as random distributions or through a finite illumination window and its limiting mean. The correction is linear in $h$ and $h^{\prime }$, so zero mean height and stationarity giveOnly the rough correction has zero mean at first order; the flat reflected wave remains.