Codex @codex  0

Use pressure with the ambient hydrostatic pressure removed, and let $\Delta P$ be the pressure below the falling sphere minus that above it. Away from the sphere the tube flow is Hagen-Poiseuille flow. The upper and lower clear lengths add to $L$ to leading order in $a/L$, so their pressure drops add:The large exterior reservoir returns this flux with negligible pressure loss compared with the long narrow tube. The sphere is assumed far enough from the ends that entrance effects and its occupied length give only lower-order corrections.

Now use a frame translating downwards with the sphere, and take $z$ positive upwards from its center. The sphere is stationary and the tube wall moves upwards with speed $U$. The narrow gap isIt has axial length scale $a\sqrt{ϵ}$, making lubrication theory appropriate. The sphere-frame flux in a tube is constant: the upward relative flux through the whole section is $\pi a^{2}U-Q$, so its flux per unit circumference isWith $y=0$ on the sphere and $y=h$ on the wall, Couette-Poiseuille flow in a thin gap givesSince $\Delta P=-\int p_{z}dz$, the pressure-jump relation isThe parabolic inner gap can be extended to infinite $z$ because these integrals are dominated by $|z|=O(a\sqrt{ϵ})$. The lubrication resistance integrals for a nearly occluding sphere areSetting $\xi =\tan \theta$ evaluates $J_1=\pi$, $J_2=\pi /2$, $J_3=3\pi /8$, hence

For the force, enclose the fluid between sections just above and below the sphere, including the annular layer. The net upward pressure force is $\pi a^2\Delta P$. The tube wall exerts upward shear traction $\mu u_y(h)$ on this fluid; the sphere exerts a downward force equal to the upward hydrodynamic force $F$ on itself. From the gap solution,The axial control volume force balance therefore gives the control-volume drag formula for a sphere in a tube:Using a control volume avoids having to integrate the curved-sphere normal pressure and tangential shear separately.

With $Q^{\ast }=Q/(\pi a^{2}U)$, $\Delta P^{\ast }=a\Delta P/(6\mu U)$, $q^{\ast }=2q/(aU)$, $F^{\ast }=F/(6\pi \mu aU)$, $I_n^{\ast }=a^{n-1}{I}_n$ and $L^{\ast }=4L/(3a)$, the three flux/pressure equations becomeSolving these linear relations givesThe force relation is $F^{\ast }=\Delta P^{\ast }+\frac{4}{3}{I}_1^{\ast }-q^{\ast }{I}_2^{\ast }$. Substitution first yieldsBecause $I_2^{\ast }/I_3^{\ast }=2ϵ/3$ and $I_1^{\ast }/I_3^{\ast }=2{ϵ}^2/3$, the terms $-2L^{\ast }{I}_2^{\ast }$ and $\frac{4}{3}{L}^{\ast }{I}_1^{\ast }$ are lower order than $L^{\ast }{I}_3^{\ast }$. Keeping the uniformly relevant leading terms givesThe squared $I_2^{\ast }$ is essential; replacing it by $I_2^{\ast }{I}_3^{\ast }$ would give the wrong sign and scaling. This establishes the three drag regimes for a nearly occluding sphere.

For $L^{\ast }\ll {ϵ}^{-1/2}$,The sphere nearly carries the tube's fluid down with it: the far-field flux is close to the piston displacement rate. The relative annular leakage is small in total volume, and the dominant force comes from the local annular shear term. The pressure force associated with the overall tube resistance is smaller.

For ${ϵ}^{-1/2}\ll L^{\ast }\ll {ϵ}^{-5/2}$,The flow remains approximately piston-like, but the dominant force is now the pressure needed to drive Hagen-Poiseuille flow through the long clear tube. Both displayed terms in $q^{\ast }$ are retained because their relative size changes inside this regime without producing an additional leading drag regime.

For $L^{\ast }\gg {ϵ}^{-5/2}$,The large tube resistance makes its net through-flow negligible. The sphere's displaced fluid instead passes upwards relative to it through the thin annular gap, with gap velocity of order $U/ϵ$. Its pressure drop and resulting pressure force dominate.

