Numerical Simulation and Optimization of Investment Casting for Turbine Strut

As a practitioner in the field of advanced manufacturing, I have extensively studied the application of investment casting for producing critical aerospace components. The turbine rear frame strut, a key structural element in jet engines, presents significant challenges due to its complex geometry, thin-walled sections, and stringent quality requirements. In this article, I will delve into the numerical simulation and optimization of the investment casting process for this component, leveraging computational tools to predict and mitigate defects such as shrinkage porosity and misruns. Investment casting, known for its ability to produce intricate shapes with high dimensional accuracy, is the focal point of this investigation. Through detailed analysis and iterative refinement, an optimized investment casting process was developed, ensuring the production of defect-free struts that meet rigorous industry standards.

The turbine rear frame strut is typically made from high-temperature alloys like K4169 and features a composite structure comprising a load-bearing lug, a concave float plate, and airfoil-shaped vanes with internal ribs. The geometry is characterized by varying wall thicknesses, from as thin as 2.5 mm in the float plate to 31 mm in the lug section, and includes undercuts and deep recesses up to 64 mm. Such complexity necessitates a meticulous investment casting approach to avoid common issues like gas entrapment, turbulent filling, and solidification-related defects. The success of investment casting hinges on controlling the molten metal flow, heat transfer, and solidification patterns, which can be effectively simulated using software like ProCAST. This study underscores the importance of numerical simulation in optimizing investment casting parameters, reducing trial-and-error cycles, and enhancing product quality.

To understand the initial investment casting process, a three-dimensional model of the strut was created, incorporating a basic gating system with feeders attached to the thick lug sections. The simulation setup involved discretizing the domain into finite elements, with dual nodes at the interface between the casting and ceramic shell to account for temperature gradients. Key parameters for the investment casting simulation are summarized in Table 1, which includes material properties, boundary conditions, and process settings. The alloy used was K4169, with a molten temperature of 1,450°C, while the shell was made of mullite, preheated to 800°C. Vacuum cooling was applied, and the interfacial heat transfer coefficient was set to 500 W/m²·K. The filling percentage was monitored up to 100%, with a focus on the initial stages where defects often originate.

Table 1: Key Parameters for the Investment Casting Simulation
Parameter Value Description
Alloy Material K4169 High-temperature nickel-based alloy
Shell Material Mullite Ceramic shell for investment casting
Cooling Method Vacuum Cooling Enhances solidification control
Shell Thickness 8 mm Uniform thickness around the pattern
Pouring Temperature 1,450°C Initial temperature of molten alloy
Shell Preheat Temperature 800°C Reduces thermal shock in investment casting
Furnace Dimensions φ2 m × 2.4 m Size of the investment casting furnace
Interfacial Heat Transfer Coefficient 500 W/m²·K Governs heat flux at casting-shell interface
Filling Percentage 100% Complete mold filling in simulation

The fluid dynamics during filling in investment casting can be described by the Navier-Stokes equations, which account for momentum, mass, and energy conservation. For incompressible flow, the continuity and momentum equations are:

$$ \nabla \cdot \mathbf{u} = 0 $$

$$ \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} = -\frac{1}{\rho} \nabla p + \nu \nabla^2 \mathbf{u} + \mathbf{g} $$

where \(\mathbf{u}\) is the velocity vector, \(p\) is pressure, \(\rho\) is density, \(\nu\) is kinematic viscosity, and \(\mathbf{g}\) is gravitational acceleration. In investment casting simulations, these equations are coupled with energy transport to model temperature changes:

$$ \rho c_p \left( \frac{\partial T}{\partial t} + \mathbf{u} \cdot \nabla T \right) = \nabla \cdot (k \nabla T) + Q $$

