alexa Casting Simulation and Prediction of Shrinkage Cavities

ISSN: 2168-9679

Journal of Applied & Computational Mathematics

  • Research Article   
  • J Appl Computat Math 6:371, Vol 6(4)
  • DOI: 10.4172/2168-9679.1000371

Casting Simulation and Prediction of Shrinkage Cavities

Jabur AS1* and Kushnaw FM2
1University of Basrah, Iraq
2University of Technology, Baghdad, Iraq
*Corresponding Author: Jabur AS, University of Basrah, Iraq, Tel: 9647802443867, Email: [email protected]

Received Date: Sep 08, 2017 / Accepted Date: Oct 13, 2017 / Published Date: Oct 23, 2017


This work aims to predict the shrinkage defects in Al-Si castings by determination the suitable parameters and techniques which can be applied in the casting simulation system. Also, it aims to specify the role of silicon content in amount, morphology, and distribution of these defects. The Numerical solution has been carried out using an explicit 3-D finite difference method for the given system of the casting and a mold. Additionally, an experimental casting of the studied samples was achieved. It was found that the shrinkage porosities increased with increasing the silicon content up to 7%, so at this peak, they spread in all cast regions and cannot be predicted. The low silicon alloys suffered from only the cavities defects that can be predicted by mapping the solidus time contours. Finally, it was concluded that the critical temperature gradient value of the porosities development in the eutectic (Al-12%Si) alloys was 1.3°C/cm.

Keywords: Simulation; Al-Si alloys; Casting structures; Numerical analysis; Shrinkage defects


Producing sound aluminum-silicon castings is one of great economic significance to the foundry industry; however, shrinkage defects are the main categories of casting defects which take many forms; shrinkage cavity, porosity and surface sink. Such defects have a negative economic impact on casting production; their consequences range from high rework costs to casting rejection as reported by investigators [1]. During recent years, the application of some popular commercial software’s as computer simulation tools has become widely accepted within the foundry industry [2]. The application of casting simulation has been most beneficial, for avoiding shrinkage scrap, by predication of shrinkage defects without having to discover them in the foundry through the usual trial and error process [3], which can be very tedious, time consuming, and expensive [4]. One of the methods to describe solidification modeling is the heat transfer model, which solves the energy equation [5]. There have been few studies to model the effect of solidification shrinkage from different perspectives. Trovant and Argyropoulos [6] proposed an algorithm to account for shrinkage and consequently determine the shrinkage profile resulting from phase and density change. Only the energy equation is solved in the domain and the effects of solidification shrinkage are imposed through the proposed algorithms. Bride [7] presented a model that accounts for shrinkage during the directional solidification or dendritic binary alloy under the assumption that the densities of liquid and solid phase are different but constant. Also, Kim and Ro [8] performed a study on shrinkage formation during solidification of a material in a twodimensional rectangular cavity. Anderaus and Delliisola [9] proposed a 2-D finite difference model to describe the influence of solidification shrinkage on the final shape of cast ingot. Paszndiden-Fard [10] develop a numerical method to simulate the solidification shrinkage of the solidifying liquid [11-13].

The aim of this work is to find the suitable criteria or methods that can be used in prediction the shrinkage defects. Also, it aims to study the effect of silicon content on the morphology, quantity, and distribution of the shrinkage defects in the Al-Si castings.

Mathematical Model

Metal region

The variation of temperature T with time could be described by solving the Fourier heat equation, which is generally written as:

Equation (1)

Where convection is neglected. The density denotes as ρ, specific heat of casting as Cp and the term Equation is the transient term where, τ is a time step. In the right side the three terms represent the heat conduction in three-direction with thermal conductivity K. A source term Equation is an internal heat generation for solidification under equilibrium condition. This term can be described as:

Equation (2)

Where L denotes latent heat distributed over the solidification range, so it is a function of the solid fraction (fs) as confirmed by Poirier and Salcudean [11].

Thermo physical properties in the mushy region will depend upon the amount of liquid and solid (i.e. solid fraction) and can be calculated as following:

Equation (3)

Where fs(T) are the solid fraction at T and Psol and Pliq are values of the property at solidus and liquidus temperatures, respectively.

If the alloy has very narrow solidification interval or it is a pure metal or a eutectic, TS will be close to or equal to TL. In this case an artificial temperature interval has to be assumed. This temperature interval should be as small as possible to be near the real material behavior.

Equation (4)

Mold region

For heat conduction in the mold, the energy equation is:

Equation (5)

the subscript m is denoted the mold and Km is the mold thermal conductivity.

Boundary conditions

On the outer surface, heat transfer is done by convection as following:

Equation (6)

where Tm is the temperature of outer surface of mold and Ta is ambient temperature.

Assumptions made

The following assumptions were made for the formulation of the problem:

1. Mold is filled instantaneously.

2. Phase transformation proceeds in one direction: liquid to solid.

3. No convection is present in the liquid metal.

4. Thermal properties of the metal casting for liquid and solid phase are different but constant.

