Open access peer-reviewed chapter

# Simulation of Liquid Flow Permeability for Dendritic Structures during Solidification Process

Written By

S. M. H. Mirbagheri, H. Baiani, M. Barzegari and S. Firoozi

Submitted: November 4th, 2010 Reviewed: March 1st, 2011 Published: July 5th, 2011

DOI: 10.5772/19647

From the Edited Volume

## Computational Fluid Dynamics

Edited by Igor V. Minin and Oleg V. Minin

Chapter metrics overview

View Full Metrics

## 1. Introduction

Alloy solidification is a two stage process, which starts by nucleation and ends by growth of solid phases. Subsequently, number, distribution and morphology (dendritic or non-dendritic) of the grains are formed during the solidification. Some critical defects such as micro/macro segregation, micro/macro porosities and micro/macro shrinkage take place in the solidification stage. The micro-defects are located in the interdendritic space, which are micro-channels that fluid flow through them in the last stage of the solidification. Herein, the region in the grain growth stage is introduced as a mushy zone (or porous media), where the solid phase is constantly progressing; and the ability of fluid to flow into the mushy zone is known as permeability of interdendritic liquid. Therefore, formation of the micro-defects depends on controlling of the permeability factor. In a great number of studies micro/macro solidification models have been simulated based on the permeability factor using Darcy’s law (Ganesan & Poirier, 1990; Nandapurkar et al., 1991; Poirier, 1987; Worster, 1991).

Interdendritic flow, in many CFD documents, is described using Darcy’s law, which relates the fluid flow rate to the pressure gradient, fluid viscosity, and permeability of the porous medium. To obtain an expression for the permeability as a function of the porosity of the porous medium, one generally considers flow through an idealized medium geometry, since it is impractical to solve the flow equations for the complex flow between the particles. Fig. 1 presents two viewpoints for investigation of the permeability in the porous media: metallurgical view and non-metallurgical (or common) view.

As shown in Fig. 1a two of the most commonly used geometries for analytical models are capillaries (Carman, 1937; Chen et al. 1995; Williams et al., 1974) and an array of spaced particles. A more realistic approach will be introduced that assumes geometry of a periodic or random array of cylinders. Since it is not possible to solve analytically for this type of flow over the full range of porosities, two limiting closed form solutions are used for lubrication and point-particle (dilute) models in low and high porosities, respectively. Analysis of permeability for Stokes flow through periodic arrays of cylinders were done by Sangani & Acrivos (1982), Sparrow & Loefler (1959) and, Larson & Higdon (1986).

The effect of fluid inertia on pressure drop required to drive the flow is a function of Reynolds number. Several authors computed the fluid flow through periodic arrays of cylinders as the function of Reynolds numbers for three ranges of low, moderate, and high Reynolds number (Cai & Berdichevsky, 1993; Edwards et al., 1990; Eidsath et al., 1983; Ghaddar, 1995; Koch & Ladd, 1997; McCartney, 1994; Nagelhout et al., 1995; Sangani & Mo, 1994; Sangani & Yao, 1988; Thom & Aplelt, 1961). Particularly noteworthy is the work of Koch & Ladd (1997) for modeling permeability and drag force based on Reynolds number using a Lattice-Boltzmann formulation. Also McCartney (1994) calculated the permeability in the range of low Reynolds numbers up to about 150 by Lattice Gas Cellular Automat (LGCA) model.

In the field of metallurgy, however, the models mentioned above had an important application in the mathematical modeling of flow through arrays of dendrites during the solidification of mushy alloys. Fig. 1b shows a schematic category for modeling of permeability during solidification process. There are many investigations for experimental measuring of permeability by variant methods during solidification process (Apelian et al., 1974; Duncan et al., 1999; Murakami et al., 1984; Murakami et al., 1983; Nielsen et al., 2001; Poirier & Ocansey, 1993; Streat & Weinberg, 1976). However, most of the experimental methods had two goals;

1. an experimental goal is to measure the permeability of an alloy at a constant temperature and to correlate the permeability with microstructures quantities. In this case, it is necessary to obtain a quantitative stereological characterization of the solid-phase morphology. This has been achieved by quenching of the mushy sample during the permeability experiment and subsequent metallography and image analysis, which distinguish between the phases that were solid and liquid during the experiment.

2. Another experimental goal is to measure the permeability and the temperature of an alloy during solidification/re-melting and to correlate the permeability with the solid fraction, which can be estimated from the temperature-time curves.

In this case, the composition of the interdendritic liquid varies throughout the experiment, which complicates the experimental design, the experiment method, and the data analysis (Nielsen et al., 2001).

The mathematical method includes two categories;

1. flow through a network of equiaxed grains,

2. flow through columnar dendrite networks with flows parallel and normal to the primary dendrites.