Here, \(T\) is temperature, \(c_p\) is specific heat, \(k\) is thermal conductivity, and \(Q\) represents heat sources such as latent heat release during solidification. The pressure evolution within the mold cavity, crucial for detecting gas entrapment in investment casting, is derived from the pressure-velocity coupling, often solved using algorithms like SIMPLE. The reference pressure \(P_{ref}\) in the cavity is computed as:

$$ P_{ref} = |P_{calc} – P_{set}| $$

where \(P_{calc}\) is the calculated pressure from fluid dynamics, and \(P_{set}\) is a preset threshold. Values exceeding 0.2 MPa indicate potential “gas suffocation,” a common issue in investment casting that leads to misruns.

Simulating the initial investment casting process revealed significant problems. During filling, at 20% completion, molten metal flowed from the lugs into the float plate and vanes, with the vane sections partially filled and pressure below 0.07 MPa. However, by 40% filling, the pressure in unfilled regions of the vanes and lugs surged above 0.22 MPa, indicating severe gas entrapment. Turbulent flow was observed, with average velocities reaching 2.2 m/s at the vane-plate junction, exacerbating defect formation. At 60% filling, a large area in the float plate remained unfilled with pressure over 0.27 MPa, confirming gas blockage. Even at 80%, high-pressure zones persisted, suggesting that the initial investment casting design was prone to misruns due to inadequate venting.

Defect prediction using ProCAST’s shrinkage module showed porosity and shrinkage cavities in critical areas. The vane sections, especially at the edges, centers, and rib reinforcements, exhibited concentrated defects, as quantified in Table 2. This resulted from rapid filling causing temperature gradients, with thin sections solidifying faster than thicker ones, leading to inadequate feeding. The float plate also had two shrinkage zones due to turbulent flow, while the lugs, with ample feeders, remained defect-free. The solidification sequence, governed by the Fourier heat equation, highlights the importance of thermal management in investment casting:

$$ \frac{\partial T}{\partial t} = \alpha \nabla^2 T $$

where \(\alpha = k/(\rho c_p)\) is thermal diffusivity. Regions with low diffusivity, like the vane centers, solidified last, creating hot spots and shrinkage defects. The Niyama criterion, often used to predict porosity in investment casting, is expressed as:

$$ G / \sqrt{\dot{T}} \leq C $$

where \(G\) is temperature gradient, \(\dot{T}\) is cooling rate, and \(C\) is a material constant. Low values indicate a high risk of microporosity, which aligned with the simulated defect locations.

Table 2: Simulated Defect Distribution in Initial Investment Casting Process
Component Region Defect Type Severity Level Probable Cause
Vane Edges Shrinkage Porosity Medium Rapid cooling and turbulent flow
Vane Centers Shrinkage Cavities High Low thermal diffusivity and poor feeding
Rib Reinforcements Microporosity Medium Geometric constraints in investment casting
Float Plate Shrinkage Zones Low Flow instability during filling
Load-Bearing Lugs None None Adequate feeder design

To optimize the investment casting process, several modifications were implemented. First, the venting system was enhanced from one to three exhaust channels in the float plate area to alleviate gas entrapment. The pouring time was increased by 1.5 times to reduce flow turbulence, and the strut was reoriented by 15° relative to the exhaust plane to promote smoother filling. Additionally, a larger sandbox measuring 440 mm × 360 mm × 600 mm was used, and 20 mm thick iron sand was added around the vane centers and ribs to modify the solidification sequence. These adjustments aimed to create a more controlled investment casting environment, minimizing defects. The optimized gating system, as modeled in ProCAST, incorporated these changes, with simulation parameters adjusted accordingly. The key improvements are summarized in Table 3, emphasizing the role of venting and thermal control in successful investment casting.

