Advanced Process Simulation of Low Pressure Die Cast A356 Aluminum Automotive Wheels—Part II Modeling Methodology and Validation

: This manuscript presents an advanced modeling methodology developed to accurately simulate the temperature ﬁeld evolution in the die and wheel in an industrial low-pressure die casting (LPDC) machine employed in the production of A356 automotive wheels. The model was developed in the commercial casting simulation platform ProCAST for a production die operating under production conditions. Key elements in the development of the model included the deﬁnition of the resistance to heat transfer across the die / casting interfaces and die / water-cooling channel interfaces. To examine the robustness of the modeling methodology, the model was applied to simulate production and non-production process conditions for a die cooled by a combination of water and air-cooling (Die-A), and to a second die for a di ﬀ erent wheel geometry (Die-B) utilizing only water cooling for production conditions. In each case, the model predictions with respect to in-die and in-wheel temperature evolution were compared to industrially derived thermocouple (TC) data, and were found to be in good agreement. Once tuned to the process conditions for Die-A operating under production conditions, no further tuning of the die / casting interface resistance was applied. Additionally, the model results, in terms of the prediction of pockets of solid encapsulated liquid, were used to compare to x-ray images of wheels. This comparison indicated that the model was able to predict clusters of porosity associated with encapsulated liquid with an equivalent radius of ~27 mm.


Introduction
Low-pressure die casting (LPDC) is the primary process adopted to manufacture load-bearing aluminum automotive parts thanks to its ability to produce parts with a good balance of cost-effectiveness and performance [1][2][3][4]. It is a complex process that requires careful control of the filling and solidification stages in order to avoid the occurrence of defects that degrade component performance. Computer-based simulation, as a complement to experience, has become a crucial aspect in the development and optimization of both die-tooling and process operational parameters [5]. Driven by the increasingly strict requirements on weight, mechanical performance, geometric tolerance, microstructure and pore defects, increased emphasis is being placed on computer-based simulation as a means of achieving the necessary quality. In recognition of this, original equipment

Model Development
The model was developed using the commercial software package ProCAST (A trademark of the ESI Group), which is one of the most commonly used packages in the industry for the simulation of casting processes.

Geometry and Mesh
The components included in the computational domain are as follows: top die, side dies, bottom die and wheel. To reduce the computational size of the problem, the transfer tube and die support structures-e.g., top and bottom plates-were not included in the domain, since they are not expected to have a significant influence on the accuracy of the analysis. For the same reason, the venting features and ejection pins were also removed from the die structure. In addition, to take advantage of the mirror symmetry of the die and wheel structures, half of the die and wheel geometry has been used as the computational domain.
The strategy for meshing was to seek a balance between model accuracy and computational efficiency. In general, in regions that are relatively distant from the wheel and/or where a small temperature gradient is expected, a relatively coarse mesh was used. In contrast, in the wheel and regions where a high temperature gradient is expected, a fine mesh was employed. The following specific approaches were adopted: • In the rim of the wheel that has a relatively thin cross section, a 2 mm mesh size was selected to have at least 3-4 elements across the rim; • At the locations where strong heat transfer exists, i.e., in proximity to the cooling channels, a mesh size of 2 mm was selected; • A gradual increase of mesh size was implemented moving away from the wheel/die interface and cooling channels; • A default mesh size of 10 mm was used for the balance of the die geometry.
The resulting mesh for the wheel and die (Die-A) has 378,608 2D (surface) elements and 2,137,114 3D (body) elements (see the wheel mesh in Figure 1). The calculations were carried out on four processors of an Intel Xeon CPU (2.4 GHz). Computation time was approximately 1/2 h to run a single thermal-only analysis and 2 h to run a coupled thermal-fluid flow analysis.
Metals 2020, 10, x FOR PEER REVIEW 3 of 25 of the mirror symmetry of the die and wheel structures, half of the die and wheel geometry has been used as the computational domain. The strategy for meshing was to seek a balance between model accuracy and computational efficiency. In general, in regions that are relatively distant from the wheel and/or where a small temperature gradient is expected, a relatively coarse mesh was used. In contrast, in the wheel and regions where a high temperature gradient is expected, a fine mesh was employed. The following specific approaches were adopted: • In the rim of the wheel that has a relatively thin cross section, a 2 mm mesh size was selected to have at least 3-4 elements across the rim; • At the locations where strong heat transfer exists, i.e., in proximity to the cooling channels, a mesh size of 2 mm was selected; • A gradual increase of mesh size was implemented moving away from the wheel/die interface and cooling channels; • A default mesh size of 10 mm was used for the balance of the die geometry.
The resulting mesh for the wheel and die (Die-A) has 378,608 2D (surface) elements and 2,137,114 3D (body) elements (see the wheel mesh in Figure 1). The calculations were carried out on four processors of an Intel Xeon CPU (2.4 GHz). Computation time was approximately 1/2 h to run a single thermal-only analysis and 2 h to run a coupled thermal-fluid flow analysis.

Governing Equations
The governing equations describing energy, mass and momentum conservation in the computational domain are the classic Navier-Stokes Equations, which are given below: where (kg/m 3 ) is the density, (s) is the time, (m/s) is the velocity of the fluid, (Pa) is the pressure, (Pa) is the shear stress, (m/s 2 ) is the gravitational acceleration, (J/(kg K)) is the

