Boundary Layer Theory, Drag Coefficients and Frictional Losses

Heat moves by conduction within a material, convection associated with moving fluid, and radiation exchanged through electromagnetic fields. Each mechanism has a distinct constitutive equation. Thermal conductivity is a material property, convection coefficient is a system response, and emissivity characterizes radiative behavior. Calling them all transport coefficients can be useful mathematically, but disguises their different physical meanings. An appropriate heat-flow network starts with a diagram indicating exactly which boundaries touch and which fluid or surface temperatures are known.

A plane wall with no internal generation and constant conductivity has a linear steady temperature profile. Its resistance is thickness divided by conductivity times area. A cylindrical insulation layer does not: the area changes with radius and the resistance becomes logarithmic. For a spherical layer, reciprocal radius appears instead. These are not cosmetic differences; using the wrong geometry can substantially misestimate losses for thick insulation. The same energy rate must cross each surface of a steady series system, even while temperature gradients change from one material to the next.

Transient cooling introduces energy storage. A uniform solid temperature requires that internal conduction equilibrate much faster than surface convection changes the body temperature. The Biot number compares those resistances; it provides a screening check before a lumped-capacitance model is used. If the Biot number is too large, center and surface temperatures diverge and a distributed heat equation is necessary. A measured single exponential may sometimes mimic a complex transient over a short interval, so it is not proof of uniform internal temperature.

Radiative transfer behaves differently from conduction or weak convection because the net exchange scales with the difference of fourth powers of absolute temperatures. At elevated surface temperature it can become an appreciable part of the energy balance. The Stefan–Boltzmann constant has units that require Kelvin, and emissivity alone cannot describe exchange between arbitrary finite surfaces without geometry and view factors. A radiation estimate calibrated for a blackened surface should not be transferred to polished metal without checking actual surface emissivity.

A useful learning experiment is to compare an uninsulated hot pipe with a multilayer-insulated alternative. First calculate conduction resistance per length. Then attach the proper outer convection resistance, using the new outer radius. Consider also radiation to the surroundings if temperatures justify it. Evaluate loss at two insulation thicknesses and check energy conservation through the layers. This procedure exposes the difference between fixed surface-temperature simulations and genuinely coupled wall–environment calculations, and clarifies why an apparently simple heat-loss number can carry meaningful uncertainty.

Unit 3, Lesson 1: Fourier law and plane-wall conduction

Learning objective

Calculate wall temperature gradients and identify where the constant-k assumption matters.

Conceptual model

Fourier law says conductive heat flux follows the negative temperature gradient. It is not Newton cooling: conductivity k is a material property measured in W/(m K), whereas convection coefficient h depends on fluid conditions and geometry.

For a wall with fixed surface temperatures and no internal source, energy conservation makes heat rate uniform along x. Integrating Fourier law gives q=kAΔT/L. The interior temperature is linear even though heat flows continuously.

Check if k changes strongly with temperature before using a constant. Layer interfaces require equal steady heat rate, but they do not demand equal temperature gradients because materials have different k values.

Governing relationship

q = kA(T_hot−T_cold)/L

Engineering interpretation

The relevant phenomena are thermal conductivity; heat flux; steady state. Start with the conservation law, select boundary conditions and state the physical range of the constitutive model.

Case 1: laboratory rig

Characteristic length: 1.950 m. Increase the driving gradient by 20% and evaluate q = kA(T_hot−T_cold)/L. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 2: pilot-scale trial

Characteristic length: 2.080 m. Increase the driving gradient by 30% and evaluate q = kA(T_hot−T_cold)/L. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 3: design review

Characteristic length: 2.210 m. Increase the driving gradient by 40% and evaluate q = kA(T_hot−T_cold)/L. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 4: field investigation

Characteristic length: 2.340 m. Increase the driving gradient by 50% and evaluate q = kA(T_hot−T_cold)/L. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Fully solved numerical illustration

Problem. A wall is 0.12 m thick with k=1.6 W/(m K) and area 2.0 m². The two face temperatures differ by 60 K. Compute heat flow.