Piwonka & Flemings (1966), Apelian et al. (1974), Streat & Weinberg (1976), Liu et al. (1989) and Murakami et al. (1984) reported the permeability in equiaxed dendritic structures. In these investigations, however, the microstructural length scales were not reported or they could not be used to estimate permeability in a solidification model designed to calculate macro segregation. Ganesan & Poirier (1990 ) and Ocansey & Poirier (1993) measured and reported the permeability with inverse of the specific area of the solid selected as the length scale in equiaxed microstructures based on Kozeny-Carman model. There are at least two length scales associated with equiaxed dendritic solidification: the secondary dendrite arm spacing (or interdendritic flow) and the grain size (or extra-dendritic flow) (de Groh et al., 1993; Wang et al., 1995). Brown et al. (2002) developed a numerical model for the simulation of 3D flow through equiaxed dendrites of an Al13Cu3Si alloy and the determined the variation in permeability of structure as solidification progressed. The model involved the evolution of an equiaxed dendrite and the application of a CFD program to calculate permeability from Darcy’s law.

Streat & Weinberg (1976), Poirior (1987), Ganesan & Poirier (1990), Nandapurkar et al. (1991) Worster (1991) reported the permeability in columnar dendritic structures for macroscopic scale. However, due to the very complex microstructure of the dendrites, permeability determination still remains a challenge. Indeed, the dendritic columnar region is characterized, first, by a strong anisotropy, which requires knowledge of the different components of the permeability tensor (Murakami et al., 1984; Murakami et al., 1983; Poirier, 1987) and, second, by the non-uniform macroscopic properties such as the liquid volume fraction, which continuously varies from unity in the melt to zero in the solid region. Ganesan et al. (1992) showed that the permeability for a flow parallel to primary dendritic arms is hardly dependent on the secondary dendrite arm spacing. In the columnar configuration, data in different dendritic structures have been summarized by Poirier (1987) and quantitative relationships for parallel and normal permeability have been derived using a regression analysis. In the range of liquid fraction considered (>0.66) the regressions were in good agreement with the classical physical models, but the extrapolation beyond the upper limit failed. Numerical experiments for parallel and normal flow to the primary dendritic arms in columnar structures with a high liquid volume fraction (>0.6) have been found to be more successful (Ganesan et al., 1992). Mirbagheri (2008; Mirbagheri & Khajeh, 2008) measured and simulated interdendritic flow for mushy alloy based on permeability factor and modeled some micro structural factors on the mushy alloys permeability. Bahat et al. (1995) calculated from digitized images, the permeability for flow normal to primary dendrite arm using a Navier-Stokes finite element solver.

As mentioned above, permeability for flow through transverse sections of columnar-dendritic alloy have modeled as flows through array of circles, rhombi, cruciform and schematic dendrites. These models used surface area to volume ratio of the solid for normalizing permeability without predicting the effect of protuberance and radius of dendrite for a fixed fraction of solid as well the effect of angle of dendrite radius via liquid streamlines. In spite of these valuable researches in this field, determination of dendritic structural permeability due to the complex microstructure of the dendrites still remains a challenge in both fields of mathematical and experimental methods.

In this chapter, a numerical model has been introduced for determination of liquid flow permeability through dendritic solid phases during growth. The model includes three stages; first, numerical simulation of nucleation and growth of the equiaxed grains using a novel Cellular Automation Finite Difference (CAFD) method, and second, numerical simulation of micro fluid flow for interdendritic liquid alloys using Computational Fluid Dynamics (CFD) technique, and third, calculation of permeability based on Darcy’s law by the pressure and velocity results of CFD code. Finally determining and modeling the permeability variations versus the cooling rate and the solidification rate during growth. This model can be linked as a module into a commercial macro fluid flow code in order to predict the micro-defects such as the micro porosities, shrinkages or micro-segregations.

## 2. Computation models

In the present work, two separate computation models of the nucleation and grain growth, and interdendritic liquid flow have been developed and coupled for modeling of permeability in mushy alloys.

This is achieved by combining sub-models for each of these processes, i.e., computation of nucleation and grain growth by using the CAFD, and the micro fluid flow by using the CFD model for calculation of the interdendritic permeability. The governing equations are described in details in the next sections.

### 2.1. Solid phase generation code

The solid phase generation code, which has been developed in this investigation is based on the CA and KGT model (Atwood & Lee, 2000; Lee et al., 2001). In other words, the solidification model is non-constrained nucleation and growth. The model comprises of:

1. Stochastic grain density based on local under-cooling, and

2. grain growth based on local thermal and under-cooling (Kurz et al., 2001; Kurz et al., 1986).

However, in the present work, for a binary alloy system (Fig. 2) in the mushy zone, it is assumed that there is no constitutional under-cooling during growth phenomenon. Number, distribution, and morphology of nuclei as well as the growth rate are controlled by only thermal under-cooling, which is produced by Newtonian’s heat transfer in a two dimensional (2D) space. In this condition, at the beginning, all the liquid have the same under- cooling (i.e., the gradient of the under-cooling or temperature is equal to zero). As shown in Fig. 3, once the nucleation takes place, temperature around of the nuclei is raised, because of the liberation of latent heat. Therefore, the temperature gradient around of the grains is negative.

