Next Article in Journal
Thermodynamic Features of the Intensive Formation of Hydrocarbon Hydrates
Previous Article in Journal
Impact of the Drying Temperature and Grinding Technique on Biomass Grindability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Propane Pre-Reforming into Methane-Rich Gas over Ni Catalyst: Experiment and Kinetics Elucidation via Genetic Algorithm

by
Sergey I. Uskov
1,2,
Dmitriy I. Potemkin
1,2,3,*,
Leniza V. Enikeeva
2,4,
Pavel V. Snytnikov
1,2,
Irek M. Gubaydullin
5,6 and
Vladimir A. Sobyanin
1
1
Boreskov Institute of Catalysis, Pr. Lavrentieva 5, Novosibirsk 630090, Russia
2
Energy-Saving Catalytic Processes Laboratory, Novosibirsk State University, Pirogova St. 2, Novosibirsk 630090, Russia
3
Department of Environmental Engineering, Novosibirsk State Technical University, Karl Marx Pr. 20, Novosibirsk 630073, Russia
4
Mathematical Department, Ufa State Petroleum Technological University, Kosmonavtov St. 1, Ufa 450062, Russia
5
Technological Department, Ufa State Petroleum Technological University, Kosmonavtov St. 1, Ufa 450062, Russia
6
Institute of Petrochemistry and Catalysis of RAS, October Pr. 141, Ufa 450075, Russia
*
Author to whom correspondence should be addressed.
Energies 2020, 13(13), 3393; https://doi.org/10.3390/en13133393
Submission received: 21 May 2020 / Revised: 16 June 2020 / Accepted: 30 June 2020 / Published: 2 July 2020
(This article belongs to the Section A5: Hydrogen Energy)

Abstract

:
Pre-reforming of propane was studied over an industrial nickel-chromium catalyst under pressures of 1 and 5 bar, at a low steam to carbon molar ratio of 1, in the temperature range of 220–380 °C and at flow rates of 4000 and 12,000 h−1. It was shown that propane conversion proceeded more efficiently at low pressure (1 atm) and temperatures above 350 °C. A genetic algorithm was applied to search for kinetic parameters better fitting experimental results in such a wide range of experimental conditions. Power law and Langmuir–Hinshelwood kinetics were considered. It was shown that only Langmuir–Hinshelwood type kinetics correctly described the experimental data and could be used to simulate the process of propane pre-reforming and predict propane conversion under the given reaction conditions. The significance of Langmuir–Hinshelwood kinetics increases under high pressure and temperatures below 350 °C.

1. Introduction

