# Prediction and Sensitivity Analysis of Bubble Dissolution Time in 3D Selective Laser Sintering Using Ensemble Decision Trees

Hai-Bang Ly, Eric Monteiro, Tien-Thinh Le, Vuong Minh Le, Morgan Dal, Gilles Régnier, Binh Thai Pham

► To cite this version:

Hai-Bang Ly, Eric Monteiro, Tien-Thinh Le, Vuong Minh Le, Morgan Dal, et al.. Prediction and Sensitivity Analysis of Bubble Dissolution Time in 3D Selective Laser Sintering Using Ensemble Decision Trees. Materials, 2019, 11 (3-4), pp.1544. 10.3390/ma12091544. hal-02164646v2

HAL Id: hal-02164646

https://hal.science/hal-02164646v2

Submitted on 28 Jun 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Prediction and Sensitivity Analysis of Bubble Dissolution Time in 3D Selective Laser Sintering Using Ensemble Decision Trees

Hai-Bang Ly <fnref n="1" />,*, Eric Monteiro <fnref n="2" />, Tien-Thinh Le <fnref n="3" />,*, Vuong Minh Le <fnref n="4" />, Morgan Dal <fnref n="2" />, Gilles Regnier <fnref n="2" /> and Binh Thai Pham <fnref n="1" />,\*

Received: 31 March 2019; Accepted: 6 May 2019; Published: 10 May 2019

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** The presence of defects like gas bubble in fabricated parts is inherent in the selective laser sintering process and the prediction of bubble shrinkage dynamics is crucial. In this paper, two artificial intelligence (AI) models based on Decision Trees algorithm were constructed in order to predict bubble dissolution time, namely the Ensemble Bagged Trees (EDT Bagged) and Ensemble Boosted Trees (EDT Boosted). A metadata including 68644 data were generated with the help of our previously developed numerical tool. The AI models used the initial bubble size, external domain size, diffusion coefficient, surface tension, viscosity, initial concentration, and chamber pressure as input parameters, whereas bubble dissolution time was considered as output variable. Evaluation of the models' performance was achieved by criteria such as Mean Absolute Error (MAE), Root Mean Squared Error (RMSE) and coefficient of determination (R<fnref n="2" />). The results showed that EDT Bagged outperformed EDT Boosted. Sensitivity analysis was then conducted thanks to the Monte Carlo approach and it was found that three most important inputs for the problem were the diffusion coefficient, initial concentration, and bubble initial size. This study might help in quick prediction of bubble dissolution time to improve the production quality from industry.

**Keywords:** 3D selective laser sintering; artificial intelligence; decision trees; bubble dissolution time; sensitivity analysis

## 1. Introduction

Selective laser sintering (SLS) is one of the most important 3D printing technique widely used in industry [1]. The main idea of this method is that powders are sintered by laser in order to produce solid objects predefined by 3D sketch models [2]. The addictive manufacturing technique exhibits many interesting behaviors not found in traditional methods. For instance, using 3D sketches combined with laser sources, SLS technology could easily handle complex geometries of objects [3–7]. Secondly, laser can melt a lot of materials such as metal, ceramic, polymer or glass powders, especially in case of polymeric particles [8–11]. In industry of polymer processing, the SLS method exposure also has great advantages, particularly in terms of saving time in fabrication [12]. However, combining particles and laser technique requires broadly information involving many multiphysical phenomena. For instance, depend on the nature of polymeric material, relationship between many parameters such as particle

<footnote n="1">¹ University of Transport Technology, Hanoi 100000, Vietnam</footnote>

<footnote n="2">² PIMM, ENSAM, CNRS, CNAM, HESAM Université, 151 Boulevard de l'Hôpital, 75013 Paris, France; Eric.MONTEIRO@ensam.eu (E.M.); Morgan.DAL@ensam.eu (M.D.); gilles.regnier@ensam.eu (G.R.)</footnote>

<footnote n="3">³ Institute of Research and Development, Duy Tan University, Da Nang 550000, Vietnam</footnote>

<footnote n="4">⁴ NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam; vuongminhle09@gmail.com</footnote>

<footnote>* Correspondence: banglh@utt.edu.vn (H.-B.L.); tienthinhle.vn@gmail.com (T.-T.L.); binhpt@utt.edu.vn (B.T.P.)</footnote>
granulometry and morphology, laser power, and scan speed should be established in order to improve
the production quality. Unfortunately, solid structures created using SLS technique are under the
presence of many imperfections such as porosities, cracks, or lacks of fusion [12,13]. The existence of
these defects could be explained by the complexity and the coupling of the different physics occurring
in the process, for instance melting, coalescence, gas diffusion, and crystallization [14–16].

One of the most important defects using SLS technology is the formation of gas bubbles in fabricated parts [8]. As an illustration, for a considered layer during the process, laser beam irradiates polymer particles spread over the previous layer. Gas bubbles are therefore trapped in the melted layer by the difference between resulting temperatures [1]. The larger the initial size of the bubbles, the more dissolution time requires. That means the understanding of dissolution time of gas bubbles has a crucial impact on the processing of solids using SLS technique. However, it is a complex problem which depends on many factors, such as bubble internal radius, external radius, diffusion coefficient, surface tension, viscosity, saturation, and pressure [17]. Until now, no known analytical work has been proposed to predict such bubble dissolution time. Only semi-analytical or numerical methods have been established for investigating the influence of parameters on dissolution time. For instance Kontopoulou et al. [8] studied the impact of viscosity, initial bubble size, and applied pressure in order to determine the kinetic of bubble shrinkage. In another work, surface tension, diffusion coefficient, and saturation of polymer melt have been integrated to establish an equation of bubble dissolution time [18]. Many other attempts have also been introduced in order to relate the dissolution of bubbles and physical constant of the problem [19–21].

In general, there are three techniques commonly used in the literature to predict the mechanical
properties of materials: computational modeling and simulation, experimental studies, or artificial
intelligence (AI) approach. The modeling and simulation of physical-mechanical phenomena is usually
complicated, requires extensive computation time, and is thus not widely used in practice. Experimental
studies are costly, time consuming, and required a lot of specific equipment. Out of these, the artificial
intelligence techniques have been recently used in many fields of science [22–24]. This technique exhibits
a lot of advantages not found in traditional computational simulation [25], especially for predicting
material behaviors, for instance, rheological properties of materials [26], permeability of polymers [27],
composite materials [28], structure under compression [29], or concrete [30]. Nevertheless, to the best
of our knowledge, no artificial intelligence approaches have been applied to predict bubble dissolution
time in the SLS process as well as quantify the relationship between physical parameters of the problem.
Moreover, works related to the modeling of bubble dissolution process only consider an infinite melted
medium [8,18], which is not the case in SLS technique due to the entrapped gas bubble are close to the
surface of powder bed. A proper approach is needed to take into account the finite medium in SLS
process. In addition, laboratory experiments could only be used to observe the growth and collapse
phenomena [31] but not the factors that affect the mechanism. The approach using AI algorithms, once
developed, could provide efficient information and help to improve the quality of SLS production.

This contribution applied two AI algorithms namely the Ensemble Bagged Trees (EDT Bagged) and Ensemble Boosted Trees (EDT Boosted) to predict bubble dissolution time. The outline of the present study is organized as follows. The Introduction (Section 1) starts with the state-of-the-art of the problem. Section 2 (Materials and Methods) is devoted to describing the mechanism of the bubble shrinkage problem and the construction of data used for AI modeling. Bubble initial radius, external domain radius, diffusion coefficient, surface tension, viscosity, saturation, and chamber pressure were considered as inputs; and bubble dissolution time is output for the prediction tools. A brief introduction of two AI techniques based on the Ensemble Decision Trees algorithms is presented next, together with the Monte Carlo approach for the input sensitivity analysis. Principal results and discussions are presented in Section 3 (Results and Discussions), including the prediction capability of each model along with the sensitivity analysis of input parameters. Finally, the paper is completed with several conclusions and perspective in Section 4 (Conclusions).
## 2. Materials and Methods

