*Article* **Influence of Moisture Diffusion on the Dynamic Compressive Behavior of Glass/Polyester Composite Joints for Marine Engineering Applications**

**Oumnia Lagdani 1, Mostapha Tarfaoui 1,2, Marwane Rouway 3,4, Houda Laaouidi 1, Sara Jamoudi Sbai 4,5, Mohamed Amine Dabachi 6, Abdelwahed Aamir <sup>7</sup> and Mourad Nachtane 8,\***


**Abstract:** Thermoset polymers offer great opportunities for mass production of fiber-reinforced composites and are being adopted across a large range of applications within the automotive, aerospace, construction and renewable energy sectors. They are usually chosen for marine engineering applications for their excellent mechanical behavior, including low density and low-cost compared to conventional materials. In the marine environment, these materials are confronted by severe conditions, thus there is the necessity to understand their mechanical behavior under critical loads. The high strain rate performance of bonded joints composite under hygrothermal aging has been studied in this paper. Initially, the bonded composite specimens were hygrothermal aged with the conditions of 50 ◦C and 80% in temperature and relative humidity, respectively. After that, gravimetric testing is used to describe the moisture diffusion properties for the adhesively bonded composite samples and exhibit lower weight gain for this material. Then, the in-plane dynamic compression experiments were carried out at different impact pressures ranging from 445 to 1240 s−<sup>1</sup> using the SHPB (Split Hopkinson Pressure Bar) technique. The experimental results demonstrated that the dynamic behavior varies with the variation of strain rate. Buckling and delamination of fiber are the dominant damage criteria observed in the sample during in-plane compression tests.

**Keywords:** moisture diffusion; compressive behavior; SHPB technique; glass/polyester composite joints

### **1. Introduction**

Glass fiber composites are usually chosen for marine engineering applications due to their excellent mechanical behavior and particularly, their low densities and cost compared to traditional materials [1–5]. However, the most structural requirement of composites is the ability to maintain a high proportion of its load-carrying performance over prolonged time under extreme environmental circumstances, such as temperature changes, humidity, oxidation, microbial attack, etc. [6–8]. Over the past few decades, the strengthening and repair of marine engineering structures using glass fiber reinforced polymer (GFRP) laminates have gained much attention [9,10]. As the science of polymers has advanced, composite bonding has developed into a major joining method for many applications, because it affords a lot of benefits compared to metallic joints, such as high fatigue resistance,

**Citation:** Lagdani, O.; Tarfaoui, M.; Rouway, M.; Laaouidi, H.; Sbai, S.J.; Dabachi, M.A.; Aamir, A.; Nachtane, M. Influence of Moisture Diffusion on the Dynamic Compressive Behavior of Glass/Polyester Composite Joints for Marine Engineering Applications. *J. Compos. Sci.* **2022**, *6*, 94. https:// doi.org/10.3390/jcs6030094

Academic Editor: Stelios K. Georgantzinos

