A New Conformal Cooling Design Procedure for Injection Molding Based on Temperature Clusters and Multidimensional Discrete Models

This paper presents a new method for the automated design of the conformal cooling system for injection molding technology based on a discrete multidimensional model of the plastic part. The algorithm surpasses the current state of the art since it uses as input variables firstly the discrete map of temperatures of the melt plastic flow at the end of the filling phase, and secondly a set of geometrical parameters extracted from the discrete mesh together with technological and functional requirements of cooling in injection molds. In the first phase, the algorithm groups and classifies the discrete temperature of the nodes at the end of the filling phase in geometrical areas called temperature clusters. The topological and rheological information of the clusters along with the geometrical and manufacturing information of the surface mesh remains stored in a multidimensional discrete model of the plastic part. Taking advantage of using genetic evolutionary algorithms and by applying a physical model linked to the cluster specifications the proposed algorithm automatically designs and dimensions all the parameters required for the conformal cooling system. The method presented improves on any conventional cooling system design model since the cooling times obtained are analogous to the cooling times of analytical models, including boundary conditions and ideal solutions not exceeding 5% of relative error in the cases analyzed. The final quality of the plastic parts after the cooling phase meets the minimum criteria and requirements established by the injection industry. As an additional advantage the proposed algorithm allows the validation and dimensioning of the injection mold cooling system automatically, without requiring experienced mold designers with extensive skills in manual computing.


Introduction
Injection molding is today the most widely used plastic material processing technology in the industry. The injection molding process is characterized by its high versatility in the manufacture of plastic components by providing parts of complex geometry with good surface finish, high precision and low cost in industrial production [1,2]. Injection molding essentially consists of four stages: injection, packing, cooling, and ejection. During the injection process the plastic material is heated until obtaining a homogeneous melt which is introduced at high pressure into the mold cavity. When the filling of the cavity is almost finished, a subsequent packing pressure is applied which allows the remaining volume of the cavity to be filled, compensating in this way the shrinkage caused by the material cooling. When the material located on the cavity gate solidifies, the process begins of cooling the whole part by thermal exchange with the coolant flow. Once the part is rigid enough to be removed from the mold the part is ejected, and the cycle starts again [3]. In every cyclic manufacturing process cycle time is a key factor in the productivity of the process. In the injection molding process, the time of each cycle depends significantly on the time spent on cooling the molded part [4], this being around 70-80 percent of the cycle time [5][6][7]. Additionally, the quality of the molded part depends to a great extent on the characteristics of the mold cooling process [8,9]. An efficient cooling system plays a vital role in the technology of injection molding, since it optimizes the process of heat exchange between the molded part and the coolant flow. The requirements of the mold cooling process depend fundamentally on the topology of the plastic part, the design and dimensioning of the cooling channels and on the thermal properties of the mold material. In this line, the cooling of the molded part must be as uniform and balanced as possible, so that it makes it possible to eliminate defects such as sink marks, differential shrinkage, residual stress, warping, etc. [10,11]. The design of the cooling channels plays a critical role in the performance of the injection mold [12]. Traditionally, straight cooling channels prevent uniform and effective cooling in free form parts. With the advancement of the additive manufacturing technology the conformal cooling channels mean that the coolant flow can follow the geometry of the molded part, transferring in this way the heat from the cavity to the cooling channel efficiently. Another additional advantage of the conformal cooling channels is the reduction in the cooling time as well as a distribution improvement in the surface temperature of the part. Efficiency in the cooling process is then achieved by designing an optimal channel system that maximizes the process of thermal exchange between the part and the coolant flow. On the other hand, the temperature distribution at the end of the filling phase depends fundamentally on the geometrical features of the part [13]. A complex and non-optimized part topology can cause shearing in the injected material and a lack of uniformity in the surface temperatures. However, despite the importance of prior knowledge of the temperature distribution at the end of the injection phase no method of designing and dimensioning of the cooling system takes this information as input to obtain an optimal and accurate cooling system. Although some researchers focused on the optimization of conformal cooling channel design, there still is no a concrete approach for designing the best cooling system for the mold [14]. Additionally, the number of studies focused on obtaining the design parameters for conformal channels for free form parts automatically is very limited.
To solve these problems, this paper presents a new algorithm for obtaining the automated design of conformal cooling channels. The algorithm surpasses the current state of the art since it uses as input variables firstly the discrete map of temperatures of the melting front in each node of the surface mesh, and secondly a set of geometrical parameters of the part together with the technological and functional requirements of cooling in injection molds. Thus, the dimensioning of the main technological and geometrical parameters of the conformal cooling system is performed via an optimal and precise procedure. In addition, the algorithm allows the grouping, using a discrete format, of geometrical areas of the plastic part in clusters of temperature, from the nodal temperature of the melting front at the end of the filling phase. Finally, the temperature of each cluster coincides with the average temperature of the cluster nodes, which minimizes the standard deviation with respect to the rest of the cluster marginal temperatures. The technological and geometrical parameters obtained after the generation of the clusters is stored in a multidimensional discrete model of the plastic part. With this information it is possible to automatically dimension, by using expert algorithms of the genetic type, all of the parameters required for the cooling system, such as the diameter of the conformal channels and the optimum distance and positioning of each channel with respect to the mold cavity surface. The new algorithm uses the advantages of the process of additive manufacturing in rapid tooling to improve the cooling process in deep zones of the part, eliminating in this way quality failures and achieving as a result a complete and uniform cooling of the plastic part. The algorithm presented here allows the validation and dimensioning of the cooling system of the injection mold without requiring the help of mold makers and expert designers.