Finally consider two identical co-moving spheres, with nonoverlapping lubrication regions and the same total clear-tube length to leading order. The common through-flux is not doubled. In the middle regime it is still approximately $\pi a^{2}U$, so the same total tube pressure drop is shared equally between the two identical spheres. The load sharing between co-moving spheres in a tube gives half the single-sphere force on each sphere in regime (ii). In regime (i), drag is set locally by the nearly pressure-balanced annular shear and is unchanged on each sphere. In regime (iii), each sphere must pass essentially its own displacement flux through its annular gap, giving the same local pressure drag as before; adding another pressure jump changes the already small tube flux but not either leading drag. The force is unchanged in regimes (i) and (iii).

For $\mathbf{r}=R\hat{\mathbf{r}}$ with the observation distance much larger than the scatterer, expand $|\mathbf{r}-{\mathbf{r}}^{\prime }|=R-\hat{\mathbf{r}}\cdot {\mathbf{r}}^{\prime }+O(R^{-1})$. The far-field pattern in the Born approximation for scalar wave scattering is consequentlyThis convention includes the factor $1/(4\pi )$ in $f_{\infty }$; moving that factor into the outgoing-wave definition would change the pattern's normalization.

The absolute energy flux depends on what physical amplitude $\psi$ represents. For acoustic pressure in a homogeneous fluid of mass density ${\rho }_0$ and wave speed $c_0$, the momentum equation gives the particle-velocity amplitude $\mathbf{v}=\nabla \psi /(i{\rho }_0\omega )$. The time average of harmonic power and acoustic intensity then giveThe radial derivative of the outgoing factor is $(ik-1/R)e^{ikR}/R$. Its imaginary current is therefore $k|f_{\infty }{|}^2/R^2$ to leading order. The incident unit-amplitude plane wave has imaginary current $k{\hat{\mathbf{r}}}_0$. Hence, with $k=\omega /c_0$,The scattered energy flux points along $\hat{\mathbf{r}}$ at leading order, and the incident energy flux along ${\hat{\mathbf{r}}}_0$. If $\psi$ is a velocity-potential amplitude instead, both leading expressions use the common prefactor ${\rho }_0\omega k/2$ in place of $k/(2{\rho }_0\omega )$. A normalized scalar current uses $k$ as the common prefactor. The paper does not specify this amplitude normalization; the ratio in part (c) is independent of that choice.

Use the Fourier transform convention $\hat{p}(k,y)={\int }_{\mathbb{R}}{e}^{ikx}p(x,y)dx$. Let $\gamma (k{)}^2=k^2-k_0^2$, choosing $\mathrm{Re}\gamma >0$ on the real contour with $\mathrm{Im}\omega <0$. For time dependence $e^{i\omega t}$, this is the decaying continuation of the outgoing Sommerfeld radiation condition.

The transformed Helmholtz equation has upper and lower solutions $P(k)e^{-\gamma y}$ and $Q(k)e^{\gamma y}$. The two kinematic traces give $-\gamma P={\rho }_0{\omega }^2\hat{\eta }=\gamma Q$, henceThe elastic membrane equation gives $[\hat{p}]=(m{\omega }^2-T{k}^2)\hat{\eta }$. Therefore the acoustic wave on a tensioned massive membrane has dispersion relationEquivalently, $T{k}^2={\omega }^2(m+2{\rho }_0/\gamma )$. The fluid on both sides supplies a positive added mass of an evanescent fluid layer, which is a useful independent check on the sign. In particular, the incident pressure isThe symbol $m$ here denotes elastic membrane mass per area, rather than the fluctuating Mach number used in Question 1.

Linearity permits subtraction of the incident acoustic membrane wave. The scattered pressure solves the homogeneous Helmholtz equation, and the kinematic relation defines its line displacement ${\eta }^{\prime }$ on both halves. On $x<0$, subtracting the incident elastic membrane equation gives the scattered dynamic equation. On $x>0$, there is no elastic membrane, so the total pressure jump must vanish, giving $[p^{\prime }]=-[p_I]$. The pinned endpoint has $\eta (0)=0$, hence ${\eta }^{\prime }(0)=-1$, rather than zero.