Governing Equations
The governing equations describing energy, mass and momentum conservation in the computational domain are the classic Navier-Stokes Equations, which are given below: where ρ (kg/m 3 ) is the density, t (s) is the time, u (m/s) is the velocity of the fluid, p (Pa) is the pressure, τ (Pa) is the shear stress, g (m/s 2 ) is the gravitational acceleration, C P (J/(kg K)) is the specific heat at constant pressure, T (K) is the temperature, k (W/(m K)) is the thermal conductivity, and . Q (J/(m 3 s)) is the volumetric source term of heat, which incorporates the latent heat associated with the in-wheel solidification in this model. Note: the calculation of . Q will be presented in Section 2.3.1. The group of equations was solved subject to the specification of appropriate material properties, initial conditions and boundary conditions. Please note that the momentum conservation equation (Equation (2)) is only activated during the filling stage in ProCAST to save computation time. This means that the model assumes no fluid flow once the die cavity is fully filled and will no longer generate fluid flow-related results, i.e., fluid velocity. The modeling methodology was initially developed using the production process conditions for Die-A. This process involved some "tuning" of the internal and external interface conditions, which are described below in Section 2.5. The model is then applied to a second set of process conditions for Die-A (non-production) without retuning the model. The only changes to the model are related to the timing of activation and deactivation of the die cooling elements, and are applied in a consistent manner with the actual die machine process timings. Finally, the model is applied to a second die geometry, Die-B, with a third set of process timings (cooling element activation/deactivation and pressure fill curve), without retuning.

Materials and Properties
The materials used in the model are listed in Table 1, and their thermophysical properties are presented graphically in the following subsections. For those components that remain solid throughout the process-e.g., the die structure-only the properties required for solving the energy conservation equation are needed, which are density, thermal conductivity and specific heat. In contrast, the wheel material initially enters the die as a liquid, which then is transformed to a solid. Therefore, additional properties related to solving the momentum and continuity conservation equations, and dealing with the latent heat of solidification in the energy conservation equation, are needed for the wheel material A356. These include viscosity, solid fraction vs. temperature in the semi-solid region, and the latent heat of solidification. The properties of the A356 used in the model are temperature-dependent and presented in Figure 2 [18,19]. Please note that, as indicated in Figure 2e, the solidus, eutectic transformation temperature and liquidus are 500 • C, 560 • C and 614 • C, respectively. Additionally, the latent heat associated with the in-wheel solidification (the heat source term, . (3)) is calculated by Equation (4).

Q, in Equation
where . Q (J/(m 3 s)) is the volumetric rate of heat release due to solidification, ρ (kg/m 3 ) is the density, L (J/kg) is the specific latent heat (390 kJ/kg was used in the model),

H13 Tool Steel-Top Die and Bottom Die
H13 tool steel is used for the top and bottom dies, and the properties are given in Figure 3 [15,18].     The side dies are made from 35CrMo tool steel-the Chinese designation, which is equivalent to the German designation, 34CrMo4. The chemical compositions of both alloys are given in Table 2. It can be seen that the two alloys have negligible variation in composition, and therefore, the temperature-dependent thermal conductivity and specific heat capacity of 34CrMo4 are used in the model-see Figure 4, as those of 35CrMo are unavailable [20]. For the density, data in the appropriate temperature range is not available for both 34CrMo4 and 35CrMo. However, the density of alloys will not be significantly affected by small variations in the composition of alloying elements, and thus density data for H13 is used in the model. Note: Table 2 gives the composition of H13. The comparison of the densities of H13 and 34CrMo4 at room temperature only shows an ignorable difference-less than 1%. The side dies are made from 35CrMo tool steel-the Chinese designation, which is equivalent to the German designation, 34CrMo4. The chemical compositions of both alloys are given in Table 2. It can be seen that the two alloys have negligible variation in composition, and therefore, the temperature-dependent thermal conductivity and specific heat capacity of 34CrMo4 are used in the model-see Figure 4, as those of 35CrMo are unavailable [20]. For the density, data in the appropriate temperature range is not available for both 34CrMo4 and 35CrMo. However, the density of alloys will not be significantly affected by small variations in the composition of alloying elements, and thus density data for H13 is used in the model. Note: Table 2 gives the composition of H13. The comparison of the densities of H13 and 34CrMo4 at room temperature only shows an ignorable difference-less than 1%.

Initial Conditions
In the commercial process, the die structure is typically pre-heated to between 300 and 500 °C using propane burners before the first casting cycle. However, due to the cyclic nature of the LPDC process, the die temperature will evolve over the first few casting cycles until eventually reaching the cyclic steady state for a given set of process conditions. The approach adopted in the present work is that the model is initially run for nine thermal-only cycles followed by six thermal and fluid flow cycles. For the initial thermal-only cycle, the die is set to 400 °C. In each subsequent cycle, the die structure is set to the temperature field calculated from the previous run. In the thermal only analyses, the calculation is executed with the die cavity full and at a uniform temperature of 700 °C at the beginning of each cycle. For the thermal-fluid flow analyses, the metal inlet temperature to the die is set to 700 °C and die filling is simulated together with the heat transfer during die filling. The approach described above greatly reduces the time required to achieve cyclic steady state as a thermal-only cycle takes ~25 min to run, whereas a coupled thermal and fluid flow cycle takes ~110 min. Cyclic steady state is assumed to have been achieved when the cycle-to-cycle thermal field in the die at the beginning of a cycle does not vary more than 5 °C.