Background and Related Work
The advanced design of conformal cooling channels attracted the attention of designers and mold makers who centered their research on the conformal cooling channels' design and layout.
An advanced design of the cooling structure for an injection mold can improve the quality of the products and reduce the cycle time. Several studies highlight the fact that a cooling design using conformal cooling channels presents a better performance than a traditional layout with straight channels for cooling free form surfaces [15]. In a first study, Sachs et al. [16] used additive manufacturing with stainless steel powder and a polymer binder to obtain the cavity inserts, integrating a conformal channel lay out into a cylindrical part. In his research he contrasted the performance of the conformal system versus a conventional refrigeration made of steel by means of Computer Numerical Control technologies. The results show a more balanced temperature arrangement in the conformal approach. Nevertheless, the thermal features of the mold material vary using conventional steel, this being an important factor not quantified in the performance of the cooling system. Conforming to these results, Xu et al. [17] developed a new model for sizing the main parameters of the cooling system. In the research presented, Xu divided the mold surface into micro areas called cooling areas in order to obtain the global cooling system. In a second step, the cooling areas were decomposed into cooling units in a space of six design boundary rules for each unit. Finally, the results for each discrete cooling area were again regrouped in order to obtain the cooling system. The results obtained by Xu indicated a reduction of 15% in cycle time and 9% in part warping. Dimla et al. [18] described a procedure for modeling and improving the cooling and heating for injection molds by means of simulated models. Their methodology included a numerical simulation of the designed conformal cooling in order to find the best position for the conformal channels. Quiao et al. [19] applied a procedure for improving the Computer Aided Design model of the conformal channels, analyzing the cooling by means of a sensitivity study established on the boundary elements method and a hybrid optimizer. The sensitivity analysis together with boundary element method presented a good performance; however, they are not applicable to real industrial parts. Kitayama et al. [20] examined the cooling design of conformal channels, analyzing the cycle time and the warpage in industrial molded parts through numerical and experimental tests. The results confirmed that conformal cooling is effective due to the short cycle time and the warpage reduction. Kitayama et al. [21] performed a multi-objective design optimization in order to maximize the minimum weld line temperature for weld lines reduction and to minimize the clamping force, achieving as a result an increase in the quality of the product. Park et al. [22] presented the development of an injection mold with conformal cooling channels for manufacturing a free form industrial component with variable thickness. Their study indicated that conformal cooling channels decrease the cycle time by about 30% in comparison to traditional cooling systems. Park et al. [13] performed a method for maximizing the efficiency of conformal cooling channels manufactured with 3D printing technology. The research results showed that the cooling and cycle time can be reduced by more than 50%. Marques et al. [23] investigated two designs of conformal cooling: the parallel circuit and the serial circuit. Both proposals were evaluated by Computer Aided Engineering simulation to produce parts of polypropylene against traditional cooling circuits (linear channels). The results showed that the deformation of the product can be reduced significantly by using conformal cooling channels. Zink et al. [24] analyzed the effect of mold designing and mold material on cooling efficiency using numerical methods. The simulation model was adjusted to the measurement results by considering the gap between the mold inserts. More recently Jahan et al. [25] proposed a modeling procedure for designing optimized layouts of conformal cooling systems in injection molding. They applied the design of experiments (DOEs) methodology in order to analyze the influence of the main design dimensions, as well as the cross-section topologies. Berger et al. [26] studied the performance of three cooling layouts for cooling a complex-shaped injection molded PA6 GF35 part that exhibited a mass accumulation. The results of the experiments indicated that the 3D printed insert core could greatly improve the uniformity of the surface temperature in the area with defects and reduce the weld line of the manufactured product. Conformal cooling designs were also Polymers 2020, 12, 154 4 of 35 applied in molds with fast cooling cycles. This technique consists of heating the cavity during filling and cooling the rest of the cycle rapidly. Recent studies can be found in [27][28][29]. Cho et al. [29] studied a 3D printed insert core of an injection molding tool for reducing weld line defects on automotive crash pads with a high standard of appearance. The heat transfer experiments were conducted in order to compare the efficiency of two different types of building up methods for insert cores, and a heat transfer simulation was conducted to design a heating and cooling channel for the conformal surface temperature of the insert core and for thermal response efficiency. For conformal cooling uses in small scale molds Clemente et al. [30] introduced the flow channels svelteness as a total geometric dimensioning value for obtaining their layout.
Unfortunately, all the designs presented had the disadvantage of a totally manual execution. To resolve this issue, several researchers studied the possibility of automating the conformal cooling design for injection molding. To decrease the progressive rise in the coolant temperature from the start point to the final point of the channels lay out, Au et al. [31] designed a methodology to adapt the positioning of the conformal cooling channels to the mold cavity. Wang et al. [32] performed a new procedure for automated designing of the conformal cooling system. Unfortunately, this new procedure required a large quantity of information with regard to the coolant, cycle time, sizing of the cooling channels etc. The values were obtained by the user, and so were partially dependent on designers. Additionally, the cooling channel design developed by Wang was made up of complex networks. To resolve this issue, Wang et al. [33] presented a new methodology for obtaining conformal cooling designs by means of spiral channels. This arrangement improved the cooling channels and decreased production costs. Recently, Li et al. [14] proposed a topology optimization method for designing the conformal cooling layout for injection molds. In the course of the modeling process the cycle-averaged procedure was employed to facilitate the study of the cooling performance, and the boundary element method (BEM) was employed to resolve the governing algorithms. Unfortunately, the channel arrangement required the performance of a self-complied code to obtain the sensitivities.
Additive manufacturing can create three-dimensional lattices and scaffoldings with specific mechanical and thermal properties [34][35][36][37][38]. Several authors developed new cooling proposals using scaffoldings or porous structures. Porous structures can achieve a larger surface-area-to-volume ratio, which may further reduce the thermal inertia of an injection mold. Compared with conventional circular channels, porous structures exhibit better mechanical performance during the injection molding process. Some authors presented new proposals using uniform porous scaffolds with a cubic shape which are used to fill the cooling area uniformly. Au et al. [39] presented a balanced accuracy lattice construction for cooling channels design. The conclusions indicated that this arrangement supplied a more balanced cooling geometry in comparison to traditional cooling layouts and enhanced the mechanical strength of the cavities in the injection mold. In Au et al. [40] they followed up their work by performing a different design procedure comprising of a conformal cooling passageway with porous multi-connected features, established on the theory of duality. The conclusions of their research presented a more balanced mold cooling distribution in comparison with available conformal layouts. Crema et al. [41] presented a new mold insert designed using scaffoldings and validated by numerical simulations for rapid heating. Several other cell topologies, including cross and N-shape cells, were also investigated. Along these lines Brooks et al. [42] presented a study based on using conformal cooling layers designed using unit cells without supports. The results presented a decrease in cooling time of 26% compared to the traditional cooling systems. Unfortunately, there are still some drawbacks to using uniform porous structures to achieve conformal cooling. Firstly, the pressure drop in uniform porous structure is usually much higher than that of conventional cooling channels. Secondly, the distance between the cooling surface and the mold is not constant. To resolve these issues, Tang et al. [43] proposed a general design method of conformal porous structure for injection mold cooling that can achieve better cooling performance by reducing the maximum average mold temperature and the pressure drop. More recently, Mercado et al. [44] proposed a new layout for a conformal cooling system with the surface of a parametric lattice. The parameters of the cooling Polymers 2020, 12, 154 5 of 35 channels, such as the diameter of the incoming lattice channels, are automatically sized requiring as input only information of the material and the part topology. Even so, the topologies of the lattice arrangements [39][40][41][42][43][44] could originate too much turbulence and pressure drops in the coolant flow as it advances through the scaffolds.
Despite all the advantages of conformal cooling channels, their benefits are maximized only when an automated design method capable of optimizing the heat exchange process between the part and the coolant flow is carried out. The knowledge of the temperature at each point on the part surface at the beginning of the cooling phase helps to increase the accuracy in the dimensioning and arrangement of the channels with respect to the part surface. However, none of the methods presented in the state of the art include as an input variable the non-uniform distribution of temperatures on the surface of the part at the end of the filling phase. Considering a constant injection temperature in the whole part as input for the layout and sizing of the mold cooling can cause important design errors. At the end of the injection phase rheological problems such as an increase in the temperature of some areas of the part due to the shearing of the material, lack of surface thermal uniformity or hot spots caused by a thickness variability can occur. All these factors have a great influence on the geometrical design and sizing of the cooling system. Additionally, for designing cooling channels in injection molding it is essential to include geometrical, functional, and technological requirements that allow us to achieve, at the end of the cooling phase, a uniformity in the part surface temperatures eliminating warpage and residual stresses as well as reducing the cooling time.
To solve these problems this paper presents a new design system for the automated design of conformal cooling channels. The methodology for sizing the geometrical and technological variables of the conformal cooling system is applied, taking into account the geometric discretization of the plastic surface in temperature clusters and the advantages of using expert optimization algorithms of an evolutionary type. In this way the average temperature of the molten plastic front is used as an input parameter to the physical model for the dimensioning of the conformal cooling system, linking the physical model to the requirements of each cluster. It is thereby possible to obtain the magnitude of the geometrical and technological variables of the conformal cooling system associated with each cluster, obtaining as a result an optimum uniformity in the cooling of the molded part. The new proposal employs the benefits of the technology of additive manufacturing (rapid tooling), enhancing the cooling of the deep areas of the geometry and removing quality failures. The new system permits the validation and sizing of the cooling mold design without requiring skills in manual computing and experience on the part of the mold designer.

Discrete Model of the Plastic Part
The algorithm proposed in this paper is based on the previous geometrical analysis of the plastic part under study P ∈ R 3 (see Figure 1). Throughout this process the 3D geometry of the plastic part is discretized with the aim of obtaining a three-dimensional mesh P' ∈ R 3 (see Figure 2). The three-dimensional discrete mesh of the plastic part P' ∈ R 3 (see Figure 2) is composed of remarkable elements of the node type P ij = {X ij , Y ij , Z ij } (with the geometrical information of their Cartesian coordinates, see Figure 2) and facets or polygonal surfaces F i = {P i1 , P i2 , P i3 } ∈ P' f (with the geometrical information about the nodes that compose it, see Figure 2). For the methodology proposed in this paper the facets or polygonal surfaces are defined as triangular. Nevertheless, the scope of the proposed methodology can be extended to facets with different geometrical definitions. The accuracy of the mesh or average length of the facets of each of the plastic parts analyzed is establish according their thickness. Figure 2 shows the graphical representation of the plastic part virtual mesh model and in detail the geometry of the triangular facets and the nodes that define it.