After the solidification, the microstructure of the solidified alloy will consist of fully equiaxed grains and no columnar grains. Therefore, in this investigation, initial and boundary conditions for the simulation of heat transfer and liquid flow are based on the non-constrained solidification.

#### 2.1.1. Governing equations

1. Heat transfer equations

In solidification process, there are two terms of heat transfer and latent heat (Hf). Interaction of these two terms affects the domain of thermal distribution. The heat transfer equation for the mushy zone may be written as:

ρCpTt=(KT)+ρHffstE1
fs=TliqTTliqTsolE2

where ρ,Cpand Kare the density, heat capacity and thermal conductivity, respectively; fsis solid fraction and Tsol, Tliqare solidus and liquidus temperature, respectively. ρHf(∂fs/∂t) is the heat source term, which is a function of temperature in the mushy zone, and is written as follows:

ρHffst=ρHf(fsT)(Tt)E3

By differentiating Eq. 2, and substitution in Eq. 3:

ρHffst=ρHffsTTt=ρHf(1TliqTsol)(Tt)E4

The fraction of solid in the mushy zone is estimated by Eq. 2. The release of latent heat between liquidus and solidus temperature is calculated by substituting Eq. 4 into the second term of Eq. 1. Therefore heat transfer equation is given by:

ρ(CPHf(1TliqTsol))Tt=(KT)E5
ρ[CPeq]Tt=(KT)E6

where CPeqcan be considered as a quasi-specific heat given by:

CPeq=[CPHf(1TliqTsol)]E7

Physical properties of the liquid and solid are assumed to be constant above Tliq and below Tsol, respectively. However, in the mushy zone, coefficients of heat conductivity and thermal capacity are presented as kmu= fLkL+ fSkS, andCpmu=fLCPl+fSCPS(Mirbagheri & Silk, 2007). In a case that the alloy composition has no solidification range (i.e., ΔT0 = TliqTsol = 0), such as eutectic composition, the right hand side in Eq. 4 approaches infinity and as a result, a virtual solidification range of 0.1-1 C is assumed.

2. Nucleation equations

When the temperature falls below the liquidus temperature, nucleation begins. In this condition the number of nuclei at each temperature and time are calculated from Eq. 7, as follows:

Ns=Ntotalexp(CmTliq2T(TliqT)2)E8

After determining the number of nuclei in each time step, by assigning a random distribution function, nuclei are distributed in the liquid domain. In the next step, the latent heat of solidification is calculated to adjust the temperature.

3. Growth equations

The final stage of solidification process is the growth simulation. As mentioned before, once a negative temperature gradient is present in liquid adjacent to the grains, equiaxed grains grow as shown in Fig 3. The direction of primary arms of the equiaxed dendrites depend on crystal structure. Here, a B.C.C. crystal structure is assumed, in which each equiaxed dendrite has four perpendicular primary arms in 2D space, where they can grow in 48 crystalline directions (An et al., 2000). In order to simulate the morphology of the growth, a simple shape was used for grains based on Eq. 8 in polar coordinates.

r=Rd+PdCos(4(θ+θ0))if0PdRd1;(R*)hR*=Rd+PdE9

Fig. 4 shows a “cloverleaf” morphology for a dendrite section created based on Eq. 8, where Pdis perturbation, Rdradius of spherical nuclei prior to perturbation, and θangle between primary arm direction and stream line.

A function (Eq. 9) needs to be defined for the dendrites radius growth rate (dr), which is added to the surface of existing grains in each time step of solidification stage.

rnew=rold+dr;dr=dfs2πrE10

Finally, results of nucleation and growth simulation at each dfs, is used in the CFD code to calculate the permeability in domain.

#### 2.1.2. CA-FD model

1. CA-FD model

The finite difference approximation of heat transfer equation is:

Ti,jn+1=Ti,jnΔt{(UTX+VTY)+(DQX+DQY)}E11

where UTX, DTXL, DTXR and DQX are defined as followed:

UTX=0.5{(1+α)Ui1,j(DTXL)+(1α)Ui,j(DTXR)}E12
DTXL=Ti,jTi1,jΔxDTXR=Ti+1,jTi,jΔxDQX=QXRQXLρCΡΔxE13

The terms of VTY and DQY are found in a similar way to UTX and DQX, respectively. If cells (i, j) and (i+1, j) are liquid:

QXR=kRTi+1,jnTi,jnΔxQXL=klTi,jnTi1,jnΔxE14

KRand KLare thermal conductivity and found through the following equations:

KR1=0.5(1ki,j+1ki+1,j)KL1=0.5(1ki,j+1ki1,j)E15

If cell (i, j, k) contains a mixture of solid and liquid metal, the thermal conductivity of this cell is:

ki,j=FSi,jkSi,j+FLi,jkLi,j,FSi,j+FLi,j=1,FLi,j=Ti,jTsTlTs;TsTTlE16

In the freezing range the specific heat and liquid fraction of the mushy metal is found through the following equation:

CpLS=ΔHfTlTs;Ts(T(TlE18

where Ts, Tlare solidus and liquidus temperatures respectively. In finite difference form, Cpis calculated as follows:

CΡ=CPL;Ti,jn)Ti,jn+1)TE19
CΡ=CPS;Ti,jn+1)Ti,jn)TSE20
CΡ=CΡLTi,jnTLTi,jnTi,jn+1+CΡLSTLTi,jn+1Ti,jnTi,jn+1;Ti,jn)TL)Ti,jn+1)TSE21
CΡ=CΡLSTi,jnTSTi,jnTi,jn+1+CΡSTSTi,jn+1Ti,jnTi,jn+1;TLTi,jn)TS)Ti,jn+1E22

It should be noted that iteration is required here, because not only does Ti,jdepends on Cp, but also Cpdepends on Ti,j.

2. Direction of crystalline growth

For a solidifying cell, as the solid fraction within the cell becomes greater than zero, the local temperature of the particles is obtained using the phase diagram and the under-cooling is calculated accordingly. In each solidifying cell, the change in solid fraction is primarily determined by KGT model (Kurz et al., 1986), which calculates the maximum growth rate based on a given under-cooling at near absolute stability limit. The solid fraction is further corrected by diffusion-controlled growth once the solid fraction reaches a critical value; it can grow into its neighboring liquid cells, providing the cells are under-cooled. Where the solid fraction of cells approaches unity, the cell is considered as fully solidified and the grain growth ends. A captured liquid cell by a growing neighboring cell is assigned the same grain orientation as its growing neighbor (Lee et al., 2001). Fig.5 shows the CA-FD algorithm of heat transfer during nucleation and growth

### 2.2. Inter-dendritic liquid flow code

#### 2.2.1. Governing equations

1. Fluid flow equations

The Navier-Stokes and continuity equations are used to simulate flow of the interdendritic liquid through the network of dendritic solids. The Navier-Stokes equation for the incompressible liquid is given by the following equation (Bahat et al., 1995):

ρLDVDt=P+ρLg+μ2VE23

Also, the mass continuity equation for the incompressible liquid is as follows:

ρt+.(ρV)=0E24

It is often desirable to reproduce large scale physical experiments in scaled-down and more manageable laboratory settings. Information on the flow is contained in the parameters which characterize it, such as dynamic viscosity, velocity and density. If these parameters are combined in a suitable way to yield dimensionless quantities, then these enable one to make the desired statements relating the flows on the large and small scales. By introducing the following dimensionless variables (Griebel et al., 1998):

x*=xL,t*=utL,u*=uu,g*=Lu2g,p*=ppρu2,andRe=ρuLμE25

Navier-Stokes equation can be rewritten as the dimensionless form:

DVDt=P+g+1Re2VE26

Thus by solving Eqs. 16a, 16b and 16c, one can get the, velocity and pressure fields for an incompressible viscous fluid:

ut+px=1Re(2ux2+2uy2)(u2)x(uv)y+gxE27
vt+px=1Re(2vx2+2vy2)(uv)x(v2)y+gxE28
ux+vy=0E29

Equation 15 can be written as two independent equations based on projection model as follows:

(δV)t+V.V=1Re2V;δV=V*VE30
(δV)t+P=0E31

Equation 17a is independent of pressure and by using the Poisson’s equation (Eq. 18) the Navier-Stokes equation is solved implicitly.

2P=1Δt.V*E32

After solution of the Navier-Stokes equations, the pressure and velocity fields can be calculated to obtain the permeability of the mushy zone in each growth sequence.

2. Permeability equations

Permeability is a measure of the ability of a porous material to transmit fluids. Darcy's law is a simple proportional relationship between the instantaneous discharge rate through a porous medium (V), the viscosity of the fluid (μ) and the pressure drop over a given distance. The law was formulated by Henry Darcy based on the results of experiments on the flow of water through beds of sand (Darcy, 1856).

V=KμPE33

where, Kis permeability. The problem with this law is that it is not valid at high Reynolds number. As the velocity increases, inertia effects appear, which complicates the calculation of permeability. To overcome this weakness, several corrections have been applied to the Darcy’s law by various researchers. In this investigation, the effective Darcy’s law equation was used for calculation of permeability (Eq. 19).

V=KμflPE34

#### 2.2.2. Numerical solution of the fluid flow governing equations

The main purpose of the permeability simulation in the solidification process is calculating the velocity profile for inter-dendritic space during equiaxed grains growth. Therefore, governing equations are solved by FDM. The numerical solution method can be considered in four steps:

1. Meshing of the system,

