In the sand casting foundry industry, the production of complex valve bodies often encounters defects such as shrinkage porosity and cavities due to improper gating system design and process parameters. To address these challenges, I employed the ProCAST software to perform a comprehensive numerical simulation of the filling and solidification process for a typical valve body. The goal was to identify the root causes of defects and to optimize the casting process, thereby improving product quality and yield in the sand casting foundry environment.
The valve body studied features a total length of 300 mm, stepped inner diameters of 130 mm, 160 mm, and 180 mm, and an average wall thickness of 10 mm. It includes four boss holes on the base, a front ear seat of 15 mm thickness, a top channel with an inner diameter of 80 mm and wall thickness of 8 mm, and a flange of 20 mm thickness. Such geometric complexity makes the sand casting foundry process prone to thermal and fluid flow issues. I created a 3D model of the valve body using Pro/E and then imported it into the ProCAST environment for mesh generation and simulation setup.
Theoretical Framework
The flow of molten metal during sand casting foundry is governed by the conservation laws of mass, momentum, and energy. For an incompressible, viscous, and unsteady flow with a free surface, the governing equations are as follows:
Continuity equation (mass conservation):
$$
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0
$$
Momentum conservation equations in the x, y, and z directions:
$$
\rho \left( \frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z} \right) = -\frac{\partial p}{\partial x} + \rho g_x + \mu \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)
$$
$$
\rho \left( \frac{\partial v}{\partial t} + u\frac{\partial v}{\partial x} + v\frac{\partial v}{\partial y} + w\frac{\partial v}{\partial z} \right) = -\frac{\partial p}{\partial y} + \rho g_y + \mu \left( \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} + \frac{\partial^2 v}{\partial z^2} \right)
$$
$$
\rho \left( \frac{\partial w}{\partial t} + u\frac{\partial w}{\partial x} + v\frac{\partial w}{\partial y} + w\frac{\partial w}{\partial z} \right) = -\frac{\partial p}{\partial z} + \rho g_z + \mu \left( \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} + \frac{\partial^2 w}{\partial z^2} \right)
$$
Energy conservation equation:
$$
\rho c \left( \frac{\partial T}{\partial t} + u\frac{\partial T}{\partial x} + v\frac{\partial T}{\partial y} + w\frac{\partial T}{\partial z} \right) = \lambda \left( \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} \right) + Q
$$
In these equations, u, v, w are velocity components, ρ is density, p is pressure, μ is dynamic viscosity, g is gravitational acceleration, c is specific heat, T is temperature, λ is thermal conductivity, and Q represents internal heat sources. These equations are solved numerically by ProCAST to predict the flow and solidification behavior in the sand casting foundry process.
| Symbol | Description | Unit |
|---|---|---|
| u, v, w | Velocity components | m/s |
| ρ | Fluid density | kg/m³ |
| p | Pressure | Pa |
| μ | Dynamic viscosity | Pa·s |
| g | Gravitational acceleration | m/s² |
| c | Specific heat capacity | J/(kg·K) |
| λ | Thermal conductivity | W/(m·K) |
| T | Temperature | °C or K |
Initial Gating System Design and Simulation Setup
The initial gating system for the sand casting foundry was designed using a two-part molding approach. A middle injection scheme was adopted where the gating system was placed at the parting line on the right side of the mold. Two ingates connected via a runner to a sprue and pouring cup. Risers were added at three locations: the front ear seat, the top flange, and the rear outer diameter where heavy sections existed. The 3D model of the initial design was saved in Parasolid format and imported into the MeshCAST module of ProCAST for meshing. A total of 25,673 elements were generated to balance accuracy and computational efficiency.
| Parameter | Value |
|---|---|
| Material | Cast steel |
| Pouring type | Gravity pouring |
| Pouring temperature | 1600 °C |
| Filling velocity | 5 m/s |
| Mold material | Dry quartz sand |
| Mold initial temperature | 25 °C |
| Ambient temperature | 25 °C |
| Total simulation steps | 5000 |
| Final temperature | 25 °C |
| Initial velocity | 1 m/s |
| Heat transfer coefficient (mold-core) | 1000 W/(m²·K) |
| Heat transfer coefficient (mold-casting) | 500 W/(m²·K) |
| Heat transfer coefficient (casting-core) | 500 W/(m²·K) |
Filling Process Analysis
At 0.647 seconds after pouring, the molten metal entered the sprue and split into two streams through the runner, then entered the cavity through the ingates. Due to gravity, the metal first flowed downward, filling the lower regions, and then spread to other parts. By 3.073 seconds, the cavity was nearly filled, but a small boss located at a certain position (where air entrapment was likely) remained incompletely filled. This could lead to a misrun defect in the sand casting foundry product. At 4.865 seconds, the mold was completely filled. The temperature at the end of filling was measured at around 1250 °C, indicating that the pouring velocity was relatively high, which might cause turbulent flow and gas entrainment.
The filling sequence showed that the metal front advanced unevenly. The central region filled faster than the extremities, and the thin sections near the bosses showed slower flow due to higher resistance. This is a common challenge in sand casting foundry where complex geometries cause non-uniform filling.
| Time (s) | Event |
|---|---|
| 0.647 | Metal enters cavity through ingates |
| 1.425 | Metal fills lower part of cavity |
| 3.073 | Nearly filled; incomplete boss region observed |
| 4.865 | Mold completely filled; temperature ~1250 °C |
Solidification Process Analysis
After filling, the temperature began to drop. At 5.741 seconds, solidification commenced. The solidification sequence was critical for predicting shrinkage defects. At 123.06 seconds, the flange riser solidified first, losing its ability to feed the adjacent heavy sections. The thin walls (e.g., the ear seat) solidified before the thicker central body, causing the risers to become ineffective. Moreover, the ingate region solidified earlier than the thick side plate, isolating it from the sprue feed. As a result, isolated liquid pools formed, leading to shrinkage porosity upon final solidification.
At 326.6 seconds, the side plate region, being the thickest, remained liquid while surrounding areas had solidified. This created a “hot spot” that could not be fed adequately. Finally, at 7326 seconds, the casting reached ambient temperature. The distribution of shrinkage defects predicted by ProCAST showed concentrated porosity in the flange, base, and side plate regions. The boss that had incomplete filling also exhibited some porosity.
| Time (s) | Event |
|---|---|
| 5.741 | Solidification begins |
| 123.06 | Flange riser solidifies; thin walls solidify |
| 326.6 | Side plate remains liquid; hot spot forms |
| 7326 | Complete solidification to ambient temperature |
| Location | Defect Type | Severity |
|---|---|---|
| Flange area | Shrinkage porosity | High |
| Base plate | Shrinkage cavity | High |
| Side thick plate | Shrinkage porosity | Moderate |
| Boss (near incomplete fill) | Misrun and porosity | Low |
Improvement of the Gating System
Based on the simulation results, I modified the gating system to address the identified shortcomings in the sand casting foundry process. The flange riser was enlarged and repositioned to sit on the flange face (which would be machined later), ensuring better feeding and easier removal. Additional risers were added on both sides of the front and rear ends to improve feeding of the ends. For the thick side plate, a riser was placed directly above it, and chills were applied on the inner wall to accelerate cooling and eliminate the hot spot. Furthermore, the cross-sectional area of the ingates was increased to delay their solidification and allow continued feeding from the runner system.
The modified design was re-meshed and simulated under the same process conditions. The filling process improved significantly, with no incomplete regions observed. The solidification sequence became more favorable: the risers solidified last, providing ample feed metal to the casting. The side plate region cooled faster due to the chills, eliminating the isolated liquid pool.
Results of the Improved Design
After implementing the improvements, the shrinkage defects nearly disappeared. The few remaining small porosities were located in areas designated for machining, thus they would be removed during post-processing. The predicted shrinkage porosity fraction was reduced to below 0.1% in all critical areas, compared to up to 0.97% in the initial design. The casting quality was significantly enhanced, meeting the requirements for high-pressure water jet applications.
| Location | Initial Defect Fraction | Improved Defect Fraction |
|---|---|---|
| Flange | 0.85% | 0.02% |
| Base | 0.92% | 0.01% |
| Side plate | 0.65% | 0.05% (in machining allowance) |
| Overall casting | 0.75% average | 0.03% average |
The simulation demonstrated that careful design of riser size, location, and the use of chills can effectively control solidification shrinkage in sand casting foundry. The numerical approach allowed me to optimize the process without costly trial-and-error experiments.
Conclusion
Through this numerical simulation study using ProCAST, I successfully identified the causes of shrinkage defects in a sand casting foundry process for a valve body. The initial design suffered from improper riser placement and insufficient feeding capacity, leading to hot spots and isolated liquid pools. By modifying the riser geometry, adding chills, and enlarging ingates, the solidification sequence was corrected, and defect levels were reduced to acceptable limits. This work highlights the power of numerical simulation in improving sand casting foundry quality, reducing development time, and lowering production costs. The optimized process has been validated for production, yielding defect-free castings that meet stringent performance requirements.