Temperature Map of the Melting Front in the Plastic Part
With the aim of designing and dimensioning a conformal cooling system optimized and adapted to the geometrical and technological variables of the plastic part, the algorithm proposed in this paper establishes as an input variable the temperature of the melting front after the filling phase. From this technological variable it is possible to access the melt plastic flow temperature in each area and, in particular, in each node Pij ∈ P'n that forms the discrete geometry of the plastic part. Hence, the magnitude of the technological and geometrical variables of the conformal cooling system are adjusted to the temperature map of the melt plastic flow of the plastic part by optimizing and adapting the cooling system to each geometrical area, improving the final surface quality, the structural behavior, and the functional requirements. It also manages to reduce the residual stress acquired by the gradient of temperatures between different areas of the plastic part as well as a reduction in the cycle time, improving the productivity of the manufacturing process.

Mold Cavity Filling Stage Modeling
Given the plastic part virtual mesh model P' ∈ ℝ 3 and the thermoplastic material of which it is composed, the temperature of the thermoplastic filling flow is defined by a numerical simulation of the filling process in the injection mold cavity. This numerical analysis was carried out using the numerical simulation software Moldex 3D R16 for injection molding [45]. Firstly, the plastic part virtual mesh model under study is imported into the numerical simulation software. For this task the discrete mesh P' ∈ ℝ 3 must be previously defined, together with the Cartesian coordinates Pfill = {Xfill, Yfill, Zfill} ∈ ℝ 3 and geometrical dimensions of the injection point or gate Rgate (see Figure 3).

Temperature Map of the Melting Front in the Plastic Part
With the aim of designing and dimensioning a conformal cooling system optimized and adapted to the geometrical and technological variables of the plastic part, the algorithm proposed in this paper establishes as an input variable the temperature of the melting front after the filling phase. From this technological variable it is possible to access the melt plastic flow temperature in each area and, in particular, in each node Pij ∈ P'n that forms the discrete geometry of the plastic part. Hence, the magnitude of the technological and geometrical variables of the conformal cooling system are adjusted to the temperature map of the melt plastic flow of the plastic part by optimizing and adapting the cooling system to each geometrical area, improving the final surface quality, the structural behavior, and the functional requirements. It also manages to reduce the residual stress acquired by the gradient of temperatures between different areas of the plastic part as well as a reduction in the cycle time, improving the productivity of the manufacturing process.

Mold Cavity Filling Stage Modeling
Given the plastic part virtual mesh model P' ∈ ℝ 3 and the thermoplastic material of which it is composed, the temperature of the thermoplastic filling flow is defined by a numerical simulation of the filling process in the injection mold cavity. This numerical analysis was carried out using the numerical simulation software Moldex 3D R16 for injection molding [45]. Firstly, the plastic part virtual mesh model under study is imported into the numerical simulation software. For this task the discrete mesh P' ∈ ℝ 3 must be previously defined, together with the Cartesian coordinates Pfill = {Xfill, Yfill, Zfill} ∈ ℝ 3 and geometrical dimensions of the injection point or gate Rgate (see Figure 3). . Discrete geometry P' ∈ R 3 of the plastic part and detail of a triangular facet F i and the nodes P i1 , P i2 , P i3.

Temperature Map of the Melting Front in the Plastic Part
With the aim of designing and dimensioning a conformal cooling system optimized and adapted to the geometrical and technological variables of the plastic part, the algorithm proposed in this paper establishes as an input variable the temperature of the melting front after the filling phase. From this technological variable it is possible to access the melt plastic flow temperature in each area and, in particular, in each node P ij ∈ P' n that forms the discrete geometry of the plastic part. Hence, the magnitude of the technological and geometrical variables of the conformal cooling system are adjusted to the temperature map of the melt plastic flow of the plastic part by optimizing and adapting the cooling system to each geometrical area, improving the final surface quality, the structural behavior, and the functional requirements. It also manages to reduce the residual stress acquired by the gradient of temperatures between different areas of the plastic part as well as a reduction in the cycle time, improving the productivity of the manufacturing process.

Mold Cavity Filling Stage Modeling
Given the plastic part virtual mesh model P' ∈ R 3 and the thermoplastic material of which it is composed, the temperature of the thermoplastic filling flow is defined by a numerical simulation of the filling process in the injection mold cavity. This numerical analysis was carried out using the numerical simulation software Moldex 3D R16 for injection molding [45]. Firstly, the plastic part virtual mesh model under study is imported into the numerical simulation software. For this task the discrete mesh As shown in Figures 3 and 4, the numerical simulation software for injection molds Moldex 3D R16 discretizes the internal volume of the plastic part from the discrete mesh P' ∈ ℝ 3 using SOLID 186 tetrahedron-type elements with average length equal to the precision ε of P' ∈ ℝ 3 . These elements have a quadratic thermal behavior, being composed of 10 nodes (four nodes on the tetrahedron vertex and six nodes on the midpoints of the edges). The tetrahedrical elements include three degrees of freedom per node, these three nodal translations being in the X, Y and Z directions. Furthermore, to improve the numerical modeling of the contact area between the injection mold surface and the melting front a set of elements with a variable and adaptable size called boundary layers was defined (see Figure 4). This set of elements covers the whole contact surface between the injection mold and the plastic part. Each element is composed of five layers of SOLID 186 prismatic elements whose layer height is established by an offset ratio, which relates the height of the prismatic element and the average size of the element. The offset ratio used for each of the plastic parts analyzed in this paper was equal to 0.100. This value corresponds to the height of the prismatic element in the boundary layers, being equal to the average length of the prismatic element multiplied by the offset ratio. These elements have a quadratic thermal behavior, being composed of 15 nodes (6 nodes in the vertex of the prism and nine nodes in the midpoints) and three degrees of freedom in each node: translation nodal in the X, Y and Z directions. Figure 4 shows the precision, the number of nodes and elements and element type of each of the plastic parts meshes analyzed in this paper. Next, the thermoplastic material of the plastic part under study is established. The numerical simulation software Moldex 3D R16 has a wide database of thermoplastic materials from different As shown in Figures 3 and 4, the numerical simulation software for injection molds Moldex 3D R16 discretizes the internal volume of the plastic part from the discrete mesh P' ∈ R 3 using SOLID 186 tetrahedron-type elements with average length equal to the precision ε of P' ∈ R 3 . These elements have a quadratic thermal behavior, being composed of 10 nodes (four nodes on the tetrahedron vertex and six nodes on the midpoints of the edges). The tetrahedrical elements include three degrees of freedom per node, these three nodal translations being in the X, Y and Z directions. Furthermore, to improve the numerical modeling of the contact area between the injection mold surface and the melting front a set of elements with a variable and adaptable size called boundary layers was defined (see Figure 4). This set of elements covers the whole contact surface between the injection mold and the plastic part. Each element is composed of five layers of SOLID 186 prismatic elements whose layer height is established by an offset ratio, which relates the height of the prismatic element and the average size of the element. The offset ratio used for each of the plastic parts analyzed in this paper was equal to 0.100. This value corresponds to the height of the prismatic element in the boundary layers, being equal to the average length of the prismatic element multiplied by the offset ratio. These elements have a quadratic thermal behavior, being composed of 15 nodes (6 nodes in the vertex of the prism and nine nodes in the midpoints) and three degrees of freedom in each node: translation nodal in the X, Y and Z directions. Figure 4 shows the precision, the number of nodes and elements and element type of each of the plastic parts meshes analyzed in this paper. As shown in Figures 3 and 4, the numerical simulation software for injection molds Moldex 3D R16 discretizes the internal volume of the plastic part from the discrete mesh P' ∈ ℝ 3 using SOLID 186 tetrahedron-type elements with average length equal to the precision ε of P' ∈ ℝ 3 . These elements have a quadratic thermal behavior, being composed of 10 nodes (four nodes on the tetrahedron vertex and six nodes on the midpoints of the edges). The tetrahedrical elements include three degrees of freedom per node, these three nodal translations being in the X, Y and Z directions. Furthermore, to improve the numerical modeling of the contact area between the injection mold surface and the melting front a set of elements with a variable and adaptable size called boundary layers was defined (see Figure 4). This set of elements covers the whole contact surface between the injection mold and the plastic part. Each element is composed of five layers of SOLID 186 prismatic elements whose layer height is established by an offset ratio, which relates the height of the prismatic element and the average size of the element. The offset ratio used for each of the plastic parts analyzed in this paper was equal to 0.100. This value corresponds to the height of the prismatic element in the boundary layers, being equal to the average length of the prismatic element multiplied by the offset ratio. These elements have a quadratic thermal behavior, being composed of 15 nodes (6 nodes in the vertex of the prism and nine nodes in the midpoints) and three degrees of freedom in each node: translation nodal in the X, Y and Z directions. Figure 4 shows the precision, the number of nodes and elements and element type of each of the plastic parts meshes analyzed in this paper. Next, the thermoplastic material of the plastic part under study is established. The numerical simulation software Moldex 3D R16 has a wide database of thermoplastic materials from different Next, the thermoplastic material of the plastic part under study is established. The numerical simulation software Moldex 3D R16 has a wide database of thermoplastic materials from different manufacturers, so the material selection is performed during the numerical simulation configuration of the injection mold cavity filling. The thermoplastic materials used for the case studies analyzed in this paper are: Cycolac MG47 (ABS), Finalloy EBP-830 (PP) and Lexan HF1110R (PC). Finally, the technological parameters that govern the numerical analysis of the filling phase are shown in Table 1. The whole set of input technological parameters corresponding to the Cartesian coordinates of the injection points P fill = {X fill , Y fill , Z fill } ∈ R 3 are determined according to [4], in order to automate the process of numerical simulation during the mold cavity filling. The input technological parameters required for the algorithm were the dimensions of the gate R gate , the filling time of the mold cavity T fill , the packing time T pack and the cooling time T cooling . The remaining technological parameters: melt temperature T melt , mold temperature T mold , plastic part ejection temperature T eject and maximum injection pressure P inj are automatically defined with the selection of the thermoplastic material, as the manufacturer recommends their magnitude for an optimal cavity filling process of the injection mold. Each of the variables is presented in Tables 1-3.
Figures 5 and 6 graphically show the temperature map distribution of the melt plastic flow TMF ij ∈ R ( • C) (see Equation (1)) in each node of the plastic part virtual mesh model under study. The result of the temperature of the thermoplastic filling flow obtained for each node P ij ∈ P' n belonging to P' ∈ R 3 is exported in ASCII format for a later use as parameter for the calculation of the geometrical and technological variables of the conformal cooling system in the proposed algorithm. In this way, by combining the result of the technological variable TMF ij ∈ R obtained from the numerical simulation and the Cartesian coordinates of each node P ij = {X ij , Y ij , Z ij } ∈ P' n , it is possible to establish a hybrid matrix PH ij ∈ R 4 (see Figure 6) of range 4·n x m that combines geometrical information of the part with technological information about its manufacture, conforming a multidimensional discrete model of the plastic part. (see Equation (1)). L p-c m Separation distance between the surface of the plastic part and the conformal cooling channels x (5) . . . x (11) Table 3. Set-up of the genetic optimization algorithm.