### 2.1. Physical Problem and Data Preparation

#### 2.1.1. Mechanical Description of the Bubble Dissolution Process

When fabricating parts in the SLS process, grains coalesce faster at the surface than deeper in the powder bed because of higher laser radiation absorption at its surface. Air porosities are then embedded in the melted layers, whose thickness can reach a few hundred of microns. Air porosity, or air bubbles, is an essential factor accounting for the final mechanical properties of the fabricated parts. Thus, understanding bubble dissolution process in the melted polymer is important, especially in SLS technique.

Bubble dissolution mechanism is a highly coupled and nonlinear physical problem. Noting that proper development of equations has been presented in our previous work, in this paper, only principal equations that have impact on the AI modeling part will be briefly recalled. The nature of the resorption process lied on the diffusion of gas inside the bubbles (the gas phase) into the melted materials, or the liquid phase, such as metal, glass, or polymer. During the resorption process, the movement of the interface between gas/liquid phases obeyed the fluid dynamics equations, represented by the equation of conservation of momentum [17]:

$$[P_b(t) - P_c(t)]R_i(t) - 2\gamma\left(1 + \frac{R_i(t)}{R_e(t)}\right) + 4\eta\left(\frac{R_i^3(t)}{R_e^3(t)} - 1\right)\dot{R}_i(t) = 0 \quad (1)$$

where $P_b(t)$ and $P_c(t)$ represent the pressure of the perfect gas inside the bubble and the chamber of the SLS machine, respectively; $R_i(t)$, $R_e(t)$ denotes the bubble initial radius and the radius of the outer liquid phase; $\dot{R}_i(t)$ is the rate of bubble dissolution; $\eta$ refers to the fluid viscosity and $\gamma$ stands for the surface tension. The liquid phase is characterized by fluid viscosity, whereas the surface tension depends on both media as it is resulted at the air-liquid interface. Besides, the rate of gas diffusion is related to the movement of the interface through the diffusion equation, or the Fick law [17]:

$$\frac{\partial c}{\partial t} = \frac{D}{r^2} \frac{\partial}{\partial r} (r^2 \frac{\partial c}{\partial r}) \quad (2)$$

where $D$ stands for the diffusion coefficient, $c$ is the concentration of gas in the fluid phase. The degree of saturation is introduced in order to take into account the boundary conditions derived from Henry's law such as [17]:

$$c(r, 0) = x_{sat} K_H P_C(0) \quad (3)$$

with $K_H$ the Henry's law constant. Finally, the mass balance equation at the bubble interface, in considering the decrease of perfect gas within the bubble is proportional to the flux across the interface, reads [17]:

$$\frac{\partial}{\partial t}\left(\frac{P_b(t)R_i^3(t)}{R_g T}\right) = 3R_i^2(t)D\left.\frac{\partial c(r,t)}{\partial r}\right|_{r=R_i} \quad (4)$$

with $R_g$ the universal gas constant and $T$ the temperature. For summary, the dissolution mechanism depends on the nature of the fluids (the molten materials), the mechanism of gas diffusion (air) and the mass balance at the gas/fluid interface. In the aim of predicting bubble dissolution time in the selective laser sintering (SLS) process, the need of experimental values of physical constants (inputs) as well as results (targets) are important. However, to the best of our knowledge, very few existing investigations clearly determining such parameters have been reported so far. These values are extremely difficult to determine. Therefore, simulation data constructed from previously developed numerical tool [17] have been used in this investigation. In this work, a finite element-based numerical tool has been developed in order to model and explain the influence of diverse parameters on the dissolution time, such as
initial bubble size ($R_i$), the external radius of the fluid phase ($R_e$), the surface tension ($\gamma$), viscosity ($\eta$), diffusion coefficient ($D$), initial melt concentration ($x_{sat}$), and the chamber pressure ($P$). The fluid domain was assumed Newtonian. It is noteworthy that results extracted from such a numerical tool have already been compared and validated with experimental studies on the bubble removal in rotational molding process [8,18], a technique similar to the SLS.

### 2.1.2. Factors that Affecting Bubble Dissolution Time

In this paper, the prediction capability and sensitivity analysis of physical parameters that may have influence on the dissolution process were performed. It is important noticed that a change of any parameter affects all others, which is governed by the dynamic, thermodynamics, or physics of the bubble shrinkage problem [17]. For example, the change of gas inside bubble or the polymer melt causes changes in viscosity, surface tension, and diffusion coefficient applied to both liquid and gas phases. Therefore, all the parameters were varied in a broad range in order to fully collect results and analyze their effect to the bubble dissolution time.

Input parameters considered for the AI simulation were bubble radius ($I_1$), bubble external radius ($I_2$), diffusion coefficient ($I_3$), surface tension ($I_4$), viscosity of the melted materials ($I_5$), initial saturation ($I_6$) and chamber pressure ($I_7$). The prediction output was bubble dissolution time (O). The values of physical parameters, i.e., $I_1$, $I_2$, $I_5$, (used in Equation (1)) were taken based on experimental data in the work of Kontopoulou et al. [8]. The diffusion coefficient ($I_3$) was varied around that in the case of air diffuse into polyethylene melt (used in Equations (2) and (4)), which was estimated previously in the work of Griskey [9]. Surface tension values ($I_4$) had been taken based on the case of polyethylene melts (used in Equation (1)), which was gathered from earlier study of Wu [32]. The viscosity of the melted polymer ($I_5$—used in Equation (1)) was collected following the works of Bird et al. [33] and Durrill et al. [34]. Being used in Equation (3), the initial saturation ($I_6$) and the chamber pressure ($I_7$) were varied from that of the analytical study of Gogos [18]. Parameters such as the Henry's law constant, air parameters or the temperature dependence needed for applying an Arrhenius equation were not considered as inputs, as they are often directly related to one or several inputs $I_1$ to $I_7$. Simulated data were taken to cover all the possibilities in the SLS process. A total number of 68,644 data were generated, lied on the fact that the corresponding dissolution time of these data varied from 20 s to 1200 s, which is of interest in the SLS process. If it is lower than 20 s, the dissolution process is quite quick and improvement of the SLS process appears less important. If it is higher than 1200 s, the dissolution process is too long so that it appears not interesting to be studied. A summary of the input and outputs values are presented in Table 1.

<table><caption>**Table 1.** Input and output quantities used in the artificial intelligence (AI) models.</caption><thead><tr><td>Input Parameter</td><td>Minimum</td><td>Maximum</td><td>Average</td><td>Standard Deviation</td></tr></thead><tbody><tr><td>I<sub>1</sub> (bubble radius, µm)</td><td>30</td><td>150</td><td>112.4</td><td>19.37</td></tr><tr><td>I<sub>2</sub> (external radius, µm)</td><td>33</td><td>150,000</td><td>333.1</td><td>422.15</td></tr><tr><td>I<sub>3</sub> (diffusion coefficient, m<sup>2</sup>·s<sup>-1</sup>)</td><td>0.1 × 10<sup>-9</sup></td><td>50 × 10<sup>-9</sup></td><td>4.1 × 10<sup>-9</sup></td><td>6.79 × 10<sup>-9</sup></td></tr><tr><td>I<sub>4</sub> (surface tension, N·m<sup>-1</sup>)</td><td>0.0100</td><td>0.0500</td><td>0.0288</td><td>0.0144</td></tr><tr><td>I<sub>5</sub> (viscosity, Pa·s)</td><td>100</td><td>10,000</td><td>3670.5</td><td>3701.3</td></tr><tr><td>I<sub>6</sub> (saturation)</td><td>0.1</td><td>1</td><td>0.54</td><td>0.37</td></tr><tr><td>I<sub>7</sub> (chamber pressure, atm)</td><td>0.25</td><td>1.5</td><td>0.88</td><td>0.41</td></tr><tr><td>O (bubble dissolution time, s)</td><td>20</td><td>1200</td><td>194.2</td><td>234.9</td></tr></tbody></table>

