## 1. Introduction

Casting process is desired because it is very versatile, flexible, and economical and happens to be the shortest and fastest way to transform raw material into finished product. Squeeze casting belongs to permanent mould casting method which offers considerable saving in cost for large production quantities when the size of the casting is not large.Squeeze casting has the advantage of producing good surface finish, close dimensional tolerance and the absence of sand inclusions on the cast surfaces of the products as opined by Das and Chatterjee, (1981).

The solidification process of the molten aluminium metal in the steel mould takes a complex form, (Hosford and Caddell, 1993) and (Potter and Easterling, 1993). During solidification all mechanisms of heat transfer are involved and the solidifying metal undergoes state and phase changes. The final structure and properties of the cast product obtained depend on the casting parameters applied i.e. applied pressures, die pre-heat temperature, delay time and period of applied pressure on the solidifying metal, (Potter and Easterling, 1993), (Bolton, 1989) and (Callister, 1997). The prediction of temperature distribution and solidification rate in metal casting is very important in modern foundry technologies. This helps to control the fundamental parameters such as the occurrence of defects, as well as, the influence on final properties of cast products and the mould wall / cast metal interface contact surface.

Heat transfer coefficients during squeeze cast of commercial aluminium were determined using the solidification temperature versus time curves obtained for varying applied pressures during squeeze casting process. The steel mould / cast aluminium metal interface temperatures versus times curve obtained through polynomial curves fitting and extrapolation was compared with the numerically obtained temperatures versus times curve. Interfacial heat transfer coefficients were determined experimentally from measured values of heating and cooling temperatures of steel mould and cast metal and compared with the numerically obtained values and found to be fairly close in values.

Aluminium is a product with unique properties, making it a natural partner for the building and other manufacturing industries.The commercially pure aluminium metal used for this research work finds extensive use in the building, manufacturing and process industries, both as a material of construction and household goods. Products of squeeze casting are of improved mechanical properties and could be given heat treatment. Heat dissipation from the squeeze cast specimen is fast thus producing products of fine grains as compared to the slow cooling of sand casting, which produces large grains. Products obtained through squeeze casting are with improved mechanical properties.

## 2. Squeeze casting procedure

A metered quantity of molten metal was poured into the steel mould cavity at a supper-heat temperature of between 40-60 ^{0}C fast but avoiding turbulence. The upper die was then released to close the mould cavity with and without applying any load on the upper die. Thermocouples were inserted into the drilled holes made in the die, which were used to monitor both the die and cast metal temperatures. The terminals of the thermocouples were connected to the chart recorder/plotter (set at the highest speed of 10mm/s and voltage 100mV) through the cold junction apparatus, maintained at 0 ^{0}C throughout the measuring period.

## 3. Assumptions made

Heat transfer in the molten metal cast zone is due to both conduction and convection while conduction heat transfer occurs in the steel mould, it is convection at the outer surface of the steel mould.

The thickness of the cast specimen is much smaller than the diameter (radial dimension), thus giving one dimensional heat transfer process.

Considering the symmetrical nature of the cast specimen, solidification process was assumed symmetrical and only lower half of the specimen’s thickness was analysed see figure 1.

The bottom of the squeeze casting rig was and the heat losses to the atmosphere was small and neglected.

Heat losses through conduction and convection to the atmosphere at the punch were neglected, as a result of short time of pressure application.

The process of analyses in the cast specimen starts only when the steel mould cavity had been filled with the required quantity of liquid molten metal (i.e. heat transfer processes during pouring of molten aluminium into the steel mould are not considered).

Density of the molten and solidified aluminium metal was assumed to be the same and independent of temperature.

Thermal conductivity and specific heat of aluminium metal were dependent on the cast temperatures.

## 4. Heat transfer governing equations

### 4.1. Without pressure application on the cast metal

A measured quantity of molten aluminium metal was poured into the steel mould cavity. The process of solidification begins from the steel mould/cast metal interface and continues inwards into the cast metal. As this process continues, there was an increase in the thickness of the solidified layer and a decrease in the liquid molten metal portion. For the situation when nopressure was applied on the solidified molten metal, the governing heat transfer equations in one dimension are given by equation (1).

From figure 1, equation (1) is defined within the region with;

Steel mould,

Solidified molten,

Liquid molten metal,

At the phase change boundary condition;

where,

#### 4.3. Casting with pressure application on the solidified molten metal

