Codex @codex  0

Biological fluid dynamics studies flows generated by and acting on organisms, cells, tissues, and biological transport networks.

For a passive scalar on a half-line, impose constant total solute flux at its inlet, with initially zero concentration and constant positive advection speed $U$ and mass diffusivity $\mathcal{D}$. A Laplace transform in time givesWith $\xi =Ux/\mathcal{D}$, $\vartheta =U^{2}t/\mathcal{D}$ and $z_{\pm }=(\xi \pm \vartheta )/(2\sqrt{\vartheta })$, its inverse isThe complementary error function describes the smeared front. This is a flux boundary condition, not the constant-concentration inlet solution. At large axial Péclet number, both have the leading advancing front $\frac{1}{2}\mathrm{erfc}(z_{-})$.

Salt lowers the freezing temperature of water. Consequently ice can ablate in saline water below the pure-water melting temperature if the ocean is warmer than its own liquidus. Solute diffusion towards the interface sustains the freezing-point depression. In the diffusive glacier-ablation interface law, small compositional Péclet number makes interfacial salinity close to ocean salinity and gives $P\simeq [T_{\infty }-T_m+\Gamma (C_{\infty }-C_s)]/[LD/(c_p\kappa )+\Gamma (C_{\infty }-C_s)]$ near the no-ablation threshold. This is diffusion-assisted dissolution rather than a supply of heat above the pure-water melting point; latent heat must still be supplied to complete the phase change.

For an almost spherical inviscid bubble in exterior Stokes flow $\mathbf{Ex}$, let $\gamma ={\gamma }_0-{\gamma }_1(C-C_0)$ and $K=C_0{a}^2/(2D_s)$. The quadrupolar surfactant distribution on a spherical interface makes the tangential velocity ${\mathbf{u}}_s=\alpha {\mathbf{I}}_s\mathbf{Ex}$ satisfyThe Marangoni stress approaches the stress needed to suppress surface motion as $M$ grows. The quadrupolar shape perturbation $r=a(1+\mathbf{n}\cdot \mathbf{Dn})$ isThese expressions require small deformation and small surface Péclet number; large $M$ alone does not justify a small concentration perturbation.

The bacterium samples its surroundings in time. A sustained increase in chemoattractant makes a favorable run persist by reducing the chance of tumbling; declining chemoattractant makes reorientation more likely. Individual tumbles need not point up the gradient. Biasing run durations can generate net migration even if reorientation itself is random. This is temporal comparison, not a direct measurement of a concentration gradient across the short cell body.

Compare the swimming time over length $\ell$, $t_a=\ell /v$, with diffusion over the same scale, $t_D\sim {\ell }^2/D$. Their ratio is the Péclet number $\mathrm{Pe}=v\ell /D$. The crossover estimate isOrder-one factors depend on the chosen diffusion-distance convention. A hypothetical stirrer of extent $d$ must instead have $Ud/D$ of order one, givingThese are transport-scale estimates, not a proof of a particular increase in uptake. Around a one-micrometre body, $vd/D\simeq 0.03$: local diffusive supply is already rapid. Actual uptake enhancement depends on the nutrient field and the cell's uptake boundary condition, neither supplied. The long run can help the organism reach a better nutrient region without implying that it greatly increases the flux from a uniform bath. E. coli normally uses bacterial flagella; the cilium in this comparison is a hypothetical generic stirring appendage, not its actual locomotor organelle.

A small swimming member of the bacteria samples too little spatial separation across its own body to compare weak chemical gradients accurately. Temporal sensing in bacterial chemotaxis instead compares recent receptor activity with an adapted memory. Straight runs carry it between locations; tumbles reorient it. When conditions improve, tumbling becomes less frequent and favorable runs last longer. An unfavorable direction is abandoned sooner. Repetition creates a biased random walk without requiring the cell to identify the best direction at each tumble. Rotational diffusion and environmental obstacles eventually destroy the benefit of maintaining a direction indefinitely.