As shown in Table 1, the dissolution time of bubble considered in this study varied from 20 to 1200 s, with a mean value of 194.2 s and a standard deviation of 234.9 s. That means the distribution of bubble dissolution time was highly asymmetric and exhibited a big amount of quantity distributed between 20 s and 200 s. This interval of dissolution time was very important in the SLS processing industry. The lower dissolution time of bubbles, the higher production quality.
### 2.1.3. Datasets

In general, the training and testing datasets are required for constructing and validation of AI models. Data were divided into 2 groups at a ratio of 70/30. A total number of 48,051 training data (70%) were used to construct the prediction models, whereas 20,593 testing data (30%) were used for evaluation and validation of the models. The 70/30 ratio for generating the datasets in this study has been done as suggested by Khorsheed and Mohammad [35], and Leema et al. [36]. Description of two AI algorithms in this study will be briefly presented in the next section.

## 2.2. Background of Models Used

### 2.2.1. Decision Trees Methods

Decision trees (DT), proposed by Quinlan [37], is a popular machine learning method that can be applied in solving many real world problems such as identifying risk factors for drug use [38], flash flood prediction [39], landslide spatial prediction [40], prediction of a short-term photovoltaic power [41]. Main principle of the DT is to use a series of rules to identify the regions with the most homogeneous output variables to input variables on which a constant is fitted to each region. Main advantages of the DT are (i) it represents information in an intuitive and easy way for visualization, (ii) it is powerful for mining nonlinear and interactions effects between dependent and independent variables, (iii) it requires no mathematical assumptions between output and input variables, and (iv) it is capable to handle missing values and outliers [42]. On the other hands, the DT also faces several problems such as (i) it is difficult to model smooth functions, (ii) the structure of tree is sensitive with the sample data as only small change in training data can give very different results, and (iii) it has high variance and low bias [43]. Therefore, many techniques have proposed to enhance the predictive capability of the DT such as Ensemble Bagged Trees and Ensemble Boosted Trees.

#### Ensemble Bagged Trees (EDT Bagged)

EDT Bagged is a hybrid model which is a combination of bagging algorithm and decision trees [44]. Out of these methods, bagging method uses the bootstrap sampling to optimize the input training data for learning the decision trees [45]. This process can be carried out in several steps such as (i) it takes a bootstrap sample from the original training dataset, (ii) it fits the decision trees with the extracted dataset, (iii) it then finds the best fit models with the optimized dataset, (iv) it thereafter uses each of the fitted models with the optimized dataset to predict the results, and (v) it finally averages the prediction results of all the models [46]. EDT Bagged can significantly improve the prediction accuracy as bootstrap aggregation used in bagging method can reduced the variance in individual prediction methods like decision trees. EDT Bagged has been applied to solve many real world problems namely travel time prediction [44], disease prediction [45], and bankruptcy prediction [47]. In this study, we have used the EDT Bagged for prediction of bubble dissolution time when manufacturing solid structures by SLS technique.

Suppose that $x = (x_1, x_2, \dots, x_n)$ is a vector of input variables ($n$ is the number of input variables), the bootstrap aggregation used in EDT Bagged to average the predictions is defined as follows [46]:

$$f_{bag}(x) = \frac{1}{N} \sum_{i=1}^{N} f_{i}^{*}(x) \qquad (5)$$

where $N$ is the number of various training datasets generated by bootstrapping, and $f_i^*(x)$ are defined as predictors.

#### Ensemble Boosted Trees (EDT Boosted)

EDT Boosted is an additive regression model in which individual terms are simple trees, fitted in a forward, stagewise estimation [48]. It is well-known as an ensemble method which combines
regression trees and boosting techniques in which the regression trees are used to connect output variables to input variables by recursive binary splits whereas the boosting technique is used to combine many single models to improve the predictive capability of the ensemble models [49]. EDT Boosted takes advantages of tree-based methods, and even overcomes the drawback of single tree models as (i) it is able to select relevant variables to fit accurate functions, (ii) it can be fitted with different amounts of data using stochastic boosting, and (ii) it can reduce both bias by forward stagewise fitting and variance by model averaging [48]. EDT Boosted has been applied to many fields such as banking [49], ecology [50], and medical [51]. In this study, EDT Boosted is a first time use for predicting the dissolution time of bubble removal when using laser to melt polymer powders in a SLS procedure.

Suppose that $x = (x_1, x_2, ..., x_n)$ is a vector of predictors ($n$ is the number of predictors) and $y$ is a response, the EDT Boosted algorithm is trained by following function [52]:

$$f(x) = \sum_{k} f_k(x) = \sum_{k} \alpha_k t(x, y_k) \qquad (6)$$

where $\alpha_k$ values represent weights given to the nodes of each tree, $y_k$ is defined as the split variables, $t(x, y_k)$ represents single decision trees.

### 2.2.2. Monte Carlo Method

In many fields of science, average values of parameters are not sufficient enough for investigating the problem and also the complex relationship between parameters. For instance, material variability and structural variability could largely impact the effective mechanical behaviors of composite structure [53–56]. Along with the development of computational algorithms and computer performance, Monte Carlo method appeared then as a highly efficient technique for accounting variability of parameters when they become relevant [57]. This method could help researchers to determine which parameters are important on the overall responses following variability of each input. Indeed, Monte Carlo method exhibits high capability to quantify numerically linear or non-linear relationship in data. The main idea of this method is to repeat random sampling of inputs using probability density distribution of each. By doing so, variability and of inputs could be fully propagated to the overall properties though constructed models. Statistical analysis of outputs is next carried out in order to quantify the impact of input variability on the responses. A schematization of using Monte Carlo method for propagating input variability is presented in Figure 1.

<figure><img src="image_3.png" /><figcaption>Figure 1. Schematization of using Monte Carlo method for the propagation of input variability.</figcaption></figure>
Monte Carlo method could also be used for analyzing the sensitivity of input parameters on the responses through constructed models [58]. For instance, the considered input has no influence on the final responses if it is removed from the input space and the final responses has statistically no differences compared to the case of simulation with all inputs. This way, the sensitivity of all inputs to the overall properties could be entirely quantified. This analysis has an important role for designer, for instance, dimension of the problem could be reduced and time consuming could be saved by eliminating insignificant parameters. To the contrary, sensitive parameters exposure important role and receive so huge attention.

In this study, the “normdiff” function of random variable P was introduced in order to: (i) determine the corresponding optimal number of Monte Carlo realizations and (ii) evaluate the degree of fluctuation of P around its average value (in % unit). The normdiff function was defined as [55]:

$$ \text{normdiff}(N_r) = \frac{1}{\bar{P}N_r} \sum_{i=1}^{N_r} P_i \qquad (7) $$

where $N_r$ is the actual number of Monte Carlo runs and $\bar{P}$ is the average value of P. In this work, the normdiff was applied for estimating the statistical convergence of several random variables such as Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and coefficient of determination (R²) respectively (see Section 2.2.3.).

