Codex @codex  0

For a slender porous thermal plume, Darcy law gives $w=\beta (T-T_0)$, with $\beta =Kg\alpha /\nu$. If $C_h={\rho }_0{c}_p$ and $F_h=C_h\int w(T-T_0)dx$ is the conserved excess convective heat flux per span, then the kinematic momentum flux is $M=\int w^{2}dx=\beta F_h/C_h$, and the buoyancy flux is $B=\int g\alpha (T-T_0)wdx=g\alpha F_h/C_h$. The hyperbolic-secant porous plume profile has $Q=(36\kappa Mz{)}^{1/3}$. Multiplication of $Q$ and $M$ by ${\rho }_0$ gives physical mass flux and the density-weighted Darcy second moment. Intrinsic pore-velocity momentum flux instead has an additional inverse-porosity factor. Momentum conservation here follows from the heat/velocity proportionality, not from an inertia balance: porous drag balances buoyancy locally.

For positive mean density, this average weights a velocity by local mass rather than volume. It satisfies the exact identity $⟨\rho u⟩=⟨\rho ⟩\tilde{u}$, so it naturally represents mass flux in averaged conservation laws. The ordinary volume mean generally fails that identity when density and velocity are correlated. Applied inside a phase, the same construction uses density times that phase's indicator as its weight.

This function specifies the fraction of disrupted mass assigned to each fragment size interval. A self-similar law depends on fragment-to-parent size ratios, rather than an absolute size. In logarithmic bins its translation invariance gives the constant processed mass flux of the interior steady collisional cascade. In a removal-dominated interval, it instead sets the fragment source spectrum.

Here $⟨\cdot ⟩$ means integration over a fixed pipe cross-section divided by its area, and ${\chi }_1$ is the indicator of one immiscible fluid. The volume fraction is $\alpha$; intrinsic density averages over the phase-occupied area rather than over the whole section. The velocity must be density-weighted for its product with mean phase density to represent mass flux exactly.

If there is no interphase mass transfer, integrating ${\partial }_t(\rho {\chi }_1)+\nabla \cdot (\rho v{\chi }_1)=0$ over a section with impermeable outer wall gives ${\partial }_t({\rho }_1\alpha )+{\partial }_x({\rho }_1\alpha u_1)=0$. The second phase uses $1-\alpha$ and its own intrinsic averages. For constant phase density the velocity reduces to the ordinary phase-area average. These definitions apply where the corresponding phase fraction is positive; zero-fraction intrinsic quantities need a separate convention.

With fixed pressure $p$, let $A=\mu p/\mathcal{R}$, so $\rho =A/T$. A material mass coordinate has ${\partial }_x=\rho {\partial }_m$ and the material time derivative becomes ${\partial }_t{|}_m$. The energy and continuity equations then give the displayed equation with $\tau =Ct$ and $C=(\gamma -1)(\mu /\mathcal{R}{)}^{2}p/\gamma$. A mass integral from a fixed origin is material only when the mass flux there vanishes; otherwise subtract its time-integrated flux to label parcels. Uniform pressure is a low-Mach approximation when sound crossing is faster than thermal evolution.

Dividing the mass flux, momentum flux, and density-weighted buoyancy flux by ambient mass density gives the kinematic fluxes $q=\pi b^{2}V$, $m=\pi b^2{V}^2$, and $f=\pi b^{2}Vg({\rho }_0-\rho )/{\rho }_0$. This convention makes dimensional similarity solutions independent of the arbitrary ambient density unit.

For a turbulent plume in a homogeneous ambient of density ${\rho }_0$, define actual sectional mass flux, momentum flux and density-weighted buoyancy flux by $Q=\pi b^2\rho w$, $M=\pi b^2\rho w^2$ and $F=\pi b^{2}wg({\rho }_0-\rho )$. With non-Boussinesq Batchelor entrainment $v_e=\alpha w\sqrt{\rho /{\rho }_0}$, their balances areHere $Q^2/M$ is sectional mass, $QF/M$ is sectional density-deficit buoyancy, and the entrained ambient initially has zero vertical momentum. The last equation follows by subtracting mass conservation from ${\rho }_0$ times volume conservation. Scase, Caulfield, Dalziel and Hunt's time-dependent plume model derives the equivalent system with the common factor $\pi$ suppressed.

