May 28, 2023

Separation of hydrocarbons from activated carbon as a porous substance in a glycol regeneration process using supercritical carbon dioxide

Scientific Reports volume 12, Article number: 19910 (2022) Cite this article

1190 Accesses

1 Citations

1 Altmetric

Metrics details

A Publisher Correction to this article was published on 14 December 2022

This article has been updated

Activated carbons are used in industrial applications; their cost is a major barrier to their more widespread application. Regeneration of activated carbons is indispensable to minimize operational costs and product wastage. Supercritical carbon dioxide (SC-CO2) as green technology was used to regenerate activated carbons. In this work, response surface methodology was employed to optimize the supercritical regeneration process and to evaluate the effect of the operational parameters including pressure (100–300 bar), temperature (313–333 K), flow rate (2–6 g/min), and dynamic time (30–150 min) on the regeneration yield. The maximum regeneration yield (93.71%) was achieved at 285 bar, 333 K, 4 g/min, and 147 min. Mathematical modeling was done using two one-parameter kinetics models, which agree well with the experimental data. The fitting parameter of the model was obtained by using a differential evolution algorithm. The chemical composition of the substances extracted from the activated carbon was identified by gas chromatography. The results showed that the regeneration of activated carbon by SC-CO2 can be an alternative method to conventional methods.

Activated carbon has been recognized as one of the widely applied adsorbents for spray coating1, food processing2, biomass3, pharmaceuticals4, chemicals5, waste water treatment6 petroleum7, and nuclear industries8, as well as the glycol treatment for removing organic and volatile organic compounds (VOCs) pollutants in natural gas industries9. Activated carbon is used in natural gas sweetening and dehydration systems with fluids such as amines and glycols10,11. The conventional techniques for the regeneration of activated carbon include thermal volatilization12, chemical extraction13, ultrasound14, microwave15, electrochemical16, and bio-regeneration. These methods have several disadvantages such as loss of carbon, damage of its porous structure , treatment of exhaust gases, chemical regenerations using solvents are not necessarily acceptable because additional separation and environmental problems and bio-regeneration requires long reaction time for regeneration. Recently, supercritical fluids (SCFs) has attracted widespread attention in many fields and regeneration of activated carbon as one of the applications of this technology has been studied17. The unique characteristics of SCFs have made these solvents attractive. Particularly, solvent density, and hence its dissolvent properties, can be controlled modifying pressure and temperature. Besides, liquid-like density and gas-like viscosity, coupled with diffusion coefficients that are at least an order of magnitude higher than those of liquids, contribute to the enhancement of the mass transfer processes18,19,20,21,22,23. Among the different substances, carbon dioxide is the best choice as it is an environment-friendly solvent providing such advantages as nontoxicity and high chemical stability.

Regeneration of activated carbon using SCFs has been studied by several researchers. DeFilippi et al. observed that the supercritical regeneration was economical even though the operating temperature and pressure were above 387 K and 150 atm, respectively. They proposed a local equilibrium model (Freundlich isotherm) fitting well with the experimental data24. Regeneration of activated carbon loaded with phenol using supercritical carbon dioxide was studied by Kander and Paulaitis25. They found that supercritical carbon dioxide offered no significant advantages for the regeneration of carbon loaded with phenol. However, they suggested that for organic compounds which are not adsorbed strongly onto activated carbon, supercritical carbon dioxide would be a powerful adsorbent. Tan and Liou investigated desorption by supercritical carbon dioxide of activated carbon loaded with either ethyl acetate or toluene. They presented that this regeneration method would give better results than the steam regeneration method and accordingly, presented a linear desorption kinetics model that was found to fit experimental data quite well17,26.

Madras et al.27 studied the desorption of several nonvolatile solids such as hexachlorobenzene and pentachlorophenol via break through experiments at a fixed carbon dioxide density. Macnaughton and Foster28 measured the adsorption equilibrium of 1,1,1-trichloro-2,2-bis (p-chlorophenyl) ethane (DDT) on an activated carbon and illustrated that one limiting factor on supercritical desorption process is adsorption equilibrium at a low CO2 flow rate. Using a three-parameter mode, the adsorption breakthrough and equilibrium of ethylbenzene on an activated carbon from supercritical carbon dioxide was investigated by Harikrishnan et al.29. Tan and Liou26 studied desorption behavior of ethyl acetate from activated carbon in CO2 at temperatures from 300 to 338 K and pressures from 8.7 to 12.9 MPa. They observed that increasing pressure and decreasing temperature increased the desorption rate. Benkhedda et al.30 also presented similar temperature and pressure dependence for m-xylene desorption behavior using SC-CO2 at temperatures from 313 to 333 K and pressures from 10.0 to 15.0 MPa. They suggested that desorption behavior of VOCs from activated carbons using SC-CO2 largely depends on temperature and pressure conditions. Salvador et al.31 studied the regeneration of activated carbons contaminated with phenol using supercritical water at 260 bar and temperature ranging between 400 and 500 °C. Regeneration temperature and time were the variables that most affected the process. Ikuo Ushiki et al.9 investigated the desorption behavior of various VOCs (toluene, acetone, n-hexane, n-octane, methanol, ethanol, 2-propanol, and propylene glycol monomethyl ether) from activated carbon using a fixed-bed method employing SC-CO2 at temperature values ranging between 313 and 353 K and pressure values ranging between 10.0 and 15.0 MPa.

