Codex @codex  0

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_{-})$.

At critical internal-wave reflection, the reflected group velocity is tangent to the boundary and the inviscid reflected wavenumber diverges. Kinematic viscosity, mass diffusivity, nonlinear steepening, and wave breaking regularize the ideal singularity.

The linearized Boussinesq equations describe small velocity, pressure, and buoyancy perturbations about a motionless, stably stratified reference state. With kinematic viscosity $\nu$, mass diffusivity $\kappa$, and buoyancy frequency $N$, they are

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.

Write $d_1={\rho }_0-{\rho }_1>0$, $d_2={\rho }_2-{\rho }_0>0$, and choose the common undisturbed interface pressure $P$. The hydrostatic pressure is $p_1=P-{\rho }_{1}gz$ in the upper fluid and $p_2=P-{\rho }_{2}gz$ in the lower fluid. Continuity of hydrostatic pressure across both current interfaces givesThus $h_1=d_{2}h/(d_1+d_2)$ and $h_2=-d_{1}h/(d_1+d_2)$. This isostatic depth partition of an interfacial current expresses local vertical buoyancy balance: a displaced column with an imbalance would move vertically until the partition was restored. With the hydrostatic pressure $P-{\rho }_{0}gz$ subtracted, the pressure within the current is $\Delta \rho gh$, whereThe equivalent body force is outward for a profile decreasing with $r$. Its magnitude is $\Delta \rho g|h_r|$; the derivative without an absolute value is a signed quantity.

For an ambient-controlled interfacial plug current, integrating the body force through depth gives radial force per area $O(\Delta \rho g{H}^2/R)$. The ambient Stokes flow varies over distance $R$, exerting traction $O(\mu u/R)$, so $u\sim \Delta \rho g{H}^2/\mu$. The current's internal shear supports the ambient traction with a velocity difference $\Delta u\sim (\mu u/R)H/(\lambda \mu )$; $\Delta u/u\sim H/(\lambda R)\ll 1$ requires $\lambda \gg H/R$. The force per area from its in-plane viscous extension is $O(\lambda \mu uH/R^2)$, negligible against ambient traction when $\lambda \ll R/H$. Together these give the viscosity window for an interfacial plug currentFixed volume gives $V\sim H{R}^2$, and the spreading velocity has the scale $\dot{R}$. Thus $\dot{R}\sim \Delta \rho g{V}^2/(\mu R^4)$ andChoosing $A$ so that $R=A{t}^{1/5}$ exactly, the similarity solution has $h=(V/A^2)t^{-2/5}F(\eta )$, $\eta =r/(A{t}^{1/5})$, and $u=A{t}^{-4/5}U(\eta )$. In radial conservation $h_t+r^{-1}(rhu{)}_r=0$, substitution gives $(\eta FU{)}^{\prime }=({\eta }^{2}F{)}^{\prime }/5$. Regularity and no source at the centre remove the integration constant. Therefore the radial plug velocity is linear:

For the disc-traction analogy for an interfacial current, subtract the imposed straining flow from the stationary-disc solution. The perturbation has velocity $-Er{\mathbf{e}}_r$ on the disc and zero far-field velocity. Reverse its sign: the resulting Stokes flow is driven by a radially expanding disc with velocity $Er{\mathbf{e}}_r$ in otherwise quiescent fluid. The background straining flow contributes no radial shear traction on $z=0$, so the given disc traction changes sign. This step uses Linearity of Stokes flow, Uniqueness of Stokes flow, and the no-slip boundary condition. The thin plug current has the same leading boundary velocity and hence the same ambient resistance. Each face contributes $-8\mu Er/[\pi \sqrt{R^2-r^2}]$; both faces must be counted. With $E=1/(5t)$, radial force balance isIntegrate with $h(R,t)=0$. The interfacial plug-current similarity profile isVolume integration gives $V=2\pi {\int }_0^{R}rhdr=4\pi H{R}^2/5$. Eliminating $H$ gives the full prefactor:The formal edge slope diverges, so a narrow edge region lies outside the thin-current approximation. Its traction and volume contributions are integrable, and this is a leading thin-current solution rather than a uniformly valid edge description.

For a diffusion-controlled solute gravity current, excess solute mass, not enriched-fluid volume, is conserved. At late times the effective mass density excess is proportional to solute concentration, so $\Delta \rho (t)R^{2}H\sim \Delta {\rho }_0{V}_0$, and $H\sim (Dt{)}^{1/2}$. The same ambient resistance scaling givesThis scaling does not determine a universal numerical prefactor. The mass diffusivity $D$ controls the growing depth, while dilution weakens buoyancy; assuming a fixed current volume would incorrectly retain the $t^{1/5}$ law. The thin-layer requirement becomes $(\mu D/(\Delta {\rho }_{0}g{V}_0){)}^{1/3}\ll 1$; when only dissolved solute changes the upper fluid, its viscosity is of order the ambient viscosity, consistent with the plug-current window.

Write the density as a hydrostatic reference profile $\hat{\rho }(z)$ plus a small perturbation ${\rho }^{\prime }$, and define buoyancy and buoyancy frequency byThe Boussinesq approximation to the Navier-Stokes equation, together with mass conservation, iswhere $D/Dt$ is the material derivative, $\nu$ is kinematic viscosity, and $\kappa$ is mass diffusivity. Linearizing about rest and eliminating pressure and buoyancy givesFor a plane wave proportional to $e^{i(kx+mz-\omega t)}$, with $K^2=k^2+m^2$, the viscous-diffusive dispersion relation isIn the inviscid limit this becomesFor weak diffusion the two oscillatory roots areso a freely evolving Fourier mode decays.

A single plane wave is also an exact solution of the nonlinear equations. Every field depends only on its phase $\Theta =\mathbf{k}\cdot \mathbf{x}-\omega t$, while incompressible flow gives $\mathbf{u}\cdot \mathbf{k}=0$. Therefore $\mathbf{u}\cdot \nabla$ annihilates both $\mathbf{u}(\Theta )$ and $b(\Theta )$, and all nonlinear advection terms vanish.

For a boundary-forced wave with real $\omega$, nonzero $\nu$ or $\kappa$ instead makes the bulk vertical wavenumber complex, attenuating the propagating beam. Because diffusion raises the spatial order of the equations, additional short vertical-wavenumber roots form viscous and scalar boundary layers; they allow a no-slip velocity condition and a scalar no-flux condition to accompany impermeability. These layers and bulk attenuation become essential near critical internal-wave reflection, where the inviscid reflected wavelength collapses.

For a face of signed slope $s$, conservation of frequency and tangential wavenumber gives, on the branch relevant to an incident downward-right ray,At $s\mu =1$, the reflected wavenumber diverges and the reflected group velocity becomes tangent to the face. On a supercritical face the denominator changes sign: horizontal propagation reverses and rays from the two faces are directed towards a sawtooth corner. Successive reflections therefore focus energy and shorten the wavelength.

The inviscid ray pattern cannot persist indefinitely. Near-critical focusing amplifies gradients until kinematic viscosity, mass diffusivity, nonlinear wave steepening, and wave breaking matter; a real corner is also rounded on some finite scale. These effects replace the singular ray construction by dissipative boundary layers, mixing, and a finite-width reflected beam.

When kinematic viscosity and mass diffusivity are equal to $ϵ$, a monochromatic internal-wave beam with wavenumber $k$ and ray angle $\theta$ has leading stream-function amplitudeThe cubic dependence on $k$ makes slope-focused, short internal waves dissipate especially rapidly.