An order-of-magnitude transport estimate compares the time $\ell /v$ to swim a distance $\ell$ with the time ${\ell }^2/D$ for an attractant to diffuse that distance. Their crossover is a Péclet number of order one:This estimates a useful run scale: over much shorter distances, nutrient diffusion rapidly replenishes the local environment compared with the organism's displacement. Runs of order a second also permit comparison along a substantial path before reorientation. The given $D$ and $v$ determine this scale, not an exact biological mean run length or a complete sensory optimization. If a one-coordinate root-mean-square diffusion length $\sqrt{2D\tau }$ is equated with $v\tau$, the resulting prefactor is two; a three-dimensional root-mean-square convention gives six. Such constants cannot be fixed by a scaling argument alone.

For a cilium of length $d$, compare its advective time $d/u$ with the local diffusion time $d^2/D$. Significant stirring relative to diffusion requires ${\mathrm{Pe}}_d=ud/D\gtrsim 1$, soThis is about $33$ times the stated swimming speed. If the tip travelled on a circular path of radius $d$, that speed would correspond to an illustrative frequency $u/(2\pi d)\simeq 160\mathrm{Hz}$. That geometric conversion is not a specified beat pattern. At microscopic scales the flow has small Reynolds number; useful transport depends on the actual nonreciprocal motion and resulting flow, and rapid motion alone does not guarantee increased uptake. The scaling explains why simply stirring with a very short appendage is demanding when molecular diffusion is already fast.

For stationary horizontally homogeneous means, incompressibility gives ${\bar{w}}_z=0$ and the impermeable boundary makes $\bar{w}=0$ throughout. With no imposed horizontal mean pressure gradient, the balances reduce toThe last equation includes the vertical Reynolds stress; omitting it requires an additional small-normal-stress approximation. The total vertical buoyancy flux is $B_0$ and the horizontal momentum flux is a signed constant $J_M$ with $|J_M|=U_{\ast }^2$. These are fluxes per unit reference density.

The local down-gradient closure defines eddy viscosity and eddy diffusivity byA parcel moving upwards typically retains the smaller momentum or scalar value from its previous height in an increasing mean profile; its fluctuation then correlates negatively with upward motion. This motivates positive coefficients, but is a closure assumption, not a universal theorem about turbulent or convective fluxes. The mean gradients obeyFor positive applied traction on a fluid above a bottom boundary, $J_M=+U_{\ast }^2$ and the mean velocity decreases away from the forcing boundary. The usual atmospheric convention specifies the stress absorbed by the boundary and has the opposite sign of $J_M$; scalar and Richardson-number results are unchanged. At large Reynolds number and Péclet number, outside the molecular wall sublayers, $K_M\gg \nu$ and $K_B\gg \kappa$, so molecular transport can be neglected in these flux-gradient relations. Large bulk numbers do not justify doing this exactly at a smooth solid wall.

The material derivative follows surfactant carried by the moving interface. Tangential surface divergence changes the local surface area and hence dilutes or concentrates insoluble material. Normal motion changes area through the curvature: an expanding sphere has positive $u_n{\nabla }_s\cdot \mathbf{n}$ and its concentration decreases. The last term is surface diffusion, which smooths concentration gradients. There is no exchange with the bulk; this is conservation of insoluble surfactant on a moving interface.

The interfacial speed is of order $a|\mathbf{F}|$. Diffusion dominates redistribution when the surface Péclet number satisfiesThen $C^{\prime }/C_0=O({\mathrm{Pe}}_s)$, and advection of $C^{\prime }$ is a second-order effect. Write $S=\mathbf{n}\cdot \mathbf{Fn}$. Since ${\mathbf{u}}_s=\mathbf{Fx}-S\mathbf{x}$, the projected divergence of the first term is $\mathrm{tr}\mathbf{F}-S=-S$, whereas the second has divergence $2S$; its gradient contracted with $\mathbf{x}$ is zero. Hence ${\nabla }_s\cdot {\mathbf{u}}_s=-3S$.

The leading steady transport equation is $D_s{\Delta }_s{C}^{\prime }=C_0{\nabla }_s\cdot {\mathbf{u}}_s=-3C_{0}S$. Trace-freeness and the supplied surface Laplacian identity give ${\Delta }_{s}S=-6S/a^2$. Fixing the conserved total surfactant removes the uniform null mode and yieldsThis is quadrupolar surfactant distribution on a spherical interface; its spherical mean is zero because $⟨n_i{n}_j⟩={\delta }_{ij}/3$.

