Optimization of Lost Foam Casting Process for Automotive Brake Hubs Based on Numerical Simulation

The lost foam casting process has established itself as a pivotal, environmentally conscious manufacturing technique within the foundry industry. Characterized by its utilization of expandable polystyrene (EPS) patterns that vaporize upon contact with molten metal, this process offers exceptional benefits for producing complex, near-net-shape components with high dimensional accuracy and minimal machining allowances. Its application is particularly advantageous for high-volume automotive parts, where consistency, cost-effectiveness, and geometric fidelity are paramount. However, the inherent nature of the lost foam casting process also presents unique challenges, primarily concerning defect formation during solidification, such as shrinkage porosity and cavities, which can critically compromise the mechanical integrity and service life of the final casting.

This article provides a comprehensive, first-person perspective on the systematic optimization of the lost foam casting process for a critical automotive component—the brake hub. Brake hubs are subjected to extreme thermal and mechanical stresses, demanding high strength, wear resistance, and fatigue performance. Any internal defects like shrinkage are unacceptable. We will delve into a detailed methodology that employs advanced numerical simulation as the core tool for diagnosis and optimization. The discussion will encompass the fundamental governing equations of the process, a step-by-step analysis of a specific case study, the formulation and validation of an optimized lost foam casting process scheme, and the resultant improvements in casting quality.

Numerical Simulation Fundamentals for Lost Foam Casting

The fidelity of optimizing the lost foam casting process via simulation hinges on accurately modeling the coupled physics of fluid flow, heat transfer, and phase change. The process is governed by a set of conservation equations, solved computationally to predict the behavior of the molten metal from the moment it enters the flask until complete solidification.

The flow of the molten metal during mold filling is described by the Navier-Stokes equations, which for an incompressible fluid (assuming negligible density variation during filling) can be expressed in their conservation form. The general transport equation for a scalar property $\phi$ (which can represent velocity components, temperature, etc.) in a fluid flow is given by:

$$
\frac{\partial (\rho \phi)}{\partial t} + \nabla \cdot (\rho \vec{v} \phi) = \nabla \cdot (D_{\phi} \nabla \phi) + S_{\phi}
$$

Where:
$\rho$ is the fluid density,
$\vec{v}$ is the velocity vector,
$t$ is time,
$D_{\phi}$ is the diffusion coefficient for $\phi$,
$S_{\phi}$ is the source term for $\phi$.

Specifically, when $\phi$ represents a velocity component, this leads to the momentum equations. Coupled with this is the continuity equation for mass conservation:

$$
\nabla \cdot (\rho \vec{v}) = 0
$$

In the context of the lost foam casting process, a critical additional phenomenon is the thermal decomposition of the EPS pattern. The heat transfer model must account for the energy consumed in vaporizing the foam. A simplified form of the energy equation considering this interaction can be written as:

$$
\rho_L c_{p,L} \frac{\partial T}{\partial t} + \rho_L c_{p,L} (\vec{v} \cdot \nabla T) = \nabla \cdot (k_L \nabla T) + \rho_P L_P \frac{\partial f_s}{\partial t}
$$

Where the subscript $L$ denotes liquid metal properties and $P$ denotes foam pattern properties:
$c_{p}$ is the specific heat capacity,
$k$ is the thermal conductivity,
$L_P$ is the latent heat of vaporization/decomposition of the foam,
$f_s$ is the fraction of foam decomposed.

Finally, the solidification process is predicted using models that track the evolution of the solid fraction $f_s^{metal}$ in the metal, often employing a Scheil-Gulliver model or a lever rule for equilibrium solidification. The release of latent heat $L_L$ during metal solidification is a crucial source term in the energy equation:

$$
S_{latent} = \rho_L L_L \frac{\partial f_s^{metal}}{\partial t}
$$

The boundary conditions are vital. Key interfacial heat transfer coefficients (HTC) must be defined, such as those between the sand mold and the coating, the coating and the decomposing foam, and the metal and the coating. A summary of typical simulation parameters is provided below.

Table 1: Key Parameters for Numerical Simulation of Lost Foam Casting
Parameter Category Specific Parameter Value / Specification
Material Properties Casting Alloy Ductile Iron QT600-3
Pattern Material Expandable Polystyrene (EPS)
Molding Sand 35-50 mesh Silica Sand
Process Conditions Pouring Temperature 1450°C
Mold Temperature 25°C
Pouring Time 26.4 s
Vacuum Pressure -0.4 bar
Interfacial Conditions Sand-Foam HTC 100 W/m²·K
Sand-Casting HTC 500 W/m²·K
Coating Thickness & Material 1.0 mm, Silica-based refractory coating

Case Study: Brake Hub Casting – Initial Process Analysis

The subject component is a ductile iron brake hub with an outer diameter of 341 mm, a height of 231 mm, and an approximate mass of 31.2 kg. Its geometry is essentially a cylindrical shell with a mounting flange. While most walls are relatively uniform (~11 mm), a significant thicker section (up to ~75 mm) exists at the top closed-end of the hub, creating a pronounced thermal mass. A one-pattern-four-casting gating system was designed for productivity. To promote smooth filling and minimize turbulence, a bottom-gating system was employed, with ingates positioned at the base of each hub cavity. A top riser was initially placed on each casting to feed shrinkage.

The simulation of the initial lost foam casting process scheme revealed valuable insights. The filling sequence, completed in 26.4 seconds, was indeed smooth and progressive from the bottom upward, confirming the gating design was adequate for defect-free mold filling. The critical analysis, however, lies in the solidification sequence and temperature field evolution.