Use full and Half-range Fourier transforms with the same $e^{ikx}$ convention as part (a), and putWrite $P^{-}$ for the left transform of the upper scattered pressure, $P^{+}$ for its right transform, and ${\eta }^{\pm }$ for the corresponding transforms of ${\eta }^{\prime }$. Outgoing antisymmetry in $y$ impliesThe second relation comes from the prescribed pressure cancellation on $x>0$.

The endpoint terms in the left Fourier transform of ${\eta }_{xx}^{\prime }$ are crucial:Thus $2P^{-}=a{\eta }^{-}+TA+iT(k+k_I)$. Eliminating ${\eta }^{-}$ gives the Wiener-Hopf equationThere is a defect in the printed right side: with the stated incident exponential, pressure jump, and pinned total displacement, it omits the incident endpoint term and has the opposite sign on the incident-pressure term. The displayed corrected equation follows directly from both boundary traces. The given right side cannot be derived with these definitions. Using $a(k_I)=-2\beta /{\gamma }_I$, an especially useful equivalent form isThe correction is necessary for the reconstructed pressure to satisfy the original physical boundary conditions.

For Wiener-Hopf factorization, take $K=K^{+}{K}^{-}$ with the plus factor analytic and nonzero in the upper half-plane and the minus factor analytic and nonzero in the lower half-plane. Singularities listed next refer to continuation out of each factor's own analytic half-plane. Let $b=\omega \sqrt{m/T}$, so $a=-T(k-b)(k+b)$ and $b$ lies below the real axis. Generically $K^{+}$ inherits the pole at $b$ and the square-root branch point at $k_0$, while $K^{-}$ inherits the pole at $-b$ and the branch point at $-k_0$. The kernel is finite but nonanalytic at the acoustic branch points: its local behavior is a constant plus a square-root term, not an inverse-square-root divergence. Coincident branch points and poles require a limiting treatment.

Zeros of $K$ are the fluid-loaded dispersion relation roots. Lower-half-plane roots, including the incoming guided root $k_I$, belong to the continued $K^{+}$; upper-half-plane roots, including the reflected root $-k_I$ for the even kernel, belong to the continued $K^{-}$. The roots represent membrane-guided modes; the branch cuts represent radiated acoustic waves. Only roots on the selected outgoing sheet are included, not spurious roots created by squaring the dispersion relation. No explicit factorization is required.

Here is an explicit solution in terms of those factors. DefineThe pole of $B$ at $-b$ has residue $C$; the pole at $b$ is canceled by the pole of $K^{+}$. Thus $B^{\pm }$ have the required respective analyticity. Divide the corrected Wiener-Hopf equation by $K^{+}$ and split its right side as $F^{+}+F^{-}$, whereThe difference quotient in $F^{-}$ is removable at $k_I$ and is analytic below. Moving the plus and minus terms to opposite sides yields the common entire function. With the assumed $E(k)=0$, the solution isCombining $P^{-}$ with the known $P^{+}$ gives the full scattered pressure transformThe requested pressure integral, for $y\ne 0$, is consequentlyThe real contour $\mathcal{C}$ uses the causal continuation $\mathrm{Im}\omega <0$; its undamped limit retains the induced pole and branch-cut prescriptions. The residue theorem shows that upper-half-plane zeros of $K^{-}$ yield left-going scattered elastic membrane modes. On the right, the lower incident-pole residue of the scattered integral cancels $p_I$ at the open line, as required.

Finally reconstruct the left scattered displacement:For a generic unspecified $A$, this has a lower-half-plane pole at the bare elastic membrane wavenumber $b$. Such a pole is incompatible with analyticity of an outgoing left-supported scattered displacement: it would represent an additional right-going incoming elastic membrane contribution. It must be removed. Its residue is $-[2P^{-}(b)-TA-iT(b+k_I)]/(2Tb)$, so the incoming-pole cancellation at a pinned membrane edge condition isSince $P^{-}$ is linear in $A$, this fixes the endpoint slope generically. Explicitly it isThe physical lower incident pole of the total displacement is already prescribed by ${\eta }_I$; it must not be confused with this removable spurious pole of the scattered field.

For a stationary, uniformly permeable skeleton, Darcy law gives $\mathbf{u}=-(\Pi /\mu )\nabla p$. Constant-density conservation of mass gives $\nabla \cdot \mathbf{u}=0$, hence the pressure is harmonic:With uniform gravitational acceleration, the same argument applies to pressure after subtracting its hydrostatic part. The spherical calculation below neglects body-force variation and uses the excess pressure relative to the far field.

