Next Article in Journal
Dual-Function Photocatalysis in the Visible Spectrum: Ag-G-TiO2 for Simultaneous Dye Wastewater Degradation and Hydrogen Production
Next Article in Special Issue
Nano-TiO2-Enhanced Surface Functionalization of Recycled Concrete Aggregates for Improved Degradation Efficiency of Low-Concentration Sulfur Dioxide
Previous Article in Journal
Defective TiO2/MIL-88B(Fe) Photocatalyst for Tetracycline Degradation: Characterization and Augmented Photocatalytic Efficiency
Previous Article in Special Issue
Advanced Technologies Conciliating Desulfurization and Denitrogenation to Prepare Clean Fuels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Production of Green Fuel Using a New Synthetic Magnetite Mesoporous Nano-Silica Composite Catalyst for Oxidative Desulfurization: Experiments and Process Modeling

by
Aysar T. Jarullah
1,
Ahmed K. Hussein
1,
Ban A. Al-Tabbakh
2,
Shymaa A. Hameed
1,
Iqbal M. Mujtaba
3,*,
Liqaa I. Saeed
4 and
Jasim I. Humadi
5
1
Chemical Engineering Department, College of Engineering, Tikrit University, Tikrit 34001, Iraq
2
Petroleum Research & Development Center, The Iraqi Ministry of Oil, Baghdad 00964, Iraq
3
Chemical Engineering Department, Faculty of Engineering & Informatics, University of Bradford, Bradford BD7 1DP, UK
4
Department of Mining Engineering, College of Petroleum and Mining Engineering, University of Mosul, Mosul 41002, Iraq
5
Department of Petroleum and Gas Refining Engineering, College of Petroleum Processes Engineering, Tikrit University, Tikrit 34001, Iraq
*
Author to whom correspondence should be addressed.
Catalysts 2024, 14(8), 529; https://doi.org/10.3390/catal14080529
Submission received: 10 June 2024 / Revised: 7 August 2024 / Accepted: 8 August 2024 / Published: 15 August 2024

Abstract

:
Producing an eco-friendly fuel with the least amount of sulfur compounds has been an ongoing issue for petroleum refineries. In this study, bentonite (which is a cheap material and is locally available in abundance) is employed to prepare a nano-silica catalyst with a high surface area to be used for the oxidative desulfurization of kerosene. Two composite catalysts of Fe/silica were supported on CAT-1 (0% HY-zeolite and 100% nano-silica) and CAT-2 (20% HY-zeolite and 80% nano-silica). The activity of the catalysts was evaluated in a batch ODS (oxidative desulfurization) process at temperatures of 30, 60, 90, and 120 °C, a pressure of 1 atm, and a reaction time of 30, 60, 90, and 120 min using 120 L/h of air as the oxidant. The results revealed that the highest total sulfur removal efficiency was 50% and 87.88% for 100% nano-silica (CAT-1) and 80% nano-silica (CAT-2), respectively. The experimental data were then used to construct and validate an accurate mathematical model of the process. The operational parameters for eliminating more than 99% of sulfur and producing eco-friendly fuel were then achieved by using the model. The testing methods for these characterizing materials included X-ray diffraction (XRD), thermal gravimetric examination (TGA), X-ray fluorescence (XRF), and surface area (BET). The outcomes indicated that the addition of HY-zeolite increased the activity of the catalyst (CAT-2 > CAT-1).

1. Introduction

Currently, the production of eco-friendly fuel to combat climate change is a major concern worldwide. Globally, 82% of energy is provided by fossil fuels, but it has been noticed that the use of fossil fuels has declined over the years [1]. The main source of energy used globally are petroleum fractions and their main usage is in the transportation sector, including jet fuel, diesel, and gasoline. The removal of toxic components such as sulfur [2] from the fuel that is widely used in transportation and in the petrochemical industry is absolutely essential from an environmental point of view. Fuels containing sulfur, like Benzothiophenes (BTs), Thiophenes (Ts), Dibenzothiophenes (DBTs), and Mercaptans (RSH), can be subjected to oxidative desulfurization, resulting in the production of sulfur-di-oxide, which can be easily removed from the fuel, thus producing clan fuel [3,4]. This process uses various types of oxidants, viz., air, H2O2, O2, and O3. Air is utilized due to its accessibility; it is inexpensive and causes no further pollution [5,6]. However, the choice of catalyst for sulfur compounds is considered the most critical issue in designing an efficient ODS process. Thus, researchers have continued to develop different types of catalysts for more rapid, economical, and efficient ODS processes in order to select the best type for the petroleum industry [7,8]. Magnetite mesoporous is an excellent catalyst, as it has a tremendous capacity to eliminate methyl mercaptan due to the robust molecular affinity of magnetite and the huge surface area and pore dimensions of mesoporous silica [9,10]. In this study, the authors have prepared a novel catalyst based on iron oxide magnetite supported on nano-silica and zeolite nanoparticles. The use of such catalysts has not yet entered the public domain. The primary objective of this study is to produce an eco-friendly fuel by eliminating sulfur compounds from kerosene in a batch reactor using the ODS technique. Another goal of this study is to create a mathematical model of the process that has excellent precision and can be optimized and simulated, providing in-depth understanding of the process. Experimental data and model predictions will be used within an optimization framework to estimate unknown model parameters (such as kinetic parameters).

2. Results and Discussion

2.1. Evaluation of Prepared Catalysts

2.1.1. Chemical Composition

X-Ray fluorescence (XRF) was used to determine the existence of metal oxides, the quantity of active metal (Fe2O3) and nano-silica support, and the structural composition. Table 1 demonstrates the synthesized nano-silica’s purity, and Table 2 represents the XRF analysis of the prepared catalyst.

2.1.2. X-ray Diffraction

The oldest and most widely used technique for characterizing catalysts is thought to be X-ray diffraction. Lattice structural parameters are used to identify the crystalline phase within catalysts and to determine the size of the particles. The X-ray diffraction instrument used in this study was capable of generating Ni-filtered cuK radiation with a wavelength of 1.5418 Å. The instrument operated at 40 kV and 30 mA during the XRD measurements of the catalyst samples. The diffraction patterns were recorded over a 2θ angular range of 3° to 60°. X-ray diffraction is a well-established technique that can distinguish materials based on their characteristic diffraction patterns, with a scanning rate of 5°/min. It offers statistics on how far a specific structure deviates from the overall structure, in addition to section identification. High-resolution X-ray diffraction was used to record the crystal structure of the particles. Additionally, the presence of an amorphous silica component was confirmed by a large diffraction peak at 2 = 20°–24° [11]. The clearly defined nano-silica diffraction peak revealed the hypothesized existence of a highly structured hexagonal mesostructure. Figure 1 below illustrates the XRD for amorphous nano-silica.

2.1.3. Fourier Transform Infrared (FTIR)

A PerkinElmer FTIR spectrometer was utilized to obtain the Fourier transform infrared spectroscopy (FTIR) results depicted in Figure 2. The peaks observed at 1087.85 and 796.90 cm−1 indicate the presence of asymmetric and symmetric stretching modes of the Si-O-Si bonds. Furthermore, a band centered at 462 cm−1 was detected, which is attributed to the Si-O-bending vibration of the sample containing iron and no other elements aside from silicon and oxygen. Additionally, features at 92 cm−1 are noted at DLA. The broad band around approximately 3481 cm−1 is attributed to the increasing vibrational frequencies of the silanol group and the remaining absorbed water.
Notably, no peaks are observed in the 2800–3000 cm−1 range in the spectra, indicating the absence of significant vibrational modes related to the water–silica matrix interactions. This actually implies that there were no actual natural mixes in the silica after controlled calcination. The top positions of the collected nano-silica showed no notable changes in the FTIR continuous flow. Silica with a large surface area is the best multipart for use as a catalyst backing. Additionally, it has been discovered that the surface area is related to the silica reactivity.
Spectroscopic investigation provided unequivocal evidence that the mesopore walls were made of a zeolite-like structure [12,13,14]. The scientific literature frequently cites bands in the range of 550 to 600 cm−1 when investigating materials containing MFI-type and zeolite β seed structures. These bands are believed to correspond to the asymmetric stretching vibration of the 5-membered ring components within the material’s framework. The tetrahedral coordination of Al is a further intriguing characteristic of the aluminosilicate phase. This is significant because it emphasizes the existence of a potent Bronsted acid center of the zeolite variety, which is useful for research as a microporous material, particularly as a possible heterogeneous catalyst [13,14,15,16,17]. Figure 3 shows the results of the produced catalyst’s FTIR examination.

2.1.4. TGA (Thermal Gravimetric Analysis)

A TGA Model SDT 2960 device, owned by the Oil Research and Development Center located in Baghdad, Iraq, was utilized to examine the silica particles. The tests were conducted at a heating rate of 5 °C/min in a nitrogen environment from ambient to 900 °C. TGA profiles for the nano-silica and composite samples are shown in Figure 4A,B. It can evidently be seen from these figures that the mass losses of both samples were the same before and after the process of loading active metals. The curves display two areas of mass loss, the first of which (between 20 and 120 °C), with mass loss values of around 10%, is caused by the evaporation of water that has been physically adsorbed on the solid. In the TGA measurements, mass loss values of 16.74% for the support’s nano-silica and 16.79% for the catalysts were noted. Such favorable outcomes demonstrated the strong thermal stability of the supports and of the manufactured catalysts.