Initial Conditions
In the commercial process, the die structure is typically pre-heated to between 300 and 500 • C using propane burners before the first casting cycle. However, due to the cyclic nature of the LPDC process, the die temperature will evolve over the first few casting cycles until eventually reaching the cyclic steady state for a given set of process conditions. The approach adopted in the present work is that the model is initially run for nine thermal-only cycles followed by six thermal and fluid flow cycles. For the initial thermal-only cycle, the die is set to 400 • C. In each subsequent cycle, the die structure is set to the temperature field calculated from the previous run. In the thermal only analyses, the calculation is executed with the die cavity full and at a uniform temperature of 700 • C at the beginning of each cycle. For the thermal-fluid flow analyses, the metal inlet temperature to the die is set to 700 • C and die filling is simulated together with the heat transfer during die filling. The approach described above greatly reduces the time required to achieve cyclic steady state as a thermal-only cycle takes~25 min to run, whereas a coupled thermal and fluid flow cycle takes~110 min. Cyclic steady state is assumed to have been achieved when the cycle-to-cycle thermal field in the die at the beginning of a cycle does not vary more than 5 • C.

Boundary Conditions
Since the model involves both heat transport and fluid flow, it requires the definition of two types of boundary conditions: (1) those associated with the thermal aspects of the analysis, and (2) those associated with the fluid flow aspects.
where q (W/m 2 ) is the heat flux across the boundary, h (W/(m 2 K)) is the heat transfer coefficient (HTC), and T surf (K) and T ref (K) are the boundary surface temperature and the reference temperature, respectively. Note: the HTC represents the inverse of the resistance to heat transport across and at the interface between two solids or through a boundary layer between a solid and a fluid. There are two types of boundaries: one type is for interfaces where the boundary surface is internal and in contact with the surface of another component of the analysis domain, and the second type are interface is where the boundary is external and there is no contact with any other components. For the internal interfaces, the reference temperature, T ref , in Equation (5) is the temperature of the adjacent surface in contact at the location being processed in the analysis, which varies with time and is calculated by the model. For the external surfaces, T ref is set to the environment temperature, which is a user input. Please note that when the die is opened and the wheel is removed close to the end of a cycle, the casting/die interface is shifted from an internal type to an external type for the duration of the analysis of the cycle.
Heat Transport Across the Casting/Die Interfaces-The HTCs applied at the casting/die interfaces were determined through a combination of trial-and-error fitting to the production process data and experience, using the following principles: • The HTCs applied at all casting/die interfaces are temperature-dependent, following the general behavior shown in Figure 5. When the temperature is above the eutectic transformation point of A356, and the metal is fully, or mainly, liquid, the HTCs are set to relatively high values since the liquid metal is in close-to-perfect contact with the die. In contrast, when the temperature is below the solidus and the wheel is fully solidified, the contact degrades, resulting in lower values of HTCs due to the formation of solid/solid interfaces. Depending on the interface location, it is possible for a gap to form between the casting and the die surface, which leads to lower HTC values. In the region between the eutectic transformation point and solidus, the HTCs are assumed to vary linearly.