Uniform distribution was applied to perform random sampling of input data for Monte Carlo simulations. Random subsets of input data were drawn uniformly between 1 and last sample index. The rectangular distribution ensures in this case no constraints for generating input realizations (i.e., all sub-spaces of the input space were accounted in the simulation).

### 2.2.3. Validation Criteria

In this study, popular validation criteria namely Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and coefficient of determination (R²) were selected for validating the predictive capability of ensemble decision trees [59,60]. RMSE measures the average squared difference between actual and predicted values, MAE measures the average magnitude of the errors, and R² measures the correlation between actual and predicted values [59,60]. Whereas R² varies from 0 (lowest models) and 1 (best models), higher RMSE and MAE shows lower predictive capability of the models. We can calculate these RMSE and MAE and R² by following equations [61,62]:

$$ RMSE = \sqrt{\frac{1}{m} \sum_{i=1}^{m} (v_i - \bar{v}_i)^2} \qquad (8) $$

$$ MAE = \frac{1}{m} \sum_{i=1}^{m} (v_i - \bar{v}_i) \qquad (9) $$

$$ R^2 = 1 - \frac{\sum_{i=1}^{m} (v_i - \bar{v}_i)^2}{\sum_{i=1}^{m} (v_i - \bar{v}_i)^2} \qquad (10) $$

where $m$ is considered as the number of the samples, $v_i$, $\bar{v}_i$, and $\bar{v}$ are the actual output, the predicted output, and the mean of the $\bar{v}_i$, respectively.

### 2.2.4. Methodology Chart

We present here different steps of the present work’s methodology, which contains five main steps as followings (as shown in Figure 2):
**Step 1.** Numerical tools for collecting data. In this step, all input parameters such as bubble initial radius ($I_1$), bubble external radius or the domain size ($I_2$), the diffusion coefficient ($I_3$), surface tension ($I_4$), viscosity of the melted materials ($I_5$), the initial saturation ($I_6$), and the chamber pressure ($I_7$) were used to compute the prediction output-bubble dissolution time (O), using the developed finite element model [17]. Fluid dynamics, mass balance and gas diffusion were taken into account in this numerical model.

**Step 2.** Dataset preparation for training AI models. In this step, all input and output parameters were used to create a complete set of data. A number of 70% data (48,051 training data) were taken from the initial dataset for training the two proposed AI models. The remaining 30% data (20,593 testing data) were used for validating step, after the AI models were trained.

**Step 3.** Training models. In this step, the AI models were trained using the training data set. Two AI algorithms were used: EDT Bagged and EDT Boosted. The algorithm of these two models were introduced in the previous section. This step was repeated until the models were successfully trained (a tolerance error criteria was reached).

**Step 4.** Validating models. In this step, after the two AI models were trained, the validation process was performed using the testing data set (30% of the initial data set). The models were validated using different statistical measures such as: Root Mean Squared Error (RMSE), Mean Absolute Error (MAE) and the coefficient of determination ($R^2$).

**Step 5.** Sensitivity analysis. After validation of the two proposed AI models using the testing data set, the sensitivity of input parameters was performed. This process was carried out thanks to the Monte Carlo approach.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Methodology of the proposed analysis.</figcaption></figure>

# 3. Results and Discussion

## 3.1. *Comparison of Ensemble Decision Trees (EDT) Algorithms*

Validation of the two proposed AI models are illustrated in Figure 3 for the training part including 48,051 data (left) and testing part with 20,593 data (right). Error frequency are plotted in function of relative errors between target and output values for EDT Bagged (blue line) and EDT Boosted (orange line) algorithms. In both training and testing parts, the EDT Boosted exhibit more error between output and target, as the values ranging from -1 to 6. The predicted outputs are ranging from 1 time lower to 6 times higher than the target ones. As regard to the training part, the over predicted values outside the range -1 to 1 were 1880 samples (or 3.9%) for EDT Boosted and 275 samples (or 0.57%) for EDT Bagged. For the testing part, the over predicted values outside the range -1 to 1 were 893 samples (or 4.3%) for EDT Boosted and 163 samples (or 0.79%) for EDT Bagged, showing reasonable performance
of the AI models. Considering the EDT Bagged algorithm, the error range between output and target is rather narrower than the EDT Boosted model, i.e., ranging from -0.5 to 2. The peaks (blue line) of both training and testing parts are centered at 0, showing excellence prediction capability of EDT Bagged. For both AI algorithms, the number of over predicted values are more than that of under predicted values, confirmed by the error curves are rather on the right side of 0. It is worth noticed that for this simulation, the error measurements are: RMSE = 24.38, MAE = 9.70, R² = 0.989 for training EDT Bagged, RMSE = 28.28, MAE = 9.91, R² = 0.986 for testing EDT Bagged, whereas these values are RMSE = 141.48, MAE = 79.64, R² = 0.852 for training EDT Boosted, RMSE = 144.80, MAE = 81.80, R² = 0.843 for testing EDT Boosted. In general, both of methods show good predictive capability based on the validation criteria R², RMSE and MAE. Comparing the prediction performance using testing parts, EDT Bagged appears a better candidate for this problem, as the RMSE, MAE, and the number of wrong predictions are smaller than those of EDT Boosted, whereas the R² values are higher than that of the EDT Boosted model.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Histogram of relative error between predicted and output values of training (**left**) and testing (**right**) part of Ensemble Bagged Trees (EDT Bagged) and Ensemble Boosted Trees (EDT Boosted) algorithms.</figcaption></figure>

In validation stage, it is difficult to conclude the performance of given AI algorithms unless a fully analysis is performed. As 70% of data are randomly taken to train the AI prediction tools, the corresponding performance measurements (RMSE, MAE, or R²) are different for each simulation [63]. Such values depend of the choice of combination of 70% input data and can be varied from one to another simulation. In order to check the robustness of EDT Bagged and EDT Boosted models, thus, 1000 different arrangements of data were generated using uniform distribution and the error measurements in each case were collected. RMSE, MAE, or R² errors for 1000 different simulations are plotted under scatter points in Figure 4 for training and testing parts. The variation of error measurements in the training part is rather smaller than that of the testing part in all cases, which clearly demonstrate the dependence of the predicted results on input combination. In the case of EDT Boosted model (Figure 4a-c), the training and testing part vary around the same values, i.e., RMSE = 141.21, MAE = 79.65, and R² = 0.855. Regarding EDT Bagged algorithm (Figure 4d-f), the training results are RMSE = 25.82, MAE = 9.83 and R² = 0.988, whereas the testing one are RMSE = 23.78, MAE = 9.14, and R² = 0.990. This indicates that EDT Bagged outperforms EDT Boosted. However, both EDT Boosted and EDT Bagged algorithms possess ability to well-predict bubble dissolution time in the SLS process.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Error distribution (Root Mean Squared Error (RMSE), Mean Absolute Error (MAE), and coefficient of determination (R²)) of different models for 1000 simulations: (a–c) EDT Boosted and (d–f) EDT Bagged algorithms.</figcaption></figure>

The normalized statistical convergence analysis of the obtained results is calculated using Equation (3) and plotted in Figure 5. It is observed that over two Monte Carlo runs for EDT Boosted model are required to achieve the convergence within the 1% range around the converged values. On the contrary, RMSE and MAE criteria require at least 30 simulations to achieve converged values. For R² criterion, both algorithms are converged in a 1% error range without any simulations. These results demonstrate that the two proposed AI algorithms are promising methods for the prediction of bubble dissolution time with an interesting convergence rate. Out of these, the EDT Bagged is better in term of precision than the EDT Boosted but it requires more simulations to achieve converged statistical values.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Statistical convergence analysis of (a) RMSE, (b) MAE, and (c) R² over 1000 Monte Carlo simulations for EDT Bagged and EDT Boosted algorithms. Blue lines representing ±1% deviation around the average value.</figcaption></figure>