Despite the growing contribution of renewable energy sources, such as solar/wind energy or biomass, and respective technologies and devices for energy generation and storage, it is believed that for many decades, fossil resources (primarily oil and natural gas) will still form the basis of global energy consumption. An increase in the efficiency of using hydrocarbon raw materials while improving operational characteristics (an important factor for the final consumer) is associated with the use of fuel cells (FCs) and hydrogen storage technologies for the needs of distributed energy [1,2]. Therefore, fuel cells, which allow for converting the energy of chemical bonds into electrical energy with high efficiency even in compact and mobile systems, are one of the most dynamically developing energy sectors.
At the moment, the infrastructure for supplying hydrogen and logistics are not well developed to meet mass consumption needs: there is a lack of hydrogen gas stations, and long-term storage of large quantities of hydrogen remains problematic. In this regard, most of the world research centers involved in the development of fuel cell based power plants show interest in the use of widely available traditional logistic fuels (gasoline, propane-butane, kerosene and diesel fuel). Besides, in comparison with other hydrogen storage technologies available today, hydrocarbon fuels contain the largest amount of hydrogen per unit volume [3]. Therefore, the use of liquid or easily liquefied hydrocarbon fuels with high energy density is relevant for fuel cells. However, such fuels cannot be directly oxidized in the anode space of the fuel cell, since they initiate the carbonization of the electrode and the failure of the fuel cell. In this case, the fuels are first converted to a hydrogen-containing gas and then oxidized in the fuel cell. The processes of steam conversion or steam reforming (SC), autothermal reforming (ATR) and partial oxidation (PO) are used to obtain synthesis gas from hydrocarbons. The most effective process is SC, but it requires high operating temperatures, consumes a significant amount of heat to evaporate water and preheats steam and needs a water purification system. Therefore, SC can be implemented only for stationary power plants based on high-power fuel cells. PO of liquid hydrocarbon fuels has not yet been implemented due to thermodynamic limitations, an undesirable side process of carbon formation and high operating temperatures.
In order to solve the problem of partial oxidation of C2–C10 hydrocarbons, reach high performance and minimize the amount of water to be vaporized, it is proposed to combine the pre-reforming of hydrocarbons and the partial oxidation of the resulting methane-containing mixture into synthesis gas. The first stage is carried out at low molar ratios of H2O/C = 0.7–1.5 and temperatures of 250–350 °C. At these conditions, the main reaction products are methane, hydrogen and carbon dioxide. The second stage (PO) provides a high performance for the overall process of synthesis gas production.
Although catalytic steam conversion of hydrocarbons has been studied quite well, in the vast majority of works, the process was studied at temperatures of 450–700 °C and a molar ratio of H2O/C > 2 to obtain synthesis gas directly from hydrocarbons. SC of C2+-hydrocarbons at H2O/C < 1 and temperatures below 350 °C are poorly studied. However, there is a limited number of works [4,5,6,7] showing the feasibility of such an approach to obtain methane-hydrogen mixtures.
Nickel-based catalysts are usually used in hydrocarbons steam reforming because of their high activity and low cost compared to platinum metals. Currently, nickel catalysts for methanation of carbon oxides, the pre-reforming of C2+-hydrocarbons and the steam conversion of natural gas to synthesis gas are available commercially. Natural gas steam reforming catalysts possess high mechanical strength, a Ni content of about 10–15 wt.% and can be used at temperatures of 700–900 °C. As a result, the catalysts for natural gas steam conversion have a low nickel surface area and low catalytic activity at low temperatures. Conventional pre-reforming catalysts contain 15–30 wt.% of Ni and can be used at temperatures of 450–600 °C. These catalysts were developed with particular attention to carbonization resistance [8] rather than catalytic activity. CO/CO2 methanation catalysts have a high nickel content and dispersion but lack thermal stability and are usually operated at temperatures of 300–320 °C. Since the pre-reforming of hydrocarbons at such low steam to carbon molar ratios is carried out at about 250–350 °C [5], the use of methanation catalysts seems to be the optimal choice.
There are a lack of data on the reaction kinetics for pre-reforming at such low temperatures and the steam to carbon molar ratio. Ethane steam conversion was studied in methane excess at molar ratios of CH4/C2H6 ≈ 20 and H2O/CC2+ = 0.6−6 over a Ni-Cr catalyst at 300–360 °C, 1 atm pressure [9]. It was shown that ethane conversion rate could be described by a Langmuir–Hinshelwood type equation:
WC2H6 = k∙PC2H6/(1 + A∙PH2O/PH2)
where k and A are constants dependable of temperature; PC2H6, PH2O and PH2 are partial pressures of the components. The equation corresponds to a scheme with ethane adsorption being the rate-determining step.
A two-stage macrokinetic scheme was suggested for the low-temperature steam conversion of methane-propane mixtures over Ni-based catalysts [10]. Conversion of the methane-pentane mixtures under 1 bar pressure was simulated using the scheme, which included irreversible reactions of the steam conversion of C2+-hydrocarbons, Equation (2), with the formation of CO2 and H2, followed by reversible CO2 methanation, Equation (3):
CnH2n+2 + 2nH2O → (3n + 1)H2 + nCO2 (n > 1),
CO2 + 4H2 ⇄ CH4 + 2H2O.
Kinetic parameters were determined: CO2 methanation was shown to occur in a quasi-equilibrium mode at T > 250 °C. The apparent reaction order with respect to C3H8 was close to one. The effect of the CH4, CO2, H2 and H2O concentrations on the C3H8 conversion was insignificant. It should be noted that the kinetics are valid for C2+-hydrocarbons partial pressure much lower that 1 bar. Thus, kinetics for a wider pressure range are required.
The kinetic study of heterogenous catalytic reactions is generally a challenging task as mechanisms are typically multistep and the majority of experimental data are obtained in a plug-flow reactor in which the reaction conditions and reaction rates are different at all points of the reactor. The kinetic parameters are determined by solving an inverse problem of chemical kinetics.
The pre-reforming of the propane reaction is a complex chemical reaction, in other words, a multi-stage one, and for each stage of the reaction, it is necessary to find the kinetic parameters (activation energies and preexponential factors). Moreover, as will be shown later, we will look for additional reaction parameters when solving the inverse problem. Thus, an ill-posed optimization problem with high dimensionality is formulated. For these reasons, an efficient multidimensional optimization technique such as genetic algorithms (GAs) are, therefore, necessary. Genetic algorithms were designed for the accelerated search of a global minimum in various optimization problems [11,12,13,14]. It was successfully applied for the kinetics elucidation of [15,16,17,18,19].
In this work, the experiments on propane pre-reforming were carried out over an industrial nickel methanation catalyst NIAP-07-05, and a mathematical simulation was applied to describe the process. The effect of temperature and pressure was studied. The experiments data were fitted by the power law (PL) and Langmuir–Hinshelwood (LH) kinetics.