Genetic Algorithm Options Structure Defined Set-up Description
Constraints Constraints restrict the value of the optimization problem variables, allowing the solution to converge on the defined domain.

Reproduction Two Points
Future generations are created from the arithmetic mean of the two random pairs of data of the previous generation population.

Mutation 2%
Definition of the number of individuals of a generation that can be mutated in future generations.
Crossover factor 80% Definition of the number of individuals that can be crossed in future generations.

Migration
Direction: Forward Migration Interval: 20 Migration factor: 30% Defines the migration direction of individuals between generations and the individuals fraction that can be migrated between the computational domains.
(1)  (1)) in each node of the plastic part virtual mesh model under study. The result of the temperature of the thermoplastic filling flow obtained for each node Pij ∈ P'n belonging to P' ∈ ℝ 3 is exported in ASCII format for a later use as parameter for the calculation of the geometrical and technological variables of the conformal cooling system in the proposed algorithm. In this way, by combining the result of the technological variable TMFij ∈ ℝ obtained from the numerical simulation and the Cartesian coordinates of each node Pij = {Xij, Yij, Zij} ∈ P'n, it is possible to establish a hybrid matrix PHij ∈ ℝ 4 (see Figure 6) of range 4·n x m that combines geometrical information of the part with technological information about its manufacture, conforming a multidimensional discrete model of the plastic part. (see Equation (1)).    (1)) in each node of the plastic part virtual mesh model under study. The result of the temperature of the thermoplastic filling flow obtained for each node Pij ∈ P'n belonging to P' ∈ ℝ 3 is exported in ASCII format for a later use as parameter for the calculation of the geometrical and technological variables of the conformal cooling system in the proposed algorithm. In this way, by combining the result of the technological variable TMFij ∈ ℝ obtained from the numerical simulation and the Cartesian coordinates of each node Pij = {Xij, Yij, Zij} ∈ P'n, it is possible to establish a hybrid matrix PHij ∈ ℝ 4 (see Figure 6) of range 4·n x m that combines geometrical information of the part with technological information about its manufacture, conforming a multidimensional discrete model of the plastic part. (see Equation (1)).   . Definition of the hybrid matrix PH ij from the discrete geometry P' n ∈ R 3 and from the temperature of the melt plastic front TMF ij in each node.

Clustering of Discrete Geometrical Areas Via Melting Front Temperature
The temperature definition of the thermoplastic filling flow establishes the temperature TMF ij ∈ R, in each node P ij ∈ P' n during the filling phase. To dimension in an optimized manner, the main technological and geometrical variables of the channels that conform to the cooling system, the nodes P ij ∈ P' n of the discrete mesh P' ∈ R 3 are grouped based on the temperature of the thermoplastic filling flow TMF ij ∈ R. Hence, it is feasible to dimension the elements of the cooling system adapted to the technological requirements of the plastic part to be molded. Thus, it is possible to perform discrete groupings of geometrical areas of the plastic part C k (see Figure 7 and Equation (2)) whose nodal temperature of the melting front is similar.

Clustering of Discrete Geometrical Areas Via Melting Front Temperature
The temperature definition of the thermoplastic filling flow establishes the temperature TMFij ∈ ℝ, in each node Pij ∈ P'n during the filling phase. To dimension in an optimized manner, the main technological and geometrical variables of the channels that conform to the cooling system, the nodes Pij ∈ P'n of the discrete mesh P' ∈ ℝ 3 are grouped based on the temperature of the thermoplastic filling flow TMFij ∈ ℝ. Hence, it is feasible to dimension the elements of the cooling system adapted to the technological requirements of the plastic part to be molded. Thus, it is possible to perform discrete groupings of geometrical areas of the plastic part Ck (see Figure 7 and Equation (2)) whose nodal temperature of the melting front is similar. Where nclusters represents the number of discrete groups of geometrical regions of the plastic part, Ck the index of the discrete groups or clusters and Tk (°C) the average temperature of the melt plastic flow for each discrete group or cluster. In this way the range of the hybrid matrix PHij ∈ ℝ 6 becomes 6·n x m, expanding the information of the discrete multimodal model to include technological and geometrical parameters obtained after the generation of the clusters. Therefore, it is possible to establish from the hybrid matrix PHij ∈ ℝ 6 the index Ck and the average temperature of the molten plastic front Tk of the discrete clusters to which each node Pij ∈ P'n belongs.
The discrete groups of geometrical areas of the plastic part are established by means of the statistical algorithm of generation of k-means clusters [46,47]. This algorithm has as a main objective the division of a set of n observations into a determined number of k subsets or clusters. So each observation belongs to the subset or clusters whose average value is closer. For the application of the k-means algorithm to the methodology proposed in this paper the temperature of the melt plastic flow TMFij ∈ ℝ is established as the initial set of observations, which in turn are related to the Cartesian coordinates of each node Pij ∈ P'n by the hybrid matrix PHij ∈ ℝ 6 .The temperature of the melt plastic flow TMFij ∈ ℝ is a vector (vector column of the hybrid matrix PHij ∈ ℝ 6 ), and from it the subsets or clusters Ck are constructed in order to minimize the sum of the squares of the average temperatures of the melt plastic flow in each of them (see Equation (3)). Where n clusters represents the number of discrete groups of geometrical regions of the plastic part, C k the index of the discrete groups or clusters and T k ( • C) the average temperature of the melt plastic flow for each discrete group or cluster. In this way the range of the hybrid matrix PH ij ∈ R 6 becomes 6·n x m, expanding the information of the discrete multimodal model to include technological and geometrical parameters obtained after the generation of the clusters. Therefore, it is possible to establish from the hybrid matrix PH ij ∈ R 6 the index C k and the average temperature of the molten plastic front T k of the discrete clusters to which each node P ij ∈ P' n belongs.
The discrete groups of geometrical areas of the plastic part are established by means of the statistical algorithm of generation of k-means clusters [46,47]. This algorithm has as a main objective the division of a set of n observations into a determined number of k subsets or clusters. So each observation belongs to the subset or clusters whose average value is closer. For the application of the k-means algorithm to the methodology proposed in this paper the temperature of the melt plastic flow TMF ij ∈ R is established as the initial set of observations, which in turn are related to the Cartesian coordinates of each node P ij ∈ P' n by the hybrid matrix PH ij ∈ R 6 .The temperature of the melt plastic flow TMF ij ∈ R is a vector (vector column of the hybrid matrix PH ij ∈ R 6 ), and from it the subsets or clusters C k are constructed in order to minimize the sum of the squares of the average temperatures of the melt plastic flow in each of them (see Equation (3)).
To establish an optimal number of discrete groups of geometrical areas of the plastic part or clusters, the Davies-Bouldin clustering evaluation criterion [48] is included. The Davies-Bouldin criterion is based on a ratio of within-cluster and between-cluster distances. The Davies-Bouldin index DB is defined as shown by Equation (4).
where D u,v is the within-to-between cluster distance ratio for the u and v clusters. Equation (5) shows this parameter in mathematical terms.
where d u is the average distance between each point in the u cluster and the centroid of the u cluster, d v is the average distance between each point in the v cluster and the centroid of the v cluster, and d u,v is the Euclidean distance between the centroids of the u and v clusters. The maximum value of D u,v represents the worst-case within-to-between cluster ratio for cluster u. The optimal clustering solution has the smallest Davies-Bouldin index value. In a similar manner, in order to improve the convergence criterion of the algorithm a vector V k ∈ R was proposed, from which the optimum number of discrete groups of geometrical areas of the plastic part or clusters will be established (see Equation (6)).
V k = n clusters min . . . n clusters max ∈ R where n clusters min represents the minimum number of optimal clusters and n clusters max represents the maximum number of optimal clusters. For the algorithm established in this manuscript the magnitude of these parameters is 2 and 6, respectively.