With $+$ outside, $-$ inside, and the normal directed out of the bubble, the interfacial stress balance with variable surface tension isNow ${\nabla }_{s}S=2I_s\mathbf{Fn}/a$, so $-{\nabla }_s\gamma =2A{\gamma }_1{I}_s\mathbf{Fn}/a$. Multiplying the curvature by the tension and retaining linear terms givesThe undeformed isotropic pressure jump is kept with its normal; the perturbation terms are evaluated on the reference sphere. Products of deformation and concentration perturbations are beyond the approximation.

The imposed tensor is symmetric and traceless, so $\mathbf{x}\cdot \mathbf{Ex}$ is a harmonic scalar of degree two. Its decaying partner is $\mathbf{E}:\nabla \nabla (1/r)$; the corresponding harmonic vector uses $\mathbf{E}\nabla (1/r)$. Linearity and rotational symmetry give precisely these tensorial modes, with the growing scalar chosen to reproduce the imposed strain. Thus in the unscaled convention one can takeThe coefficients multiply distinct decaying modes. No other angular harmonic is forced at this order.

No penetration in the steady state gives $P=3Q-1$. The given tangential surface velocity implies $\alpha =1+2Q$. The tangential stress balance, with $M=A{\gamma }_1/(\mu a)$, gives $1+P-8Q=M\alpha$, or $-5Q=M(1+2Q)$. ThereforeThe nonuniform normal stress is unaffected by the uniform inviscid interior pressure. Its balance isSubstitution and rearrangement giveAs $M\to \infty$, $\alpha \to 0$: the concentration-induced Marangoni stress suppresses tangential surface motion, making the exterior boundary behave like an immobile surface despite the inviscid interior. The deformation coefficient tends to $5\mu a/(2{\gamma }_0)$; the clean limit is $2\mu a/{\gamma }_0$. Both small surface Péclet number and small deformation, of order $\mu a|\mathbf{E}|/{\gamma }_0$, are required. This is Marangoni immobilization of a bubble in straining flow.

Let $\mathbf{n}$ point from fluid $-$ to fluid $+$, $P_s=I-\mathbf{nn}$, ${\nabla }_s=P_s\nabla$, and $\mathcal{K}={\nabla }_s\cdot \mathbf{n}$. A clean material interface without phase change has continuous velocity, with normal interface speed $V_n={\mathbf{u}}_{\pm }\cdot \mathbf{n}$. Its surface stress is $\gamma P_s$. The surface divergence of that tensor is ${\nabla }_s\gamma -\gamma \mathcal{K}\mathbf{n}$. Balancing forces on a small interfacial patch therefore gives the interfacial stress balance with variable surface tensionThe normal condition is $\mathbf{n}\cdot ({\mathbf{\sigma }}_{+}-{\mathbf{\sigma }}_{-})\mathbf{n}=\gamma \mathcal{K}$; the tangential condition is $P_s({\mathbf{\sigma }}_{+}-{\mathbf{\sigma }}_{-})\mathbf{n}=-{\nabla }_s\gamma$. On a stationary sphere with outward normal these reduce to $p_{-}-p_{+}=2\gamma /a$ when viscous normal stresses vanish. This fixes the sign convention for the Marangoni stress rather than leaving the direction ambiguous.

For the bubble take negligible internal dynamic viscosity, no contaminants, and initially the standard effectively insulating gas-bubble approximation. With heat advection neglected, the exterior temperature satisfies Laplace's equation, approaches $T_0+\mathbf{H}\cdot \mathbf{r}$ remotely, and has zero normal heat flux at $r=a$. The temperature dipole around an insulating sphere givesThe zero-flux condition follows by differentiating the radial factor. Let ${\gamma }_T=(d\gamma /dT{)}_{T_0}$. To first order in the weak gradient,

