Chemical engineering transport & equilibrium foundations of Transient Lumped Capacitance.
In chemical process engineering, metallurgical heat treatment, polymer extrusion pelletizing, and food blanching, unsteady-state (transient) heat transfer dictates processing cycle times, microstructure evolution, and thermal stress generation. When a hot solid body is abruptly plunged into a cooler agitated fluid medium, heat conducts from the interior of the solid to its surface and is subsequently removed by convective boundary layer dissipation into the surrounding fluid. The governing thermal behavior depends critically on the relative magnitude of internal conduction resistance to external surface convective resistance, quantified by the dimensionless Biot number: Bi = h · L_c / k, where h is the convective heat transfer coefficient, k is solid thermal conductivity, and L_c = V / A_s is the characteristic body dimension. When Bi < 0.1, internal conductive resistance is negligible compared to surface convection, meaning temperature gradients within the solid remain virtually flat (spatially uniform temperature). Under this lumped thermal capacitance condition, the transient energy balance yields pure exponential decay: (T(t) - T_inf) / (T_0 - T_inf) = exp(-t / τ), where thermal time constant τ = (ρ · V · c_p) / (h · A_s) = 1 / k_rate. In this laboratory investigation, a brass test cylinder instrumented with centerline thermocouples is quenched from 85.0 °C into an agitated cooling bath to confirm lumped capacitance dynamics and evaluate the cooling decay rate k_rate.
1. Unsteady-State Heat Conduction Fundamentals
Transient conduction occurs when boundary conditions or internal heat generation vary with time, causing solid temperatures to evolve dynamically. In industrial operations such as metal quenching, batch autoclave sterilization, and catalyst regeneration, predicting transient thermal trajectories is critical.
2. The Biot Number & Lumped Capacitance Formulation
The Biot number compares internal conductive resistance to external surface convective resistance:
\[ Bi = \frac{R_{cond}}{R_{conv}} = \frac{L_c / (k A_s)}{1 / (h A_s)} = \frac{h L_c}{k} \]where $L_c = V / A_s$ is the characteristic body dimension. When $Bi < 0.1$, temperature gradients inside the solid are negligible, and the body can be treated as a single lumped thermal mass at uniform temperature $T(t)$.
3. Derivation of the Exponential Cooling Law
Applying the first law of thermodynamics to a solid body of volume $V$, density $\rho$, and specific heat $c_p$ cooling in a bath at temperature $T_\infty$:
\[ -\dot{E}_{out} = \dot{E}_{stored} \implies -h A_s (T - T_\infty) = \rho V c_p \frac{dT}{dt} \]Defining temperature excess $\theta(t) = T(t) - T_\infty$ and initial excess $\theta_0 = T_0 - T_\infty$:
\[ \frac{d\theta}{\theta} = -\left( \frac{h A_s}{\rho V c_p} \right) dt = -\frac{dt}{\tau} = -k_{rate} dt \]where $\tau = \frac{\rho V c_p}{h A_s}$ is the thermal time constant and $k_{rate} = 1/\tau$ is the cooling decay rate constant. Integrating with $\theta(0) = \theta_0$ yields Newton's cooling law:
\[ \frac{\theta(t)}{\theta_0} = \exp(-k_{rate} t) \implies \ln\left( \frac{\theta(t)}{\theta_0} \right) = -k_{rate} \cdot t \]Plotting $\ln(\theta/\theta_0)$ versus time $t$ yields a straight line with slope $-k_{rate}$, from which convective coefficient $h$ is evaluated:
\[ h = \frac{k_{rate} \rho V c_p}{A_s} \]Viva voce questions & analytical theory
What is the physical criterion for the validity of the lumped thermal capacitance model?
The Biot number must be less than 0.1 (Bi = h L_c / k < 0.1); this ensures internal conduction resistance is less than 10% of external convective resistance, keeping spatial temperature variations within the solid below 5%.
What is the characteristic length L_c for a solid body in lumped capacitance analysis?
Characteristic length is defined as body volume divided by external surface area (L_c = V / A_s); for a sphere L_c = r_o / 3, and for an infinitely long cylinder L_c = r_o / 2.
How does the thermal time constant τ relate to thermal capacitance and thermal resistance?
The time constant is the product of thermal resistance R_th = 1 / (h A_s) and lumped thermal capacitance C_th = ρ V c_p: τ = R_th C_th = (ρ V c_p) / (h A_s), representing the time required to complete 63.2% of the total temperature change.
What happens when Bi > 0.1 (e.g., in low-conductivity polymers or ceramics)?
Lumped capacitance fails because significant internal temperature gradients develop; transient behavior must be solved via Fourier series partial differential equations (Heisler charts / Bessel functions).
How does fluid agitation in the quench bath affect the cooling time constant τ?
Higher agitation increases the convective heat transfer coefficient h, decreasing convective resistance R_th and shortening the time constant τ, accelerating the cooling rate.
Why is brass (or copper) an ideal material for demonstrating lumped capacitance in the laboratory?
Brass has high thermal conductivity (k ≈ 115 W/(m·K)), which keeps the Biot number small (Bi ≈ 0.025) even in strongly agitated water baths, ensuring textbook exponential lumped behavior.
Experimental precautions and process safety
- Ensure the transfer from heater to quench bath occurs in < 0.2 s to prevent convective pre-cooling in ambient air.
- Maintain vigorous bath agitation throughout quenching to keep the convective coefficient h constant.
- Ensure thermocouple leads are insulated to prevent conduction wicking along the thermocouple wire.
- Avoid touching the hot specimen; use automated pneumatic transfer mechanisms.
Theoretical assumptions & boundary condition limits
- The transport and thermodynamic models apply within the specified calibrated operating window (10.0 to 60.0 s).
- Assumes ideal continuum flow and homogeneous phase mixing without unmodeled parasitic wall losses.
- Continuous sensor readings reflect ordinary least-squares regression across successive observations.
- High-temperature or extreme concentration excursions require non-ideal equation-of-state corrections.