Several factors such as pressure, temperature, CO2 flow rate, time, particle size, among others, simultaneous influenced the supercritical regeneration of activated carbons. In most of the previously published studies, the process conditions have just been optimized by conducting one-factor- at-a-time experiments. The results of one-factor-at-a-time experiments do not reflect real changes in the environment as they ignore interactions between the different variables of the process. Thus, to avoid this problem, the design, optimization and assessment of the process using design of experiment (DOE) seems to be essential32. Response surface methodology (RSM) is a mathematical tool to specify the effect of each variable along with their interactions, on the yield of process, and to predict its behavior under a given sets of conditions. The main advantage of the RSM is the reduced number of experiments required to optimize a process. RSM including Central Composite Design (CCD) and Box–Behnken Design (BBD) have been applied to optimize a supercritical fluid extraction (SFE) process. CCD is usually more expensive than BBD due to requires a higher number of runs than BBD33.

However, despite the fact that there can be various contaminants that need to be removed to regenerate activated carbons, most of these works investigated polluting substances intentionally added to the activated carbon. To our knowledge this is the first study that propose the use of SC-CO2 to regenerate activated carbon used in the gas industry, developing a comprehensive study to optimize and modeling the process. CCD design as a powerful experimental tool was applied to evaluate the influence of operating parameters on the regeneration process. In addition, the process was modeled via mathematical modeling by correlating the experimental data. Moreover, the extract obtained from the SC-CO2 method was characterized by GC and the physical properties of activated carbon were studied by SEM and iodine number methods, respectively.

The activated carbon samples used in this study were taken, at summer 2017 during three months, from a gas dehydration unit of the South Zagros Oil and Gas Company (Iran). It should be remarked that, in gas the gas dehydration units activated carbon is used to remove hydrocarbons and aromatic compounds in triethylene glycol. CO2 (99.90% purity) was purchased from Aboughadareh Gas Factory.

Raw natural gas is fully saturated with water vapor when produced from an underground reservoir. Because most of the water vapor has to be removed from natural gas before it can be commercially marketed, natural gas is subjected to a dehydration process. One of the most substances used for removing the water from produced gas is glycol. Figure 1 shows a scheme of the typical equipment used for the dehydration process using glycol. While the overall process equipment is similar for all glycol dehydration units, there can be considerable variation among installations. The gas flows through a separator to remove condensed liquids or any solids that might be in the gas. Some absorbers incorporate the separator at its bottom section, in which case the gas then flows upward through a chimney tray into the glycol absorber portion of the vessel. The glycol contactor provides the close contact between the gas and the glycol. The glycol is highly hygroscopic, and most of the water vapor in the gas is absorbed by the glycol. The water-rich glycol is withdrawn from the contactor near the bottom of the vessel above the chimney tray through a liquid level control valve and passes to the regeneration section. The treated gas leaves the contactor at the top through a mist eliminator and usually meets the specified water content. The water-rich glycol can be routed through a heat exchange coil in the top of the reboiler column called the steel column. The heat exchange generates some reflux for the separation of the water from the glycol in the top of the steel and heats the rich glycol somewhat. In some installations, the rich solution passes to a flash tank operating pressure at about 15 to 50 psig, which allows to the absorbed hydrocarbon gas to separate from the glycol. The glycol then flows into the still through an activated carbon filter and a heat exchanger, exchanging heat with the regenerated glycol. Then, it drops through a packed section in the still into the glycol reboiler vessel where it is heated to the necessary high regeneration temperature to near atmospheric pressure. At the high temperature, the glycol loses its ability to hold water; the water is vaporized and leaves through the top of the still. The regenerated glycol flows to the surge tank, from which it is routed through the lean/rich heat exchanger to the glycol pump. The pump boosts the pressure of the lean glycol to the contactor pressure. Prior to entering the contactor, it exchanges heat with the dry gas leaving the contactor or some other heat exchange medium34.

Schematic of unit process.

