Visual walkthrough — Ablative cooling — charring, blowing

Every symbol on this page is earned before it is used. If you have never seen a derivative or a logarithm, you will still be able to follow from line one.

---

Step 1 — The picture we start from: gas sliding past a wall

WHAT. Imagine hot gas flowing horizontally, and a solid wall underneath it. Right at the wall the gas is stuck (glued by friction), so it moves slowly; far above, it moves fast. In between there is a thin sheet where the speed climbs from zero to full. That sheet is the boundary layer (studied in Boundary Layer Theory).

WHY start here. Heat only reaches the wall by crossing this thin sheet. If we understand the sheet, we understand the heating. Nothing else matters yet.

PICTURE. In the figure, height above the wall is the vertical axis $y$. The horizontal arrows show gas speed: short near the wall, long up top. The dashed line marks the top of the boundary layer at height $\delta$.

---

Step 2 — Temperature crosses the same thin sheet

WHAT. The gas up top is hot (call its energy content $H_{gas}$). The wall is cooler ($H_{wall}$). So as we walk from the wall upward, the temperature — really the enthalpy — rises from wall value to gas value across the very same sheet.

WHY enthalpy $H$, not temperature $T$? At thousands of kelvin the gas molecules break apart and recombine, storing energy in chemistry, not just in "how fast they jiggle". Temperature only measures the jiggle. ==Enthalpy $H$ counts all the stored energy per kilogram== — so it is the honest measure of "how much heat is available to dump into the wall." (This is the fourth mistake warned about in the parent note.)

PICTURE. A smooth curve of enthalpy $H$ versus height $y$: low at the wall, rising to $H_{gas}$ at the top. The steepness of that curve right at the wall is what drives heat in.

---

Step 3 — The slope at the wall is the heat flux (why we need a derivative)

WHAT. Heat flows from hot to cold, and it flows faster where the temperature changes more sharply. The sharpness of the enthalpy curve at the wall is measured by its slope there. In symbols we write this slope as ${\frac{dH}{dy}|}_{y=0}$.

WHY a derivative — and why this tool and no other? We have a curve (enthalpy vs. height) and we need one specific number from it: how steep is it at exactly the wall? The derivative is the mathematical machine whose entire job is "give me the steepness of a curve at one chosen point." No algebra or ratio can do that for a bending curve — only the derivative. So it enters here by necessity.