In the laboratory frame the force-free bubble can be described by the harmonic Papkovich–Neuber representation potentialsThey give constant exterior pressure and the decaying source dipole fieldAt the surface $\mathbf{u}\cdot \mathbf{n}=\mathbf{U}\cdot \mathbf{n}$, so the normal kinematic boundary condition holds for translation. In coordinates aligned with $\mathbf{U}$, $u_r=a^{3}U\cos \theta /r^3$ and $u_{\theta }=a^{3}U\sin \theta /(2r^3)$. Direct differentiation givesEquivalently, the surface tangential traction is $(3\mu /a)P_s\mathbf{U}$. The tangential stress condition therefore requiresFor the usual negative ${\gamma }_T$, the bubble moves toward hotter liquid. The normal condition is consistent as well: its viscous dipole is $-6\mu \mathbf{U}\cdot \mathbf{n}/a=3{\gamma }_T\mathbf{H}\cdot \mathbf{n}$, exactly the dipolar part of $2{\gamma }_s/a$. The constant part fixes $p_{\mathrm{gas}}-p_{\infty }=2{\gamma }_0/a$. There is no Stokeslet term, hence no net applied force or drag imbalance.

The thermal assumption matters. The problem gives the exterior thermal diffusivity but no gas thermal conductivity. For an inviscid bubble with conductivity ratio $K=k_{\mathrm{gas}}/k_{\mathrm{liquid}}$, temperature and normal heat-flux continuity instead giveso the same stress calculation givesThe insulating-bubble result uses $K\ll 1$. If instead an undisturbed surface temperature gradient is imposed, corresponding to $K=1$, the coefficient is $1/3$ rather than $1/2$. Exterior thermal diffusivity alone does not select between those thermal boundary conditions. This dependence is explicit rather than silently treating a missing gas property as a prescribed value.

The small parameters can now be stated precisely. With $U_{\ast }=a|{\gamma }_T||\mathbf{H}|/(2\mu )$ and $\nu =\mu /\rho$, requireSmall Reynolds number justifies Stokes flow; small Péclet number allows quasi-static heat conduction without advective distortion; small capillary number and fractional surface-tension variation support a nearly spherical interface and a weak perturbation of its base surface tension. These last two conditions are proportional for this velocity scale but describe different approximations. Linearizing $\gamma (T)$ also needs $a|\mathbf{H}||{\gamma }_{TT}/{\gamma }_T|\ll 1$ if its derivative varies appreciably. Negligible internal viscosity and the insulating specialization require the corresponding gas-to-liquid viscosity and conductivity ratios to be small. Gravity has been excluded; if restored, its buoyancy-induced migration must be small compared with $U_{\ast }$, not merely small enough to avoid deformation.

Finally, surfactants change both the constitutive surface tension and the surface stress dynamics. The surface flow driven by the thermal Marangoni effect carries insoluble surfactant toward the high-tension side, lowering its tension and producing an opposing concentration-driven gradient. When surface diffusion or exchange is slow, this tends to immobilize the interface, suppressing the clean-bubble thermocapillary migration of an insulating bubble, sometimes nearly stopping it. If the surfactant stays uniformly equilibrated, it can instead change ${\gamma }_T$ without a large opposing gradient. Its concentration, surface diffusivity and adsorption kinetics must be supplied to compute a replacement speed; no universal numerical reduction factor follows merely from its presence.

The material derivative follows surfactant carried by the moving interface. The term $-C{\nabla }_s\cdot u_s$ is dilution or concentration by tangential area expansion or contraction. The term $-C(u\cdot n){\nabla }_s\cdot n$ is the additional area change produced by normal motion of a curved interface. The term $D_s{\nabla }_s^{2}C$ represents surface diffusion, while $-k(C-C_0)$ relaxes the concentration towards a reservoir value through exchange. All these terms refer to concentration per surface area, as in surfactant transport with exchange relaxation.

In the steady bubble frame, the spherical interface has $u\cdot n=0$. Write $C=C_0+C^{\prime }$ with uniform $C_0$ and regard both $u$ and $C^{\prime }$ as first order in the rising speed. Advection of $C^{\prime }$ and its product with surface dilation are second order. Hence the linear steady equation isSpherical symmetry and linearity leave only the tangential projection of the imposed translation vector: $u_s=A{I}_{s}U$. An azimuthal term $n\times U$ would change sign under reflection and is excluded for this non-chiral rising bubble. The degree-one forcing does not excite other spherical-harmonic degrees in the linear problem.

