Abstract
The partial oxidation of ethane to ethylene and acetic acid on supported MoVNbPd/TiO2 (P25 of Degussa) has been investigated. Pd was added in a nano-metallic form. The catalyst composition was also different from similar studied catalysts. This results in a better selectivity towards acetic acid formation. The reaction was carried out in a tubular reactor at temperature range 225–275 °C, total pressure range 0–200 psig and oxygen percentage in the feed gas of 10–40%. The feed gas contains ethane and oxygen. In this work, we develop a kinetic model for the reaction for the developed catalyst. In this model, we assume that oxidation reactions take place on different sites; ethane oxidation takes place on one site, ethylene oxidation on another site, and CO is oxidized to CO2 on a third site. The model exhibits good agreement with the experimental data.
Similar content being viewed by others
Introduction
Acetic acid is mainly produced industrially on large scale from methanol carbonylation which accounts for 65% of world capacity, but gives rise to many problems, such as corrosion and waste disposal. Other processes include acetaldehyde oxidation and ethane partial oxidation. A plant for partial oxidation of ethane to acetic acid is operated by SABIC of Saudi Arabia since 2005. Vinegar which is about 4–6% diluted form of acetic acid is obtained by fermentation of renewable carbon sources containing ethanol. Acetic acid is primarily used in the production of vinyl acetate monomer (VAM), acetic anhydride and as a solvent in the production of terephthalic acid (TPA).
MoVNb catalysts were shown [1] to be able to effect the partial oxidation of ethane to ethylene and acetic acid. The addition of Pd makes the catalyst more selective towards acetic acid formation.
The design of industrial reactors for partial oxidation of ethane over MoVNbPd catalyst requires a reliable kinetic model. In this paper we attempt to obtain kinetic expression for ethane partial oxidation based on known mechanisms that are usually used in partial oxidation reactions.
The following three mechanisms are usually used in partial oxidation reactions;
-
1.
Langmuir–Hinshelwood (LH) in which the adsorbed oxygen and adsorbed ethane on the active sites of the catalyst react in the adsorbed states,
-
2.
Eley–Rideal (ER) in which one of the reactants is adsorbed and react in the adsorbed state with the other reactant in the gas phase.
-
3.
Mars–Van Krevelen redox (MvK) in which the oxidized catalyst react with ethane and thus is reduced and is reoxidized with molecular oxygen.
In the pioneering work of Thorsteinson et al. Mo–V–Nb catalysts [1], have been reported as active for the partial oxidation of ethane to ethylene and acetic acid using oxygen. In their work the kinetics for Mo16V8Nb2.29 was derived assuming LH mechanism. Two schemes have been suggested to explain the mechanism. In one scheme, ethoxy–hydroxyl species are formed and these decompose to ethylene. In the second scheme a hydroxy–alkoxy is formed and this decomposes to ethylene and water. Bare metal sites M reacts with water to form MOH which reacts with ethylene to give an ethoxide MOC2H5. The ethoxide is oxidized to acetic acid.
For acetic acid production, a pressure of about 300 psig gives acetic acid selectivity about 20%, and ethylene selectivity about 70%. In addition, the use of high pressure is advantageous in reducing contact time for a specified ethane conversion and thus reducing the catalytic reactor bed.
To avoid the formation of inflammable mixture, the feed to the reactor contains low percentage of oxygen limiting ethane conversion. This requires the recycle of unreacted ethane.
The kinetics of Thoresteinson et al. has been used for simulation of fixed bed reactor [2], and fluidized bed reactor [3]. To avoid the formation of inflammable mixture, El-Sherehy et al. [4] studied the use of distributed oxygen along the reactor bed.
Karim et al. [5] showed that the addition of Pd to the MoVNb catalyst greatly increases the selectivity to acetic acid to about 80% and completely oxidizes CO to CO2.
Fakeeha et al. [6] obtained the kinetics for MoVPO catalyst. A model based on Eley–Rideal and Mars and Van Krevelen redox mechanisms were found adequate.
The mechanism and kinetics of the reaction with a Mo1V0.25Nb0.12Pd0.0005O x catalyst were studied by Linke et al. [7, 8]. It has been shown that ethylene inhibits ethane oxidation through depletion of lattice oxygen (O*). The consecutive oxidation of ethylene to acetic acid is itself catalyzed by the palladium oxide in a heterogeneous analogue of the Wacker process.
Their kinetic model contains some negative activation energy which is physically wrong. They have indicated that the C–H bond activation in ethane by the oxygen-saturated catalyst surfaces is the rate limiting step [7]. Water is also believed to increase acetic acid selectivity by promoting the desorption of acetate species as acetic acid [7].
The kinetics of MoV type catalyst was re-investigated by Rahman et al. [9] but the catalyst composition was not mentioned. Most probably it contains Nb and Pd. Rahman et al. [9] developed a two-site Eley–Rideal-Redox (ERR) model to predict the partial oxidation of ethane to ethylene and the partial ethane oxidation to acetic acid over the catalyst. To be used in acetic acid plant design, we need to know how much CO is produced. However, their model assumes that the oxidation of ethylene and acetic acid to CO is insignificant. Also the oxidation of ethane to CO and CO2 is negligible.
Li and Iglesia [10, 11] found that precipitation of Mo, V and Nb salts solution in the presence of colloidal TiO2 (titania P25 from Degussa) led to a tenfold increase in ethylene and acetic acid rates (per active oxide) without significant changes in selectivity relative to unsupported samples. The introduction of trace amounts of Pd (0.0025–0.01 wt%) as a physical mixture of separate 0.3 (wt%) Pd/SiO2 led to the near complete depletion of ethylene and to a significant increase in acetic acid synthesis rate.
A suitable kinetic model for the catalytic oxidation of ethane to acetic acid on MoVNb catalyst must be able to explain the following key observations:
-
1.
Rate of depletion of ethane is first order with respect to ethane as noted by Thorsteinson et al. [1]. Burch and Swarnakar [12] determined the rate of reaction to be 0.8–1.0 with respect to ethane depending on temperature.
-
2.
Rate of depletion of ethane is a fraction close to zero with respect to oxygen as noted by Thorsteinson et al. [1]. Burch and Swarnakar [12] determined the rate of reaction to be 0.07–0.5 with respect to oxygen.
-
3.
High pressure favors the formation of acetic acid.
-
4.
Temperature increase makes the reaction more selective to ethylene with respect to acetic acid.
-
5.
The strong acceleration of the rate of the ethylene oxidation to acetic acid by the presence of water (Linke et al. [7, 8]; Rahman et al. [9]). Water, whether formed as a byproduct or added with C2H6–O2 reactants, increases acetic acid selectivity by promoting the desorption of adsorbed acetate species as acetic acid. Thus, the promoting effect of water on acetic acid selectivity reflects a specific increase in the rate at which ethane converts to acetic acid via direct pathways. The absence of water led to larger (ethylene/acetic acid) ratios in products and to slightly higher CO x selectivity.
-
6.
Palladium catalyzes the oxidation of CO to CO2. Berlowitz et al. [13] have shown that for some range of operating conditions the rate of CO oxidation is first order in oxygen and negative first order in CO.
The catalyst used in our experiments has the composition Mo16V6.37Nb2.05O x Pd 00.0037 /TiO2 (P25).The catalyst loading on Titania is 30%. Pd was added. The method of preparation, catalyst characterization and the effect of operating parameters for the partial oxidation reaction of ethane such as temperature, space time and feed composition are described in references [15, 16]. The experimental tubular reactor used to carry out the experiments is also described in these two references. This catalyst composition was optimized for maximum acetic acid yield, and thus it is required to develop a kinetic model for its use in the partial oxidation reaction of ethane to acetic acid.
In the next section, the experiments carried out and the results are presented. This is followed by the kinetic model suggested. Parameter estimation results are then presented. The results are then discussed and reaction mechanism is suggested. Finally, conclusions are presented.
Experimental
The reaction system set up is described in Ref. [16]. Briefly 0.3 g of the catalyst diluted with silica particles of the same size is placed in the reactor tube of overall length 150 mm with inside diameter of made of 6.35 mm stainless steel and surrounded by brass block. The block is surrounded by mica band heater.
At the beginning of the experiments, we tested the effect of reducing the catalyst particle size on ethane conversion. No noticeable change happened. We concluded that intraparticle mass transfer is negligibly small.
For the present kinetic study sixty-three experiments were carried out for ethane oxidation. Temperatures varied between 225 and 275 °C, pressures from atmospheric to 200 psig, total gas inlet flow-rate containing ethane and oxygen from 10 to 40 ml/min, and O2% in the feed between 10 and 40. In all experiments, catalyst weight was 0.3 g. Specifically, the following experiments were carried out:
-
1.
At pressure 200 psig, and flow-rate 10 ml/min;
(a) One experiment at 225 °C and 10% O2
(b) Two experiments at 225, 240 °C and 20% O2
(c) One experiment at 225 °C and 40% O2
-
2.
At pressure 200 psig, and flow-rate 20 ml/min;
(a) Two experiments at 225, 240 °C and 10% O2
(b) Three experiments at 225, 240, 250 °C and 20% O2
(c) Three experiments at 225, 250, 275 °C and 40% O2
-
3.
At pressure 200 psig and flow-rate 40 ml/min
(a) Two experiments at 225, 250 °C and 10% O2
(b) Four experiments at 225, 250 °C and 20 and 40% O2
-
4.
At pressure 100 psig and flow-rate 10 ml/min
(a) One experiment at 225 °C and 10% O2
(b) Two experiments at 225, 250 °C and 20% O2
(c) Three experiments at 225, 250, 275 °C and 40% O2
-
5.
At pressure 100 psig and flow-rate 20 ml/min,
(a) One experiment at 225 °C and 10% O2
(b) Two experiments at 225, 250 °C and 40% O2.
-
6.
At pressure 100 psig and flow-rate 40 ml/min, nine experiments at 225, 250, 275 °C and 10, 20 and 40% O2.
-
7.
At pressure 0 psig, twenty-seven experiments at 225, 250, 275 °C and 10, 20 and 40% O2, flow-rates of 10, 20 and 40 ml/min.
The results are shown in Table 1.
Some observations can be made from testing results:
-
1.
Acetic acid is not formed at atmospheric pressure.
-
2.
Ethylene is the primary product at low temperature and low pressure.
-
3.
Ethylene selectivity decreases with the increase of temperature.
-
4.
Acetic acid is the primary product at high pressure and low temperature.
-
5.
In all experiments the selectivity of CO was very low.
-
6.
The selectivity of CO2 increases with the increase of temperature and the increase in O2% in the feed.
-
7.
The pressure range for the experiments is very wide with high selectivity to ethylene and CO2 at low pressure and high selectivity to acetic acid at high pressure. It would be a challenge to find a model that is suitable for wide variation in the pressure.
Kinetic modeling of partial oxidation of ethane
Here we attempt to obtain kinetic expression based on known mechanisms that are usually used in partial oxidation reactions. These are Langmuir–Hinshelwood (LH) kinetics assuming reaction taking place between adsorbed species, Eley–Rideal kinetics assuming reaction taking place between adsorbed species and gaseous reactant, and Mars–Van Krevelen redox mechanism in which the oxidized catalyst react with ethane and thus is reduced and is reoxidized with molecular oxygen.
For the reaction
where A could be a hydrocarbon and B is the corresponding dehydrogenated hydrocarbon.
Langmuir–Hinshelwood kinetics takes the form
While the redox kinetics take the form
The reaction network is assumed in its most general form of Fig. 1 (ri’s denote reaction rates of different species).
This assumes that acetic acid can be obtained from ethane and ethylene and that all reactants and products are oxidized to CO and CO2 and CO is converted to CO2.
We assume the conversion of ethane is X1, the conversion to CO is X2, CO2 is X3, and CH3COOH is X4. In addition, we assume that the reactor operates in the integral mode. Thus, we have
With the initial conditions at
where W = catalyst weight FC2H6 = volumetric flow-rate of ethane in the feed.
The rate of depletion of ethane \((r_{{{\text{C}}_{2} {\text{H}}_{6} }} )\) is given by
where r1 is the rate of oxidation of ethane to ethylene; r2 is the rate of oxidation of ethane to acetic acid; r3 is the rate of formation of CO by ethane oxidation and r4 is the rate of formation of CO2 by ethane oxidation.
The rate of formation of CO \((r_{\text{CO}} )\) is given by
where \(r_{5}\) is the rate of formation of CO from ethylene oxidation, \(r_{8}\) is the rate of formation of CO from acetic acid oxidation, and \(r_{10}\) is the rate of oxidation of CO to CO2.
The rate of formation of CO2 \((r_{{{\text{CO}}_{2} }} )\) is given by
where \(r_{6}\) is the rate of formation of CO2 from ethylene oxidation, and \(r_{9}\) is the rate of formation of CO2 from acetic acid oxidation, and \(r_{10}\) is the rate of oxidation of CO to CO2.
The rate of formation of CH3COOH (rCH3OOH) is given by
where \(r_{7}\) is the rate of formation of acetic acid from ethylene.
The expressions for r1–r10 are given by
with the k’s are defined as
where: all pressures are in atmosphere, the rates are in Ncm3/min g, k i ’s units can be deduced. K i ’s are dimensionless, \(E_{i}\) is activation energy in cal/mol, \(T\) is in K and \(R = 1.98\) cal/mol K. The dependence of the adsorption constants (K’s) on temperature was insignificant.
The calculated conversions are obtained by integrating equations subject to the initial condition. The pressures of the different components in terms of conversions are obtained as follows.
Let \(X_{1}\) be the conversion of ethane, \(X_{2}\) number of moles of CO produced per mole of ethane, \(X_{3}\) number of moles of CO2 produced per mole of ethane X4 number of mole of CH3COOH produced per mole of ethane.,
Y1 moles of O2 per mole of ethane and \(nt\) total number of moles of products per mole of ethane. Thus
Thus, the partial pressures of the different compounds are as follows:
where \(P_{t} = {\text{total pressure}}\)
Parameter estimation was done by the minimization of weighted residual squares of the conversions. The resulting non-linear equations are solved by Marquardt method. The method used is described in details by Hosten et al. [14].
The computer package for carrying out estimation calculations provides statistical data. They are F value for judging goodness of fit. t values for significance of parameters and 95% confidence regions of the parameter estimation and some others.
If the F value is small, the fit of the model is poor. At fixed number of parameters and experiments the higher is the F value the better is the fit. If a parameter is found to have very small t value, it means that it does not have a significant effect in determining the reaction rates and should be eliminated from the rate equations.
Parameter estimation results
Preliminary screening led us to conclude that the rate of oxidation of ethane to acetic acid and CO is negligibly small, i.e., \(r_{2} = 0\) and \(r_{3} = 0\). In addition, the rate of oxidation of ethylene to CO and CO2 is negligibly small, i.e., \(r_{5} = 0\) and \(r_{6} = 0\). We obtain the following parameter estimates (Table 2).
Discussion
We were not able to fit our experimental results with the model of Rahman et al. [9]. Thus, the model needs modifications for our case since we are using a different catalyst. Parameter estimation results indicate that the main reactions taking place are the oxidation of ethane to ethylene and CO2 and the oxidation of ethylene to acetic acid. Some of acetic acid is also oxidized to CO and CO2. CO is oxidized to CO2. The final reaction network takes the form of Fig. 2.
The results indicate that oxidation reactions take place on different sites. Ethane oxidation takes place on one site. Ethylene oxidation in presence of water to acetic acid takes place on another site. CO is oxidized to CO2 on a third site.
The estimation results indicate that a redox mechanism is more suitable to explain the reactions taking place for the oxidation of ethane to ethylene. Equation (14) for ethane oxidation to ethylene takes the familiar form of redox kinetics. The parameters associated with the oxidation of ethane to ethylene are determined with large t values giving high confidence in their values. For the oxidation of ethane to ethylene, our results agree with Linke et al. [7]. For the oxidation of ethylene to acetic acid our results also agree with that of Linke et al. [7] that acetic acid is related to the formation of a hydroxyl group from water on the active site. There are some controversy of whether ethane or ethylene oxidizes to CO and CO2. Burch and Swanakar [12] suggested that CO and CO2.are formed to a large extent from ethane. On the contrary Thorsteinson et al. [1] suggested that ethylene oxidizes into CO and CO2. Our results suggest that CO2 is mainly obtained from ethane and acetic acid whereas CO is obtained from acetic acid. This CO is then oxidized to CO2 on sites containing Pd. For reactions involving ethylene, acetic acid and CO a LH kinetic model is suggested.
For the oxidation of ethylene to acetic acid, Rahman et al. [9] suggested two routes one involves water and the second does not involve water. We found that acetic acid formation from the second route is not significant. The first route involving water needs modification. Let Y be a catalytic center for ethylene oxidation to acetic acid. Consider the following scheme;
This leads to the following rate equation
For large k this equation simplifies to Eq. (20).
For the oxidation of acetic acid to CO and CO2, we suggest the following scheme
which leads to the following rate equations
For large k, these two equations can be approximated by Eqs. (21) and (22).
For the oxidation of CO to CO2 this occurs through adsorbed oxygen and adsorbed CO on a site denoted by Z such that
For large k this equation can be approximated by Eq. (23).
The order of magnitude of the activation energies is similar to the values reported in the literature [1, 8, 9]. Large activation energy means that the reaction is favored at high temperature. Thus, it is noticed from the values of activation energy and from experimental results that higher temperature favors the formation of CO2 from oxidation of ethane.
Figures 3, 4 and 5 show the goodness of fit between the observed ethane conversion, CO2 yield and acetic acid yield and that calculated from the kinetic model. Reasonable fit is obtained for high conversion which would be of practical importance. For low conversion, the fit is slightly distorted possibly because inaccuracy in the experimental results.
Conclusion
Our kinetic model has a phenomenological and statistical significance. We assumed a general network, but found through model development and parameter estimation of model parameters that the significant reactions are ethane oxidation to ethylene and CO2, ethylene to acetic acid, CO and CO2, and CO to CO2. The analysis showed that ethane oxidation to acetic acid and CO is not significant. The reaction takes place on three sites, one for ethane oxidation to ethylene, the second for ethylene oxidation to acetic acid, CO and CO2, and the third for the oxidation of CO to CO2 In spite of reasonable acetic acid selectivity, there is still scope for catalyst improvement because of the high CO2 selectivity.
Future work can be directed towards deep understanding of all steps, such as ethylene oxidation, acetic acid oxidation, CO oxidation with these experiments carried out in the presence and absence of water vapor. The species which lead to the partial oxidation products should be identified.
References
Thorsteinson EM, Wilson TP, Young FG, Kasai PH (1978) The oxidative dehydrogenation of ethane over catalysts containing mixed oxides of molybdenum and vanadium. J Catal 52:116–132
Fakeeha AH, Soliman MA, Al-Rashidy MO, Ibrahim AA (1999) Sensitivity analysis of oxydehydrogenation of ethane in fixed bed reactor. Afinidad 56(484):369–376
Fakeeha AH, Al-Zeghayer YS, Soliman MA (2001) Sensitivity analysis for partial oxidation of ethane in a fluidized bed reactor. JEAS 48(2):385–399
Al-Sherehy FA, Adris AM, Soliman MA, Hughes R (1998) Avoidance of flammability and temperature runaway during oxidative dehydrogenation using a distributed feed. Chem Eng Sci 53(23):3965–3976
Karim K, Mamedov E, Al-Hazmi MH, Fakeeha AH, Soliman MA, Al-Zeghayer YS, Al-Fatish AS, Al-Arify AA (2000) Catalysts for producing acetic acid from ethane oxidation, processes for making same and method of using same. US Patent No 6, 030, 920, Feb 29
Fakeeha AH, Fahmy YM, Soliman MA, Alwahabi SM (2000) A kinetic model for partial oxidation of ethane to acetic acid on promoted VPO catalyst. J Chem Technol Biotechnol 75:1160
Linke D, Wolf D, Dingerdissen U, Zeyβ BS, Baerns M (2002) Catalytic partial oxidation of ethane to acetic acid over Mo1 V0.25 Nb0.12 Pd0.0005 O x part I: catalyst performance and reaction mechanism. J Catal 205:16–31
Linke D, Wolf D, Dingerdissen U, Zeyβ BS, Baerns M (2002) Catalytic partial oxidation of ethane to acetic acid over Mo1 V0.25 Nb0.12 Pd0.0005 O x part II: kinetic modeling. J Catal 205:32–43
Rahman F, Loughlin KF, Al-Saleh MA, Saeed MR, Tukur NM, Hossain MM, Karim K, Mamedov A (2010) Kinetics and mechanism of partial oxidation of ethane to ethylene and acetic acid over MoV type catalysts. Appl Catal A Gen 375(1):17–25
Li X, Iglesia E (2008) Support and promoter effects in the selective oxidation of ethane to acetic acid catalyzed by Mo–V–Nb oxides. Appl Catal A Gen 334(1–2):339–347
Li X, Iglesia E (2008) Kinetics and mechanism of ethane oxidation to acetic acid on catalysts based on Mo–V–Nb oxides. J Phys Chem 112(38):15001–15008
Burch R, Swarnakar R (1991) Oxidative dehydrogenation of ethane on vanadium-molybdenum oxide and vanadium-niobium-molybdenum oxide catalysts. Appl Catal A Gen 70:129–148
Berlowitz PJ, Peden CHF, Goodman DW (1988) Kinetics of carbon monoxide oxidation on single-crystal palladium, platinum, and iridium. J Phys Chem 92:5213–5221
Hosten LH, Froment GF (1975) Parameter estimation in multiresponse models. Period Polytech Chem Eng 19:123–136
Al-Zeghayer YS, Al-Awadi AS, Jibril BY, Soliman MA, Al-Mayman S (2013) Partial oxidation of ethane to acetic acid on titania supported MoVNbPd catalyst. Asian J Chem 25:7979–7985
Al-Mayman S, Al-awadi AS, Al-Zeghayer YS, Soliman MA (2018) Partial oxidation of ethane to acetic acid using metallic Pd promoted MoVNb catalyst supported on titania. J Chem Soc Pak 40(1):76–88
Acknowledgements
The authors acknowledge the financial support provided by King Abdul Aziz City for Science and Technology (KACST) for this research under Grant number AR-29-256.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Al-Mayman, S.I., Soliman, M.A., Al-Awadi, A.S. et al. Reaction kinetics of ethane partial oxidation to acetic acid. Appl Petrochem Res 8, 29–38 (2018). https://doi.org/10.1007/s13203-018-0195-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13203-018-0195-8