Adaptive computation for boundary control of radiative heat transfer in glass

This paper is concerned with an optimal boundary control of the cooling down process of glass, an important step in glass manufacturing. Since the computation of the complete radiative heat transfer equations is too complex for optimization purposes, we use simplified approximations of spherical harmonics including a practically relevant frequency bands model. The optimal control problem is considered as a constrained optimization problem. A first-order optimality system is derived and decoupled with the help of a gradient method based on the solution to the adjoint equations. The arising partial differential-algebraic equations of mixed parabolic-elliptic type are numerically solved by a self-adaptive method of lines approach of Rothe type. Adaptive finite elements in space and one-step methods of Rosenbrock-type with variable step sizes in time are applied. We present numerical results for a two-dimensional glass cooling problem.