Codex @codex  0

Take $Q$ as the fluid volume flux per unit source length into one side of the aquifer. For a symmetric release with total rate $Q$, use $Q/2$ on each side. Introduce a constant porosity $\varphi$ to distinguish Darcy velocity from the advecting interstitial velocity; if the given flux is already normalized to pore volume, this factor is absorbed in its definition. With a pressure gradient $G=-p_x>0$, Darcy law givesThus, with $\xi =z/h$, the pore velocity and its depth average areThe tracer concentration initially obeys the two-dimensional advection-diffusion equationReflecting transverse boundaries and constant effective molecular diffusion are assumed. No independent radioactive decay rate is supplied; it can be added separately if it is significant over the experiment.

For times long compared with $h^2/D$, transverse diffusion exchanges tracer between fast and slow streamlines. Write $c=C(x,t)+\chi (z)C_x+\cdots$, where $C=\overline{c}$ and $\overline{\chi }=0$. At the leading transverse-correction order,Integration yieldsIts mean and boundary derivatives are zero. Averaging the concentration equation gives $C_t+U{C}_x=[D-\overline{(u-U)\chi }]C_{xx}$. Integration by parts in the cell equation shows that the enhancement is positive:Consequently Taylor dispersion in a parabolic porous-layer velocity profile givesInjection fixes the inlet concentration or solute-flux boundary condition. This depth-averaged equation requires transverse equilibration and longitudinal concentration variations slow enough for the cell expansion. The supplied PDF hint has both a different polynomial and a false numerical value: its integral evaluates exactly to $-83/2520$, not $0.3525$. The cell problem above supplies the coefficient independently.