The main components of the apparatus were a CO2 cylinder (E-1), needle valve (E-2), molecular sieve filter (E-3), refrigerator unit (E-4), air compressor (E-5), high-pressure pump (air-driven liquid pump, Haskel, USA) (E-6), oven (Memert) (E-7), Surge tank (E-8), stainless-steel impregnation cell (E-9), back-pressure valve (Xi'an Shelok Instrument Technology Co., Shaanxi, China) (E-10), Flowmeter (E-11), Sampler (E-12), and automation system (E-13). At first, CO2 from the reservoir passed through a filter that is rated as a 1µ nominal size to avoid passing any impurity to the process path. Then, CO2 was guided towards a refrigerator unit, in which its temperature was drop to -10 °C and 60 Bar. The liquefied CO2 was pressurized using a reciprocating pump, and translated towards shell and tube configuration surge Drum. At this drum the fluid temperature was set to proper temperature to reach the supercritical fluid condition. This drum has been designed to reduce the fluctuations and the heat exchange happens by circulating water using 1000 W elements with a temperature accuracy setting of 1 K. A pressure gauge was installed at the output of surge drum to illustrate the outlet fluid pressure. Afterward the fluid was guided towards an extraction cell loaded with 10 g of activated carbon along with glass beads. At the input and output of the extraction cell two filters having Porosity of 1 micron have been installed to prevent the escaping of particles. Existence of glass beads in cell allowed the homogenous distribution of CO2 inside the extraction cell. After loading the extraction cell with SC-CO2, the cell was set at appropriate temperature and pressure for about 45 min. Following the static stage, the dynamic one begins by opening the sample valve. To prevent the valve from freezing, it was wrapped in a heating element. During the dynamic time, sample collection temperature was kept below 0 °C by using ice booth. Each experiment was done in triplicate35. Figure 2 shows the details of supercritical fluid apparatus. The yield of supercritical recovery was calculated through dividing the solute weight per the initial activated carbon. Weight of initial solute in Eq. (1) was determined as follows, first, a known mass of activated carbon was loaded into the filter. Then, after three months, the active carbon is going out from the filter. The difference in the primary and secondary weight of activated carbon indicates the initial value in Eq. (1).

Schematic of Supercritical carbon dioxide extraction.

Response surface methodology (RSM) was employed to evaluate the effects of pressure (A), temperature (B), CO2 flow rate (C) and time (D) on the regeneration yield. The coded and uncoded independent parameters applied in the RSM and their respective levels were tabulated in Table 1. The remaining independent parameters (e.g., particle size and volume of cell) were kept constant during the experimental procedures. The experimental design was based on the central composite design (CCD) using a 30 factorial and star design with six central points as shown in Table 2. Although the cost of this method is more than the Box Behnken method, the authors used this method to evaluate the effect of their work (According to the request of oil industry).

All experimental runs, which included 8 factorial points, 6 axial points and 5 center points, were based on run order performed. Furthermore, second-order polynomial equation was used to express the yield as a function of the independent parameters.

Design expert software (Trial version 7.1) was applied for statistical treatment of the results. Analysis of variance (ANOVA) was employed to determine the statistically significant parameters and interactions using Fisher’s test and its associated probability p(F). The determination coefficients, R2, and their adjusted values, R2,adj, were used to evaluate the goodness of fit of the regression models.

A Varian CP-3800 gas chromatograph (Varian instruments) with an FID detector and CP-sil 9 CB capillary column (100 m 0.25 mm, 0.25 μm film thickness) was used to accomplish chemical analysis in this work. The oven temperature was maintained at 35 °C for 7 min before being raised up to 250 °C at 3 °C/min. Injector and detector temperatures were set at 275 °C and 275 °C, respectively. Helium was used as carrier gas at a flow rate of 1.6 ml/min with the samples being injected manually under a split ratio of 1:80. Peaks’ area percentages were used to obtain quantitative data by using of the DHA software.

The iodine number (IN) was determined according to the ASTM D4607-94 method36,37. The iodine number is defined as the milligrams of iodine adsorbed by 1.0 g of carbon when the iodine concentration of the filtrate is 0.02 N (0.02 mol/L). This method is based upon a three-point isotherm. A standard iodine solution is treated with three different weights of activated carbon under specified conditions. The experiment consists of treating the activated carbon sample with 10.0 mL of 5% HCl. This mixture is boiled for 30 s and then cooled. Soon afterwards, 100.0 mL of 0.1 N (0.1 mol L-1) iodine solution is added to the mixture and stirred for 30 s. The resulting solution is filtered and 50.0 mL of the filtrate is titrated with 0.1 N (0.1 mol L-1) sodium thiosulfate, using starch as indicator. The iodine amount adsorbed per gram of carbon (X/M) is plotted against the iodine concentration in the filtrate (C), using logarithmic axes.

The iodine number is the X/M value when the residual concentration (C) is 0.02 N (0.02 mol.L-1). The X/M and C values are calculated by the Eqs. 3 and 4 respectively.

where N1 is the iodine solution normality, V1 is the added volume of iodine solution, VHCl is the added volume of 5% HCl, VF is the filtrate volume used in titration, NNa2S2O3 is the sodium thiosulfate solution normality, VNa2S2O3 is the consumed volume of sodium thiosulfate solution and MC is the mass of activated carbon.

Differential evolution algorithm (DEA) is a population-based algorithm like genetic algorithms by similar operators; crossover, mutation, and selection. The main difference in constructing better solutions is that genetic algorithms depend on crossover while DE relies on mutation operation. This main operation is found on the differences of randomly sampled pairs of solutions in the population. The algorithm uses mutation operation as a seek mechanism and selection operation to direct the search toward the probable regions in the search space. The DE algorithm also uses a non-uniform crossover that can take child vector parameters from one parent more often than it does from others. Among the DE’s advantages its simple structure, ease of use, speed, and robustness can be mentioned38,39,40.

CCD was used for the design of experiments and process optimization. In this way, four coded operating parameters were defined in five levels, namely − 2, − 1, 0, 1 and 2, to obtain the optimal regeneration yield. Pressure (100–300) bar, temperature (313–333 K), flow rate (2–6 g/min), and dynamic time (30–150 min) were considered as the operating parameters of the process to be studied. Table 2 reports the experiments responses at the different values of operating conditions maintaining constant the particle size at 1 mm. Among models presented in CCD, the quadratic model represents the best model for the prediction and optimization of the regeneration yield. Equation (5) represents the obtained quadratic polynomial model through which the yield was correlated, as a response to the coded independent parameters. Therefore, the RSM predictive model was obtained as follows:

This model was found to be significant at 95% level of confidence while the corresponding F-value and p-value were calculated to be 92.58 and 0.0001, respectively. Also, accuracy of the model can be tested by investigating the determination coefficient (R2). The values of the coefficient of determination (R2) and adjusted coefficient of determination (R2adj.) were calculated to be 98.86 and 97.79%, respectively. The value of R2 indicates a good agreement between the experimental and predicted response values. The value of adjusted R2 revealed that only 2.21% of total variations failed to be explained by the model. The lack-of-fit measures the failure of the model to represent data in the experimental domain at points which are not included in the regression. The p-value of the lack-of-fit was higher than 0.05 indicating an excellent fit. Furthermore, Signal-to-noise ratio (SNR) is known as an adequate precision measure. In fact, SNR compares the range of predicted values at design points to the average prediction error. As a SNR greater than 4 represents desirable results, the achieved ratio (35.30) served as an adequate signal as it was much smaller than the actual effect size. Low C.V. (3.32%) indicates the reliability of the carried-out experiments. The ANOVA of quadratic model was carried out using the data tabulated in Table 3 with the purpose of examining significance of the variables as linear, quadratic and interaction coefficients of the RSM. Those variables and their interactions with higher regression coefficient and smaller p-value (p < 0.05) have a significant influence on the regeneration yield41,42. Analysis results reported in Table 3, as well as, the parameters which statistically indicated highly significant impact on the yield were linear term of pressure, CO2 flow rate and time (p < 0.0001), followed using the quadratic term of flow rate(p < 0.0001). The quadratic term of time (p = 0.0003), the linear term of temperature (p = 0.0037) as well as the quadratic term of pressure and temperature (p = 0.0005), represent significant influence on the yield.

In this section the effect of the different operational parameters on the supercritical regeneration of activated carbon were examined with three dimensional graphs, while the other variables were maintained at its respective fixed middle level(Pressure, 200 bar, Temperature, 318 K, flow rate,4 g/min and time 90 min), corresponding to zero code.

Figure 3 shows the 3D plot for the influence of pressure and temperature. It can be observed that the regeneration yield increased with increasing pressure and temperature from 100 to 300 bar and from 313 to 333 K, respectively. Based on CCD, pressure had the most important influence on the regeneration yield. Particularly, the increase of pressure leads to increase solvent power of CO2 and therefore solubility of the solutes increased as well9,43. On another hand, the increase in temperature had a slight positive effect on the regeneration yield (Fig. 4). At a pressure of 100 to 300 bar, with the increase in temperature, the solubility of solute increases, and it is due to the increment of the solute vapor pressure effect. therefore, the solubility of solute increases with the increase in temperature. On the other hand, the extract vapor pressure is raised with increasing temperature, which leads to an enhancement in the supercritical fluid diffusivity. Figure 5 denotes the 3D plot of the yield as a function of pressure and flow rate at 323 K and 90 min. Increasing the CO2 flow rate from 1 (code =  − 2) to 4 g/min (code = 1) reduces the film thickness around the particles leading to lower resistance for mass transfer around the particles and accordingly enhancing the regeneration of the activated carbon. Whereas at more than 4 g/min with reduction of residence time (SC-CO2- particles contact time) a negative effect of flow rate was observed on the regeneration yield. This contrast effect was also observed in Figs. 6, 8 and 9.

3D graphic surface for the effects of temperature and pressure on the yield.

3D graphic surface for the effects of temperature and time on the yield.

3D graphic surface for the effects of flow and pressure on the yield.

3D graphic surface for the effects of temperature and flow on the yield.

Figures 4, 7, and 8 illustrate the effect of dynamic time and flow rate on the regeneration yield. As can be seen in Figs. 4, 7, and 8 the regeneration yields gradually increase with increasing the dynamic time, achieving its highest value about 150 min. This behavior has been reported by other researchers and explained in terms of the increase of the ratio between SC-CO2 and the solute as the dynamic time increase35. On another hand, the residence time decrease as flow rate increase, while the external mass transfer coefficient increased, so these opposite phenomena canceled their effects out leading to the almost constant yield. Moreover, the increase of the flow rate could reduce the contact time between the solvent and the activated carbon. This phenomenon could be related to the channeling effect, where SC-CO2 at high flow rates would just flow around the samples with no ability to diffuse through the pores within the samples. Furthermore, the increased flow rate could cause the sample compaction limiting the amount of CO2 able to being in direct contact with the sample matrix. Similar results were reported by some researchers33,44,45,46,47,48. In addition, the perturbation plot (Fig. 9) revealed the significant effect of all process variables on the extraction. A perturbation plot does not show the effect of interactions and it is like one factor-at-a-time experimentation. The perturbation plot helps to compare the effect of all independent variables at a particular point in the design space. The response is plotted by changing only one factor over its range while holding the other factors constant.

3D graphic surface for the effects of pressure and time on the yield.

3D graphic surface for the effects of flow and time on the yield.

Perturbation plot of operational parameters obtained through CCD on the yield.

The optimal conditions to obtain the highest yield of regeneration from activated carbon were determined at 285 bar, 333 K, 4 g/min, and 147 min and the predicted yield was 93.75%. Average actual yield (94.25%) was in good agreement with the estimated value, revealing the ability of the developed model in terms of extraction process prediction and optimization.

The solute composition determined for the extract, at the optimum condition, by gaseous chromatography (GC) were shown in Table 4. Thirty-two components were identified in the extract obtained using SC-CO2 which comprised 95.2% of the extract. The main components extracted by SC-CO2 were N-C6 (4.90), N-C7 (9.45%), N-C8 (9.04%), N-C9 (11.90%), N-C10 (12.40%) and N-C11 (4.73%).

Iodine number is a measure of iodine molecules adsorbed in the pores of a particle, which indicates the pore volume capacity and extent of micro pore distribution in the activated carbon. Iodine number can be correlated with the adsorbent ability to adsorb low molecular weight substances. Iodine adsorbed results in extracted and non-extracted activated carbon samples were shown in Table 5. Results of physical and adsorptive characteristics denote that SC-CO2 appears to be a suitable way for the regeneration of activated carbon.

As was previously mentioned the structure and surface morphology of non-extracted and extracted activated carbon samples was evaluated by scanning electron microscopy (SEM). Figure 10a,b presents two micrographs of activated carbon before (a) and after supercritical regeneration (b). The micrograph of the unprocessed sample shows an uneven and rough surface covered uniformly with a layer of solute. The structure of the treated sample is more porous and the surface is clearly deflated and depleted of solute due to the high recovery efficiency obtained by the extraction with supercritical carbon dioxide.

(a,b) SEM image of (a) un-process, (b) process of activated carbon.

The experimental desorption was correlated with a kinetic model proposed by Tan and Liou17,49. This model assumes linear desorption kinetics in the adsorbed phase, and has a high applicability for correlating desorption behavior using SC-CO2 with only one fitting parameter. The following assumptions were considered in the formulation of the model:

The system was isothermal and isobaric.

The physical properties of the SC-CO2 were constant during the extraction process.

The extraction model was expressed as an irreversible desorption process.

All particles were considered to be spherical and the solutes were uniformly distributed in their structures.

The volume of the solid matrix (particles) was not changed during the extraction process.

The solvent flow rate was constant along the bed and uniformly distributed without radial dispersion.

Axial dispersion was neglected

A schematic diagram of the particles and bed was presented in Fig. 11. According to the above assumptions, the material balances for the solute in the solid and bulk phases are as follows:

Bulk phase:

A schematic diagram of the particles and bed.

Solid phase:

The concentration at the exit of the packed bed can then be obtained by:

K is the adjustable parameter of the kinetic model.

The adjustable parameter was computed by minimizing the errors between experimental and calculated yield values. The average absolute relative deviation (AARD), described using the following equation was applied to evaluate the adjustable parameter:

The mathematical modeling was applied to study the impact of pressure, temperature and flow rate as on the regeneration process. Figure 12 and Table 6 show the modeling results by using the kinetics model. Fig 12a shows a positive effect of the pressure on the regeneration process at fixed temperature, flow rate, particle size and dynamic time. This effect was related to the increase of the SC-CO2 density as pressure increase. As was expected, the increase of density leads to an increase in the solubility of the solutes in SC-CO2. On the other side, increasing pressure decreased the diffusion coefficient of CO2. The decreasing effect of mass transfer coefficient may lead decreasing the yield. Nevertheless, the effect of increasing density and solubility overcomes the decreasing effect of diffusivity on the final yield of extraction22,50. The adjustable parameter of the model (K) and AARD values for the different extraction runs were reported in Table 6. The effect of temperature on the regeneration yield of activated carbon was shown in Fig. 12b. The yield increased with increasing the temperature at fixed pressure and flow rate. As previously mentioned, this may be due to the increase of the vapor pressure of the solute. In the present section, the effect of temperature on the regeneration yield at points of 318, 328 K and optimum temperature were modeled. The flow rate of carbon dioxide was another parameter that was investigated. The experimental data and results of model in three flow rates including 1, 2 and 4 (optimum) were shown in Fig. 12. Based on the results indicated in Fig. 12c, the yield was higher at optimal points. As can be seen in Fig. 12 the experimental data were well described with the kinetics model. Values of average absolute relative deviation are in the range 5.91-9.04%.

Effect of process parameters on the yield of extraction versus time.

Single sphere model (SSM) proposed by Crank51 with an assumption that particle size is one of the main factors on the diffusivity step in the extraction process. In addition, diffusion is assumed to occur in the sphere surface area of a particle solute. Dissemination of SC-CO2 in sphere equations is conducted to determine the diffusiveness of solvent to dissolve in matrices. Other assumptions made for an SSM model are: (1) The resistance of mass transfer is zero, (2) all of the particle sizes are homogenous, (3) the main factor in the extraction process is intraparticle of mass transfer, (4) the solute is in the inert porous sphere, and lastly, (5) all of the solutes in the bed will be extracted, and the extracted component will be dissolved in the particles by a process similar to diffusion. Equation 13 exhibits the diffusion equation for a constant coefficient.

where Mt is the total amount of the diffusing substances entered on the sheet at a specific time, M∞ is the particular quantity after countless time, Deff is the diffusivity, R is the radius of particle solute, and t is the time.

The single sphere model is usually applied to determine the diffusivity coefficient and mass transfer between solvent and solute. Compared to another kinetic model, the single sphere model is easily applied, particularly the shrinking core model52 and the broken intact cell53 due to one adjustable parameter. This is because one adjustable parameter is fit enough to determine the mass transfer process of experimental data. Furthermore, a single-sphere model can assess the effect of the parameter on the diffusivity process between solvent and solute54.

The single sphere model was fitted to the experimental data with effective diffusivity as a fitting parameter, using average absolute deviation (AARD). Table 7 shows the AARD and effective diffusivity. The minimum %AARD was 2.08% at 285 bar and 333 K, while the maximum %AARD was 19.04% at 200 bar and 323 K. As shown in Fig. 13 the model successfully fitted the exponential trends of the extraction process. Putra et al.1 applied single sphere model for fitting the experimental data of modified supercritical carbon dioxide with error below 5% with the highest diffusivity coefficient was 6.794 × 10‐12 m2/s at operating condition was 10 MPa, 40 °C and particle size 425 µm. Moreover, Aris et al.2 determined Momordica charantia extract yield with different mean particle size as well as diffusion coefficient, De, in the extraction process with and without co-extractant. Based on the results, mean particle size of 0.3 mm gave the highest extract yield, 3.32% and 1.34% with and without co-extractant respectively. Whereas, the value of De at 0.3 mm mean particle size, with and without co-extractant are 8.820 × 10−12 and 7.920 × 10−12 m2/s respectively.

Effect of process parameters on the yield of extraction versus time.

Response surface methodology (RSM) was used in order to find out the effect of pressure, temperature, flow rate and dynamic time on the regeneration yield of activated carbon using SC-CO2. The optimal conditions for regeneration were found at 285 bar, 333 K, 5 g/min and 147 min which resulted in a regeneration yield equal to 93.71%. The regeneration effectiveness was revealed by the iodine number and visualized by SEM micrographs. The kinetic models proposed by Crank and Tan et al. successfully described the regeneration process using SC-CO2 and the estimated regeneration yields using the model adequately agreed with the experimental data for all studied conditions. The model parameters (K) and effective diffusivity was obtained by best fitting procedure applying differential evolution (DE) algorithm. Gas chromatography analysis was performed to identify the composition of solutes at the optimal extraction condition.

All data are available within the published paper.

A Correction to this paper has been published:

Refers to calculated

Concentration of the solute in Co2 (mol/m3)

Exit concentration of solute (mol/m3)

Effective diffusivity of solute (m2/s)

Particle size (mm)

Refers to experimental

Adjustable parameter

Length of extractor (m)

Number of data

Pressure (bar)

Mass flow rate of CO2 (kg/s)

Radius of particle

Time (s)

Temperature (K)

Solvent velocity (m/s)


Axial position in column

Density (kg/m3)

Bed porosity

Chen, L. et al. Optimization of a Fe–Al–Ce nano-adsorbent granulation process that used spray coating in a fluidized bed for fluoride removal from drinking water. Powder Technol. 206, 291–296 (2011).

Article CAS Google Scholar

Lu, J. et al. Removal of imidacloprid and acetamiprid in tea (Camellia sinensis) infusion by activated carbon and determination by HPLC. Food Control 131, 108395 (2022).

Article CAS Google Scholar

Ouyang, J., Zhou, L., Liu, Z., Heng, J. Y. & Chen, W. Biomass-derived activated carbons for the removal of pharmaceutical mircopollutants from wastewater: A review. Sep. Purif. Technol. 253, 117536 (2020).

Article CAS Google Scholar

Vedenyapina, M., Kurmysheva, A. Y., Rakishev, A. & Kryazhev, Y. G. Activated carbon as sorbents for treatment of pharmaceutical wastewater, Solid Fuel. Chemistry 53, 382–394 (2019).

CAS Google Scholar

Jeirani, Z., Niu, C. H. & Soltan, J. Adsorption of emerging pollutants on activated carbon. Rev. Chem. Eng. 33, 491–522 (2017).

Article CAS Google Scholar

Jjagwe, J., Olupot, P. W., Menya, E. & Kalibbala, H. M. Synthesis and application of Granular activated carbon from biomass waste materials for water treatment: A review. J. Bioresour. Bioprod. 6, 292–322 (2021).

Article CAS Google Scholar

Saleh, T. A. Characterization, determination and elimination technologies for sulfur from petroleum: Toward cleaner fuel and a safe environment. Trends Environ. Anal. Chem. 25, e00080 (2020).

Article CAS Google Scholar

Ickecan, D., Turkan, M. N. & Gulbicim, H. Investigation of shielding properties of impregnated activated carbon for gamma-rays. Appl. Radiat. Isotopes 172, 109687 (2021).

Article CAS Google Scholar

Ushiki, I. et al. Desorption behavior of various volatile organic compounds from activated carbon in supercritical carbon dioxide: Measurement and kinetic modeling. J. Supercrit. Fluids 121, 41–51 (2017).

Article CAS Google Scholar

Sayed, A.E.-R., Ashour, I. & Gadalla, M. Integrated process development for an optimum gas processing plant. Chem. Eng. Res. Des. 124, 114–123 (2017).

Article CAS Google Scholar

Rahimpour, M., Saidi, M. & Seifi, M. Improvement of natural gas dehydration performance by optimization of operating conditions: A case study in Sarkhun gas processing plant. J. Nat. Gas Sci. Eng. 15, 118–126 (2013).

Article CAS Google Scholar

San Miguel, G., Lambert, S. & Graham, N. Thermal regeneration of granular activated carbons using inert atmospheric conditions. Environ. Technol. 23, 1337–1346 (2002).

Article CAS Google Scholar

Dimpe, K. M., Ngila, J. & Nomngongo, P. N. Preparation and application of a tyre-based activated carbon solid phase extraction of heavy metals in wastewater samples. Phys. Chem. Earth Parts A/B/C 105, 161–169 (2018).

Article ADS Google Scholar

Lu, J. Y., Wang, S. F. & Zhu, J. F. Regeneration of granular activated carbon using dual-frequency ultrasound. Adv. Mater. Res 12, 1707–1710 (2013).

Google Scholar

Yang, D. & Du, X. A review about microwave regeneration technology of waste activated carbon. In IOP Conference Series: Earth and Environmental Science, 012101. IOP Publishing (2022)

McQuillan, R. V., Stevens, G. W. & Mumford, K. A. The electrochemical regeneration of granular activated carbons: a review. J. Hazard. Mater. 355, 34–49 (2018).

Article CAS Google Scholar

Yesilada, E., Tanaka, S., Tabata, M. & Sezik, E. The antiinflammatory activity of the fractions from Eryngium billardieri in mice. Phytother. Res. 3, 38–40 (1989).

Article Google Scholar

Sodeifian, G., Ardestani, N. S. & Sajadian, S. A. Extraction of seed oil from Diospyros lotus optimized using response surface methodology. J. For. Res. 30, 709–719 (2019).

Article CAS Google Scholar

Sodeifian, G., Ardestani, N. S., Sajadian, S. A. & Moghadamian, K. Properties of Portulaca oleracea seed oil via supercritical fluid extraction: experimental and optimization. J. Supercrit. Fluids 135, 34–44 (2018).

Article CAS Google Scholar

Sodeifian, G., Sajadian, S. A. & Honarvar, B. Mathematical modelling for extraction of oil from Dracocephalum kotschyi seeds in supercritical carbon dioxide. Nat. Prod. Res. 32, 795–803 (2018).

Article CAS Google Scholar

Sodeifian, G. & Sajadian, S. A. Investigation of essential oil extraction and antioxidant activity of Echinophora platyloba DC. using supercritical carbon dioxide. J. Supercrit. Fluids 121, 52–62 (2017).

Article CAS Google Scholar

Esfandiari, N. & Sajadian, S. A. CO2 utilization as gas antisolvent for the pharmaceutical micro and nanoparticle production: A review. Arab. J. Chem., 104164 (2022).

Sodeifian, G., Sajadian, S. A. & Ardestani, N. S. Experimental optimization and mathematical modeling of the supercritical fluid extraction of essential oil from Eryngium billardieri: Application of simulated annealing (SA) algorithm. J. Supercrit. Fluids 127, 146–157 (2017).

Article CAS Google Scholar

DeFilippi, R., Kyukonis, V., Robey, R. & Modell, M. Supercritical fluid regeneration of activated carbon for adsorption of pesticides (1980).

Kander, R., Paulaitis, M., Paulaitis, M., Penninger, J., Gray, R., & Davidson, R. Chemical Engineering and Supercritical Conditions, Ann Arbor, Michigan (1983).

Tan, C. S. & Liou, D. C. Desorption of ethyl acetate from activated carbon by supercritical carbon dioxide. Ind. Eng. Chem. Res. 27, 988–991 (1988).

Article CAS Google Scholar

Madras, G., Erkey, C. & Akgerman, A. Supercritical fluid regeneration of activated carbon loaded with heavy molecular weight organics. Ind. Eng. Chem. Res. 32, 1163–1168 (1993).

Article CAS Google Scholar

Macnaughton, S. J. & Foster, N. R. Supercritical adsorption and desorption behavior of DDT on activated carbon using carbon dioxide. Ind. Eng. Chem. Res. 34, 275–282 (1995).

Article CAS Google Scholar

Harikrishnan, R., Srinivasan, M. & Ching, C. Adsorption of ethyl benzene on activated carbon from supercritical CO2. Am. Inst. Chem. Eng. AIChE J. 44, 2620 (1998).

Article CAS Google Scholar

Lucas, S., Cocero, M. J., Zetzl, C. & Brunner, G. Adsorption isotherms for ethylacetate and furfural on activated carbon from supercritical carbon dioxide. Fluid Phase Equilib. 219, 171–179 (2004).

Article CAS Google Scholar

Salvador, F., Martin-Sanchez, N., Sanchez-Montero, M. J., Montero, J. & Izquierdo, C. Regeneration of activated carbons contaminated by phenol using supercritical water. J. Supercrit. Fluids 74, 1–7 (2013).

Article CAS Google Scholar

Khajeh, M. Optimization of process variables for essential oil components from Satureja hortensis by supercritical fluid extraction using Box-Behnken experimental design. J. Supercrit. Fluids 55, 944–948 (2011).

Article CAS Google Scholar

Sodeifian, G., Sajadian, S. A. & Ardestani, N. S. Extraction of Dracocephalum kotschyi Boiss using supercritical carbon dioxide: Experimental and optimization. J. Supercrit. Fluids. 107, 137–144 (2016).

Article CAS Google Scholar

Ballard, D. How to operate a glycol plant. Hydrocarb. Process. 45, 171–180 (1966).

CAS Google Scholar

Honarvar, B., Sajadian, S., Khorram, M. & Samimi, A. Mathematical modeling of supercritical fluid extraction of oil from canola and sesame seeds. Braz. J. Chem. Eng. 30, 159–166 (2013).

Article CAS Google Scholar

D. Astm, 4607-94, Standard test method for determination of iodine number of activated carbon (2006).

Nunes, C. A. & Guerreiro, M. C. Estimation of surface area and pore volume of activated carbons by methylene blue and iodine numbers. Quim. Nova 34, 472–476 (2011).

Article CAS Google Scholar

Aboosadi, Z. A., Jahanmiri, A. & Rahimpour, M. Optimization of tri-reformer reactor to produce synthesis gas for methanol production using differential evolution (DE) method. Appl. Energy 88, 2691–2701 (2011).

Article Google Scholar

Khajeh, S., Aboosadi, Z. A. & Honarvar, B. A comparative study between operability of fluidized-bed and fixed-bed reactors to produce synthesis gas through tri-reforming. J. Nat. Gas Sci. Eng. 19, 152–160 (2014).

Article CAS Google Scholar

Rahimpour, M. R., Aboosadi, Z. A. & Jahanmiri, A. H. Differential evolution (DE) strategy for optimization of methane steam reforming and hydrogenation of nitrobenzene in a hydrogen perm-selective membrane thermally coupled reactor. Int. J. Energy Res. 37, 868–878 (2013).

Article CAS Google Scholar

Rai, A., Mohanty, B. & Bhargava, R. Modeling and response surface analysis of supercritical extraction of watermelon seed oil using carbon dioxide. Sep. Purif. Technol. 141, 354–365 (2015).

Article CAS Google Scholar

Maran, J. P. & Priya, B. Supercritical fluid extraction of oil from muskmelon (Cucumis melo) seeds. J. Taiwan Inst. Chem. Eng. 47, 71–78 (2015).

Article CAS Google Scholar

Takahashi, N. et al. Measurement and prediction of desorption behavior of five volatile organic compounds (acetone, n-hexane, methanol, toluene, and n-decane) from activated carbon for supercritical carbon dioxide regeneration. J. Supercrit. Fluids 107, 226–233 (2016).

Article CAS Google Scholar

Salimi, A., Fatemi, S., Zakizadeh Nei Nei, H. & Safaralie, A. Mathematical Modeling of Supercritical Extraction of Valerenic Acid from Valeriana officinalis L. Chem. Eng. Technol. 31, 1470–1480 (2008).

Article CAS Google Scholar

Döker, O., Salgin, U., Yildiz, N., Aydoğmuş, M. & Çalimli, A. Extraction of sesame seed oil using supercritical CO 2 and mathematical modeling. J. Food Eng. 97, 360–366 (2010).

Article Google Scholar

Sodeifian, G., Ghorbandoost, S., Sajadian, S. A. & Ardestani, N. S. Extraction of oil from Pistacia khinjuk using supercritical carbon dioxide: Experimental and modeling. J. Supercrit. Fluids 110, 265–274 (2016).

Article CAS Google Scholar

Han, X., Cheng, L., Zhang, R. & Bi, J. Extraction of safflower seed oil by supercritical CO 2. J. Food Eng. 92, 370–376 (2009).

Article CAS Google Scholar

Danielski, L. et al. Marigold (Calendula officinalis L.) oleoresin: solubility in SC-CO 2 and composition profile. Chem. Eng. Process. Process Intensif. 46, 99–106 (2007).

Article CAS Google Scholar

Tan, C.-S. & Liou, D.-C. Modeling of desorption at supercritical conditions. AIChE J. 35, 1029–1031 (1989).

Article CAS Google Scholar

Sodeifian, G., Sajadian, S. A. & Ardestani, N. S. Evaluation of the response surface and hybrid artificial neural network-genetic algorithm methodologies to determine extraction yield of Ferulago angulata through supercritical fluid. J. Taiwan Inst. Chem. Eng. 60, 165–173 (2016).

Article CAS Google Scholar

Crank, J. The Mathematics of Diffusion (Oxford University Press, 1979).

MATH Google Scholar

Goto, M., Roy, B. C., Kodama, A. & Hirose, T. Modeling supercritical fluid extraction process involving solute-solid interaction. J. Chem. Eng. Jpn. 31, 171–177 (1998).

Article CAS Google Scholar

Sovová, H. Mathematical model for supercritical fluid extraction of natural products and extraction curve evaluation. J. Supercrit. Fluids 33, 35–52 (2005).

Article Google Scholar

Putra, N., Yian, L., Nasir, H., Idham, Z.B. & Yunus, M. Effects of process parameters on peanut skins extract and CO2 diffusivity by supercritical fluid extraction. In IOP Conference Series: Materials Science and Engineering, p. 012057. IOP Publishing (2018).

Download references

The authors gratefully acknowledge support from the Research Council of the University of Islamic Azad Marvdasht. The authors also gratefully acknowledge the financial support from South Zagros Oil & Gas Production Co. Furthermore, the authors would like to appreciate Senior Engineer Mr. R. Ghobadi, CEO of Central workshop, for technical supports.

This research received no external funding.

Department of Chemical Engineering, Islamic Azad University of Marvdasht, Marvdasht, Iran

Mohammad Najafi, Zahra Arab Abousadi & Bizhan Honarvar

South Zagros Oil and Gas Production, National Iranian Oil Company, Shiraz, 7135717991, Iran

Mohammad Najafi & Seyed Ali Sajadian

Department of Chemical Engineering, Faculty of Engineering, University of Kashan, Kashan, 87317-53153, Iran

Seyed Ali Sajadian

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Conceptualization, S.A.S., and M.N.; methodology, M.N., S.A.S., B.H., Z.A.A., validation, S.A.S., Z.A.A., B.H., and M.N.; writing—original draft preparation, M.N.., S.A.S., and Z.A.A.; writing—review and editing, S.A.S., Z.A.A., and B.H.; visualization, S.A.S., and B. H.; supervision, B.H., Z.A.A.; project administration, S.A.S. and Z.A.A.; funding acquisition, and B.H. All authors have read and agreed to the published version of the manuscript.

Correspondence to Seyed Ali Sajadian.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original online version of this Article was revised: The original version of this Article contained an error in Reference 22, which was incorrectly given as “Sodeifian, G. & Sajadian, S. A. CO2 utilization as gas antisolvent for the pharmaceutical micro and nanoparticle production: A review. Arab. J. Chem., 104164 (2016).” The correct Reference is: “Esfandiari, N. & Sajadian, S. A. CO2 utilization as gas antisolvent for the pharmaceutical micro and nanoparticle production: A review. Arab. J. Chem., 104164 (2022).”

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

Najafi, M., Abousadi, Z.A., Honarvar, B. et al. Separation of hydrocarbons from activated carbon as a porous substance in a glycol regeneration process using supercritical carbon dioxide. Sci Rep 12, 19910 (2022).

Download citation

Received: 12 September 2022

Accepted: 03 November 2022

Published: 19 November 2022


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.