## 3.2. Sensitivity Analysis of Input Parameters

Prediction of bubble dissolution time is a complex problem in which physical and mechanical equations are highly coupled and nonlinear. Thus, the sensitivity analysis is then carried out in order to evaluate the impact of each single input parameter to the predicted output. This could be an attempt
to reduce the input space if one input variable is found not affect the final prediction results. All the input parameters were successively excluded from the input space by setting the column to zero values. Monte Carlo simulations were then performed in order to quantify the influence of each parameter, thanks to the prediction performance, i.e., RMSE, MAE, and R². In excluding successively I₁, ..., I₇, a total number of 7 groups of analysis were generated and simulated with a number of 1000 runs for each group. The statistical performances (RMSE, MAE, and R²) of two AI models are highlighted in Table 2. Based on mean values of RMSE, MAE and R² over 14,000 simulations, the EDT Bagged outperforms the EDT Boosted. Moreover, standard deviation values of EDT Bagged model are smaller than that of the EDT Boosted, indicating that the EDT Bagged algorithm is a more stable prediction technique.

<table><caption>**Table 2.** Sensitivity analysis of input parameters with respect to RMSE, MAE, and R² obtained after 1000 Monte Carlo simulations for each case: full simulation with 7 inputs (No excl.) and excluding successively I₁ to I₇ (I₁, ..., I₇ excl.).</caption><thead><tr><th rowspan="2">Criteria</th><th rowspan="2">Input Excl.</th><th colspan="4">EDT Bagged</th><th colspan="4">EDT Boosted</th></tr><tr><th>Min</th><th>Max</th><th>Mean</th><th>SD</th><th>Min</th><th>Max</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td rowspan="7">RMSE</td><td>No excl.</td><td>19.20</td><td>28.27</td><td>23.78</td><td>1.26</td><td>135.86</td><td>146.70</td><td>141.16</td><td>1.79</td></tr><tr><td>I<sub>1</sub> excl.</td><td>97.48</td><td>104.97</td><td>100.80</td><td>1.15</td><td>153.04</td><td>161.99</td><td>157.11</td><td>1.49</td></tr><tr><td>I<sub>2</sub> excl.</td><td>56.13</td><td>63.04</td><td>58.84</td><td>1.02</td><td>138.32</td><td>150.13</td><td>142.84</td><td>1.68</td></tr><tr><td>I<sub>3</sub> excl.</td><td>215.67</td><td>224.28</td><td>220.44</td><td>1.18</td><td>222.84</td><td>223.37</td><td>228.18</td><td>1.61</td></tr><tr><td>I<sub>4</sub> excl.</td><td>40.06</td><td>48.21</td><td>43.71</td><td>1.12</td><td>136.29</td><td>147.91</td><td>141.28</td><td>1.69</td></tr><tr><td>I<sub>5</sub> excl.</td><td>23.50</td><td>36.27</td><td>27.83</td><td>1.70</td><td>136.29</td><td>147.91</td><td>141.28</td><td>1.69</td></tr><tr><td>I<sub>6</sub> excl.</td><td>120.65</td><td>128.45</td><td>124.63</td><td>1.12</td><td>178.98</td><td>189.90</td><td>184.26</td><td>1.67</td></tr><tr><td>I<sub>7</sub> excl.</td><td>33.37</td><td>41.85</td><td>37.24</td><td>1.28</td><td>136.29</td><td>147.91</td><td>141.28</td><td>1.69</td></tr><tr><td rowspan="7">MAE</td><td>No excl.</td><td>7.93</td><td>10.82</td><td>9.14</td><td>0.47</td><td>77.30</td><td>82.33</td><td>79.62</td><td>0.91</td></tr><tr><td>I<sub>1</sub> excl.</td><td>46.87</td><td>51.16</td><td>48.86</td><td>0.69</td><td>85.76</td><td>90.81</td><td>87.85</td><td>0.82</td></tr><tr><td>I<sub>2</sub> excl.</td><td>24.95</td><td>31.25</td><td>26.57</td><td>0.63</td><td>78.39</td><td>83.60</td><td>80.54</td><td>0.85</td></tr><tr><td>I<sub>3</sub> excl.</td><td>163.07</td><td>167.72</td><td>165.46</td><td>0.73</td><td>145.98</td><td>152.37</td><td>149.14</td><td>0.88</td></tr><tr><td>I<sub>4</sub> excl.</td><td>15.13</td><td>20.76</td><td>17.15</td><td>0.75</td><td>77.41</td><td>82.77</td><td>79.65</td><td>0.86</td></tr><tr><td>I<sub>5</sub> excl.</td><td>10.15</td><td>15.99</td><td>12.54</td><td>0.92</td><td>77.41</td><td>82.77</td><td>79.65</td><td>0.86</td></tr><tr><td>I<sub>6</sub> excl.</td><td>69.15</td><td>73.83</td><td>71.64</td><td>0.71</td><td>98.73</td><td>104.90</td><td>101.63</td><td>0.90</td></tr><tr><td>I<sub>7</sub> excl.</td><td>13.24</td><td>19.13</td><td>15.18</td><td>0.83</td><td>77.41</td><td>82.77</td><td>79.65</td><td>0.86</td></tr><tr><td rowspan="7">R<sup>2</sup></td><td>No excl.</td><td>0.986</td><td>0.993</td><td>0.990</td><td>0.001</td><td>0.832</td><td>0.870</td><td>0.856</td><td>0.006</td></tr><tr><td>I<sub>1</sub> excl.</td><td>0.802</td><td>0.827</td><td>0.816</td><td>0.004</td><td>0.721</td><td>0.758</td><td>0.742</td><td>0.005</td></tr><tr><td>I<sub>2</sub> excl.</td><td>0.929</td><td>0.943</td><td>0.938</td><td>0.002</td><td>0.808</td><td>0.858</td><td>0.843</td><td>0.007</td></tr><tr><td>I<sub>3</sub> excl.</td><td>0.100</td><td>0.144</td><td>0.122</td><td>0.006</td><td>0.153</td><td>0.183</td><td>0.168</td><td>0.005</td></tr><tr><td>I<sub>4</sub> excl.</td><td>0.960</td><td>0.971</td><td>0.966</td><td>0.002</td><td>0.831</td><td>0.872</td><td>0.855</td><td>0.006</td></tr><tr><td>I<sub>5</sub> excl.</td><td>0.978</td><td>0.990</td><td>0.987</td><td>0.002</td><td>0.831</td><td>0.872</td><td>0.855</td><td>0.006</td></tr><tr><td>I<sub>6</sub> excl.</td><td>0.704</td><td>0.731</td><td>0.719</td><td>0.004</td><td>0.533</td><td>0.573</td><td>0.554</td><td>0.007</td></tr><tr><td>I<sub>7</sub> excl.</td><td>0.970</td><td>0.980</td><td>0.975</td><td>0.002</td><td>0.831</td><td>0.872</td><td>0.855</td><td>0.006</td></tr></tbody></table>

Let us consider the case where input parameters I₁ to I₇ are successively excluded from the input space. It is observed that when I₃ is excluded from the prediction process, highest mean values of RMSE and MAE are obtained, i.e., RMSE = 220.44, MAE = 165.46 for EDT Bagged algorithm and RMSE = 228.18, MAE = 149.14 for EDT Boosted model. Regarding R², lowest mean values are obtained: R² = 0.122 for EDT Bagged algorithm and R² = 0.168 for EDT Boosted model. This means that without the diffusion coefficient (I₃) in the input space, it is difficult to achieve acceptable performance of the AI prediction tools. Therefore, the diffusion coefficient (I₃) is the most important input variable in the prediction process of bubble dissolution time. Although there is no known analytical solution to predict the bubble dissolution process in Selective Laser Sintering until now, several observations in good agreement with the results can be found in the literature. Kontopoulou et al. [8], Gogos [18] reported the kinetic of bubble shrinkage and concluded that such mechanism is mostly controlled by the diffusion of air from the bubble to the polymer melt. The distribution of the resulted dissolution
time was highly distributed in the range between 20 s and 200 s. This in good agreement with this finding, as the mean value of $I_3$ was rather high, i.e., $4.1 \times 10^{-9} \text{ m}^2\cdot\text{s}^{-1}$.

Based on the RMSE, MAE, and $R^2$ values, the second and the third important input factors in the prediction process are the initial concentration (saturation rate-$I_6$) and the initial bubble size ($I_1$). Indeed, without $I_6$, the as-obtained errors are: RMSE = 124.63, MAE = 71.64, $R^2$ = 0.719 for EDT Bagged algorithm, whereas RMSE = 184.26, MAE = 101.63, $R^2$ = 0.554 in case of EDT Boosted model. Besides, in excluding $I_1$ out of the input space, the corresponding errors are: RMSE = 100.80, MAE = 48.86, $R^2$ = 0.816 for EDT Bagged algorithm, whereas RMSE = 157.11, MAE = 87.85, $R^2$ = 0.742 for EDT Boosted model. The higher values of RMSE, MAE and the lower values of $R^2$ compared to that of the simulation using 7 inputs indicates the necessity of these variables in the prediction of bubble dissolution time which is in good agreement with the work of Kontopoulou et al. [8] where the authors mentioned that the air concentration and bubble initial size were found to be of great importance in the bubble shrinkage process. Again, the air initial concentration distribution was taken to achieve mean value of $I_6$ = 0.54, SD = 0.37, and the initial bubble size mean value $I_1$ = 112.4 µm, SD = 19.37 µm. The skewness in the input space affected the bubble dissolution time to be in a range between 20 s–200 s, which was in good agreement with observations in the literature.

In our previous work [17], the contribution lies on the simulation tool for the prediction of bubble dissolution process. The polymer domain is generally considered an infinite domain in the literature, whereas we have introduced the parameter that accounting for the size of the polymer domain ($I_2$). In the sensitivity analysis of this study, such parameter is found to be more important than the others ($I_4$, $I_5$ and $I_7$). The error when excluding $I_2$ are RMSE = 58.84, MAE = 26.57, $R^2$ = 0.938 for EDT Bagged algorithm and RMSE = 142.84, MAE = 80.54, $R^2$ = 0.843 for EDT Boosted model. The surface tension ($I_4$), applied pressure ($I_7$), and viscosity ($I_5$) are found to be the most insensitive input parameters for the prediction problem (Table 2). The role of these parameters has also been discussed in several works relating the SLS process, such as the surface tension was neglected to derive an analytical solution [64], or its effect is negligible under certain condition, i.e., small degree of saturation [18]. The bubble shrinkage process is not influenced considerably by melt viscosity, as reported in the work of Kontopoulou et al. [8]. Last but not least, the effect of pressure is important only in case when it is applied after the formation of the bubble [8,18]. In our developed numerical tool, the pressure is considered constant during the whole SLS process [17]. Therefore, no conclusion on the effect of pressure can be deduced at this moment.

The histogram of RMSE, MAE, and $R^2$ values of 14,000 simulations are plotted in Figure 6. All the peaks are rather narrow, demonstrating small values of standard deviation. The three most influenced parameters are easily observed, as they are separated to the rest. However, depending on the AI models, other input parameters can be detected or not. In case of EDT Boosted algorithm, the curves relating $I_2$, $I_4$, $I_5$, and $I_7$ are superimposed for RMSE, MAE and $R^2$ criteria. Regarding the EDT Bagged algorithm, these input parameters can be slightly identified by RMSE and $R^2$ criteria, whereas it is difficult for MAE criterion. As a conclusion, sensitivity analysis requires not only adapted error criteria but also a good prediction model that can capture the difference influence between variables.

The normalized statistical convergence analysis of the results are plotted in Figure 7 (see also Equation (3)), for the cases of three most influenced input parameters. Again, it is observed that over two Monte Carlo simulations for both EDT Boosted and EDT Bagged models are required to achieve the convergence within the 1% range around the converged values. On the contrary, for $R^2$ criterion, both algorithms are converged in a 1% error range within 10 Monte Carlo simulations. These results demonstrate that the two proposed AI algorithms are potential candidates for the prediction of bubble dissolution time with an interesting convergence rate. Out of these, the EDT Bagged is better in terms of precision than the EDT Boosted but it requires more simulations to achieve converged statistical values.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Histogram of RMSE, MAE, and $R^2$ for 14,000 simulations in case of: (a–c) EDT Boosted model, and (d–f) EDT Bagged algorithm.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Statistical convergence analysis of RMSE, MAE, and $R^2$ in case of: (a–c) diffusion coefficient excluded, (d–f) initial saturation excluded, and (g–i) initial bubble size excluded. Blue lines representing ±1% deviation around the average value.</figcaption></figure>

4. Conclusions

In this study, two ensemble AI models namely EDT Bagged and EDT Boosted were proposed and compared for predicting bubble dissolution time of the SLS process. A metadata of 68,644
simulated ones were created with the help of our previously developed numerical tool. They were used for generating datasets, which included input parameters (initial bubble size, external domain size, diffusion coefficient, surface tension, viscosity, initial concentration, the SLS chamber pressure) and output parameter (bubble dissolution time). Validation of the models was achieved using several criteria such as MAE, RMSE, and R². Monte Carlo approach was finally used to carry out the sensitivity analysis for determination of the most important input parameter for the prediction of bubble dissolution time.

The results showed that both AI models performed well for prediction of the bubble dissolution time of the SLS process but EDT Bagged (MAE = 9.14, RMSE = 23.78, and R² = 0.990) outperforms EDT Boosted (MAE = 79.62, RMSE = 141.16, and R² = 0.856). The sensitivity analysis results showed that three most important input of the prediction problem were diffusion coefficient, initial concentration and bubble initial size. The present study simplified the mechanical problem of the SLS process in considering the dissolution of a single spherical gas. Therefore, a perspective of the study is to consider 2D and 3D problem to mimic the real bubble shrinkage problem. Other hybrid AI models should also be studied to improve the predictive capability. In the case of simplified approach, the computation time was estimated about 2 min for each dataset. However, the AI approaches would be more significant and useful when considering the prediction of 2D and 3D problem, where computation time would be extensively high for the generation of one dataset. Investigation of other sensitivity analysis method such as Global Sensitivity Analysis (GSA) [65] will also be a perspective of this study.

**Author Contributions:** Conceptualization, H.-B.L., T.-T.L. and E.M.; Methodology, H.-B.L. and B.T.P.; Software, H.-B.L. and T.-T.L.; Validation, H.-B.L., V.M.L. and T.-T.L.; Formal Analysis, H.-B.L. and E.M.; Investigation, H.-B.L. and T.-T.L.; Resources, H.-B.L., E.M. and G.R.; Data Curation, H.-B.L., E.M., M.D. and G.R.; Writing—Original Draft Preparation, H.-B.L., T.-T.L., V.M.L. and B.T.P.; Writing—Review & Editing, H.-B.L. and E.M.; Visualization, H.-B.L., T.-T.L. and V.M.L.; Supervision, H.-B.L., E.M. and G.R.; Project Administration, H.-B.L. and E.M.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Schmid, M. *Laser Sintering with Plastics: Technology, Processes, and Materials*; Carl Hanser Verlag GmbH & Co. KG: München, Germany, 2018; ISBN 978-1-56990-683-5.