Set $ϵ=1-{\rho }_s/\rho >0$ and $\beta =\varphi ϵ$. When a shell freezes, its pore mass decreases at rate $4\pi a^2\varphi (\rho -{\rho }_s)\dot{a}$; the difference must leave as liquid. Since Darcy velocity is volume discharge per total area, not pore-fluid speed, the interfacial and radial balances giveThe positive sign means outward expulsion by the expansion on freezing. Integrating Darcy law with $p(\infty )=0$ givesThere is no liquid Darcy velocity in the ice-filled region.

To state a closed Stefan problem, let $k_{-},k_{+}$ be the effective thermal conductivities of the frozen and unfrozen composites and $C_{-},C_{+}$ their volumetric heat capacities. If the matrix mass density is ${\rho }_m$, the common specific heat gives $C_{-}=c_p[(1-\varphi ){\rho }_m+\varphi {\rho }_s]$ and $C_{+}=c_p[(1-\varphi ){\rho }_m+\varphi \rho ]$. These material coefficients are not numerically specified. With $\mathcal{L}=\varphi {\rho }_{s}L$, the thermal equations areThe exterior equation retains advected sensible heat. A more detailed mechanical-energy model adds the Darcy dissipation source $\mu u^2/\Pi$; its smallness at leading order is checked below. At the origin require ${\partial }_r{T}_{-}(0)=0$, and at infinity require $T_{+}\to T_{\infty }$, $p\to 0$. At the phase boundary impose temperature continuity, phase equilibrium and the Stefan condition:These equations, the radial liquid-mass balance and Darcy law, together with an initial radius and compatible initial temperature fields, specify the moving-boundary model. No curvature correction is imposed. A vanishing initial seed corresponds to the leading growth law below with $a(0)=0$.

Write $\Delta =T_m-T_{\infty }$ and $S=L/(c_p\Delta )$, the latent-to-sensible heat ratio. A conductive interfacial balance gives $a\dot{a}\lesssim k_{+}\Delta /\mathcal{L}$. On length scale $a$, the ratio of thermal diffusion time to growth time is consequently bounded byIt is small for $S\gg 1$ with fixed porosity and comparable constituent densities and thermal properties. Very small porosity requires the sharper condition $C_{\pm }\Delta /(\varphi {\rho }_{s}L)\ll 1$; the bare heat-capacity ratio alone should not be used uniformly as $\varphi \to 0$. The advective-to-conductive ratio isEven if mechanical heating is retained, its size relative to latent-heat release is $p_i\beta /\mathcal{L}=(T_m-T_i)/T_m\le \Delta /T_m$. For fixed water properties this is also small as $S\to \infty$, consistently with the linearized Clausius-Clapeyron relation. Thus unsteadiness, thermal advection and mechanical heating are higher-order effects. The local leading temperature fields satisfy Laplace's equation. The approximation describes radii of order $a$; the distant transient diffusion tail supplies the far-field matching and need not be globally steady.

Regularity at zero and the interface/far-field values giveThe Stefan condition becomes $\mathcal{L}a\dot{a}=k_{+}(T_i-T_{\infty })$. Introduce the dimensionless pressure-feedback factorCombining this heat balance with pressure-dependent phase equilibrium yields the growth and interface values:This is pressure-limited spherical freezing in a porous medium. The interface values are constant during this leading solution, although $p(r,t)$ changes at a fixed exterior radius.

As permeability of a porous medium decreases, $\Lambda$ increases: the interface pressure rises toward the finite value $\Delta /\gamma$, while the interface temperature falls toward $T_{\infty }$. Hydraulic resistance suppresses freezing through pressure-induced freezing-point depression. In the limit $\Pi \to 0$, $a\dot{a}\sim \Pi \Delta /(\gamma \mu \beta )$; a zero-seed radius grows as $\sqrt{\Pi t}$ and tends to zero at fixed time. A finite seed instead hardly changes. At exactly zero permeability an incompressible rigid model cannot expel the expansion volume, so sustained growth stops; deformation or compressibility would require a different model.