In the field of metallurgical engineering, the production of high-quality sand casting parts is critical for applications such as blast furnace cooling components. As a researcher focused on casting processes and finite element simulation, I have investigated the solidification behavior of copper-based sand casting parts, specifically water-cooled crucibles, using a combined metal mold and sand mold approach. This study aims to provide insights into optimizing casting processes for sand casting parts, ensuring superior mechanical properties and minimizing defects through controlled solidification. The importance of sand casting parts lies in their widespread use in industrial settings, where factors like temperature gradients and sequential solidification directly impact performance and longevity. Here, I present a comprehensive analysis based on finite element simulations, incorporating theoretical models and empirical data to elucidate the solidification dynamics of these sand casting parts.
The foundation of this work rests on the application of ANSYS finite element software to simulate the transient thermal behavior during solidification. Sand casting parts often involve complex geometries and material interactions, making numerical simulation an invaluable tool for predicting outcomes. My approach involves modeling a pure copper casting with dimensions of 2 m × 1 m × 0.12 m, produced using a sand mold on the top and a tapered ductile iron metal mold on the bottom, with elliptical cooling channels formed by sand cores. This configuration is representative of typical sand casting parts used in high-temperature environments. To streamline the simulation, several assumptions were made: the liquid metal fills the mold cavity instantaneously, no macroscopic convection occurs within the melt, the material is homogeneous and isotropic in both solid and liquid phases, and the gating system is omitted to reduce computational time. These simplifications are common in studies of sand casting parts, allowing for efficient analysis while capturing essential physics.
Initial conditions were set to reflect practical casting scenarios. The molds were preheated to 200–300°C using hot air, establishing a steady-state temperature distribution prior to pouring. The initial temperature of the casting was set at the pouring temperature of 1150°C, which is typical for copper sand casting parts. Boundary conditions accounted for various heat transfer mechanisms, including conduction at the casting-mold interface, convective heat transfer with the surrounding air, and radiation losses. Specifically, the interfaces between the casting and both the sand and metal molds were modeled with appropriate heat transfer coefficients, while the external surfaces of the molds and the cooling channel interiors were subjected to natural convection. These details are crucial for accurately simulating the thermal history of sand casting parts, as they influence cooling rates and solidification patterns.
The thermal properties of the materials play a pivotal role in the simulation. For pure copper, the thermal conductivity varies with temperature, as summarized in Table 1. The density and specific heat capacity are temperature-dependent, described by the following equations:
$$ \rho = 8900 – 0.2667 (T – 25) \, \text{kg} \cdot \text{m}^{-3} $$
$$ c_p = 385 + 0.0998 T \, \text{J} \cdot \text{kg}^{-1} \cdot \text{K}^{-1} $$
where \( T \) is the temperature in °C. For the mold materials, ductile iron has temperature-dependent thermal conductivity, while the sand mold properties are considered constant. Table 2 provides these parameters. The release of latent heat during solidification is handled by defining enthalpy as a function of temperature:
$$ \Delta H(T) = \int_0^T \rho \cdot c_p(t) \, dt $$
This ensures that the phase change effects are accurately captured in the simulation of sand casting parts.
| Temperature (°C) | Thermal Conductivity (W·m⁻¹·°C⁻¹) |
|---|---|
| 20 | 399 |
| 100 | 387 |
| 200 | 379 |
| 400 | 374 |
| 600 | 363 |
| 800 | 353 |
| 1083 | 321 |
| Material | Temperature (°C) | Thermal Conductivity (W·m⁻¹·°C⁻¹) | Density (kg·m⁻³) | Specific Heat Capacity (J·kg⁻¹·°C⁻¹) |
|---|---|---|---|---|
| Ductile Iron | 20 | 42.3 | 7100 | 500 |
| 200 | 36.5 | |||
| 400 | 30.0 | |||
| 800 | 21.2 | |||
| 1000 | 17.0 | |||
| Sand Mold | 20 | 0.58 | 1700 | 1220 |
The simulation was performed using ANSYS with a transient analysis type. A steady-state analysis was first conducted with a time step of 0.01 s to establish the initial temperature field, which was then used as the starting point for the transient analysis. The total simulation time was set to 1200 s to cover the entire solidification and cooling process. This methodology allows for a detailed examination of temperature evolution in sand casting parts, providing data on solidification fronts and thermal gradients.
Results from the simulation reveal that the solidification of these sand casting parts follows a progressive pattern. The temperature distribution on longitudinal sections shows that solidification initiates at the end away from the gating system and progresses toward the feeder, while simultaneously moving from the metal mold side to the sand mold side. This sequential solidification is advantageous for sand casting parts, as it promotes effective feeding and reduces shrinkage defects. The sand mold side exhibits delayed solidification due to its lower thermal conductivity, allowing liquid metal to compensate for solidification shrinkage on the metal mold side. Additionally, tilting the mold by 30° after pouring enhances this effect, ensuring that the gating system acts as a riser for optimal feeding. Such insights are vital for designing robust processes for sand casting parts.

To quantify the solidification behavior, I analyzed solidification curves at various cross-sections of the sand casting parts. Three transverse sections were selected: at x = 0.35 m, x = 1.1 m, and x = 1.83 m, corresponding to locations along the casting length. At each section, temperatures were recorded at different thickness positions: the metal mold-casting interface (z = 0.15 m), near the cooling channels (z = 0.21 m), and the sand mold-casting interface (z = 0.27 m). The curves display distinct phases, with the sand mold interface showing a clear solidification plateau at 1083°C, while the metal mold interface does not due to rapid cooling. This underscores the heterogeneous cooling inherent in sand casting parts when using composite molds. Furthermore, on a longitudinal section at y = 0.826 m and z = 0.213 m, temperatures at different lengths (x = 0.172 m, 1.08 m, 1.92 m) show that regions closer to the gating system remain hotter, increasing temperature gradients over time and reinforcing sequential solidification in sand casting parts.
For a theoretical interpretation, I applied the Schwarz model to describe these solidification curves. The model is expressed as:
$$ T_n(x, t) = A_n + B_n \, \text{erf} \left( \frac{x}{2 \sqrt{a t}} \right) $$
where \( T_n \) is the temperature at position \( x \) and time \( t \), \( a \) is the thermal diffusivity, and \( A_n \) and \( B_n \) are constants determined by boundary conditions. This model is particularly useful for analyzing sand casting parts, as it separates the solidification into two stages: liquid undercooling and solid cooling. For the liquid undercooling stage, at \( x = 0 \), the temperature is constant and equal to the pouring temperature:
$$ T_1(0, t) = A_1 = \text{constant} $$
For the solid cooling stage, at \( x = 0 \), the temperature represents the average interface temperature:
$$ T_2(0, t) = A_2 = \text{constant} $$
Fitting this model to the simulated data yielded excellent agreement, confirming its applicability to sand casting parts. The parameter \( A_1 \) ranged from 1112 to 1146°C, closely matching the pouring temperature of 1150°C, while \( A_2 \) varied with location, reflecting lower interface temperatures farther from the gating system and on the metal mold side. This analysis demonstrates that the Schwarz model effectively captures the thermal dynamics of sand casting parts during solidification.
Temperature gradients are a key metric for assessing solidification quality in sand casting parts. I calculated gradients on both longitudinal and transverse sections. On the longitudinal section at y = 0.826 m and z = 0.213 m, the temperature gradient increases with time, reaching a maximum of approximately 73°C/m during cooling. This value surpasses the minimum required for progressive solidification in copper castings, estimated at 10.1°C/m based on literature, indicating that sequential solidification is achievable in these sand casting parts. On transverse sections, the gradients exhibit more complex behavior due to latent heat release. For instance, at sections x = 0.35 m, 1.1 m, and 1.83 m, maximum gradients of 272°C/m, 352°C/m, and 369°C/m were observed around 74 s, respectively. These peaks occur as initial solidification releases latent heat, creating negative gradients that later turn positive. As solidification progresses, latent heat accumulation reduces gradients to minima of 142°C/m, 208°C/m, and 238°C/m, before gradients rise again post-solidification. Notably, gradients are smaller near the gating system, consistent with sequential solidification in sand casting parts. Compared to previous studies using steel metal molds, the switch to ductile iron molds here enhances temperature gradients, further benefiting the integrity of sand casting parts.
The implications of these findings are significant for the manufacturing of sand casting parts. The progressive solidification mode ensures that the hot face (working surface) of the casting solidifies first, with the sand mold side providing adequate feeding to prevent defects. This is crucial for sand casting parts used in demanding applications like blast furnace cooling, where reliability is paramount. The use of composite molds with different thermal properties allows for tailored cooling rates, optimizing the microstructure and mechanical properties of sand casting parts. Moreover, the Schwarz model provides a robust framework for predicting solidification curves, aiding in process design for various sand casting parts. By controlling parameters such as mold materials and pouring temperatures, manufacturers can achieve desired temperature gradients and solidification sequences, reducing trial-and-error in production.
In terms of practical applications, this research highlights the importance of finite element simulation in advancing sand casting technology. For instance, simulating different mold configurations or cooling channel layouts can help optimize the design of sand casting parts before physical prototyping. Additionally, the analysis of temperature gradients can guide the placement of chills or insulation to control solidification in critical regions of sand casting parts. These approaches contribute to higher yield rates and improved performance of sand casting parts in industries ranging from metallurgy to automotive.
To further elaborate, let’s consider the mathematical modeling aspect. The heat transfer during solidification of sand casting parts can be described by the Fourier equation with a phase change term:
$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + L \frac{\partial f_s}{\partial t} $$
where \( k \) is thermal conductivity, \( L \) is latent heat, and \( f_s \) is the solid fraction. This equation underpins the finite element simulation, accounting for the complex interactions in sand casting parts. By discretizing the domain and solving numerically, we obtain detailed temperature fields that inform solidification analysis. For sand casting parts with irregular geometries, adaptive meshing techniques can be employed to improve accuracy, though this was not necessary in the current study due to the relatively simple shape.
Another aspect is the effect of cooling channels on solidification. In these sand casting parts, the elliptical channels act as internal heat sinks, influencing local cooling rates. The simulation shows that regions near the channels solidify faster, contributing to the overall temperature gradient. This is beneficial for achieving directional solidification in sand casting parts, as it helps align the solidification front with the desired direction. However, care must be taken to avoid excessive thermal stresses that could lead to cracking in sand casting parts, a topic for future investigation.
In comparison to other casting methods, such as die casting or investment casting, sand casting parts offer advantages in terms of cost and flexibility for large components. However, controlling solidification is more challenging due to the lower thermal conductivity of sand molds. This study demonstrates that combining sand with metal molds can mitigate this issue, making it a viable strategy for high-performance sand casting parts. The key is to balance the cooling rates to avoid defects like porosity or hot tears, which are common concerns in sand casting parts.
From a materials science perspective, the solidification kinetics of copper in sand casting parts can be analyzed using growth models. For example, the interface velocity during solidification can be estimated from temperature gradients and cooling rates. Assuming a planar solidification front, the growth velocity \( v \) is related to the temperature gradient \( G \) and the cooling rate \( \dot{T} \) by:
$$ v = \frac{\dot{T}}{G} $$
For the sand casting parts in this study, the cooling rates near the metal mold interface are higher, leading to faster growth and finer microstructures. This can enhance the mechanical properties of sand casting parts, such as strength and ductility. Microstructural analysis of actual castings would complement these simulations, but the numerical results provide a strong foundation for understanding the behavior of sand casting parts.
In conclusion, this work underscores the value of finite element simulation in optimizing the solidification of sand casting parts. The sequential solidification achieved with metal-sand composite molds ensures high-quality castings with minimal defects. The Schwarz model offers a reliable way to interpret solidification curves, while temperature gradient calculations confirm the feasibility of progressive solidification in sand casting parts. Future research could explore other materials or more complex geometries for sand casting parts, as well as incorporate fluid flow simulations to account for filling effects. By continuing to refine these models, we can advance the production of sand casting parts for critical engineering applications, ensuring durability and efficiency in their end-use environments.
To summarize key points in a tabular format, Table 3 lists the maximum temperature gradients observed in different sections of the sand casting parts, highlighting the influence of mold design. This data can guide engineers in selecting appropriate parameters for similar sand casting parts.
| Section | Position | Maximum Temperature Gradient (°C/m) |
|---|---|---|
| Longitudinal | y = 0.826 m, z = 0.213 m | 73 |
| Transverse | x = 0.35 m | 272 |
| Transverse | x = 1.1 m | 352 |
| Transverse | x = 1.83 m | 369 |
Overall, the integration of simulation and theoretical analysis provides a comprehensive framework for improving the manufacturing of sand casting parts. As demand for reliable and efficient components grows, such studies will play an increasingly important role in the foundry industry, driving innovations in sand casting technology.