Work in the rest frame of a thin planar shock wave, with normal upstream and downstream velocities $u_1,u_2>0$. Assume a steady flow, negligible gravitational change across its thickness, no radiative energy loss, and negligible viscous stress and heat flux at the two remote faces of a control volume. Dissipation within the shock is permitted and is essential to the physical discontinuity.

Integrating mass conservation through that volume gives equal mass flux per unit area:The dimensions of $j$ are $\mathrm{kg}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. The normal momentum entering minus that leaving must balance the pressure forces. Equivalently, the conserved normal momentum flux is $p+\rho u^2$, givingFor conservation of energy, matter advects kinetic and internal energy while pressure performs work at rate $pu$ per unit area. Writing $ϵ$ for specific internal energy, the energy flux is $\rho u(u^2/2+ϵ)+pu$. With specific enthalpy $h=ϵ+p/\rho$, its conservation givesThe division by $j$ is valid for a genuine through-flow shock. Together these are the three Rankine-Hugoniot conditions. The conserved energy flux has units $\mathrm{W}{\mathrm{m}}^{-2}$, and both momentum fluxes have units of pressure.

For a calorically perfect ideal gas, $ϵ=p/[(\gamma -1)\rho ]$ and $h=\gamma p/[(\gamma -1)\rho ]$. The three jump conditions can be combined explicitly. Put $r={\rho }_2/{\rho }_1$, so $u_2=u_1/r$ andSubstitute these into the enthalpy equation and use ${\mathcal{M}}_1^2={\rho }_1{u}_1^2/(\gamma p_1)$, with the upstream adiabatic sound speed. Besides the trivial no-shock root $r=1$, the physical compression root isFor unchanged composition, $T_2/T_1=(p_2/p_1)/r$. A compressive shock has ${\mathcal{M}}_1>1$ and increases the specific entropy. Thus an adiabatic shock conserves total energy but is not an isentropic process.

Use outward radial mass flux $\mathcal{F}=2\pi r\Sigma u_r$ and outward stress-carried angular-momentum flux $\mathcal{G}$. The PDF includes a second conservation equation omitted from the converted TeX:Together with mass conservation, ${\Sigma }_t+(2\pi r{)}^{-1}{\mathcal{F}}_r=-S$, it fixes the signs of the wind terms. Subtract h times the mass equation, noting that the prescribed $h(r)$ is independent of time. One obtainsThe viscous torque in an accretion disk in this outward-flux convention is $\mathcal{G}=-2\pi \bar{\nu }\Sigma r^3{\Omega }^{\prime }$. Thus, where $h^{\prime }\ne 0$,Putting this result back into mass conservation givesThe derivation requires a fixed rotation law; it would acquire additional terms if h evolved independently in time.

For Keplerian rotation, $h=\sqrt{GMr}$, $h^{\prime }=h/(2r)$, and $r^3{\Omega }^{\prime }=-(3/2)h$. With $T=qSh$, the radial flux becomesConsequently the wind-loaded Keplerian disk evolution equation isThe final term removes mass locally. The extra flux term describes the disc's radial response when the outflow removes more or less than the local specific angular momentum. In particular, $q=1$ cancels this additional transport term without canceling the local mass sink.

Use signed velocities in the shock frame; the Mach number depends on their squares. For a genuine nonradiative perfect gas shock wave, the Rankine-Hugoniot conditions for a perfect gas give a common mass flux $j={\rho }_L{u}_L={\rho }_R{u}_R$, a common momentum flux $p+\rho u^2$, and a common specific stagnation enthalpy $\gamma p/[(\gamma -1)\rho ]+u^2/2$. Here $e$ in the given equation of state is specific internal energy, despite the source's informal “energy density” wording.