2. Converting differential equations to finite difference approximation,

3. Solution of the finite difference approximations of momentum in order to calculate velocity profile, and

4. Calculation of permeability of the mushy zone by adding Darcy’s law in each growth sequences.

1. Meshing

The computational domain is divided into a number of cells with ΔX, ΔYdimensions, and all cell dimensions are equal for all calculations. Curved boundaries of system are approximated by stepped boundaries. For example, if the solid phase fraction occupies more than 0.9 volume of a cell, the cell is considered to be a complete solid cell. As shown in Fig. 6, domain is discretized as a staggered grid, i.e. the scalar and the vector variables are located at the centre and sides of the computational cell, respectively.

The diffusive terms of the momentum equation is discretized using central differencing. The discretization of the convective terms of momentum equation has been done using the donor-cell scheme (Mirbagheri et al., 2003)

2. Finite difference approximations

The finite difference approximation of momentum equations are:

Ui,jn+1=Ui,jn+Δt[gxFUXFUY+1Re(VISX)]+Δt(Pi,jn+1Pi+1,jn+1Δx)Ui,jn+1=Axn+Δt(Pi,jn+1Pi+1,jn+1Δx)E35
Vi,jn+1=Vi,jn+Δt[gYFVXFVY+1Re(VISY)]+Δt(Pi,jn+1Pi,j+1n+1Δy)Vi,jn+1=Ayn+Δt(Pi,jn+1Pi,j+1n+1Δy)E36
Axn=Ui,jn+Δt[gxFUXFUY+1Re(VISX)]Ayn=Vi,jn+Δt[gYFVXFVY+1Re(VISY)]E37
FUX=Ui,jn2Δx[DUL+DUR+αSign(Uα)(DULDUR)];DUL=Ui,jnUi1,jn;DUR=Ui+1,jnUi,jn;Uα=(DUR+DUL)/2;VISX=(Ui+1,jn2Ui,jn+Ui1,jnΔx2+Ui,j+1n2Ui,jn+Ui,j1nΔy2)E38

Other terms of the flux and the viscosity in x and y directions such as FUY, FVX, FVY, and VISY are obtained in a similar way as FUX and VISX terms, respectively. The finite difference approximation of Poisson’s equation (Eq. 23) is

pi+1,jn+12pi,jn+1+pi1,jn+1(δx)2+pi,j+1n+12pi,jn+1+pi,j1n+1(δy)2=1δt(Ax,i,jnAx,i1,jnδx+Ay,i,jnAy,i,j1nδy)E39
3. Solving and computing procedures

Various methods are available to adjust the pressure term in the Navier-Stokes equation (Hong, 2004; Versteeg & Malalasekra, 1995). Methods like SMAC (Amsden & Harlow, 1970) and SOLA (Hirt et al., 1975) use an explicit scheme to solve the pressure adjustment equation, which is based on divergence of pressure in each computational cell. The computational time are relatively long in these methods and are seldom used today due their low performance. Instead, other methods such as projection method are used to adjust the pressure term and compute new velocities that satisfy the continuity equation. Semi-implicit methods such as SIMPLE, SIMPLER and PISO are also frequently used (Versteeg & Malalasekra, 1995).

To solve the governing equations of fluid flow in a mushy zone for binary alloys, a new CFD code has been developed based on fundamentals presented by (Griebel, 2011). Therefore, to correct the pressures calculated from the Navier-Stokes momentum equation, the Chorin’s projection methodis used, and the resulting Poisson’s equation is solved using Successive Over-Relaxation (SOR) iterative method.

Projection method uses an auxiliary velocity V* to obtain a Poisson’s equation for pressure. This equation can be solved using any solution algorithm such as the SOR or Gauss-Seidel methods. The momentum and continuity equations in time derivative form are:

Vn+1VnΔt+Vn.Vn+Pn+1=1Re2VnE40
Vn+1=0E41

The superscripts (n) and (n+1) denote old and new time level, respectively. In projection method Vn+1domain is calculated at each new time step. Using the auxiliary velocity, the momentum equation can be split to two independent equations with Eq. 26 with no pressure term and Eq. 27 with pressure term.

V*VnΔt+Vn.Vn=1Re2VnE42
Vn+1V*Δt+Pn+1=0E43

The divergence of Eq. 26 takes the form:

.Vn+1.V*Δt+2Pn+1=0E44

Continuity equation (Eq. 25), requires that Vn+1to be zero, thus

2Pn+1=1Δt.V*E45

The overall solution process is using Eq. 26 to obtain V*, then to solving the Poisson’s equation (Eq. 29) to find the adjusted pressure values and finally solving Eq. 27 to obtain Vn+1. Fig. 7 describes the algorithm used in the present code.

4. Calculation of permeability

After obtaining the pressure and velocity fields, which is called herein as the temporary permeability subroutine, and using Darcy’s law, coefficient of temporary permeability is calculated for domain in each solid fraction and their changes are saved until the end of solidification. This subroutine is showed at the lower part of the flowchart in Fig 7.

