An adjoint method for the inverse design of solidification processes with natural convection

Abstract

This paper presents a finite element algorithm based on the adjoint method for the design of a certain class of solidification processes. In particular, the paper addresses the design of directional solidification processes for pure materials such that a desired freezing front heat flux and growth velocity are achieved. This is the first time that an infinite-dimensional continuum adjoint formulation is obtained and implemented for the solution of such inverse/design problems with moving boundaries and Boussinesq incompressible flow. The present design problem belongs to a category of inverse problems in which one is looking for the unknown conditions in part of the boundary, while overspecified boundary conditions are supplied in another part of the boundary (here the freezing interface). The solidification design problem is mathematically posed as a whole time-domain optimization problem. The gradient of the cost functional is calculated using the solution of an appropriately defined continuous adjoint problem. The minimization process is realized by the conjugate gradient method via the solutions of the direct, adjoint and sensitivity sub-problems. The proposed methodology is demonstrated with the solidification of an initially superheated liquid aluminum confined in a square mold. The non-uniformity in the casting product in the direction of gravity due to the existence of natural convection in the melt is emphasized. The inverse design problem is then posed as finding the appropriate spatial-temporal variations of the boundary heat flux on the vertical mold walls that can eliminate or reduce the effects of convection on the freezing interface heat fluxes and growth velocity. The numerical example demonstrates the accuracy and convergence of the adjoint formulation. Finally, open related research design problems are discussed. © 1998 John Wiley & Sons, Ltd.