As the cast aluminium metal solidifies, pressure is applied on the specimen, observing lapse or delay time, t_{l,} while varying the values of pressure applied. The time between the end of pouring of molten metal and pressure application known as lapse time, is recorded. This is necessary such that the cast specimen will not stick to the upper punch or cause the cast metal to tear with pressure application. Due to the applied pressure, an internal energy ∆q is generated within the solidified molten metal, (see figure 2).

Inserting the internal energy into the heat transfer equation (1) for solidified molten metal, it becomes equation (8),

where,

∆q-internal energy generated by applied pressure,

∆q_{P}-energy due to plastic strain within the solidified molten metal material,

∆q_{f}-frictional energy generated during pressure application,

∆q_{fP}-frictional energy due to punch / solidified molten metal interface,

∆q_{fm}-frictional energy due to steel mould cylindrical surface / solidified molten metal interface.

## 5. Formulation of heat transfer equations

Finite difference expressions for the nodal temperatures are obtained by either energy balance within an elemental volume around the node or by substitutions into the governing partial differential equations. Partial time derivative of temperature contained in the equations can be written in terms of a moving gradient as in White (1991) at a velocity of dx/dt. In the heat transfer equation, temperature is defined as a function of distance, r with time, t and represented in equation (11).

### 5.1. In the steel mould

where,

_{st}-thermal diffusivity of steel mould material.

### 5.2. In the solidified molten metal

where,

_{S}-thermal diffusivity of solidified molten metal material.

### 5.3. In the liquid molten metal

where,

_{L} -thermal diffusivity of liquid molten metal material.

## 6. Nodal divisions

The steel mould, the solidified metal and the molten metal regions were discretized separately. Each of these regions was divided into a fixed number of gridal points as in figure 1.

### 6.1. In the steel mould

I = 1, 2, 3,…, G-1

### 6.2. In the solidified molten metal

I = G+1, G+2, G+3, …, M-1

### 6.3. In the liquid molten metal portion

I = M+1, M+2, M+3 … N-1

### 6.4. In the phase change boundaries

The phase change is represented by the equations;

### 6.5. At the completion of solidification

As solidification time progresses, the boundary locations change, and the thickness of the solidified molten metal in the radial direction increases. The rate of change of boundary location with time is represented, mathematically in equation (25),

Where,

-time interval.

## 7. Boundary conditions

The problem of phase change during solidification is that the location of the solidifying molten metal / liquid molten metal interface is not known and this is determined continuously by appropriate mathematical analysis. This moving interface is normally expressed mathematically by the energy balance equations at the interfaces. In the numerical analysis, as solidification of molten aluminium metal progresses, three boundary interfaces occurred as:

### 7.1. Steel mould-atmosphere interface (I = 1)

The heat conducted to the steel mould material (from I = 2 to 1 = 1) equals the sum of the change in the internal energy and heat convected from the surface of the steel metal mould material into the atmosphere, mathematically represented in equation (26).

I = 1

#### 7.2. In the solidified molten metal-steel mould interface (I= G)

The sum of the heat conducted from the solidifying molten metal/steel mould interface and the change (decrease) in the internal energy at the boundary equal the sum of the heat conducted to the steel mould and the change (increase) in the internal energy at the interface.

I * =* G

### 7.3. In the liquid molten metal – Solidified molten metal interface, (I= M)

The sum of the heat conducted from the liquid molten metal and the internal energy generated equal the sum of heat conducted to the solidified molten metal and the internal energy generated at the interface as in equation (28).

I = M

### 7.4. At the corners

The usual one dimensional heat transfer analysis does not take into consideration heat loss at the corners as represented in figure 3. Heat loss at both the external and internal corner nodes of the steel mould have been considered and analysed while the die is lagged at the bottom surface.

### 7.5. External corner effect

The heat conducted from point R_{i+1}to point R_{i}added to the change in the internal energy equal the amount of heat convected out to the atmosphere at point R_{I}., Fig. (3a). The finite form of equation (29) is represented by equation (30),

### 7.6. Internal corner effect

The heat conducted from point P_{I+1} to point P_{I}added to the change in the internal energy amount to the heat conducted into the steel mould at point I, fig. (3b) thus becoming equation (31);

### 7.7. First time analysis

Solidification takes place only in the radial direction, a one-dimensional heat solidification problem was assumed numerically to take place in the radial direction only after filling the steel mould cavity with the liquid molten metal.

For the first time analysis, the specimen is considered to be in the molten stage and therefore equation (32) for the liquid molten metal is used for computation. Thus;

This equation is subjected to the boundary conditions with;

The instantaneous radius r_{i} in the first time analysis is given by equation (34);