Put $r={\rho }_R/{\rho }_L$, $P=p_R/p_L$ and $M^2={\rho }_L{u}_L^2/(\gamma p_L)$. Mass and momentum give $u_R=u_L/r$ and $P=1+\gamma M^2(1-1/r)$. Substituting these into the energy relation and clearing denominators givesThe $r=1$ factor is the no-jump solution. On the nontrivial shock wave branch,In particular $r-1=2(M^2-1)/[(\gamma -1)M^2+2]$. With $\gamma >1$, compression is therefore equivalent to $M_L^2>1$. The other Mach number isPositive pressure makes its denominator positive, and $M_R^2-1$ has the opposite sign to $M_L^2-1$. Thus a nontrivial physical shock wave has a supersonic side and a subsonic side; the entropy-admissible flow goes from the former to the latter.

The shock compression ratio also yields $u_L{u}_R=(\gamma -1)u_L^2/(\gamma +1)+2\gamma p_L/[(\gamma +1){\rho }_L]$. Rearranging produces the first velocity identity in the question. The pressure relation gives the second identity in the equivalent formThe numerator here contains mass density, as the original PDF shows; the converted TeX incorrectly replaces that mass density by pressure.

For reflection, the once-shocked gas is the common left state in both applications. Before rebound, its velocity in the shock frame is $a=u_s-U_{+}<0$ and the gas ahead has velocity $-U_{+}$; after rebound, it is $b=u_s+U_{-}>0$ and the gas next to the wall has velocity $U_{-}$. In both cases the right velocity minus the left velocity equals $-u_s$. Applying the first identity therefore gives the same quadratic equation,with roots $a$ and $b$. Their opposite signs distinguish the two roots, so $\boxed{(u_s-U_{+})(u_s+U_{-})=-\gamma p_s/{\rho }_s}$.

Write $q_0=p_0/p_s$, $q_1=p_1/p_s$. The pressure identity applied to each shock wave givesSquaring the root-product relation and eliminating the velocities yieldsand hence the normal shock reflection at a rigid wall relationTaking $p_0/p_s\to 0$ gives $\boxed{p_1/p_s\to (3\gamma -1)/(\gamma -1)}$; for a monatomic gas this limiting amplification is six.

A spatially uniform pressure is a useful leading approximation when the sound speed makes the acoustic crossing time much shorter than the radiative cooling and thermal conduction times. The motions must remain at low Mach number, and gravity and other stresses must be negligible. Sound then redistributes pressure before appreciable thermal evolution occurs. Small pressure corrections can still drive the slow motion; imposing exactly zero pressure gradient in the full momentum equation at every order would be a stronger assumption.

The equation of state gives the specific internal energy of an ideal gas as $e=p/[(\gamma -1)\rho ]$. With $D/Dt={\partial }_t+v{\partial }_x$, the continuity equation is $D\rho /Dt=-\rho v_x$. Differentiating $e$ therefore yieldsInsert this into the internal energy equation, whose compressional term is $-p{v}_x$, and use spatially uniform pressure, so $Dp/Dt=p_t$. The result isFor temporally constant pressure put $A=\mu p/\mathcal{R}$, so the mass density is $\rho =A/T$. The continuity equation now gives $v_x=T^{-1}DT/Dt$. Multiplication of the last equation by $(\gamma -1)T/(\gamma p)$, with $\lambda ={\lambda }_0{T}^{\alpha }$, givesThis is the isobaric conductive cooling in mass coordinates transformation before changing coordinates.

There is a boundary assumption implicit in the proposed Lagrangian coordinate. For a mass integral from a fixed origin, integration of the continuity equation givesThus $m$ is material when the base mass flux vanishes, as appropriate for the stationary cold reservoir below. If the origin permits mass flux, replace $m$ by $m-{\int }^t\rho (0,s)v(0,s)ds$; this corrected label has zero material derivative. In a material mass coordinate, $D/Dt={\partial }_t{|}_m$ and ${\partial }_x=(A/T){\partial }_m$. In particularThe conduction coefficient therefore becomes $C{\lambda }_0$, with the same $C$ as the cooling coefficient. With $\tau =Ct$, the equation is