•
Values of the HTCs for the different casting/die interfaces are given in Table 3. For the casting/bottom die interface, the HTC for the liquid/solid regime of behavior is higher compared to the other two interfaces. The rationale is that the bottom die has geometrically complicated features resulting in more complex flow during die filling, and therefore enhanced heat transfer. On the other hand, for the casting/top die interface, the HTC for the solid/solid regime of behavior is the highest as the wheel contracts onto the top die, placing the interface under pressure, and minimizing the development of resistance associated with gap formation. below the solidus and the wheel is fully solidified, the contact degrades, resulting in lower values of HTCs due to the formation of solid/solid interfaces. Depending on the interface location, it is possible for a gap to form between the casting and the die surface, which leads to lower HTC values. In the region between the eutectic transformation point and solidus, the HTCs are assumed to vary linearly. • Values of the HTCs for the different casting/die interfaces are given in Table 3. For the casting/bottom die interface, the HTC for the liquid/solid regime of behavior is higher compared to the other two interfaces. The rationale is that the bottom die has geometrically complicated features resulting in more complex flow during die filling, and therefore enhanced heat transfer.
On the other hand, for the casting/top die interface, the HTC for the solid/solid regime of Figure 5. Schematic profile of the HTCs applied at the casting/die interface.
As previously described, the casting/die interface changes from an internal interface to an external interface once the die is opened and the wheel is ejected. Correspondingly, these boundaries will be in contact with environmental air, and an HTC of 20 W/(m 2 K) is used together with a reference temperature of 50 • C-see Equation (5).
Heat Transport Across the Die/Die Interfaces-The die/die internal interfaces are solid/solid interfaces throughout the process, and therefore, the HTCs of these interfaces are assumed to be invariant with temperature. Table 4 gives the HTCs used for the die/die interfaces. The values are again determined by a combination of trial-and-error and experience. Generally, the HTCs are high, and reflective of close contact following die assembly and closure of the die components. Similar to the casting/die interfaces, when the various die component interfaces lose contact due to die operations (e.g., separations), then 20 W/(m 2 K) is used for the HTC and 50 • C for the reference temperature. Water Cooling-Heat transfer associated with the water cooling in the cooling channels is one of the most important, but most complicated and dynamic, phenomena occurring during the LPDC process. Production-derived data indicates that boiling of the cooling water occurs at both the water-on and -off stages-refer to Part I of this study [17]. To incorporate boiling and its effect on heat transfer, the method described by Moayedinia [11] is used to calculate the cooling channel HTCs. An example of the relationship used in the model is presented in Figure 6. When the temperature is relatively low, no boiling occurs, and heat is extracted by forced convection. The HTC varies in a range between 4200 and 5000 W/(m 2 K), dependent on the flow rate, diameter of the channel and temperature-dependent properties of the water. At interface temperatures above~120 • C, the water in the cooling channels at the interface starts to boil. Note: the channel pressure is >1 atm. These bubbles agitate the boundary layer between the fluid and die, resulting in enhanced heat transfer-i.e., increases in the HTC. Peak heat transfer is predicted to occur at approximately~160 • C. Above~160 • C, interface contact gradually becomes reduced due to the formation and expansion of large vapor bubbles at the interface, and the HTC drops. Above~240 • C, the HTC is assumed to be constant at 5000 W/(m 2 K). The reference temperature for the water-cooling boundary conditions, T ref , is set to 50 • C, which is based on the average temperature change measured in the cooling water. In addition, as indicated by the in-channel temperature and pressure data, boiling routinely occurs after the cooling water is switched off [17]. Note that, in the commercial process, even though the channels are purged with air for 10 s, there is still a certain amount of residual water left in the cooling channels, which then begins to boil, extracting significant amounts of heat from the die. The amount of water in a given channel was found to be dependent on the location of the channel in the die and the arrangement of the supply/drain lines affixed to the die. In the model, to represent the heat extracted by the post-water-off boiling, a constant heat flux of −1.5 × 10 5 W/m 2 is applied to the channels during the water-off periods. It is acknowledged that using a constant heat flux is a simplified approach. However, this approach was necessary to overcome limitations in the version of the software used in the study-i.e., only one temperature-and time-dependent HTC can be applied per cooling channel.
Air Cooling-Forced air cooling was employed at a number of locations in the dies studied in this work. In the model, an HTC of 200 W/(m 2 K) and reference temperature of 50 °C are used at the locations where air cooling was applied.
Heat Transfer to Environment-Heat transfer on the outer surfaces of the die to the surrounding plant environment can occur through convection and radiation. The convective heat transfer is described in Equation (5), while the radiative transport has been described using the Stefan-Boltzmann equation-Equation (6). The parameters related to the environmental heat transfer are given in Table 5.
where (W/m 2 ) is the heat flux by radiation, is the emissivity and (W/(m 2 K 4 )) is the Stefan-Boltzmann constant. * The duration of air cooling under the non-production process condition is significantly longer than the production process condition which decreases the environment temperature, , and thus the cooling efficiency increases. To standardize the model parameters, instead of changing , proportionally greater HTCs are used for the non-production process condition. ** Depending on the reference die components, different are applied at different locations of the surface based on measured data from the production process condition. In addition, as indicated by the in-channel temperature and pressure data, boiling routinely occurs after the cooling water is switched off [17]. Note that, in the commercial process, even though the channels are purged with air for 10 s, there is still a certain amount of residual water left in the cooling channels, which then begins to boil, extracting significant amounts of heat from the die. The amount of water in a given channel was found to be dependent on the location of the channel in the die and the arrangement of the supply/drain lines affixed to the die. In the model, to represent the heat extracted by the post-water-off boiling, a constant heat flux of −1.5 × 10 5 W/m 2 is applied to the channels during the water-off periods. It is acknowledged that using a constant heat flux is a simplified approach. However, this approach was necessary to overcome limitations in the version of the software used in the study-i.e., only one temperature-and time-dependent HTC can be applied per cooling channel.
Air Cooling-Forced air cooling was employed at a number of locations in the dies studied in this work. In the model, an HTC of 200 W/(m 2 K) and reference temperature of 50 • C are used at the locations where air cooling was applied.
Heat Transfer to Environment-Heat transfer on the outer surfaces of the die to the surrounding plant environment can occur through convection and radiation. The convective heat transfer is described in Equation (5), while the radiative transport has been described using the Stefan-Boltzmann equation-Equation (6). The parameters related to the environmental heat transfer are given in Table 5.
where q rad (W/m 2 ) is the heat flux by radiation, ε is the emissivity and σ (W/(m 2 K 4 )) is the Stefan-Boltzmann constant. * The duration of air cooling under the non-production process condition is significantly longer than the production process condition which decreases the environment temperature, T env , and thus the cooling efficiency increases. To standardize the model parameters, instead of changing T env , proportionally greater HTCs are used for the non-production process condition. ** Depending on the reference die components, different T env are applied at different locations of the surface based on measured data from the production process condition.

