2.1. Preparation of fibers and biocomposites
In this study, WF fibers were used. This fiber is easily grown locally in the region of Skikda (Algeria). WF fibers are extracted manually and then immersed in distilled water to wash off surface impurities. To remove moisture, the fibers are dried under natural conditions (at room temperature, 20°C for 10 days) before the fabrication of the biocomposite. The polymer used as a matrix is polyethylene (HDPE) reference PPC10642. It was supplied by SABIC Petrochemicals (HDPE blow molding, 95% purity and above, an average MFI of 0.4 and an average density of 0.96).
The biocomposite is manufactured on two Thermotron-C.W. Brabender type rollers (Model T303). First, 20% by mass HDPE is melted on the rollers at 170° C [22]. The WF fibers are cut (between 5 and 10 mm) and the rest of the HDPE is added and mixed for 7 min at 60 rpm. Then, the biocomposite is mixed 5 to 6 times for 5 minutes to improve the homogeneity of the material. After this, the processed biocomposite is removed from the rollers and cut to fit the mold size. The samples are prepared in a mold that is held at 190 ± 3°C using a Dake brand press for 20 minutes under a pressure of 20 MPa. Finally, the mold is cooled to 60°C using cold water. Table 1 shows the mass compositions of the different formulations produced.
Table 1
Mass composition of the different formulations elaborated in this work.
Designation
|
Biocomposites
|
HDPE/10%WF
|
High density polyethylene HDPE with10% of WF by weight
|
HDPE /20%WF
|
High density polyethylene HDPE with 20% of WF by weight
|
HDPE /30%WF
|
High density polyethylene HDPE with 30% of WF by weight
|
2.2. Water absorption
The study of the effect of WF fiber content on the diffusion kinetics of HDPE/WF biocomposites by total immersion in distilled water at room temperature consists of following the mass evolution of the samples over time, measured at regular intervals over a total period of 21 days (500h). At the time of the measurement of the mass, the samples were taken out of the bath and lightly wiped with an absorbent paper to eliminate the film of water present on the surface. Weighing was done on an electronic balance (accuracy of measurement = 0.0001 g). Water absorption tests were performed according to the ASTM standard method [23]. During aging, the mass gain of a biocomposite sample at time t expressed as a percentage, Mt (%), was determined according to Eq. 1.
$$Mt\left(\%\right)= \frac{\left(mt-mo\right)}{mo}\times 100$$
1
with mt: weight of the composite at the time t; m0: initial mass of the dried sample.
2.3. Fick’s model
Mathematical modeling of the water diffusion process in biocomposites is very important to control the diffusion mechanism. The numerous diffusion phenomena present in nature are described by Fick’s laws and can be characterized by their diffusion coefficient D. We have the first Fick’s law given by the following equation:
where D: diffusion coefficient of the medium inmm2/s;
C: concentration of the solvent in the medium
Considering a thin plate of thickness h, in which the solvent diffuses, initially at the concentration C0, and whose surfaces are kept at the uniform concentration C1, the spatial and temporal evolution of the solvent concentration is given by
$$\frac{C-{C}_{0}}{C-{C}_{0}}=1-\frac{4}{\pi }\sum _{n=0}^{\infty }\frac{{\left(-1\right)}^{n}}{\left(2n+1\right)}exp\left(-D\frac{{\left(2n+1\right)}^{2}}{{h}^{2}}{\pi }^{2}\bullet t\right)\text{cos}\left(\frac{\left(2n+1\right)\pi }{h}\right)$$
3
where D: diffusion coefficient; t: aging time; h: thickness of the plate
If we note Ms as the mass of water absorbed after an infinite time, then Eq. 3 is written as [24]:
$$\frac{{M}_{t}}{{M}_{s}}=1-\frac{8}{\pi }\sum _{n=1}^{\infty }\frac{1}{{\left(2n+1\right)}^{2}}exp\left(-D\frac{{\left(2n+1\right)}^{2}}{{h}^{2}}{\pi }^{2}\bullet t\right)$$
4
where Mt: the water content at time t; Mm: the maximum mass of water in equilibrium.
When Mt / Mm is less than 0.6, Eq. 4 becomes approximately [14]
$$\frac{{M}_{t}}{{M}_{s}}=\frac{4}{h}\sqrt{\frac{Dt}{\pi }}$$
5
When Mt / Mm is greater than 0.6, the equation describing the moisture absorption curve is [14]:
$$\frac{{M}_{t}}{{M}_{s}}=1-exp\left[-7.3{\left(\frac{Dt}{{h}^{2}}\right)}^{0.75}\right]$$
6
The diffusion coefficient D can be deduced from Eq. 5 as
$$D=\pi {\left(\frac{k}{4{M}_{m}}\right)}^{2}$$
7
where k is the slope of the linear part of the curve 𝑀𝑡 = (√𝑡⁄ℎ)
2.4. Response surface methodology modeling
A statistical design of experiment, called Response Surface Methodology (RSM), allows process variables to be varied simultaneously, as compared to traditional testing, in order to derive the relationship among these variables. RSM offers a faster and more cost-effective method of collecting research findings than conventional experimentation with one variable as well as with a full variable [25–28]. In the current research, the effect of WF fiber content on the diffusion kinetics of HDPE/WF biocomposites by full immersion in distilled water at ambient conditions of the samples was investigated. The studied variables were fiber content (X1) and immersion time (X2). In Table 2, the variables and their corresponding values are reported. A two-factor central composite design (CCD) was used, involving 16 series of experiments in the RSM analysis, computed according to Eq. 8:
Table 2
Design of experiments for RSM model.
N°
|
Factors
|
Notation
|
Units
|
Levels
|
|
|
|
|
Low level
|
Intermediate level
|
High level
|
1
|
Immersion time
|
Time
|
hour
|
2
|
408
|
800
|
2
|
Content fibers
|
W
|
%
|
10
|
20
|
30
|
$$N= {2}^{k}+2k+{k}_{c}$$
8
In which N stands for the sum of experimental of tests to be performed with k standing for the number of parameters used and kc stands for center-repeated series.. An experimental model was developed to correlate the response with the two factors in the process based on a quadratic polynomial model as indicated in Eq. 9.
$$Y= {B}_{0 }+ \sum {B}_{i}{X}_{i}+ \sum {B}_{ii}{X}^{2}+ \sum {B}_{ij}{X}_{i}{X}_{j}+E$$
9
In which Y represents the value of the predicted response, where B0 is the constant value and Bi, Bii and Bij are the coefficients of the linear, quadratic and interactive terms, respectively [29].
The experimental plan and the measured response (% water absorption) for the two parameters and the 16 experimental runs produced are presented in Table 3. To perform the RSM regression analysis, the Design Expert software was used to optimize the biocomposite composition data generated with the input data. An ANOVA analysis including quadratic, linear and coefficient of interaction is carried out for the model’s statistical test using F-test for the empirical interrelationship of the input parameters and the output model. For testing the model fit performance, each individual term of the model was statistically tested, confirming F-value significance with p < 0.05. Also, in order to verify the goodness of the proposed polynomial, the R2, the adequate accuracy and predicted R2 and the adjusted R2 values of the models were acquired. Both the response surface map and the contour map were constructed to display the input-output interactions.
Table 3
Experimental data with RSM and ANN model results for water absorption of HDPE/WF biocomposite.
Experiment number
|
Input variables
|
Output variables
|
WF
(%)
|
Time
(h)
|
EXP
(%)
|
RSM
(%)
|
ANN
(%)
|
1
|
0
|
0
|
0.0000
|
0.0818
|
0.0002
|
2
|
10
|
2
|
0.4325
|
0.4064
|
0.4371
|
3
|
20
|
2
|
0.9013
|
0.9566
|
0.8985
|
4
|
30
|
2
|
1.6516
|
1.7378
|
1.6511
|
5
|
10
|
24
|
0.5627
|
0.4789
|
0.5552
|
6
|
20
|
24
|
1.0141
|
1.0457
|
1.0240
|
7
|
30
|
24
|
1.8308
|
1.8434
|
1.7763
|
8
|
10
|
48
|
0.6145
|
0.5543
|
0.6598
|
9
|
20
|
48
|
1.1251
|
1.1392
|
1.1659
|
10
|
30
|
48
|
1.9661
|
1.9550
|
1.9649
|
11
|
10
|
120
|
0.8504
|
0.7578
|
0.8471
|
12
|
20
|
120
|
1.4063
|
1.3970
|
1.4831
|
13
|
30
|
120
|
2.2632
|
2.2671
|
2.2359
|
14
|
10
|
288
|
1.1432
|
1.1000
|
1.1385
|
15
|
20
|
288
|
1.8832
|
1.8660
|
1.8694
|
16
|
30
|
288
|
2.9807
|
2.8628
|
3.0051
|
17
|
10
|
408
|
1.1431
|
1.2308
|
1.1030
|
18
|
20
|
408
|
2.0111
|
2.0873
|
2.0191
|
19
|
30
|
408
|
3.2515
|
3.1746
|
3.2532
|
20
|
10
|
572
|
1.1404
|
1.2565
|
1.1481
|
21
|
20
|
572
|
2.1448
|
2.2366
|
2.1438
|
22
|
30
|
572
|
3.4793
|
3.4476
|
3.4796
|
23
|
10
|
800
|
1.1408
|
0.9982
|
1.1395
|
24
|
20
|
800
|
2.1480
|
2.1503
|
2.1473
|
25
|
30
|
800
|
3.4801
|
3.5333
|
3.4417
|
2.5. Artificial neural network
To train the network, a multilayer perceptron (MPL) was used, and the algorithm of back propagation training was used for modelling. It consisted of an input layer, a hidden layer, and an output layer. In the input layer, the variables were the fiber content and the immersion time of the biocomposite in water, while, in the output layer, it was the water absorption. To determine the ideal neural number in the hidden layer, a test series of approaches was used to obtain the number of neurons that would yield the least root mean square error (RMSE) and the most correlation coefficient (R2) [30–32]. This approach was intended to provide a minimum discrepancy in predicted and experimental outputs as well as reduce the opportunity for excessive model fitting. High as well as low neural numbers were omitted, as they cause complication in fitting and reduced rate of convergence, respectively [27, 33]. To train the network, 70% of the datasets were utilized. Of the remaining, 15% were utilized to test the network, and the other 15% were utilized for result validation.
The accuracy of the model given by ANN was validated by the values of correlation coefficient (R2) and the mean squared error (MSE). The lower the values of MSE, the more precise the ANN model is in the prediction. The three parameters were given by equations 10–12 [34].
$$MSE=\frac{1}{n}{\sum }_{i=1}^{n}\left|{\left({Y}_{Predicted}-{Y}_{Experiment}\right)}^{2}\right|$$
10
$${R}^{2}=1-\frac{{\sum }_{i=1}^{n}{\left({Y}_{Predicted}-{Y}_{Experiment}\right)}^{2}}{{\sum }_{i=1}^{n}{\left({Y}_{Predicted}-{Y}_{Mean}\right)}^{2}}$$
11
$$RMSE= \sum _{1}^{{n}}\sqrt{\frac{{({X}_{Predicted}-{X}_{Experiment})}^{2}}{{n}}}$$
12