2.1.5. Particle Size Distribution

The connection between an assumed particle dimension and the frequency of that size or diameter is known as “particle size distribution”. Atomic Force Microscopy (AFM) is conducted by ethanol dispersion to measure the particle size distribution. The basis of this discussion rests on the perceived notion that the majority of particles exhibit remarkably similar shapes. It is crucial, when evaluating the catalytic properties of metal-backed catalysts, to account for the actual sizes and distributions of the particles. The physical dimensions of the metal particles determine their surface area, and consequently the potential number of active sites, as well as the electrical characteristics of the metal, which are dependent on particle size. Even though such catalysts may share similar average particle sizes, the distributions of particle extents can vary dramatically, and the characteristics of the compared catalysts may differ significantly as a result. Figure 5 explains the particle size distributions of nano-silica and the produced catalysts.

2.1.6. FESEM-EDX

The SiO2 nanocatalyst’s surface shape and real particle size were visible in SEM, in large part due to the catalyst’s enormous surface area and mesoporous frequencies functioning as beds and supports. This technique was utilized using a Zeiss-EM10C-100KV FESEM device (Oberkochen, Germany). Figure 6 below displays a SEM picture of the nano-silica.
Figure 7 displays the SiO2 nanocatalyst’s surface elemental composition as determined by EDX analysis. Si and O are present in every sample. The map displays the relative dispersion density and the EDX element map. The distribution and consistency of the dispersed particles can be examined through an EDX (energy-dispersive X-ray) image of the catalytic converter’s external surface. The EDX mapping analysis revealed the presence of the following elements: silicon (37.6% by weight), oxygen (32% by weight), carbon (8.7% by weight), sodium (5.7% by weight), iron (5.4% by weight), and aluminum (10.6% by weight). Figure 7, which depicts the EDX image for the nano-silica generated in this work, illustrates these details.
The bimetal nanocatalyst’s (Fe3O4/SiO2) surface shape and particle size were demonstrated by FESEM. SiO2 acted as a couch support for Fe3O4 and zeolite due to its variable surface and mesoporous frequencies. As seen in the image, several forms of iron were immobilized on a mesoporous silica surface. The SEM image of the catalyst (Fe/nano-silica) created in our work is given in Figure 8 below. Additionally, the micrographs (magnification of 1 m) demonstrate a uniform 100 distribution of the active metal (Fe) on the catalyst’s surface along with the presence of spherical particles and well-developed pores.
The scanning electron micrograph in Figure 9 depicts the Fe/SiO2 nanocatalyst, and the corresponding energy-dispersive X-ray (EDX) analysis shows its surface component composition. The EDX results presented in Figure 9 characterize the sample compositions, which contain Si, O, and Fe components. As the metal oxide loading increased, the surface Fe concentration rose, while the Si concentration decreased, indicating a more uniform distribution of Fe on the SiO2 nano-silica surface. The relative dispersion density and EDX elemental map are shown in the provided figure. When a sample was taken from the exterior surface of the catalytic converters and analyzed using EDX, the detected elements were found to be distributed uniformly across the converter’s outer surface. The EDX mapping reveals a composition of C (5% wt), O (34.2% wt), Na (4.7% wt), Si (34.1% wt), Fe (9.6% wt), and Al (12.4% wt).

2.2. Experimental Section

A lab-synthesized catalytic system with nanoparticles on zeolite was tested in a batch reactor. The influence of reaction temperature, reaction duration, and initial sulfur species concentration was investigated.

2.2.1. The Influence of Reaction Temperature

The reaction temperature investigated in the transformation of ODS reactions was 303 K, 333 K, 363 K, and 393 K. The results are illustrated in Figure 10.
The observed increase in reaction temperature appears to coincide with an enhanced conversion of sulfur-containing compounds. This observation can be attributed to the following factors:
  • The higher temperature favorably influenced the reaction rate, as described by the Arrhenius equation. This leads to a notable rise in the conversion of sulfur-based compounds, as reported in the literature [18].
K = K o e ( E A R T   )
In addition to the activation energy, the reactant temperature affects both the reaction rate constant and the Arrhenius equation in an inverse and direct proportion, respectively.
Higher temperatures increase the number of molecules with sufficient activation energy, facilitating the reaction between the oxidant and organic compound, thereby improving conversion. This behavior is consistent with the literature [19], suggesting that temperature can influence the oxidation reaction rate based on the activation energy.
  • Furthermore, temperature can affect various transport and physical properties of the gases and liquids involved, such as density, viscosity, gas diffusivity, and Henry’s constant. Specifically, an increase in temperature leads to an increase in both Henry’s constant and the diffusivity of reactants (sulfur and oxidant in the presence of catalyst) in the reaction media. Additionally, temperature impacts properties like viscosity, which decrease with increasing temperature. However, the oxidation reaction is favored with higher temperatures, as the absorption of molecular oxygen in the liquid, the dispersion of sulfur compounds, and the rate of dissolved oxygen reaching active sites inside the catalyst pores all increase with rising temperatures.
  • The higher temperatures can be explained by the fact that faster-moving reactants have a higher probability of interacting. Moreover, for the catalyst used, the finest sulfones may have been readily desorbed from the surface at high temperatures [20].
  • The catalysts’ oxidation temperature has a substantial effect on their activity, and the density of acidic centers, which determine the catalytic properties, depends on the temperature at which the catalyst has been calcined. This has led to enhanced removal of sulfur compounds from diesel fuel [21,22].

2.2.2. The Impact of Reaction Time

This study investigated an oxidation reaction at various time intervals, including 30 min, 60 min, 90 min, and 120 min. The results were presented graphically, and it was observed that the experimental findings for different temperatures were unique for each catalyst. In all cases, the desulfurization efficiency increased as the reaction time was extended, as a longer reaction time allowed for enhanced contact between the reactants. The interpretation of the figures can be summarized as follows: Figure 11A indicates that an increase in reaction time led to an increase in the conversion rate. With extended reaction time, the reactants had more opportunity to interact with the active sites of the catalyst, thereby enhancing the overall interaction process [23,24,25]. For instance, by simply increasing the reaction time from 30 min to 120 min at 303 K, the conversion of sulfur compounds rose from 33% to 79.83%. Images presenting the obtained results are provided in Figure 11 below.

2.3. Estimation of Kinetic Parameters

Kinetic parameter values are obtained when the model equations produce the closest correspondence to the experimental data [26,27,28]. In this study, the mathematical model was leveraged to estimate the kinetic parameters of the ODS process under various operational conditions and sulfur concentrations, where the maximum correlation with the measured data was achieved. Two methodologies were considered: linear growth and non-linear growth. The model takes into account the constant parameters detailed in Table 3 below.

2.3.1. Linear Approach

The optimal kinetic parameter values from the linear regression approach, presented in Table 4 and Table 5, enable absolute errors of less than 5%.
As can be seen in Table 4 and Table 5, the sequence of the reaction varies depending on the type of catalyst used. For CAT-1 and CAT-2, the ODS reaction order for sulfur compounds in kerosene was found to be 2 and 1.9860986, respectively.

Activation Energy (E)

Figure 12 and Figure 13 depict the Arrhenius equation, which was used to determine the activation energy for each catalyst by plotting (lnk) vs (1/T) and finding that the slope of the line is equal to (−E/R). Table 6 displays the calculated activation energies.

2.3.2. Approach with Non-linear Regression Method

This study simultaneously assessed the reaction order, activation energy, and pre-exponential factor, with the results shown in Table 7 and Table 8 below.
These tables show that, under the same conditions, CAT-2 provides a higher conversion rate than the other catalysts. Under these conditions (reaction temperature = 393 K, reaction time = 120 min), CAT-1 converted sulfur compounds at a rate of 74.77%, while CAT-2 achieved a rate of 85.35%. The alternative catalysts demonstrate greater activation energy requirements for the reaction compared to CAT-2. This suggests that the use of CAT-2 leads to a reduction in the activation energy, which in turn correlates with an increase in the reaction rate. These parameters were estimated using both linear and non-linear methods, with the latter yielding lower sums of squared errors (SSEs). Given this, it follows that the second method yields more precise parameter estimation, at least to some extent, based on (SSE).
In comparison to the concurrent determination of kinetic parameters through non-linear regression, the magnitude of error introduced during the computation of activation energy and pre-exponential factor via the linearization of the Arrhenius equation is notably higher. The final sulfur concentration and reaction rate models obtained by non-linear methods for CAT-1 and CAT-2 are as follows:

For CAT-1

  • Reaction rate:
= r s u l f u r = 578.64 ƞ 0   e ( 4126.2184 T )   C s u l f u r 2
  • Concentration profile:
  C s u l f u r = [ 1.33155 + 578.64 ƞ 0 e ( 4126.2184 T ) . t ] 1

For CAT-2

  • Reaction rate:
r s u l f u r = 1681.534 ƞ 0   e ( 4182.68994 T )   C s u l f u r 1.975313
  • Concentration profile:
  C s u l f u r = [ 1.32217 + 1640.021971 ƞ 0 e ( 4182.68994 T ) . t ] 1.025311

2.4. Experimental Findings and Simulation Outcomes