Fluid Flow Boundary Conditions
Pressure Inlet-In the actual process, the flow rate of liquid into the die is controlled via the pressure in the holding furnace. In the model, the flow is controlled by the pressure applied at the inlet of the die cavity (sprue). The pressure applied in the model is a modified form of the machine-specified pressure curve following Equation (7), since the transfer tube is excluded from the computational domain.
where P model (Pa) is the pressure used in the model and P machine (Pa) is the process-specified pressure variation with time programmed into the die casting machine. ρ A356 (kg/m 3 ) is the density of A356, g (m/s 2 ) is the gravitational acceleration and h tube (m) is the height of the space between the sprue and liquid metal surface in the holding furnace. The term ρ A356 gh tube represents the pressure required to raise the liquid metal in the transfer tube to the sprue.

Time Steps
The model starts with an initial time step of 0.001 s. Subject to solution convergence, the time step is automatically adjusted but limited to the maximum time steps of 0.1 and 1 s for the filling stage (solution of fluid flow is activated) and the solidification stage (thermal only), respectively. Note: a sensitivity analysis was conducted with these parameters to assess their impact on the model results vs. execution time. These parameters were selected because they offer a balance between computational accuracy and efficiency.

Results and Discussion
Part I of this work reported on the method used for the process characterization of Die-A for two process conditions, and presented and discussed a representative selection of the data obtained-see reference [17]. The model validations for the two process conditions for Die-A are presented in Section 3.1, Section 3.2 and Section 3.3, followed by a comparison of the model predictions with a selection of data obtained from the production process for Die-B (Section 3.4).

In-Die Temperature
The location of the nine in-die TCs used for model verification are shown in Figure 7a. Each die component has three TCs, which are distributed widely-see BD61, 62 and 63, SD44, 43 and 42, and TD11, 13 and 17. Figure 7a also shows the location of the cooling elements in each of the die components of Die-A-see B1, 2, 3, 5, S2, and T1, 2, 3, 4, 5, 7. B1, B5, S2, T1, T3, T5 and T7 are all water-cooling channels/elements. The remainder are air-cooling elements. The approach to die design is to have a generic number of elements, some of which are then activated for a given wheel geometry, thus providing flexibility in process control.  Production Process Condition-For the production process condition, Figure 8a-c present the indie temperature comparison between the model and measured data at the TC locations in the top, side and bottom dies, respectively. Note: the cooling timings of the active cooling elements in each die component are also illustrated in the figures by the grey bars. In addition, the averaged error Production Process Condition-For the production process condition, Figure 8a-c present the in-die temperature comparison between the model and measured data at the TC locations in the top, side and bottom dies, respectively. Note: the cooling timings of the active cooling elements in each die component are also illustrated in the figures by the grey bars. In addition, the averaged error between the modeled and measured temperature data has been calculated by Equation (8) and displayed in Figure 9.
where T error (K) is the error in temperature between the model, T i, model (K), and experimental measurements, T i, exp (K), at time index i, and N is a large number (300), representing the number of sub-increments the casting cycle time has been divided into. measurements, , (K), at time index i, and N is a large number (300), representing the number of sub-increments the casting cycle time has been divided into. Overall, it can be seen that the model agrees very well with the measured data at all the locations in each die component with respect to both the magnitude and variation of temperature with time. The error ( ) is in the range of 6 to 27 °C. More specifically, the model is capable of capturing all three stages in the temperature evolution at all of the locations. The three stages and phenomena that give rise to them are discussed in detail in Part I of this work [17]. Briefly: • Stage-1 is characterized by a small drop in temperature until such time as the pulse of heat associated with the presence of liquid metal in the die cavity is able to diffuse through the die to the TC location; • Stage-2 is characterized by a rise in temperature associated with the transfer of heat from the aluminum into the die; • Stage-3 is characterized by a period of slow decline in die temperature associated with the balance between removal of heat from the die and the supply of heat from the solidifying wheel.      A closer look at the comparison shows that relatively large discrepancies appear at two locations, TD17 and BD61, and are both approximately 50 °C in magnitude at their maximum. They also both occur only in the third stage of the temperature evolution. As the third stage, characterized by cooling, is dependent on the activation of water-cooling elements, it is therefore speculated that the boundary conditions in the model describing the water cooling are responsible. Firstly, the effect of boiling has been included; however, its description is still simplified, and secondly, there is a delay between the time when the cooling water valve opens and when the water actually arrives at the in-die channelssee the details in Sections 4.1.2 and 4.1.4 of reference [17]. However, considering that one major aim of this work is to standardize the model parameters for extensive industrial use, and given that the delay varies from channel-to-channel, the cooling timings in the model were not adjusted relative to the programmable logic controller (PLC) set values to improve the fitting.
Non-Production Process Condition-In the non-production process condition, the cooling parameters were adjusted while the rest of the process parameters and boundary conditions were kept identical to the production process condition for Die-A. Figure 10a-c present the in-die temperature comparison at the selected TC locations in the top, side and bottom dies, respectively. Additionally, the error in temperature between model and measured data is presented in Figure 11. In general, at most locations (except for TD17), the model results match closely to the measured data Overall, it can be seen that the model agrees very well with the measured data at all the locations in each die component with respect to both the magnitude and variation of temperature with time. The error (T error ) is in the range of 6 to 27 • C. More specifically, the model is capable of capturing all three stages in the temperature evolution at all of the locations. The three stages and phenomena that give rise to them are discussed in detail in Part I of this work [17]. Briefly: • Stage-1 is characterized by a small drop in temperature until such time as the pulse of heat associated with the presence of liquid metal in the die cavity is able to diffuse through the die to the TC location; • Stage-2 is characterized by a rise in temperature associated with the transfer of heat from the aluminum into the die; • Stage-3 is characterized by a period of slow decline in die temperature associated with the balance between removal of heat from the die and the supply of heat from the solidifying wheel.
A closer look at the comparison shows that relatively large discrepancies appear at two locations, TD17 and BD61, and are both approximately 50 • C in magnitude at their maximum. They also both occur only in the third stage of the temperature evolution. As the third stage, characterized by cooling, is dependent on the activation of water-cooling elements, it is therefore speculated that the boundary conditions in the model describing the water cooling are responsible. Firstly, the effect of boiling has been included; however, its description is still simplified, and secondly, there is a delay between the time when the cooling water valve opens and when the water actually arrives at the in-die channels-see the details in Sections 4.1.2 and 4.1.4 of reference [17]. However, considering that one major aim of this work is to standardize the model parameters for extensive industrial use, and given that the delay varies from channel-to-channel, the cooling timings in the model were not adjusted relative to the programmable logic controller (PLC) set values to improve the fitting.
Non-Production Process Condition-In the non-production process condition, the cooling parameters were adjusted while the rest of the process parameters and boundary conditions were kept identical to the production process condition for Die-A. Figure 10a-c present the in-die temperature comparison at the selected TC locations in the top, side and bottom dies, respectively. Additionally, the error in temperature between model and measured data is presented in Figure 11. In general, at most locations (except for TD17), the model results match closely to the measured data (T error < 23 • C). This demonstrates that the modeling approach is able to accurately capture a broad range of process cooling timings.   Temperature comparison between the model and measured data for the non-production process condition in Die-A. Figure 11. Error between the model and measured data. Figure 11. Error between the model and measured data.
To assess the extent to which the changes in the process conditions (production vs. non-production) have changed the die temperature field at the steady state, and also the ability of the model to capture this change, the average change in temperature ∆T avg has been calculated for the top, bottom and side TC locations according to the expression given in Equation (9).
where ∆T i (K) is the difference in temperature between the production and non-production temperatures at time index i and N is a large number (300), representing the number of sub-increments the casting cycle time has been divided into. This has been done for both the TC data and for the model output at the TC locations. The results are presented in Figure 12 and indicate that the temperature difference between the two process conditions examined ranges from 30 to 90 • C, for most of the TC locations. At these locations, the differences based on the model predictions are in relatively good agreement. At TD17, however, a much larger average temperature difference in the process data is observed (~185 • C), whereas the model predicts this difference to be~120 • C-a mismatch of approximately 65 • C. This mismatch is large and warrants a closer look.