I = G, G+1, G+2, G+3,…., N

The boundary velocity of moving coordinate of equation (32) is given as equation (35);

### 7.8. At the completion of solidification

At the completion of solidification, the whole molten region becomes solidified molten metal, and just before pressure is applied, the governing heat transfer equation becomes equation (36) representing the solidified portion;

The equation is applicable within the region defined by equation (37);

The boundary motion at the completion of solidification is as in equation (38);

I = G, G+1, G+2, G+3,…,.N

The boundary velocity at the completion of solidification is expressed as equation (39);

where,

L -constant value

## 8. Finite difference of governing heat transfer equations

The heat transfer equations generated in the cast metal and interfaces are written in the finite difference forms. These equations are presented in the various regions thus;

### 8.1. In the steel mould region

where,

I = 1,2,3,……………G

### 8.2. In the solidified molten metal portion

where,

I = G+1, G+2, G+3,………., M

### 8.3. In the liquid molten metal region

where,

I = M, M+1, M+2, M+3, …, N

### 8.4. In the phase change boundary condition (I = M)

where,

### 8.5. In the steel mould / atmosphere interface (I = 1)

I = 1

where;

I = 1, 2, 3, …, G

### 8.6. In the solidified molten metal / steel mould interface (I = G)

where,

### 8.7. In the liquid molten metal / solidified molten metal interface (I = M)

where,

I = M, M+1, M+2, M+3,…,N

### 8.8. External corner effect (I = 1)

where,

### 8.9. Internal corner effect (I = G)

where,

### 8.10. First time analysis

Finite difference form of equation (32) therefore becomes equation (58);

where,

I = G+1, G+2, G+3, ……….,.N

### 8.11. At the completion of solidification

Finite difference form of equation (36) becomes equation (60);

where;

I = G+1, G+2, G+3,…, (N-1)

## 9. Stability criteria

For stability criteria to be achieved, the values of temperature

## 10. Casting with pressure application and die heating

Pressure was applied only when the cast specimen was solidified, the governing heat transfer equation therefore, takes the form of solidified molten metal (completion of solidification). The finite difference of equation (8) is written as equation (62);

where,

T -temperature change resulting from pressure application

The cast specimen height, h_{c}, is pressure dependent and the relationship is expressed with the equation (63) below after performing series of experiment with various applied pressure;

where coefficient of correlation r = 0.996

The plastic flow stress, σ_{(T)}, is dependent on both the applied pressure, P, and die temperature, TM,(White, 1991), and expressed with the equation (64);

where coefficient of correlation r = 0.9508

## 11. Casting with die heating

Aluminium cast specimens were produced with die pre- heating temperatures of between 100- 300^{0}C without applying pressure on the solidifying aluminium metal. The die heating process was carried out, using three electric heater rods (100Watts each) that were connected to a.c supply. The required die temperatures were set and controlled, using a bimetallic thermostat.

## 12. Heat transfer coefficient evaluations

The method of calculating heat transfer coefficients as reported by Santos et al (2001) and Maleki et al, (2006) is based on the knowledge of known temperature histories at the interior points of the casting or mould together with the numerical models of heat flow during solidification. These temperatures are difficult to measure due to the difficulty in locating accurate position of thermocouple at the interface. Therefore, the inverse heat conduction problem based on non-linear estimation technique of Chattopadhyay,(2007) and Hu and Yu,(2002), has been adopted to determine the values of interface heat transfer coefficients, as a function of time during solidification of squeeze casting. Solidification of squeeze casting of aluminium involves phase change and therefore thermal properties of aluminium are temperature dependent, making the inverse heat conduction problem non-linear.

The governing heat transfer equation in one-dimensional cylindrical coordinates is given by equation (65):

Equation (65) holds within the boundary condition as expressed in equation (66):

The thermal conductivity K_{al}(Reed-Hill and Abbaschian 1973) and (Elliot, 1988) of aluminium is dependent upon casting temperature, T_{al}, and expressed in equation (67):

The heat flow across the casting/mould interface can be characterized by an average interfacial heat transfer coefficient, h_{al}(T) as obtained by Gafur et al (2003) and Santos et al (2004). This is expressed mathematically in equation (68):

The heat transfer coefficient, h, at the interface is estimated by minimizing the errors between numerically estimated and measured temperatures defined by equation (69):

where,

T_{est} and T_{exp} -are the estimated and experimentally measured temperatures at various thermocouples location and times,

n -iteration stage

## 13. Numerical simulations of differential equations