The density considered variable with temperature.

Numerical solution

After discretization of the whole domain and using the forward difference for the time derivative and the central difference for the spatial second derivative, the approximation resulting in :

Equation (7)

Equation (14) is called the Forward Time, Centered Space or FTCS approximation to the heat equation. It can be reduced as:

Equation (8)

The FTCS scheme is easy to implement because the values can be updated independently of each other. The entire solution is contained in two loops: an outer loop over all time steps and an inner loop over all interior nodes. Notice that The FTCS can yield unstable solutions that oscillate and grow if Δτ is too large. According to Jaluria [12], stable solutions with the FTCS scheme are only obtained if;

Equation (9)

Computer programming

A computer program has been built in using FORTRAN 90, to carry out the numerical solution, on a personal computer. The basic structure of this program is shown in Figure 1. Dimensions of casting are illustrated in Figures 2 and 3.


Figure 1: Basic Structure of FD- program.


Figure 2: Straight bar samples.


Figure 3: Rectangular Samples.

Experimental Procedure

After melting the aluminum alloys in a gas furnace and degassing with a C2Cl6 tablets, the molten metal were poured in sand molds. Seven straight bar samples (Figure 2) with different dimensions were prepared to study the effect of silicon percentage on the shrinkage defects. Also, a rectangular sample (Figure 3) with fin was prepared to evaluate the efficiency of the program in prediction the cavities location. Details of these samples were contained in Table 1.

Sample code Composition Cast length mm Cast width mm Cast thickness mm Feeder length mm Feeder Diameter mm Pouring temp. °C Solidus temp. °C Liquius temp. °C
1 Al-12% Si 500 50 30 145 45 610 542 573
2 Al-12% Si 5 520 35 30 145 45 610 542 573
3 Al-12% Si 480 70 30 145 45 610 542 573
4 Al-7% Si 4 470 45 30 145 45 650 567 614
5 Al-0.85% Si 480 45 30 145 45 690 643 648
6 Pure Al 500 50 30 145 45 705 660 660
7 Pure Al 100 60 30 - 45 705 660 660

Table 1: Dimension of the casting samples.

Results and Discussion

Verification of the model

In order to check the program, numerical simulation results were compared with experimental data. It conducted by inserting thermocouple in the cast Sample (1). Figure 4 shows the simulated and experimental cooling curves at cast center position. It is obvious that the slope of the simulated curve after solidification is smaller than experimental one which may be due to the effect of solid contraction and consequently the air gap development. But, at the interval between the liquidus and solidus temperature (local solidification time), It represents a quite similar between them.


Figure 4: Comparison Between Simulated And Experimental Results At The Centre Of The Cast (1).

The effect of Si percentage on the shrinkage morphology

Four samples with different silicon percentages are arranged in Table 2. It contains the shrinkage morphologies extracted from Figure 5 shows that, the pure aluminum (sample 6), suffered from only one type of shrinkage defects. Feeder tends to pipe deeply; the pipe is smooth and geometric in shape, with the absence of any type of shrinkage porosities as shown in Figures 6 and 7. This can be attributed to the plane front solidification behavior of pure alloys. Two reasons to the severe extent of pipe in these samples. First, the high solidification contraction of pure aluminum (7.14%). Second, the absence of the other types of shrinkages leads to accumulation in the pipe only.

Sample Code % Si Pipe Centerline porosities Dispersed porosities Surface sink
1 12 Medium Present Nil Nil
4 7 Light Nil Severe Nil
5 0.85 severe Nil Nil Present
6 0 severe Nil Nil Nil

Table 2: Shrinkage Defect Morphologies Of The Casting Samples.


Figure 5: Photographs of A-Sample (1), B-Sample (4), C-Sample (5), And D-Sample (6) Sections, Which Illustrate The Shrinkage Defect Morphologies.


Figure 6: A Photograph of Polished Section of Sample (4) Shows The Dispersed Porosities.


Figure 7: A Photograph of Polished Section of Sample (6) Shows The Porosity Free Surface.

Sample (5) of low silicon content (0.85%) had approximately the same shrinkage behavior as the pure aluminum samples. Here the pipe is less regular and the end is ragged as shown in Figure 5c.

Sample (4) of 7% silicon exhibited light pipe and a dispersed porosity in all cast regions, as shown in Figure 5b and Figure 6. This is explained by the long freezing range of this composition which offered a high constitutional super cooling and an effective dendritic solidification behavior and consequently a limited feeding to the interdendritic shrinkage regions so that a nucleation of porosities enhanced. On the other hand, the reason of the presence of small pipe in comparison with the previous alloys is firstly; that the total contraction percentage of Al-Si alloys decreases with increasing the silicon content. Secondly; the total contraction here shared between the pipe and wide dispersed porosities.

