Injection Molding

1. Create Case

This tutorial will require high CPU usage and might take several hours to finish the calculation depending on the machine used.

After opening SimFlow, we will create a new case named housing_injection

1. Go to New panel
2. Provide name housing_injection
3. Click Create Case

2. Import Geometry

Firstly we need to Download Geometry: housing.stl

1. Click Import Geometry
2. Select geometry file housing.stl
3. Click Open

3. Geometry - Housing

After importing geometry, it will appear in the 3D view.

1. Click Fit View to fit geometry in the view

4. Scale Geometry - Housing (I)

The geometry is defined in millimeters unit. We need to scale it to meters.

1. Expand the Options list for housing geometry
2. Select Scale

5. Scale Geometry - Housing (II)

1. Check Uniform option
2. Type Scaling Factors as 1e-03 and press Enter
3. Click Scale

6. Split Geometry - Housing (I)

The imported geometry is made up of a single surface. We need to split it into separate faces for further processing.

1. Expand the Options list
2. Select Split

7. Split Geometry - Housing (II)

1. Click Split

8. Create Face Groups (I)

We will add 3 face groups to the original geometry. SimFlow treats each face group as a mesh boundary. By using face groups, we can create separate boundaries for inlet and outlets.

1. Select the inlet face by holding CTRL button and clicking on the geometry face
2. Press Create face group from selection button

9. Create Face Groups (II)

New component should be appeared on the Face Groups list.

1. Change face group name from group_1 to inlet
(double click to edit name and press Enter to confirm)
2. To uncheck the selected face from graphics click on Clear Selection or press the Esc key

10. Create Face Groups (III)

Repeat this operation for the outlets. Additionally, we need to add outlet tube walls as separate groups to set the mesh refinement for these regions. Each group should contain a single face.
1. Create and rename new groups as follows:

11. Create Face Groups (IV)

In order to be able to set different settings per each face, we will convert the single geometry into multiple geometry objects.

1. Select Geometry Faces
2. Click on the inlet group
3. Click on Extract group
4. Repeat this operation for the rest of created face groups

12. Meshing Parameters - Housing

Before we start meshing we will specify geometry options used during the meshing. We will also increase the refinement level for better accuracy.

1. Go to Hex Meshing panel
2. Select the housing component
3. Enable Mesh Geometry
4. Set Refinement to Min 2 Max 2

13. Copy Meshing Parameters

We will use the copy option to apply common settings.

1. Select housing
2. Click on Copy Meshing Properties
3. Select all remaining faces as a target
4. Click on Copy

14. Meshing Parameters - Housing Outlet 1 Tube Wall

For outlet tubes, we will increase copied refinement levels.

1. Select the housing_outlet1_tube_wall
2. Increase the Refinement to Min 3 Max 3

15. Meshing Parameters - Housing Outlet 2 Tube Wall

1. Select the housing_outlet2_tube_wall
2. Increase the Refinement to Min 3 Max 3

16. Base Mesh - Domain

Under the Base tab we will define the background mesh for the housing.

1. Go to Base tab
2. Click on Autosize
3. Set the division of the base mesh
Division 41 20 20
4. Change Max [m] Z to 0.0205

17. Material Point

Material Point tells the meshing algorithm on which side of the geometry the mesh is to be retained. Since we are considering flow inside the housing we need to place the material point inside the geometry.

1. Go to Point tab
2. Specify location inside housing geometry
Material Point 0.019 1e-02 1e-02

You can specify the point location from the 3D view. Hold the Ctrl key and drag the arrows to the destination.

18. Start Meshing

1. Switch to Mesh tab
2. Start the meshing process with Mesh button

19. Check Mesh - Section View (I)

After the meshing process is finished the mesh should appear in the graphics window. In order to check the quality of the mesh, we view the element’s shape using section view.

1. Select Section View from panel
2. Check the Section View
3. Set the normal vector to the section plane
Normal [-] -1 0 0
4. Set the point to set the location of section plane
Point [m] 2e-03 0 0

20. Check Mesh - Section View (II)

After the mesh check, we will turn off the section view to return to the full model view.

1. Uncheck the Section View
2. Click Section View options or press ESC key to close panel

21. Mesh - Boundaries

To be able to assign inlet and outlet boundary conditions we have to adjust the boundary types.

1. Set the boundary type to Patch for inlet and outlets boundaries
inlet Patch
outlet1 Patch
outlet2 Patch

22. Select Solver