Squeeze casting consists of two stages, the first of which is mould filling: - the mould is filled with the required quantity of liquid molten metal; the second is cooling, this continues until the part has solidified completely. Controlling both stages is of major importance for obtaining sound casts with the required geometry and mechanical properties as observed by (Kobryn and Semiatin, (2000), Browne and O’Mahoney,(2001) and Martorano and Capocchi,(2000). When molten metal is poured into the mould cavity, it is initially in the liquid state with a high fluidity. It quickly becomes very viscous, in the early stage of solidification, and later completely solidifies (Gafur et al, 2003). For the numerical analysis of heat transfer problem, the appropriate set of equations were determined that described the heat transfer behaviour in the cast metal (Hearn, 1992). With the boundary conditions, initial conditions, and thermo-physical properties of the materials being known, it is possible to obtain the temperature and variation of the whole casting system (Ozisik, 1985) and (Liu et al (1993). Finite number at discrete points (Adams and Rogers,(1973), Shampire,(1994) and Bayazitoglu and Ozisik,(1988)) within the cast specimen was employed as the numerical method of solution. This method provides the temperature at a discrete number of points in the cast region. In the numerical method, the cast region is defined and divided into discrete number of points. As temperature difference is imposed in the system, heat flows from the high-temperature region to the low-temperature region as shown in figure 1.

To determine the temperature distribution, energy conservation equations were used for each of the nodal points of the unknown temperature at the interfaces and the cast regions((Incropera and Dewitt, 1985) and (Janna, 1988)). Temperatures were monitored at distance 2mm into the cast metal, represented by grid point M, and at the steel mould/cast metal interface (see figure 4). By using measured temperatures in both the casting and the steel mould, together with the numerical solutions of the solidification problem, heat transfer coefficients were determined based on Beck (1970)solution of the inverse heat conduction problem. The estimation of the surface heat transfer coefficients or heat flux density utilizing a measured temperature history inside a heat-conducting solid is called the inverse heat conduction problem (Cho and Hong 1996). This problem becomes non-linear, as the thermal properties (thermal conductivity, specific heat) are temperature dependent.

## 14. Experimental procedure

Chromel-Alumel thermocouples TC2, TC3, TC4, TC5 and TC6 were positioned on the sides of the cylindrical steel container, while TC1 and TC7 were positioned in the cast aluminium metal in the cylindrical and bottom flat surfaces respectively as shown in figure 4below.

Thermocouples of chromel-alumel type, 3mm in diameter were used to determine the solidifying temperatures of the cast molten metal and heating temperatures of the steel mould at the various positions in the cylindrical steel container of figure 4. The solidifying temperatures at both the cylindrical and flat bottom surfaces of the cast molten aluminium metal were monitored at a position 2mm (from the surface of the steel mould –cast aluminium metal interface) into the cast molten aluminium metal.

At the steel mould wall in the cylindrical surface, thermocouples were positioned at X2 = 4mm, X3 = 8mm, X4 = 12mm, X5 = 16mm and X6 = 20mm measured from the cast aluminium metal / steel mould interface to monitor the heating temperatures at these positions of the steel mould wall as shown in fig. 4. From the temperatures versus time curves obtained for each position in the steel mould, the interface heating temperature versus time curve at the cast aluminium metal / steel mould, for position when X = 0 was obtained by using the polynomial curve fitting method.

This was done by selecting a particular time of heating of steel mould, say t = 10sec. and drawing vertical lines cutting across the heating temperatures versus time curves at various thermocouples’ distances within the steel mould. At the point of intersection with each curve, the value of temperature was read against distance X, for the chosen time, t = 10sec. The value of interface steel mould / cast temperature at time say, t = 10sec.was determined at the steel mould / cast metal interface by substituting the value of X = 0 in the polynomial curve fitting equation(70) obtained from the values of temperatures at various distances in the steel mould at a chosen time, t = 10sec.

The temperature obtained by this method corresponds to the interface steel mould / cast metal temperature at a distance X = 0 for the chosen time t = 10sec. If this procedure is repeated for a number of time increments, the temperatures obtained with corresponding times represent the temperature at X=0, for such time increments. The graph of extrapolated temperatures versus time is drawn for position when X = 0 to represent the heating temperatures versus time curve at the steel mould / cast aluminium metal interface is shown in figure 5.

## 15. Interface heat transfer coefficients determination

Extrapolated temperature versus time curve of figure 5 for position when X = 0 (i.e. cast aluminium metal / steel mould interface) was used to determine the heat transfer coefficients of solidifying molten aluminium metal. It was used to determine the interface heat transfer coefficients in the cast aluminium metal / steel mould for no pressure and with pressure applications at both the cylindrical and bottom flat surfaces of the steel mould.

The interface heat transfer coefficients between the steel mould and cast aluminium metal at the cylindrical and bottom flat surfaces were determined from the extrapolated experimental heating temperature versus time curve obtained for position X = 0 and aluminium cast solidification temperature versus time curves obtained for the cylindrical and bottom flat surfaces, using equations (66) and (67).

The interface heat transfer coefficients were determined also numerically by the inverse method using the Finite Difference Method (FDM) and the obtained results were compared with the experimentally derived values.

## 16. Discussions of results

### 16.1. Temperature-time curves

Figure 6 shows typical temperature versus time curves for solidifying molten aluminium metaland steel mould respectively without the application of pressure on solidifying metal.

This figure shows the comparison of the numerical method usually applied by Cho and Hong (1996) to determine interface steel mould / cast metal temperature versus time curve with the extrapolated experimental method of this present work. The heating curve, as obtained through extrapolations of polynomial curves fitting equations and numerical methods are in close agreement and the deviations from the values obtained numerically varied from between 1.26- 19.31%.

Typical result obtained under pressure is also shown in figure 7, indicating the solidification and heating curves generated for solidifying molten aluminium metal and steel mould which follow the same patterns to the curves in figure 6.

With the application of pressure, the peak temperatures recorded are about the same 649^{o}C and 648^{0}C for a pressure of 85.86 MPa at the bottom flat and cylindrical surfaces of the steel mould respectively (see figure 7). The peak temperature (649^{o}C) obtained at the bottom flat surface of the steel mould under applied pressure is found to be higher than that temperature (607^{o}C) without pressure application. This effect may be associated to additional internal heat generated, resulting to higher temperature during pressure application on the solidifying molten aluminium.

### 16.2. Interface heat transfer coefficients with time

From the temperature with time curves of figure 6, the heat transfer coefficients for both cylindrical and bottom flat surfaces were determined for both numerical and calculated values and shown in figure 8. The maximum heat transfer coefficients of 2927.92 W/m^{2}K and 2975.14 W/m^{2}K are obtained at the cylindrical and bottom flat surfaces respectively for no pressure application, which is close to 2900 W/m^{2}K as obtained for pure aluminium by Kim and Lee (1997). The values of heat transfer coefficients decrease rapidly for both the cylindrical and bottom flat surfaces to a level of 866.70 W/m^{2}K and 969.50 W/m^{2}K respectively in 90 seconds. These values further decrease to 361.80 W/m^{2}K and 478.80 W/m^{2}K at these surfaces in another 150 seconds and further decrease then becomes not so noticeable.

From figure 8, the peak values of interface heat transfer coefficients are 2927.92 W/m^{2}K and 2956.73 W/m^{2}K as obtained by experimental and numerical determinations respectively at the cylindrical surface for no pressure application. For times within 40 seconds to 120 seconds the values of the interfacial heat transfer coefficients obtained numerically and experimentally are found to show higher values of about 19.83 % for numerical results to experimental results.

With the application of pressure on the solidifying aluminium metal, the heat transfer coefficients reach maximum values of 3085.34 W/m^{2}K and 3351.08 W/m^{2}K in the cylindrical and bottom flat surfaces respectively (see figure 9). These values also decrease to 847.80 W/m^{2}K and 783.63 W/m^{2}K in 240 seconds in the cylindrical and bottom flat surfaces respectively, while further decrease with time of solidification is no longer pronounced.

### 16.3. Interface heat transfer coefficients with solidification temperatures

Figures 10 and 12 show the calculated experimental interface heat transfer coefficients for solidifying molten aluminium metal as a function of solidification temperatures of the solidifying molten aluminium metal. Figure 10 shows the variation of heat transfer coefficient with solidification temperatures of aluminium at the cylindrical surface, while figure 12 is the interface heat transfer coefficients with solidifying temperature at the bottom flat surface with and without the application of pressure on the solidifying metal. From the two graphs, the maximum interface heat transfer coefficients obtained without pressure and with pressure application in the bottom flat surface of the steel mould are 2975.14 W/m^{2}K and 3351.08 W/m^{2}K respectively.

Figure 11 shows the numerical values of the variation of interface heat transfer coefficients with the solidification temperature of aluminium metal at the cylindrical surface with the application of pressure. The maximum value of heat transfer coefficient of 3397.29 W/m^{2}K at applied pressure of 85.86Mpa as compared to 3351.08 W/m^{2}K obtained through the experimental procedure.

At solidifying temperature above 600^{o}C, a sharp reduction in the interface heat transfer coefficients is noticed at both surfaces as is shown on figures 10, 11 and 12.

For temperatures below 500^{0}C, the interface heat transfer coefficients for both no pressure and pressure applications are close in values. This shows that at temperature below 500^{0}C, the effect of applied pressure is no longer significant on the interface heat transfer coefficient values. The drop in temperature results in solidification of the molten aluminium, which in turn leads to a drop in the heat transfer coefficient values. The effect of applied pressure on the heat transfer coefficients of aluminium becomes more pronounced at solidifying temperatures above 500^{0}C which was also reported by Cho and Hong (1996). Below this temperature, the effect of applied pressure on interface heat transfer coefficient values becomes less pronounced.

Therefore, from figures 10, 11 and 12, it is observed that the effect of applied pressure becomes more significant at temperature close to the liquidus temperature of aluminium as measured along the bottom flat surface of the steel mould (see figure 12). The maximum value of 3397.29 W/m^{2}K is obtained for pressure level of 85.86MPa as compared to 2975.14 W/m^{2}K for no pressure at the bottom flat surface of the steel mould.

### 16.4. Peak interface heat transfer coefficients with applied pressures

Figure 13 shows the variation of peak values of interface heat transfer coefficients with and without pressure applications. Higher experimental values of heat transfer coefficients are obtained at the bottom flat surface than at the cylindrical surface of the steel mould (see figure 13). This can be associated to greater effect of pressure application experienced at the bottom flat surface than at the side cylindrical surface, thus leading possibly to greater additional internal heat energy generated, and hence obtained higher values of heat transfer coefficients. The results of numerical determination of heat transfer coefficients as in figure 13 shows higher values as compared to the values obtained by experimental method. At applied pressure of 85.86MPa, the obtained heat transfer coefficients are 3397.29 W/m^{2}K and 3351.08 W/m^{2}K by numerical and experimental procedures respectively. From the curves of heat transfer coefficients obtained with temperatures of figures 10, 11 and 12, three distinct portions are noticed. These portions are easily differentiated by the aluminium solidification temperature. These temperatures are below 500^{0}C, solidus phase, 500 to 660^{0}C, liquidus-solidus phase, and above 660^{0}C liquidus phase of solidification of molten aluminium.

The liquidus-solidus phase, which occurs over a solidification temperature range of 500^{0}C to 660^{0}C, instead of a constant solidification temperature of 660^{0}C can be attributed to the presence of impurities such as silicon, magnesium and manganese in the commercially pure aluminium a fact supported by Higgins (1983).

The empirical equations for each of the distinct phase changes as a function of applied pressures and solidification temperatures are determined for both the experimental and numerical methods. The empirical equations obtained are for mean average values of heat transfer coefficients based on the experimental method at the cylindrical and flat bottom surfaces. These are equations (71-73):

Temperatures below 500^{0}C (solidus phase),

with coefficient of correlation r = 0.9545.

Temperatures between 500^{0}C and 660^{0}C (liquidus-solidus phase),

with coefficient of correlation r = 0.9884.

Temperatures above 660^{0}C (super heat, liquidus phase),

with coefficient of correlation r = 0.7825.

The empirical equations (74-76) obtained through the numerical methods are from the results of the computer simulations of heat transfer coefficients at the cylindrical cast metal / steel mould interface by the application of various applied pressures. These empirical equations are:

Temperatures below 500^{0}C (solidus phase),

with coefficient of correlation r = 0.969

Temperatures between 500^{0}C and 660^{0}C (liquidus-solidus phase),

with coefficient of correlation r = 0.964

Temperatures above 660^{0}C (super heat, liquidus phase),

with coefficient of correlation r = 0.772

### 16.5. Die heating effect

Figure 14 is the effect of die pre-heat temperatures on the values of heat transfer coefficients of aluminium metal without the application of pressure. From the figure, the heat transfer coefficients become lower with increase in die pre-heat temperatures. At die temperature of 95^{0}C, the heat transfer coefficient is 3185.34 W/m^{2}K and drop to a value of 2476.73 W/m^{2}K at die temperature of 300^{0}C. For all the die temperatures, there is a fall in the heat transfer coefficient’s values as solidification temperature decreases.

### 16.6. Comparison of heat transfer coefficient with semi-empirical method

The values of heat transfer coefficients determined using experimental, heat differential equations (numerical), and methods of semi-empirical equations are shown in figure 15 under a pressure application of 85.86MPa. From this graph, the peak values of interface heat transfer coefficient are 3358.19 W/m^{2}K and 3198.79 W/m^{2}K as obtained by heat differential and method of semi-empirical equations respectively for a pressure application of 85.86MPa. The heat transfer coefficients’ values for the three methods drop with time and are found to be 1708.03, 1976.81 and 1838.72 W/m^{2}K in 100seconds for experimental, differential and methods of semi-empirical equations respectively.

With die temperature of TM=150^{0}C, the peak heat transfer coefficients of 3088.99 W/m^{2}K, 3249.84 W/m^{2}K and 2982.60 W/m^{2}K are obtained for experimental, heat differential and method of semi-empirical equations as shown on figure 16 following the same pattern as in figure 15.

## 17. Conclusions

The following conclusions can be made from the present investigation:

The graph of temperature against time curves obtained by extrapolating to steel mould / cast metal interface by polynomial curve fitting to heating temperatures graphs with times at various steel mould locations are found to agree in values to the usual numerical methods obtained by previous authors. The interface heat transfer coefficients obtained by the numerical and experimental methods without the application of pressure are found to have values close to that of the numerical methods. The values of the numerical methods were higher by about 19.83%.

Effect of pressure application on the solidifying molten aluminium is more pronounced at casting temperatures above 500^{0}C of the cast aluminium specimen on the values of the interface heat transfer coefficients obtained. Interface heat transfer coefficients are found to decrease with increase in solidification time in both the cylindrical and bottom flat surfaces of the steel mould and thereafter remain fairly constant at temperature below 500^{0}C

Values of experimental peak heat transfer coefficients at the bottom flat surface are found to be higher with pressure application on the solidifying aluminium metal than at the cylindrical surface.

The empirical equations, relating the values of interface heat transfer coefficients with the applied pressures and solidification temperatures at three distinct stages of solidifying molten aluminium are determined and can be applied to determine the heat transfer coefficients.

The values of heat transfer coefficients obtained by heat differential equations incorporating internal heat energy and methods of semi-empirical equations are very close in values. The values as obtained by semi-empirical equations were higher by about 1.7% within the first 5 seconds of solidification.

The semi-empirical equations generated are flexible and could be used to predict the casting temperatures of other metals if the heat transfer coefficient values at the three phase changes are known.

### Nomenclature

C_{L- }Specific heat of liquid molten aluminium metal,C _{S}-Specific heat of solidified molten aluminium metal,C _{st}- Specific heat of steel mould,d _{sc}-Spatial step size in the solidified molten metal at completion of solidification,d _{L}-Spatial step sizes at the Liquid portion,d _{PS}-Spatial step size at the solidified molten metal phase change,d _{PL}-Spatial step size at the solidified molten metal phase change,d _{s}-Spatial step size in the solidified molten metal,h _{0}-Initial height of cast specimen metal (no pressure application),h _{PL}-Spatial step size at the liquid molten metal phase change,h _{PS}-Spatial step size at the solidified molten metal phase change,H ^{*}-Convective heat transfer of the surrounding atmosphere around the steel mould,I-grid point, G-grid point at the steel mould-solidified molten metal interface, M-grid point at the solidified molten metal-liquid molten metal interface, N-grid point at the centre of the solidifying cast specimen metal, K _{L}-Thermal conductivity of liquid molten metal,K _{st}-Thermal conductivity of steel mould material,K _{S}-Thermal conductivity of solidified molten metal,R _{o}-external radius of steel mould,L-internal radius of steel mould, r-radial axis/coordinate, T-Temperature, T _{∞}-ambient temperature,T _{i}-temperature at a specified grid point(i, j),T _{L}-temperature of liquid molten metal,T _{st,} T_{M}-temperature of steel mould,T _{MM}-melting temperature of cast (aluminium) metal,T _{p}-pouring temperature of liquid molten metal (super heat temperature),T _{S}-temperature of solidified molten metal,X ^{i}_{r}-thickness of solidified molten metal at a specified grid point-radial direction,Q _{}-Thickness of the cylindrical steel mould. |

### Greek symbols

_{L}-thermal diffusivity of the liquid molten metal,_{st}-thermal diffusivity of the steel mould material,_{S}-thermal diffusivity of solidified molten metal,-time interval step, _{st}-density of steel mould material,_{S}-density of solidified molten metal,_{L}-density of liquid molten metal,_{(T)}-Plastic flow stress (temperature dependent). |

## References

- 1.
Adams J. Alan Rogers David. F. 1973 Computer Aided Heat Transfer Analysis”, McGraw-Hill Publishing Company, Tokyo. - 2.
Bayazitoglu Yildiz. Ozisik M. Necati 1988 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium McGraw-Hill Book Company, New York. - 3.
Beck J. V. 1970 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Int. J. Heat Mass Transfer,13 703 716 - 4.
Bolton W. 1989 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Butterworths-Heinemann Limited, UK. - 5.
Browne David. J. O’Mahoney D. 2001 Interface Heat Transfer in Investment casting of Aluminium”, Metallurgical and Materials Transactions A, Dec.32A 3055 3063 - 6.
Callister William. D. Jr 1997 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium ^{th}Edition, John Wiley and Sons Inc. - 7.
Chattopadhyay Himadri. 2007 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium J. Materials Processing Technology, 186,174 178 - 8.
Cho I. S. Hong C. P. 1996 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Int. J. Cast Metals Res.,9 227 232 - 9.
Das A. Chatterjee S. 1981 Squeeze Casting of an Aluminium Alloy Containing Small Amount of Silicon Carbide Whiskers”, The Metallurgist and Materials Technologist,137 142 - 10.
Elliott R. 1988 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Butterworths, London, U.K. - 11.
Narayan Prabhu, (Gafur M. A. Nasrul Haque. K. 2003 Effects of chill Thickness and Superheat on Casting/Chill Interfacial Heat Transfer During Solidification of Commercially Pure Aluminium”, J Materials Processing Technology, 133,257 265 - 12.
(Hearn E. J. 1992 ), “Mechanics of Materials”, 2^{nd}Edition, Pergamon Press, UK. - 13.
Higgins Raymon. A. 1983 Engineering Metallurgy Part I:Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium ^{th}Edition, ELBS with Edward Arnold, UK. - 14.
Hosford William. F. Caddell Robert. M. 1993 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium ^{nd}edition, PTR Prentice-Hall Englewood, NJ. - 15.
Hu H. Yu A. 2002 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Modelling Simul. Mater Sci. Eng.,10 1 11 - 16.
Incropera Frank. P. Dewitt David. P. 1985 Fundamental of Heat and Mass Transfer”, 3^{rd}Edition, John Wiley and Sons, NY. - 17.
Janna William. S. 1988 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium SI Edition, Van Nostrand Reinhold (International0, U.K. - 18.
Kim T. G. Lee Z. H. 1997 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Int J, heat mass transfer, 40(15),3513 3525 - 19.
Kobryn P. A. Semiatin S. L. 2000 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Metallurgical and materials Transactions, August,32B 685 695 - 20.
Liu A. Voth T. E. Bergman T. L. 1993 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Int. J. Heat Mass Transfer,36 2 441 442 - 21.
Maleki A. Niroumand B. Shafyei A. 2006 Experimental Determination of Heat Transfer Coefficients During Squeeze Casting of Aluminium Materials Science and Engineering A, 428135 140 - 22.
Martorano M. A. Capocchi J. D. T. 2000 Heat Transfer Coefficient at the Metal-Mold Interface in the Unidirectional Solidification of Cu-8%Sn alloys”, Intl J. Heat Mass Transfer, 43,2541 2552 - 23.
Ozisik M. Necati 1985 Heat Transfer:- A Basic Approach”, McGraw-Hill Publishing, Company U.K. - 24.
Potter D. A. Easterling K. E. 1993 Phase Transformations in Metals”, 2^{nd}Edition, Chapman & Hall, London. - 25.
Reed-Hill R. E. Abbaschian R. 1973 Physical Metallurgy Principles”, 3^{rd}Edition, PWS-KENT Publishing Company, Boston. - 26.
Santos C. A. Quaresma J. M. V. Garcia A. 2001 Determination of Transient Heat Transfer Coefficients in Chill Mold Castings”, Journal of Alloys and Compounds, 139,174 186 - 27.
Santos C. A. Garcia A. Frick C. R. Spim J. A. 2004 Evaluation of heat transfer coefficients along the secondary cooling zones in the continuous casting of steel billets, Inverse problems, Design and Optimization symposium”, Rio de Janeiro,1 8 - 28.
Shampire Lawrence. F. 1994 Numerical Solution of Ordinary Differential Equations”, Chapman & Hall, New York. - 29.
White Frank. M. . 1991 Heat and Mass Transfer”, Addison Wesley Publishing Company, Massachusetts.