### 2.3. Validation of the present CFD code

The developed CFD code was validated by comparing the pressure gradient for an equiaxed grain with the Fluent code at the same initial and boundary condition. Fig. 8 shows four pressure fields adjacent to a single dendrite at solid fractions of 0.02, 0.08, 0.19, and 0.34, for and Al-6%wt.Si alloy. Fig. 9 shows comparison of the pressure gradient results between the present CFD code and the Fluent code based on data of Table 1. The CFD code predictions at high and low solid fractions are in good agreement with the corresponding Fluent code. Therefore, the present code could be valid for simulation of other domains.

## 3. Temporary permeability and effective parameters

It is clear that the permeability is a function of pressure gradient based on Darcy’s law; therefore, the simulation of the pressure field in the mushy zone and the factors that affect it is our target. These factors such as, cooling rate, nucleation and growth rate, distribution of nucleus and grains, and morphology of grains, subsequently affects the permeability.

Since, in this investigation a Newtonian thermal condition is assumed; nucleation phenomena are entirely random. Therefore, the distribution of nuclei location could affect the behavior of the pressure field due to the drag force of nucleus on fluid flow. Figs. 10 to 12 show the effect of distribution of nuclei and grains on the pressure field for three constant solid fractions, 0.02, 0.10, 0.24, and 0.43 at a fixed cooling rate. Stream lines are substantially different in these figures, especial at high solid fractions (i.e. fs=0.43), which are affected from the number and size of the grains as well as the distribution of nuclei and grains.

Fig 13 shows the effect of cooling rate on the pressure field for a fixed solid fraction (fs=0.43%). Results show that by increasing the rate of heat extraction from the domain (dQ), some parameter such as size, numbers, and distribution of nucleus and grains have changed and subsequently the pressure field has changed. For example at dQ= 70, grains are coarse and large, however at dQ= 300, grains are fine and small.

Fig 14 shows the temporary permeability that was calculated from the simulated pressure fields. In low solid fraction (fs< 0.06) and high solid fraction (fs>0.78), the behavior of permeability is as asymptotic function, because of flfactor in Eq. 19a. However between two solid fractions, behavior of the temporary permeability is almost liner. Markers, ■, ●,▲, and ◊, demonstrate the effect of different grains distributions on the permeability coefficient, which do not show a significant effect, especially at low and high solid fractions.

## 4. Conclusions

It is possible, especially in commercial codes, to obtain the thermal history in casting or solidification codes. In fact all numerical fluid flow and heat transfer software can save the cooling rate(Tt), temperature gradient (G=Qk=Tx), and solidification rate (v=fst) for each location of the meshed domain. Therefore, in the present code if the interdendritic permeability is defined as a function of cooling rate, the micro-structure will be related to the thermal history in a macro scale, and the formation of micro defects could be predicted. In other words, the present code is capable of predicting the micro-defects based on thermal history. However, the cooling rate is not an independent variable and should be determined by the temperature gradient and the solidification rate, as follows:

Tt=(Tx)(xt)= G.v, for a single phase, andE46
Tt=(k.G)/(1+ΔHm.v), during phase changeE47

It means that increasing the cooling rate results in an increase in the number of nuclei and the change of dendrite morphology from faceted surface to dendritic surface (cloverleaf grain) and that in turn affects the interdendritic permeability. Therefore, the effect of temporary solid fraction (dfs) and heat extraction rate (dQ) on the temporary permeability is modeled as an asymptotic function and plotted in Fig. 15 and 16. Results shown in Fig. 16 show that after formation of 0.72% solid, grains are inter-connected and join together. Consequently, there is no fluid stream into the mushy zone, but only some micro fluid flow as a vortex of the entrapped liquid between the grains. It seems that the critical permeability for the formation of micro-porosity, in a mushy alloy such as Al-6%wt.Si, is less than 1E-10 to 1E-11 (m2). This is true, even at very high cooling rates of over 3000 J/s and different distributions of the equiaxed grains.

In this investigation, the temporary permeability results, which was obtained based on the cooling and solidification rate can be modeled as a 3-dimensional surface. The resulting equation of this 3D-surface as given in Eq. 30 can be used in all casting simulation software.

Log(Permeability)=7.175fs+0.002dQ+5.185E48

Finally, the significant findings of this investigation can be phrased as follows:

1. In this investigation an algorithm was developed to calculate the permeability of mushy alloys. To simulate the permeability, CFD and CA-FD codes were coupled. The CA-FD code was used for nucleation and growth of the equiaxed dendrites and the CFD code was used for simulation of the pressure and velocity fields adjacent to the nuclei and grains.

2. The temporary permeability was calculated based on the present codes (CFD and CA-FD)

3. Permeability simulation results showed that apart from low and high solid fractions, fS<0.06 and fS>0.72, respectively, the temporary permeability versus the solid fraction during solidification has a linier behavior.