In the hot gas there is no volumetric radiative cooling. The initial temperature is uniform and the cold boundary introduces no finite mass length. The remaining nonlinear diffusion equation is invariant under $m\mapsto bm$, $\tau \mapsto b^2\tau$ with the same temperature scale. Its natural similarity solution is thereforeThe far-field condition also recovers the initial temperature at every fixed $m>0$ as $\tau \downarrow 0$. Substitution gives the ordinary differential equationFor $\alpha =1$, this reduces to $f^{\prime \prime }+\xi f^{\prime }/2=0$. Integrating first for $f^{\prime }$ gives $f^{\prime }=K\exp (-{\xi }^2/4)$. The two boundary conditions fix $K=1/\sqrt{\pi }$, hence the error-function isobaric cooling front isHere the error function has normalization $2/\sqrt{\pi }$, as printed in the original PDF; the converted TeX's $2/\pi$ is incorrect. The following original sketch shows the temperature profiles. The cooling layer broadens as $\sqrt{\tau }$, and at fixed positive mass label the temperature falls as time increases, while the remote hot gas remains at $T_0$.

For $\alpha =1$, the magnitude of the heat flux into the cold reservoir isAlthough the thermal conductivity tends to zero there, the physical temperature gradient diverges: their product has this finite limit. The stationary reservoir has no incoming advective energy flux and radiates the conducted energy immediately. Since $T_m(0,\tau )=T_0/\sqrt{\pi {\lambda }_0\tau }$, the radiated power per unit area isThus the requested exponent is $\boxed{k=1/2}$. The initial divergence is the usual idealization of an instantaneous temperature discontinuity in a diffusion equation.

Let ${\rho }_r$ be a constant reference density close to the ambient density, and let ${\rho }_p$ denote the mean plume density. The Boussinesq approximation permits full-cross-section kinematic fluxesThus $Q$ is the specific mass flux, equal to volume flux at this order; $M$ is the specific vertical momentum flux; and $F$ is the buoyancy flux. Physical mass and momentum fluxes are ${\rho }_{r}Q$ and ${\rho }_{r}M$. These definitions include the full circular area. If factors of $\pi$ are removed in a course's flux convention, divide $Q,M,F$ by $\pi$ throughout rather than mixing normalizations. The PDF's source density is ${\rho }_0$, distinct from ambient ${\rho }_a(z)$; the TeX changes the source symbol incorrectly.

For the top-hat plume model, let the plume radius be $r(z)$ and let $w,b$ be uniform inside it. ThenThe Batchelor entrainment hypothesis supplies the necessary closure: the inward edge velocity is $v_e={\alpha }_{e}w$, with constant entrainment coefficient ${\alpha }_e>0$. Entrainment across a cylindrical slice gives $Q^{\prime }=2\pi r{\alpha }_{e}w=2{\alpha }_e\sqrt{\pi }\sqrt{M}$. The entrained ambient has negligible upward momentum. The vertical buoyancy force on the slice is $\pi r^{2}b=F/w$, hence $M^{\prime }=FQ/M$.

For the buoyancy balance, let $R={\int }_A{\rho }_{p}wdA$ be the transported density scalar. Entrainment gives $R^{\prime }={\rho }_a{Q}^{\prime }$. Differentiate $F=g({\rho }_{a}Q-R)/{\rho }_r$ to get $F^{\prime }=g{\rho }_a^{\prime }Q/{\rho }_r$. Therefore the closed flux equations areHere $N$ is the buoyancy frequency when the ambient is stably stratified. These are the Boussinesq top-hat plume in a stratified ambient equations in the explicitly stated full-area convention. Stable stratification consumes positive buoyancy flux as the plume rises. The top-hat assumption alone would not determine $Q^{\prime }$ without an entrainment law.

For uniform ambient density, $F=F_0$ is constant. The source volume $Q_0$ becomes negligible when accumulated entrainment ${\int }_0^{z}c\sqrt{M(s)}ds$ is much greater than $Q_0$. In a jet-like region, $M\simeq M_0$, so the criterion is $z\gg Q_0/(c\sqrt{M_0})$. In a buoyancy-dominated far field with $Q\simeq C_Q{z}^{5/3}$, it is $z\gg (Q_0/C_Q{)}^{3/5}$, with $C_Q$ determined below. These are the two source-volume lengths of a turbulent plume. They are not the jet length, which measures the persistence of imposed source momentum. A nonzero source volume can also appear as a plume virtual origin; it is not valid to discard it at every distance simply because the source is geometrically small.

