A Finite Element Approximation of an Optimal Insulation Problem with Convective Heat Transfer

Abstract

A finite element discretization of an optimal insulation problem with convective heat transfer is considered. The model is formulated as a non-smooth, two-variable convex minimization problem. It accounts for the temperature distribution in a thermally conducting body $\hat{I}\text{c◯}\subseteq {\mathbb{R}}^d$, with $d\in \{2,3\}$, and the distribution of a given amount of insulation material on an insulated boundary part $\hat{I}{“}_I\subseteq \partial \hat{I}\text{c◯}$. The surface integral over the insulated boundary $\hat{I}{“}_I$ is approximated by a mass-lumping quadrature that preserves the structure of the continuous setting and, in particular, yields discrete optimality conditions mirroring their continuous counterparts. Well-posedness, stability, and weak convergence of discrete solutions to the continuous ones are established. Furthermore, a block coordinate descent algorithm for the computation of the discrete solutions is formulated and its linear convergence is derived. Under suitable regularity assumptions, uniform $L^{\infty }(\hat{I}{“}_I)$-bounds and $\text{a priori}$ error estimates for both the temperature distribution and the distribution of a given amount of insulation material are obtained. Numerical experiments are carried out that confirm the predicted error decay rates and demonstrate the method in a qualitative three-dimensional test on a realistic spacecraft crew module capsule geometry with idealized reentry-heating Robin data.

26 pages, 3 figures