For injection molding, we will use Inter (interFoam) solver. This solver is designed to model two-phase flow with interface capturing capabilities.

1. Go to SETUP panel
2. Select Transient filter
3. Check Multiphase model filter
4. Pick Inter (interFoam) solver
5. Select solver

23. Transport Properties - HDPE

In the transport properties panel, we will define fluid parameters. Firstly, we will define High-Density Polyethylene (HDPE) phase using Power Low transport model.

1. Go to Transport Properties panel
2. Set the surface tension between fluids
σ[N/m] 0.026
3. Change phase name from phase1 to HDPE
4. Set the density
ρ[kg/m3] 952
5. Extend the list for Transport Model and choose Power Law
6. Set the parameters accordingly
νmin[m2/s 6e-03
νmax[m2/s 1000
k[m3/s 21.0084
n[-] 0.41

24. Transport Properties - Air

For the second phase, we will use predefined air properties from the material database.

1. Change phase name from phase2 to air
2. Click on Material Database
3. Select the air
4. Click Apply

25. Solution - PIMPLE

For the purpose of the simulation, we will change the default PIMPLE algorithm parameters.

1. Go to Solution panel
2. Switch to the PIMPLE tab
3. Set Non-Orthogonal Correctors to 0
4. Uncheck Momentum Predictor for stability

26. Solution - Limits

During computations sometimes velocity can spontaneously approach unphysical values due to numerical errors. In order to prevent computations from divergence, we can limit the velocity. As the result, we should get more stable computations.

1. Switch to the Limits tab
2. Check Velocity Damping
3. Set the maximum velocity Umax to 10

27. Boundary Conditions - Inlet (Flow I)

In order to define injection, we will use pressure inlet boundary. We will define the pressure by time history table.

1. Go to Boundary Conditions panel
2. Select housing_inlet
3. Change the pressure type to Uniform Fixed Value
4. Set the Value [Pa] as table
5. Press Edit Table

28. Boundary Conditions - Inlet (Flow II)

Fill in inlet pressure values history.

1. Press Add Row two times
2. Fill the table by values
0 1600000
0.5 30000000
1 50000000
1.5 30000000
3. Close the window

29. Boundary Conditions - Inlet (Phases)

To assign injected material, the phase fraction parameter phase is used. HDPE equal 1 will stand for HDPE phase.

1. Select housing_inlet
2. Switch to the Phases tab
3. Set the Type to Fixed Value
4. Set the Value [-] to 1

30. Boundary Conditions - Outlet (Flow I)

Outlets’ boundary conditions simulate outflow to the atmosphere. The default phase setup allows both fluids to exit freely while in case of backflow only air can be sucked into the housing.

1. Select housing_outlet1
2. Switch the tab to Flow
3. Set the pressure type and its value accordingly
p-ρgh Type Total Pressure
p-ρgh p0[Pa] 100000
4. Change the velocity Type to Zero Gradient

31. Boundary Conditions - Outlet (Flow II)

We will copy the boundary condition settings to the second outlet.

1. Select housing_outlet1
2. Press Copy Boundary Conditions
3. Check the housing_outlet2
4. Press Copy

32. Create Geometry - Patch

Before we will set the initial condition we need to define the region representing the initial location of HDPE. For this purpose, we will create a box.

1. Go to GEOMETRYpanel
2. Select Create Box
3. Change geometry name from box_1 to patch
4. Set the origin and box dimensions
Origin[m] 0.012 5e-03 6e-03
Dimensions [m] 1e-02 1e-02 0.02

33. Initial Conditions - Patch

By using the previously created geometry we will tell the solver where should be localized HDPE at the initial state.

1. Go to Initial Conditions panel
2. Switch to Patch tab
3. Check the patch geometry
4. Expand the Fields
5. Check αHDPE
6. Override αHDPE with the value of 1

34. Monitors – Sampling

During calculation, we can observe intermediate results on a section plane or probe plots. Note that runtime post-processing can only be defined before starting calculations and can not be changed later on. Firstly, we need to specify variables that will be sampled. After that, we will define the section plane and probes.

1. Go to Monitors panel
2. Switch to Sampling tab
3. ExpandFields list
4. Select the velocity U and the HDPE phaseαHDPE

35. Monitors – Create Slice (I)

To add sampling data on a plane we need to create a slice. To define a plane we need to set the normal vector and then select the point. In the current step, we will choose the normal vector and in the next step, we will choose a point that belongs to the plane.

1. Select Create Slice
2. Set the normal along Y axis
Normal [-] 0 1 0

36. Monitors – Create Slice (II)

Now we will place the slice to the destination by setting the point coordinates. To achieve this we will use arrows in the graphic window. Before that, we will hide unnecessary geometries.

1. Extend Graphic Objecyt List
2. Uncheck Geometries
3. Set YZ View orientation
4. Click Fit View
5. Holding CTRL button, click on the red arrow and drag it to the middle of the model
6. The slice point coordinates should be like this
Point 0 0.01 0

37. Monitors – Create Probe

To track sampled parameters in a specific point in space, we can add the probe location monitor.

1. Add New Probe
2. Press the viewfinder icon next to Location [m] coordinates
3. Set the probe location on the 3D view. Holding CTRL button, move the probe to the specific location
4. The probe point location should be like this
Location [m] 2e-03 0.017 1e-02

38. Run - Time Control

For any simulation, it is very convenient to let the solver automatically determine the proper time step value. To use this option we need to define time step constraints by providing the initial time step (adjusted by the solver during computations), maximal time step value and Courant number.

1. Go to RUN panel
2. Set the Simulation Time [s] to 1.5
3. Change Time Stepping to Automatic
4. Set the parameters:
Initial Δt [s] 1e-06
(solver will start computation with this value and adjust it in the next iterations)
Max Δt [s] 1e-03
Max Co [-] 1
Max Alpha Co [-] 0.5

This calculation will require high CPU usage and might take several hours to finish the calculation depending on the machine used. If you want to shorten the computation time, you can decrease the Simulation Time to 0.4 s – this is enough for purpose of further postprocessing.

39. Run - Output

We can control how often results should be saved on the hard drive. Only this data will be available for postprocessing.

1. Switch to Output tab
2. Set the Interval [s] to 1e-02

40. Run - CPU

Due to the non-newtonian fluid and 2 phase model, the calculation is very time-consuming. To speed up the calculation process increase the number of CPUs basing on your PC capability. We recommend using at least 6 cores for this tutorial. If you are using a free version you can use the contact form to Request 30-day Trial

Estimated computation time for 2 processors: over 5 hours

1. Switch to CPU tab
2. Use parallel mode
3. Increase the Number of processors
4. Click Run Simulation button

41. Monitor Solution - Residuals

When the calculation is finished we should see a similar residual plot.

42. Results - Slice

Slices tab appears next to Residuals. Under this tab, we can preview results on the defined slice planes. The results preview is available during the calculation and we can track it on a regular basis. The newly calculated time step will be actualized automatically as long as the time selector points to the latest time step.

1. Switch to Slices tab
2. Click ViewXY
3. Click Fit View to zoom the geometry
4. Choose alpha.HDPE field
5. Click Adjust range to data
6. Play with an animation buttons to track the results of analysis
7. Probes outputs are available to preview on the sheets next to Slice

43. Calculate - Mass Flow

After computations are finished, we will investigate the mass flow rate during the injection. We will compute it by using already existing data (this feature is not available as a monitor).

1. Go to Calculate panel
2. Switch to Boundary tab
3. Select Full Time Range
4. Select the boundary as housing_inlet
5. Select mass φ
6. Click on Calculate

44. Results - Mass Flow

1. Computed variables are displayed on the chart in the tab Integral of phi

45. Postprocessing – ParaView

We can do the complex visualization of our results with ParaView.

2. Click on Run ParaView

46. ParaView - Load the results

Load the results into the program.

1. Mark housing_injection.foam
2. Click Apply to load results into ParaView
3. Select the phase alpha.HDPE from the list
4. After loading results they will be shown in the 3D graphic window

47. ParaView – Results (I)

We will focus on the HPDE phase behavior and observe form filling during the injection.

1. Select Clip from top menu
2. Change the clip type to Scalar
3. Select alpha.HDPE from the scalar list
4. Click Refresh
5. Uncheck Invert option
6. Click Apply
7. Select the phase alpha.HDPE to preview
8. Play with the animation buttons to track the HDPE phase movement

48. ParaView – Results (II)

To show the housing together with the fluids we will follow these steps

1. Unhide housing_injection.foam
2. Set the coloring by the Solid Color
3. Set the opacity to 0.3

49. Final Results

Now we can see the final results with the original geometry.

This website uses cookies to offer you the best experience online. By continuing to use our website, you agree to the use of cookies.