The oxidation process was simulated after the optimal kinetic parameters had been estimated. Table 9, Table 10, Table 11 and Table 12 below show the predicted results compared to the experimental findings under a variety of operational scenarios.
The results in Table 13 show that for the same number of training data points and depending on the target function, the estimations of the parameters provided by the second approach were definitely more accurate than the ones given by the first approach. It was because of this characteristic that, for the same number of training data points, the second approach actually gave rise to a smaller sum of squared errors (SSEs). The experimental errors obtained through a linearization (first method) of the Arrhenius equation to find activation energy (E) and pre-exponential factor were far above those for the non-linear fit of data (second method).
The predicted sulfur conversion outcomes are depicted in Figure 14 and Figure 15, alongside the experimental results. The figure displays experimental values (x-axis) versus observed values (y-axis) under the same conditions, demonstrating reasonable agreement between the experimental data and the simulations. As can be observed from the figures, the predicted and experimental data exhibit a high degree of conformity, with the ratio between them being close to a linear relationship with a slope value near one.

2.5. Optimal Conversion to Eco-Friendly Fuel

Using best-fit kinetic parameters, Table 14 and Table 15 show the best operating conditions for each catalyst.
The utilization of CAT-2 under particular conditions, which are 550 K reaction temperature, 193.63318 min contact time, and an initial sulfur compound concentration of 0.2443 (wt%), meet environmental standards and attain a low sulfur content, resulting in high-quality fuel oil, as shown in Table 12 and Table 13.
Maximum conversion (almost 99%) was achieved using the optimal operating conditions, as demonstrated in this section. The following is a brief synopsis of the prerequisites that need to be met in order to have a product with a low level of sulfur dioxide, which is better for the environment.
Temperature: As temperature rises, the sulfur compound is more likely to be converted. This is because the reaction rate constant increases with temperature, as described by the Arrhenius equation. Elevating the temperature also enhances other critical factors like diffusivity and viscosity. Consequently, the sulfur compound’s diffusion rate, air molecule penetration through catalyst pores, and air absorption into kerosene all increase with higher temperature.
Time: As the reaction goes on for longer, there will be more time for the reactants to make contact with one another on the active sites of the catalyst. As a direct consequence of the prolonged residence time, a significant amount of reaction product is produced.
Table 16 presents a comparison between the present study and some previous studies.

3. Experimental Methodology

This study investigated the oxidation of sulfur in kerosene using a novel composite nanocatalyst and atmospheric oxygen as the oxidizing agent. Four distinct catalysts were prepared following a series of steps.
(a)
Nano-silica was made using bentonite.
(b)
Using the precipitation process, two composite supports (nano-silica + HY-zeolite) were prepared with varying percentages of HY-zeolite: support-1: nanoparticles of silica only; support-2: nanoparticle composed of 80% nano-silica and 20% HY-zeolite.
(c)
Using the impregnation process, two types of nanocatalysts were created:
  • Fe2O3/nano-silica nanoparticle (CAT-1).
  • Fe2O3/(80% nano-silica + 20% HY-zeolite) nanoparticle (CAT-2).
The efficiency of the two prepared catalysts was evaluated in a batch reactor undergoing oxidative desulfurization at different operational conditions, shown in Table 17.

3.1. Materials and Equipment

3.1.1. Kerosene Feedstock (Fuel Feed)

Kerosene was utilized as the fuel feedstock here. The consumption of this type of fuel is significantly high in daily life, and in aircraft, its consumption is rapidly increasing. Commercial kerosene, with a total sulfur content equal to 0.1966 wt. %, was provided by the Baiji refinery, North Refineries Company. The detailed specification of the feedstock is shown in Table 18.

3.1.2. Oxidant

Utilizing a compressor, air was delivered primarily as an oxidizing mediator. Sulfones were created by oxidizing the sulfur compound with oxygen found in the air.

3.1.3. Active Components for Catalyst Preparation

The chemical compounds that were utilized to provide catalysts for the active components in the present research study are listed in Table 19.

3.1.4. HY-Zeolite

HY-zeolite was synthesized and characterized by the Iraqi Ministry of Oil’s Laboratory. The particle properties and XRF data are shown in Table 20 and Table 21, respectively.

3.1.5. Equipment

Table 22 shows the required equipment used in the present study.

3.2. Nano-Silica Preparation

Mesoporous silicate/aluminosilicate materials (M41S phase) were synthesized in 1992 with a pore size range of 2 to 30 nm. These materials belong to the mesoporous silica molecular sieves (MCM) series—MCM-41, featuring a hexagonal distribution of mesopores; MCM-48, cubic; and MCM-50, layered [32]. The complicated steric arrangements of these materials, which are well-organized mesoporous solids, have encouraged scientists to develop new methods for synthesizing new types that have unconventional physical and chemical characteristics. Therefore, the synthesis of zeolite from these groups maintains the mesopores’ construction integrity.
In other words, the optimal porous material comprises small, clean pores surrounded with crystalline walls, in such a way that the mass transfer is carried out with the maximum possible efficiency, and thermal and hydrothermal stability are coupled with better control of its active properties. The compounds and materials required for the precipitation process—which yields nano-silica—are listed in Table 23.
The sand was mixed with sodium hydroxide, which led to the creation of very small (around 45 nanometer) particles. This mixture was then stirred together with bentonite sand and heated in a programmable electric furnace at 550 degrees Celsius for 30 min. Water was added to the heated mixture, forming a consistent, transparent liquid. Concentrated sulfuric acid (98%) was slowly introduced into this liquid until the pH reached 3. The gel that formed was separated by filtering, and then it was rinsed once with filter paper. The filtered gel was then washed with deionized water (pH 7) to remove sulfate ions. Finally, the gel was dried overnight at 110 degrees Celsius, producing the material. The overall process of preparing this nano-silica support material is summarized in Figure 16.

3.3. Preparation of Iron-Supported Nano-Silica Catalyst

Incorporating different elements into microporous and mesoporous materials has been an active area of research. It is driven by the need to develop more stable and effective materials for applications like catalysis, separation, coatings, and sensing [33,34]. The discovery of the useful catalytic properties of titanium silicate molecular sieves sparked interest in incorporating other metals into similar porous materials [34]. Iron-containing zeolites and related microporous substances were found to have unique catalytic abilities applicable to gas-phase reactions like hydrocarbon oxidation, ammonia decomposition, nitrous oxide reduction, and selective catalytic reduction of nitric oxide and nitrous oxide [35,36]. In this study, the researchers thoroughly mixed 3 g of ferric nitrate nonahydrate with 75 milliliters of water for 30 min using an ultrasonic mixer. Then, 22 g of nano-silica was slowly added to the mixture, which was further stirred in an ultrasonic water bath for 60 min to ensure thorough mixing. The resulting gel was filtered using a Buchner funnel and a vacuum pump, and the dry sample was dried overnight at 110 °C. To prepare the catalyst, a calcination process was carried out at 550 °C for 3 h. The catalyst preparation process and the loading of the active metal are illustrated in Figure 17.
The catalyst was tested at the Petroleum Research Center in Baghdad, and its physical properties were examined and are shown in Table 24.

3.4. Oxidative Desulfurization Reactor

3.4.1. Batch Reactor

Oxidative desulfurization was carried out in a 250 mL three-necked round-bottom flask. One neck connects to a condenser that separates kerosene vapor from air. Another has an air inlet connected to a compressor to force air into the flask. The third neck allows for temperature monitoring and sampling. The contents are heated and stirred using a magnetic stirrer. A diagram is shown in Figure 18.

3.4.2. The ODS Reaction

The following section presents the likely mechanism behind the oxidation of sulfur compounds during the heterogeneous desulfurization of fuels using a catalyst. The complete reaction of converting ethyl mercaptan to a disulfide is achieved with the aid of 5% iron core/(silica + zeolite) nanoparticles. The reactions start with the metals reacting with oxygen in the air to form peroxide species [37,38].
The peroxide compounds incorporate the sulfur atom and react with the sulfur atoms in ethyl mercaptan. This causes the sulfur atoms to be oxidized, forming a disulfide compound. At the same time, the primary metal oxide is produced as a byproduct of this reaction. The resulting disulfide can then be easily extracted using a polar solvent, which is a suitable method for driving the reaction [37,39,40]. The overall reaction can be summarized as follows:
2RSH + O2 → 2RSSR + H2O

3.5. Experimental Testing

The experimental runs for the ODS process of kerosene fuel were conducted under different operating conditions, as shown in Table 17.
A controlled airflow of 120 L/h was adjusted to obtain a pressure of 1 atm. The points taken in the experiment were carefully selected so that all maxima and minima were considered while also taking into account safe and moderate conditions in between.

Experimental Procedure for ODS Reactions

To prepare for the experiments, sulfur-supplying feedstock oil, specifically kerosene, was used. This included the preparation of the necessary catalysts and other equipment required. During each experimental run, the following steps were carried out:
  • A round-bottom flask was charged with 80 mL of the feedstock before proceeding with other steps.
  • The flask was placed in a heating mantle stirrer, and an air hose was connected to a condenser as well.
  • Measures were taken to ensure that heat was not transferred to the condenser through the kerosene vapor, and the condenser was supplied with cooled water.
  • A thermometer measured the temperature at which the overall reaction occurred.