2. Materials and Methods

2.1. Catalytic Experiments

Propane pre-reforming was studied using an NIAP-07-05 industrial catalyst of the following composition (wt.%): 42 NiO, 12 Cr2O3, 46 Al2O3 and 4 graphite. This catalyst was previously used in the studies on the low-temperature steam conversion of light hydrocarbons in methane excess [5,10]. The catalyst advantages are high activity and stability in propane reforming and resistance to carbon formation, which was confirmed by experiments during more than 40 h.
The experiments were carried out in a fixed-bed U-shaped quartz reactor (i.e., 8 mm, catalyst bed volume 1.7 mL). Experimental conditions were as follows: 1–5 bar pressure, 220–380 °C, gas hourly space velocity (GHSV) is 4000 and 12,000 h−1. A measure of 2 g of the catalyst (fraction 0.5–1 mm) was previously reduced at 420 °C in a flow of 10 mL/min H2 and 100 mL/min Ar. The composition of the reaction mixture was 25 vol.% C3H8 and 75 vol.% H2O, which corresponded to H2O/C molar ratio of 1.
The temperature of the catalyst was measured directly in the middle of the catalyst bed at the centerline point with the use of a K-type thermocouple. The GC analysis was applied to determine the composition of the gas phase before and after the reactor. Components were determined using a Chromos GC-1000 chromatograph equipped with two thermal conductivity detectors (TCD) and one flame ionization detector (FID) detector with a methanator. CH4, CO2 and H2 were registered as main reaction products. The concentration of carbon monoxide was negligible (below 0.1 vol.%) and, therefore, CO was not taken into account when processing the results. The carbon imbalance in all the experiments did not exceed ±1 rel. %. The obtained results were compared with equilibrium concentrations, calculated using the HSC Chemistry 7.1 software package (Outotec, Finland).

2.2. Kinetic Modeling