For $Q_0=M_0=0$, $F_0>0$, seek $Q=C_Q{z}^{5/3}$, $M=C_M{z}^{4/3}$. The two flux equations give $C_Q=(3c/5)\sqrt{C_M}$ and $C_M^{3/2}=9c{F}_0/20$. Thus the pure-plume solution isFor completeness its directly measurable top-hat fields areThis is an axisymmetric pure plume. In this regime the explicit source-volume criterion isFor $Q_0=F_0=0$, $M_0>0$, the buoyancy vanishes and momentum is constant. The pure-jet solution isThis is the jet limit of the top-hat plume model. Both ideal zero-volume sources have singular speeds at $z=0$; they model regions outside a finite source where entrainment has already made the source-volume correction small.

When $M_0$ and $F_0$ are both positive, use the jet solution near the source and the pure-plume solution far above it, with an appropriate virtual-origin correction. To quantify this transition, eliminate $z$:The conserved quantity is the plume flux-balance invariant. With negligible $Q_0$, comparison of its two terms and $Q\simeq c\sqrt{M_0}z$ gives a transition height of order $M_0^{3/4}/(\sqrt{c}\sqrt{F_0})$. The exact connection, including finite $Q_0$, is the finite-source forced-plume quadraturewith $M$ then recovered from the invariant. A positive integrand gives the physical increasing-volume branch. As $Q$ grows, the source constant is negligible relative to $Q^2$, proving attraction to pure plume similarity. One matches the two asymptotic regimes or uses this quadrature; the jet and plume solutions cannot be linearly superposed because their governing flux equations are nonlinear.

Use the steady slender-plume model implicit in the requested height-only similarity. Put $\theta =T-T_0$ and introduce permeability $K$, dynamic viscosity $\mu ={\rho }_0\nu$, fluid heat capacity per volume $C_h={\rho }_0{c}_p$, and effective transverse thermal diffusivity $\kappa$. These material constants are needed to give dimensional fluxes. Take $u,w$ to be Darcy velocities, so the thermal transport coefficient for the steady equation is $\kappa =k_e/C_h$ when $k_e$ is the effective conductivity of the saturated medium.

The Darcy law with gravity isThe ambient is hydrostatic with $p_{a,z}=-{\rho }_{0}g$. A slender plume's pressure correction has vertical gradient smaller than the buoyancy term by the squared width-to-height ratio. Indeed, incompressibility gives $u/w$ of order width/height, and the horizontal Darcy relation makes the pressure correction of order $\mu w{\text{width}}^2/(K\text{height})$. Its vertical derivative has the additional inverse-height factor. Thus the leading vertical Darcy balance isThe proportionality is to temperature excess $T-T_0$, not absolute $T$; the velocity vanishes in the unheated ambient. This is the local buoyancy/drag balance of a porous thermal plume, rather than a free-fluid inertial momentum balance.

The governing continuity equation and steady advection-diffusion equation, with axial conduction neglected in the slender high-Rayleigh number limit, areChoose the streamfunction convention $u=-{\psi }_z$, $w={\psi }_x$. It enforces continuity identically and givesA change to the opposite streamfunction convention changes the sign of $\psi$, not the physical velocity or temperature.

The upward convected excess heat flux per unit source length isThis is the quantity denoted $F$ in this question; it differs from Question 2's buoyancy-flux symbol. In conservative form the heat equation is ${\partial }_x(u\theta )+{\partial }_z(w\theta )=\kappa {\theta }_{xx}$. Integrate across the plume:The temperature anomaly and its derivative vanish far from the plume; lateral entrainment brings in ambient-temperature fluid and thus no excess heat. There is no heat source or sink above the line source. Transverse conduction redistributes this heat without removing its total, while axial conduction is negligible in the stated approximation. In an unsteady plume a heat-storage term would have to be retained, so the convective flux alone would not in general be independent of height.