The catalyst was added to the reactor once the required temperature was reached, and the time was recorded. The sulfur content of the final product was then analyzed to determine the reaction rate.

3.6. Analysis of Liquid Samples

The ASTM-D4294 method combined with X-ray fluorescence was used to analyze the overall content of sulfur in the oil feedstock and the final products [40]. All the available analytical methods were used to characterize the properties of the feedstock. In every operational condition, the analysis of products was duplicated for each sample so that the consistency and validation of the results can be demonstrated. The results obtained from each run were taken as average and a maximum permissible variance of 0.5 percent was assumed for all runs.

4. Mathematical Modeling of ODS Process

A mathematical model was developed for the ODS process based on the following assumptions:
  • The reactor is isothermal at a constant pressure.
  • A reactant that is in excess is in the gaseous state.
  • Reactor design produces a homogeneous environment throughout.

4.1. Mass Balance Equation

The general mass balance equation for the sulfur component of the catalytic batch reactor at an increment of time (dt) is as follows:
[input] = [output] + [consumption by reaction] + [accumulation]
  • Input of sulfur, mass/time = 0.
  • Output of sulfur, mass/time = 0.
  • Consumption of sulfur by reaction, mass = −(rRS)V dt.
Accumulation of sulfur, mass = V d C R S d t dt.
Introducing these terms into Equation (6) will give:
0 = 0 + ( r RS )   Vdt + V d C R S d t dt
Then,
( r RS ) = d C R S d t
So:
dt = d C R S ( r R S )
where:
  • at t = 0, CRS = CRSO
  • t = t, CRS = CRS
Thus:
0 t d t   = C R S O C R S d C R S ( r R S )
Finally, the design correlation for the ideal batch reactor is written in terms of batch time to be
t = C R S O C R S d C R S ( r R S )
Equation (9) can be integrated to obtain the concentration profile of sulfur compounds in the catalytic reactor, although the term (−rRS) is primarily dependent on the sulfur concentration.

4.2. Reaction Rate

The oxidation of sulfur by air is a complex multi-step process. The reaction kinetics can be described using an nth-order rate equation to capture the complexity.
(−rRS) = kapp CRSn
The apparent kinetic constant kapp is controlled by the intrinsic kinetic constant kin, modified by the catalyst effectiveness factor ηo [41]:
kapp = ƞo kin
Substituting this expression into the main Equation (6) and integrating it yields the final expression for the catalytic oxidation reaction of sulfur with nth-order kinetics.
C R S = [ C R S O 1 n + n 1 . t . k a p p ] ( 1 1 n )

4.3. Efficiency of Reactors

Equation (6) is rewritten as follows (using ƞ o , kin instead of k a p p ) [41] to account for more physical factors:
C R S = [ C R S O 1 n + n 1 . t . k i n . ƞ o ] ( 1 1 n )
The Arrhenius equation may be used to calculate the reaction rate constant (Kin) for each ODS reaction:
k i n = K 0   e E A R T
K0 represents the frequency or pre-exponential factor, while (E) denotes the activation energy of the reaction. This expression fits experiments conducted over extensive temperature ranges and is widely advocated from various perspectives as a highly accurate approximation of the actual temperature dependence. The chemical reaction rate can thus be formulated as follows:
C R S = [ C R S O 1 n + n 1 . t . K 0 e E A R T . ƞ o ] ( 1 1 n )
Many factors, including diffusivities, oil viscosity, and efficacy factors, are involved in the oxidation process of sulfur compounds. These parameters are obtained using the following formulae.

4.3.1. Effectiveness Ratio ( ƞ o )

The effectiveness factor is dependent on the Thiele modulus (Φ) and is determined by the following correlation, which is applicable for spherical particles [41,42]:
ƞ o = 3 ( Φ coth Φ 1 ) Φ 2
The Thiele modulus for an nth-order irreversible reaction is calculated using the following generalized equation [43,44]:
Φ = V p S p ( n + 1 ) 2 K i n C R S ( n 1 ) ρ p D e i

4.3.2. Catalyst Surface (Sp) and External Volume (Vp)

For a spherical particle, the external surface area and volume of the catalyst can be determined as follows:
V p = 4 3 π · ( r p ) 3
Sp = 4π (rp)2

4.3.3. Effective Diffusivity (Dei)

Effective diffusivity characterizes the catalyst structure (porosity and tortuosity) by considering the pore network within the particle, as described below [41,43]:
D e i = ϵ B τ 1 1 D m i + 1 D k i
The catalyst porosity (ϵB) is determined using the following two equations based on experimental data:
ϵ B = ρ p V g
ρ p = ρ B 1 ϵ B
The values of the tortuosity factor (τ) of the pore network are from 2 to 7 [44]. τ values in studies in the literature have been reported to be 4 [41,44].
Within the catalyst particle, Dei comprises two diffusivity types: molecular diffusivity (Dmi) and Knudsen diffusivity (Dki).
Molecular diffusivity is computed using the Tyn–Calus equation [45,46]:
D m i = 8.93 10 8 υ L 0.267 T υ R S 0.433 μ L
Knudsen diffusivity is determined according to the following procedure [42,47]:
D k i = 9700 r g ( T M i ) 0.5

4.3.4. Viscosity

The viscosity of light naphtha can be calculated using the equation below [48]:
μ L = 3.141 10 10 T 460 3.444 [ l o g 10 A P I ] a
a = 10.313 l o g 10 T 460 36.447
This equation is used to estimate the API.
API = 141.5 s p . g r 15.6 131.5

4.3.5. Molar Volume

The molar volume of the model sulfur compound is estimated according to the following correlation [49]:
υ R S = 0.285 ( υ c R S ) 1.048
The critical volume of liquid (light naphtha) is obtained using a Riazi–Daubert equation [49]:
υ L = 0.285   ( υ c L ) 1.048
υ c L = ( 7.5214 10 3 ( T m e A B P ) 0.2896 ( ρ L , 15.6 ) 0.7666 ) M L
The mathematical model equations for the ODS reaction outlined in this section were implemented and solved concurrently using the gPROMS (general Process Modeling System—Version 5.6) package.

5. Technique Used to Estimate Kinetic Parameters

Model predictions can be guaranteed only if the unknown parameters are estimated correctly. The design and control strategies of the process can be compromised if an accurate process model is not used.
In this work, two methodologies were applied to find the optimum kinetic parameters that depend on the amount of sulfur in the ODS process under diverse conditions of operation. The methods are given below according to references [49,50,51]:
In general, the reaction order (n) and the reaction rate constant (k) are first determined using linear regression. Subsequently, the activation energy (E) and the pre-exponential factor (Ko) are determined by linear regression using the Arrhenius equation.
The simultaneous evaluation of the reaction order (n), activation energy (E), and pre-exponential factor (Ko) was performed using non-linear regression.
The following function was minimized in order to find the best value for the kinetic parameter:
O B J = n = 1 N t ( C R S e x p C R S p r e d ) 2  
In Equation (33), Nt is the number of runs, C R S e x p is the concentration of sulfur obtained experimentally, and C R S p r e d is the value predicted by the model.

Optimization Problem Formulation for Parameter Estimation

The optimization problem formulation to estimate the parameters is as follows:
GivenThe reactor configuration, the type of catalyst used, and the process conditions of a specific process
DetermineThe order of the oxidation reaction (n1) and the reaction rate constant (k) at different temperatures (303, 333, 363, 393 K) for each catalyst and calculate the activation energy (E) and pre-exponential factor using linear regression of the Arrhenius equation.
So as to minimize The sum of squared error (SSE).
Subjected toAll optimization variables are subject to conversion rate limits and linear limitations.
The optimization problem may be expressed mathematically using linear regression as follows:
MinSSE (OBJ)
n j , k i j , ( i = 1 4 , j = CAT-1, CAT-2, CAT-3 and CAT-4),S.t. f(z, x(z),ẋ(z), u(z), v) = 0 (batch reactor)
The problem is formulated as follows when employing the non-linear regression technique. Here, reaction order (n), activation energy (E), and the pre-exponential factor (ko) are estimated simultaneously for each catalyst.
MinSSE
n j , E A j , k i j , j = CAT-1, CAT-2, CAT-3 and CAT-4),S.t.f (z, x(z),ẋ(z), u(z), v) = 0 (batch reactor)
The equation f(z, x(z), x’(z), u(z), v) = 0 represents a model where u(z) is the decision variable, z is the independent variable, x(z) are the variables, and x’(z) are their derivatives. Variable v is the design variable, and C, CL, CU, L, U are related to concentration and bounds.
The optimization technique used is a two-step method described in [52]. Step 1 simulates the process for a decision variable and calculates the objective and constraints. Step 2 updates the optimization variables.

6. Conclusions