Table 3: Optimization Measures for Investment Casting Process
Optimization Aspect Initial Design Optimized Design Impact on Investment Casting
Exhaust Channels 1 channel 3 channels Reduces gas pressure and prevents suffocation
Pouring Time Base time T 1.5T Decreases flow velocity and turbulence
Strut Orientation 0° rotation 15° rotation Enhances metal flow and venting efficiency
Sandbox Size Standard size 440 mm × 360 mm × 600 mm Provides better insulation and cooling control
Thermal Modifiers None Iron sand added Alters solidification pattern to reduce hot spots
Gating System Simple feeders Revised with multiple vents Improves feeding and reduces shrinkage in investment casting

The optimized investment casting process was simulated under the same conditions as the initial run. At 20% filling, molten metal flowed steadily from one lug into the float plate and vanes, with unfilled vane areas maintaining external connectivity and pressures below 0.7 MPa. The average flow velocity in the vanes was reduced to 0.9 m/s, 56% lower than before, indicating laminar flow. By 40% filling, the vanes were completely filled, with only minor areas in the float plate and one lug remaining, and pressures stayed under 0.7 MPa, showing no gas entrapment. At 60% filling, the entire strut was filled, with pressures ranging from 0.7 to 1.3 MPa only in exhaust channels and feeders, which are non-critical. The filling process remained stable throughout, demonstrating the effectiveness of the optimization in investment casting.

Defect analysis of the optimized investment casting simulation revealed no shrinkage porosity or cavities in the strut body. Minor defects were confined to the exhaust channels and feeder sections, which are typically removed during post-processing. The solidification sequence, now more uniform, prevented hot spot formation. This can be quantified using the solidification time \(t_s\), derived from Chvorinov’s rule, commonly applied in investment casting:

$$ t_s = B \left( \frac{V}{A} \right)^n $$

where \(V\) is volume, \(A\) is surface area, \(B\) is a mold constant, and \(n\) is an exponent (typically 2). For the optimized design, the modulus \(V/A\) was balanced across regions, reducing variations in \(t_s\) and promoting directional solidification. The improved thermal profile is reflected in the Fourier number \(Fo\), a dimensionless parameter for heat conduction:

$$ Fo = \frac{\alpha t}{L^2} $$

where \(L\) is characteristic length. Higher \(Fo\) values in critical areas indicated better heat dissipation, minimizing defects. The success of this investment casting optimization underscores the value of numerical simulation in predicting and enhancing process outcomes.

To validate the optimized investment casting process, a production trial was conducted. Five prototype struts were cast using the revised parameters, and all were inspected via digital radiography according to ASTM E192 standards. The results confirmed the absence of internal defects such as gas pores, inclusions, or shrinkage, meeting the stringent requirements for aerospace components. Subsequently, a batch of 40 struts was produced using the same investment casting methodology, with all units passing quality checks. This real-world verification demonstrates the accuracy of the numerical simulations and the robustness of the optimized investment casting process. The ability to consistently produce defect-free struts highlights the transformative impact of simulation-driven design in investment casting, reducing scrap rates and accelerating development cycles.

In conclusion, this study exemplifies the integration of numerical simulation into the investment casting workflow for complex turbine components. By analyzing fluid flow, pressure evolution, and solidification patterns, initial design flaws were identified and rectified through targeted optimizations. The enhanced venting, adjusted pouring parameters, and modified thermal management collectively eliminated defects, yielding a reliable investment casting process. The methodologies discussed here—from mathematical modeling to practical validation—provide a framework for advancing investment casting in high-stakes industries. Future work could explore multi-scale simulations or machine learning for further refinement, but the present results affirm that investment casting, when coupled with computational tools, can achieve unparalleled precision and quality in manufacturing critical aerospace parts.

The broader implications for investment casting are significant. As industries demand lighter, stronger, and more intricate components, the role of simulation in investment casting will only grow. By preemptively addressing issues like gas entrapment and shrinkage, manufacturers can reduce costs and improve sustainability. This case study on the turbine strut serves as a testament to the power of investment casting innovation, driven by data and technology. I encourage practitioners to embrace these approaches, fostering a new era of excellence in investment casting for aerospace and beyond.

Scroll to Top