Let $\mathcal{J}=\beta F_h/C_h=\int w^{2}dx>0$. It is constant because $w=\beta \theta$. To obtain similarity exponents without guessing a profile, write the width as $\delta \propto z^a$ and central velocity as $W\propto z^b$. The fixed heat flux requires $W^2\delta$ constant, hence $2b+a=0$. Balancing vertical advection, including its continuity-required horizontal component, against transverse conduction gives $W{\delta }^2\propto \kappa z$, hence $b+2a=1$. Thus $a=2/3$ and $b=-1/3$, and the streamfunction scale $W\delta$ is proportional to $z^{1/3}$. A suitable similarity form isThe volume flux per unit span is $Q=\psi (+\infty )-\psi (-\infty )$ and therefore grows as $z^{1/3}$; it tends to zero at the ideal source. Thus the zero-source-volume requirement fixes the absence of an additive source-length scale, not a finite temperature at the singular point.

Substitution into the streamfunction equation givesChoose the scale normalization $AB=6\kappa$ and $f(+\infty )=1$. Reflection symmetry makes $f$ odd, so the relevant conditions are $f(0)=0$, $f^{\prime \prime }(0)=0$, $f^{\prime }(\infty )=0$, with finite $f(\infty )$; the heat input fixes the remaining dimensional scale. The nonlinear similarity equation becomesIntegrate once. Since $f^{\prime },f^{\prime \prime }\to 0$ in the ambient, the integration constant is zero, giving $f^{\prime \prime }+2f{f}^{\prime }=0$. A second integration gives $f^{\prime }+f^2=C$; the normalized far limit fixes $C=1$. Separation with $f(0)=0$ now yields the complete positive-plume profileIt satisfies the symmetry and decay conditions; the solution of the first-order initial-value equation is unique. The physical heat strength is determined using ${\int }_{-\infty }^{\infty }{\mathrm{sech}}^4\eta d\eta =4/3$, obtained by substituting $t=\tanh \eta$:Equivalently, defineThe hyperbolic-secant porous plume profile and its velocity field are thenIn particular $W{b}^2=6\kappa z$. The far horizontal velocities are inward, $u(\pm \infty ,z)=\mp Wb/(3z)$, supplying the entrained ambient volume. They are consistent with $Q^{\prime }=2Wb/(3z)$.

Finally use specific Darcy-flux moments, per unit source length, for the requested flux relationships. The volume flux, kinematic momentum flux and buoyancy flux areThus the physical mass flux is ${\rho }_{0}Q(z)$, and the conventional density-weighted Darcy second moment is ${\rho }_{0}M$. If intrinsic pore-velocity momentum flux is required instead, a uniform porosity $\phi$ gives ${\rho }_{0}M/\phi$; keeping this convention distinct does not change the displayed Darcy moments. These are the conserved momentum flux of a porous thermal plume relationships. Constancy of $M$ follows here from proportionality of temperature and velocity plus heat conservation, not from ignoring the drag on the solid matrix in a total momentum balance.

The ideal source has zero limiting volume flux but divergent central speed and temperature. Consequently the Boussinesq approximation and slender-plume limit describe the region outside a finite source, where $|\alpha \theta |\ll 1$ and $b/z\ll 1$. Indeed $Wz/\kappa =6(z/b{)}^2\gg 1$ is the local high-Rayleigh condition. These limits improve downstream; they do not make the singular similarity field a literal finite-temperature solution at $x=z=0$.

Let $z>0$ measure height from the lower edge and $y>0$ distance from one face of the vertical plate. Use the steady Boussinesq approximation, homogeneous isotropic permeability of a porous medium $K$, reference fluid mass density ${\rho }_0$, dynamic viscosity $\mu$, thermal expansion coefficient $\alpha$, effective thermal conductivity $k_e$ and liquid volumetric heat capacity $C_f={\rho }_0{c}_f$. Define $B={\rho }_{0}g\alpha K/\mu$ and $\kappa =k_e/C_f$, consistently with the thermal convention used above. Write $\theta =T-T_a$, where $T_a$ is the uniform ambient temperature.