Dimensioning of the Conformal Cooling Channels via Genetic Algorithms
The calculation process of the technological and geometrical variables magnitude of conformal cooling channels via genetic optimization algorithms is described in the present section. The use of this type of algorithms allows that the geometrical design of the conformal cooling system to be adapted to a set of boundary conditions of cooling in injection molds. In this way it is possible to optimize in turn the heat exchange process between the coolant flow and the injection mold cavity. Table 2 describes the technological and geometrical variables of the conformal cooling channels that the genetic algorithm proceeds to optimize. It should be noted that this dimensioning process is carried out according to the physical model presented and, in turn, it is applied taking into account each of the geometrical discretizations or clusters (see Figure 7) previously defined, i.e., since each geometrical discretization or cluster C k has an associated average temperature of the thermoplastic filling flow T k , the present dimensioning process is adjusted to the technological variables and conditions of each geometrical discretization or cluster C k .
The magnitude of the thermal and technological properties of the injection mold material as well as the plastic part thermoplastic material were established from the manufacturer´s information. In addition, it should be noted that the magnitude of the variable temperature of the thermoplastic filling flow T melt ( • C) is determined from the magnitude of the average temperature of the thermoplastic filling flow T k ( • C), for each geometrical discretization or previously defined C k cluster.
As shown in Table 2, the number of variables of the genetic optimization algorithm for dimensioning the geometrical and technological variables of the cooling channels is not constant. The number of variables is a function of the number of previously defined geometrical clusters (see Figure 7) (minimum 2, maximum 6). In this way, it is established that the variable inlet temperature of the coolant flow T coolant ( • C), the separation distance between cooling conformal channels L conformal (m) and the diameter of the conformal cooling channels D conformal (m) are constant for the whole domain of the plastic part regardless of the geometrical cluster discretization to which they belong (see Figure 7). However, the variable of separation distance between the surface of the plastic part and the conformal cooling channels varies depending on the geometrical discretization or clusters (see Figure 7) to which they belong. In this way it is possible to develop a conformal cooling system that allows a uniform temperatures map to be obtained after the cooling phase around the plastic part surface.
Genetic optimization algorithms demand the implementation of a set of fitness functions to be optimized. These functions are related to the boundary conditions defined during the conformal cooling channels design process. Furthermore, constraint limits are defined to meet the additive manufacturing requirements. Table 3 describes the set-up of the genetic optimization algorithm for the theoretical model proposed.
Where T mold min ( • C) and T mold max ( • C) are defined as suggested values of the thermoplastic material manufacturer, L min (mm), L max (mm), D max (mm) and D min (mm) are the maximum and minimum separation distance between any element of the conformal cooling system (including the plastic part) and the maximum and minimum diameter of the conformal cooling channels, respectively. The magnitude of these geometrical parameters is recommended by the metallic material manufacturer in order to achieve the additive manufacturing criteria without using supports [43]. Table 4 describes the constraint values of the optimization problem variables for each material used in this paper. To finalize the genetic optimization algorithm definition, the three objective functions that compose it are defined. The first fitness function (see Equation (7)) reduces the sum of differences between the temperature of the mold cavity surface T mold ( • C), x(1) (optimization problem variable) and its final magnitude according to the proposed theoretical model (see Appendix A) for each geometrical discretization or cluster of the plastic part (see Figure 7).
where k represents the index of each geometrical discretization or clusters and T k ( • C) the average temperature of the thermoplastic filling flow for each geometrical discretization or clusters. The second fitness function (see Equation (8)) minimizes the technological variable cooling time t cooling ( • C), allowing the production process to be more effective and decreasing costs related to it. Furthermore, residual stress on plastic parts after the cooling phase is an important issue that usually is associated with the final quality of the product. Thus, the established constraints of the technological variables associated with the cooling process temperatures avert this type of defects from happening in the final product.
Fitness Function 2 : Min t cooling = Min The third fitness function (see Equation (9)) minimizes the sum of the differences between the average temperatures obtained on the surface of the injection mold, according to the proposed physical model (see Appendix A), for each geometrical discretization or clusters of the plastic part (see Figure 7). The definition of this objective function is to obtain a conformal cooling system whose temperature map after the cooling phase and cycle time is uniform for the whole domain of the plastic part.
Fitness Function 3 : min Finally, the results achieved by the genetic optimization algorithm in dimensioning the geometrical and technological variables of the conformal cooling channels are analyzed by the numerical simulation software Moldex 3D R16. In this way it is possible to validate the design of the resulting conformal cooling system.