Based on the experimental dependencies, the inverse chemical kinetics task was solved to determine the reaction rate constants and activation energies. These parameters were obtained by solving the direct task numerous times by enumerating the rate constants of the stages (or a set of pre-exponents and activation energies). The model of an iso-thermal plug-flow rector was used for solving the direct problem. The task was to minimize the functional Equation (4) of the deviation of the experimental data from the calculated data:
F = i = 1 M j = 1 N | x i j c a l c x i j e x p | m i n ,
where xijcalc are the calculated values, xijexp are the experimental data, M is the number of experimental points (temperatures) and N is the number of substances involved in the reaction.
When solving the direct task, the concentrations of the reaction components were calculated. The results were compared with experimental data, which mathematically meant minimizing the difference between the experimental and calculated data. A genetic algorithm was chosen as a method for solving the inverse task. The authors developed a program in the Matlab environment using the Optimization toolbox (see Supplementary Materials). Population sizes between 50 and 250 were used. As a rule, the higher the population size, the faster the algorithm converges to a solution, but this greatly affects the calculation time. The algorithm terminates when there is no further improvement in fitness function F (see Equation (4)) after several hundred generations.
Four kinetics constants were optimized when solving the inverse task: Eref and Emet (kJ/mol), and kref and kmet. These parameters were limited according to the concepts of the kinetics of low-temperature steam conversion of C2+-hydrocarbons in methane excess [8]: 90 ≤ Eref ≤ 150; 30 ≤ Emet ≤ 60; 1010 ≤ kref ≤ 1013; 104 ≤ kmet ≤ 107. Also, the apparent reaction order with respect to propane was varied within the range −1 < n < 1 for PL kinetics. Constants m and B (see Equation (10)) were varied for LH kinetics.

3. Results and Discussion

Propane pre-reforming experiments were carried out at 220–380 °C. The main reaction products were methane, hydrogen and carbon dioxide, as well as residual propane. Figure 1a–d shows the temperature dependences of the outlet concentrations of propane, methane, carbon dioxide and hydrogen during propane pre-reforming under pressures 1 bar and 5 bar. The flow rates (GHSV) for both pressures were 4000 and 12,000 h−1. It is seen that the conversion of propane proceeds more efficiently under 1 bar pressure. Complete propane conversion can be reached at temperatures above 350 °C.
It can be seen that at 1 bar pressure and a flow rate of 12,000 h−1, there is a maximum of hydrogen concentration at 290 °C, followed by a decrease to a minimum at 320 °C, and then a further increase. This can be explained that at low temperatures, hydrogen and CO2 were accumulating during the steam conversion reaction. With increasing temperature, hydrogen is consumed by the hydrogenolysis of propane. At high temperatures, the hydrogen concentration is determined by the equilibrium distribution of the products and increases according to the equilibrium of the CO2 methanation reaction.
The inverse problem was then solved for the experiments on propane pre-reforming. The reaction scheme consisted of two reactions: steam conversion of propane and methanation of CO2.
C3H8 + 6H2O → 10H2 + 3CO2,
CO2 + 4H2 = CH4 + 2H2O.
Earlier, it was shown that this approach can be quite effective in the case of the low-temperature steam conversion of small concentrations of propane in methane excess [10]. In this work, we use the following form of the kinetic equations:
W ref = k ref · exp ( E ref RT ) · C C 3 H 8 n ,
W met = k met · exp ( E met RT ) · C H 2 · [ 1 P CH 4 P H 2 O 2 K eq P CO 2 P H 2 4 ] ,
where Wref and Wmet are the reaction rates (mol∙m−3∙sec−1); Eref and Emet are the observed activation energies (J/mol); kref and kmet are pre-exponential factors ((mol∙m3)n−1∙sec−1 and sec−1). The indices ref and met refer to the pre-reforming and methanation reactions, respectively. CC3H8 and CH2 are concentrations of propane and hydrogen (mol/m3); n is effective reaction order with respect to propane, which ranged from −2 to 2; Keq is the equilibrium constant for CO2 methanation reaction and PCH4, PH2O, PCO2, PH2 are partial pressures of the corresponding substances (bar).
As a result of solving the inverse task, the values of E ref , E met , k ref , k met and n were optimized. The obtained values are listed in Table 1.
These parameters were substituted into Equations (7) and (8). Next, within the framework of the proposed reaction scheme and the model of a quasi-homogeneous isothermal plug flow reactor, the system of the material balance equations was solved to calculate the temperature dependences of the reaction product concentrations. The system of differential equations is written in the following form:
G dy i dz = ν i W ref ρ i .
Here, G is the mass flow of the gas mixture (kg∙m−2∙sec−1); yi is the mass fraction of the i-th component, i = { C 3 H 8 , CH 4 , H 2 O , H 2 , CO 2 } ; z is the coordinate along the catalyst bed (0 ≤ z ≤ 0.030 m); νi are stoichiometric coefficients for steam reforming reactions; Wref is the steam conversion rate of propane (mol∙m−3∙s−1) and ρi is the density of the i-th component (kg∙m−3). The boundary conditions were expressed as: yi = yi0 at z = 0.
The calculation results for the experiments under 1 and 5 bar and flow rates of 4000 and 12,000 h−1 are shown in Figure 2a–d. It is seen that power low type kinetics correctly describes the experimental data at low temperatures for both pressures, while for high pressure, it predicts the reaching the equilibrium too fast. Note that the effective reaction order with respect to propane (n = 0.15) in the case of the pre-reforming of pure propane is less than or close to zero, in contrast to the conversion of diluted propane-methane mixtures [10]. The reason for this is the high coverage of the catalyst surface with propane molecules (θC3H8), which, according to the Langmuir–Hinshelwood model, reduces the effective reaction order.
Aside from the power law dependence, attempts were made to find the solution of the inverse task in the form of the Langmuir–Hinshelwood (LH) model:
Wref = k ref · exp ( E ref RT ) · C C 3 H 8 ( 1 + B · C C 3 H 8 ) m .
The value of B constant varied from 0 to 5, and the value of m varied from 0 to 2. When solving the inverse problem, the values of E ref , E met , k ref and k met also varied (the upper and lower restrictions were the same as for the power law model). The results are shown in Table 2.
These values were used to solve the direct task. Figure 3a–d shows the results of the calculations for experiments on propane pre-reforming under 1 bar and 5 bar.
It is seen that the model based on the Langmuir–Hinshelwood mechanism correctly describes the experimental data. It is seen that apparent activation energies for both fits are quite close. The main difference is in the nature of the dependence of the rate of conversion of propane on its concentration.
Figure 4 shows the temperature dependence of the denominator of Equation (10), which shows how the value of [1 + B∙CC3H8]m varies with temperature and pressure (with parameters B and m selected). It is seen that the denominator varies within the range of 1–3 under the conditions studied, and decreases with pressure decrease and temperature increase. The results obtained explain the inconsistency of the power law kinetics: at a low temperature, the apparent reaction order with respect to propane is close to zero, providing good agreement of experiment and fitting; however, with temperature increase, the apparent propane reaction order becomes closer to one, but power law kinetics cannot reflect this, therefore making it inappropriate.

4. Conclusions

The temperature dependences of the concentrations of the products for propane pre-reforming were obtained under 1 and 5 bar, H2O/C = 1 ratio, in the temperature range of 220–380 °C and at flow rates of 4000 and 12,000 h−1. It was shown that propane conversion over the industrial NIAP-07-05 catalyst (42% NiO, 12% Cr2O3, 46% Al2O3, 4% C) proceeds more efficiently at low pressure (1 atm) and temperatures above 350 °C. Mathematical modeling of the experimental dependencies was carried out according to two models: a power law model and the Langmuir–Hinshelwood mechanism. A genetic algorithm was used to obtain the best fit within the models used. It was shown that the LH model correctly describes the experimental data and can be used to simulate the process of propane pre-reforming and predict propane conversion under given reaction conditions. The significance of the LH model increases under high pressure and temperatures below 350 °C.

Supplementary Materials

The following are available online at https://www.mdpi.com/1996-1073/13/13/3393/s1, Supplementary materials for Propane Pre-Reforming into Methane-Rich Gas over Ni Catalyst: Experiment and Kinetics Elucidation via Genetic Algorithm.

Author Contributions

Conceptualization, D.I.P.; methodology, P.V.S.; software, L.V.E.; validation, I.M.G.; formal analysis, D.I.P.; investigation, D.I.P.; resources, D.I.P.; data curation, S.I.U.; writing—original draft preparation, D.I.P., S.I.U. and L.V.E.; supervision, V.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

In the part of experimental study the work is supported by the RFBR project 19-33-60008 “Perspectiva” (D.I.P.). In the part of mathematical modelling the work is supported by the RFBR project 19-37-60014 “Perspectiva” (L.V.E.). S.I.U. appreciates support by the Ministry of Science and Higher Education of the Russian Federation (Project AAAA-A17-117041710088-0).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ren, J.; Musyoka, N.M.; Langmi, H.W.; Mathe, M.; Liao, S. Current research trends and perspectives on materials-based hydrogen storage solutions: A critical review. Int. J. Hydrog. Energy 2017, 42, 289–311. [Google Scholar] [CrossRef]
  2. Zhang, F.; Zhao, P.; Niu, M.; Maddy, J. The survey of key technologies in hydrogen energy storage. Int. J. Hydrog. Energy 2016, 41, 14535–14552. [Google Scholar] [CrossRef]
  3. Shekhawat, D.; Berry, D.A.; Spivey, J.J. Introduction to Fuel Processing. In Fuel Cells: Technologies for Fuel Processing; Elsevier: Amsterdam, The Netherlands, 2011; pp. 1–9. [Google Scholar] [CrossRef]
  4. Christensen, T.S. Adiabatic prereforming of hydrocarbons–an important step in syngas production. Appl. Catal. A Gen. 1996, 138, 285–309. [Google Scholar] [CrossRef]
  5. Zyryanova, M.M.; Snytnikov, P.V.; Amosov, Y.I.; Belyaev, V.D.; Kireenkov, V.V.; Kuzin, N.A.; Vernikovskaya, M.V.; Kirillov, V.A.; Sobyanin, V.A. Upgrading of associated petroleum gas into methane-rich gas for power plant feeding applications. Technological and economic benefits. Fuel 2013, 108, 2822–2891. [Google Scholar] [CrossRef]
  6. Rostrup-Nielsen, J.R.; Christiansen, L.J. Concepts in Syngas Manufacture; Imperial College Press: London UK, 2011; p. 379. [Google Scholar]
  7. Rostrup-Nielsen, J.R. Catalytic Steam Reforming; Springer: Berlin, Germany, 1984; p. 119. [Google Scholar]
  8. Rostrup-Nielsen, J.R.; Christensen, T.; Dybkjaer, I. Steam reforming of liquid hydrocarbons. Recent. Adv. Basic Appl. Asp. Ind. Catal. 1998, 113, 81–95. [Google Scholar]
  9. Meshenko, N.T.; Veselov, V.V.; Shub, F.S.; Temkin, M.I. Kinetics of the low-temperature steam reforming of ethane on a nickel–chromium catalyst. Kinet. Catal. 1977, 18, 962–967. [Google Scholar]
  10. Uskov, S.I.; Potemkin, D.I.; Shigarov, A.B.; Snytnikov, P.V.; Kirillov, V.A.; Sobyanin, V.A. Low-temperature steam conversion of flare gases for various applications. Chem. Eng. J. 2019, 368, 5335–5340. [Google Scholar] [CrossRef]
  11. Bhoskar, T.; Kulkarni, N.K.; Patekar, S.; Bhoskar, T.; Kulkarni, O.; Kakandikar, G.M.; Nandedkar, V.M. Genetic Algorithm and its Applications to Mechanical Engineering: A Review. Mater. Today Proc. 2015, 2, 1624–2630. [Google Scholar] [CrossRef]
  12. Wiśniewski, J. Optimal Design of Reinforcing Fibres in Multilayer Composites using Genetic Algorithms. Fibres Text. East. Eur. 2004, 12, 90–543. [Google Scholar]
  13. Lam RatnaRaju, T.K. Design Optimization of Perforated Plate Heat Exchangers using Genetic Algorithm. Adv. Mater. Manuf. Charact. 2013, 3, 89–93. [Google Scholar] [CrossRef] [Green Version]
  14. Chakraborti, N.; Kumar, R.; Jain, D. A study of the continuous casting mold using a pareto-converging genetic algorithm. Appl. Math. Model. 2001, 25, 2872–2897. [Google Scholar] [CrossRef]
  15. Elliott, L.; Ingham, D.B.; Kyne, A.G.; Mera, N.S.; Pourkashanian, M.; Wilson, C.W. Genetic algorithms for optimisation of chemical kinetics reaction mechanisms. Prog. Energy Combust. Sci. 2004, 30, 297–328. [Google Scholar] [CrossRef]
  16. Wolf, D.; Moros, R. Estimating rate constants of heterogeneous catalytic reactions without supposition of rate determining surface steps—An application of a genetic algorithm. Chem. Eng. Sci. 1997, 52, 1189–1199. [Google Scholar] [CrossRef]
  17. Garcia, S. Experimental Design Optimization and Thermophysical Parameter Estimation of Composite Materials Using Genetic Algorithms. Ph.D. Thesis, Virginia Polytechnic Institute (and Universite de Nantes), Nantes, France, 1999. [Google Scholar]
  18. Lautenberger, C.; Rein, G.; Fernandez-Pello, C. The application of a genetic algorithm to estimate material properties for fire modeling from bench-scale fire test data. Fire Safe. J. 2006, 41, 204–214. [Google Scholar] [CrossRef] [Green Version]
  19. Rein, G.; Lautenberger, C.; Fernandez-Pello, A.C.; Torero, J.L.; Urban, D.L. Application of genetic algorithms and thermogravimetry to determine the kinetics of polyurethane foam in smoldering combustion. Combust. Flame 2006, 146, 95–108. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Temperature dependences of the outlet concentrations on a dry basis of propane (a), methane (b), carbon dioxide (c) and hydrogen (d) during propane pre-reforming. Experimental conditions: 220–380 °C, GHSV 4000 and 12,000 h−1, 1 and 5 bar pressure and reaction mixture: 25 vol.% C3H8, 75 vol.% H2O. Concentrations of the gas components on the figure are given on the dry basis.
Figure 1. Temperature dependences of the outlet concentrations on a dry basis of propane (a), methane (b), carbon dioxide (c) and hydrogen (d) during propane pre-reforming. Experimental conditions: 220–380 °C, GHSV 4000 and 12,000 h−1, 1 and 5 bar pressure and reaction mixture: 25 vol.% C3H8, 75 vol.% H2O. Concentrations of the gas components on the figure are given on the dry basis.
Energies 13 03393 g001
Figure 2. Temperature dependences of the outlet concentrations of propane, methane, CO2 and hydrogen in propane pre-reforming. Points are the experiment, and lines are the modeling. Experimental conditions: 25 vol.% C3H8, 75% H2O and pressure 1 bar. The solution of the inverse problem is power law dependence with a propane effective reaction order of 0.15. Concentrations of the gas components on the figure are given on the dry basis. Figure 2a shows the results for the experiments under 1 bar and 4000 h−1, Figure 2b—1 bar and 12,000 h−1, Figure 2c—5 bar and 4000 h−1, Figure 2d—5 bar and 12,000 h−1.
Figure 2. Temperature dependences of the outlet concentrations of propane, methane, CO2 and hydrogen in propane pre-reforming. Points are the experiment, and lines are the modeling. Experimental conditions: 25 vol.% C3H8, 75% H2O and pressure 1 bar. The solution of the inverse problem is power law dependence with a propane effective reaction order of 0.15. Concentrations of the gas components on the figure are given on the dry basis. Figure 2a shows the results for the experiments under 1 bar and 4000 h−1, Figure 2b—1 bar and 12,000 h−1, Figure 2c—5 bar and 4000 h−1, Figure 2d—5 bar and 12,000 h−1.
Energies 13 03393 g002aEnergies 13 03393 g002b
Figure 3. Temperature dependences of the outlet concentrations of propane, methane, CO2 and hydrogen in propane pre-reforming. Points are the experiment, and lines are the modeling. Experimental conditions: 25 vol.% C3H8, 75% H2O and pressure 1 bar. The solution of the inverse problem is according to the Langmuir–Hinshelwood (LH) model. Concentrations of the gas components on the figure are given on the dry basis. Figure 3a shows the results for the experiments under 1 bar and 4000 h−1, Figure 3b—1 bar and 12,000 h−1, Figure 3c—5 bar and 4000 h−1, Figure 3d—5 bar and 12,000 h−1.
Figure 3. Temperature dependences of the outlet concentrations of propane, methane, CO2 and hydrogen in propane pre-reforming. Points are the experiment, and lines are the modeling. Experimental conditions: 25 vol.% C3H8, 75% H2O and pressure 1 bar. The solution of the inverse problem is according to the Langmuir–Hinshelwood (LH) model. Concentrations of the gas components on the figure are given on the dry basis. Figure 3a shows the results for the experiments under 1 bar and 4000 h−1, Figure 3b—1 bar and 12,000 h−1, Figure 3c—5 bar and 4000 h−1, Figure 3d—5 bar and 12,000 h−1.
Energies 13 03393 g003
Figure 4. Temperature dependence of the denominator of the Langmuir–Hinshelwood model for the experiments on propane pre-reforming under 1 bar and 5 bar and at flow rates of 4000 and 12,000 h−1.
Figure 4. Temperature dependence of the denominator of the Langmuir–Hinshelwood model for the experiments on propane pre-reforming under 1 bar and 5 bar and at flow rates of 4000 and 12,000 h−1.
Energies 13 03393 g004
Table 1. The values of E ref , E met , k ref , k met and n obtained by solving the inverse task.
Table 1. The values of E ref , E met , k ref , k met and n obtained by solving the inverse task.
Eref, kJ/molEmet, kJ/molkref, (mol∙m−3)n−1∙sec−1kmet, sec−1n
107.546.74.0∙10103.2∙1050.15
Table 2. The values of B, m, E ref , E met , k ref and k met optimized in the framework of the Langmuir–Hinshelwood model.
Table 2. The values of B, m, E ref , E met , k ref and k met optimized in the framework of the Langmuir–Hinshelwood model.
Eref, kJ/molEmet, kJ/molkref, (mol∙m−3)n−1∙sec−1kmet, sec−1mB
105.144.21.3∙10101.6∙1051.0670.204

Share and Cite

MDPI and ACS Style

Uskov, S.I.; Potemkin, D.I.; Enikeeva, L.V.; Snytnikov, P.V.; Gubaydullin, I.M.; Sobyanin, V.A. Propane Pre-Reforming into Methane-Rich Gas over Ni Catalyst: Experiment and Kinetics Elucidation via Genetic Algorithm. Energies 2020, 13, 3393. https://doi.org/10.3390/en13133393

AMA Style

Uskov SI, Potemkin DI, Enikeeva LV, Snytnikov PV, Gubaydullin IM, Sobyanin VA. Propane Pre-Reforming into Methane-Rich Gas over Ni Catalyst: Experiment and Kinetics Elucidation via Genetic Algorithm. Energies. 2020; 13(13):3393. https://doi.org/10.3390/en13133393

Chicago/Turabian Style

Uskov, Sergey I., Dmitriy I. Potemkin, Leniza V. Enikeeva, Pavel V. Snytnikov, Irek M. Gubaydullin, and Vladimir A. Sobyanin. 2020. "Propane Pre-Reforming into Methane-Rich Gas over Ni Catalyst: Experiment and Kinetics Elucidation via Genetic Algorithm" Energies 13, no. 13: 3393. https://doi.org/10.3390/en13133393

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