In-Wheel Temperature
The temperature evolution with time within the wheel at multiple locations was also examined in this work, which provides an opportunity for a direct validation of the model with respect to the thermal phenomena occurring during the solidification process. Figure 14a,b present a comparison between the in-wheel temperature measurements and the model predictions for the production and non-production process conditions, respectively. The in-wheel TC locations are shown in Figure 15. Those locations labeled with a P are for the production process condition and those labeled with an TD17 was located very close to the water-cooling channel T3, which was instrumented with a TC (TD21)-refer to reference [17] Section 4.1.4 for the details related to the instrumentation of water-cooling channels. The results for TD21 for the production and non-production process conditions are presented in Figure 13a,b, respectively. Note: the water-on and air-purging periods are marked by the light purple and green bars, respectively, in the two figures. The major difference between these two sets of data is that, in the production process condition, the residual water remaining after air purging is fully vaporized and the channel dries out (after~280 s). The "dry out" is indicated by the rise in temperature to well in excess of 100 • C, which continues through into the next casting cycle until the water-cooling is reactivated, whereas in the non-production process, the residual water in channel T3 is not fully vaporized (the temperature within the cooling channel remains just below 100 • C). Thus, the extended cooling time applied to channel T3 has had a secondary effect, in that by cooling the die to a larger extent, residual water is retained in the channel. Consequently, the period over which boiling occurs during the water-off period for the channel is extended from~40 to~90 s, leading to the large ∆T avg observed in Figure 12. This effect is not currently captured within the model, hence the model under-predicts the effect of extending the cooling timing in channel T3. It also points to the need to fully understand the cooling/boiling behavior in the water channel during and after air purging.

In-Wheel Temperature
The temperature evolution with time within the wheel at multiple locations was also examined in this work, which provides an opportunity for a direct validation of the model with respect to the thermal phenomena occurring during the solidification process. Figure 14a,b present a comparison between the in-wheel temperature measurements and the model predictions for the production and non-production process conditions, respectively. The in-wheel TC locations are shown in Figure 15. Those locations labeled with a P are for the production process condition and those labeled with an NP are for the non-production process condition. Note that there was only one successful insertion of a TC into the die cavity for the non-production process condition. It should be mentioned that these are challenging measurements to make, and that the TCs have to be inserted by hand into the die cavity and are prone to be cut when the die closes. There are two periods of time in which the data from the TCs should be disregarded: the first is at the beginning of the cycle prior to the liquid contacting the TC, and the second is at the end of the cycle when the TC is cut prior to the ejection of the wheel from the top die. In the figures, these periods are highlighted in light gray for clarity.