The synthesis of new efficient catalysts through oxidative desulfurization will change oil pollutants into clean fuel. In view of the environmental concerns of high sulfur content and the difficulties in designing potential catalysts to develop better-quality fuel, achieving this synthesis is one of the most essential tasks. Oxidation reactions are very promising when catalyzed by nanoparticles. Silica can possibly be safely, economically, and environmentally extracted from bentonite. Therefore, when compared to findings obtained in the literature, this study’s preparation of nano-silica from the cheapest raw material (bentonite) produced an excellent nano material with a high surface area. Here, two composite catalysts (CAT-1 (0% HY-zeolite and 100% nano-silica) and CAT-2 (20% HY-zeolite and 80% nano-silica)) of Fe/silica were experimentally designed. The results revealed that the highest total sulfur removal efficiency was 50% and 87.88% for CAT-1 and 80% CAT-2, respectively. CAT-2 showed remarkable enhancement in the performance of the ODS process due to the positive impact of HY-zeolite on sulfur removal efficiency and the improved activity and selectivity of the catalyst in sulfur compound removal. Compared to previous studies, the obtained conversion provides a promising technology for efficient ODS processes.
The novel nanocatalyst designated as CAT-2 exhibited a markedly enhanced capacity to convert sulfur-containing compounds when air was utilized as the oxidizing agent under the predetermined conditions of a temperature of 393 Kelvin and a duration of 120 min. Both linear and non-linear optimization methods were employed in order to obtain the optimal fit of kinetic parameters, with the non-linear approach proving to be more accurate. The final mean absolute error of the model remained below 5% across all conditions tested. The efficiencies achieved with this catalyst reached over 99% of the process at a temperature of 550 Kelvin and a batch time of 193.63318 min. Ultimately, the new nanocatalyst enabled a reduction in fuel sulfur content with high quality while also lowering environmental pollution.

Author Contributions

Conceptualization, A.T.J. and A.K.H.; Methodology, A.T.J.; Software, S.A.H.; Validation, A.K.H., B.A.A.-T. and J.I.H.; Formal Analysis, I.M.M.; Investigation, L.I.S.; Resources, J.I.H.; Data Curation, B.A.A.-T.; Writing—Original Draft Preparation, A.T.J.; Writing—Review and Editing, L.I.S. and A.K.H.; Proofreading, L.I.S.; Visualization, S.A.H.; Supervision, I.M.M.; Project Administration, L.I.S.; Funding Acquisition, I.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding, and the APC was not funded by any external source.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest to disclose related to this article.

Nomenclature

D m i Molecular diffusivity (cm2/s)
S p . g r 15.6 Specific gravity of kerosene at 15.6 °C
M w L Liquid molecular weight of kerosene (g/gmol)
M w i Molecular weight of sulfur (g/gmol)
R Gas constant
r R S Reaction rate of sulfur
r g Pore radius (nm)
r p Particle radius (nm)
S p External surface area of catalyst particle (cm2)
S g Specific surface area of particle (cm2/g)
V p External volume of catalyst particle (cm3)
V g Pore volume (cm3)
T m e A B P Mean average boiling point (°C)
Greek Symbols
η 0 Effectiveness factor
Φ Thiel modulus
E B Porosity
T Tortuosity
ρ B Bulk density
ρ p Particle density
ρ L 15.6 Density of kerosene at 15.6 °C ( g m / c m 3 ) .
μ l Viscosity of liquid
v l Liquid molar volume
v c l Critical molar volume of liquid
v R S Critical volume of sulfur
o Initial (at time = 0)