Geometrical Design of the Conformal Cooling Channels Pathing
The algorithm for obtaining the automated geometrical design of the conformal cooling channels takes as a basis the results of the dimensioning expert algorithms and the advantages of an arrangement of channels in zig-zag. The conformal cooling channels design begins with the generation of a set of flat sections S conformal ∈ R 2 (see Equation (10)) obtained by intersecting the set of planes P conformal ∈ R 3 (see Figure 8) with the discrete geometry of the plastic part P' ∈ R 3 (see Figure 2). The planar sections S conformal ∈ R 2 are separated from each other by the parametrized distance L conformal obtained from the dimensioning genetic optimization algorithms. Similarly, each planar section S conformal (i) ∈ R 2 is defined by a set of n conformal (i) ∈ R 3 nodes generated by the intersection of the set of P conformal ∈ R 3 and the edges of the facets that conform to the discrete geometry of the plastic part P' ∈ R 3 (see Equation (10)). The normal direction that defines the set of planes P conformal ∈ R 3 is coincident with the longitudinal direction of the plastic part (see Figure 8) established from the edge of the bounding box that presents a larger dimension related to the layout of the zig-zag conformal cooling channels.
The geometrical design of the main plates, cavity, and core, which conforms the injection mold, is defined using the analysis of the plastic part demoldability. The main components of the injection mold cooling system have to be distributed uniformly along the main plates of the mold in order to optimize the cooling process. With the objective of distributing the conformal cooling channel paths in a uniform way along the cavity and core plates, the nodes n conformal ∈ R 3 belonging to the flat sections S conformal ∈ R 2 are classified according to their demoldability. The classification of the nodes n conformal ∈ R 3 is established by evaluating the demoldability of the facets P' f ∈ R 3 in which the nodes are located. The demoldability of the facets P' f ∈ R 3 of the discrete geometry P' ∈ R 3 is obtained from [1] (see Figure 9). As shown in Figure 10, with this classification two subset of nodes are determined, n conformal cav ∈ R 3 and n conformal cor ∈ R 3 (see Equation (11)), each of them defining the geometry of the cooling channel paths housed in the cavity and core plates respectively.
∀ P conformal (i) ∈ R 2 ∃ n conformal (i) ∈ R 2 n conformal (i) = n conformal cav (i) n conformal cor (i) (11) coincident with the longitudinal direction of the plastic part (see Figure 8) established from the edge of the bounding box that presents a larger dimension related to the layout of the zig-zag conformal cooling channels. The geometrical design of the main plates, cavity, and core, which conforms the injection mold, is defined using the analysis of the plastic part demoldability. The main components of the injection Figure 8. Definition of S conformal (i) ∈ R 2 , n conformal (i) ∈ R 2 and P conformal (i) ∈ R 3 for the generation of the layout of the conformal cooling channels.
Polymers 2020, 12, 154 14 of 34 mold cooling system have to be distributed uniformly along the main plates of the mold in order to optimize the cooling process. With the objective of distributing the conformal cooling channel paths in a uniform way along the cavity and core plates, the nodes nconformal ∈ ℝ 3 belonging to the flat sections Sconformal ∈ ℝ 2 are classified according to their demoldability. The classification of the nodes nconformal ∈ ℝ 3 is established by evaluating the demoldability of the facets P'f ∈ ℝ 3 in which the nodes are located. The demoldability of the facets P'f ∈ ℝ 3 of the discrete geometry P' ∈ ℝ 3 is obtained from [1] (see Figure 9). As shown in Figure 10, with this classification two subset of nodes are determined, nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 (see Equation (11)), each of them defining the geometry of the cooling channel paths housed in the cavity and core plates respectively.  The set of elements that make up the cooling system of an injection mold are located at a certain distance from the surface of the plastic part. This value is defined on the basis of the geometrical and technological requirements of the plastic part as well as of the injection mold. The methodology presented in this paper focuses on adapting the design of the conformal cooling channels by generating parameterized paths located to variable distances regarding the part surface geometry, taking into account the requirements of the manufacturing process. Therefore, the separation distance between each of the nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 with respect to the surface of the plastic part mold cooling system have to be distributed uniformly along the main plates of the mold in order to optimize the cooling process. With the objective of distributing the conformal cooling channel paths in a uniform way along the cavity and core plates, the nodes nconformal ∈ ℝ 3 belonging to the flat sections Sconformal ∈ ℝ 2 are classified according to their demoldability. The classification of the nodes nconformal ∈ ℝ 3 is established by evaluating the demoldability of the facets P'f ∈ ℝ 3 in which the nodes are located. The demoldability of the facets P'f ∈ ℝ 3 of the discrete geometry P' ∈ ℝ 3 is obtained from [1] (see Figure 9). As shown in Figure 10, with this classification two subset of nodes are determined, nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 (see Equation (11)), each of them defining the geometry of the cooling channel paths housed in the cavity and core plates respectively.  The set of elements that make up the cooling system of an injection mold are located at a certain distance from the surface of the plastic part. This value is defined on the basis of the geometrical and technological requirements of the plastic part as well as of the injection mold. The methodology presented in this paper focuses on adapting the design of the conformal cooling channels by generating parameterized paths located to variable distances regarding the part surface geometry, taking into account the requirements of the manufacturing process. Therefore, the separation distance between each of the nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 with respect to the surface of the plastic part Figure 10. Definition of n conformal cav ∈ R 3 y n conformal cor ∈ R 3 for the generation of the pathing of the conformal cooling channels.
The set of elements that make up the cooling system of an injection mold are located at a certain distance from the surface of the plastic part. This value is defined on the basis of the geometrical and technological requirements of the plastic part as well as of the injection mold. The methodology presented in this paper focuses on adapting the design of the conformal cooling channels by generating parameterized paths located to variable distances regarding the part surface geometry, taking into account the requirements of the manufacturing process. Therefore, the separation distance between each of the n conformal cav ∈ R 3 and n conformal cor ∈ R 3 with respect to the surface of the plastic part is defined by its proximity to the nodes P' n ∈ R 3 belonging to the hybrid matrix PH ij ∈ R 6 (see Figure 11). In turn, each cluster C k (see Figure 7) of the hybrid matrix PH ij ∈ R 6 has an associated separation distance between the cooling channels and the L P-C plastic part. The value of L P-C is defined during the dimensioning process by using genetic optimization algorithms. In this way, once the separation distance between the nodes n conformal cav ∈ R 3 and n conformal cor ∈ R 3 and the discrete geometry of the plastic part P' ∈ R 3 are determined, a geometrical bidirectional offset is designed following the normal direction to the plastic part surface in an outward sense (see Figure 11).
Polymers 2020, 12,154 15 of 34 distance between the cooling channels and the LP-C plastic part. The value of LP-C is defined during the dimensioning process by using genetic optimization algorithms. In this way, once the separation distance between the nodes nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 and the discrete geometry of the plastic part P' ∈ ℝ 3 are determined, a geometrical bidirectional offset is designed following the normal direction to the plastic part surface in an outward sense (see Figure 11). Subsequently, as a result of the 2D offset geometrical operation the nodes nconformal cav ∈ ℝ 3 and nconformal cor ∈ ℝ 3 are defined and displaced a distance LP-C with respect to the surface of the plastic part. The nodes nconformal cav ∈ ℝ 3 y nconformal cor ∈ ℝ 3 define the axes of the main paths of the cooling channels. To achieve this, a set of smoothed lines are modeled by using rounding radius with a diameter equal of those to the conformal cooling channels Dconformal. Finally, the algorithm performs a join operation in order to link all the individual lines designing the zig-zag path. The path starts from the polylines that are sited at the end of the plastic part that coincides with the inlet and outlet of the conformal cooling system (see Figure 12). A sweeping algorithm using a circular profile predimensioned by the expert system performs the final layout of the optimal conformal cooling distribution (see Figure 12). Subsequently, as a result of the 2D offset geometrical operation the nodes n conformal cav ∈ R 3 and n conformal cor ∈ R 3 are defined and displaced a distance L P-C with respect to the surface of the plastic part. The nodes n conformal cav ∈ R 3 y n conformal cor ∈ R 3 define the axes of the main paths of the cooling channels. To achieve this, a set of smoothed lines are modeled by using rounding radius with a diameter equal of those to the conformal cooling channels D conformal . Finally, the algorithm performs a join operation in order to link all the individual lines designing the zig-zag path. The path starts from the polylines that are sited at the end of the plastic part that coincides with the inlet and outlet of the conformal cooling system (see Figure 12). A sweeping algorithm using a circular profile predimensioned by the expert system performs the final layout of the optimal conformal cooling distribution (see Figure 12).

Numerical Analysis of the Conformal Cooling System
This section describes the configuring process for performing the numerical analysis of the study cases. The main objective of the analysis is to evaluate and validate the thermal performance of the results obtained according to the technological requirements demanded in the injection molding industry. The thermal and rheological numerical simulations analyze the behavior of the coolant flow along the conformal cooling system and the thermal exchange between it and the plastic part during the cooling phase. The set of thermal and rheological simulations are performed in the commercial numerical software Moldex 3D R16 [45]. As shown in Figure 13, the set of solid elements analyzed was divided into three main domains: conformal cooling channels (coolant domain), injection mold (solid domain alloy steel 1.2709), filling system (thermoplastic solid domain) and plastic part (thermoplastic solid domain). The hypotheses defined in the numerical simulations performed are as follows: • The materials used were: Water (Pure) for the coolant domain, steel alloy 1.2709 for the injection mold domain and the thermoplastics Cycolac MG47 (ABS), Finalloy EBP-830 (PP) and Lexan HF1110R (PC) for the domain of the plastic part and filling system. Table 5 shows the magnitude of the main properties of these materials.