Solution with units. Integrate steady one-dimensional Fourier law to get q=kAΔT/L =1.6×2×60/0.12=1600 W. The thermal resistance is L/(kA)=0.0375 K/W, and q=60/0.0375 agrees. If each face is measured independently, their temperatures—not ambient conditions—belong in this expression.

Model audit. A convective film should be added as a resistance only if the boundary temperatures are fluid temperatures rather than known wall-face temperatures.

Questions

1. Explain the difference between the physical driving force and the constitutive coefficient in fourier law and plane-wall conduction, identifying dimensions for every quantity in q = kA(T_hot−T_cold)/L.
2. What observable field or flux would you measure to check the claimed relationship for fourier law and plane-wall conduction? Design a practical before-and-after comparison.
3. Give a limiting case in which the equation for fourier law and plane-wall conduction simplifies, and another in which the equation becomes unreliable. Explain why.

Unit 3, Lesson 2: Cylindrical and composite conduction

Learning objective

Model a hot pipe with insulation and identify the correct area for a surface film.

Conceptual model

Radial conduction is not equivalent to a planar slab with an arbitrarily selected area. In cylindrical geometry the heat-transfer area grows as 2πrL, so a fixed radial heat rate produces a gradient that changes with r.

Integrating radial Fourier law yields R=ln(r_o/r_i)/(2πkL). Composite layers have resistances in series, and external convection adds 1/(h 2πr_oL), based on the exposed outer radius.

Insulation can occasionally increase heat loss on a very small cylinder because it also increases convection area. This critical-radius effect requires a convective boundary; it cannot appear in the simple fixed-surface-temperature cylinder model.

Governing relationship

R_cyl = ln(r₂/r₁)/(2πkL)

Engineering interpretation

The relevant phenomena are logarithmic resistance; multilayer insulation. Start with the conservation law, select boundary conditions and state the physical range of the constitutive model.

Case 1: laboratory rig

Characteristic length: 2.020 m. Increase the driving gradient by 20% and evaluate R_cyl = ln(r₂/r₁)/(2πkL). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 2: pilot-scale trial

Characteristic length: 2.150 m. Increase the driving gradient by 30% and evaluate R_cyl = ln(r₂/r₁)/(2πkL). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 3: design review

Characteristic length: 2.280 m. Increase the driving gradient by 40% and evaluate R_cyl = ln(r₂/r₁)/(2πkL). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 4: field investigation

Characteristic length: 2.410 m. Increase the driving gradient by 50% and evaluate R_cyl = ln(r₂/r₁)/(2πkL). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Fully solved numerical illustration

Problem. A cylinder with inner radius 0.050 m, outer radius 0.100 m, length 2 m and insulation k=0.20 W/(m K) has a 80 K surface-temperature difference.

Solution with units. The logarithmic conduction resistance is R=ln(0.1/0.05)/[2π(0.2)(2)]=0.2758 K/W. Heat rate is q=80/R=290.0 W. Dividing by length gives 145.0 W/m. A planar calculation using a fixed area would not preserve the correct radial geometry.

Model audit. The temperature profile changes with ln(r), not linearly with r, because radial area is not constant.

Questions

1. Explain the difference between the physical driving force and the constitutive coefficient in cylindrical and composite conduction, identifying dimensions for every quantity in R_cyl = ln(r₂/r₁)/(2πkL).
2. What observable field or flux would you measure to check the claimed relationship for cylindrical and composite conduction? Design a practical before-and-after comparison.
3. Give a limiting case in which the equation for cylindrical and composite conduction simplifies, and another in which the equation becomes unreliable. Explain why.

Unit 3, Lesson 3: Convection and transient cooling

Learning objective

Check Bi < 0.1 before assuming the interior temperature is spatially uniform.

Conceptual model

The lumped approximation replaces a spatial PDE by one energy-storage ODE. Under this assumption, ρVc_p dT/dt=−hA(T−T∞). Temperature approaches ambient exponentially but never crosses it without another energy source.

A small Biot number hL_c/k, often judged against 0.1 as a rough screening limit, supports the assumption that internal resistance is small compared with surface resistance. The time constant τ=ρVc_p/(hA) controls thermal response.

Fitting an exponential to measured temperature does not by itself validate constant properties, negligible radiation or a uniform body temperature. A thermometer at the center and another at the surface can reveal an internal gradient.

Governing relationship

T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]

Engineering interpretation

The relevant phenomena are film coefficient; Biot number; lumped capacitance. Start with the conservation law, select boundary conditions and state the physical range of the constitutive model.

Case 1: laboratory rig

Characteristic length: 2.090 m. Increase the driving gradient by 20% and evaluate T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 2: pilot-scale trial

Characteristic length: 2.220 m. Increase the driving gradient by 30% and evaluate T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 3: design review

Characteristic length: 2.350 m. Increase the driving gradient by 40% and evaluate T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 4: field investigation

Characteristic length: 2.480 m. Increase the driving gradient by 50% and evaluate T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]. Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Fully solved numerical illustration

Problem. A solid body has mass 27 kg and c_p=900 J/(kg K). Its exposed area is 0.50 m² and h=15 W/(m² K). Find the lumped cooling time constant.

Solution with units. The total thermal capacitance is m c_p=27×900=24,300 J/K. The film conductance is hA=7.50 W/K, yielding τ=3240 s or 54 minutes. After one time constant the excess temperature drops to e⁻¹≈36.8% of its initial excess.

Model audit. This is a candidate model only: measure body conductivity and characteristic dimension to check Bi=hL_c/k before declaring it valid.

Questions

1. Explain the difference between the physical driving force and the constitutive coefficient in convection and transient cooling, identifying dimensions for every quantity in T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)].
2. What observable field or flux would you measure to check the claimed relationship for convection and transient cooling? Design a practical before-and-after comparison.
3. Give a limiting case in which the equation for convection and transient cooling simplifies, and another in which the equation becomes unreliable. Explain why.

Unit 3, Lesson 4: Thermal radiation

Learning objective

Compare radiative and convective losses across a realistic temperature interval.

Conceptual model

Thermal radiation exchanges electromagnetic energy even across a vacuum. A gray diffuse surface in a large isothermal enclosure approximately obeys q=εσA(T_s⁴−T_sur⁴); all temperatures inside the fourth power must be absolute Kelvin.

The sensitivity to temperature is steep, making radiation increasingly important at high temperatures. Emissivity varies by material, finish and wavelength; a value measured under one condition cannot always be carried into another.

For two finite bodies, view factors and radiative resistances matter. Radiation can be parallel with convection, and a high-temperature equipment balance must account for both without double-counting surface areas.

Governing relationship

q = εσA(T_s⁴−T_sur⁴)

Engineering interpretation

The relevant phenomena are emissivity; absolute temperature; view factor. Start with the conservation law, select boundary conditions and state the physical range of the constitutive model.

Case 1: laboratory rig

Characteristic length: 2.160 m. Increase the driving gradient by 20% and evaluate q = εσA(T_s⁴−T_sur⁴). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 2: pilot-scale trial

Characteristic length: 2.290 m. Increase the driving gradient by 30% and evaluate q = εσA(T_s⁴−T_sur⁴). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 3: design review

Characteristic length: 2.420 m. Increase the driving gradient by 40% and evaluate q = εσA(T_s⁴−T_sur⁴). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Case 4: field investigation

Characteristic length: 2.550 m. Increase the driving gradient by 50% and evaluate q = εσA(T_s⁴−T_sur⁴). Establish the control volume, boundary conditions and independent balance. In a linear constant-property regime the output scales with the driving force; outside it recalculate transport coefficients. Identify a limiting case, a validation measurement and an uncertainty source.

Fully solved numerical illustration

Problem. An approximately gray surface of emissivity 0.80 and area 1.0 m² is held at 500 K in a large 300 K enclosure. Calculate net radiative transfer.

Solution with units. Using σ=5.670374419×10⁻⁸ W/(m² K⁴), q=0.8σ(500⁴−300⁴)=2.468×10³ W (rounded). The temperature difference in Celsius must not be inserted into the fourth-power formula; only absolute temperatures are correct.

Model audit. A polished metal can have much lower emissivity than a black coating. Changes in surface finish may dominate the uncertainty.

Questions

1. Explain the difference between the physical driving force and the constitutive coefficient in thermal radiation, identifying dimensions for every quantity in q = εσA(T_s⁴−T_sur⁴).
2. What observable field or flux would you measure to check the claimed relationship for thermal radiation? Design a practical before-and-after comparison.
3. Give a limiting case in which the equation for thermal radiation simplifies, and another in which the equation becomes unreliable. Explain why.

Applied Project 3 — Heat transport

Develop an engineering model. Supply geometry, conservation equations, units, two operating points, simulation plots, uncertainty, and a validation plan.

Detailed project brief

Evaluate insulation on a hot vessel jacket: the inner and outer radii are 50 and 130 mm, jacket length is 1 m, insulation conductivity is 0.08 W/(m K), and the hot/cold interface temperatures are 160 and 35 °C. First compute radial conductive loss from known boundary temperatures. Next add a realistic outer convection resistance, explicitly noting that the outer surface temperature becomes an unknown and therefore the original boundary data cannot be imposed simultaneously. Compare two insulation thicknesses and identify any geometry-related nonlinear behavior. Consider a radiation component only after choosing an emissivity and enclosure temperature. Show why a planar slab estimate is biased for a thick radial layer and summarize the dominant uncertainty.

Unit 3 Assessment — Heat transport

1. Explain the physical principle supporting Fourier law and plane-wall conduction and state all assumptions in q = kA(T_hot−T_cold)/L.
2. Derive a dimensionless or limiting-form check of q = kA(T_hot−T_cold)/L; what happens when the driving gradient tends to zero?
3. Design an experiment to distinguish a boundary-condition error from an incorrect material coefficient for Fourier law and plane-wall conduction.
4. Explain the physical principle supporting Cylindrical and composite conduction and state all assumptions in R_cyl = ln(r₂/r₁)/(2πkL).
5. Derive a dimensionless or limiting-form check of R_cyl = ln(r₂/r₁)/(2πkL); what happens when the driving gradient tends to zero?
6. Design an experiment to distinguish a boundary-condition error from an incorrect material coefficient for Cylindrical and composite conduction.
7. Explain the physical principle supporting Convection and transient cooling and state all assumptions in T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)].
8. Derive a dimensionless or limiting-form check of T(t)−T∞ = (T₀−T∞)exp[−hAt/(ρVcₚ)]; what happens when the driving gradient tends to zero?
9. Design an experiment to distinguish a boundary-condition error from an incorrect material coefficient for Convection and transient cooling.
10. Explain the physical principle supporting Thermal radiation and state all assumptions in q = εσA(T_s⁴−T_sur⁴).
11. Derive a dimensionless or limiting-form check of q = εσA(T_s⁴−T_sur⁴); what happens when the driving gradient tends to zero?
12. Design an experiment to distinguish a boundary-condition error from an incorrect material coefficient for Thermal radiation.

Laboratory manual The laboratory environment is fully local: no tracking, paid API or network dependency. Each model lists assumptions and solves an analytic expression or an explicitly labeled numerical approximation. Graphs show a physical profile and parameter sensitivity. How to investigate Copy the equation and state variable units. Derive a limiting or conservation check. Change only one input, then record both plots and the key numerical output. Export the CSV and independently reproduce a calculation. Identify an assumption that could fail in practice and an experiment to test it. Safety and scope All property and equipment data in this book are illustrative and intended solely for learning. None of the models selects relief equipment, predicts hazardous releases, sizes industrial units or replaces measured properties and qualified process-safety review.

Engineering measurement, uncertainty and model validation

Why equations need measured inputs

An equation can be mathematically exact under its assumptions and still predict a real process poorly if its inputs were measured incorrectly. In transport calculations, viscosity, effective diffusivity, thermal conductivity and boundary conditions can all carry uncertainty. Before comparing model and experiment, list each measured quantity, its units, calibration source and expected accuracy. Keep instrument resolution separate from repeatability: repeated readings may agree closely while every reading shares the same calibration bias.

Case A: pipe radius controls flow