References

  1. Fichman, B.T. Annual Energy Review 2011; USDOE Energy Information Administration (EIA): Washington, DC, USA, 2012. [Google Scholar]
  2. Speight, J.G. The Desulfurization of Heavy Oils and Residua; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  3. Yengejeh, S.M.; Allahyari, S.; Rahemi, N. Efficient oxidative desulfurization of model fuel by visible-light-driven MoS2-CeO2/SiO2-Al2O3 nano photocatalyst coating. Process Saf. Environ. Prot. 2020, 143, 25–35. [Google Scholar] [CrossRef]
  4. Humadi, J.I.; Nawaf, A.T.; Jarullah, A.T.; Ahmed, M.A.; Hameed, S.A.; Mujtaba, I.M. Design of new nano-catalysts and digital basket reactor for oxidative desulfurization of fuel: Experiments and modelling. Chem. Eng. Res. Des. 2023, 190, 634–650. [Google Scholar] [CrossRef]
  5. Wan, M.-W. Development of a Portable, Modular Unit for the Optimization of Ultrasound-Assisted Oxidative Desulfurization of Diesel; University of Southern California: Los Angeles, CA, USA, 2006. [Google Scholar]
  6. Shi, M.; Zhang, D.; Yu, X.; Li, Y.; Wang, X.; Yang, W. Deep oxidative desulfurization catalyzed by (NH4) 5H6PV8Mo4O40 using molecular oxygen as an oxidant. Fuel Process. Technol. 2017, 160, 136–142. [Google Scholar] [CrossRef]
  7. Hassan, S.U.; Ahmad, S.; Farid, M.A.; Nadeem, S.; Ali, Z.; Abro, R.; Mohyuddin, A.; Nazir, M.S.; Hussain, M.; Park, Y.K. Designing polyoxometalate based hybrid catalysts for efficient removal of hazardous sulfur from fuel via heterogeneous oxidative desulfurization. J. Ind. Eng. Chem. 2022, 116, 438–446. [Google Scholar] [CrossRef]
  8. Kargar, H.; Ghahramaninezhad, M.; Shahrak, M.N.; Balula, S.S. An effective magnetic catalyst for oxidative desulfurization of model and real fuels: Fe3O4/ZIF-8/TiO2. Microporous Mesoporous Mater. 2021, 317, 10992. [Google Scholar] [CrossRef]
  9. Zhang, T.; Huang, B.; Elzatahry, A.A.; Alghamdi, A.; Yue, Q.; Deng, Y. Synthesis of podlike magnetic mesoporous silica nanochains for use as enzyme support and nanostirrer in biocatalysis. ACS Appl. Mater. Interfaces 2020, 12, 17901–17908. [Google Scholar] [CrossRef] [PubMed]
  10. Saleh, T.A.; Sulaiman, K.O.; Al-Hammadi, S.A.; Dafalla, H.; Danmaliki, G.I. Adsorptive desulfurization of thiophene, benzothiophene and dibenzothiophene over activated carbon manganese oxide nanocomposite: With column system evaluation. J. Clean. Prod. 2017, 154, 401–412. [Google Scholar] [CrossRef]
  11. Yan, F.; Jiang, J.; Tian, S.; Liu, Z.; Shi, J.; Li, K.; Chen, X.; Xu, Y. A green and facile synthesis of ordered mesoporous nanosilica using coal fly ash. ACS Sustain. Chem. Eng. 2016, 4, 4654–4661. [Google Scholar] [CrossRef]
  12. Wan Abu Bakar, W.; Rusmidah, A.; Abdul Kadir, A.; Wan Mokhtar, W. The Role of Molybdenum Oxide Based Catalysts on Oxidative Desulfurization of Diesel Fuel. Mod. Chem. Appl. 2015, 3, 1–3. [Google Scholar] [CrossRef]
  13. Kamel, L.; Anbia, M. Composite beads of molybdenum oxide supported on textured silicon as an oxidative desulfurization nanocatalyst. Silicon 2020, 12, 1325–1336. [Google Scholar] [CrossRef]
  14. Zhu, Z.; Ma, H.; Liao, W.; Tang, P.; Yang, K.; Su, T.; Ren, W.; Lü, H. Insight into tri-coordinated aluminum dependent catalytic properties of dealuminated Y zeolites in oxidative desulfurization. Appl. Catal. B Environ. 2021, 288, 120022. [Google Scholar] [CrossRef]
  15. Aguado, J.; Serrano, D.; Rodríguez, J. Zeolite Beta with hierarchical porosity prepared from organofunctionalized seeds. Microporous Mesoporous Mater. 2008, 115, 504–513. [Google Scholar] [CrossRef]
  16. Zhang, S.; Lv, T.; Feng, Z.; Wang, F.; Liu, X.; Wang, W.; Wang, Y.; Meng, C. OSDA-free synthesis of zeolite beta by magadiite hydrothermal conversion method and an insight into the changes of medium-range structure during crystallization. Microporous Mesoporous Mater. 2019, 278, 81–90. [Google Scholar] [CrossRef]
  17. Hamidzadeh, M.; Komeili, S.; Saeidi, M. Seed-induced synthesis of ZSM-5 aggregates using the Silicate-1 as a seed: Characterization and effect of the Silicate-1 composition. Microporous Mesoporous Mater. 2018, 268, 153–161. [Google Scholar] [CrossRef]
  18. Dionigi, C.; Ivanovska, T.; Liscio, F.; Milita, S.; Corticelli, F.; Ruani, G. Fabrication and properties of non-isolating γ-alumina meso-foam. J. Alloys Compd. 2016, 666, 101–107. [Google Scholar] [CrossRef]
  19. Huang, P.; Luo, G.; Kang, L.; Zhu, M.; Dai, B. Preparation, characterization and catalytic performance of HPW/aEVM catalyst on oxidative desulfurization. RSC Adv. 2017, 7, 4681–4687. [Google Scholar] [CrossRef]
  20. Paraskos, J.; Frayer, J.; Shah, Y. Effect of holdup incomplete catalyst wetting and backmixing during hydroprocessing in trickle bed reactors. Ind. Eng. Chem. Process Des. Dev. 1975, 14, 315–322. [Google Scholar] [CrossRef]
  21. Alheety, M.A.; Al-Jibori, S.A.; Karadağ, A.; Akbaş, H.; Ahmed, M.H. A novel synthesis of MnO2, nanoflowers as an efficient heterogeneous catalyst for oxidative desulfurization of thiophenes. Nano-Struct. Nano-Objects 2019, 20, 100392. [Google Scholar] [CrossRef]
  22. Salman, M.S.; Jafar, S.A.; Abdullah, G.H.; Humadi, J.I.; Ahmed, M.A.; Mohammed, A.M. Development of an electrochemical reactor with rotating anode for fast and ultra-deep catalytic desulfurization of diesel: Experimental and modeling. Chem. Eng. Commun. 2024, 211, 1508–1523. [Google Scholar] [CrossRef]
  23. Zulfiqar, U.; Subhani, T.; Husain, S. Synthesis and characterization of silica nanoparticles from clay. J. Asian Ceram. Soc. 2016, 4, 91–96. [Google Scholar] [CrossRef]
  24. Al-Ghouti, M.A.; Al-Degs, Y.S.; Khalili, F.I. Minimisation of organosulphur compounds by activated carbon from Diesel fuel: Mechanistic study. Chem. Eng. J. 2010, 162, 669–676. [Google Scholar] [CrossRef]
  25. Wei, S.; He, H.; Cheng, Y.; Yang, C.; Zeng, G.; Qiu, L. Performances, kinetics and mechanisms of catalytic oxidative desulfurization from oils. RSC Adv. 2016, 6, 103253–103269. [Google Scholar] [CrossRef]
  26. Poyton, A.; Varziri, M.S.; McAuley, K.B.; McLellan, P.J.; Ramsay, J.O. Parameter estimation in continuous-time dynamic models using principal differential analysis. Comput. Chem. Eng. 2006, 30, 698–708. [Google Scholar] [CrossRef]
  27. Nawaf, A.T.; Humadi, J.I.; Jarullah, A.T.; Ahmed, M.A.; Hameed, S.A.; Mujtaba, I.M. Design of nano-catalyst for removal of phenolic compounds from wastewater by oxidation using modified digital basket baffle batch reactor: Experiments and modeling. Processes 2023, 11, 1990. [Google Scholar] [CrossRef]
  28. Maslova, O.; Senko, O.; Akopyan, A.; Lysenko, S.; Anisimov, A.; Efremenko, E. Nanocatalysts for oxidative desulfurization of liquid fuel: Modern solutions and the perspectives of application in hybrid chemical-biocatalytic processes. Catalysts 2021, 11, 1131. [Google Scholar] [CrossRef]
  29. Santos, M.R.F.D.; Pedrosa, A.M.G.; Souza, M.J.B.D. Oxidative desulfurization of thiophene on TiO2/ZSM-12 zeolite. Mater. Res. 2016, 19, 24–30. [Google Scholar] [CrossRef]
  30. Ugal, J.R.; Jima’a, R.B.; Al-Jubori, W.M.K.; Bayader, F.A.; Al-Jubori, N.M. Oxidative desulfurization of hydrotreated gas oil using Fe2O3 and Pd loaded over activated carbon as catalysts. Orient. J. Chem. 2018, 34, 1091. [Google Scholar] [CrossRef]
  31. Dehkordi, A.M.; Kiaei, Z.; Sobati, M.A. Oxidative desulfurization of simulated light fuel oil and untreated kerosene. Fuel Process. Technol. 2009, 90, 435–445. [Google Scholar] [CrossRef]
  32. Perego, C.; Millini, R. Porous materials in catalysis: Challenges for mesoporous materials. Chem. Soc. Rev. 2013, 42, 3956–3976. [Google Scholar] [CrossRef]
  33. Goldberg, M.A.; Akopyan, A.V.; Gafurov, M.R.; Makshakova, O.N.; Donskaya, N.O.; Fomin, A.S.; Polikarpova, P.P.; Anisimov, A.V.; Murzakhanov, F.F.; Leonov, A.V.; et al. Iron-doped mesoporous powders of hydroxyapatite as molybdenum-impregnated catalysts for deep oxidative desulfurization of model fuel: Synthesis and experimental and theoretical studies. J. Phys. Chem. C 2021, 125, 11604–11619. [Google Scholar] [CrossRef]
  34. Corma, A. From microporous to mesoporous molecular sieve materials and their use in catalysis. Chem. Rev. 1997, 97, 2373–2420. [Google Scholar] [CrossRef]
  35. Karre, A.V.; Dadyburjor, D.B. Review of iron-based catalysts with and without zeolite supports used in fischer-tropsch processes. Chem. Eng. Commun. 2022, 209, 967–987. [Google Scholar] [CrossRef]
  36. Zhang, J.; Tang, X.; Yi, H.; Yu, Q.; Zhang, Y.; Wei, J.; Yuan, Y. Synthesis, characterization and application of Fe-zeolite: A review. Appl. Catal. A Gen. 2022, 630, 118467. [Google Scholar] [CrossRef]
  37. Abd Zaid, S.M.; AbdulRazak, A.A.; Abid, M.F. A review of nano-catalyst applications in kerosene desulfurization techniques. J. Appl. Sci. Nanotechnol. 2022, 2, 86–102. [Google Scholar] [CrossRef]
  38. Srivastava, V.C. An evaluation of desulfurization technologies for sulfur removal from liquid fuels. RSC Adv. 2012, 2, 759–783. [Google Scholar] [CrossRef]
  39. Khan, N.A.; Hasan, Z.; Jhung, S.H. Adsorptive removal of hazardous materials using metal–organic frameworks (MOFs): A review. J. Hazard. Mater. 2013, 244–245, 444–456. [Google Scholar] [CrossRef] [PubMed]
  40. Nadkarni, R.K. Spectroscopic Methods for the Determination of Sulfur in Petroleum Products; ASTM International: West Conshohocken, PA, USA, 2011. [Google Scholar]
  41. Tam, P.S.; Kittrell, J.R.; Eldridge, J.W. Desulfurization of fuel oil by oxidation and extraction. 1. Enhancement of extraction oil yield. Ind. Eng. Chem. Res. 1990, 29, 321–324. [Google Scholar] [CrossRef]
  42. Nehlsen, J.; Benziger, J.; Kevrekidis, I. Oxidation of aliphatic and aromatic sulfides using sulfuric acid. Ind. Eng. Chem. Res. 2006, 45, 518–524. [Google Scholar] [CrossRef]
  43. Bagherzadeh, M.; Amini, M. Synthesis, characterization and catalytic study of a novel iron (III)-tridentate Schiff base complex in sulfide oxidation by UHP. Inorg. Chem. Commun. 2009, 12, 21–25. [Google Scholar] [CrossRef]
  44. Jiang, X.; Li, H.; Zhu, W.; He, L.; Shu, H.; Lu, J. Deep desulfurization of fuels catalyzed by surfactant-type decatungstates using H2O2 as oxidant. Fuel 2009, 88, 431–436. [Google Scholar] [CrossRef]
  45. Di Giuseppe, A.; Crucianelli, M.; De Angelis, F.; Crestini, C.; Saladino, R. Efficient oxidation of thiophene derivatives with homogeneous and heterogeneous MTO/H2O2 systems: A novel approach for, oxidative desulfurization (ODS) of diesel fuel. Appl. Catal. B Environ. 2009, 89, 239–245. [Google Scholar] [CrossRef]
  46. Wallace, T.; Schriesheim, A.; Hurwitz, H.; Glaser, M. Base-catalyzed oxidation of mercaptans in presence of inorganic transition metal complexes. Ind. Eng. Chem. Process Des. Dev. 1964, 3, 237–241. [Google Scholar] [CrossRef]
  47. Zannikos, F.; Lois, E.; Stournas, S. Desulfurization of petroleum fractions by oxidation and solvent extraction. Fuel Process. Technol. 1995, 42, 35–45. [Google Scholar] [CrossRef]
  48. Sengupta, A.; Kamble, P.D.; Basu, J.K.; Sengupta, S. Kinetic study and optimization of oxidative desulfurization of benzothiophene using mesoporous titanium silicate-1 catalyst. Ind. Eng. Chem. Res. 2012, 51, 147–157. [Google Scholar] [CrossRef]
  49. Al-Malki, A. Desulfurization of Gasoline and Diesel Fuels, Using Non-Hydrogen Consuming Techniques; King Fahd University of Petroleum and Minerals: Dhahran, Saudi Arabia, 2004. [Google Scholar]
  50. Nawaf, A.T.; Gheni, S.A.; Jarullah, A.T.; Mujtaba, I.M. Optimal design of a trickle bed reactor for light fuel oxidative desulfurization based on experiments and modeling. Energy Fuels 2015, 29, 3366–3376. [Google Scholar] [CrossRef]
  51. Ion, C.S.; Bombos, M.; Doukeh, R.; Vasilievici, G.; Matei, V. Kinetics of 1-dodecanethiol Desulfurization by Reactive Adsorbtion on MgO/dolomite. Rev. Chim. 2018, 69, 3439–3444. [Google Scholar] [CrossRef]
  52. Jarullah, A.T.; Mujtaba, I.M.; Wood, A.S. Kinetic model development and simulation of simultaneous hydrodenitrogenation and hydrodemetallization of crude oil in trickle bed reactor. Fuel 2011, 90, 2165–2181. [Google Scholar] [CrossRef]