2. Nelson, J.C. *Selective Laser Sintering: A Definition of the Process and an Empirical Sintering Model*; University of Texas at Austin: Austin, TX, USA, 1993.

3. Oter, Z.C.; Coskun, M.; Akca, Y.; Surmen, O.; Yilmaz, M.S.; Ozer, G.; Tarakci, G.; Khan, H.M.; Koc, E. Support optimization for overhanging parts in direct metal laser sintering. *Optik* **2019**, *181*, 575–581. [CrossRef]

4. Cheng, B.; Chou, K. Geometric consideration of support structures in part overhang fabrications by electron beam additive manufacturing. *Comput. Aided Des.* **2015**, *69*, 102–111. [CrossRef]

5. Calignano, F. Design optimization of supports for overhanging structures in aluminum and titanium alloys by selective laser melting. *Mater. Des.* **2014**, *64*, 203–213. [CrossRef]

6. Yuan, S.; Zheng, Y.; Chua, C.K.; Yan, Q.; Zhou, K. Electrical and thermal conductivities of MWCNT/polymer composites fabricated by selective laser sintering. *Compos. Part A Appl. Sci. Manuf.* **2018**, *105*, 203–213. [CrossRef]

7. Josupeit, S.; Ordia, L.; Schmid, H.-J. Modelling of temperatures and heat flow within laser sintered part cakes. *Addit. Manuf.* **2016**, *12*, 189–196. [CrossRef]

8. Kontopoulou, M.; Vlachopoulos, J. Bubble dissolution in molten polymers and its role in rotational molding. *Polym. Eng. Sci.* **1999**, *39*, 1189–1198. [CrossRef]

9. Griskey, R. *Polymer Process Engineering*; Springer Science & Business Media: Dordrecht, The Netherlands, 2012.

10. Doremus, R.H. Diffusion of Oxygen from Contracting Bubbles in Molten Glass. *J. Am. Ceram. Soc.* **1960**, *43*, 655–661. [CrossRef]
11. Wei, Q.; Zhong, J.; Xu, Z.; Xu, Q.; Liu, B. Microstructure evolution and mechanical properties of ceramic shell moulds for investment casting of turbine blades by selective laser sintering. *Ceram. Int.* **2018**, *44*, 12088-12097. [CrossRef]

12. Mokrane, A.; Boutaous, M.; Xin, S. Process of selective laser sintering of polymer powders: Modeling, simulation, and validation. *C. R. Méc.* **2018**, *346*, 1087-1103. [CrossRef]

13. Venuvinod, P.K.; Ma, W. Selective Laser Sintering (SLS). In *Rapid Prototyping: Laser-based and Other Technologies*; Venuvinod, P.K., Ma, W., Eds.; Springer: Boston, MA, USA, 2004; pp. 245-277. ISBN 978-1-4757-6361-4.

14. Bourell, D.L.; Watt, T.J.; Leigh, D.K.; Fulcher, B. Performance Limitations in Polymer Laser Sintering. *Phys. Procedia* **2014**, *56*, 147-156. [CrossRef]

15. Schmid, M.; Wegener, K. Additive Manufacturing: Polymers Applicable for Laser Sintering (LS). *Procedia Eng.* **2016**, *149*, 457-464. [CrossRef]

16. Wudy, K.; Lanzl, L.; Drummer, D. Selective Laser Sintering of Filled Polymer Systems: Bulk Properties and Laser Beam Material Interaction. *Phys. Procedia* **2016**, *83*, 991-1002. [CrossRef]

17. Ly, H.-B.; Monteiro, E.; Dal, M.; Regnier, G. On the factors affecting porosity dissolution in selective laser sintering process. *AIP Conf. Proc.* **2018**, *1960*, 120014.

18. Gogos, G. Bubble removal in rotational molding. *Polym. Eng. Sci.* **2004**, *44*, 388-394. [CrossRef]

19. Weinberg, M.C. Surface tension effects in gas bubble dissolution and growth. *Chem. Eng. Sci.* **1981**, *36*, 137-141. [CrossRef]

20. Duda, J.; Vrentas, J. Mathematical analysis of bubble dissolution. *AIChE J.* **1969**, *15*, 351-356. [CrossRef]

21. Naji Meidani, A.R.; Hasan, M. Mathematical and physical modelling of bubble growth due to ultrasound. *Appl. Math. Model.* **2004**, *28*, 333-351. [CrossRef]

22. Haykin, S. *Neural Networks: A Comprehensive Foundation.*, 2nd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1998; ISBN 978-0-13-273350-2.

23. Jang, J.-S.R. *Neuro-Fuzzy and Soft Computing: A Computational Approach to Learning and Machine Intelligence*; Prentice Hall: Upper Saddle River, NJ, USA, 1997; ISBN 978-0-13-261066-7.

24. Witten, I.H.; Frank, E.; Hall, M.A.; Pal, C.J. *Data Mining: Practical Machine Learning Tools and Techniques*, 4th ed.; Morgan Kaufmann: San Francisco, CA, USA, 2016; ISBN 978-0-12-804357-8.

25. Cosmin, A.; Elena, A.; Naif, A.; Timon, R. Artificial Neural Network Methods for the Solution of Second Order Boundary Value Problems. *Comput. Mater. Contin.* **2019**, *59*, 345-359.

26. Yusoff, N.I.M.; Ibrahim Alhamali, D.; Ibrahim, A.N.H.; Rosyidi, S.A.P.; Abdul Hassan, N. Engineering characteristics of nanosilica/polymer-modified bitumen and predicting their rheological properties using multilayer perceptron neural network model. *Constr. Build. Mater.* **2019**, *204*, 781-799. [CrossRef]

27. Hasnaoui, H.; Krea, M.; Roizard, D. Neural networks for the prediction of polymer permeability to gases. *J. Membr. Sci.* **2017**, *541*, 541-549. [CrossRef]

28. Diaconescu, R.-M.; Barbuta, M.; Harja, M. Prediction of properties of polymer concrete composite with tire rubber using neural networks. *Mater. Sci. Eng. B* **2013**, *178*, 1259-1267. [CrossRef]

29. Khader, M.H.; Hamid, G.; Xiaoying, Z.; Naif, A.; Timon, R. Computational Machine Learning Representation for the Flexoelectricity Effect in Truncated Pyramid Structures. *Comput. Mater. Contin.* **2019**, *59*, 79-87.

30. Dao, D.V.; Trinh, S.H.; Ly, H.-B.; Pham, B.T. Prediction of Compressive Strength of Geopolymer Concrete Using Entirely Steel Slag Aggregates: Novel Hybrid Artificial Intelligence Approaches. *Appl. Sci.* **2019**, *9*, 1113. [CrossRef]

31. Xu, D.; Pop-Iliev, R.; Park, C.B.; Fenton, R. Fundamental study of CBA-blown bubble growth and collapse under atmospheric pressure. *J. Cell. Plast.* **2005**, *41*, 519-538. [CrossRef]

32. Wu, S. *Polymer Interface and Adhesion*; CRC Press: Boca Raton, FL, USA, 1982.

33. Bird, R.B.; Stewart, W.E.; Lightfoot, E.N. *Transport Phenomena*; John Wiley & Sons: Hoboken, NJ, USA, 2007.