Sample (1) of 12% silicon (Figure 5a and Figure 8) exhibited medium size pipe and localized porosities in the centerline. The short freezing range of this alloy produced a simple constitutional super cooling resulted in cellular or light dendritic freezing mode for the eutectic layers or for the aluminum grains which contain the porosities in between. The pipe of this sample is smaller than that of sample (6) because the total contraction percentage of this alloy is 3.8%; also it is larger than that of sample (4) because the porosities here are localized in the centerline only.


Figure 8: A Photograph of Polished Section of Sample (1) Shows The Centerline Porosities.

Simulation results

Aluminum - 12% Silicon

Figures 9-11 shows the distribution of the temperature gradient through (x-y) mid surface of samples (1, 2 and 3) respectively. The temperature gradients near the wall area are high due to rapid wall chilling.


Figure 9: The Gradient variation over the mid surface of sample (1).


Figure 10: The Gradient variation over the mid surface of sample (2).


Figure 11: The Gradient variation over the mid surface of sample (3).

Figure 12 contains the values of the porosity length extensions in these samples in addition to the gradient values at the starting and ending points. The region of shrinkage coincides with the region of shallow temperature gradient (center line shrinkage). In fact, all of the shrinkage porosities were found within the region of calculated temperature gradient of lower than 1.3°C/cm. Therefore, 1.3°C/cm was tentatively taken as the critical value of shrinkage occurrence. The equivalence of these values at the two points for a one sample and between the three different dimensions samples gives a good support for this proposal. Niyama et al. [13] reports the values of the critical temperature gradient; 2 to 3°C/cm for the steel bars. The difference between Niyama values and ours is not surprising because the values must depend on casting type, alloy composition. Niyama did not take into account the change in thermal properties for liquid and solid. While Niyama’s criterion works well for low carbon steel that has very small difference between solid and liquid properties that makes its application by many nonferrous foundries is questionable because they have large diverge of thermal properties between solid and liquid metal.


Figure 12: Values of Critical Temperature Gradients Which Is Gained From Simulation.

Aluminum - 7% silicon

Figure 5b and Figure 6 include photographs of the polished section of sample (4), which show a wide dispersed of shrinkage porosities through the sample section without any concentration in a specified region. Figure 13 illustrates the distribution of the temperature gradient through the mid surface (x-y). It clarifies that; the shrinkage porosities distribution in this silicon concentration does not depend on this parameter values, and they cannot be predicted with this tool.


Figure 13: The Gradient variation over the mid surface of sample (4).

Aluminum - 0.85% silicon

Figure 14 shows the distribution of the temperature gradient through the mid surface (x-y) for sample (5). Also, the gradient decreased for very low values and did not cause a development of porosities. Therefore, their solidification behavior resemble to the pure aluminum.


Figure 14: The Gradient variation over the mid surface of sample (5).

The pure aluminum

The first extension of the modeling is the evaluation of the temperature gradient over the mid surface (x-y) of the pure aluminum (sample 6) as in Figure 15. In general, the gradients are highest at the casting edges, but decrease with moving toward the casting center. These samples are free from the porosities in a different manner than steel, i.e. the development of the porosities depends on the decreasing of the temperature gradient under a critical value as the work of Niyama et al. [13]. Consequently, the development of the porosities in pure aluminum did not depend on the gradient values because the solidification front type (planar) has the predominant role in this manner and the gradient has no role.


Figure 15: The Gradient variation over the mid surface of sample (6).

The solidification time can be used in the cavities prediction .The appearance of a closed loop in the solidification time contours is the suitable criterion for this type of shrinkage. For checking the efficiency of the present simulation program with employing the solidification time criterion in prediction the cavities position, a rectangular sample designed to produce an internal cavity was tested (sample 7). Figure 16 includes a comparison between the sliced sample (to identify the position of the real cavity) and a simulated graph of the solidus time contours over the same sliced surface. There is a good agreement between them that make this criterion a good tool in cavities prediction and a successful of the present computer program in the simulation function to predict the hot spot and consequently the cavity.


Figure 16: The solidus time contours over the mid surface of sample (7).


From the present work, the main conclusions can be summarized as:

1. The shrinkage porosities increase with increasing the Si content up to 7%, so at this peak, they spread in all cast regions and cannot be predicted.

2. The low Si and pure Al-alloys castings suffered from the shrinkage cavities defects only which can be predicted by mapping the solidus time contours.

3. The critical temperature gradient value of the shrinkage porosities development in the eutectic (Al-12%Si) alloys was 1.3°C/cm.


Citation: Jabur AS, Kushnaw FM (2017) Casting Simulation and Prediction of Shrinkage Cavities. J Appl Computat Math 6: 371. Doi: 10.4172/2168-9679.1000371

Copyright: © 2017 Jabur AS, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Select your language of interest to view the total content in your interested language

Post Your Comment Citation
Share This Article
Article Usage
  • Total views: 1254
  • [From(publication date): 0-2017 - Jun 19, 2018]
  • Breakdown by view type
  • HTML page views: 1186
  • PDF downloads: 68

Post your comment

captcha   Reload  Can't read the image? click here to refresh
Leave Your Message 24x7