On $r=a$, differentiating $n=x/r$ givesEach component of $n$ is a degree-one spherical harmonic, soOne may also obtain the coefficient from $n\cdot n=1$: its surface Laplacian gives $n\cdot {\nabla }_s^{2}n=-|{\nabla }_{s}n{|}^2=-2/a^2$, while rotational symmetry makes the vector Laplacian parallel to $n$. Expanding $I_{s}U=U-n(U\cdot n)$ then yieldsThe candidate $C^{\prime }=BU\cdot n$ therefore solves the linear transport equation ifIts spherical mean is zero. This is the dipolar surfactant distribution. Assuming $k{a}^2+2D_s>0$, its consistency condition isDiffusion or exchange must restore uniform concentration faster than the interfacial flow redistributes it: this holds with a small surface Péclet number, fast exchange, or sufficiently immobilized surface velocity. These same conditions justify neglecting the perturbation's advection. Small capillary number and sufficiently small buoyancy-induced shape distortion are separate requirements for a spherical interface.

Take $+$ to mean the outer liquid, $-$ the bubble interior and $n$ outward from the bubble. The general interfacial stress balance with variable surface tension isFor $\gamma ={\gamma }_0-{\gamma }_1{C}^{\prime }$, its tangential part is $-{\nabla }_s\gamma ={\gamma }_{1}B{I}_{s}U/a$. Thus the requested form isThis positive parameter measures the strength of the restoring Marangoni stress relative to viscous stress, including the relaxation of concentration by diffusion and exchange.

The far velocity $-U$ is produced by the constant harmonic vector potential $\Phi =U$. The two independent decaying degree-one modes are the force-carrying vector potential $\alpha aU/r$ and the scalar dipole $\chi =\beta a^{3}U\cdot \nabla (1/r)$. They are harmonic outside the bubble and have precisely the symmetry of uniform translation. The Unscaled Papkovich–Neuber representation gives, at arbitrary radius,The velocity at $r=a$ agrees with the printed expression. However, the original PDF's traction formula has a factor-of-two error. Direct differentiation with $\sigma =-pI+\mu (\nabla u+\nabla u^T)$ givesnot the printed prefactor $12\mu /a$. This is traction of translating-sphere Papkovich–Neuber potentials. The $\alpha$ term follows immediately from the Stokeslet stress, and the scalar-potential part gives $6\mu \beta [U-3(U\cdot n)n]/a$; neither contribution contains the extra factor two.

With this consistent traction, no penetration requires $\beta -\alpha =1/2$, tangential velocity requires $A=-(1+\alpha +\beta )$, and tangential stress balance requires $\beta =A\lambda$. Solving,For comparison, taking the printed $12\mu /a$ literally would instead giveThose solve the supplied algebraic pair but not the Newtonian stress calculated from the supplied potentials. The factor-six result is the physically consistent translating surfactant-coated bubble solution.

As $\lambda \to 0$, $A\to -1/2$ and $\beta \to 0$: the clean bubble has a mobile interface and zero tangential traction. As $\lambda \to \infty$, $A\to 0$, $\alpha \to -3/4$ and $\beta \to -1/4$: the interface becomes immobilized and the exterior flow approaches that past a no-slip rigid sphere. The limiting tangential traction is $-3\mu I_{s}U/(2a)$, a useful check on the correct normalization.

The normal dynamic traction isapart from an arbitrary constant pressure. In the normal stress balance on a translating bubble, its degree-one part is balanced by the buoyant hydrostatic pressure difference and by the variation $-2{\gamma }_1{C}^{\prime }/a$ of normal capillary stress. A uniform internal gas pressure balances the constant Laplace-pressure term but cannot by itself balance this angular variation. If $U$ is upwards and $\Delta \rho$ is liquid minus bubble mass density, the degree-one balance can be writtenIt fixes $U=\Delta \rho g{a}^2(1+2\lambda )/[3\mu (1+3\lambda )]$, tending to the clean-bubble and rigid-sphere rise speeds. A degree-one radial displacement translates a sphere rather than changing its curvature, so a supposed dipolar shape distortion cannot replace this force balance.

