Dissolution Kinetics of R-Glass Fibres: Influence of Water Acidity, Temperature, and Stress Corrosion

Glass fibres slowly degrade due to dissolution when exposed to water. Such environmental aging results in the deterioration of the mechanical properties. In structural offshore and marine applications, as well as in the wind energy sector, R-glass fibre composites are continuously exposed to water and humid environments for decades, with a typical design lifetime being around 25 years or more. During this lifetime, these materials are affected by various temperatures, acidity levels, and mechanical loads. A Dissolving Cylinder Zero-Order Kinetic (DCZOK) model was able to explain the long-term dissolution of R-glass fibres, considering the influence of the pH, temperature, and stress corrosion. The effects of these environmental conditions on the dissolution rate constants and activation energies of dissolution were obtained. Experimentally, dissolution was measured using High Resolution Inductively Coupled Plasma Mass Spectrometry (HR-ICP-MS). For stress corrosion, a custom rig was designed and used. The temperature showed an Arrhenius-type influence on the kinetics, increasing the rate of dissolution exponentially with increasing temperature. In comparison with neutral conditions, basic and acidic aqueous environments showed an increase in the dissolution rates, affecting the lifetime of glass fibres negatively. External loads also increased glass dissolution rates due to stress corrosion. The model was able to capture all of these effects.


