COMSOL 6.4 - Cooling of an Injection Mold

Cooling is an important process in the production of injection molded plastics. First of all, the cooling time may well represent more than half of the production cycle time. Second, a homogeneous cooling process is desired to avoid defects in the manufactured parts. If plastic materials in the injection molding die are cooled down uniformly and slowly, residual stresses can be avoided, and thereby the risk of warps and cracks in the end product can be minimized.

As a consequence, the positioning and properties of the cooling channels become important aspects when designing the mold.

The simulation of heat transfer in molds of relatively complex geometries requires a 3D representation. Simulation of 3D flow and heat transfer inside the cooling channels are computationally expensive. An efficient short-cut alternative is to model the flow and heat transfer in the cooling channels with 1D pipe flow equations, and still model the surrounding mold and product in 3D.

This example shows how you can use the interface together with the interface to model a mold cooling process. The equations describing the cooling channels are fully coupled to the heat transfer equations of the mold and the polyurethane part using multiphysics coupling.

Figure 1: The steering wheel of a car, made from polyurethane.

Model Definition

Model Geometry and Process conditions

The polyurethane material used for a steering wheel is produced by several different molds. The part considered in this model is the top half of the wheel grip, shown in gray in Figure 2.

Figure 2: Polyurethane parts for a steering wheel. The top half of the grip is modeled in this example.

The mold consists of a 50-by-50-by-15 cm steel block. Two cooling channels, 1 cm in diameter, are machined into the block as illustrated in Figure 3.

Figure 3: Mold block and cooling channels.

After the injection of the polyurethane, the average temperature of the mold and the plastic material is 473 K. Water at room temperature is used as cooling fluid and flows through the channels at a rate of 10 liters/min. The model simulates a 10 min cooling process.

For numerical stability reasons, the model is set up with an initial water temperature of 473 K, which is ramped down to 288 K during the first few seconds.

Results and Discussion

The mold and polyurethane part, initially at 473 K, are cooled for 10 minutes by water at room temperature. Figures below show sample results when flow rate of the cooling water is 10 liters/minute and the surface roughness of the channels is 46 μm. After two minutes of cooling, the hottest and coldest parts of the polyurethane part differ by approximately 50 K (Figure 4).

Figure 4: Temperature distribution in the polyurethane part and the cooling channels after 2 minutes of cooling.

Figure 5 shows the temperature distribution in the steel mold after 2 minutes. The temperature footprint of the cooling channels is clearly visible.

Figure 5: Temperature distribution in the steel mold block after 2 minutes of cooling.

After 10 minutes of cooling, the temperature in the mold block is more uniform, with a temperature at the center of approximately 333 K (Figure 6). Still, the faces with cooling channel inlets and outlets are more than 20 K hotter.

Figure 6: Temperature distribution in the steel mold block after 10 minutes of cooling.

To evaluate the influence of factors affecting the cooling time, use Material and Parametric sweeps. Figure 7 shows the average temperature of the polyurethane part as function of the cooling time for the several flow rates of the cooling water, the surface roughness of the cooling channels, and the mold materials.

Figure 7: Average temperature of the polyurethane part as function of time and cooling conditions.

Pipe_Flow_Module/Heat_Transfer/mold_cooling

Modeling Instructions

From the menu, choose .

New

In the window, click

Model Wizard

In the window, click

In the tree, select > > .

Click .

In the tree, select > .

Click .

Click

In the tree, select > .

Click

Global Definitions

Parameters 1

In the window, under click .

In the window for , locate the section.

In the table, enter the following settings:


Name	Expression	Value	Description
T_init_mold	473.15[K]	473.15 K	Initial temperature, mold
T_coolant	288.15[K]	288.15 K	Steady-state inlet temperature, coolant
Qw	10[l/min]	1.6667E-4 m³/s	Coolant flow rate
e	0.046[mm]	4.6E-5 m	Surface roughness

Step 1 (step1)

Create a smooth step function to decrease the coolant temperature at the beginning of the process.

In the toolbar, click

and choose > .

In the window for , locate the section.

In the text field, type 2.5.

In the text field, type 1.

In the text field, type 0.

Click to expand the section. In the text field, type 5.

Optionally, you can inspect the shape of the step function:

Click

Variables 1

In the toolbar, click

and choose .

In the window for , locate the section.

In the table, enter the following settings:


Name	Expression	Unit	Description
T_inlet	T_coolant+(T_init_mold-T_coolant)*step1(t[1/s])	K	Ramped inlet temperature, coolant

Geometry 1

Import 1 (imp1)

First, import the steering wheel part from a CAD design file.

In the toolbar, click