FROM temperature-conduction TO enthalpy-conduction (why $c_p$ appears). The physical law of heat conduction (Fourier's law) is written in temperature: the conductive flux is $\dot{q}=-k\frac{dT}{dy}$ — heat flows down the temperature slope, hence the minus sign. But we chose to work in enthalpy $H$, not $T$. For a gas of constant heat capacity, a change in enthalpy and a change in temperature are locked together by $dH=c_{p}dT$, i.e. $\frac{dT}{dy}=\frac{1}{c_p}\frac{dH}{dy}$. Substituting turns Fourier's law into its enthalpy form: $\dot{q}=-k\frac{dT}{dy}=-\frac{k}{c_p}\frac{dH}{dy}.$ The $c_p$ in the denominator is nothing mysterious — it is just the exchange rate that converts a temperature slope into an enthalpy slope.

PICTURE. Zoom into the wall. Draw the straight line that just kisses the curve at $y=0$ (the tangent). Its tilt is the derivative. A steep tilt = lots of heat pouring in; a gentle tilt = little heat.

---

Step 4 — Turn on the blowing: gas seeps up out of the wall

WHAT. The charring material releases pyrolysis gas that leaves the surface moving straight up, with a small vertical velocity $v_w$. This is a steady creep of gas away from the wall, everywhere across it.

WHY it changes the curve. That rising gas is cool. As it climbs it physically carries the near-wall region upward, stretching the low-enthalpy zone taller. The enthalpy curve gets pushed outward — the same rise from $H_{wall}$ to $H_{gas}$ now happens over a taller distance. A taller rise over which to spread means a gentler slope at the wall.

PICTURE. Two enthalpy curves overlaid: the old "no-blowing" curve (steep at the wall) and the new "blowing" curve (pushed up, gentle at the wall). Small upward arrows all along the wall show the injected gas $v_w$.

---

Step 5 — Write the energy balance in the thin sheet (the 1-D model)

WHAT. Take a razor-thin horizontal slice of gas inside the boundary layer, between height $y$ and $y+dy$. Energy can enter/leave two ways: carried up by the blowing (advection) and conducted by temperature difference (diffusion). In steady state, whatever energy crosses the bottom face of the slice must also cross the top face — nothing piles up.

WHY this is the honest 1-D model. We ignore sideways changes and keep only up–down transport. That is the Couette picture the parent note named: a clean 1-D duel between "gas physically hauling heat up" and "heat leaking down the gradient."

Sign convention, and how it survives into the PDE. Count everything as positive when it moves upward ($+y$). The total upward energy flux at any height is the sum of two upward-counted pieces: $F(y)=\underset{\text{advection, rising gas carries }H\text{ up}}{\underbrace{{\dot{m}}^{\prime \prime }H}}+\underset{\text{conduction (Step 3), points DOWN so it is negative here}}{\underbrace{\left(-\frac{k}{c_p}\frac{dH}{dy}\right)}}.$ The minus sign on the conduction term is exactly Fourier's law from Step 3: because $H$ rises with $y$, conduction actually points downward, so counted in the $+y$ convention it enters as a negative number. Steady state with no source means this total upward flux is the same at every height: $\frac{dF}{dy}=0$. Differentiating, ${\dot{m}}^{\prime \prime }\frac{dH}{dy}-\frac{d}{dy}\left(\frac{k}{c_p}\frac{dH}{dy}\right)=0.$ Moving the conduction term to the other side flips its sign, giving the tidy balance below. So the minus sign was never dropped — it became the two terms sitting on opposite sides of the equals sign.

PICTURE. The slice, with an up-arrow (mass ${\dot{m}}^{\prime \prime }$ carrying enthalpy $H$) and a down-arrow (conduction, minus sign) drawn on it, their net change across the slice cancelling.

---

Step 6 — Solve it: the no-blowing baseline, then why an exponential appears

WHAT. First name the constant group. The conduction coefficient $\frac{k}{c_p}$ is fixed (all of $k$, $c_p$ are constants by our small-blowing ground rule), so call it $\Gamma \equiv \frac{k}{c_p}$ — an enthalpy conductivity with units that make $\Gamma dH/dy$ a heat flux. (Note: this is not the thermal diffusivity $\alpha =k/(\rho c_p)$; those differ by a factor of density $\rho$. We keep $\Gamma =k/c_p$ because it is exactly the coefficient that multiplies the enthalpy slope in the flux, which is what our equation contains.) The balance becomes ${\dot{m}}^{\prime \prime }\frac{dH}{dy}=\Gamma \frac{d^{2}H}{d{y}^2}.$

FIRST the baseline: no blowing. Set ${\dot{m}}^{\prime \prime }=0$. The equation collapses to $\Gamma \frac{d^{2}H}{d{y}^2}=0$, i.e. $\frac{d^{2}H}{d{y}^2}=0$: the curvature is zero everywhere. A curve with zero curvature is a straight line. Integrating $H^{\prime \prime }=0$ twice gives $H=a+by$; applying $H(0)=H_{wall}$ and $H(\delta )=H_{gas}$ fixes it to $H_{\text{no blow}}(y)=H_{wall}+(H_{gas}-H_{wall})\frac{y}{\delta }.$ Its slope is constant, $\frac{dH}{dy}=\frac{H_{gas}-H_{wall}}{\delta }$ — this is where the baseline wall slope we quote in Step 7 comes from, now derived, not asserted.

NOW turn blowing back on — why the exponential? With ${\dot{m}}^{\prime \prime }\ne 0$ the equation says: the slope of $H$ is proportional to the slope-of-the-slope of $H$. The one function whose derivative is a copy of itself is the exponential $e^{(\cdot )}$. So the solution must be built from $e^{{\dot{m}}^{\prime \prime }y/\Gamma }$.

WHY the exponential — this exact tool? We asked: "what shape has its own steepness proportional to itself?" That question only has exponential answers. It is not a guess; it is forced by the equation's structure.

HOW we integrate it (step by step). Let $g=\frac{dH}{dy}$ (the slope itself). Then $\frac{dg}{dy}=\frac{d^{2}H}{d{y}^2}$, and the equation becomes ${\dot{m}}^{\prime \prime }g=\Gamma \frac{dg}{dy}⟹\frac{dg}{g}=\frac{{\dot{m}}^{\prime \prime }}{\Gamma }dy.$ Integrating both sides (the left gives a logarithm, the right a straight line) and exponentiating: $g=\frac{dH}{dy}=C_1{e}^{{\dot{m}}^{\prime \prime }y/\Gamma }.$ Integrate once more in $y$ (an exponential integrates to itself, plus a constant): $H(y)=A+B{e}^{{\dot{m}}^{\prime \prime }y/\Gamma },$ where $A$ and $B$ are two unknown constants left over from the two integrations.

Fix the constants with the two boundary conditions. We know the enthalpy at both ends of the sheet:

• at the wall, $y=0$: $H(0)=H_{wall}$,
• at the top, $y=\delta$: $H(\delta )=H_{gas}$.

Plugging in $y=0$: $H_{wall}=A+B$. Plugging in $y=\delta$: $H_{gas}=A+B{e}^{{\dot{m}}^{\prime \prime }\delta /\Gamma }$. Subtracting the first from the second gives $B(e^{{\dot{m}}^{\prime \prime }\delta /\Gamma }-1)=H_{gas}-H_{wall}$, so $B=\frac{H_{gas}-H_{wall}}{e^{{\dot{m}}^{\prime \prime }\delta /\Gamma }-1},A=H_{wall}-B.$ Substituting $A$ and $B$ back and simplifying gives the tidy normalized profile.

PICTURE. Put the two profiles side by side: the derived straight no-blowing line, and the curved blowing profile that hugs the wall flat before shooting up — exactly the "pushed-out, gentle-at-the-wall" shape blowing predicted in Step 4.

---

Step 7 — Read off the wall slope and complete the logarithm

WHAT. Take the derivative of $H(y)$ at $y=0$ — the wall slope from Step 3 — and compare it to the no-blowing straight-line slope we just derived. Define the blowing parameter that bundles all the constants.

Define $B^{\prime }$ once, cleanly. The exponent group is $\frac{{\dot{m}}^{\prime \prime }\delta }{\Gamma }=\frac{{\dot{m}}^{\prime \prime }\delta c_p}{k}$. Recall from Step 3's baseline that $h_0=\frac{\Gamma }{\delta }\cdot \frac{(H_{gas}-H_{wall})}{(H_{gas}-H_{wall})}=\frac{\Gamma }{\delta }=\frac{k}{c_p\delta }$ (the no-blow coefficient is conduction coefficient over thickness). Therefore $\frac{{\dot{m}}^{\prime \prime }\delta }{\Gamma }=\frac{{\dot{m}}^{\prime \prime }}{\Gamma /\delta }=\frac{{\dot{m}}^{\prime \prime }}{h_0}.$ So the single "how hard we blow" dial is $\boxed{B^{\prime }\equiv \frac{{\dot{m}}^{\prime \prime }}{h_0}=\frac{{\dot{m}}^{\prime \prime }\delta }{\Gamma }=\frac{{\dot{m}}^{\prime \prime }\delta c_p}{k}}$ Note on units: here $h_0=k/(c_p\delta )$ carries units of kg/m²·s (an enthalpy-based coefficient, flux per unit enthalpy), so $B^{\prime }={\dot{m}}^{\prime \prime }/h_0$ is dimensionless — the two forms match with no missing $c_p$. (If one instead writes $h_0$ in the temperature-based convention W/m²·K, the same dimensionless group reads $B^{\prime }={\dot{m}}^{\prime \prime }{c}_p/h_0$; both conventions give the identical number.)

HOW the ratio forms. Differentiating the profile once and setting $y=0$: ${\frac{dH}{dy}|}_0=(H_{gas}-H_{wall})\frac{{\dot{m}}^{\prime \prime }/\Gamma }{e^{{\dot{m}}^{\prime \prime }\delta /\Gamma }-1}=(H_{gas}-H_{wall})\frac{B^{\prime }/\delta }{e^{B^{\prime }}-1}.$ The no-blowing baseline slope (Step 6) was $(H_{gas}-H_{wall})/\delta$. Since $h/h_0$ equals the ratio of these two wall slopes (the gap and $\Gamma$ cancel), $\frac{h}{h_0}=\frac{B^{\prime }/\delta }{(e^{B^{\prime }}-1)/\delta }\cdot \frac{1}{1}=\frac{B^{\prime }}{e^{B^{\prime }}-1}.$

COMPLETE the classic form. The result $\frac{B^{\prime }}{e^{B^{\prime }}-1}$ is the exact Couette answer. The parent note quoted the equivalent film-model form $\frac{\ln (1+B^{\prime })}{B^{\prime }}$. They are the same law written with two different definitions of the blowing parameter: if you name the parameter by the enthalpy ratio across the film, $B_{\text{film}}^{\prime }=e^{B^{\prime }}-1$, then $B^{\prime }=\ln (1+B_{\text{film}}^{\prime })$ and substituting gives $\frac{h}{h_0}=\frac{B^{\prime }}{e^{B^{\prime }}-1}=\frac{\ln (1+B_{\text{film}}^{\prime })}{B_{\text{film}}^{\prime }}.$ Dropping the subscript, this is the boxed formula of the parent note. Either way the logarithm is the fingerprint of the exponential profile: $\ln$ is exactly the question "$e$ to what gives $1+B^{\prime }$?", and that is precisely what inverting the wall-slope ratio asks.

PICTURE. Slopes at the wall for several $B^{\prime }$ values, all starting at $H_{wall}$: bigger $B^{\prime }$ ⇒ visibly flatter start.

Recall Asymptotic sanity checks (do the edges behave?)

• $B^{\prime }\to 0$ (no blowing). For tiny $B^{\prime }$, $\ln (1+B^{\prime })\approx B^{\prime }-\frac{1}{2}{B}^{\prime 2}$, so $\frac{\ln (1+B^{\prime })}{B^{\prime }}\to 1$. Then $h\to h_0$ — with no injected gas nothing changes. ✔ Equivalently the Couette form $\frac{B^{\prime }}{e^{B^{\prime }}-1}\to 1$ as $B^{\prime }\to 0$.
• $B^{\prime }\to \infty$ (blow-off). $\ln (1+B^{\prime })$ grows far slower than $B^{\prime }$, so the ratio $\to 0$: heating shuts off. ✔

---

Step 8 — Every case: $B^{\prime }\to 0$, $B^{\prime }$ large, and $B^{\prime }$ between

WHAT & WHY. A formula you cannot test at its edges is a formula you do not yet trust. We check all regimes.

• No blowing, $B^{\prime }\to 0$. $\frac{\ln (1+B^{\prime })}{B^{\prime }}\to 1$, so $h\to h_0$: with no injected gas, nothing changes. ✔ The blowing curve in Step 6 collapses to the straight baseline line.
• Gentle blowing, $B^{\prime }=2$. $\frac{\ln 3}{2}=\frac{1.0986}{2}=0.549$ — heating cut nearly in half.
• Strong blowing, $B^{\prime }=4$. $\frac{\ln 5}{4}=0.402$ — doubling the injection barely helped: from $0.549$ down to $0.402$. This diminishing return is the log flattening out.
• Blow-off, $B^{\prime }\to \infty$. $\frac{\ln (1+B^{\prime })}{B^{\prime }}\to 0$: the wall slope goes flat, heating shuts off — but you are shedding mass furiously. The parent's "blow-off" limit. (Here the small-blowing assumption is stretched, so read this as a trend.)

PICTURE. The full curve of $\ln (1+B^{\prime })/B^{\prime }$ versus $B^{\prime }$, with the four points marked: it starts at $1$, sags through $0.549$ and $0.402$, and creeps toward $0$ — never negative, never above $1$.

---

The one-picture summary

Everything above, on one canvas: the enthalpy curve bending flatter as blowing turns on (left), which feeds the saturating $\ln (1+B^{\prime })/B^{\prime }$ law (right). Read left-to-right: inject gas → curve pushes out → wall slope drops → $h/h_0$ falls along the log curve.

Recall Feynman retelling — say it in plain words

Hot gas can only touch the wall by crossing a thin sticky sheet of slow gas. The heat that gets in depends on how steeply the temperature climbs across that sheet, right at the wall — that steepness is a derivative. Fourier's law of conduction is written in temperature, $\dot{q}=-kdT/dy$, but because enthalpy and temperature are locked by $dH=c_{p}dT$ we rewrite it as $\dot{q}=-(k/c_p)dH/dy$ — the $c_p$ is just an exchange rate. We package that wall steepness (divided by the fixed top-to-bottom enthalpy gap) into one number, $h$. Now let the burning shield seep cool gas upward out of the wall, at density $\rho$ and speed $v_w$, so mass ${\dot{m}}^{\prime \prime }=\rho v_w$ leaves per second. That rising gas stretches the sheet taller, so the same temperature rise spreads over more distance and the wall slope gets gentle — less heat gets in. Writing the steady balance between "gas hauling heat up" (positive) and "heat conducting back down" (negative in the same up-is-positive convention) puts the two terms on opposite sides of one equation with a consistent sign. With no blowing that equation has zero curvature, so the profile is a straight line and the baseline slope is just $(H_{gas}-H_{wall})/\delta$. Turn blowing on and the equation says the slope is proportional to its own curvature — and the only shape that behaves that way is the exponential. Integrating twice and pinning the two constants with the wall value $H_{wall}$ and the top value $H_{gas}$ gives the exact profile. Reading its wall slope and comparing to the straight baseline gives $B^{\prime }/(e^{B^{\prime }}-1)$, which is the very same law as $\ln (1+B^{\prime })/B^{\prime }$ once you name the blowing parameter by the film enthalpy ratio. The logarithm is just the exponential seen backwards. Its shape tells the whole story: at zero blowing it equals one (nothing changes), and as you blow harder it sags toward zero — but slower and slower, so past a point extra gas is wasted. All of this assumes a gentle, laminar, one-dimensional flow with $k$, $c_p$, $\rho$ and $\delta$ held constant; turbulence and heavy blowing shift the numbers but keep the same story.

Recall Quick self-test

• Why must a derivative appear in the heat-flux expression? ::: heat flux depends on the steepness of the enthalpy curve at one point, and only a derivative measures the slope of a bending curve at a chosen point.
• Where does the $c_p$ in $\dot{q}=(k/c_p)dH/dy$ come from? ::: from converting Fourier's law $\dot{q}=-kdT/dy$ into enthalpy using $dH=c_{p}dT$.
• Why enthalpy $H$ instead of temperature $T$? ::: at high temperature the gas stores energy in chemistry (dissociation), which $T$ misses but $H$ counts.
• What is $h$ in one sentence? ::: the heat flux per unit enthalpy gap — the wall slope packaged into a single coefficient.
• What is the no-blowing profile and why? ::: a straight line $H_{wall}+(H_{gas}-H_{wall})y/\delta$, because with ${\dot{m}}^{\prime \prime }=0$ the balance gives zero curvature $H^{\prime \prime }=0$.
• How are the two integration constants fixed? ::: by the boundary conditions $H(0)=H_{wall}$ and $H(\delta )=H_{gas}$.
• Where does the exponential come from? ::: the 1-D energy balance makes the slope proportional to the curvature, and the exponential is the unique self-derivative function.
• Is $k/c_p$ the thermal diffusivity? ::: no — the thermal diffusivity is $k/(\rho c_p)$; $k/c_p$ is the enthalpy-conduction coefficient that multiplies $dH/dy$.
• Why does the log saturate? ::: $\ln (1+B^{\prime })$ grows ever more slowly, so dividing by $B^{\prime }$ drives the ratio toward zero with diminishing returns.
• Value of $h/h_0$ at $B^{\prime }=2$? ::: $\ln 3/2=0.549$.
• What is assumed throughout? ::: laminar, 1-D Couette-type flow with $\delta$, $k$, $\rho$, $c_p$ held constant (small blowing).