Figure 1. XRD for amorphous nano-silica.
Figure 1. XRD for amorphous nano-silica.
Catalysts 14 00529 g001
Figure 2. Analysis of nano-silica using FTIR.
Figure 2. Analysis of nano-silica using FTIR.
Catalysts 14 00529 g002aCatalysts 14 00529 g002b
Figure 3. FTIR spectrum for CAT-2.
Figure 3. FTIR spectrum for CAT-2.
Catalysts 14 00529 g003
Figure 4. TGA analysis of (A) nano-silica and (B) composite catalyst.
Figure 4. TGA analysis of (A) nano-silica and (B) composite catalyst.
Catalysts 14 00529 g004
Figure 5. The size distributions of nano-silica and the produced catalysts.
Figure 5. The size distributions of nano-silica and the produced catalysts.
Catalysts 14 00529 g005
Figure 6. SEM image of nano-silica.
Figure 6. SEM image of nano-silica.
Catalysts 14 00529 g006
Figure 7. EDX image of nano-silica.
Figure 7. EDX image of nano-silica.
Catalysts 14 00529 g007
Figure 8. FESEM of the catalyst.
Figure 8. FESEM of the catalyst.
Catalysts 14 00529 g008
Figure 9. EDX mapping of the catalyst.
Figure 9. EDX mapping of the catalyst.
Catalysts 14 00529 g009
Figure 10. Effect of reaction temperature on sulfur conversion at several reaction times. Effect of reaction temperature on sulfur conversion for different catalyst types at various times: (A) 30 min, (B) 60 min, (C) 90 min, (D) 120 min.
Figure 10. Effect of reaction temperature on sulfur conversion at several reaction times. Effect of reaction temperature on sulfur conversion for different catalyst types at various times: (A) 30 min, (B) 60 min, (C) 90 min, (D) 120 min.
Catalysts 14 00529 g010
Figure 11. Effect of reaction time for different catalyst types on the sulfur conversion at various reaction temperatures. Effect of reaction time on sulfur conversion at different temperatures: (A) 303 K, (B) 333 K, (C) 363 K, (D) 393 K.
Figure 11. Effect of reaction time for different catalyst types on the sulfur conversion at various reaction temperatures. Effect of reaction time on sulfur conversion at different temperatures: (A) 303 K, (B) 333 K, (C) 363 K, (D) 393 K.
Catalysts 14 00529 g011
Figure 12. CAT-1 oxidation kinetics: ln(k) vs. 1/T.
Figure 12. CAT-1 oxidation kinetics: ln(k) vs. 1/T.
Catalysts 14 00529 g012
Figure 13. Kinetics of oxidation with CAT-2, displaying lnk against 1/T.
Figure 13. Kinetics of oxidation with CAT-2, displaying lnk against 1/T.
Catalysts 14 00529 g013
Figure 14. Comparison of CAT-1 data from experiments and simulations.
Figure 14. Comparison of CAT-1 data from experiments and simulations.
Catalysts 14 00529 g014
Figure 15. Comparison of CAT-2 data from experiments and simulations.
Figure 15. Comparison of CAT-2 data from experiments and simulations.
Catalysts 14 00529 g015
Figure 16. Steps of nano-silica preparation.
Figure 16. Steps of nano-silica preparation.
Catalysts 14 00529 g016
Figure 17. Catalyst preparation using the IWI method.
Figure 17. Catalyst preparation using the IWI method.
Catalysts 14 00529 g017
Figure 18. Schematic layout of a batch ODS reactor.
Figure 18. Schematic layout of a batch ODS reactor.
Catalysts 14 00529 g018
Table 1. Results of the produced nano-silica’s XRF analysis.
Table 1. Results of the produced nano-silica’s XRF analysis.
ElementSymbolConcentration %
SiliconSiO271.14
MagnesiumMgO0.0552
AluminumAl2O30.4179
PhosphorusP2O50.5621
SulfurSO35.729
ChlorineCl0.03876
PotassiumK2O<0.0012
CalciumCaO0.1619
TitaniumTiO23.010
ChromiumCr2O30.01824
ManganeseMnO0.00367
IronFe2O30.6977
CobaltCoO<0.00039
NickelNiO0.00648
CopperCuO0.00075
ZincZnO0.00274
GalliumGa0.00045
Another element-0.11952
Sum of concentration-81.97
Table 2. The prepared catalysts’ XRF results.
Table 2. The prepared catalysts’ XRF results.
ElementSymbolConcentration %
SiliconSiO251.56
MagnesiumMgO0.0863
AluminumAl2O311.24
PhosphorusP2O50.6316
SulfurSO33.652
PotassiumK2O0.0927
CalciumCaO0.1326
TitaniumTiO22.995
VanadiumV2O50.0314
ChromiumCr2O30.00472
ManganeseMnO0.01570
IronFe2O310.88
CopperCuO0.00603
ZincZnO0.00510
Another component 0.091201
Sum of concentration 81.45
Table 3. Parameters of the ODS model and their values.
Table 3. Parameters of the ODS model and their values.
ParameterNotationMeasurement UnitNumerical Value
Initial amountCtwt%0.1966
Durationtime1, time2, time3, time4mintime1 = 30, time2 = 60, time3 = 90, time4 = 120
TemperatureT1, T2, T3, T4KT1 = 303, T2 = 333, T3 = 363, T4 = 393
Kerosene density at 15.5 °C ρ l gm/cm30.8205
Average boiling pointTmeABP°R787
Acceleration gravitygm/sec29.81
Gas constantRJ/mole.°K8.314
Pore volume per unit catalyst massVgcm3/gmVgCAT-1 = 0.606
VgCAT-2 = 0.625
Particle surface areaSgcm2/gmSgCAT-1 = 2,426,700
SgCAT-2 = 2,501,380
Catalyst particle volumeVpcm3VpCAT-1 = 5.2158 × 10−19
VpCAT-2 = 5.2245 × 10−19
Particle’s external surface areaSpcm2SpCAT-1 = 3.13297 × 10−12
SpCAT-2 = 3.1365 × 10−12
DensityρBgm/cm3ρB.CAT-1 = 0.667
ρB.CAT-2 = 0.687
Molecular weight of keroseneMwLgm/mole170
Molecular weight of sulfurMwigm/mole32.06
Pore radius averagergnmrg,CAT-1 = 4.9944
rg,CAT-2 = 4.9972
Table 4. Parameters for CAT-1 obtained using linear regression.
Table 4. Parameters for CAT-1 obtained using linear regression.
ParameterValueUnit
n2.000_
k10.0106 ( w t 1 ) · m i n 1
k20.0233 ( w t 1 ) · m i n 1
k30.0544 ( w t 1 ) · m i n 1
k40.1296 ( w t 1 ) · m i n 1
SSE2.9922 × 10−4_
Table 5. CAT-2 model parameters obtained using linear regression.
Table 5. CAT-2 model parameters obtained using linear regression.
ParameterValueUnit
n1.9861_
k10.0207 ( w t 0.9861 ) · m i n 1
k20.0538 ( w t 0.9861 ) · m i n 1
k30.1165 ( wt 0.9861 ) · min 1
k40.2278 ( wt 0.9861 ) · min 1
SSE2.1999 × 10−4_
Table 6. Activation energy calculations using the linear regression technique.
Table 6. Activation energy calculations using the linear regression technique.
CatalystE (kJ/mol)Frequency Factor
CAT-125.8108399.2151
CAT-226.2407710.1669
Table 7. Parameters for the kinetic model of CAT-1 obtained using non-linear regression.
Table 7. Parameters for the kinetic model of CAT-1 obtained using non-linear regression.
ParameterValueUnit
n1.9955_
E30.7646KJ/mol.
ko1455.7723 w t 0.9955 · m i n 1
SSE1.4621 × 10−4_
Table 8. Parameters for the kinetic model of CAT-2 obtained using non-linear regression.
Table 8. Parameters for the kinetic model of CAT-2 obtained using non-linear regression.
ParameterValueUnit
n1.9869_
E25.8635KJ/mol.
ko1681.5340 w t 0.9869 · m i n 1
SSE1.5799 × 10−4_
Table 9. Experimental and simulation results for CAT-1 obtained using linear regression.
Table 9. Experimental and simulation results for CAT-1 obtained using linear regression.
Temperature
(K)
Time
(min)
Sulfur Content (wt. %)ConversionError
%
ExperimentalPredictedExperimentalPredicted
303300.19020.18510.03250.05902.6940
303600.18150.17480.07680.11083.6760
303900.17090.16560.13070.15763.0688
3031200.16410.15740.16530.19984.0846
333300.17490.17290.08740.12101.1594
333600.15780.15420.19730.21562.2465
333900.14570.13920.25890.29194.4216
3331200.13340.12690.32140.35454.8582
363300.15340.14880.21970.24312.9650
363600.12270.11970.37580.39112.3924
363900.09910.10020.49590.49081.0974
3631200.08940.08610.54520.56253.6793
393300.10990.11140.44090.43331.3928
393600.07990.07770.59530.60472.6920
393900.06220.05970.68360.69644.0148
3931200.04960.04840.74770.75362.3081
Table 10. Experimental and simulation results for CAT-2 obtained using linear regression.
Table 10. Experimental and simulation results for CAT-2 obtained using linear regression.
Temperature
(K)
Time
(min)
Sulfur Content (wt. %)ConversionError %
ExperimentalPredictedExperimentalPredicted
303300.16940.17470.13830.11133.1557
303600.14990.15720.23750.20044.8965
303900.13980.14290.28890.27312.2202
3031200.12950.13090.34130.33411.1193
333300.14160.14830.27970.24564.7459
333600.11520.11890.41400.39523.2956
333900.09880.09930.49750.49490.5145
3331200.08760.08520.55440.56672.7597
363300.10990.11530.44090.41404.8970
363600.07840.08140.60120.58593.8378
363900.06240.06290.68260.68050.7395
3631200.04970.05120.74720.74002.9511
393300.07950.08250.59560.58043.7622
393600.04960.05200.76140.73544.8996
393900.03990.03790.79700.80704.9001
3931200.02880.02980.85350.84823.5993
Table 11. Second-approach experimental and simulation results for CAT-1.
Table 11. Second-approach experimental and simulation results for CAT-1.
Temperature
(K)
Time
(min)
Sulfur Content (wt. %)ConversionError
%
ExperimentalPredictedExperimentalPredicted
303300.19020.18850.03250.04170.8947
303600.18150.18100.07680.07930.2555
303900.17090.17410.13070.11441.8965
3031200.16410.16780.16530.14692.2248
333300.17490.17410.08740.11440.4521
333600.15780.15620.19730.20540.9967
333900.14570.14170.25890.27942.7661
3331200.13340.12960.32140.34132.8569
363300.15340.14850.21970.24463.1823
363600.12270.11930.37580.39312.7662
363900.09910.09970.49590.49330.5862
3631200.08940.08560.54520.56514.2606
393300.10990.11530.44090.41404.8952
393600.07990.08150.59530.58542.0034
393900.06220.06300.68360.67951.3054
3931200.04960.05140.74770.73883.5296
Table 12. Second-approach experimental and simulation results for CAT-2.
Table 12. Second-approach experimental and simulation results for CAT-2.
Temperature
(K)
Time
(min)
Sulfur Content (wt. %)ConversionError %
ExperimentalPredictedExperimentalPredicted
303300.16940.17380.13830.11582.6066
303600.14990.15570.23750.20783.8948
303900.13980.14100.28890.28250.8931
3031200.12950.12890.34130.34440.4810
333300.14160.14770.27970.24864.3185
333600.11520.11820.41400.39862.6218
333900.09880.09850.49750.49890.3038
3331200.08760.08440.55440.57073.6631
363300.10990.11450.44090.41784.1502
363600.07840.08060.60120.59002.7996
363900.06240.06210.68260.68390.4212
3631200.04970.05050.74720.74291.6689
393300.07950.08250.59560.58053.7277
393600.04960.05200.76140.73544.8722
393900.03990.03790.79700.80704.9152
3931200.02880.02980.85350.84823.6126
Table 13. Comparison of mean squared errors.
Table 13. Comparison of mean squared errors.
ParametersThe Summation of Squared Errors (SSE)
Linear ApproachNon-Linear Approach
CAT-1299.2197 × 10−6146.213 × 10−6
CAT-2219.9992 × 10−6157.995 × 10−6
Table 14. Optimal conditions for applying the ODS process with CAT-1.
Table 14. Optimal conditions for applying the ODS process with CAT-1.
ParameterValueUnit
Csulfur,t0.1911wt.%
T546K
Time199min
Conversion99.0432%
Table 15. Optimal conditions for applying the ODS process with CAT-2.
Table 15. Optimal conditions for applying the ODS process with CAT-2.
ParameterValueUnit
Csulfur,t0.2443wt.%
T550K
Time193.6332min
Conversion99.1042%
Table 16. Comparison of results.
Table 16. Comparison of results.
Sulfur CompoundOxidant/CatalystSulfur RemovalReference
THn-heptane,
H2O2/TiO2/ZSM-12
60%[29]
DBTDiesel, H2O2/Fe2O3/
AC
71%[30]
S-compoundsKerosene, H2O2/formic acid87%[31]
S-compoundsKerosene, air/20% HY-zeolite and 80% nano-silica87.88%Present study
Table 17. Variables and operating conditions for this work.
Table 17. Variables and operating conditions for this work.
VariablesConditions
CatalystCAT-1, CAT-2
Temperature (°C)30, 60, 90, 120
Reaction time (min)30, 60, 90, 120
Table 18. Specification of kerosene.
Table 18. Specification of kerosene.
PropertyValue
Flash point, °C40.5
API47.5
Density at 15.6 °C, g/cm30.7904
Color+20
IBP, °C147
EP, °C265
Sulfur, ppm1966
Table 19. Materials and supplies needed to prepare catalysts.
Table 19. Materials and supplies needed to prepare catalysts.
Materials Purity %CompanyRole
Deionized water-SDISolvent of active materials
Fe(NO3)3·9H2O98HimediaActive metal
Table 20. HY-zeolite nanoparticle specifications.
Table 20. HY-zeolite nanoparticle specifications.
CharacteristicValue
Crystallinity percentage92.3
Pore volume (cm3/g)0.388
Surface area (m2/g)499.4
Particle size (nm)94.31
Bulk density (g/cm3)0.648
Pore size (nm)31.14
Table 21. XRF results for HY-zeolite nanoparticles.
Table 21. XRF results for HY-zeolite nanoparticles.
HY-Zeolite ComponentWt. %
SiO257.96
Al2O319.01
Na2O1.954
SiO2/Al2O3 molar ratio5.183
Table 22. The main equipment used.
Table 22. The main equipment used.
EquipmentSpecificationSupplier
Tubular FurnaceHeating length 300 mm, Max. temperature 1200 °CSAFTherm, Luoyang,
China
Magnetic heater stirrerStirring speed 100–1500 rpm,
capacity 100–3000 mL, Max. temperature 380 °C
JISICO, Seoul, Republic of Korea
Vacuum pump1/4 hp, model VE125NManufacturer is Value, DEIRA DUBAI,
Dubai, UAE
X-ray diffraction instrument XRD-6000 modelShimadzu company, Tokyo, Japan
Atomic Force Microscope device SPM-AT 3000,
Scanning Probe Microscope, AA 3000
Angstrom-Advance Inc., contact model Stoughton, MA, USA
FESEM deviceZeiss-EM10C-100KVZEISS Sigma company, Oberkochen, Germany
Table 23. Nano-silica production materials.
Table 23. Nano-silica production materials.
Materials and ComponentsPurity%Supplier
Sodium hydroxide (NaOH)98%Sigma Aldrich (Burlington, MA, USA)
Ammonium hydroxide (NH4OH)28%GCC Ltd. Company (London, UK)
Sulfuric acid (H2SO4)98%GCC Ltd. company
Bentonite clay97%GCC Ltd. company
Deionized water-CHD Chemicals Ltd. Company (Chandigarh, India)
Ethanol (C2H5OH)98%Analytical laboratory (Telford, PA, USA)
Table 24. The surface area, pore volume, and bulk density of the catalysts.
Table 24. The surface area, pore volume, and bulk density of the catalysts.
Catalyst TypeSurface Area, m2/gPore Volume, cm3/gBulk Density, g/cm3
SiO2 (support)341.310.8530.932
CAT-1 (0% HY-zeolite and 100% nano-silica)242.670.6060.667
CAT-2 (20% HY-zeolite and 80% nano-silica)250.1380.6250.687
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jarullah, A.T.; Hussein, A.K.; Al-Tabbakh, B.A.; Hameed, S.A.; Mujtaba, I.M.; Saeed, L.I.; Humadi, J.I. Production of Green Fuel Using a New Synthetic Magnetite Mesoporous Nano-Silica Composite Catalyst for Oxidative Desulfurization: Experiments and Process Modeling. Catalysts 2024, 14, 529. https://doi.org/10.3390/catal14080529