Use a one-sided top-hat plume model with Boussinesq approximation and neglect the solid/liquid density difference in the interfacial volume source. Define the positive ablation speed $U=-\dot{a}$, volume flux per unit wall span $Q=bw$, and kinematic momentum flux $J=b{w}^2$. Let ambient entrainment per height be $E=\epsilon w$, with entrainment coefficient $\epsilon$, and take $\rho ={\rho }_{\infty }[1+{\beta }_C(C-C_{\infty })]$ for the salinity contribution to density. Conservation of volume and vertical momentum, with negligible wall stress and negligible vertical momentum of entering fluid, givesOnly the ocean-facing edge entrains. Using two edges would change the entrainment coefficient in these formulas.

To derive the heat balance, choose liquid specific enthalpy $h_l=c_{p}T+L$ and solid specific enthalpy $h_s=c_p{T}_s$. A horizontal slab receives ambient volume $Edz$ with enthalpy $h_{\infty }=c_p{T}_{\infty }+L$, and solid mass-equivalent volume $Udz$ with enthalpy $h_s$. ThusSubtract $(c_p{T}_{\infty }+L)Q^{\prime }$ and divide by $c_p$:Similarly, the incoming salt flux is $E{C}_{\infty }+U{C}_s$, so $(QC{)}^{\prime }=E{C}_{\infty }+U{C}_s$. Subtracting $C_{\infty }{Q}^{\prime }$ provesBoth anomalies are negative for fresh, cold melt entering the ocean. These balances account for the heat needed to melt the entering ice, rather than treating it as already-liquid water at $T_s$.

Now use $L/c_p\gg T_{\infty }-T_s$, the large latent-to-sensible heat ratio convention of the paper. Put $\Delta C=C_{\infty }-C_s>0$. The stated small plume-deficit limits allow us to replace plume $T,C$ by their ambient values in the diffusive interfacial fluxes. The Stefan condition and interfacial salt balance are thenFor the compositional Péclet number $P\equiv P{e}_c=U\delta /D$, solve the second equation and impose the liquidus:Substitution in the heat balance gives the diffusive glacier-ablation interface lawHere $\Gamma$ is the liquidus slope, unrelated to the surface-energy parameter of the previous question. If $A_T=LD/(c_p\kappa )$, $\Theta =T_{\infty }-T_m$ and $B_T=\Gamma \Delta C$, the physical root isIt is positive when $\Theta +B_T>0$, exactly the condition that the ocean is warmer than its own liquidus. Since $\delta$ and ambient properties are fixed, $U$ is constant with height in this approximation.

For the melting-driven saline wall plume, measure height $Z=z-z_0$ from a virtual origin. The integrated anomaly balances areIn the specified small-deficit regime the meltwater volume is negligible compared with entrainment. Indeed, the solution below gives $(C_{\infty }-C)/\Delta C=4U/(3\epsilon w)$, so a small compositional plume deficit implies $U\ll \epsilon w$. Thus retain $Q^{\prime }=\epsilon w$ for the leading power laws, although $U$ remains essential in the salt and heat sources. Try $b=BZ$, $w=W{Z}^q$. The salt deficit is proportional to $Z^{-q}$, and the momentum balance has powers $Z^{2q}=Z^{1-q}$, hence $q=1/3$. The volume balance gives $(4/3)BW=\epsilon W$, and the momentum coefficient gives $(5/3)B{W}^2=g{\beta }_{C}U\Delta C/W$. ThereforeThese are entrainment-dominated power laws, not an exact solution with the $U$ term retained in the volume equation. They fail close to the virtual origin, where the formal temperature and salinity deficits diverge. Their buoyancy source is distributed melting along the wall, unlike a line plume with a constant buoyancy flux from a localized source.

The profiles rise from $T_i,C_i$ across the laminar boundary layer to the almost-ambient plume values $T,C$, then approach $T_{\infty },C_{\infty }$ across the entraining edge. The sketch uses the linear-gradient approximation implicit in the prescribed diffusive fluxes; it is not an exact normal advection–diffusion profile at large $P$.