For a pressure-controlled laminar circular pipe, the Hagen–Poiseuille relation has Q proportional to R to the fourth power. Taking logarithmic differentials gives approximately δQ/Q = 4δR/R + δ(Δp)/(Δp) − δμ/μ − δL/L for small signed changes. If radius is uncertain by ±2% while all else is treated as exact, its first-order contribution alone is roughly ±8% in Q. This is why reporting a very precise flow from a loosely measured pipe diameter creates false confidence. In a real laboratory, measure diameter at multiple locations and inspect noncircularity; the physical pipe may not have a single well-defined radius.

Case B: heat through a layered wall

For known surface temperatures, total steady heat rate through resistances in series is q=ΔT/ΣR. If insulation conductivity is temperature dependent, evaluating k at one arbitrary temperature can systematically bias q. A physically appropriate estimate may integrate k(T) or iterate layer temperatures with measured property data. Contact resistances introduce additional drops that are absent from a perfect-contact wall model. Measure surface temperatures on both sides of each contact where possible; otherwise infer the combined thermal resistance and state that it may not be uniquely separable into component properties.

Case C: diffusivity from a release experiment

A one-dimensional Gaussian tracer pulse offers an experimental estimate of effective diffusivity because its variance grows linearly as 2Dt. Suppose measured variance increases from 0.002 m² at 10 s to 0.014 m² at 70 s. The estimated slope is (0.014−0.002)/(70−10)=0.0002 m²/s, hence D=0.0001 m²/s. This conclusion assumes a uniform infinite medium without advection, reaction or unmodeled mixing. If the mean position also shifts, convective transport must be considered. A single measurement of peak concentration cannot reliably identify both injected mass and diffusivity.

Case D: boundary-layer coefficient identification

A reported convective coefficient h is often computed as measured heat rate divided by exposed area and temperature difference. If the actual wall temperature varies across the surface, substituting a single thermocouple value may be inaccurate. In an external-flow heated-plate trial, place sensors at multiple downstream distances. Repeat at varied mean velocity and compare Nusselt trends with a correlation whose Reynolds and Prandtl limits match the equipment. Surface contamination, radiant exchange and conduction through the support can all masquerade as convection. Correcting for these parasitic paths is necessary before interpreting h as a clean fluid transfer coefficient.

Case E: numerical uncertainty versus model uncertainty

A computational solver can exhibit excellent grid convergence while its material-property model is wrong. Suppose a refined model predicts 12.0 W, a still finer mesh predicts 11.9 W, but the only experimentally measured heat flow is 9.5 W. The mesh difference is 0.1 W, yet the model–experiment discrepancy is 2.4 W. Refining the mesh further is unlikely to resolve a missing convection resistance or a conductivity error. Check residuals, boundary settings and energy conservation first, then challenge physical assumptions and measurements separately. State which uncertainty sources were quantified and which remain qualitative.

Practical uncertainty propagation

For a prediction y=f(x₁,…,xₙ) with small independent input uncertainties, one common first-order method uses u_y²≈Σ(∂f/∂x_i)²u_i². For correlated inputs there are additional covariance terms. The approximation may fail with strongly nonlinear responses, broad uncertainties, discontinuous regimes or parameters near physical limits. Monte Carlo sampling can help if plausible input distributions are available, but arbitrary normal distributions must not be treated as evidence. In all uncertainty analyses, specify whether quoted intervals are standard uncertainties, expanded uncertainties or observed ranges.

Engineering reporting checklist

A reproducible transport calculation should make it possible for a separate reader to reconstruct its main result. Record the domain geometry; draw the direction of positive flux; write all governing laws and their units; identify every boundary and initial condition; document material properties and their sources; and state any empirical correlation and validity limits. Include at least one independent limiting-case calculation, such as zero gradient causing zero pure conductive flux or a vanishing reaction rate giving nearly complete catalyst utilization. Report numerical solver settings if applicable, including grid resolution, timestep and residual definition. Finally, distinguish what was directly measured, what was fitted, what was computed and what was inferred. This separation prevents a simplified teaching simulation from being mistaken for a validated industrial design tool.