Certified Lumped Approximations for the Conduction Dunking Problem

TL;DR

This paper develops and validates asymptotic approximations for the heat conduction problem with heterogeneous materials in the small Biot number limit, providing error estimates and improving lumped model criteria.

Contribution

It introduces first- and second-order asymptotic approximations for the conduction problem, with error bounds and a new criterion for domain suitability beyond the standard Biot number.

Findings

Asymptotic error estimates are validated numerically for small Biot number.

Second-order approximation improves accuracy over the standard lumped model.

A new domain criterion identifies when the lumped model is insufficient.

Abstract

We consider the dunking problem: a solid body at uniform temperature $T_\text{i}$ is placed in a environment characterized by farfield temperature $T_\infty$ and time-independent spatially uniform heat transfer coefficient; we permit heterogeneous material composition. The problem is described by a heat equation with Robin boundary conditions. The crucial parameter is the Biot number, a nondimensional heat transfer coefficient; we consider the limit of small Biot number. We introduce first-order and second-order asymptotic approximations (in Biot number) for the spatial domain average temperature as a function of time; the first-order approximation is the standard `lumped model'. We provide asymptotic error estimates for the first-order and second-order approximations for small Biot number, and also, for the first-order approximation, non-asymptotic bounds valid for all Biot number. We also develop a second-order approximation and associated asymptotic error estimate for the normalized difference in the domain average and boundary average temperatures. Companion numerical solutions of the heat equation confirm the effectiveness of the error estimates for small Biot number. The second-order approximation and the first-order and second-order error estimates depend on several functional outputs associated with an elliptic partial differential equation; the latter can be derived from Biot-sensitivity analysis of the heat equation eigenproblem in the limit of small Biot number. Most important is the functional output $\phi$, the only functional output required for the first-order error estimate and also the second-order approximation; $\phi$ admits a simple physical interpretation in terms of conduction length scale. We characterize a class of spatial domains for which the standard lumped-model criterion -- Biot number (based on volume-to-area length scale) small -- is deficient.