# New Die Spray Cooling Model

In this last blog on FLOW-3D CAST v4.2’s developments, I discuss the new Die Spray Cooling Model for thermal die cycling simulations.

In order to accurately predict the temperature distribution in the die in thermal die cycling simulations, the spatial variation in the spray cooling needs to be modeled. The Die Spray Cooling model was developed for this purpose, and is now available in the recent release of FLOW-3D CAST v4.2. This model explicitly computes the cooling from each spray, instead of assuming a constant heat transfer coefficient across the entire die cavity. The spray area on the die surface is computed and updated constantly due to the movement of the spray nozzles. The model also takes into account the blocking of the fluid being sprayed, and thus the effect on cooling, due to the spray angle and the shape of the die surface. The new model provides accurate temperature distribution on the die surfaces with reliable and realistic input parameters, which helps engineers to better design and optimize the cooling process to eliminate hot spots. A sample thermal die cycling simulation demonstrates the capability of the new model, and the importance of explicitly simulating the die spraying process.

Spray area during spray cooling. Simulation courtesy of Audi AG.

Die surface temperature at the fifth cycle of spray cooling. Simulation courtesy of Audi AG.

This simulation is based on the production of a vehicle structural part with large dimension and thin wall thickness. For this example, the simulation is limited only on the ejector part of the die. The two animations show the spray area on the die surface during spray cooling and the die surface temperature at the fifth cycle of spray cooling. The effect of global spraying and hot spot spraying can be clearly seen and identified.

## Spray Area Computation

The influence of the shape of the die surface and the position of the spraying nozzles on the spray cooling is taken into account in the new model. As shown in the figure below, some areas under the spray on the die surface are blocked, and some areas are overlapped by two or more sprays. These areas are computed and identified using a ray tracing algorithm to distinguish the different spray cooling effects. The spray area can be visualized in FlowSight, along with other properties, such as the total spray time and total heat removed per unit area by the spray cooling.

## Determining Heat Transfer Coefficients

Spray cooling mechanisms are complicated, and the spray cooling heat transfer coefficient (HTC) depends on many variables such as spray shape, coolant flow rate, spray pressure, mold temperature, spray angles, and spray distance. To simplify the spray cooling HTC calculation, for every spray surface element, the HTC is computed by multiplying the base HTC by dependence factors, for example for a conical spray,

$\displaystyle HTC=HT{{C}_{0}}(T)\cdot {{f}_{d}}(d)\cdot {{f}_{b}}(\beta )\cdot {{f}_{e}}(\varepsilon )$, where

• HTC0 is the base spray heat transfer coefficient when the nozzle sprays on the mold from a specified distance. The base heat transfer coefficient is dependent on the spray cone characteristics, spray medium, and spray pressure, etc., and is a function of the mold surface temperature.
• fd is the distance d dependent factor function.
• fb is the spray angle β dependent factor function.
• fe is the spray angle (between the surface normal and the spray direction) ε dependent factor function

The meaning of the spray distance d, and spray angles β and ε are shown in the picture below.

The base heat transfer coefficient and the dependence factor functions can be curve fitted from experimental measurements, derived from theory, or from experience. If the sprays are not conical in shape, the dependence factors could be different.