Introduction
Glass fibres (GFs) are the most common fibrous reinforcement material used in fibre-reinforced composites (fibre-reinforced polymers; FRPs) [1]. The GFs possess desirable reinforcement material properties, such as high hardness, strength, and stiffness [2,3]. Various types of GFs exist, such as E, ECR, R, and S-glass, listed in the order of their increasing mechanical strength. Glass fibre-reinforced composites (GFRPs) are often used in structural applications in marine and offshore industries, as well as in the wind energy sector [4][5][6][7][8][9][10][11]. In these applications, GFRPs are continuously exposed to water and humid environments for decades, with a typical design lifetime being around 25 years or more [8]. When exposed to such environments, hydrolytic degradation of the glass material occurs [12,13], leading to a reduction of the mechanical strength and the initiation of corrosion-induced defects in fibres and composites [14]. The degradation of glass fibres due to environmental attacks can severely affect the performance of GFRPs [1]. Aqueous environments clearly have a negative effect on the fatigue properties of GFRPs and such environments act primarily to reduce the fibre strength [15]. The fact that even water may corrode glass fibres has been known for many years [15]. In the long-term, hydrolytic degradation occurs mainly via the glass dissolution mechanism [12]. Thus, understanding the dissolution behaviour of glass is necessary to predict the deterioration of the composite's physical and mechanical properties [16].
The long-term environmental degradation of glass materials has been previously studied, mainly with respect to nuclear waste applications [17,18]. Most of the existing works on glass dissolution have been performed with bulk silicate glass and fibres are not studied often [12]. Recently, there has been an increasing interest in the environmental aging of FRPs, where R-glass and E-glass are often used as reinforcement [4,12,16,[19][20][21]. Even so, very few studies exist on the kinetics of GF dissolution (mostly on E-glass) [22,23] and even less on R-glass [12,16].
The degradation of GFs in water can be divided into two kinetic regions (Phases I and II). Phase I is the short-term non-steady-state. Phase II is the long-term steady-state [12,16]. Phase I is complex and involves such processes as ion exchange, gel formation, and dissolution [12,17]. However, once the long-term kinetic region is reached, glass dissolution becomes a dominant process [12,17,24]. For the studied GFs, this transition to steady-state occurs in about a week (166 h) [12,16]. This time is very negligible compared to the typical design lifetime of about 25 years or more for the structures in which these materials are implemented [8]. Thus, as in the case of structural applications, when the long-term degradation is mainly of interest, dissolution is the governing process, which has to be understood and well-predictable.
During glass-water interactions, several chemical reactions may occur, shown in the chemical reactions (1)- (11) [12,17,24,25]: MeCl x The chemical reaction (10) can also be written as a combination of the subsequent reactions (12) and (13), meaning that initially H 2 SiO 3 is formed, which dissociates weakly and further reacts with water to form silicic acid: The elements that are released during the degradation of R-glass are Na, K, Ca, Mg, Fe, Al, Si, and Cl [12]. The cumulative mass loss is the summary mass of all of the ions released [12]. The total mass loss is what causes the radius reduction [12]. As shown in the chemical reactions (1)- (11), various competing reactions happen in parallel. During Phase I, these reactions happen at independent rates, then one process later becomes limiting and dominates the behaviour (Phase II). Thus, the degradation process should be divided into two aforementioned stages. Furthermore, the sub-processes may follow individual kinetics and may be interdependent during Phase I [12]. During Phase I, the process is rather complex. It involves competing leaching mechanisms in addition to the dissolution reactions, i.e., ion exchange and gel formation. Therefore, it cannot be fully described with just one kinetic equation [12,17,24]. The complexity and varying rates of the competing processes during Phase I turn this stage of the glass degradation into a non-linear and non-steady-state process.
During Phase II, however, dissolution becomes dominant [17,24] and the process becomes steady-state, governed by the dissolution kinetics [12,16]. For R-glass, the contribution from Si to the total mass loss is the largest (56.1% by mass) and governs the long-term steady-state dissolution [12]. Once the steady-state is reached, the reactions proceed at constant independent rates [12]. It may be speculated that, during Phase II, an equilibrium composition of glass (likely different from the bulk composition) is obtained in the outer layers in contact with water. This allows elements to dissolve at a limiting rate and proportionally to their content in the outer layers of the glass, being indicative of the glass outer layer composition [12]. In other words, the long-term dissolution is dominated by the slowest dissolving element and the other elements following proportionally, at some equilibrium composition [12].
The non-linear non-steady-state stage of degradation can be approximated using the Dissolving Cylinder Zero-Order Kinetic (DCZOK) model, when the exact shape of the non-linear short-term part is not important, i.e., when the long-term degradation is of the most interest, such as in cases of the long service life of structures implementing GFs [8,12].
Glass dissolution kinetics depend on the surface area of the glass in contact with water. This surface area is proportional to the glass fibre radius. As the dissolution proceeds, the radius decreases and results in the deceleration of the mass loss; the DCZOK model accounts for this effect [12,16]. Glass dissolution kinetics are affected by the pH, temperature, glass composition (i.e., E, ECR, R, S-glass), solution composition (distilled water or seawater), surface area of the glass, the protective effect of the sizing, and the ionic strength of the solution, to a lesser extent. The temperature and pH are considered to be key parameters for the long-term glass dissolution [16,17,20,24,26,27]. Elevated temperatures accelerate the dissolution since the dissolution is exothermic [17,20,26]. The activity of the water is dependent on the composition of the aqueous environment, notably, on the, as well as ionic concentration of some species in the glass composition [28]. Interestingly enough, for R-glass, the sizing was found to slow down the dissolution kinetics of the GFs by almost an order of magnitude [16]. In addition, external loads may have an influence on glass dissolution kinetics, due to stress corrosion [27,28].
Humid and water environments reduce the strength of silicate glasses due to the growth of flaws, such as surface cracks, under the combination of stress and chemical attacks, known as stress corrosion [1,15]. The glass corrosion (dissolution) rate is accelerated by this stress [1].
The classical theory of stress corrosion involves the chemical reaction of a water molecule with silica, and was described elsewhere [27]. A model and the details of the chemical interaction between the adsorbed water and the Si-O bond was described by Michalske and Freiman [29]. Their model showed that water is efficient in the stress-activated corrosion of glass [27]. The basic mechanism of the stress corrosion reaction was related to the stress-enhanced thermal activation of a dissociative hydrolysis reaction, presented in Figure 1 [28,29]. Any molecule that is able to donate a proton (hydrogen ion and has a lone pair of electrons, such as water, can react with the strained Si-O bond at the tip of a crack, provided that the molecular diameter is smaller than 0.5 nm [30]. The first property is necessary for the harmonized reaction of adsorption and the chemical scission of the siloxane bridges of the glass network. The small diameter is necessary to avoid steric problems, so that the reactive water molecules can reach the strained bonds at the crack tip [28]. The constrained Si-O bonds have an enhanced chemical reactivity with water due to the intensity of the strain, as was shown via hydrolysis experiments for different siloxane ring structures, after [31]. Michalske and Bunker [32] also supported the effect of stress on the stress corrosion reaction rate. They proved this hypothesis by conducting molecular-orbital simulations of the stress corrosion reaction on distorted siloxane bonds. Several other reactions are possible for other glass compositions, but they all share the feature of being stress-enhanced and thermally activated, therefore explaining the exponential (Arrhenius type) relationship [28]. In addition, it is believed that the interaction of water with cations in the flaws leads to hydrolysis and a local increase in the concentration of hydroxyl ions, leading to an increase in the pH locally [33].
Various approaches have been proposed on how to model the kinetics of glass degradation. Grambow has developed a general rate equation for the reaction of glass in nuclear wastewater solutions with different pH levels, surface areas, times, temperatures, and compositions of the reagents [18]. Delage et al. described the dissolution of borosilicate glass in aqueous media using a combination of surface reaction kinetics and diffusion [34]. Mišíková et al. have studied the kinetics of the leaching of E-glass fibres in distilled water at various temperatures, linking the glass fibre dissolution to the temperature using statistical methods [22]. Tournié et al. presented an extensive experimental study of the glass corrosion mechanism in both acidic and basic solutions. They reported that a strong acid attack was dependent on the composition of the glass [13]. Ma et al. presented a model for the glass dissolution, based on the combination of the chemical affinity and diffusion [35]. Bashir et al. studied the kinetics of the dissolution of E-glass fibres in hot alkaline solutions using zero-order and shrinking cylinder models [23,36]. They reported that both an increased temperature and increased alkalinity (higher pH) accelerated the dissolution reaction rates [23]. Bashir et al. concluded that the rate-limiting step was either the glass dissolution itself or the diffusion of hydroxide ions [23]. Krauklis and Echtermeyer developed and presented a Dissolving Cylinder Zero-Order Kinetic (DCZOK) model that explained the long-term dissolution of R-glass in water for both fibres [12,16] and GFRPs [16]. The DCZOK model is able to predict both the mass loss and the radius reduction kinetics due to dissolution [12].
Glasses formulated with significant amounts of various metal oxides, such as R-glass, have rarely been studied, but they are becoming popular for structural applications [12,18]. Thus, this work focuses on the dissolution kinetics of R-glass in water at various environmental conditions (pH, temperature, stress corrosion). However, the model should be applicable to other types of glass as well, as SiO 2 is the major component in virtually all types of glass [23].
The aim of this work is to study the effects of the pH, temperature, and stress corrosion on the R-glass fibre dissolution kinetics, in terms of the rate constants and activation energies, and to capture these effects using the Dissolving Cylinder Zero-Order Kinetics (DCZOK) model.

Materials
A typical glass fibre used for marine, oil, and gas applications was selected for this study. Boron-free and fluorine-free high strength, high modulus 3B HiPer-Tex™ W2020 R-glass fibre bundles were used. These were classified as high-strength, high modulus R-glass (defined by an international standard [37]). The average fibre diameter was 17 ± 2 µm [12,38]. The density of the glass (ρ f ) was 2.54 g/cm 3 [12,38]. The authors estimated that a bundle had about 4098 fibres [12]. The specific surface area of the glass fibres was determined to be 0.09 m 2 /g from the geometrical considerations, as a product of the number, circumference, and length of the fibres [12]. All of the fibres used during this work were sized, whereas the dissolution of unsized R-glass fibres was previously studied in another study by Krauklis and Echtermeyer [16]. A micrograph of R-glass fibre bundles is shown in Figure 2. Distilled water (0.5-1.0 MΩ·cm) was used for the conditioning of the GF bundles, produced using a water purification system, the Aquatron A4000 (Cole-Parmer, Vernon Hills, IL, USA). The pH of the distilled water was 5.650 ± 0.010, being lower than neutral due to the dissolved CO 2 from the atmosphere in equilibrium. IUPAC standard buffer solutions (Radiometer analytical, Lyon, France) were used for studying the effect of the pH on the kinetics of the GF dissolution. The solutions of pH 1.679 ± 0.010, 4.005 ± 0.010, 5.650 ± 0.010, 7.000 ± 0.010, and 10.012 ± 0.010 were used.

Glass Dissolution Experiments:
Aqueous dissolution experiments of GF samples were performed using a batch system. The samples for the dissolution study were weighed using AG204 analytical scales (±0.1 mg; Mettler Toledo, OH, USA) before the experiments. The samples were placed in inert closed vessels filled with 50 mL of distilled water or pH buffer solutions. The tight sealing of samples was ensured. The water-tight vessels, with the samples and water solutions in them, were placed into the water bath.
The water's temperature (25,40,60, 80 • C) in the bath was controlled via proportional-integralderivative (PID)-controlled heating, giving an accuracy of ±1 • C. A two-stage heating system was used in order to make sure that there is no contact of the sample water with potential contaminating ion release sources, i.e., the heating element.
pH measurements were performed using a standard MeterLab PHM210 pH-meter (Radiometer analytical, France) calibrated with IUPAC standard buffer solutions, providing an accuracy of ±0.010 pH.
The concentration of the dissolved ions in the water from the vessels was analysed over time via high resolution inductively coupled plasma mass spectrometry (HR-ICP-MS), providing glass dissolution kinetics. Experimentally, the total mass loss of glass material was measured as a sum of all of the ions' release, quantified with HR-ICP-MS. Analyses were performed using a double focusing magnetic sector field ICP-MS Finnigan ELEMENT 2 (Thermo-Scientific, Waltham, MA, USA), equipped with a PrepFAST sample introduction system (ESI/Elemental Scientific) and an UltraClave pre-treatment/digestion machine (Milestone, Brøndby Kommune, Denmark). The acidification of samples was performed using an ultra-pure grade HNO 3 SubPur (Milestone) in order to avoid the adsorption of ions into the wall of the sample vials that would otherwise give erroneous results.
The benefit of HR-ICP-MS versus gravimetric analysis is that it allows the measurement of the dissolution kinetics of each separate ion, as well as the total mass loss [12]. The data obtained from the HR-ICP-MS experiments are in the form of the non-cumulative mass concentration at each point in time (g/L) and need to be converted to the cumulative m dissolved form by using Equation (14): where V water is the volume of a water sample using the HR-ICP-MS measurement. The V water used for the experiments was 50 mL. Equation (14) is valid for each individual ion release and for the total mass loss.
For studying the influence of the external loads on the dissolution kinetics, a stress corrosion rig was designed, built, and used during the experiments. The experimental rig is shown in Figure 3 and the design schematics of the stress corrosion setup are attached in Appendix A. The principle is that glass fibres are inside an inert cylinder pushed by an inert rod, which transfers the stress from the weights to the fibres. The water samples for the HR-ICP-MS analyses were taken from the main cylinder, where the distilled water was in contact with the stressed glass fibres. The temperature was PID-controlled by using a two-stage heating system, ensuring a temperature of 60 ± 1 • C. Polymeric spheres, seen in Figure 3, were used to reduce the evaporation rate of the heating water.

The Model and Its Assumptions
An approach used to describe the glass dissolution kinetics in the fibres was proposed and explained in works by Krauklis and Echtermeyer [12,16]. This approach is used as a basis for this study. In this work, the kinetic constants are obtained for the Dissolving Cylinder Zero-Order Kinetic (DCZOK) model for various pH and temperature conditions, as well as for different stress levels. The protective effect of sizing is accounted for in this model [16].
The rate of the dissolution is dependent on the glass dissolution rate constant (K 0 ), the glass surface area exposed to water (S), and the effect of sizing (ξ sizing ) [12]. For fibres, in infinite water availability conditions, the dissolution, being a surface reaction, can be well-described with zero-order kinetics [23], but the decrease in the fibre radius, and thus in the surface area over time, should be accounted for [12]. For sizeless fibres, Equation (15) should be used [16]: For sized fibres, the effect of sizing, ξ sizing , should also be accounted for, as shown in Equation (16) [16]: where m is the total cumulative mass dissolved after time t, K 0 is a material/environment interaction property, and ξ sizing is the protective effect of sizing. K 0 is affected by environmental effects [12].
In this Equation (16), K 0 and ξ sizing are the time-independent parameters, whereas S changes with time. As the dissolution continues, the fibre's radius is reduced and the total surface area (S) is decreased. This leads to a decrease in the dissolution rate. The ion release rate slows down linearly with a decrease in the fibre's radius [12]. The glass fibre's radius reduction is linear with time. The proportionality is given by the dissolution rate constant (K 0 (g/(m 2 ·s))) and the density of the glass (ρ glass (g/m 3 )) [12]. A linear radius decrease of the glass fibres over time was previously experimentally observed elsewhere [23] and is accounted for in the DCZOK model [12].
The effect of sizing on the glass dissolution (ξ sizing ) for the studied R-glass is 0.165, protecting the fibres from water by almost an order of magnitude, in terms of dissolution rates. The effect of sizing (ξ sizing ) is assumed to be time-independent [16]. ξ sizing goes into the model linearly, since fibres are covered in "islands" or patches of sizing, meaning not the whole surface is coated [39]. These areas of the bare unsized glass dominate the dissolution behaviour (the difference being roughly an order of magnitude) [16].
Taking into account the geometry of the fibres, the mass loss kinetic DCZOK model equation in its differential form is the following, as shown in Equation (17) [12]: where n is the number of fibres (-), l is the length of the fibres (m), r 0 is the initial fibre radius (m), and ρ glass is the density of the glass (g/m 3 ). The model shows the rate of mass loss due to dissolution. It can be used for the total mass loss ∂m ∂t , as well as for Si ∂m Si ∂t or another ion's release. As can be seen from the differential expression of the DCZOK model in Equation (17), the radius reduction over time is accounted for in the model. The pH, temperature, and stress corrosion affect the material-environment energy-activated interactions, thus affecting the dissolution rate constants (K 0 ), as represented by Equation (18).
where pH is the acidity of the environment (-), T is the temperature (K), and σ is the stress (MPa). The model involves the following assumptions: the model is deterministic and all of the glass fibres are assumed to have the same initial fibre radius (r 0 ); the cross-sectional surface area at the ends of the fibres is assumed negligible in the surface area calculations; the length of the fibres (l) is assumed to be constant during the whole dissolution process. During the whole degradation process, the density of the glass material is assumed to be constant (ρ glass ). The effect of sizing (ξ sizing ) is assumed to be independent of the environmental conditions and time [12,16]. For free fibre bundles (not embedded in the composite), the condition of the infinite availability of water is ensured by using large volumes of water, thus making the rate of the reaction independent of the water concentration [12,16].

The Effect of the Environment on the Dissolution Rate Constant
It has been reported that the temperature dependence of the rate of the glass fibre radius reduction follows an exponential trend and the diameter reduction is linearly related to the loss of material in fibres due to dissolution [12,23]. The fact that the temperature dependence of a dissolution rate constant is described by the Arrhenius equation is also confirmed by the evidence reported further in this manuscript. The temperature dependence of a dissolution rate constant can be described using the Arrhenius equation (Equation (19)), being an exponential function: where A is the pre-exponential factor (g/(m 2 ·s)), R is the gas constant being 8.314 J/(mol·K), T is the absolute temperature (K), and E A is the activation energy (J/mol). Both the pH and stress corrosion are thought to affect the activation energy term in the Arrhenius equation [31,40].
The activation energy can be obtained through the graphing of the constant pH and σ, using Equation (20) [18]: The rate constant can be calculated for each temperature from the general differential mass loss DCZOK model equation. To obtain the E A , the lnK 0 is plotted against 1 T . The pre-exponential factor (A) is assumed to be constant. The DCZOK model can be expanded (Equation (21)):

Modelling of the Two Distinct Stages of Glass Degradation
Mass loss curves made during the aging of glass over time can be divided into two parts: Short-term non-steady-state degradation and long-term steady-state dissolution [12]. When the long-term dissolution (Phase II) is of the most interest, the shape of the non-linear short-term stage (Phase I) is not important. In such case, an approximate solution can then be used to describe the Phase I of degradation as well, using the DCZOK model. In the case when the short-term degradation (Phase i) is of the most interest, the kinetic models for each of the individual sub-processes, i.e., ion exchange, gel formation, and dissolution, should be implemented.
The model can also predict the fibre radius reduction kinetics during the dissolution, using the same four parameters that can be determined experimentally [12]. The DCZOK model can be used to describe the dissolution for each element separately, as well as for the total mass loss. The glass mass loss (all elements included) is described by four parameters (r 0 , t st , K I 0 , K I I 0 ). For each ion, the individual K I 0 and K I I 0 values have to be obtained, while the t st and r 0 for each individual ion release are the same and are equivalent to that for the total mass loss [12]. K I 0 and K I I 0 are the dissolution rate constants (g/(m 2 ·s)) for the Phase I and Phase II regions, respectively.
For each material, the activation energy function (E A (pH, σ)) of the pH and stress has to be established experimentally for both the short-term and long-term stages of glass degradation.
The total material loss and the release of Si at various environmental conditions were modelled using the DCZOK model Equation (21) and the rate constants were obtained and reported.

Results
From the measured cumulative concentration data of the ions released during the dissolution of glass, the dissolution rate constants were obtained using non-linear regression for Si (K I 0 Si and K I I 0 Si ) and for the total mass loss (K I 0 total and K I I 0 total ). The Generalized Reduced Gradient (GRG) non-linear regression method was used. It involved minimization of the sum of the squares of the differences between the modelled values and the experimental values. The obtained parameters are systematized in Tables 1-3 for the effect of the pH, temperature, and stress, respectively. The values obtained are comparable with the dissolution rates of other glass fibres, reported in the literature [12,16,41].
In Figure 4, and further also for the glass dissolution kinetics at various stress levels and temperatures, the points represent the experimental data, whereas the lines represent the DCZOK modelled curves using the rate constants reported in Tables 1-3, respectively. In figures of the glass dissolution kinetics at various environmental conditions, the glass mass loss is reported as normalized per the initial surface area of the fibres (S 0 ), as is common for surface reaction kinetics [12,16,42]. The effect of sizing on the glass dissolution (ξ sizing ) for the studied R-glass was 0.165 [16].

Effect of pH on the Glass Dissolution Kinetics
The dissolution rate constants for the glass dissolution from fibres in water at various pH levels at 60 ± 1 • C are shown in Table 1. This steady-state was achieved after about a week. Table 1. The glass dissolution rate constants obtained via the regression of the experimental data of the R-glass fibre bundles using the Dissolving Cylinder Zero-Order Kinetic (DCZOK) model for Si and the total mass loss at 60 • C and various pH levels. K I 0 Si , K II 0 Si , K I 0 total and K II 0 total are dissolution rate constants for Phase I and II for Si and total glass dissolution, as indicated by indices. pH K I 0Si (g/(m 2 ·s)) K II 0Si (g/(m 2 ·s)) K I 0total (g/(m 2 ·s)) K II 0total (g/(m 2 ·s)) Reference The experimental glass dissolution data and the modelled DCZOK dissolution curves for various pH levels are shown in Figure 4.
It is possible to obtain a pH influence on the activation energies of dissolution by rearranging Equation (19) into the following form (Equation (22)): A pH function of the activation energy of dissolution E I I A = f (pH) is shown in Figure 5. It is fine to approximate the kinetic pH dependency as a polynomial function, being a parabolic function in this case [43][44][45][46].  Both the glass and the Si dissolution rate constant trends are similar, with an exception at low pH levels, i.e., at a pH of 1, the rates of Si dissolution increase more dramatically than the total glass dissolution rates, indicating a higher contribution of Si to the material loss during the dissolution in strongly acidic environments. Thus, GFRPs should not be used in strongly acidic conditions. This is in agreement with an observation in another study, stating that many GFRPs fail catastrophically after a critical time when exposed to acids [1].

The Effect of Temperature on Glass Dissolution Kinetics
The dissolution rate constants for glass dissolution from the fibres in water at various temperatures at a pH of 5.650 ± 0.010 are shown in Table 2. The steady-state was achieved after about a week. Table 2. The glass dissolution rate constants obtained via the regression of the experimental data using the Dissolving Cylinder Zero-Order Kinetic (DCZOK) model for Si and the total mass loss at a pH of 5.65 and at various temperatures. K I 0 Si , K II 0 Si , K I 0 total and K II 0 total are dissolution rate constants for Phase I and II for Si and total glass dissolution, as indicated by indices.
Temperature ( • C) K I 0Si (g/(m 2 ·s)) K II 0Si (g/(m 2 ·s)) K I 0total (g/(m 2 ·s)) K II 0total (g/(m 2 ·s)) Reference The Arrhenius approach was used to obtain the activation energy of the steady-state Si and glass dissolution. The activation energy was obtained through graphing at a constant pH and σ, similar to what was done in another study, using Equation (23) [23]: The graphing approach is shown in Figure 6. The obtained activation energy (E A ) of Si dissolution is 53.46 kJ/mol (using K I I 0 Si values). The pre-exponential factor (A) for Si dissolution is 6.82·10 −1 g/(m 2 ·s). The obtained activation energy (E A ) of the total glass dissolution is 34.84 kJ/mol (using K I I 0 total values). The pre-exponential factor for the glass dissolution (A) is 1.67·10 −3 g/(m 2 ·s). The obtained values are consistent with the values reported in the literature, being slightly lower than for E-glass (58-79 kJ/mol in alkaline solutions [23]). The activation energy for the Si dissolution of various silica-based materials may be around 60 kJ/mol [47].

The Effect of Stress Corrosion on Glass Dissolution Kinetics
The dissolution rate constants for glass dissolution from fibres in water loaded with various weights (0.03, 5, 8, and 10 kg) at a pH of 5.650 ± 0.010 and at 60 ± 1 • C are shown in Table 3. The stress was calculated as half of the weight divided by the cross-sectional area of the glass fibre bundles (0.93 mm 2 ). The division by two is due to the fact that the weight is carried by a bundle which goes from each side of the rig, meaning that the total area is that of two bundles. In order to visualize the setup and the fibre arrangement, the experimental rig in action is shown in Figure 3. The steady-state was achieved after about 5 days (120 h), sooner than in all of the unstressed cases, indicating that stress may slightly accelerate the transition towards the steady-state dissolution. Table 3. The glass dissolution rate constants obtained via the regression of the experimental data using the Dissolving Cylinder Zero-Order Kinetic (DCZOK) model for Si and the total mass loss at 60 • C and at a pH of 5.65 at various stress levels. σ stands for stress; K I 0 Si , K II 0 Si , K I 0 total and K II 0 total are dissolution rate constants for Phase I and II for Si and total glass dissolution, as indicated by indices.

Weight (kg)
σ K I 0Si (g/(m 2 ·s)) K II 0Si (g/(m 2 ·s)) K I 0total (g/(m 2 ·s)) K II 0total (g/(m 2 ·s)) Reference The experimental glass dissolution data and the modelled DCZOK dissolution curves for various stress levels are shown in Figure 8. A stress function of the activation energy of dissolution E I I A = f (σ) for R-glass fibre bundles is shown in Figure 9. An increase in the stress seems to accelerate glass dissolution, showing a similar trend for both Si and the total mass loss kinetics. The increase in the stress reduces the activation energy of dissolution linearly and accelerates the glass dissolution rates exponentially. Additional studies involving higher stresses are recommended to establish this influence in more detail. Figure 9. The activation energy of dissolution (yellow) and the dissolution rate constants (red), functions of the stress (σ) for the glass fibre bundles at the steady-state. Determination coefficients R 2 are shown. 4E-9 stands for "4 × 10 −9 ".

Discussion
R-glass fibres were studied in this work. However, there should be no limitations to applying the DCZOK model to other types of glass fibres. Validation of the model with various types of glass fibres, such as E, ECR, and S, is advised. The processes for all of these fibres are expected to follow similar trends, whereas only the values of the terms in the model have to be assessed for each material individually.
The industry is interested in the pH, temperature, and stress corrosion for environmental effects. Temperatures are also interesting to the industry for accelerated testing purposes. As the industry is concerned with lowering the testing time for the fibre-dominated property deterioration in GFRPs, the model allows the prediction of the loss of glass material. While it takes time to obtain the parameters, the most significant time saving comes from using the kinetic constants and the model.
The influence of the pH, temperature, and stress corrosion on the rate constants and dissolution activation energies was obtained, however, the influence of each parameter was studied once at a time. Thus, a suggestion for future work includes a cross-parametric study of the pH, temperature, and stress corrosion, in order to deduce a general analytical solution for the environmental influence on the activation energy of dissolution and to study whether there is a coupled effect. The DCZOK model should be extendable to include the effect of the ionic strength and, finally, to be further used in linking the reduction in the mechanical strength of the fibres with the dissolution kinetics.
The authors think that validating the DCZOK model for seawater conditions (about 1.84-12.62 mg SiO 2 /kg water; pH of seawater at 7.8) would be highly beneficial, especially for the marine and offshore industries, since the real-life structures most often operate in the seawater environment. When GFs are used in seawater, the dissolution of glass occurs at a slower rate, due to the presence of silica in the seawater, which is in seawater from the contact with sand and minerals [48][49][50][51]. The approach in distilled water is conservative with regards to seawater, meaning that structures implementing glass fibres designed for distilled water conditions should not encounter penalties to their service time in the seawater.
The model should be applicable to other types of glass as well, as SiO 2 is the major component in virtually all types of glass [23], but it would be beneficial to validate this model experimentally with other types of glass fibres. Information on the exact composition of the studied R-glass was not available from the manufacturer, supplier, or technical data sheets. However, it may be possible to speculate that the composition of the dissolving ions during steady-state Phase II is somewhat representative of the glass composition, since all ions follow the slowest dissolving element, as was reported elsewhere for the same material [12]. The ion composition during the steady-state was the following: 2.7 wt% Na, 1.4 wt% K, 14.3 wt% Ca, 15.1 wt% Mg, 0.1 wt% Fe, 6.9 wt% Al, 56.1 wt% Si, and 3.4 wt% Cl [12].
Another aspect that should be investigated is the effect of the sizing on the glass fibre dissolution, i.e., its chemical composition, amount, and distribution on the fibres. For the studied R-glass, it is known that the sizing is epoxy-compatible and its loss on ignition (LOI) value is 0.64 wt%. The protective effect of sizing for the studied R-glass fibres was reported in another study [16].
It was previously observed by Scheffler et al. [52] that strong alkaline solutions can significantly damage glass fibres due to the chemical attack [52]. This observation is consistent with the results in this work, indicating much higher dissolution rates at a pH of 10 (at least an order of magnitude higher) compared to the neutral conditions.

Conclusions
The analytical model, termed the Dissolving Cylinder Zero-Order Kinetic model, was successfully used to explain the long-term glass dissolution experiments of R-glass fibres with various environmental conditions: The model accounts for the influence of the pH, temperature, and stress, as well as the effects of sizing protection (ξ sizing for the studied R-glass is 0.165). The glass dissolution rate constants were obtained and reported for various pH levels and temperatures, as well as for various stress corrosion conditions.
The temperature showed an Arrhenius-type influence on the kinetics, increasing the rate of dissolution exponentially with an increasing temperature. The activation energy of the steady-state glass dissolution was obtained and reported at a pH of 5.65 and with no stress (53.46 and 34.84 kJ/mol for Si and the total glass dissolution, respectively). The activation energy of dissolution is affected by the pH and stress. The activation energy decreases linearly as the stress increases. The influence of the pH is more complicated and may be described by a parabolic polynomial function; the activation energy peaks at a pH of 5.65 and decreases towards more basic, as well as more acidic, conditions.
In comparison with neutral conditions, basic and acidic aqueous environments showed an increase in the dissolution rates, affecting the lifetime of the glass fibres negatively.
External loads also increased the glass dissolution rates due to stress corrosion. The stress seems to accelerate the glass dissolution rates exponentially. The model was able to capture all of these effects. Acknowledgments: This work is part of the DNV GL-led Joint Industry Project "Affordable Composites" with 19 industrial partners and the Norwegian University of Science and Technology (NTNU). The authors would like to express their thanks for the financial support from The Research Council of Norway (Project 245606/E30 in the Petromaks 2 programme). The authors are thankful to Børge Holen, Carl-Magnus Midtbø, Cristian Torres Rodriguez, Melanie Shebel, Syverin Lierhagen, and Konstantīns Viligurs. Andrey is especially grateful to Oksana V. Golubova.

Conflicts of Interest:
The authors declare no conflicts of interest.

Abbreviations r
Fibre radius (m) r 0 Initial fibre radius (m) K 0 Dissolution rate constant (g/(m 2 ·s)) K I

Si
Dissolution rate constant for Si only (non-steady-state) (g/(m 2 ·s)) K II

Si
Dissolution rate constant for Si only (steady-state) (g/(m 2 ·s)) K I

total
Dissolution rate constant for glass (non-steady state) (g/(m 2 ·s)) K II 0 total Dissolution rate constant for glass (steady-state) (g/(m 2 ·s)) t st Time when long-term steady state is reached (s) n Number of fibres (-) l Length of fibres (m) S Glass fibre surface area (m 2 ) S 0 Initial glass fibre surface area (m 2 ) pH Acidity of the environment (-) A Pre-exponential factor (g/(m 2 ·s)) R Universal gas constant (8.314 J/(mol·K)) T Absolute temperature (K)