4. Results showed that the temporary permeability during nucleation and growth of the equiaxed grains depends on cooling and solidification rate.

5. At high cooling and solidification rate fluctuations on temporary permeability was observed, which was presumably due to the perturbation of the grain surfaces or the evolution of semi-spherical grains to cloverleaf dendrites.

6. Permeability results of present code can be utilized in commercial casting and solidification software.

## References

1. 1. AmsdenA. A.HarlowF. H.1970The SMAC method: A Numerical Technique for Calculating Incompressible Fluid Flows.Los Alamos Scientific Laboratory Report LA-4370New Mexico.
2. 2. AnS. U.LarionovV.MonastyrskiV.MonastyrskaiaE.GrafasI.OhJ. M.LimO. D.KimS. H.LeeJ. H.SeoD. Y.2000The thermal analysis of the mushy zone and grain structure changes during directional solidification of superalloys. In:Superalloys 2000. Green, K. A., Pollock, T. M. & Kissinger, R. D., TMS.
3. 3. ApelianD.FlemingsM. C.MehrabianR.1974Specific permeability of partially solidified dendritic networks of Al-Si alloys.Metallurgical and Materials Transactions B.525332537.
4. 4. AtwoodR. C.LeeP. D.2000A Combined cellular automaton and diffusion model for the prediction of porosity formation during solidification.Modelling of casting welding and advanced solidification processing IX, Aachen, Germany.
5. 5. BahatM. S.PoirierD. R.HeinrichJ. C.1995Permeability for cross flow through columnar-dendritic alloys.Metallurgical and Materials Transactions B.265):10491056.
6. 6. BrownS. G. R.SpittleJ. A.Walden-BevanR.2002Numerical Determination of Liquid Flow Permeabilities for Equiaxed Dendritic Structures.Acta Materialia.5015591569.
7. 7. CaiZ.BerdichevskyA. L.1993Numerical simulation on the permeability variations of a fiber assembly.Polymer Composites.146):529539.
8. 8. CarmanP. C.1937Fluid flow through a granular bed.Transactions of the Institution of Chemical Engineers.15150166.
9. 9. ChenY. T.DavisH. T.MacoskoC. W.1995Wetting of fiber mats for composites manufacturing: I. Visualization experiments.4110):22612273.
10. 10. DarcyH.1856Les Fontaines Publiques De Ville de Dijon, Paris.
11. 11. de GrohH. C.WeidmanP. D.ZakhemR.AhujaS.BeckermannC.1993Calculation of dendrite settling velocities using a porous envelopeMetallurgical and Materials Transactions B.245:749753.
12. 12. DuncanA. J.HanQ.ViswanathanS.1999Measurement of liquid permeability in the mushy zones of aluminum-copper alloysMetallurgical and Materials Transactions B.304):745750.
13. 13. EdwardsD. A.ShapiroM.Bar-YosephP.ShapiraM.1990The influence of Reynolds number upon the apparent permeability of spatially periodic arrays of cylinders.Physics of Fluids A.24560.
14. 14. EidsathA.R. G.C.WhitakerS.HerrmannL. R.1983Dispersion in pulsed systems- III : Comparison between theory and experiments for packed beds.Chemical Engineering Science.3818031816.
15. 15. GanesanS.ChanC. L.PoirierD. R.1992Permeability for flow parallel to primary dendrite arms.Materials Science and Engineering A.15197105.
16. 16. GanesanS.PoirierD. R.1990Metallurgical and Materials Transactions.21B):173181.
17. 17. GanesanS.PoirierD. R.1990Conservation of Mass and Momentum for the Flow of Interdendritic Liquid during Solidification.Metallurgical and Materials Transactions B.21173181.
18. 18. GhaddarC. K.1995On the permeability of unidirectional fibrous media: A parallel computational approach.Physics of FluidsA.725632585.
19. 19. GriebelM.2011Retrieved 2011, fromhttp://www.ins.uni-bonn.de/institut/.
20. 20. HirtC. W.NicholsB. D.RomeroN. C.1975SOLA- A Numerical Solution Algorithm for Transient Fluid Flows.Los Alamos Scientific Laboratory Report LA-5852New Mexico.
21. 21. HongC. P.2004Computer modelling of heat and fluid flow in materials processing, Taylor & Francis.
22. 22. KochD. L.LaddA. J. C.1997Moderate Reynolds number flows through periodic and random arrays of aligned cylinders.J. Fluid Mech.3493166.
23. 23. KurzW.BezenconC.GäumannM.2001Columnar to equiaxed transition in solidification processing.Science Technology and Advance Materials.21):185191.
24. 24. KurzW.GiovanolaB.TrivediR.1986Theory of Microstructrual Development during Rapid Solidification.Acta Metallurgica.345):823
25. 25. LarsonR. E.HigdonJ. J. L.1986Microscopic flow near the surface of two-dimensional porous media, Part 1. Axial flow.Journal of Fluid Mechanics.166449472.
26. 26. LeeP. D.AtwoodR. C.DashwoodR. J.NagaumiH.2001Modeling of Porosity Formation in Direct Chill Cast Aluminium-Magnesium Alloys.Materials Science and Engineering. Vol.: pp. Submitted.
27. 27. LeeP. D.ChiraziA.SeeD.2001Modeling Microporosity in Aluminum-Silicon Alloys: a Review.Journal of Light Metals.11):15
28. 28. LiuC. Y.MurakamiK.OkamotoT.1989Permeability of dendrite network of cubic alloys.Materials Science and technology.511):11481152.
29. 29. Mc CartneyJ. F.1994Flow through arrays of cylinders: Lattice gas cellular automata simulations.Physics of Fluids.62):435437.
30. 30. MirbagheriS. M. H.2008Modeling of the equiaxed dendrite coarsening based on the interdendritic liquid permeability during alloy solidification.Metallurgical and Materials Transactions B-Process Metallurgy and Materials Processing Science.393):469483.
31. 31. MirbagheriS. M. H.AshuriH.VarahramN.DavamiP.2003Simulation of mould filling in lost foam casting process.International Journal of Cast Metals Research.166):554565.
32. 32. MirbagheriS. M. H.KhajehE.2008Modelling and simulation of equiaxed dendritic structures permeability for Pb-Sn alloys.Materials Science and Technology.2412):14441451.
33. 33. MirbagheriS. M. H.SilkJ.2007Simulation of Si concentration effect on the permeability for columnar dendrite structures during solidification of Al-Si alloy.Materials & Design.281):356361.
34. 34. MurakamiK.ShiraishiA.OkamotoT.1983Interdendritic fluid flow normal to primary dendrite-arms in cubic alloys.Acta Metallurgica.3114171424.
35. 35. MurakamiK.ShiraishiA.OkamotoT.1984Fluid flow in interdendritic space in cubic alloys.Acta Metallurgica.3214231428.
36. 36. NagelhoutD.BhatM. S.HeinrichaJ. C.PoirierD. R.1995Permeability for flow normal to a sparse array of fibres.Materials Science and Engineering A.1911-2):203208.
37. 37. NandapurkarP.PoirierD. R.HeinrichJ. C.1991Momentum equation for dendritic solidification.Numerical Heat Transfer.19A297311.
38. 38. NielsenO.MoA.ApplolaireB.CombeauH.2001Measurements and modeling of the microstructural morphology during equiaxed solidification of Al-Cu alloys.Metallurgical and Materials Transactions A.328):1204912060.
39. 39. PiwonkaT. S.FlemingsM. C.1966Part VIII- Pore Formation in Solidification.Transactions of the Metallurgical Society of AIME.23611571165.
40. 40. PoirierD. R.1987Permeability for Flow of interdendritic liquid in columnar-dendritic alloys.Metallurgical and Materials Transactions B.18245255.
41. 41. PoirierD. R.OcanseyP.1993Permeability for flow of liquid through equiaxial mushy zones.Materials Sciens and Engineering A.171231240.
42. 42. SanganiA. S.AcrivosA.1982Slow flow past periodic arrays of cylinders with application to heat transfer.International Journal of Multiphase Flow.8193206.
43. 43. SanganiA. S.MoG.1994Inclusion of lubrication forces in dynamic simulations.Physics of Fluids.616531662.
44. 44. SanganiA. S.YaoA.1988Transport processes in random arrays of cylinders. II. Viscous flow.Physics of Fluids.319):24352442.
45. 45. SparrowE. M.LoeflerA. L.1959Longitudinal laminar flow between cylinders arranged in regular array.AIChE Journal.5325330.
46. 46. StreatN.WeinbergF.1976Interdendritic fluid flow in a lead-tin alloy.Metallurgical and Materials Transactions B.73):417423.
47. 47. ThomA.ApleltC. J.1961Field Computations in Engineering and PhysicsVan Nostrand, London.
48. 48. VersteegH.MalalasekraW.1995An Introduction to Computational Fluid Dynamics: The Finite Volume Method Approach,Longman Group Ltd., Harlow.
49. 49. WangC. Y.AhujaS.BeckermanaC.de GrohH. C.1995Multiparticle interfacial drag in equiaxed solidification.Metallurgical and Materials Transactions B.261):111119.
50. 50. WilliamsJ. G.MorrisC. E. M.EnnisB. C.1974Liquid flow through aligned fiber beds.Polymer Engineering & Science.146):413419.
51. 51. WorsterM. G.1991Natural convection in a mushy layer.Journal of fluid mechanics.224325359.

Written By

S. M. H. Mirbagheri, H. Baiani, M. Barzegari and S. Firoozi

Submitted: November 4th, 2010 Reviewed: March 1st, 2011 Published: July 5th, 2011