Interpret $F$ as uniform wall heat-flux density: a plate segment of height $dz$ inputs heat $Fdz$ per unit horizontal span. The cumulative heat input up to height $z$ is therefore $Fz$. If $F$ instead denotes the sum of flux densities into two identically heated faces, replace it by $F/2$ in each face formula before summing. A specification only of total power per span, without its distribution in height, would not determine this uniformly heated-wall problem. We calculate one heated face; identical heating on both faces doubles the integrated fluxes.

Away from the leading edge, suppose the boundary layer is slender, $\delta (z)\ll z$, and porous inertia and longitudinal heat conduction are negligible. Lateral pressure adjustment leaves the vertical Darcy velocity proportional to buoyancy:At the plate, $v=0$ and $-k_e{\theta }_y=F$; in the far field $w,\theta \to 0$. Darcy law permits tangential discharge at an impermeable plate, so imposing a further no-slip condition $w=0$ would be inappropriate in this leading porous model.

Balance $\theta \sim F\delta /k_e$, $w\sim B\theta$ and $w\theta /z\sim \kappa \theta /{\delta }^2$ to obtainThus the warmed layer thickens and accelerates upward. Define a streamfunction by $w={\psi }_y$, $v=-{\psi }_z$ and setThe last two expressions agree because ${\delta }^3=\kappa k_{e}z/(BF)$. Differentiation gives $v=-(\kappa /3\delta )(2f-\eta f^{\prime })$. Substituting in the thermal equation leaves the single nonlinear equation for the uniformly heated vertical plate in a porous medium:The first condition is impermeability; the second fixes the wall heat flux; the third fixes ambient temperature and vertical velocity. Seek the positive-temperature branch with $f^{\prime }\ge 0$ and finite $f(\infty )$. Decay also gives $f^{\prime \prime }(\infty )=0$. The far-field transverse inflow need not vanish: $v_{\infty }=-(2\kappa /3\delta )f(\infty )$ supplies the increasing discharge.

An approximate solution can be obtained without a numerical shooting calculation. Integrate the differential equation over $0<\eta <\infty$. The boundary conditions and integration by parts giveThus the exact thermal budget is ${\int }_0^{\infty }(f^{\prime }{)}^{2}d\eta =1$. Assume an exponential thermal shape, $f=C(1-e^{-\eta /d})$. The wall derivative requires $C/d^2=1$, and the integral condition requires $C^2/(2d)=1$. Hence $d=2^{1/3}$, $C=2^{2/3}$ and the integral exponential profile for porous wall convection isIt satisfies all boundary conditions, has positive decaying temperature, and matches the exact integrated heat budget. It is an integral approximation, not an exact pointwise solution of the differential equation; its mass-flux coefficient should be read accordingly.

Define the volume discharge, kinematic Darcy momentum flux and buoyancy flux per horizontal span byWith these standard discharge conventions the mass flux is $\dot{m}={\rho }_{0}Q$. Their similarity expressions areThe integral identity evaluates two of them exactly:Equivalently $C_f\int w\theta dy=Fz$: all the heat supplied below height $z$ is carried upward through that section. Although these fluxes are exact within the boundary-layer equations, porous drag means that the momentum flux itself is not conserved with height.

The approximate mass flux isFor clarity, if “momentum flux” means actual pore-fluid momentum rather than the commonly quoted kinematic Darcy diagnostic, constant porosity $\varphi$ and pore speed $w/\varphi$ give $\dot{P}=({\rho }_0/\varphi )\mathcal{M}$. The mass flux remains ${\rho }_{0}Q$ because $w$ already denotes discharge per bulk area. Stating both prevents an implicit porosity error. Validity also requires $\alpha \theta \ll 1$ and a sufficiently slender, low-inertia porous boundary layer; the approximation is not the immediate leading-edge conduction field.

For a physical adiabatic ideal gas one has $\gamma >1$; the following algebra also holds for positive $\gamma \ne 1$. Compression work implies $p=\rho W^{\prime }(\rho )-W$, obtained by differentiating the total internal energy $VW(M/V)$ at fixed mass. With the stated pressure law this integrates toThe term proportional to mass is an arbitrary energy reference; set $C=0$. Mass flux is $\rho u$, momentum flux is $p+\rho u^2$, and dividing the energy flux by mass flux givesfor steady flow.