In-Wheel Temperature
The temperature evolution with time within the wheel at multiple locations was also examined in this work, which provides an opportunity for a direct validation of the model with respect to the thermal phenomena occurring during the solidification process. Figure 14a,b present a comparison between the in-wheel temperature measurements and the model predictions for the production and non-production process conditions, respectively. The in-wheel TC locations are shown in Figure 15. Those locations labeled with a P are for the production process condition and those labeled with an NP are for the non-production process condition. Note that there was only one successful insertion of a TC into the die cavity for the non-production process condition. It should be mentioned that these are challenging measurements to make, and that the TCs have to be inserted by hand into the die cavity and are prone to be cut when the die closes. There are two periods of time in which the data from the TCs should be disregarded: the first is at the beginning of the cycle prior to the liquid contacting the TC, and the second is at the end of the cycle when the TC is cut prior to the ejection of the wheel from the top die. In the figures, these periods are highlighted in light gray for clarity.    From the comparison, it can be seen that, overall, the model-predicted in-wheel temperatures follow the profiles of the measured data relatively well. The poorest agreement is observed for locations P_cycle21B and P_cycle19B, both of which are located in the spoke near the rim-spoke junction. In P_cycle21B, the error is approximately 50 °C at the end of the cycle. As can be seen, the model and experimental data both reflect the latent heat release associated with primary α-aluminum solidification at approximately ~600 °C and eutectic formation at ~560 °C, as evidenced by periods of reduced rates in temperature decline-see reference for solidification sequence [19].

Solidification Direction
The solidification path predicted by the model has been used to assess the locations that would be prone to form porosities. The solidification front should ideally follow a path that continuously provides access to liquid to offset the density increase associated with the liquid-to-solid transformation. The failure to achieve this condition may lead to the formation of regions containing From the comparison, it can be seen that, overall, the model-predicted in-wheel temperatures follow the profiles of the measured data relatively well. The poorest agreement is observed for locations P_cycle21B and P_cycle19B, both of which are located in the spoke near the rim-spoke junction. In P_cycle21B, the error is approximately 50 • C at the end of the cycle. As can be seen, the model and experimental data both reflect the latent heat release associated with primary α-aluminum solidification at approximately~600 • C and eutectic formation at~560 • C, as evidenced by periods of reduced rates in temperature decline-see reference for solidification sequence [19].

Solidification Direction
The solidification path predicted by the model has been used to assess the locations that would be prone to form porosities. The solidification front should ideally follow a path that continuously provides access to liquid to offset the density increase associated with the liquid-to-solid transformation. The failure to achieve this condition may lead to the formation of regions containing partially solidified material that is "shut-off" from a supply of liquid and encapsulated in solid. Depending on the structural integrity of the surrounding solid, the cleanliness of the liquid metal (e.g., oxide films) and the hydrogen content in solution in the liquid metal, this could result in the formation of distributed porosity (e.g., shrinkage-or hydrogen-based) or surface depressions/deformation.
Focusing first on the production process condition, Figure 16 presents a series of images (a-f) of the wheel within the die at increasing times within the cast cycle. In each image, the portion of the wheel with a solid fraction, f s , greater than 0.7 is removed or "cut-off" in the image, thereby allowing the progress of the solidification front within the die to be observed. As can be seen in the series of images, in general the solidification front follows a path from the upper die (in-board rim flange) to the sprue. Small volumes of liquid encapsulated (~18 mm in equivalent radius) can be seen at the spoke/hub junction in image (f), at PTS (percentage of total solidification) = 93%. A close-up view of the bottom die structure in the vicinity of the area of liquid encapsulation is shown in Figure 17. As can be seen, the structure of the die is complex, resulting in a large surface area for heat transfer per volume of material undergoing solidification, likely leading to the formation of the region of encapsulation. As presented in Section 4.1.5 in reference [17], no noticeable shrinkage-based porosity was observed in the area in question from the line-of-sight X-ray images taken of example wheels cast under the production process condition in a cyclic steady state. This indicates that either the model falsely predicted the area of liquid encapsulation, or alternatively the size of the porosity was below the threshold for detection in the X-ray imaging system used.
can be seen, the structure of the die is complex, resulting in a large surface area for heat transfer per volume of material undergoing solidification, likely leading to the formation of the region of encapsulation. As presented in Section 4.1.5 in reference [17], no noticeable shrinkage-based porosity was observed in the area in question from the line-of-sight X-ray images taken of example wheels cast under the production process condition in a cyclic steady state. This indicates that either the model falsely predicted the area of liquid encapsulation, or alternatively the size of the porosity was below the threshold for detection in the X-ray imaging system used.     Turning to the non-production process condition, only the stages with the formation of liquid encapsulation are presented in Figure 18 for brevity. It can be seen that liquid encapsulation formed in two regions: the rim/spoke junction at PTS = 85%, and the spoke/hub junction at PTS = 94%. The former is ~27 mm in equivalent radius, whereas the latter has a similar location and size (~16 mm in equivalent radius) to that predicted in the production process condition wheel. Compared to the Xray radiography results [17], Figure 19, with the larger area of liquid encapsulation predicted by the model, indicates consistency. This lends some support to the conclusion that the other smaller areas of liquid encapsulation yielded porosities below the detection limit. Turning to the non-production process condition, only the stages with the formation of liquid encapsulation are presented in Figure 18 for brevity. It can be seen that liquid encapsulation formed in two regions: the rim/spoke junction at PTS = 85%, and the spoke/hub junction at PTS = 94%. The former is~27 mm in equivalent radius, whereas the latter has a similar location and size (~16 mm in equivalent radius) to that predicted in the production process condition wheel. Compared to the X-ray radiography results [17], Figure 19, with the larger area of liquid encapsulation predicted by the model, indicates consistency. This lends some support to the conclusion that the other smaller areas of liquid encapsulation yielded porosities below the detection limit.

