On the solution of an ill‐posed design solidification problem using minimization techniques in finite‐ and infinite‐dimensional function spaces

This paper provides a comparative study of two alternative methodologies for the solution of an inverse design solidification problem. It is the one-dimensional solidification problem of calculating the boundary heat flux history that achieves a desired freezing front velocity and desired heat fluxes at the freezing front. The front velocity h(t) and flux history qmS(t) on the solid side of the front control the obtained cast structure. As such, the potential applications of the proposed methods to the control of casting processes are enormous. The first technique utilizes a finite-dimensional approximation of the unknown boundary heat flux function q0(t). The second technique uses the adjoint method to calculate in L2 the derivative of the cost functional, ‖Tm – T(h(t), t;q0)‖, that expresses the square error between the calculated T(h(t), t; q0) and the given freezing front temperature Tm. Both steepest descent (SDM) and conjugate gradient methods (CGM) are examined. A front tracking FEM technique is used for the discretization of the state space. A detailed numerical analysis of the space and time discretization of the ‘parameter’ and state spaces, of the effect of the end condition of the adjoint problem and of other parameters in the solution are examined.

[1]  J. Baumeister On the treatment of free boundary problems with the heat equation via optimal control , 1979 .

[2]  Owen Richmond,et al.  Design of Two-Dimensional Stefan Processes with Desired Freezing Front Motions , 1992 .

[3]  Jürgen Sprekels,et al.  Real-time control of the free boundary in a two-phase stefan problem , 1982 .

[4]  Nicholas Zabaras,et al.  A deforming finite element method analysis of inverse stefan problems , 1989 .

[5]  D. Colton,et al.  The Numerical Solution of the Inverse Stefan Problem in Two Space Variables , 1984 .

[6]  S. Mukherjee,et al.  An analysis of inverse heat transfer problems with phase changes using an integral method , 1988 .

[7]  Nicholas Zabaras,et al.  An inverse finite‐element technique to determine the change of phase interface location in two‐dimensional melting problems , 1991 .

[8]  Prediction of the phase change front using the optimization technique , 1985 .

[9]  L. Schumaker Spline Functions: Basic Theory , 1981 .

[10]  P. Knabner Control of stefan problems by means of linear-quadratic defect minimization , 1985 .

[11]  Chung Moon-Hyun,et al.  Phase-change front prediction by measuring the wall temperature on which solidification occurs , 1987 .

[12]  Christian Saguez Simulation and Optimal Control of Free-Boundary Problems , 1982 .

[13]  Boris Rubinsky,et al.  An Inverse Finite Element Method for the Analysis of Stationary Arc Welding Processes , 1986 .

[14]  Boris Rubinsky,et al.  AN INVERSE FINITE-ELEMENT TECHNIQUE TO DETERMINE THE CHANGE OF PHASE INTERFACE LOCATION IN ONE-DIMENSIONAL MELTING PROBLEMS , 1984 .

[15]  R. Fletcher Practical Methods of Optimization , 1988 .

[16]  W. Kurz,et al.  Fundamentals of Solidification , 1990 .

[17]  Nicholas Zabaras,et al.  Inverse finite element techniques for the analysis of solidification processes , 1990 .

[18]  D. Luenberger Optimization by Vector Space Methods , 1968 .

[19]  B. Blackwell,et al.  Inverse Heat Conduction: Ill-Posed Problems , 1985 .

[20]  Nicholas Zabaras,et al.  Moving and deforming finite‐element simulation of two‐dimensional Stefan problems , 1990 .