•
Since the complete cooling process of the plastic part is analyzed over time, the type of numerical analysis used is "Cooling transient".

•
The total time of the numerical analysis in each simulation covers 30 s, with a magnitude of multiple time steps output equal to 10. For each time step the numerical calculation software stores the solution obtained. In this way, the complete evolution of the temperature map of the plastic part can be evaluated throughout the cooling phase.

•
The physical and thermal properties of the coolant flow along the conformal cooling system are analyzed using the Run 3D cooling channels analysis option • The type of solver used for the numerical analysis of the cooling phase is of the maximum variation of mold temperature type, the temperature difference being equal to 1° and the maximum number of cycles C being 10.

Numerical Analysis of the Conformal Cooling System
This section describes the configuring process for performing the numerical analysis of the study cases. The main objective of the analysis is to evaluate and validate the thermal performance of the results obtained according to the technological requirements demanded in the injection molding industry. The thermal and rheological numerical simulations analyze the behavior of the coolant flow along the conformal cooling system and the thermal exchange between it and the plastic part during the cooling phase. The set of thermal and rheological simulations are performed in the commercial numerical software Moldex 3D R16 [45]. As shown in Figure 13, the set of solid elements analyzed was divided into three main domains: conformal cooling channels (coolant domain), injection mold (solid domain alloy steel 1.2709), filling system (thermoplastic solid domain) and plastic part (thermoplastic solid domain). The hypotheses defined in the numerical simulations performed are as follows: • The materials used were: Water (Pure) for the coolant domain, steel alloy 1.2709 for the injection mold domain and the thermoplastics Cycolac MG47 (ABS), Finalloy EBP-830 (PP) and Lexan HF1110R (PC) for the domain of the plastic part and filling system. Table 5 shows the magnitude of the main properties of these materials.

•
Since the complete cooling process of the plastic part is analyzed over time, the type of numerical analysis used is "Cooling transient".

•
The total time of the numerical analysis in each simulation covers 30 s, with a magnitude of multiple time steps output equal to 10. For each time step the numerical calculation software stores the solution obtained. In this way, the complete evolution of the temperature map of the plastic part can be evaluated throughout the cooling phase.

•
The physical and thermal properties of the coolant flow along the conformal cooling system are analyzed using the Run 3D cooling channels analysis option • The type of solver used for the numerical analysis of the cooling phase is of the maximum variation of mold temperature type, the temperature difference being equal to 1 • and the maximum number of cycles C being 10.

•
The turbulence model of the numerical software is activated to analyze the contact of the coolant flow with the walls of the conformal cooling channels. The magnitude of the roughness defined on this contact surface is equal to 0.02 mm Also, Table 6 shows the configuration parameters of the numerical simulations and the magnitude of the geometrical and technological variables obtained for each case study analyzed in this paper. As shown in Figure 13, each case study presents boundary conditions of the type inlet and outlet of coolant flow. On the one hand, the coolant flow inlet temperatures are shown in Table 6, and on the other hand the coolant flow inlet velocity is set so that the dimensionless Reynold number is maintained above 4000 throughout the entire system path cooling that is in turbulent regime.
To discretize all the domains shown in Figure 13, three-dimensional elements of the tetrahedron and prismatic type (Solid 186, 10 and 15 nodes, see Figure 4) were used. To define the mesh a dimensioning operation that determines the approximate size of each element of the mesh is used, according to an average length equal to the precision ε of discrete geometry P' ∈ R 3 (see Table 7). Furthermore, in order to improve the quality of the mesh in the two interfaces formed by the contact areas between the plastic part and mold surface and the mold surface and the conformal cooling channels, five layers of SOLID 186 prismatic elements were defined (see Figure 4). The offset ratio used in this paper, defined as the ratio between tetrahedron elements size and prismatic elements size is equal to 0.1. To discretize all the domains shown in Figure 13, three-dimensional elements of the tetrahedron and prismatic type (Solid 186, 10 and 15 nodes, see Figure 4) were used. To define the mesh a dimensioning operation that determines the approximate size of each element of the mesh is used, according to an average length equal to the precision ε of discrete geometry P' ∈ ℝ 3 (see Table 7). Furthermore, in order to improve the quality of the mesh in the two interfaces formed by the contact areas between the plastic part and mold surface and the mold surface and the conformal cooling channels, five layers of SOLID 186 prismatic elements were defined (see Figure 4). The offset ratio used in this paper, defined as the ratio between tetrahedron elements size and prismatic elements size is equal to 0.1.

Implementation and Results
To demonstrate the benefits of the proposed methodology, it was applied to four industrial parts with complex shapes, manufactured by means of plastic injection technology. The accuracy of the nodal mesh (ε) for each geometry is set according to the dimension of the smallest detail of the part. The algorithms were developed in the software Moldex3D R16 [45] (to perform the cooling numerical simulations) and Matlab R2013a [49], Catia V5-6R2016 [50] (to define the geometries of the cooling channels) and with an MSI notebook with an Intel (R) Core (TM) i-77700HQ CPU @ 2.80 GHz. The 3D modeling of the conformal cooling channels was automated by creating an application in the programming environment of the CAD software Catia V5-6R2016 [50].
The proposed methodology establishes the automated design of conformal cooling, based firstly on the topology and geometrical features of the plastic part and secondly on the temperature map during the injection phase, without applying post-processing operations. The location of the cooling elements and their technological parameters are automatically defined. To validate the proposed methodology, several numerical simulations were performed using the injection molding simulation software Moldex3D R16 [45]. The proposed methodology provides a discrete multidimensional model that includes the X, Y and Z coordinates of the axes of the conformal cooling channels, and a set of geometrical and technological variables obtained from the genetic optimization algorithm in order to dimension the cooling system. All this information in discrete format is imported into the software Catia V5-6R2016 [50] to generate the CAD geometry of the cooling system.
Cycolac MG47 (Acrylonitrile butadiene styrene, ABS), Lexan HF1110R (Polycarbonate, PC) and The final design of the conformal cooling system for the mold that will manufacture the four plastic parts under study is shown graphically in Case A, B, C, D. Table 8 shows the magnitude of the geometrical features required for the analysis of the injection process of the plastic parts. To validate the proposed methodology, the variables resulting from the numerical analyses were coolant temperature, pressure map of conformal cooling channels, surface temperature of the part after the cooling phase and surface temperature gradient of the injection mold. The results of the geometrical and technological variables obtained with the genetic optimization algorithm are included for each plastic part analyzed in this manuscript in Tables 9-12. Tables 9-12 include the cooling time and the ejection time of the numerical results established for each plastic part analyzed. The cooling time and ejection time value is defined as the time after which the temperature of the mold surface is less than 10 • C. This requirement is established by the industrial sector in the injection molding area as a criterion to guarantee an optimum level of quality in the injected plastic part. As can be seen in Tables 9-12, the cooling time obtained for the numerical analyses is close to the analytical and ideal cooling time. This value is determined from the geometrical features of the plastic part and the thermoplastic material. The relative error for the plastic geometries analyzed is: 3.118%, 5.000%, 3.225% and 3.776%, respectively. In this way, it is possible to validate the thermal efficiency of the proposed conformal cooling system based on the results obtained against the ideal cooling solution of the numerical software, resulting in a minimum relative error between both solutions. Conventional cooling systems have large differences with respect to analytical or ideal cooling times, particularly in plastic parts with deep concavities, slender details, internal bosses, etc.

Case A
In this section the results obtained for Case A are presented in Figures 14-19.
Polymers 2020, 12, 154 21 of 34 resulting in a minimum relative error between both solutions. Conventional cooling systems have large differences with respect to analytical or ideal cooling times, particularly in plastic parts with deep concavities, slender details, internal bosses, etc.

Case A
In this section the results obtained for Case A are presented in Figures 14-19.    resulting in a minimum relative error between both solutions. Conventional cooling systems have large differences with respect to analytical or ideal cooling times, particularly in plastic parts with deep concavities, slender details, internal bosses, etc.

Case A
In this section the results obtained for Case A are presented in Figures 14-19.    resulting in a minimum relative error between both solutions. Conventional cooling systems have large differences with respect to analytical or ideal cooling times, particularly in plastic parts with deep concavities, slender details, internal bosses, etc.