In the window for , locate the section.

Click

Browse to the model’s Application Libraries folder and double-click the file mold_cooling_top.mphbin.

Click

Click to expand the section. Find the subsection. Click .

In the dialog, type Wheel in the text field.

Click .

Second, draw the mold and cooling channels. To simplify this step, insert a prepared geometry sequence from file. After insertion you can study each geometry step in the sequence.

Click the

button in the toolbar.

In the toolbar, click and choose .

Browse to the model’s Application Libraries folder and double-click the file mold_cooling_geom_sequence.mph.

In the toolbar, click

Click the

button in the toolbar.

Click the

button in the toolbar.

Work Plane 1 (wp1)

Create the selections to simplify the model specification.

In the window, click .

In the window for , locate the section.

Find the subsection. Click .

In the dialog, type Cooling channels in the text field.

Click .

Disable the analysis of the geometry as the remaining small geometric details can be kept.

In the window, click .

In the window for , locate the section.

Clear the checkbox.

Materials

The next step is to specify material properties for the model. First, select water from the built-in materials database.

Add Material

In the toolbar, click

to open the window.

Go to the window.

In the tree, select > .

Click the button in the window toolbar.

Materials

Water, liquid (mat1)

In the window for , locate the section.

From the list, choose .

Material Switch 1 (sw1)

Define the mold materials to switch between during a solver sweep.

In the toolbar, click

and choose > .

Select Domain 1 only.

Add Material

Go to the window.

In the tree, select > .

Right-click and choose .

In the tree, select > .

Right-click and choose .

In the toolbar, click

to close the window.

Materials

Polyurethane

Next, create a material with the properties of polyurethane.

In the window, under right-click and choose .

In the window for , type Polyurethane in the text field.

Locate the section. From the list, choose .

Locate the section. In the table, enter the following settings:


Property	Variable	Value	Unit	Property group
Thermal conductivity	k_iso ; kii = k_iso, kij = 0	0.32	W/(m·K)	Basic
Density	rho	1250	kg/m³	Basic
Heat capacity at constant pressure	Cp	1540	J/(kg·K)	Basic

Nonisothermal Pipe Flow (nipfl)

In the window, under click .

In the window for , locate the section.

From the list, choose .

Pipe Properties 1

In the window, under > click .

In the window for , locate the section.

From the list, choose .

In the di text field, type 1[cm].

Locate the section. From the list, choose .

In the text field, type e.

Initial Values 1

In the window, click .

In the window for , locate the section.

In the u text field, type 0.1.

In the T text field, type T_init_mold.

Temperature 1

In the window, click .

In the window for , locate the section.

In the Tin text field, type T_inlet.

Inlet 1

In the toolbar, click

and choose .

Select Points 3 and 4 only.

In the window for , locate the section.

From the list, choose .

In the qv,0 text field, type Qw.

Heat Outflow 1

In the toolbar, click

and choose .

Select Points 269 and 270 only.

Wall Heat Transfer 1

In the toolbar, click

and choose .

In the window for , locate the section.

From the list, choose .

Internal Film Resistance 1

In the toolbar, click

and choose .

Heat Transfer in Solids (ht)

Initial Values 1

In the window, under > click .

In the window for , locate the section.

In the T2 text field, type T_init_mold.

Heat Flux 1

In the toolbar, click

and choose .

In the window for , locate the section.

From the list, choose .

Locate the section. From the list, choose .

In the h text field, type 2.

Multiphysics

Pipe Wall Heat Transfer 1 (pwhtc1)

In the toolbar, click

and choose > .

Mesh 1

Edge 1

In the toolbar, click

and choose .

In the window for , locate the section.

From the list, choose .

Size 1

Right-click and choose .

In the window for , locate the section.

From the list, choose .

Free Tetrahedral 1

In the toolbar, click

Even with a maximum element size of 0.003 m, the mesh contains some collapsed elements, resulting in solver errors. Trial and error gives that lowering the curvature factor somewhat will create a mesh with good quality.

Size

In the window, click .

In the window for , locate the section.

Click the button.

Locate the section. In the text field, type 0.003.

In the text field, type 0.55.

In the window, right-click and choose .

Study 1

Step 1: Time Dependent

In the window, under click .

In the window for , locate the section.

In the text field, type range(0,2,28) range(30,30,600).

Definitions

To evaluate the average temperature of the polyurethane part for different conditions and mold materials (Figure 7), perform parametric and material sweeps. To avoid accumulating a lot of data while solving, keep only last solution and save the average wheel temperature in a table. For this purpose, add a global probe to the model.

Domain Probe 1 (dom1)

In the toolbar, click