Simplex space-time meshes in thermally coupled two-phase flow simulations of mold filling

Violeta Karyofylli,Loic Wendling,Michel Make,Norbert Hosters,Marek Behr

DOI: https://doi.org/10.48550/arXiv.1903.08710

2019-03-21

Abstract:The quality of plastic parts produced through injection molding depends on many factors. Especially during the filling stage, defects such as weld lines, burrs, or insufficient filling can occur. Numerical methods need to be employed to improve product quality by means of predicting and simulating the injection molding process. In the current work, a highly viscous incompressible non-isothermal two-phase flow is simulated, which takes place during the cavity filling. The injected melt exhibits a shear-thinning behavior, which is described by the Carreau-WLF model. Besides that, a novel discretization method is used in the context of 4D simplex space-time grids [2]. This method allows for local temporal refinement in the vicinity of, e.g., the evolving front of the melt [10]. Utilizing such an adaptive refinement can lead to locally improved numerical accuracy while maintaining the highest possible computational efficiency in the remaining of the domain. For demonstration purposes, a set of 2D and 3D benchmark cases, that involve the filling of various cavities with a distributor, are presented.

Computational Engineering, Finance, and Science,Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the complex flow and heat transfer problems in the mold - filling stage during the injection molding process. Specifically, the paper focuses on the two - phase flow simulation between highly viscous non - Newtonian fluids (such as molten plastics) and air during the mold - filling process. This process involves multiple physical phenomena, including but not limited to:

1. **Shear - thinning behavior**: The molten material exhibits shear - thinning characteristics, that is, its viscosity decreases as the shear rate increases.
2. **Interface capture**: It is necessary to accurately capture the interface between the molten material and air, especially in dynamically changing situations.
3. **Temperature coupling**: The viscosity of the molten material is affected by temperature, so the change of the temperature field needs to be considered.
4. **Wetting effect**: The wetting behavior of the molten material on the mold wall has an important impact on the filling effect.

In order to better understand and predict these complex physical phenomena, thereby improving product quality and reducing defects (such as weld lines, burrs or under - filling), the author proposes a new discretization method based on 4D simplex space - time grids. This method allows for local time refinement in the region of interest, thereby improving numerical accuracy while maintaining computational efficiency. Specifically, the main contributions of the paper include:

- **Introduction of the Carreau - WLF model**: Used to describe the shear - thinning behavior of the molten material and correct the temperature dependence through the WLF formula.
- **Use of the level - set method**: Used to capture the interface between the molten material and air.
- **4D simplex space - time grids**: Using this grid allows for adaptive refinement in time and space, especially near the molten front.
- **Strong - coupling partition method**: Ensure that the interactions between the Navier - Stokes equations, the level - set equations and the heat equations are accurately handled.

Through these methods, the paper aims to provide a more accurate and efficient numerical simulation tool for the injection molding process, thereby helping to optimize the production process and improve product quality.

### Formula summary

1. **Carreau - WLF model**:
\[
(\mu_{\text{eff}})_{\text{PZT}} = a_T \left[ \mu_\infty + (\mu_0 - \mu_\infty) \left( 1 + (a_T \lambda G_f)^a \right)^{\frac{n - 1}{a}} \right]
\]
where,
\[
G_f = \sqrt{2\varepsilon(u):\varepsilon(u)}
\]
and
\[
a_T = \exp \left[ \frac{c_1 (T_{\text{ref}} - T)}{c_2 + T - T_{\text{ref}}} \right]
\]

2. **Relationship between thermal conductivity and Prandtl number**:
\[
k_{\text{PZT}} = \frac{C_p \mu_\infty}{\text{Pr}}
\]
\[
\text{Pr} = \frac{C_p \mu_\infty}{k_0}
\]

3. **Navier - Stokes equations**:
\[
\rho_i \left( \frac{\partial u}{\partial t} + u\cdot\nabla u - f \right) - \nabla\cdot\sigma_i = 0
\]
\[
\nabla\cdot u = 0
\]
\[
\sigma_i(u, p) = -pI + 2 (\mu_{\text{eff}})_i \varepsilon(u)
\]
\[
\varepsilon(u) = \frac{1}{2} \left( \nabla u+(\nabla u)^T \right)
\]