Codex @codex  0

On $z>0$, consider $p_t+a{p}_z=D(z{p}_z{)}_z$ with $a,D>0$, zero scalar flux $ap-Dz{p}_z$ at the endpoints for $t>0$, and unit initial impulse at the origin. Put $r=a/D$ and $\eta =z/(Dt)$. The normalized similarity solution isThe ansatz $p=(Dt{)}^{-1}f(\eta )$ gives $[\eta f^{\prime }+(\eta -r)f{]}^{\prime }=0$. Endpoint decay sets this integrated constant to zero and hence $f\propto {\eta }^r{e}^{-\eta }$; the gamma function normalizes its integral to one. Its scale is $Dt$, so it converges weakly to a unit impulse as $t\downarrow 0$. Its gamma distribution shape is $1+r$, its maximum occurs at $z=at$, and its expected value is $(a+D)t$. For width proportional to $z$, the associated horizontally averaged plume concentration is proportional to $z^{r-1}{e}^{-z/(Dt)}$: its interior maximum is $(a-D)t$ when $a>D$, its supremum occurs at the origin when $a=D$, and it is singular there when $0<a<D$. The physical finite source regularizes this ideal-origin behaviour.

The horizontally integrated plume concentration is the integral of a local concentration $C$ across the plume. It measures scalar amount per unit vertical length and unit span, and its vertical integral is total scalar amount per unit span. Unlike a local concentration, its dimensions include one horizontal length. It must be distinguished from horizontally averaged plume concentration when a plume's width varies with height.

Use a two-sided top-hat line plume with full width $w(z)$, upward speed $U(z)$, and reduced gravity $b(z)$ uniform across the plume. Fluxes are measured per unit span along the line source: the volume flux is $V=wU$, the kinematic momentum flux is $M=w{U}^2$, and the buoyancy flux is $B=wUb$. If the inward edge speed is $eU$, where $e>0$ is the entrainment coefficient, the two exposed edges giveThese are the volume conservation, momentum conservation, and buoyancy flux balances for a Boussinesq approximation plume in an unstratified ambient. Entrained ambient fluid supplies neither vertical momentum nor reference buoyancy.

A pure plume has no persistent source length scale. Since a line-source $B$ has dimensions $L^3{T}^{-3}$, dimensional analysis gives constant $U\propto B^{1/3}$, $w\propto z$, and $V\propto z$. Substituting constant $U$ into $M=UV$ and the momentum balance determines the coefficients:Thus $\lambda =(2e{)}^{2/3}$ and $\mu =2e$ when width is full width and velocity and buoyancy have top-hat profiles. If width means half-width, $\mu =e$ instead. Other prescribed profile shapes change these numerical factors; the linear growth laws remain the same. The point-source $B$ in Question 1 has different dimensions, so its $z^{5/3}$ law must not be used here.

Let $C(x,z,t)$ be local contaminant concentration and let $c(z,t)=\int C(x,z,t)dx$ be horizontally integrated plume concentration. Its conserved amount per unit span is ${\int }_0^{\infty }cdz=V_c$. The dilute contaminant is a passive scalar; the maintained plume is unaffected by its impulsive release. A one-dimensional effective transport closure takes the integrated scalar flux to beThe constant advective speed and the eddy diffusivity scaling $K\sim Uw\propto B^{1/3}z$ follow from the line-plume scales. The entrainment coefficient alone does not determine the scalar dispersion coefficient: this Fickian closure for the integrated variable is an additional modelling assumption. Mass conservation $c_t+J_z=0$ then gives the displayed transport equation in the PDF. With this effective closure and advective speed chosen as $U$, $\alpha =(2e{)}^{-1/3}$; if $K=\kappa Uw$, then $\beta =\kappa (2e{)}^{2/3}$, with independent dimensionless mixing coefficient $\kappa$.

There is a second possible closure convention. If one instead applies local diffusive flux to horizontally averaged plume concentration $\overline{C}=c/w$, its integrated diffusive flux is $-Kw{\overline{C}}_z=-K(c_z-c/z)$. For $K=Dz$, the total scalar flux is then $(U+D)c-Dz{c}_z$. The same printed equation is obtained with $a=U+D$, rather than $a=U$. Thus the printed $\alpha$ and $\beta$ should be regarded as effective coefficients unless the averaging and closure conventions are specified; no equality between $\alpha$ and the velocity prefactor is universal.

Assume $a,D>0$, $z>0$, an initial impulse $V_c{\delta }_0$ at the origin, no further contaminant input, and zero endpoint scalar flux for $t>0$. For a similarity solution, putSubstituting into the advection-diffusion equation givesDecay at infinity makes the integrated constant zero. Hence $f^{\prime }/f=r/\eta -1$, and normalization with the gamma function gives $f(\eta )={\eta }^r{e}^{-\eta }/\Gamma (1+r)$. The resulting gamma impulse solution for linearly increasing diffusivity isIts integral over $z>0$ is $V_c$. Its flux is $J=(z/t)c$, which vanishes at both endpoints for $t>0$. The scale $Dt$ shrinks to zero as $t\downarrow 0$, so the normalized solution has weak convergence of probability measures to the required unit source impulse. The mass at the boundary is a full unit impulse on the half-line, not half of a whole-line impulse.

For the integrated quantity, ${\partial }_z\log c=r/z-1/(Dt)$, and consequentlyThis verifies the printed location for horizontally integrated plume concentration. It is a maximum because the derivative changes from positive to negative there. The normalized profile is a gamma distribution with shape $1+r$: its expected value is $(a+D)t$ and its variance is $(a+D)Dt$, so neither the mean position nor the spread should be mistaken for the modal position.

The PDF changes from “integral” to “averaged” concentration in its final request. A genuine horizontally averaged plume concentration divides by the growing width $w=\mu z$, and is thereforeIt is normalized by ${\int }_0^{\infty }w\overline{C}dz=V_c$, not by $\int \overline{C}dz$. When $a>D$, its interior maximum isWhen $a=D$ it decreases from a finite boundary supremum; when $0<a<D$ it is singular at the ideal point source and has no positive interior maximum. A finite source regularizes that singularity. Thus the final printed maximum is correct for the integrated variable $c$, or for an “average” using a fixed reference width, but is not generally the maximum of the local mean $c/w(z)$. Both quantities have been given rather than silently identifying them. If $D=0$, the smooth similarity formula is replaced by the advected impulse $V_c\delta (z-at)$; positive eddy diffusivity is essential to the gamma profile.