For $P\ll 1$,The last formula retains both first-order terms and is appropriate near the no-ablation threshold. Salt diffusion keeps the interface saline, and freezing-point depression allows ice dissolution in saline water even when $T_{\infty }<T_m$. Latent heat is still needed, but the thermodynamic permission to ablate comes from the depressed liquidus; small $P$ does not mean a negligible energy cost.

For $P\gg 1$, $C_i-C_s\sim \Delta C/P$ and $T_i\sim T_m-B_T/P$. Rapid production of fresh melt dilutes the interface, so its temperature approaches the pure-ice melting point. In the usual thermally driven regime $\Theta >0$ and $\Theta \gg B_T/P$,This is predominantly thermal melting. More generally the retained large-$P$ balance is $A_{T}P\simeq \Theta +B_T/P$. For example, $\Theta =0$ gives $P\sim \sqrt{B_T/A_T}$, which can be large while the driving term is still saline freezing-point depression. Thus the melting/dissolution interpretation must also consider the ocean's temperature relative to $T_m$, not the compositional Péclet number alone.

Use a constant porosity $\varphi$ and let $v$ denote pore velocity. Darcy law gives the Darcy velocity $q=\varphi v$. If the velocity convention has already absorbed porosity, set $\varphi =1$ in the following formulas. Assume positive permeability of a porous medium throughout the layer, so $k_0>0$ and $k_0+k_1>0$.Here $\mu$ is dynamic viscosity. The flux-weighted residence time in a porous layer is the pore volume divided by throughput. Equivalently, weighting each streamline's transit time $L/v(y)$ by its inlet flux $v(y)dy$ givesThis is also the advective travel scale after cross-layer mixing has homogenized the concentration. For resident-weighted transit time in a layered flow, sampling the initial fluid uniformly by resident volume, while suppressing cross-stream exchange, is a different experiment: its mean iswith the continuous $k_1=0$ limit $\mu \varphi L^2/(\Delta P{k}_0)$. Specifying the sampling convention avoids confusing the reciprocal of the mean speed with the mean reciprocal speed.

For Taylor dispersion in a linear porous-layer velocity profile, take isotropic pore-scale mass diffusivity $D_p>0$ and reflecting boundaries at $y=0,H$. Write $v-U=\gamma (y-H/2)$, where $\gamma =\Delta P{k}_1/(\mu \varphi LH)$. Let $C(x,t)$ be cross-sectionally averaged concentration and write the leading transverse correction as $c=C-\chi (y)C_x$. Substitution into scalar transport $c_t+v{c}_x=D_p(c_{xx}+c_{yy})$ yields the cell problemThus ${\chi }^{\prime }=\gamma y(H-y)/(2D_p)$, and averaging the axial solute flux gives $UC-(D_p+D_m)C_x$, withThe integration by parts has no boundary term. It also proves nonnegativity, even when $k_1<0$. The total effective mass diffusivity is $\mathcal{D}=D_p+D_m$. For anisotropic pore dispersion, the denominator uses the transverse coefficient and the added molecular term uses the axial coefficient.

There are two separate tests for significance. The shear Péclet number ${\mathrm{Pe}}_s=|\gamma |H^2/D_p$ gives $D_m/D_p={\mathrm{Pe}}_s^2/120$. A large value implies substantial enhancement of axial spreading. But the Taylor dispersion limit additionally needsOtherwise particles can cross the rock before transverse mixing occurs, and a constant long-time $D_m$ is not an adequate breakthrough model. The relative front width is of order $\sqrt{\mathcal{D}/(UL)}$ once the long-time model applies.

