Ensemble averaged multiphase Eulerian model for columnar/equiaxed solidification of a binary alloy: I. The mathematical model

A new multiphase Eulerian model for columnar and equiaxed dendritic solidification has been developed. The mean conservation equations are derived by means of a statistical phase averaging technique, and the mathematical formulation of the model can be used for both columnar and equiaxed solidification. The model uses three different phases, respectively, the columnar, the equiaxed solid and the liquid. The new set of equations enables us to simulate the columnar-to-equiaxed transition (CET) during the directional solidification of a binary alloy. Owing to the statistical nature of the model, we are able to treat rigorously the coexistence of equiaxed and columnar structures and consequently the CET phenomena. The averaged equations are closed by means of the cell model approximation. This technique can be successfully used to model the various interactions between the liquid and the solid. It may also incorporate the effects of the inhomogeneities of the various scalar fields, e.g. the solute and temperature gradients. An envelope model is used to parametrize the small scales, i.e. the dendrite scale. This leads us to distinguish two types of liquid, namely, the extra-dendritic and the inter-dendritic liquids. In part I the equations are rigorously derived in the purely diffusive case, whilst in part II we will present one-dimensional simulations of Sn?Pb and Al?Cu directional solidification experiments involving CET phenomena. Quasi-steady state CET maps are also computed.

[1]  Doru M. Stefanescu,et al.  Heat transfer-solidification kinetics modeling of solidification of castings , 1990 .

[2]  S. Felicelli,et al.  Simulation of freckles during vertical solidification of binary alloys , 1991 .

[3]  David R Poirier,et al.  The energy and solute conservation equations for dendritic solidification , 1991 .

[4]  A. Bejan Shape and Structure, from Engineering to Nature , 2000 .

[5]  Isao Kataoka,et al.  Local instant formulation of two-phase flow , 1986 .

[6]  C. Gandin,et al.  Boundary layer correlation for dendrite tip growth with fluid flow , 2003 .

[7]  D. R. Poirier,et al.  Conservation of mass and momentum for the flow of interdendritic liquid during solidification , 1990 .

[8]  Wilfried Kurz,et al.  Dendrite growth at the limit of stability: tip radius and spacing , 1981 .

[9]  M. Rappaz,et al.  Solute diffusion model for equiaxed dendritic growth: analytical solution , 1987 .

[10]  Xiaoping Wang,et al.  Measurement of friction under sheet forming conditions , 1996 .

[11]  J. Ni,et al.  A volume-averaged two-phase model for transport phenomena during solidification , 1991 .

[12]  M. Rappaz,et al.  Solute diffusion model for equiaxed dendritic growth , 1987 .

[13]  Martin E. Glicksman,et al.  Dendritic growth into undercooled alloy metals , 1984 .

[14]  D. Drew Mathematical Modeling of Two-Phase Flow , 1983 .

[15]  F. Incropera,et al.  A continuum model for momentum, heat and species transport in binary solid-liquid phase change systems—I. Model formulation , 1987 .

[16]  D. Drew,et al.  Theory of Multicomponent Fluids , 1998 .

[17]  Chaoyang Wang,et al.  A multiphase solute diffusion model for dendritic alloy solidification , 1993, Metallurgical and Materials Transactions A.

[18]  Marcelo de Aquino Martorano,et al.  A solutal interaction mechanism for the columnar-to-equiaxed transition in alloy solidification , 2003 .

[19]  L. Burmeister Convective heat transfer , 1983 .