34. Durrill, P.L.; Griskey, R.G. Diffusion and solution of gases into thermally softened or molten polymers: Part II. Relation of diffusivities and solubilities with temperature pressure and structural characteristics. *AIChE J.* **1969**, *15*, 106-110. [CrossRef]

35. Khorsheed, M.S.; Al-Thubaity, A.O. Comparative evaluation of text classification techniques using a large diverse Arabic dataset. *Lang Resour. Eval.* **2013**, *47*, 513-538. [CrossRef]

36. Leema, N.; Nehemiah, H.K.; Kannan, A. Neural network classifier optimization using Differential Evolution with Global Information and Back Propagation algorithm for clinical datasets. *Appl. Soft Comput.* **2016**, *49*, 834-844. [CrossRef]
37. Quinlan, J.R. Induction of decision trees. *Mach. Learn.* **1986**, *1*, 81–106. [CrossRef]

38. Amirabadizadeh, A.; Nezami, H.; Vaughn, M.G.; Nakhaee, S.; Mehrpour, O. Identifying risk factors for drug use in an Iranian treatment sample: A prediction approach using decision trees. *Subst. Use Misuse* **2018**, *53*, 1030–1040. [CrossRef]

39. Khosravi, K.; Pham, B.T.; Chapi, K.; Shirzadi, A.; Shahabi, H.; Revhaug, I.; Prakash, I.; Tien Bui, D. A comparative assessment of decision trees algorithms for flash flood susceptibility modeling at Haraz watershed, northern Iran. *Sci. Total Environ.* **2018**, *627*, 744–755. [CrossRef]

40. Thai Pham, B.; Bui, D.T.; Prakash, I. Landslide susceptibility modelling using different advanced decision trees methods. *Civ. Eng. Environ. Syst.* **2019**, *1–19*. [CrossRef]

41. Wang, J.; Li, P.; Ran, R.; Che, Y.; Zhou, Y. A Short-Term Photovoltaic Power Prediction Model Based on the Gradient Boost Decision Tree. *Appl. Sci.* **2018**, *8*, 689. [CrossRef]

42. Akin, M.; Hand, C.; Eyduran, E.; Reed, B.M. Predicting minor nutrient requirements of hazelnut shoot cultures using regression trees. *Plant Cell Tissue Organ Cult.* (PCTOC) **2018**, *132*, 545–559. [CrossRef]

43. Hastie, T.; Tibshirani, R.; Friedman, J. *The Elements of Statistical Learning: Data Mining, Inference, and Prediction*, 2nd ed.; Springer Series in Statistics; Springer: New York, NY, USA, 2009; ISBN 978-0-387-84857-0.

44. Zhang, Y.; Haghani, A. A gradient boosting method to improve travel time prediction. *Transp. Res. Part C Emerg. Technol.* **2015**, *58*, 308–324. [CrossRef]

45. Austin, P.C.; Tu, J.V.; Ho, J.E.; Levy, D.; Lee, D.S. Using methods from the data-mining and machine-learning literature for disease classification and prediction: A case study examining classification of heart failure subtypes. *J. Clin. Epidemiol.* **2013**, *66*, 398–407. [CrossRef] [PubMed]

46. Winkler, D.; Haltmeier, M.; Kleidorfer, M.; Rauch, W.; Tscheikner-Gratl, F. Pipe failure modelling for water distribution networks using boosted decision trees. *Struct. Infrastruct. Eng.* **2018**, *14*, 1402–1411. [CrossRef]

47. Serrano-Cinca, C.; Gutiérrez-Nieto, B. Partial Least Square Discriminant Analysis for bankruptcy prediction. *Decis. Support Syst.* **2013**, *54*, 1245–1255. [CrossRef]

48. Elith, J.; Leathwick, J.R.; Hastie, T. A working guide to boosted regression trees. *J. Anim. Ecol.* **2008**, *77*, 802–813. [CrossRef]

49. Zięba, M.; Tomczak, S.K.; Tomczak, J.M. Ensemble boosted trees with synthetic features generation in application to bankruptcy prediction. *Expert Syst. Appl.* **2016**, *58*, 93–101. [CrossRef]

50. De'ath, G. Boosted Trees for Ecological Modeling and Prediction. *Ecology* **2007**, *88*, 243–251. [CrossRef]

51. Abdar, M.; Yen, N.Y.; Hung, J.C.-S. Improving the diagnosis of liver disease using multilayer perceptron neural network and boosted decision trees. *J. Med. Boil. Eng.* **2017**, *38*, 953–965. [CrossRef]

52. Chung, Y.-S. Factor complexity of crash occurrence: An empirical demonstration using boosted regression trees. *Accid. Anal. Prev.* **2013**, *61*, 107–118. [CrossRef]

53. Ghasemi, H.; Rafiee, R.; Zhuang, X.; Muthu, J.; Rabczuk, T. Uncertainties propagation in metamodel-based probabilistic optimization of CNT/polymer composite structure using stochastic multi-scale modeling. *Comput. Mater. Sci.* **2014**, *85*, 295–305. [CrossRef]

54. Vu-Bac, N.; Rafiee, R.; Zhuang, X.; Lahmer, T.; Rabczuk, T. Uncertainty quantification for multiscale modeling of polymer nanocomposites with correlated parameters. *Compos. Part B Eng.* **2015**, *68*, 446–464. [CrossRef]

55. Le, T.T.; Guilleminot, J.; Soize, C. Stochastic continuum modeling of random interphases from atomistic simulations. Application to a polymer nanocomposite. *Comput. Methods Appl. Mech. Eng.* **2016**, *303*, 430–449. [CrossRef]

56. Capillon, R.; Descliers, C.; Soize, C. Uncertainty quantification in computational linear structural dynamics for viscoelastic composite structures. *Comput. Methods Appl. Mech. Eng.* **2016**, *305*, 154–172. [CrossRef]

57. Robert, C.; Casella, G. *Monte Carlo Statistical Methods*, 2nd ed.; Springer Texts in Statistics; Springer: New York, NY, USA, 2004; ISBN 978-0-387-21239-5.

58. Mordechai, S. *Applications of Monte Carlo Method in Science and Engineering*; InTech: Rijeka, Croatia, 2011; ISBN 978-953-307-691-1.

59. Chai, T.; Draxler, R.R. Root mean square error (RMSE) or mean absolute error (MAE)?—Arguments against avoiding RMSE in the literature. *Geosci. Model Dev.* **2014**, *7*, 1247–1250. [CrossRef]

60. Menard, S. Coefficients of Determination for Multiple Logistic Regression Analysis. *Am. Stat.* **2000**, *54*, 17–24.

61. Nagelkerke, N. A note on a general definition of the coefficient of determination. *Biometrika* **1991**, *78*, 691–692. [CrossRef]
62. Willmott, C.J.; Matsuura, K. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance. *Clim. Res.* **2005**, *30*, 79–82. [CrossRef]

63. Dao, D.V.; Ly, H.-B.; Trinh, S.H.; Le, T.-T.; Pham, B.T. Artificial Intelligence Approaches for Prediction of Compressive Strength of Geopolymer Concrete. *Materials* **2019**, *12*, 983. [CrossRef] [PubMed]

64. Epstein, P.S.; Plesset, M.S. On the stability of gas bubbles in liquid-gas solutions. *J. Chem. Phys.* **1950**, *18*, 1505–1509. [CrossRef]

65. Hamdia, K.M.; Msekh, M.A.; Silani, M.; Thai, T.Q.; Budarapu, P.R.; Rabczuk, T. Assessment of computational fracture models using Bayesian method. *Eng. Fract. Mech.* **2019**, *205*, 387–398. [CrossRef]

<figure><img src="image_10.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).