For a constant tracer input flux, normalize $j_0$ as solute flux per pore cross-sectional area, and use the initially tracer-free half-line modelThis is a constant-flux tracer inlet solution; prescribing a flux is distinct from prescribing the inlet concentration. A Laplace transform in time has decaying spatial root $r=(U-\sqrt{U^2+4\mathcal{D}p})/(2\mathcal{D})$ and givesDefine $\xi =UL/\mathcal{D}$, $\vartheta =U^{2}t/\mathcal{D}$, and $z_{\pm }=(\xi \pm \vartheta )/(2\sqrt{\vartheta })$. Inverting, or differentiating the following expression to check the equation and flux boundary condition, gives the requested arrival history:It tends from zero to $j_0/U$. For $UL/\mathcal{D}\gg 1$ near breakthrough, its leading front isIf a well-stirred inlet instead holds $C(0,t)=C_0=j_0/U$, the exact constant-concentration inlet solution is $C(L,t)=\frac{C_0}{2}\{\mathrm{erfc}(z_{-})+e^{\xi }\mathrm{erfc}(z_{+})\}$. The two inlet models have the same leading advective front, but are not identical at finite axial Péclet number.

The material derivative follows the interfacial surfactant. The term $-C{\nabla }_s\cdot {\mathbf{u}}_s$ describes dilution by tangential expansion and concentration by compression. The term $-C(\mathbf{u}\cdot \mathbf{n}){\nabla }_s\cdot \mathbf{n}$ accounts for changing surface area due to normal motion of a curved interface. Finally, $D_s{\nabla }_s^{2}C$ describes surface diffusion. Together these give conservation of insoluble surfactant on a moving interface.

The relevant surface Péclet number isWrite $S=\mathbf{n}\cdot \mathbf{A}\cdot \mathbf{n}$ and ${\mathbf{u}}_s=\mathbf{Ax}-S\mathbf{x}$. Since $\mathbf{A}$ has zero trace and ${\nabla }_s\mathbf{n}={\mathbf{I}}_s/a$, the surface divergence isThere is no normal motion. To first order in ${\mathrm{Pe}}_s$, the steady surfactant balance becomes $D_s{\nabla }_s^2{C}^{\prime }=C_0{\nabla }_s\cdot {\mathbf{u}}_s=-3C_{0}S$. The traceless quadratic $S$ is a degree-two spherical harmonic, so ${\nabla }_s^{2}S=-6S/a^2$. The mean of $C^{\prime }$ is zero by total surfactant conservation; consequentlyThe discarded advective term ${\mathbf{u}}_s\cdot {\nabla }_s{C}^{\prime }$ is smaller by ${\mathrm{Pe}}_s$. This is the quadrupolar surfactant distribution on a spherical interface.

With the unit normal directed from the bubble into the exterior, the interfacial stress balance with variable surface tension isUsing the specified first-order curvature and $\gamma ={\gamma }_0-K{\gamma }_{1}S$, together with ${\nabla }_{s}S=2{\mathbf{I}}_s\mathbf{An}/a$, givesThe first term is the spherical capillary pressure; the last contains the normal tension correction and the tangential Marangoni stress.

The ambient Stokes flow has the symmetry of the symmetric traceless tensor $\mathbf{E}$. In the Unscaled Papkovich–Neuber representation, ${\chi }_{\infty }=\mathbf{x}\cdot \mathbf{E}\cdot \mathbf{x}/2$ is harmonic and yields $\mathbf{Ex}$. The decaying vector potential must have the form $\mathbf{E}\nabla (1/r)$, and the scalar disturbance must have the degree-two form $\mathbf{E}:\nabla \nabla (1/r)$. They are harmonic outside the bubble and have precisely the required rotational covariance. Dimensionless amplitudes may therefore be writtenFor a steady bubble, the no-penetration boundary condition gives $1+P-3Q=0$. Its tangential velocity gives $\alpha =1+2Q$. The tangential stress boundary condition then yieldsHence $\alpha =5/(5+2M)$, $Q=-M/(5+2M)$ and $P=-(5+5M)/(5+2M)$. The inviscid interior supplies only the constant pressure balancing $2{\gamma }_0/a$. Matching the remaining normal stress givesEliminating $P,Q,\alpha$ givesAs $M\to \infty$, $\alpha \to 0$: the Marangoni stress suppresses tangential motion and effectively immobilizes the interface. This limit must retain the small surface Péclet number and small-deformation assumptions; $M$ can grow through increasing ${\gamma }_1$ without invalidating the linear concentration approximation.