In-Depth Kinetic Modeling and Chemical Analysis for the Epoxidation of Vegetable Oils in a Liquid–Liquid–Solid System

: A heterogeneous catalyst for producing epoxidized vegetable oils, an important intermediate in the production of non-isocyanate polyurethanes, is essential for product separation and for decreasing the side-reaction, i.e., ring-opening reaction, via the Prileschajew method. The development of reliable kinetic models considering key variables for both phases and the mass transfer phenomena is missing in the literature. The reaction pathway for the ring-opening reaction is also under debate. Therefore, we studied the kinetics of epoxidation of cottonseed oil by perpropionic acid over the solid acid catalyst amberlite IR-120. An in-depth kinetic model was developed by using Bayesian inference. The reaction pathway for the ring opening was investigated. Propionic acid, a weak acid, allows for a decrease in the oxirane ring-opening side reaction.


Introduction
Valorization of biomass can aid the chemical industry to use less fossil raw materials, thus increasing its sustainability.Vegetable oils in the chemical industry are mainly used for biodiesel production via transesterification [1][2][3][4][5].The other promising use of valorizing vegetable oil is for producing polymers or materials [6][7][8].Among these polymers, one can cite the production of polyurethanes, which has attracted increasing attention in recent years [9].One way to produce non-isocyanate polyurethanes is via the aminolysis of carbonated vegetable oils [10,11].
One of the first steps in valorizing vegetable oils into polymers is the functionalization of these oils, and the epoxidation of the unsaturated groups is a good way to achieve this [12].From a previous review [12], we have discussed several possible ways of vegetable oil epoxidation: − Use of gaseous oxygen with the risk of over-combustion, but one needs a carrier molecule such as cumene [13]; − Direct use of hydrogen peroxide [14], but this technique was shown to work on fatty acid methyl ester and requires the use of a solvent, such as ethyl acetate; − Use of in situ-generated percarboxylic acid, and produce the epoxidized vegetable oil via the Prileschajew method [15][16][17], which is the most common technique.
The Prileschajew method produces in situ percarboxylic acid from the perhydrolysis of carboxylic acid in the aqueous phase.Then, this percarboxylic acid, also known as an oxygen carrier, diffuses in the organic phase to epoxidize the unsaturated groups.This method works well because the solubility of hydrogen peroxide in the organic phase is low.Usually, an acid catalyst is used to accelerate the perhydrolysis reaction.
The number of articles concerning the use of the Prileschajew method for the epoxidation of vegetable oils is enormous.Research groups have focused on the development of catalysts (i.e., Al 2 O 3, Amberlite, CH 3 ReO 3 , etc.) [18][19][20], the intensification of the liquid-liquid − Two-phase system [27,[34][35][36]: mass transfer between the aqueous and organic phases is included.In this model, the mass transfer coefficient of different species must be evaluated or calculated; − Pseudo-homogeneous system [30,[37][38][39][40][41]: concentrations of species in each phase are considered, and the mass transfer phenomenon is assumed to be faster than the chemical reaction.In this model, there is no need to consider mass transfer coefficients; − Homogeneous system [15,16,[42][43][44][45][46][47][48][49]: the liquid-liquid system simplified into a homogeneous one.
To the best of the authors' knowledge, kinetic models considering the mass transfer between organic and aqueous phases, and the non-linear regression by considering the variables of the two liquid phases have not been developed and validated.To fill this gap, this study proposed such a kinetic model for the epoxidation of cottonseed oil by perpropionic acid over amberlite IR-120 in an isothermal batch reactor.An in-depth analytical study was carried out to unravel the ring-opening pathway.The development of a reliable kinetic model considering liquid-liquid mass transfer is essential for the development of industrial processes.

Analysis, NMR and FTIR Results, and GPC
It is fundamental to determine the different species present during the epoxidation.This knowledge is vital for elaborating reaction pathways.For this, we used NMR, FTIR, and GPC analytical methods.

NMR Spectra
1 H NMR could be used to determine the molecule changes, and to confirm the reactions of the epoxidation and ring opening.As displayed in Figure 1, the corresponding signal of the carbon-carbon double bond (CH=CH) from 5.2 to 5.5 ppm disappeared with the presence of a signal of the epoxidized group at 3.00 ppm.

NMR Spectra
1 H NMR could be used to determine the molecule changes, and to confirm the reactions of the epoxidation and ring opening.As displayed in Figure 1, the corresponding signal of the carbon-carbon double bond (CH=CH) from 5.2 to 5.5 ppm disappeared with the presence of a signal of the epoxidized group at 3.00 ppm.The NMR spectra of intermediate samples taken during the reaction are shown in Figure 2.With the epoxidation programming, the peak intensity of the CH=CH double bond was clearly decreased, while the signal of the oxirane ring conversely appeared at 3.00 ppm.
The absence of peaks between 3.6 and 3.8 ppm in the 1 H NMR spectra signifies the absence of hydroxyl groups resulting from hydrolytic degradation of the triglyceride, which is confirmed by the integrations of the signals of the protons in the alpha position to the oxygen of the ester groups, compared to the integration of the methyl groups.The absence of these peaks also signifies the absence of hydroxyl groups resulting from the epoxide ring opening.However, this does not mean that ring opening by water, leading to vicinal diols, did not take place, but that those elimination reactions giving access to new carbon-carbon double bonds could have taken place.This could explain the presence of a new broad peak between 5.4 and 5.6 ppm.On the other hand, the elimination reaction would also give access to enol functions, which can then be partially or totally converted into ketone functions, which are difficult to characterize from the 1 H NMR spectra of our products.The NMR spectra of intermediate samples taken during the reaction are shown in Figure 2.With the epoxidation programming, the peak intensity of the CH=CH double bond was clearly decreased, while the signal of the oxirane ring conversely appeared at 3.00 ppm.

FTIR Spectra
The FTIR spectra of CSO (cottonseed oil) and ECSO (epoxidized cottonseed oil) are displayed in Figure 3.The conversion of the CH=CH double bond in CSO as the epoxidation reaction proceeds could be verified by the disappearance of the corresponding band at about 3000-3050 cm −1 and the presence of a band at about 850 cm −1 , corresponding to C-O-C stretching vibration of the epoxide ring in CSO.The absence of peaks between 3.6 and 3.8 ppm in the 1 H NMR spectra signifies the absence of hydroxyl groups resulting from hydrolytic degradation of the triglyceride, which is confirmed by the integrations of the signals of the protons in the alpha position to the oxygen of the ester groups, compared to the integration of the methyl groups.The absence of these peaks also signifies the absence of hydroxyl groups resulting from the epoxide ring opening.However, this does not mean that ring opening by water, leading to vicinal diols, did not take place, but that those elimination reactions giving access to new carbon-carbon double bonds could have taken place.This could explain the presence of a new broad peak between 5.4 and 5.6 ppm.On the other hand, the elimination reaction would also give access to enol functions, which can then be partially or totally converted into ketone functions, which are difficult to characterize from the 1 H NMR spectra of our products.

FTIR Spectra
The FTIR spectra of CSO (cottonseed oil) and ECSO (epoxidized cottonseed oil) are displayed in Figure 3.The conversion of the CH=CH double bond in CSO as the epoxidation reaction proceeds could be verified by the disappearance of the corresponding band at about 3000-3050 cm −1 and the presence of a band at about 850 cm −1 , corresponding to C-O-C stretching vibration of the epoxide ring in CSO.

FTIR Spectra
The FTIR spectra of CSO (cottonseed oil) and ECSO (epoxidized cottonseed oil) are displayed in Figure 3.The conversion of the CH=CH double bond in CSO as the epoxidation reaction proceeds could be verified by the disappearance of the corresponding band at about 3000-3050 cm −1 and the presence of a band at about 850 cm −1 , corresponding to C-O-C stretching vibration of the epoxide ring in CSO.Furthermore, the intermediate samples were also monitored using FTIR (Figure 4), and the signal changes of the epoxidized group could be evidently observed with the extension of the reaction time.Furthermore, the intermediate samples were also monitored using FTIR (Figure 4), and the signal changes of the epoxidized group could be evidently observed with the extension of the reaction time.

GPC
The CSO GPC chromatogram shows a narrow main peak corresponding to the triglyceride molecule (Figure 5).A very small peak of higher molecular weight can be observed, which can be attributed to a siccative minority portion of the triglyceride molecule.The ECSO GPC chromatograms show a gradual shift over time of the main peak to the higher molecular weight region, suggesting that the ring-opening products resulting from the elimination of one or two water molecules are negligible.However, the appearance of a secondary peak that increases in intensity with time suggests the presence of chain extension reactions leading to oligomers.The oligomers started to be significant after 65 min of reaction (less than 15% intensity at the end of the reaction compared to the first peak).In fact, it is well established that epoxides can undergo cationic self-polymerization in an acidic environment, leading to polyethers.

Effect of Reaction Temperature on the Kinetics of Epoxidation
The effect of reaction temperature on the rate of epoxidation was observed by comparing Runs 1, 2, and 3 (Figure 6).One can notice that the rate of epoxidation increases with a reaction temperature increase from 50 °C to 70 °C.Figure 6 also shows that the ring-opening reaction occurs in conjunction with the epoxidation reaction.

Effect of Reaction Temperature on the Kinetics of Epoxidation
The effect of reaction temperature on the rate of epoxidation was observed by comparing Runs 1, 2, and 3 (Figure 6).One can notice that the rate of epoxidation increases with a reaction temperature increase from 50 • C to 70 • C. Figure 6 also shows that the ring-opening reaction occurs in conjunction with the epoxidation reaction.

Effect of HP/PA Ratio on the Kinetics of Epoxidation
The ratio of initial concentrations of hydrogen peroxide to propionic acid ( [HP] 0,as PA 0,as ) in the aqueous phase can play a role in the rate of perhydrolysis, and, thus, in the rate of epoxidation.Figure 7 shows this effect on the rate of epoxidation for Runs 2, 5, and 6. Figure 7 shows that the initial concentration ratio does not significantly impact the initial reaction rate of epoxidation, even if one can notice that an increase in [HP] 0,as PA 0,as has a slight tendency to increase the rate of epoxidation.The kinetic of Run 5 is slightly faster than those of Runs 2 and 6 (lowest ratio [HP] 0,as PA 0,as ).However, after 60 min of reaction, one can notice that the kinetic of Run 5 is faster, and the final concentration of epoxide groups is higher than those of Runs 2 and 6.This phenomenon could be explained by the fact that a higher amount of hydrogen peroxide ensures a faster production of perpropionic acid, and perpropionic acid is more stable than peracetic acid [48], thus, the thermal decomposition of perpropionic acid can be negligible.

Effect of HP/PA Ratio on the Kinetics of Epoxidation
The ratio of initial concentrations of hydrogen peroxide to propionic acid ( , , ) in the aqueous phase can play a role in the rate of perhydrolysis, and, thus, in the rate of epoxidation.Figure 7 shows this effect on the rate of epoxidation for Runs 2, 5, and 6. Figure 7 shows that the initial concentration ratio does not significantly impact the initial reaction rate of epoxidation, even if one can notice that an increase in  and perpropionic acid is more stable than peracetic acid [48], thus, the thermal decomposition of perpropionic acid can be negligible.on the rate of epoxidation.

Catalyst Loading Effect on the Kinetics of Epoxidation
The catalyst Amberlite IR-120 catalyzes the perhydrolysis reaction [48,50], in other words, the rate of perpropionic acid production.Runs 2 and 4 were compared to evaluate the effect of catalyst loading on the rate of epoxidation.Figure 8 shows that the rate of epoxidation increases with the increase in catalyst loading from 5 to 10 g.The presence of ring-opening reactions can be also noticed based on comparing the concentrations of DB and EP at 185 min.

Catalyst Loading Effect on the Kinetics of Epoxidation
The catalyst Amberlite IR-120 catalyzes the perhydrolysis reaction [48,50], in other words, the rate of perpropionic acid production.Runs 2 and 4 were compared to evaluate the effect of catalyst loading on the rate of epoxidation.Figure 8 shows that the rate of epoxidation increases with the increase in catalyst loading from 5 to 10 g.The presence of ring-opening reactions can be also noticed based on comparing the concentrations of DB and EP at 185 min.

Kinetic Modeling
The modeling consisted of two stages: regression to estimate the kinetic constants, and validation to evaluate the reliability of the models.We used a hold-out method, i.e., 10 experiments were used for the regression stage, and two experiments were used for the validation stage [51,52].Tables 1 and 2 show the selected experiments for both stages.
For the epoxidation experiments, we found that by fixing the rotating speed at 600 rpm, there is no effect of the rotation velocities on the kinetics.
Athena visual studio was used for the modeling [53].

Kinetic Modeling
The modeling consisted of two stages: regression to estimate the kinetic constants, and validation to evaluate the reliability of the models.We used a hold-out method, i.e., 10 experiments were used for the regression stage, and two experiments were used for the validation stage [51,52].Tables 1 and 2 show the selected experiments for both stages.
For the epoxidation experiments, we found that by fixing the rotating speed at 600 rpm, there is no effect of the rotation velocities on the kinetics.
Athena visual studio was used for the modeling [53].

Kinetics
Based on the preliminary experimental studies, three main reactions are involved in the epoxidation of cottonseed oil: perhydrolysis of propionic acid, epoxidation of the unsaturated group (-CH=CH-) by perpropionic acid, and oxirane ring-opening side reaction.Figure 9 shows the different reaction steps.The perhydrolysis step, i.e., the production of perpropionic acid, is catalyzed by the cation exchange resins.We assumed that the proton motion from the Amberlite IR-120 is free [50,54].The solubility of hydrogen peroxide in the organic phase is low [34,37,40,41].Thus, the perhydrolysis reaction mainly occurs in the aqueous phase.
The perhydrolysis expression rate was expressed as: where K eq is the equilibrium constant and m dried cat is the Amberlite IR-120 based on a dried basis.The expression rate for the epoxidation was expressed as: where [DB] org is the concentration of the unsaturated group from the cottonseed oil.As from a previous article of our group [40], the ring-opening reaction was derived by assuming that the rate-determining step is the protonation of the oxirane group.Thus, the rate is expressed as: In Equation ( 3), the source of protons comes from the dissociation of propionic acid in the aqueous phase.The proton transfer from the aqueous to the organic phase is supposed to be very fast [37,40].Protons from the Amberlite IR-120 are mainly located inside the resins [50].Thus, one can neglect the protonation of the oxirane group by these protons for steric hindrance reasons.
the resins [50].Thus, one can neglect the protonation of the oxirane group by these protons for steric hindrance reasons.
To decrease the estimation number of parameters, the perhydrolysis kinetics of propionic acid were studied and modeled separately following the same procedure as described in previous articles [55,56].To decrease the estimation number of parameters, the perhydrolysis kinetics of propionic acid were studied and modeled separately following the same procedure as described in previous articles [55,56].

Material Balances
Due to the vigorous stirring, no macroscopic gradient was considered in this study.Furthermore, reaction temperatures were recorded to be constant during a reaction.Thus, isothermal mode was used.
The solubilities of cottonseed oil and cottonseed oil-derivatives in the aqueous phase were neglected.The solubilities of hydrogen peroxide and water in the organic phase were not considered.The aqueous phase was the continuous phase, and the catalyst was assumed to be mainly wetted by this phase during the epoxidation.Volumes of the organic and aqueous phases were assumed to be constant during the reaction.
The double-film theory was used to consider the mass transfer phenomenon, and mass transfer resistance from the aqueous phase, i.e., the continuous phase, was supposed to be negligible.Thus, the concentration of species i at the aqueous interface, noted as [i] * aq , is the same as the concentration of species i in the bulk aqueous phase, noted as [i] aq .
The distribution coefficient or equilibrium molar ratio of a compound i, K i , is defined as: where [i] * org is the concentration of compound i at the organic interface.The concentration [i] * org can be expressed as: . Material balances for the different species are: where k PPA−org and k PA−org are mass transfer coefficients of PA and PPA in the organic phase, a is interfacial area between the organic and aqueous phases, and K PPA and K PA are the distribution coefficients for PA and PPA.In Equations ( 5)-( 13), we used bulk concentrations in the aqueous and organic phases.
In the modeling stage, the coefficients k PPA−org and k PA−org were assumed to be similar due to the similarity between both molecules.We also assumed that these coefficients were temperature-independent.Thus, we estimated the parameter k PA−org •a.
For the estimation of K PPA and K PA , we used the same approach as Wu et al. [35].
In their study, they demonstrated that these coefficients depend on the concentration of double bonds, epoxide groups, and ring-opening products: They also showed that the coefficient distribution of the carboxylic acid and percarboxylic acids are linked by the following relationships: We estimated A, B, and C. We applied Equations ( 14) and ( 15) to express these distribution coefficients.
The volume of the dispersed organic phase, V org , was calculated using the density calculated by Cai et al. [57,58].

Regression
Regression was carried out using the software Athena Visual Studio, based on the Bayesian framework.
Ordinary differential equations from the material balances were solved using a solver named DDAPLUS [53] and were based on a modified Newton algorithm [59].
The non-linear regression was carried out using the GREGPLUS subroutine.This subroutine minimizes the objective function (OF), determines the credible intervals for each estimated parameter represented by the marginal highest posterior density (HPD), and calculates the normalized parameter covariance.
The GREGPLUS subroutine minimizes the objective function via successive quadratic programming [60]: where |υ| is the determinant of the covariance matrix of the responses, b is the number of responses, and a is the number of events in response.Each element of this matrix is: where C iu is the experimental concentration, Ĉiu is the estimated value for response i and event u, C ju is the experimental concentration, and Ĉju is the estimated value for response j and event u.
One should remember that the kinetic modeling for the perhydrolysis and the epoxidation is a multiresponse system, i.e., there are several observables.Van Boekel [58] stated that the Bayesian framework is more adequate than the non-linear least square approach.
It is vital to use a modified Arrhenius equation to express the temperature dependency of the rate constants [61]: where k rx (T) is the rate constant for reaction rx at reaction temperature T, E arx is the activation energy of reaction rx, R is the ideal gas constant, and T re f is the reference temperature.

Kinetic Modeling of Propionic Acid Perhydrolysis
Epoxidation of vegetable oils through the Prileschajew method [17] requires the estimation of several parameters.We have developed a kinetic model only for the perhydrolysis of propionic acid, and for that we only performed experiments in the presence of propi-onic acid, water, hydrogen peroxide, and Amberlite IR-120 (Table 3).Table 3 shows the experimental matrix used for this modeling, and one can notice that the initial operating conditions (concentrations, temperature, catalyst loading) varied for these seven runs.Table 3. Experimental matrix for the modeling of propionic acid perhydrolysis.(HP: hydrogen peroxide, PA: propionic acid, and PPA: perpropionic acid).

Run
Dried Amberlite IR-120 (g) The concentrations of PA and PPA were tracked.Experiments were carried out in isothermal conditions and ideal batch reactors, thus, the material balances led to the following ODEs: A van't Hoff equation was used to express the dependency of K eq towards the reaction temperature.
In this modeling, K eq (T = 303.15)was fixed to 2.051 and ∆H Perh to −4.17 kJ/mol, as in the reference of Leveneur et al. [55,62].
We have estimated the following kinetic constants: ln k Perh T re f and E a Perh R•T re f , displayed in Table 4.The correlation between these two kinetic constants was found to be 0.706, lower than 0.95 [63].Thus, both kinetic constants are not correlated.Figure 10 shows the fit of the model to experimental concentrations for Run 6.One can notice that the developed model fits well with the experimental concentrations.
Figure 10 shows the fit of the model to experimental concentrations for Run 6.One can notice that the developed model fits well with the experimental concentrations.

Kinetic Modeling of Epoxidation For this model, six observables ([Ep]org, [DB]org, [PA]org, [PA]aq, [PPA]aq, and [HP]aq)
were used to estimate the kinetic constants and mass transfer parameters.The perhydrolysis kinetic and thermodynamic constants estimated in the previous section were used.Table 1 shows the experimental matrix used to develop the kinetic model for epoxidation.The ratio minitial,org/minitial,aq was kept constant.
ODEs ( 5)-( 13) were solved, and the following parameters were estimated:    Figure 11 shows the parity plots for PA and PPA, showing the excellent fit of the model.
Figure 10 shows the fit of the model to experimental concentrations for Run 6.One can notice that the developed model fits well with the experimental concentrations.

Kinetic Modeling of Epoxidation For this model, six observables ([Ep]org, [DB]org, [PA]org, [PA]aq, [PPA]aq, and [HP]aq)
were used to estimate the kinetic constants and mass transfer parameters.The perhydrolysis kinetic and thermodynamic constants estimated in the previous section were used.Table 1 shows the experimental matrix used to develop the kinetic model for epoxidation.The ratio minitial,org/minitial,aq was kept constant.
ODEs ( 5)-( 13 ) were used to estimate the kinetic constants and mass transfer parameters.The perhydrolysis kinetic and thermodynamic constants estimated in the previous section were used.Table 1 shows the experimental matrix used to develop the kinetic model for epoxidation.The ratio m initial,org /m initial,aq was kept constant.ODEs ( 5)-( 13) were solved, and the following parameters were estimated: ln k Ep T re f , A, B and C. Table 5 shows the values of the estimated parameters and their credible intervals.Credible intervals are generally low, indicating that the parameters are well identified.The credible interval for k PA−org •a is large, since the temperature effect on this value was not considered.Table 6 shows that all estimated parameters are not correlated, meaning they were well identified.In general, the model fits the experimental data quite well, as shown by the parity plot (Figure 12). Figure 13 displays the fit of the model to experimental data for Run 10. Figure 13 shows that the mass transfer of PA from the aqueous phase to the organic phase is very fast.The model exhibits a slight tendency to underestimate the experimental concentrations of DB and PPA, and slightly overestimates the experimental concentration of epoxide group (Ep).

Validation
A hold-out method was used to validate the model.Table 2 shows the experimental matrix used for the validation.Figure 14 represents the parity plot for the validation stage.One can notice that the model can fit the experimental concentrations for validation.

Validation
A hold-out method was used to validate the model.Table 2 shows the experimental matrix used for the validation.Figure 14 represents the parity plot for the validation stage.One can notice that the model can fit the experimental concentrations for validation.
Figure 15 shows the fit of the model to the validated experimental concentrations of DB and EP for Runs 1V and 2V.For the sake of clarity, Figure 15 shows the concentrations of DB and EP.One can notice that the model slightly underestimates the experimental concentrations of EP. Figure 15 shows the fit of the model to the validated experimental concentrations of DB and EP for Runs 1V and 2V.For the sake of clarity, Figure 15 shows the concentrations of DB and EP.One can notice that the model slightly underestimates the experimental concentrations of EP.

Materials and Chemicals
The following chemicals were used: refined CSO was purchased from ThermoFisher Scientific (Schwerte, Germany).Hydrogen peroxide (33 wt% in water), propionic acid (purity > 99 wt%), tetraethylammonium bromide (TEAB; 99 wt%), and perchloric acid standardized solution in acetic acid (0.1 mol/L) were obtained from VWR International (Fontenay-sous-Bois, France).Chloroform (purity > 99 wt%), iodine solution (0.1 mol/L), ferroin indicator solution, sodium thiosulfate (0.1 mol/L) solution, ammonium cerium (IV) sulfate solution (0.1 mol/L), and amberlite IR-120 were obtained from Sigma-Aldrich (Sigma-Aldrich Chemical Co., USA).Figure 15 shows the fit of the model to the validated experimental concentrations of DB and EP for Runs 1V and 2V.For the sake of clarity, Figure 15 shows the concentrations of DB and EP.One can notice that the model slightly underestimates the experimental concentrations of EP.

Reaction
Experiments were performed in a 300 mL glass water-jacketed reactor equipped with mechanical agitation (pitched blade turbine impeller with a diameter of 3.8 cm and four blades), a reflux condenser, and two temperature probes to measure the temperatures of the reaction mixture and outlet of the jacket.A mixture of cottonseed oil, hydrogen peroxide, distilled water, and amberlite IR-120 was added to the reactor for the epoxidation procedure.When the mixture temperature reached the desired value, a preheated propionic acid solution was added to start the reaction.During the reaction, samples were collected at different times (approximately at 0, 5, 35, 65, 95, 125, 155, 185, 245, and 305 min).After being taken from the reactor, the organic and aqueous phases were separated using a centrifuge (Rotofix 32 A, Hettich, Tuttlingen, Germany) at 6000 rpm/min for 4 min.Then, both phases were analyzed and kept in a fridge at 4 • C. At the end of the reaction, the agitation system was stopped to stratify the mixture into two phases.The same experimental approach was used for the perhydrolysis of propionic acid.
Tables 1-3 display the experimental matrix for the modeling of propionic acid perhydrolysis over amberlite IR-120 for the regression stage, the modeling of cottonseed oil epoxidation by perpropionic acid produced in situ, and the validation stage for the epoxidation model, respectively.

Analytical Methods
After collecting the samples in the reaction stage and separating by centrifuging, the two phases were then analyzed using titration and spectroscopy methods to determine the concentration of double bonds, epoxidized groups, and carboxylic acid in the organic phase.The concentration of acids and hydrogen peroxide in the aqueous phase was measured using titration.NMR, FTIR, and GPC were used to characterize the organic phase further and get inside for the polymerization and ring-opening reactions.

NMR
1 H NMR analysis was performed on a Bruker 300 Fourier Transform spectrometer (Bruker Corporation, Billerica, MA, USA) at 300 MHz in CDCl 3 solutions containing tetramethylsilane (TMS) as an internal standard.

FTIR
FTIR spectra were acquired with a Perkin-Elmer Spectrum 2000 FTIR (Perkin-Elmer, Waltham, MA, USA) equipped with a diamond ATR device (attenuated total reflection).Spectra were recorded from 10 scans in the 650 to 4000 cm −1 range.

GPC
The average molecular weights of polymers (Mn and Mw) and dispersity (Ð = Mw/Mn) were assessed through size exclusion chromatography (sEc).Polymers were dissolved in dichloromethane, filtered (0.45 µm), and analyzed at 25 • C using a Varian PL-GPC50 system (Agilent Technologies, Santa Clara, CA, USA) equipped with two mixed packed columns (PL gel mixed type C).Dichloromethane was used as the mobile phase and PMMA standards (from 875 to 62,000 g mol −1 ) were used for calibration.

Double Bond Content
The concentration of the double bonds was measured using the iodine value, which was determined by the Wijs method [64].Briefly, 0.20 g of oil sample or organic phase was dissolved in 20 mL of chloroform.Then, 25 mL of Wijs solution (0.1 mol/L iodine monochloride in glacial acetic acid) was added and kept in a dark place for 1 h.After a onehour reaction, 20 mL of 10% potassium iodide solution was added to stop the reaction.After adding 100 mL water, the mixture solution was titrated using sodium thiosulphate solution (0.1 mol/L) with an automatic titrator (916 Ti-Touch, Methrohm, Herisau, Switzerland).

Epoxide Content
The concentration of epoxide groups was determined by the method of Jay [65].Briefly, 0.10 g of oil sample or organic phase was dissolved in 10 mL of chloroform.After adding 10 mL of 20% TEAB solution in acetic acid, the mixture solution was titrated using standard perchloric acid in acetic acid solution (0.1 mol/L) with an automatic titrator (916 Ti-Touch, Methrohm Herisau, Switzerland)).

Hydrogen Peroxide Content
The concentration of hydrogen peroxide was determined by the method of Greenspan and Mackellar [66].Briefly, 0.20 g of aqueous phase was dissolved in 50 mL of 10% sulfuric acid in water solution.Then, three drops of Ferroin indicator solution were added and titrated using standard ammonium cerium sulfate solution (0.1 mol/L) until the color of the mixture solution reached a light blue endpoint.

Concentration of Acid
The concentration of acid was determined by the titration method [30].For the concentration of acid in the aqueous phase, a 0.30 g sample was dissolved in 50 mL distilled water and titrated using sodium hydroxide solution (0.1 mol/L) with an automatic titrator (916 Ti-Touch).For the organic phase, a 0.50 g sample was dissolved in 50 mL 2-propanol and titrated in the same way as above.

Conclusions
The production of epoxidized vegetable oils by the Prileschajew method is still the most used because it circumvents the low hydrogen peroxide solubility by producing a percarboxylic acid in situ.Nevertheless, this method has two drawbacks: thermal instability of the produced percarboxylic acid and the oxirane ring-opening side reaction by protons.A trade-off between the percarboxylic acid thermal instability and its reactivity was found by using perpropionic acid.The use of a solid acid catalyst accelerates the production of percarboxylic acid, and lowers the ring-opening side reaction because the accessibility of the acid sites to the oxirane groups due to steric hindrance is reduced.For this reason, amberlite IR-120 was used in this study.The kinetics of cottonseed oil epoxidation by perpropionic acid produced in situ was studied in an isothermal batch reactor.
A biphasic kinetic model taking into account the mass transfer of propionic acid and perpropionic acid from the aqueous continuous phase to the organic dispersed phase was developed and assessed.A Bayesian approach was used to estimate and evaluate kinetic constants and mass transfer coefficients.We developed the following kinetic modeling strategy to decrease the number of parameters to estimate: first, only model the kinetics of propionic acid perhydrolysis, which produces perpropionic acid; second, model the kinetics of the epoxidation of cottonseed oil by perpropionic acid produced in situ.A hold-out method was used to validate the epoxidation kinetic model.An in-depth analysis (NMR, GPC, FTIR, titration) was carried out to understand the ring-opening reaction, and we found that water and diol were responsible for this side reaction.
The cottonseed oil epoxidation model was found to be reliable for the following initial operating conditions: dried catalyst amount from 2 to 15 g, reaction temperature from 323.15 to 353.15 K, double bond initial concentrations from 3.33 to 3.85 mol/L, initial propionic acid concentration in the aqueous phase from 3 to 7 mol/L, and initial hydrogen peroxide concentrations from 3 to 7.01 mol/L.
A continuation of this work could be the comparison with other carboxylic acids, and an investigation on the catalyst interaction between the aqueous and organic phases.
the rate of epoxidation.The kinetic of Run 5 is slightly faster than those of Runs 2 and 6 (lowest ratio , ,).However, after 60 min of reaction, one can notice that the kinetic of Run 5 is faster, and the final concentration of epoxide groups is

Figure 6 .
Figure 6.Effect of reaction temperature on the rate of epoxidation.

Figure 7 .
Figure 7. Effect of[HP] 0,as PA 0,as on the rate of epoxidation.

Figure 8 .
Figure 8.Effect of catalyst loading on the rate of epoxidation.

4 Figure 8 .
Figure 8.Effect of catalyst loading on the rate of epoxidation.

Figure 10 .
Figure 10.The fit of the model to experimental concentrations for Run 6. (PPA: perpropionic acid and PA: propionic acid).

Figure 11
Figure11shows the parity plots for PA and PPA, showing the excellent fit of the model.

Figure 10 .
Figure 10.The fit of the model to experimental concentrations for Run 6. (PPA: perpropionic acid and PA: propionic acid).

Figure 11
Figure11shows the parity plots for PA and PPA, showing the excellent fit of the model.

Figure 15 .
Figure 15.The fit of the model to Runs 1V and 2V for the validation stage.(DB: double bond and EP: epoxide).

Figure 15 .
Figure 15.The fit of the model to Runs 1V and 2V for the validation stage.(DB: double bond and EP: epoxide).

Table 1 .
Experimental matrix for the modeling of epoxidation of cottonseed oil.(DB: double bond, EP: epoxide, PA: propionic acid, and HP: hydrogen peroxide).

Table 2 .
Experimental matrix for the validation of epoxidation of cottonseed oil model.(DB: double bond, EP: epoxide, PA: propionic acid, and HP: hydrogen peroxide).

Table 4 .
Kinetic constants and credible intervals at T ref = 343.15K.

Table 5 .
Estimated parameters and their credible intervals for the epoxidation of cottonseed oil at T ref = 343.15K.

Table 6 .
Correlation matrix of estimated parameters.