Further Validation on Die-B
In order to further evaluate the robustness of this modeling methodology, a second die structure (Die-B) was simulated. Temperature data have also been acquired from the die in a production setting and compared to the model predictions. Figure 20 displays the geometry of the second wheel. This is a conventional LPDC 17 × 7 in, 10-spoke wheel, and the die for producing this wheel was cooled by water only. The standard production parameters were used for the plant trial. Note: in comparison, the wheel produced by Die-A is a conventional 17 × 7.5 in, 6-spoke wheel, and the die was cooled by both water and air.

Further Validation on Die-B
In order to further evaluate the robustness of this modeling methodology, a second die structure (Die-B) was simulated. Temperature data have also been acquired from the die in a production setting and compared to the model predictions. Figure 20 displays the geometry of the second wheel. This is a conventional LPDC 17 × 7 in, 10-spoke wheel, and the die for producing this wheel was cooled by water only. The standard production parameters were used for the plant trial. Note: in comparison, the wheel produced by Die-A is a conventional 17 × 7.5 in, 6-spoke wheel, and the die was cooled by both water and air.   Figure 21 presents the comparison between the measured data and model predictions for the in-die temperatures (the error, T error , in Figure 22), and the comparison of in-wheel temperatures in Figure 23. The locations of the TCs are shown in Figure 24. As can be seen from the comparison, the model is able to reproduce the behavior observed at the various locations within the die and wheel selected for the comparison accurately. By comparing the temperature error of this case (Figure 22) to the other two process conditions (Figures 9 and 11), it can be seen that the accuracy does not degrade (the maximum T error is only~23 • C at BD55). It is worthwhile to mention again that the model parameters were not re-tuned for the analysis of Die-B. This result provides evidence to support the robustness of the modeling methodology, insofar as the model is able to capture the key phenomena governing the transport of heat independently of die geometry and the cooling conditions within the range encompassed by the two wheel/die combinations examined. The result also bodes well for the model being able to accurately predict a broader range of die/wheel geometries.  (c) TCs in bottom die Figure 21. In-die temperature comparison between the model and measured data for the second die/wheel geometry. Figure 22. Error between the model and measured data. Figure 22. Error between the model and measured data.

Summary
An advanced modeling methodology has been developed using the commercial software package ProCAST for the industrial LPDC-based production of aluminum alloy A356 automotive wheels. The methodology incorporates the relevant material properties, the boundary conditions to

Summary
An advanced modeling methodology has been developed using the commercial software package ProCAST for the industrial LPDC-based production of aluminum alloy A356 automotive wheels. The methodology incorporates the relevant material properties, the boundary conditions to characterize interfacial heat transfer associated with the various interfaces in the process (including

Summary
An advanced modeling methodology has been developed using the commercial software package ProCAST for the industrial LPDC-based production of aluminum alloy A356 automotive wheels. The methodology incorporates the relevant material properties, the boundary conditions to characterize interfacial heat transfer associated with the various interfaces in the process (including die/wheel, die/die, die/cooling water and die/cooling air), and the latent heat associated with solidification. The model has been comprehensively validated using in-die and in-wheel temperature data collected from a commercial production casting machine operating at steady state. Aspects of the methodology developed to describe the die/cooling water interface boundary condition benefited significantly from in-channel water temperature data, which indicated the presence of periods of boiling water heat transfer, both when water is first introduced to the channel and after the application of purging air.
The validation was carried out over a range of process conditions, including (1) production process conditions for Die-A, (2) non-production process conditions for Die-A, and (3) production process conditions for a second die structure, Die-B. The model was originally developed and some interface conditions were "tuned" based on analysis of the production conditions for Die-A. The model was then applied without re-tuning to a second set of process conditions (non-production), in which the cooling channel timing was significantly altered. Finally, the model was applied to a second die structure, Die-B, also without retuning. In addition to changes in geometry, the process conditions for Die-B were different in two ways: (1) the cooling timing was different, and (2) the die contained only water-cooling elements, unlike Die-A which utilized both air-and water-cooling elements.
The comparisons presented between the in-process-derived TC data and the model predictions indicated that the model is able to reproduce the temperature evolution at the locations examined in the wheel and die accurately for the range of process conditions/die structures examined in the study, without the need for re-tuning. Comparisons to X-ray images taken of wheels produced from Die-A for the non-production process condition indicate that the model is able to predict the location of shrinkage-based porosity associated with pockets of encapsulated liquid of~2 m equivalent radius. Other areas of encapsulated liquid predicted by the model, for both the production and non-production process conditions, that were smaller in size (~16 to 18 mm in equivalent radius) did not appear to have shrinkage porosity associated with them, according to the X-ray images. It was suggested that this may be due to limits in the resolution of the line-of-sight X-ray images used to detect smaller clusters of shrinkage porosity. Additional work is required to investigate this further.