Received: 13 February 2022 Accepted: 11 March 2022 Published: 16 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 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 (https:// creativecommons.org/licenses/by/ 4.0/).

galvanic corrosion elimination and uniform stress distribution across the joint [11–14]. The exploring of the hygrothermal aging effect on the mechanical performance of bonded joints has increasingly become recognized in the last years. The water aging influence on interface debonding has been studied in the literature. Komai et al. [15] demonstrated that water absorption leads to strength loss of the fiber/matrix interfaces of aligned reinforced carbon/epoxy composites. Ashik [16] conducted an experimental study of moisture uptake and mechanical characteristics of hybrid natural and glass fiber reinforced composites. The result shows that the addition of coconut fiber can improve strength and serve as an alternative material to glass fiber. Selzer and Friedrich [17] observed that the bonding of fiber/resin deteriorates with increasing moisture absorption. Due to the individual set of properties that synthetic polymers offer, their applications in the outdoor environment are steadily growing.

Polymers have replaced conventional materials in the marine industry. Many research papers have shown the impact of seawater on polymeric composites which leads to a degradation of performance [18]. Buehler and Seferis [19] noted interface delamination in carbon-glass reinforced epoxy due to water absorption. However, there are limited works on damage mechanisms that can change fatigue resistance, even in the absence of important seawater deterioration. So, it is substantial to better comprehend the extent and mechanisms of seawater deterioration. Schutte [20] demonstrated that hygrothermal decomposition of glass fiber composites is primarily due to the decomposition of the glass fibers, the polymer matrix and the interface. However, Mijovi´c et al. [21] showed that vinylester reinforced composites display excellent durability over epoxy-based composites. Khalilullah [22] characterized the effect of moisture and hygroscopic bulking of silicone/phosphorus composite film, also the moisture uptake increases with temperature for the film. Therefore, vinylester and epoxy resins are usually chosen due to their low cost, thermal aging resistance and high sustainability [23]. Jiang et al. [24,25] experimentally and numerically investigated the impact of humidity and hygrothermal degradation in GFRP adhesives. The results revealed that samples lost mass under 40 ◦C conditions. Grace et al. [26] presented a new approach to characterize anisotropic moisture absorption in polymer composites using gravimetric absorption. Benyahia et al. [27] investigated the mechanical behavior of composite pipes at different temperatures between −40 ◦C and 80 ◦C where the pipes degrade when temperature increases. Sassi and his colleagues [28–32] investigated the dynamic behavior of bonded composite joints at elevated strain rates using an SHPB machine, the experimental results have shown a strong material sensitivity to strain rates. Moreover, damage investigations have revealed that the failure mainly occurred in the adhesive/adherent interface because of the brittle nature of the polymeric adhesive. Results have shown good agreement with the dependency of the dynamic parameters on strain rates. Despite these efforts, the influence of moisture diffusion on the dynamic compressive behavior of glass/polyester composite joints remains immature. Hence, the importance of understanding the mechanical behavior of these materials that are subjected to critical loads for an extended period. In this context, this paper highlights the impact of hygrothermal aging on the compressive responses of bonded composite joints at an elevated strain rate using the Split Hopkinson Pressure Bar (SHPB) machine.

In summary, the outline of this article is as follows: Section 2 provides the material description and methods, followed by the experimental methodology utilized for the mechanical characterization in Section 3. Finally, Section 4 draws pertinent conclusions as well as some prospects for future developments.

### **2. Materials and Methods**

### *2.1. Materials*

The composite used in this work was made of a 45◦ biaxial glass fiber mat with 0.29 mm in thickness reinforced polyester polymer. The composite part was joined with an adhesive of polyvinylester with 1 mm in thickness purchased from NORPOL. The mechanical properties of materials are listed in Table 1 [28]. The polyester polymer is used in the maritime industry as an equivalent to other types of polymers, due to their cost being slightly less than epoxy resin and also it has an excellent performance in corrosive environments and at extreme temperatures; therefore, it is perfectly adapted to naval applications.



### *2.2. Methods*

### 2.2.1. Hygrothermal Aging Test

The samples were tested and performed by hygrothermal aging to investigate the moisture absorption properties and to measure the gain in weight of the samples in the function of time. The dimension of samples is 13 mm × 13 mm × 9 mm as shown in Figure 1. To perform the SHPB test, the samples must be prepared where both sides of the specimen are recommended to be planned and oriented with a high level of accuracy. The mass of each sample was calculated every 96 h in a regular period using a Precisa XT220A analytical balance with a tolerance of 0.0001 g. To record the weight variation after a specified interval of time, each specimen is taken out of the conditioned chamber, weighed quickly and returned into the chamber. The moisture uptake *M* absorbed by each sample is measured as a function of its initial weight *wi* and final weight *wf* [25,29] as described below:

$$M = \frac{w\_f - w\_i}{w\_i} \times 100\tag{1}$$

**Figure 1.** Adhesively bonded composite joints sample.

Before studying the aging of an adhesive, it is necessary to ensure that the polymerization is complete to avoid crosslinking phenomena during aging [2]. The DSC (Differential Scanning Calorimetry) test can determine an adequate curing cycle (polymerization of the glue). Likewise, these tests carried out on polymerized samples make it possible to determine a glass transition temperature Tg. From these data, it is then possible to set aging conditions. In general, it is necessary, for the aging temperature, to be lower than the Tg of an aged material saturated with water to overcome complex aging phenomena in a rubbery state. In addition, it is recommended that the aging temperature be much lower than the post-cure temperature as shown in Figure 2. The bonded composite specimen was

tested under the moisture diffusion effect with 50 ◦C and 80% of temperature and relative humidity, respectively, for different periods.

**Figure 2.** Glass transition temperature of the specimen.

### 2.2.2. Moisture Absorption

Many researchers investigated the performance degradation of marine water on polymeric composites [30,31]. However, the marine environment could lead to chemical bond hydrolysis and mechanical degradation and also a loss in interfacial stress transmission due to matrix plastification, in addition, the interface changes between the fiber and the matrix [18,32].

Various diffusion models have been developed to comprehend the behavior of homogeneous and heterogeneous materials, the best known and most commonly used of which is the one-dimensional isotropic Fick model [33] as seen in Figure 3, but divergence from this model are frequently remarked. On the other hand, the Langmuir-type model can accurately describe the water absorption of fiber-reinforced epoxy composite [34].

**Figure 3.** Comparative diagram of diffusion phenomena. (**a**) Homogeneous material, (**b**) Heterogeneous material.

Experimental data were expressed as the moisture uptake *M*(*t*). However, the study was conducted using average concentration *C*(*t*) as the main parameter. To compare the experimental and numerical results, the *M*(*t*) was converted to *C*(*t*). The mean concentration was evaluated using Equation (2), where the densities of the material and water at 60 ◦C are *ρmat* = 1.9 g/cm<sup>3</sup> and *ρwater* = 0.9 g/cm3, respectively.

$$\mathcal{C}(t) = \frac{\rho\_{mat}}{\rho\_{cau}} M(t) \tag{2}$$

The average concentration *Ci* is computed by using the finite element model at each integration point. So, to compare the results, *Ci* is evaluated at each increment point for the entire test tube using Equation (3), where *Ci* and *Vi* represent the concentration and volume at the integration point, respectively. The equation was executed using a Python script that extracts the appropriate information from the Abaqus results database.

$$\mathbf{C}\_{i}(t) = \frac{\sum\_{i=1}^{n} \mathbf{C}\_{i} V\_{i}}{\sum\_{i=1}^{n} V\_{i}} \tag{3}$$

Figure 4 presents the variation of the moisture uptake curve in the function of the square root of time, which usually concurs with the Fick law. The numerical results remain below the experimental and analytical values which confirm the model correlates satisfactorily with the theoretical and experimental results. The experimental data presented in Figure 4 represents the average of 10 tested samples. The moisture process of bonded joints may be defined into three period phases. In the first phase, the moisture uptake rate grows rapidly with time, achieving to moisture absorption of 0.239 wt% in 49 h. The rate of water absorption in bonded joints was changed linearly with aging time. The moisture uptake in this phase is principal because of composite material defects, such as voids and micro-cracks in the matrix. Into the second phase, the moisture decelerates until the saturation point at 0.2875 wt% in 400 h. The moisture absorption mechanisms include delamination of the fiber/resin interface and hydrolysis of the resin. At the third phase, the moisture diffusion of the composite has achieved equilibrium after approximately 420 h of aging.

**Figure 4.** Bonded joint moisture absorption, a comparison of analytical (Fick) and numerical models with experimental data.

Numerous authors have described the volume variation of the specimens under hygrothermal loading. Gazit [35] showed that the increase in sample weight is linearly proportional to the dimensions of all reinforced samples. According to the polarity of water, it can create hydrogen bonds with the hydroxyl groups. Consequently, the hydrogen bonds between chains may be interrupted to raise the length of the hydrogen bond between segments. These mechanisms reduce the glass transition temperature Tg of the wetted specimen and is known as plasticization and swelling of the polymer matrix, which leads to microstructural damage, such as fibers delamination and matrix cracking [36,37]. In contrast, matrix swelling is important in the region with high resin uptake. This induces decoherence of the fiber/resin interface. In another work, Lee [38] demonstrates that the absorption of water produces plasticization and swelling of the resin and a reduction of the glass transition temperature. These generally affect the modulus of the composite material and can be accelerated by increasing the temperature. As observed in Figure 5, the numerical results of moisture diffusion at different aging times (0 h, 216 h, 648 h and 864 h).

**Figure 5.** Illustration of numerical results at different aging times.

### **3. Dynamic Response in Moisture Absorption**

### *3.1. Dynamic Compression Test*

Many techniques were employed to evaluate the mechanical characteristics of the composites at high strain rates. The most utilized machines were the Split Hopkinson Pressure Bar (SHPB) for strain rates between 500 and 104 s−<sup>1</sup> [39–41]. Figure 6 illustrates a photograph of the used SHPB machine with different parts of the system including striker bar, input, and output bars. To register the experiment signals, incident and transmitted bars attached with two strain gauges were installed. The samples were embedded between two bars with a diameter of 20 mm in each bar. The striker, input and output bars have a length of 0.8 m, 3 m and 2 m, respectively. The bars are correctly aligned and can move easily on the base. The striker is projected onto the input bar at a specific speed.

**Figure 6.** SHPB machine.

### *3.2. Mechanical Behavior*

The thermomechanical properties of a composite are affected by moisture and temperature. However, an increase in moisture and temperature can augment the molecular displacement of the material and can significantly alter the shape and volume of composites. Such deformation is further sophisticated, especially with the dynamic compression of the material. The impact pressure effect on the stress/strain behavior of bonded composite under in-plane dynamic compression is shown in Figure 7.

**Figure 7.** Stress/strain curves of the sample under in-plane loading for different pressure and aging times. (**a**) 1 bar. (**b**) 2 bar. (**c**) 3 bar. (**d**) 4 bar.

Figure 7 presents the stress–strain responses for various impact pressures from 1 bar to 4 bar with various aging times of 0 h, 216 h, 648 h, and 864 h. It can be noted that the dynamic compressive behavior of the composite is strongly affected by strain rate and aging time. Moreover, the effect of moisture on the stress–strain curves is more significant at a low strain rate than at a high strain rate. This clearly shows that the strain rate response is sensitive to the entry pressure P in the chamber of compressed air. At the ultimate stress, the sample damage appears and loses its loadbearing capacity. It is apparent that for all impact tests, the stress–strain tendency was approximately the same during the linear elastic phase, with no damage at small strains. In addition, it can be shown that in the undamaged case, the samples regain their original state with negligible plastic deformation. On the other hand, in the damaged samples, the first peak is observed in the elastic response which shows the beginning of microscopic damage modes, such as matrix cracking. The matrix cracking of samples with in-plan direction encourages micro-buckling and fiber twisting, and leads to debonding, fiber separation, delamination and fiber breaking towards the final failure observed with the appearance of the second peak for the strain rate curves. It can be concluded that the first region of nonlinearity in the stress–strain curves are mainly due to the viscoelastic nature of the polyester polymer, while the nonlinearity noticed before the final failure is caused by the cracking of the matrix, and the elastic modulus reduced with aging time due to the rising of moisture uptake in the composite.

Figure 8 shows the evolution of strain rate as a function of time for different aging times from 0 h to 864 h and with various impact pressures on the bonded joints. It can be observed that the damage becomes significant only for high-pressure impact due to the second peak presented in the signal and it becomes increasingly important as the impact pressure increases, indicating the accumulated failure modes in the material subjected to dynamic loading. In contrast, for low impact pressure, only the residual plastic deformation occurs due to matrix cracks. Thus, the noticeable fact of increasing the strain rate is manifested by changes in the damage modes. The samples tested were damaged by fiber twisting at low strain rates, and by interfacial debonding and delamination at high strain rates.

**Figure 8.** Strain rate versus time curves of the sample at different aging times and pressure under in-plane loading. (**a**) 1 bar. (**b**) 2 bar. (**c**) 3 bar. (**d**) 4 bar.

According to our knowledge, there is no empirical study proposing a constitutive model that reflects the impact of strain rate on bonded joints under dynamic conditions with the effect of aging due to the limited availability of reliable experimental data of dynamic tests. However, the effect of strain rate on the damage of metals and composites has been investigated in many works [42,43]. Many research works have recently been proposed to model the dynamic behavior of the material subjected to impact testing [44,45].

In this paper, empirical laws have been proposed relating to the dependence of the dynamic properties of bonded joints to the impact pressure and strain rate under dynamic in-plane compression loading for different aging times. To optimize the graphical representations, the deformation rate and the maximum stress are given in relation to the

impact pressure. Figure 9 shows the trends obtained from the strain rate dependencies of impact pressure and aging time. The evolution of the strain rate is approximated by a quadratic equation and presents two phases. In the first phase (1 bar ≤ P ≤ 2 bar) the increase is less pronounced and in the second phase the increase is accelerated. On the other hand, we can see that the stiffness of the material is sensitive to the deformation rate and the aging time. It is to be noticed that these nonlinear equations were derived from empirical plots of the dynamic characterization of specimens for different aging times under dynamic compression. The obtained equations are similar to the quadratic equations of glass/epoxy composite in dynamic compression [42,43]. In addition, Tarfaoui et al. [46] employed the non-linear evaluation of damage versus strain rate for dynamic modeling under pressure impact.

**Figure 9.** Variation of maximum stress/strain rate in function of aging time. (**a**) *t* = 0 h. (**b**) *t* = 216 h. (**c**) *t* = 648 h. (**d**) *t* = 864 h.

### *3.3. Theoretical Characterization of Absorbed Energy*

The incident impact energy is the complete energy available at the start and represents the kinetic energy contributed by the impactor. At the interface between the bar and sample, a portion of this energy is absorbed by the specimen and can produce damage or plastic deformation in various forms and also can generate heat, which is associated with the occurrence of microscopic/macroscopic damage. The rest of the energy consists of reflected and transmitted energy and can be measured from the deformation profile. The absorbed energy *Wabs* is expressed as follows:

$$\mathcal{W}\_{\text{abs}} = \mathcal{W}\_{\text{inc}} - \left(\mathcal{W}\_{\text{trans}} - \mathcal{W}\_{\text{ref}}\right) \tag{4}$$

With *Winc*, *Wref* and *Wtrans* representing the incident, reflected and transmitted energy, respectively, and their expressions are given by:

$$\mathcal{W}\_{\text{inc}} = \frac{A}{\rho c} \int\_0^t \sigma\_i^2(t) dt \tag{5}$$

$$\mathcal{W}\_{ref} = \frac{A}{\rho \mathcal{C}} \int\_0^t \sigma\_r^2(t) dt \tag{6}$$

$$\mathcal{W}\_{trans} = \frac{A}{\rho c} \int\_0^t \sigma\_t^2(t) dt \tag{7}$$

Figure 10 presents typical absorbed energy by the sample at an impact pressure of 4 bar for an aging time of *t* = 648 h. The energy versus time curves representative of the aged bonded joint sample indicates that the incident and reflected energy increases quickly compared to the transmitted and absorbed energy during the propagation of the deformation wave and their values are maintained stable after achieving a required value, and the apparition of the maximum absorption energy indicating the presence of damage to the joint. Much research indicates that there is a relationship between damage modes and the characteristics of the absorbed energy [47]. The results demonstrated that the impact energy can influence the absorbed energy and can be described as follows:


**Figure 10.** Typical profile of absorbed energy *Wabs*.

Moreover, the incident energy in the initial part increases quickly with the stress wave propagation, whereas the absorbed energy of the bonded joint has a small increase with the propagation of the stress wave. This is because the majority of the incident energy is absorbed by the specimens and just a minor fraction of the incident energy is transferred into the transmitted bar. In Figure 10, the presence of different zones can be identified:


In this way, the absorbed energy is decomposed into an unrecovered elastic part and an inelastic part, and can be expressed by the equation below:

$$\mathcal{W}\_{\text{abs}} = \mathcal{W}\_{\text{clas}} + \mathcal{W}\_{\text{diss}} \tag{8}$$

*Wabs* is the energy absorbed by the sample other than the stored deformation energy that is exhausted to create damage. The inelastic part is related to the damage generation in different forms, such as matrix fracture, fiber rupture, and delamination or decohesion between the fiber and matrix. The curve of the absorbed energy tends towards a fixed value at the end of the first cycle where the transmitted signal is reflected by the incident, and permanently corresponds to the energy dissipated by the damage in the sample. Figure 11 shows the absorption energy for a compression test for an aging time of 216 h at different impact pressures. The fluctuating profiles represent the release and storage of deformation energy throughout the experimental process. It can be easily observed that at low impact energy, a considerable amount of the input energy of *Welas* is stored in the loading phase and released in the discharge phase. On the other hand, *Wabs* has an increase as the impact pressure is raised, which in turn influences the release of *Welas* through the discharge phase. The last one (*Welas*) is the largest fraction of *Wabs* at low impact energy, while *Wdiss* takes a larger fraction of *Wabs* as the impact energy increases. It should be noted that *Welas* is zero when the sample has macroscopic damage. One can also note that the loading portions of the curves for the several impact pressures are coherent and overlapped. Additionally, an increase is observed in slope with the increase of impact pressure, with the suggestion that the reaction is conditioned by the deformation rate.

**Figure 11.** Absorbed energy for an aging time *t* = 216 h.

### *3.4. Damage Mechanisms and Failure Modes*

A camera with a high-speed system is employed to visualize the damage mode of the samples. Figure 12 shows the damage pattern of the bonded composite loaded in an in-plane direction. It can be noted that the damage is observed as cracks across the diagonal of the cube. In addition, matrix fractures and intense delamination are produced in the preferred inter-laminar planes. To understand the damage behavior of the bonded joined under dynamic loading, a camera of high-speed is used to capture the loading of the sample at different impact pressures. The damage is observed only at high impact pressures. For low impact pressure, only residual plastic deformations due to matrix fractures were present. So as the damage mechanism is developed, it is expressed in the form of matrix/fiber failure, fiber pullout, and delamination of the ply pack.

**Figure 12.** Optical micrographs of bonded composite under dynamic compression for *t* = 216 h.

The fiber and matrix can be observed to be highly bonded before hygrothermal aging. The beginning of the damage includes a "V" shaped shear band and the development of delamination at its tip. Microcracks were observed at the interface of matrix/fiber after 216 h, as shown in Figure 13. By increasing the hygrothermal aging time, the interface disbonding grows and propagates as long cracks at the 864 h aging time. The obvious interfacial degradation effects of moisture absorbed in the bonded composites on the interfacial area between the adhesive and the composite are detected. In these results, the hygrothermal influence causes chemical changes in the resin matrix, plasticization and dimensional changes, such as swelling, and deterioration of the strength of the adhesive/composite interface. The bond between the fiber and matrix weakens and the interface weakens as the moisture absorption increases. All these factors reduce the impact resistance of the composite.

**Figure 13.** Optical micrographs of adhesively bonded composite loaded through at pressure of 3 bar.

### *3.5. Microstructural Analysis of Adhesively Bonded Composite*

Figure 14 illustrates the microstructure of adhesive bonded composite as remarked using SEM. There were no significant morphological differences. The materials contain the elements of C, O, Mg, Al, Si, K and Ca. It is observed that the crack propagation rate of bonded composite joints becomes more prominent with the increase in strain rate. The results of reaching such raised strength and initial modulus can be justified by the fact that the composite materials become comparatively brittle and harder, which leads to high strength and modulus [36]. The results show also that the failure stresses degrade for various aging times for the samples subjected to in-plane compression because temperature and humidity considerably modify and weaken the composite. The strain rate is seen to be increased for various aging times. These results can be explained by the plasticization of the matrix being the principal factor in high strain rate properties when materials are not totally immersed in corrosion enclosure. Therefore, the matrix is more responsive to the strain rate and thus, plays a more important role than the fiber in identifying the material performance response to the high strain rate of polymer matrix composites. From the above discussion, it is further revealed that the crack propagation model changes as the deformation rate increases, and the higher the deformation rate, the more complex the crack propagation model appears.

EDS spectra of damage (spectrum 1) and undamaged (spectrum 2) adhesive bonded composite

Virgin sample

Damaged sample at 3 bar

Damaged sample at 4 bar

**Figure 14.** SEM micrographs of the composite of bonded adhesive at 864 h.

### **4. Conclusions**

In the present paper, the influence of moisture diffusion on the dynamic compressive behavior of glass/polyester composite was investigated. The failure modes were viewed by SEM and scanning electron microscopy for different aged specimens. The results obtained show that the hygrothermal influence induces chemical changes in the resin matrix, plasticization and dimensional changes, such as swelling, and also the weakened adhesive/composite interfaces. In addition, the bond between the fiber and the matrix is weakened as the humidity uptake increases. All of these elements would affect the dynamic behavior and the impact strength of an assembly of composite materials. The conclusions from this investigation can be summarized as follows:


**Author Contributions:** O.L.: Investigation, Methodology, Writing—original draft; M.T.: Methodology, Investigation, Project administration, Conceptualization, Writing—original draft, Writing review & editing, Funding acquisition; M.R. Methodology, Investigation, Writing—original draft, Writing—review & editing; H.L., S.J.S., M.A.D., A.A.: Validation, Visualization; M.N.: Methodology, Investigation, Writing—original draft, Writing—review & editing, Conceptualization, Data curation, Supervision, Validation, Visualization. All authors have read and agreed to the published version of the manuscript.

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

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

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

### **References**

