CONTROL OF THE FREEZING INTERFACE MORPHOLOGY IN SOLIDIFICATION PROCESSES IN THE PRESENCE OF NATURAL-CONVECTION

Citation
N. Zabaras et Th. Nguyen, CONTROL OF THE FREEZING INTERFACE MORPHOLOGY IN SOLIDIFICATION PROCESSES IN THE PRESENCE OF NATURAL-CONVECTION, International journal for numerical methods in engineering, 38(9), 1995, pp. 1555-1578
Citations number
26
Categorie Soggetti
Computer Application, Chemistry & Engineering",Engineering,Mathematics
ISSN journal
00295981
Volume
38
Issue
9
Year of publication
1995
Pages
1555 - 1578
Database
ISI
SICI code
0029-5981(1995)38:9<1555:COTFIM>2.0.ZU;2-7
Abstract
This paper presents a methodology for the solution of an inverse solid ification design problem in the presence of natural convection. In par ticular, the boundary heat flux q(0) in the fixed mold wall, partial d erivative Omega(0), is calculated such that a desired freezing front v elocity and shape are obtained. As the front velocity together with th e flux history q(ms) on the solid side of the freezing front play a de terminant role in the obtained cast structure, the potential applicati ons of the proposed methods to the control of casting processes are en ormous. The proposed technique consists of first solving a direct natu ral convection problem of the liquid phase in an a priori known shrink ing cavity, Omega(L)(t), before solving an ill-posed inverse design co nduction problem in the solid phase in an a priori known growing regio n, Omega(S)(t). The direct convection problem is used to evaluate the flux q(ml) in the liquid side of the freezing Front. A front tracking deforming finite element technique is employed. The flux q(ml) can be used together with the Stefan condition to provide the freezing interf ace flux q(ms) in the solid side of the front. As such, two boundary c onditions (flux q(ms) and freezing temperature theta(m)) are especifie d along the (known) freezing interface partial derivative Omega(l)(t). The developed design technique uses the adjoint method to calculate i n L(2) the derivative of the cost functional, \\theta m - theta(x,t;q( 0))\\(2)(L2), that expresses the square error between the calculated t emperature theta(x,t;q(0)) in the solid phase along partial derivative Omega(I)(t) and the given melting temperature. The minimization of th is cost functional is performed by the conjugate gradient method via t he solutions of the direct, sensitivity and adjoint problems. A front tracking finite element technique is employed in this inverse analysis . Finally, an example is presented for the solidification of a superhe ated incompressible liquid aluminium, where the effects of natural con vection in the moving interface shape are controlled with a proper adj ustment of the cooling boundary conditions.