The simulated solidification progression clearly identified the thick section at the closed-end of the hub as the last region to solidify. As thinner sections and even parts of the gating system solidified earlier, this thick region became an isolated liquid pool. According to fundamental solidification theory, shrinkage defects form in regions that solidify last if they are not adequately fed by liquid metal. The Niyama criterion, a widely used indicator for predicting shrinkage porosity in castings, is often calculated as:

$$
N_y = \frac{G}{\sqrt{\dot{T}}}
$$

Where $G$ is the temperature gradient and $\dot{T}$ is the cooling rate at the solidus front. Low values of the Niyama criterion correlate with a high probability of microporosity formation. In the initial design, simulation results predicted a high concentration of shrinkage porosity precisely in this thermally isolated thick section, as the riser proved insufficient to maintain a feeding path throughout its solidification.

Table 2: Simulated Results for Initial Casting Design
Aspect Analyzed Simulation Observation Implication for Defect Formation
Filling Pattern Smooth, bottom-up progression. No misruns or cold shuts predicted. Gating system design is functionally sound for filling.
Solidification Sequence Thin walls solidify first. The thick top section solidifies last, becoming isolated from the riser. Creates a classical “hot spot” condition prone to shrinkage.
Shrinkage Prediction (Niyama) High defect probability (low Niyama values) concentrated in the thick top section of the hub. The existing riser is ineffective in feeding this region due to premature feeding path interruption.
Temperature Gradient (G) Low gradient in the thick section relative to surrounding areas. Promotes equiaxed, pasty solidification unfavorable for directional feeding.

Process Optimization Strategy and Redesign

The analysis of the initial lost foam casting process simulation pinpointed the root cause: an unfavorable temperature field leading to an isolated hot spot. The optimization goal was therefore to alter the solidification sequence to achieve directional solidification, progressing from the casting extremities toward a dedicated feeding source (riser). Two key modifications were implemented to the lost foam casting process scheme.

1. Modification of Casting Geometry (Functional Weight Addition): Recognizing that the top inner wall of the hub underwent machining, its thickness was intentionally increased in the pattern. This added “functional” mass serves as a chiller in reverse; it actually makes the region directly below the riser solidify slightly later than the problematic thick section of the hub body. This subtle change helps reorient the thermal gradients.

2. Redesign of the Feeding System: The small top riser was replaced with a substantial annular riser that encircles the top of the hub casting. This provides a much larger reservoir of hot metal with significant thermal mass. The combined effect of the geometry change and the annular riser is designed to establish a clear thermal hierarchy: the thickest part of the hub body should now solidify before the junction between the hub and the riser, ensuring the riser remains liquid to feed the shrinkage in the casting.

The principle can be conceptualized through Chvorinov’s rule, which states the solidification time $t_s$ is proportional to the square of the volume-to-surface-area ratio $(V/A)^2$:

$$
t_s = C \left( \frac{V}{A} \right)^2
$$

Where $C$ is the mold constant. By increasing the $(V/A)$ ratio of the riser (making it larger and more voluminous) and carefully adjusting the $(V/A)$ of the hub’s top section, the solidification time order can be controlled: $t_{s,riser} > t_{s,hub-thick-section}$.

Simulation and Validation of the Optimized Process

The redesigned lost foam casting process was subjected to the same rigorous numerical simulation. The results demonstrated a fundamental improvement in the thermal field. The solidification sequence now showed a more progressive pattern, with the annular riser remaining liquid longest. Most critically, the area of predicted shrinkage porosity was dramatically reduced within the brake hub itself and was largely confined to the central sprue and the riser—regions intended to be sacrificial and subsequently removed.

The efficacy of the optimization can be quantified by comparing key simulation metrics before and after the changes. The improved temperature gradient $G$ in the critical hub section and the shift in the shrinkage location confirm the success of the strategy.

Table 3: Comparison of Initial vs. Optimized Process Simulation Results
Evaluation Metric Initial Process Scheme Optimized Process Scheme Improvement
Last-to-Solidify Region Thick section of hub body. Annular riser and central sprue. Shrinkage moved from casting to feeder.
Shrinkage Prediction in Casting High probability in hub top. Negligible probability. Elimination of internal casting defect.
Directionality of Solidification Pasty, isolated hot spot. Directional from hub base to top riser. Established thermally favorable gradient.
Riser Efficiency Low (solidified before casting hot spot). High (remained liquid to feed casting). Feeding mechanism became functional.

The ultimate validation of any lost foam casting process optimization is physical production. Castings produced using the optimized process parameters and pattern design were inspected using non-destructive X-ray testing. The results confirmed the simulation predictions: the brake hubs were free from internal shrinkage porosity and cavities. The defect was successfully relegated to the riser and sprue, which are removed during finishing. This led to a significant increase in the casting yield and合格率.

Conclusion

This detailed exploration underscores the transformative power of numerical simulation in advancing the lost foam casting process. By applying computational tools to model the complex interplay of fluid dynamics, heat transfer, and phase changes, it is possible to move from a trial-and-error approach to a precise, science-based optimization methodology. The case of the automotive brake hub clearly illustrates the process: identifying a solidification-related defect (shrinkage) via simulation, diagnosing its root cause in the temperature field, and implementing targeted design changes—such as functional weight addition and riser redesign—to alter the thermal gradients and solidification sequence.

The optimized lost foam casting process scheme achieved directional solidification, effectively transferring shrinkage from the critical casting into the sacrificial feeding system. This was conclusively verified both in simulation and in physical production, resulting in sound, high-quality castings. The methodology presented is universally applicable, providing a robust framework for optimizing the lost foam casting process for a wide array of complex components, thereby enhancing reliability, reducing scrap, and improving the overall efficiency and sustainability of foundry operations.

Scroll to Top