AMA Style

Jarullah AT, Hussein AK, Al-Tabbakh BA, Hameed SA, Mujtaba IM, Saeed LI, Humadi JI. Production of Green Fuel Using a New Synthetic Magnetite Mesoporous Nano-Silica Composite Catalyst for Oxidative Desulfurization: Experiments and Process Modeling. Catalysts. 2024; 14(8):529. https://doi.org/10.3390/catal14080529

Chicago/Turabian Style

Jarullah, Aysar T., Ahmed K. Hussein, Ban A. Al-Tabbakh, Shymaa A. Hameed, Iqbal M. Mujtaba, Liqaa I. Saeed, and Jasim I. Humadi. 2024. "Production of Green Fuel Using a New Synthetic Magnetite Mesoporous Nano-Silica Composite Catalyst for Oxidative Desulfurization: Experiments and Process Modeling" Catalysts 14, no. 8: 529. https://doi.org/10.3390/catal14080529

APA Style

Jarullah, A. T., Hussein, A. K., Al-Tabbakh, B. A., Hameed, S. A., Mujtaba, I. M., Saeed, L. I., & Humadi, J. I. (2024). Production of Green Fuel Using a New Synthetic Magnetite Mesoporous Nano-Silica Composite Catalyst for Oxidative Desulfurization: Experiments and Process Modeling. Catalysts, 14(8), 529. https://doi.org/10.3390/catal14080529

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop