Numerical Simulation of Hot Tearing in Investment Casting for Automotive Turbocharger Turbines

The relentless pursuit of higher engine efficiency, power density, and reduced emissions has cemented turbocharging as a cornerstone of modern internal combustion engine technology. At the heart of this system lies the turbocharger turbine wheel, a component subjected to extreme temperatures and rotational stresses. To meet the demanding requirements of low inertia and high-temperature strength, nickel-based superalloys like IN713C are often employed, and the complex, thin-walled geometry of the turbine wheel is typically manufactured via the investment casting process. This process, while capable of producing intricate net-shape components, is prone to the formation of hot tears—intergranular cracks that occur in the late stages of solidification when the material’s strength and ductility are minimal. This article delves into a comprehensive numerical simulation study of the investment casting process for an automotive IN713C turbine, focusing on predicting and analyzing hot tearing defects. We will explore the governing thermomechanical models, present detailed simulation results, and systematically evaluate the influence of key process parameters on the propensity for hot tearing, concluding with actionable insights for process optimization.

A collection of complex, high-precision metal castings produced via investment casting, showcasing fine details and smooth surfaces.

Fundamentals of Hot Tearing and the Investment Casting Process

The investment casting process, also known as the lost-wax process, involves creating a ceramic shell mold around a wax pattern of the desired part. The wax is melted out, and molten metal is poured into the resulting cavity. During solidification, non-uniform cooling and geometrical constraints give rise to internal stresses. Hot tearing occurs in the mushy zone when the coherent solid skeleton is unable to accommodate thermally induced tensile strains, leading to failure through the liquid films remaining at grain boundaries. The susceptibility to hot tearing is influenced by alloy composition, casting geometry (which creates hotspots and stress concentrators), and the thermal conditions dictated by the investment casting process parameters such as pouring temperature and mold preheat.

Material: IN713C Superalloy

IN713C is a precipitation-strengthened nickel-based superalloy widely used for turbocharger wheels and gas turbine components due to its excellent high-temperature strength, creep resistance, and oxidation resistance. Its cast microstructure consists of a γ (gamma) nickel matrix strengthened by a high volume fraction of ordered γ’ (gamma prime) Ni3(Al, Ti) precipitates. Minor additions of carbon and boron form carbides and borides at grain boundaries, which can influence hot tearing behavior. The alloy’s thermal and mechanical properties, especially in the critical solidification range, are paramount for accurate simulation. Key thermophysical properties, such as fraction solid versus temperature and specific heat, were calculated using the Scheil-Gulliver model. The solidus and liquidus temperatures are approximately 1,196°C and 1,345°C, respectively. The mechanical properties required for stress simulation, including temperature-dependent Young’s modulus, yield strength, and Poisson’s ratio, were derived from material databases and literature.

Table 1: Nominal Chemical Composition of IN713C Superalloy (wt.%)
Ni Cr Mo Nb+Ta Al Ti C B Zr Fe Si Mn Co
Bal. 13.0 4.2 2.3 6.0 0.6 0.12 0.01 0.1 ≤1.8 0.4 0.23 ≤1.0

Numerical Modeling Framework for the Investment Casting Process

The simulation of the investment casting process requires coupled analysis of fluid flow, heat transfer, and stress/strain evolution. For this study, a commercial finite element software (ProCAST) was employed. The stress-strain response of the solidifying alloy was modeled using a thermal-elastoplastic constitutive model.

Thermal-Elastoplastic Stress Model

In this model, the material behavior is elastic until the yield stress is reached, after which it exhibits plastic deformation. The total strain increment ${d\varepsilon}$ is decomposed into elastic ${d\varepsilon^e}$, plastic ${d\varepsilon^p}$, and thermal ${d\varepsilon^T}$ components:
$$ {d\varepsilon} = {d\varepsilon^e} + {d\varepsilon^p} + {d\varepsilon^T} $$
The stress increment is related to the elastic strain increment via Hooke’s Law:
$$ {d\sigma} = [D^e] {d\varepsilon^e} $$
where $[D^e]$ is the elastic stiffness matrix. Combining these, the relationship in the elastoplastic regime is given by:
$$ {d\sigma} = [D^{ep}] ( {d\varepsilon} – {d\varepsilon^p} – {d\varepsilon^T} ) $$
where $[D^{ep}]$ is the elastoplastic stiffness matrix. A linear hardening rule is often assumed after yield:
$$ \sigma = \sigma_0 + H \bar{\varepsilon}^{pl} $$
where $\sigma_0$ is the initial yield stress, $H$ is the plastic modulus, and $\bar{\varepsilon}^{pl}$ is the equivalent plastic strain.

Hot Tearing Criterion: Hot Tearing Indicator (HTI)

To predict the location and severity of hot tearing, a strain-based failure criterion is commonly used. The Hot Tearing Indicator (HTI) model calculates the accumulated plastic strain in the vulnerable temperature range where the material has low strength (typically between a coherency temperature, $t_{coh}$, and the solidus temperature, $t_s$). The HTI is defined as:
$$ HTI = \int_{t_{coh}}^{t_s} \sqrt{\frac{2}{3} \, \dot{\varepsilon}^{p} : \dot{\varepsilon}^{p}} \, d\tau $$
where $\dot{\varepsilon}^{p}$ is the plastic strain rate tensor. A higher HTI value at a specific location indicates a greater propensity for hot tearing. This model effectively integrates the effects of thermal stress and the duration of material vulnerability during the investment casting process.

Case Study: Automotive Turbocharger Turbine

The subject of this simulation is a radial-inflow turbine wheel for an automotive turbocharger. The geometry features a central hub (approx. 28 mm thick) connected to ten thin, highly curved blades (approx. 0.7 mm thick). The significant difference in section thickness between the hub and the blade tips creates an inherent risk for differential cooling and stress concentration. A cluster of three turbines was simulated with a gating system consisting of a central downsprue, a pouring cup, and three ingates.

Table 2: Simulation Process Parameters for the Investment Casting Process
Parameter Values
Alloy IN713C
Liquidus Temperature 1,345 °C
Solidus Temperature 1,196 °C
Pouring Temperature (T_pour) 1,400, 1,450, 1,500, 1,550 °C
Mold Preheat Temperature (T_mold) 800, 850, 900 °C
Mold-Metal Interfacial Heat Transfer Coefficient 900 W/(m²·K)

Simulation Results: Filling, Solidification, and Stress

The filling stage of the investment casting process was rapid (≈1 second) due to the small mold size and large gating cross-section. The temperature field analysis revealed a steep thermal gradient. The thin blade tips cooled and solidified significantly faster than the massive central hub. Figure 4b (referring to the original paper’s data) shows that solidification progressed from the blade tips towards the hub, creating a large time gap between the solidification of these two regions. This sequential solidification, inherent to this geometry in the investment casting process, is a primary driver for the development of thermally induced tensile stresses in the early-solidifying, constrained blades.

The simulated stress distribution at a late stage of solidification (150 s) clearly showed stress concentration at the edges of the turbine blades. This location corresponded perfectly with the region where actual hot tears were observed in physically cast prototypes, validating the accuracy of the modeling approach for this investment casting process. The HTI map further confirmed that the highest propensity for hot tearing was localized precisely at these thin blade edges, where high tensile stress coincided with the material’s vulnerable state.

Parametric Analysis of the Investment Casting Process

A systematic parametric study was conducted to understand how key variables in the investment casting process influence hot tearing. A critical point (Point A) at the tip of a turbine blade was monitored for stress evolution and HTI under different conditions.

Effect of Pouring Temperature

With the mold preheat temperature fixed at 800°C, the impact of pouring temperature was investigated. The results revealed a non-monotonic relationship between pouring temperature and hot tearing risk.

Table 3: Effect of Pouring Temperature on Simulated Metrics at Point A (T_mold = 800°C)
Pouring Temp. (°C) Peak Thermal Stress (MPa) Time in Vulnerable Range (s) Hot Tearing Indicator (HTI) Trend
1400 ~47 Shortest Lower Baseline
1450 ~35 Longer Highest Increase
1500 ~25 Long Intermediate Decrease
1550 ~23 Longest Lowest Decrease

Analysis: Increasing the pouring temperature generally reduces the thermal gradient and thus the final residual thermal stress. However, it also prolongs the solidification time, extending the period the material spends in the brittle mushy zone (the “vulnerable range”). At 1,450°C, the extended vulnerable time dominates, leading to the highest HTI despite lower stress. At 1,500°C and 1,550°C, the significant reduction in thermal stress becomes the dominant factor, overriding the effect of longer vulnerable time and resulting in a lower HTI. This highlights a complex trade-off in the investment casting process.

Effect of Mold Preheat Temperature

With the pouring temperature fixed at 1,500°C, the effect of mold preheat was studied. Unlike pouring temperature, the effect was more straightforward.

Table 4: Effect of Mold Preheat Temperature on Simulated Metrics at Point A (T_pour = 1500°C)
Mold Temp. (°C) Peak Thermal Stress (MPa) Time in Vulnerable Range (s) Hot Tearing Indicator (HTI) Trend
800 ~25 Shorter Higher Baseline
850 ~20 Moderate Intermediate Decrease
900 ~18 Longer Lowest Decrease

Analysis: A higher mold preheat temperature reduces the initial chilling effect, lowering the cooling rate and thermal gradients throughout the casting. This leads to a direct and significant reduction in thermally induced stress. While the vulnerable time increases slightly, the stress reduction effect is overwhelmingly dominant in lowering the HTI. Therefore, within the studied range, increasing the mold preheat is a consistently beneficial strategy for mitigating hot tearing in this investment casting process.

Interaction of Process Parameters and Optimization

The interaction between pouring temperature and mold preheat is crucial for holistic process design. The following conclusions can be drawn from the simulation data:

  1. Thermal Stress: Thermal stress decreases with an increase in either pouring temperature or mold preheat temperature. The effect is most pronounced when one parameter is low; increasing the other then yields significant benefits. When both parameters are already high, further increases have diminishing returns on stress reduction.
  2. Hot Tearing Indicator (HTI): The HTI, which is the ultimate measure of risk, shows that increasing mold preheat consistently reduces hot tearing propensity. For pouring temperature, there exists an optimal window. Excessively low pouring temperatures cause high stress, while moderately high temperatures (e.g., 1,450°C) may maximize the vulnerable time. Very high pouring temperatures (e.g., 1,500-1,550°C) coupled with high mold preheat provide the best conditions for minimizing HTI.

Based on the numerical analysis of this investment casting process, an optimized set of parameters can be proposed. While a pouring temperature of 1,550°C with a 900°C mold gives the lowest theoretical HTI, practical considerations like grain growth, gas absorption, and mold-metal reaction may limit the maximum pouring temperature. Therefore, a robust optimized condition is a pouring temperature of 1,500°C with a mold preheat temperature of 900°C. This combination achieves a very low HTI while maintaining a safer margin below excessively high superheats.

Conclusion

This comprehensive numerical simulation study has successfully demonstrated the capability of thermomechanical modeling to predict hot tearing defects in the investment casting process of complex IN713C superalloy components. The simulation accurately identified the thin edges of turbine blades as critical failure sites, which was confirmed by experimental casting trials. The analysis revealed that hot tearing is governed by the interplay between thermally induced stress and the duration the material spends in a low-strength, brittle state. For the specific investment casting process of the automotive turbine:

  1. Hot tearing susceptibility at blade edges decreases monotonically with increasing mold shell preheat temperature due to dominant thermal stress reduction.
  2. The relationship with pouring temperature is non-linear, showing an initial increase in hot tearing risk (due to prolonged vulnerable time) followed by a decrease (due to dominant stress reduction) as temperature rises.
  3. Process optimization should target a combination of high mold preheat (e.g., 900°C) and a suitably high pouring temperature (e.g., 1,500°C) to minimize the integrated Hot Tearing Indicator.

This work underscores the value of numerical simulation as an indispensable tool for de-risking and optimizing the investment casting process, enabling the production of high-integrity, defect-free critical components like turbocharger turbines.

Scroll to Top