Case A
In this section the results obtained for Case A are presented in Figures 14-19.

Case B
In this section the results obtained for Case B are presented in Figures 20-25.

Case B
In this section the results obtained for Case B are presented in Figures 20-25.

Case B
In this section the results obtained for Case B are presented in Figures 20-25.

Case B
In this section the results obtained for Case B are presented in Figures 20-25.

Case C
In this section the results obtained for Case C are presented in Figures 26-31.

Case C
In this section the results obtained for Case C are presented in Figures 26-31.

Case C
In this section the results obtained for Case C are presented in Figures 26-31.

Case D
In this section the results obtained for Case D are presented in Figures 32-37.

Case D
In this section the results obtained for Case D are presented in Figures 32-37.

Case D
In this section the results obtained for Case D are presented in Figures 32-37.

Conclusions
The methodology presented in the present manuscript develops a new algorithm for the automated designing and dimensioning of the main elements of a conformal cooling system in injection molding established via the plastic part virtual mesh model and the map of temperatures of the thermoplastic filling flow during the filling phase. The algorithm groups and classifies the discrete temperature of the nodes at the end of the filling phase into geometrical areas called temperature clusters. The topological and rheological information of the clusters along with the geometrical information of the surface mesh is stored in a multidimensional discrete model of the

Conclusions
The methodology presented in the present manuscript develops a new algorithm for the automated designing and dimensioning of the main elements of a conformal cooling system in injection molding established via the plastic part virtual mesh model and the map of temperatures of the thermoplastic filling flow during the filling phase. The algorithm groups and classifies the discrete temperature of the nodes at the end of the filling phase into geometrical areas called temperature clusters. The topological and rheological information of the clusters along with the geometrical information of the surface mesh is stored in a multidimensional discrete model of the

Conclusions
The methodology presented in the present manuscript develops a new algorithm for the automated designing and dimensioning of the main elements of a conformal cooling system in injection molding established via the plastic part virtual mesh model and the map of temperatures of the thermoplastic filling flow during the filling phase. The algorithm groups and classifies the discrete temperature of the nodes at the end of the filling phase into geometrical areas called temperature clusters. The topological and rheological information of the clusters along with the geometrical information of the surface mesh is stored in a multidimensional discrete model of the

Conclusions
The methodology presented in the present manuscript develops a new algorithm for the automated designing and dimensioning of the main elements of a conformal cooling system in injection molding established via the plastic part virtual mesh model and the map of temperatures of the thermoplastic filling flow during the filling phase. The algorithm groups and classifies the discrete temperature of the nodes at the end of the filling phase into geometrical areas called temperature clusters. The topological and rheological information of the clusters along with the geometrical information of the surface mesh is stored in a multidimensional discrete model of the plastic part. The methodology was evaluated by means of its implementation on four plastic parts with diverse topology and measurements. The main geometrical and technological parameters of the conformal cooling system, such as mold surface temperature, coolant inlet temperature, cooling time, distance between conformal cooling channels and plastic part, separation between conformal cooling channels, and diameters of the conformal cooling channels were dimensioned by means of expert algorithms. The layout of the conformal cooling channels for the four cases studied were obtained in an automatic way as well as detailed data reports on the results of the application. These reports included information related to conformal channel diameters D conformal (mm), coordinates X, Y and Z of each conformal channel pathings (mm), the temperature of the injection mold surface T mold ( • C), coolant inlet temperature T coolant ( • C) and cooling time t cooling (s).
The method presented above improves on any conventional cooling system design model since the cooling times obtained are analogous to the cooling times of analytical models including boundary conditions and ideal solutions not exceeding 5% relative error in the cases analyzed. The final quality of the plastic parts after the cooling phase meets the minimum criteria and requirements established by the injection molding companies, obtaining surface temperature differences after the cooling phase of less than 10 • C in all parts analyzed. This new methodology enhances the techniques defined until now as it does not demand heuristic methods to reach the result and does not require the feature recognition method for identifying the plastic part geometry, preventing the problems related with the CAD software and the recognition of complex features. The algorithm is appropriate for every geometry since it operates regardless of the CAD software in which the geometry was modeled. Additionally, the new algorithm applies the utility of the additive manufacturing process for rapid tooling, improving the cooling in remote zones and removing malfunction characteristics, such as warping and residual stress, thus optimizing the cycle time. The proposed algorithm admits the validation and dimensioning of the injection mold cooling system automatically and without demanding high manual computing skills and experience on the part of the injection mold designer.  The theoretical model for the analytical study of the thermal exchange between the plastic part and the coolant flow which runs across the conformal cooling channels is centered on the analysis of the cooling unity. As shown in Figure A1, the cooling cell consists of a rectangular prism that contains the cross-section of the conformal cooling channel as well as the geometrical area of the corresponding plastic part. Each cooling cell is geometrically parameterized according to these variables: diameter of the conformal cooling channels D conformal (m), separation between conformal cooling channels L conformal (m) and the separation between the plastic part and the conformal cooling channels L p-c (m) (see Figure A1 and Table A1). On the other hand, the technological variables that determine the thermal performance of the coolant flow along the conformal cooling channels are: the surface temperature of the mold cavity T mold ( • C) and the coolant flow temperature in the conformal cooling channels T coolant ( • C). The present theoretical model defines the main criteria for the automated calculation of the technological and geometrical variables of the conformal cooling channels.

•
The mechanical and thermal characteristics of the injection mold metallic material where the conformal cooling channels are located are constant and isotropic.

•
The thermal exchange of the thermoplastic filling flow is treated as unidirectional because the mold thickness is lesser than the dimensions related with XY plane.

•
The thermal exchange by convection and radiation between the main elements of the injection mold and the external conditions is discounted because it corresponds to less than 5% of the total heat flow exchanged. That is to say, the thermal analysis between the plastic part and the coolant flow focuses on the thermal exchange by conduction.

•
It is assumed that the coolant flow develops in a turbulent manner (Re > 4000) along the conformal cooling channels According to Xu et al. [28], the steady temperature of the cavity mold surface, defined as Tmold (°C), can be formulated as Equation (A1). The set of geometric and technological variables used in Equation (A1) is shown in Tables A1 and A2.

•
The mechanical and thermal characteristics of the injection mold metallic material where the conformal cooling channels are located are constant and isotropic.

•
The thermal exchange of the thermoplastic filling flow is treated as unidirectional because the mold thickness is lesser than the dimensions related with XY plane.

•
The thermal exchange by convection and radiation between the main elements of the injection mold and the external conditions is discounted because it corresponds to less than 5% of the total heat flow exchanged. That is to say, the thermal analysis between the plastic part and the coolant flow focuses on the thermal exchange by conduction.

•
It is assumed that the coolant flow develops in a turbulent manner (Re > 4000) along the conformal cooling channels According to Xu et al. [28], the steady temperature of the cavity mold surface, defined as T mold ( • C), can be formulated as Equation (A1). The set of geometric and technological variables used in Equation (A1) is shown in Tables A1 and A2. T mold = T coolant + ρ p · C p · T p · (2 · δ s · L conformal + h · π · D conformal · L p−c ) · (T melt − T eject ) 2 · h · π · D cooling · δs · t cooling (A1) Table A2. Technological variables for the thermal exchange between the thermoplastic material and the injection mold. In this way, according to Equations (A1) and (A4), it is known that the equation that establishes the cooling time of the manufacturing process is related to the temperature of the mold cavity surface and vice versa. Therefore, the paper proposes genetic optimization algorithms in order to determine an optimal solution to the physical model as well as to measure the geometrical and technological variables (see Table A1) of the conformal cooling system. It should be noted that the methodology for dimensioning the geometrical and technological variables of the conformal cooling system is applied according to the geometrical discretization or clusters of the plastic part, based on the temperature of the melt plastic flow TMF ij ∈ R for each node P ij ∈ P' n , i.e., since each geometrical discretization or cluster C k (see Figure 7) is associated with an average temperature of the thermoplastic filling flow T k , T k is used as an input to the physical model for dimensioning the conformal cooling system. In this way, the physical model presented is applied to each geometrical discretization or cluster C k in order to obtain the magnitude of the geometrical and technological variables of the conformal cooling system associated with them. Figure A2. Algorithm for the design and dimensioning of the injection mold conformal cooling system.