Experimental investigation of Methane Partial Oxidation for Hydrogen Production

Similar documents
A FLAMELET MODEL FOR DIFFUSION FLAMES IN POROUS MEDIA. Max Akira Endo Kokubun. Fernando Filho Fachini

Direct pore level simulation of premixed gas combustion in porous inert media using detailed chemical kinetics

NUMERICAL SOLUTION FOR THE COMBUSTION OF METHANE IN POROUS MEDIA

Fundamentals of Combustion

Reacting Gas Mixtures

INTRODUCTION TO CATALYTIC COMBUSTION

CFD study of gas mixing efficiency and comparisons with experimental data

COMBUSTION OF FUEL 12:57:42

Methane Oxidation Reactions

Interactions between oxygen permeation and homogeneous-phase fuel conversion on the sweep side of an ion transport membrane

Premixed filtration combustion of micron and sub-micron particles in inert porous media: A theoretical analysis

Simulation of Methanol Production Process and Determination of Optimum Conditions

Hydrogen addition to the Andrussow process for HCN synthesis

Introduction. Mathematical model and experimental

Combustion. Indian Institute of Science Bangalore

Laminar Premixed Flames: Flame Structure

KINETIC MODELING OF OXIDATION METHANE CONVERSION IN REGIME OF FILTRATION COMBUSTION WITH SUPERADIABATIC HEATING

Lecture 7 Flame Extinction and Flamability Limits

Partial Oxidation of Methane to Methanol in a Catalytc Packed Bed Reactor: Kinetic Modeling and Experimental Study

Lecture (9) Reactor Sizing. Figure (1). Information needed to predict what a reactor can do.

Chemical Reaction Engineering - Part 16 - more reactors Richard K. Herz,

DETAILED MODELLING OF SHORT-CONTACT-TIME REACTORS

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore

Exploring The Fundamentals In Catalytic Partial Oxidation Of Methane: The Interaction Between Diffusion And Reaction In A Packed Bed Reactor

AAE COMBUSTION AND THERMOCHEMISTRY

Combustion: Flame Theory and Heat Produced. Arthur Anconetani Oscar Castillo Everett Henderson

Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane

The Seeding of Methane Oxidation

MICRO SOLID OXIDE FUEL CELL ON THE CHIP

Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

Chemical Kinetics Computations of Fuel Decomposition to Aldehydes for NOx Reduction in Engine Exhaust Gas

Modeling, Simulation and Control of a Tubular Fixed-bed Dimethyl Ether Reactor

Effects of Different Processing Parameters on Divinylbenzene (DVB) Production Rate

Overview of Reacting Flow

MODELING AND SIMULATION OF REFORMER AUTO- THERMAL REACTOR IN AMMONIA UNIT

Laminar flame speed (burning velocity) reactants. propagating flame front. products. reactants

Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system

CFD Simulation of Catalytic Combustion of Benzene

3.2 Alkanes. Refining crude oil. N Goalby chemrevise.org 40 C 110 C 180 C. 250 C fuel oil 300 C 340 C. Fractional Distillation: Industrially

Overview of Turbulent Reacting Flows

Modeling of the Unsteady State Methanol Synthesis at the Level of Catalyst Pellet

Combustion Theory and Applications in CFD

Analysis of homogeneous combustion in Monolithic structures

C H A P T E R 3 FUELS AND COMBUSTION

CFD Modeling Ensures Safe and Environmentally Friendly Performance in Shell Claus Off-Gas Treating (SCOT) Unit

Reactor Modeling and Process Analysis for Partial Oxidation of Natural Gas

Department of Mechanical Engineering BM 7103 FUELS AND COMBUSTION QUESTION BANK UNIT-1-FUELS

While entry is at the discretion of the centre, candidates would normally be expected to have attained one of the following, or equivalent:

and mol of Cl 2 was heated in a vessel of fixed volume to a constant temperature, the following reaction reached equilibrium.

Contents. Preface... xvii

Thermal Energy Final Exam Fall 2002

Fuel, Air, and Combustion Thermodynamics

Fundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur

Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory

Well Stirred Reactor Stabilization of flames

Fundamentals Of Combustion (Part 1) Dr. D.P. Mishra Department of Aerospace Engineering Indian Institute of Technology, Kanpur

Chapter 8: Quantities in Chemical Reactions

Two-dimensional mathematical modeling of oxidative coupling of methane in a membrane reactor

Chemical Engineering 140. Chemical Process Analysis C.J. Radke Tentative Schedule Fall 2013

a 16 It involves a change of laminar burning velocity, widening or narrowing combustion limits for

The coupling effect of gas-phase chemistry and surface reactions on oxygen permeation and fuel conversion in ITM reactors

S. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives

Influence of a Dispersed Ignition in the Explosion of Two-Phase Mixtures

Analytical and Numerical Analysis of Micro Combustor for Gas Turbine Engine

AE 205 Materials and Energy Balances Asst. Prof. Dr. Tippabust Eksangsri. Chapter 4 Stoichiometry and MB with Reactions

TABLE OF CONTENT. Chapter 4 Multiple Reaction Systems 61 Parallel Reactions 61 Quantitative Treatment of Product Distribution 63 Series Reactions 65

Inhomogeneous Mixing Behavior of Recirculated Exhaust Gas in a Lean Premixed Flame

Effects of Ethane Partial Pressure on the Apparent Rate Expressions of Oxidative Couplin... Page 1 of 16

Thermochemistry X.S. Bai Thermochemistry

SIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR)

Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions

Modeling of Packed Bed Reactors: Hydrogen Production By the Steam Reforming of Methane and Glycerol

Developing Carbon Tolerance Catalyst for Dry Methane Reforming

SECOND ENGINEER REG. III/2 APPLIED HEAT

TRITIUM RECOVERY FROM WASTE USING A PALLADIUM MEMBRANE REACTOR

Steam Reforming in a Packed Bed Membrane Reactor

Fischer-Tropsch Synthesis over Co/ɣ-Al 2 O 3 Catalyst: Activation by Synthesis Gas

CFD Flow and Heat Transfer Simulation for Empty and Packed Fixed Bed Reactor in Catalytic Cracking of Naphtha

Kinetics of palladium catalysed partial oxidation of methane

Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati

A Project for Thermodynamics II. Entitled

Chemical Kinetics of Ethane Oxidation and Methane Oxidation with Platinum

Chapter 15. In the preceding chapters we limited our consideration to CHEMICAL REACTIONS. Objectives

AP Chemistry A. Allan Chapter Six Notes - Thermochemistry

Chemistry: The Central Science. Chapter 5: Thermochemistry

11B, 11E Temperature and heat are related but not identical.

Effect of the Reaction Bath Temperature in a Fixed-bed Reactor for Oxidation of o-xylene over V 2 O 5 /TiO 2 Catalysts

Modeling heterogeneous and homogeneous reactions in the high-temperature catalytic combustion of methane

By Rogéria Amaral and Sébastien Thomas

Experimental study of the combustion properties of methane/hydrogen mixtures Gersen, Sander

REPORT ON WORKSHOP: DEFINING RTE SOLUTION METHODS FOR MULTI-MODE TRANSFER. June 3-4, 2016, Istanbul Jack Howell and Pinar Mengüç

Numerical Modeling of Anode Baking Process with COMSOL Multiphysics

O 2 CH 4 CO 2 H

EFFECT OF CARBON DIOXIDE, ARGON AND HYDROCARBON FUELS ON THE STABILITY OF HYDROGEN JET FLAMES

NAME INDEX NUMBER CANDIDATE SIGNATURE DATE.. 233/1/ CHEMISTRY PAPER 1 (THEORY) JULY/AUGUST HOURS

HEAT TRANSFER 1 INTRODUCTION AND BASIC CONCEPTS 5 2 CONDUCTION

Vacuum Gas Carburizing with Acetylene - Gas Phase Modeling of a Bench Scale Reactor

Gas Chromatography. Presented By Mr. Venkateswarlu Mpharm KTPC

APPENDIX A: LAMINAR AND TURBULENT FLAME PROPAGATION IN HYDROGEN AIR STEAM MIXTURES*

Simple Kinetic Modeling of Selective Reduction of Nitric Oxide in Diesel Exhaust Over Cu-Zn/ZSM-5 Monolithic Catalyst

Transcription:

Research Article Journal of Energy Management and Technology (JEMT) Vol. 2, Issue 1 20 Experimental investigation of Methane Partial Oxidation for Hydrogen Production HAMID REZA LARI 1 AND MOHAMMAD REZA SHAHNAZARI 2 1,2 K.N.Toosi University of Technology, Tehran, Iran * Corresponding author:hlari@nri.ac.ir Manuscript received December 12, 2017; accepted January 18,2018. Paper no. JEMT-1712-1051 The thermal partial oxidation process of methane was investigated experimentally. Thermodynamic calculations and kinetic simulation were performed to determine the practical operating conditions. Experimentally, a porous material-based reactor was built to perform the partial oxidation process. Al2O3 granular with different size selected as porous structures and installed in the reaction zone. Physical parameters include diameter of the rector is considered for other experimental test. As the result the temperature profiles along the reactor central axis and concentration profiles of CO and H2 were measured. Also the effects of equivalence air ratio on the reforming process were investigated. Different equivalence air ratio 1.5 and 2 was considered to perform the partial oxidation process. 2018 Journal of Energy Management and Technology keywords: methane reforming, hydrogen, partial Oxidation, porous media, combustion. http://dx.doi.org/10.22109/jemt.2018.110940.1051 1. INTRODUCTION Fuel cells are considered to be the green power sources for the 21st century, and may make the hydrogen economy reality. Solid oxide fuel cell is a technology that H2 as well as CO can be utilized as fuel. Thus, synthesis gas (syngas), which can be produced from hydrocarbon fuels, can be directly utilized in such high temperature fuel cells. The syngas can be produced from natural gas, by steam reforming, catalytic or thermal partial oxidation (POX) and auto-thermal oxidation [1]. Compared with the other reforming technologies partial oxidation process has several advantageous. There is no need for external heat sources and additional feeds like water as in steam reforming. Moreover, the process is catalyst free and thus catalyst deactivation problems are completely out of consideration. It has a good dynamic response, a very simple system design and can be applied to almost all hydrocarbons [1]. The process of partial oxidation is considered in porous media rector. Gas mixture combustion in porous media is substantially different from free flame homogeneous oxidation. In the free flame propagation, the dominant mechanism for heat transfer is the convection of hot combustion products. Due to poor heat transfer properties of gas mixture, combustion occurs inside a thin layer of reaction. This leads to a poor upstream heat transfer. A solid porous material placed in the combustion media can transfer the total heat more by radiation and thermal conductivity resulting from the thermal properties of the solid matrix (thermal conductivity and thermal radiation) [2]. Therefore, the upstream flame enthalpy is partially transferred through radiation and conduction, improving the stability of oxidation. Fig. 1 displays the variations in temperature of gas and porous solid along the flame length. In this paper, the thermal partial oxidation process of methane in an inert porous reactor is studied experimentally. The objective of this paper is to find the best performance of reactor for optimum hydrogen production. In numerical model for comparing to experimental results the kinetic and thermodynamic data sets GRI Mech. 3.0 were used. In the experimental part, tabular reactor with different diameters, based on porous media, was built. In the diverging conical section the reaction front (oxidation zone) is stabilized. An insulation layer was used in order to reach near adiabatic conditions. Temperature profiles along the reactor central axis and concentration profiles of CO, H2, CO2, and CH4 were measured in three sections. Different sizes of granule for Al2O3 as porous structures were installed in the reaction zone. The effects of mas flow rate of fuel and equivalence ratio on the reforming process were investigated. 2. RELEVANT LITERATURE Numerous experimental and numerical studies have been conducted on hydrogen production in porous media. In discussing hydrogen production in porous media, the literature revolves around the effects of reactor design, structure of porous media material, equivalence ratio, flow rate and partial oxidation in porous media. Using the Péclet number, the incomplete heat

Research Article Journal of Energy Management and Technology (JEMT) Vol. 2, Issue 1 21 the effects of different parameters such as inlet speed (heat load), equivalence ratio and porous material (thermal conductivity and specific heat of porous material). In this research, the output products from the reactor were cooled by an exchanger and then directed as dry gas into the gas chromatograph or gas analyzer system [14]. The product compositions were analyzed for H2, CO, CO2, N2, O2, CH4 and C2H2. Zhdanok conducted several tests for different porous materials, and geometries of porous media [15]. Fig. 1. Variations in temperature of gas and porous solid along the flame length [3] combustion of liquid and gas hydrocarbon fuels took place in the inert base reformer of a porous medium [4]. Babkin investigated the incomplete combustion of methane in two porous media with different thermal conductivity and radiation properties, including silicon carbide foam and aluminum oxide AL2O3 [5]. In an experiment using the stability parameter of Peclet number, Pedersen et al. examined the rich combustion of methanol, methane, octane and gasoline inside an inert porous medium in a quartz or porcelain tube where hydrogen production could possibly take place [6,7]. Furthermore, they employed various types of porous media with different properties including aluminum oxide foam and aluminum oxide grains. This experiment also used solid granules of aluminum oxide and silicon carbide foam in different sizes. Hsu et al. found that the heat transfer rate and radiation coefficients significantly affect the prediction of gas and solid phase temperature in porous media [8]. To achieve correct results in the numerical simulation, it is essential to apply the porous media properties based on accurate empirical data. Howell et al. characterized porous media with several properties such as effective conductivity, permeability coefficients and vanishing coefficient based on empirical tests [9]. Miguel et al. investigated numerically methane thermal partial oxidation (TPOX) within a small scale inert porous media based reactor in order to explore the operating conditions and possible procedures for maximizing the reforming efficiency and minimizing the soot formation [10]. Han et al. optimized detailed reaction mechanism for methane combustion for methane partial oxidation, based on experimental results from a flow reactor. So that the optimized mechanism can predict CO and H2 production more accurately in methane partial oxidation [11]. Using several experiments and numerical results, Pan et al. obtained the effective conductivity of porous media for a few types of ceramics [12]. Malico and Pereira explored the effect of radiation properties of porous media on the function of cylindrical porous torches in a two-dimensional scenario [13]. They found that temperature distribution is strongly dependent on radiation parameters, particularly scattering phase function. Regardless of the radiation, the results will be in good agreement with the available experimental data. In other empirical studies on incomplete combustion in porous media, Al-Hamamre et al. examined 3. POROUS MEDIA COMBUSTION It is essential to provide a few definitions prior to exploring the combustion in a porous medium. A porous medium is composed of a heterogeneous system made of solid matrix with empty space filled by fluids. It contains a number of pores between a few phase particles inside a duct or control volume. Porous media are characterized by certain parameters constituting the properties of porous material. Oxidation in porous media can be classified into two stationary and non-stationary. The former can be achieved through stabilizing the flame inside the porous media at a specific zone. On the other hand, the combustion wave in the non-stationary mode travels alternatively forward and backward in a porous solid matrix. 4. FUEL-AIR RATIO AND EQUIVALENCE RATION The mass ratio of fuel to air is called fuel-air ratio represented by FA. FA = m f m a (1) m f and m a are the mass of fuel and mass of air, respectively. Air-fuel ratio plays a crucial role in the quality of combustion. Combustion will not occur if the air-fuel ratio is too large or too small. In between the two boundaries, however, this ratio should be selected depending on the condition until the highest efficiency is achieved. 5. EXPERIMENTAL APPARATUS Experiments on partial oxidation were conducted using the setup schematically is shown in Fig. 2. The apparatus consisted of a combustion tube filled with a porous medium, fuel and air supply system, temperature measurement system. For experimental tests, combustion tube with three internal diameters (Fig. 3), wall thickness of 2 mm were made of quartz. Also for investigating the effect of the length of the tubes were made in three size of length. The inner surface of the combustion tube was covered with a 2mm layer of ceramic paper insulation. Also to prevent heat losses and achieve quasi-uniform temperature profiles, additional 20mm thick high-temperature insulation was applied on the external diameter of the reactor. The packed bed of 1, 3 and 5 mm solid Al2O3 spheres were used as a porous medium having a different porosity (Fig. 4). The combustible mixtures of methane with air were prepared by a continuous flow method where the fuel and air flows were metered using a set of mass flow controllers. To ensure uniform gas composition the reactants were premixed in a mixing chamber. The outlet gas of the reactor after passing through the gas cooling system was gathered into a reservoir tank to make a possible for testing sample. In the exhaust, the online oxygen, CO, CO2 and methane analyzers were used with measurement instruments. Fig. 5 shows the test set up for experimental investigation.

Research Article Journal of Energy Management and Technology (JEMT) Vol. 2, Issue 1 22 Fig. 5. Experimental set up for partial oxidation of methane filled with porous material. With an emphasis on higher hydrogen production, this study focused on the effects of porous media, including structural properties of materials such as particle diameter as well as physical changes in the reactor such as diameter of the reactor, and finally the equivalence ratio of fuel and air, and inlet flow rate of fuel and oxygen. A. Temperature distribution Fig. 2. Schematic of the Set up In Fig.6 the modeling and experimental result of the temperature distribution inside of rector is shown. The model shows the results in different section. In experimental test as explained 3 thermometer were installed for measuring the temperature at different length of reactor. It should be mentioned that the test results are for the flow rate for methane and Oxygen at 100 lit/hour. The equivalence ratio in this model is fixed ( ϕ = 2). Fig. 3. Quartz tubes as combustion reactor with different length and diameter Fig. 6. Variation of temperature distribution in rector, Qmethane = 100lit/hour, QOxygen = 100lit/hour ϕ = 2 B. Effect of the mass flow rate Fig. 4. Packed bed of 1, 3 and 5 mm solid Al2O3 6. RESULTS AND DISCUSSION As noted earlier, this paper attempted to make an experimental test set for partial oxidation of methane for a tubular reactor From various aspects, the mass flow rate affects the thermal incomplete combustion in porous media: stationary time, temperature distribution and heat load. It has been found that high mass flow rate results in a shorter stationary time in the reactor. This implies that the reactant has shorter time to complete the reaction. Furthermore, the high mass flow rate of fuel/air was coupled with a high degree of distribution, covering a heat load higher than low mass flows. However, if the reactor is built for a wide range of fluctuations, such as porous-medium reactors, the effect of stationary time on the degree of conversion will be smaller than the effect of turbulence and heat released in porous

Research Article Journal of Energy Management and Technology (JEMT) Vol. 2, Issue 1 23 structure. It is therefore expected that hydrogen concentrations in high thermal loads be greater than those in low thermal loads, despite shorter stationary time. This is due to the fact that relative heat loss occurs lower in high thermal loads than lower ones. This implies that a small portion of the heat is transferred to preheat the incoming gases toward the upstream flow. Consequently, in the low power scenario, the hydrogen product is low and low temperature is expected in the reactor. Fig.7 displays the model and experimental results of temperature variations at the reactor outlet. Fig.8 displays the mole percentage of hydrogen, while Fig.9 displays the mole percentage of carbon monoxide in the reactor output for different flow rates. The diameter of the reactor is 30 mm, the length is 30 cm and equivalence ratio of 2 for different flow rates. As the figures show at flow rate 50 lit/h we have high difference between the model and the experimental test. It is because the stability of the combustion cannot be controlled and stabilized. Fig. 9. CO mole fraction variation with fuel flow rate Equivalence ratio 2, D=30mm & L= 30cm Fig. 7. Temperature variation with fuel flow rate Equivalence ratio 2, D=30mm & L= 30cm In this modeling, the particle diameter was considered to be between 1 mm and 5 mm given the thermo-physical properties of aluminum oxide. The experimental test was done for 1, 3, and 5 mm of granules. Fig. 10 provides the temperature variations in the reactor outlet with aluminum oxide granules at different diameters in a reactor with a diameter of 30 mm, a length of 30 cm and equivalence ratio of 1.5 and 2. In Fig. 11 also displays the changes in the mole percentage of hydrogen along with changes in granule diameter. As can be observed, the lower granule diameter leads to greater mole percentage of hydrogen. This is because of the surface area of porous structure affecting the fuel conversion in the partial oxidation process. Lower diameters of aluminum oxide granules create larger surface area and a more tortuous paths for the gas flow inside the porous matrix and intensifying the prorogation. This in turn leads to greater heat transfer between the gas and solid porous matrix, which can subsequently achieve greater fuel conversion. Fig. 12 displays the mole percentage distribution of CO mole fraction produced in the reactor. Fig. 8. Hydrogen mole fraction variation with fuel flow rate Equivalence ratio 2, D=30mm & L= 30cm Fig. 10. Temperature variation with diameter of granule, D=30mm & L= 30cm C. Effect of material diameter in porous matrix As stated earlier, high rate of heat transfer in porous bodies is not only due to high thermal conductivity, but also due to numerous pores in porous bodies, which increases the heat transfer surface and subsequently the heat transfer. One parameter effective in this process is the diameter of granules in the porous medium. 7. CONCLUSION Partial oxidation of methane for hydrogen production in an inert porous media is studied numerically and experimentally. A tabular pack bed rector considered and methane oxygen mixtures is used for feeding in the reactor. Al2O3 as inert porous

Research Article Journal of Energy Management and Technology (JEMT) Vol. 2, Issue 1 24 5. Babkin V., Filtrational combustion of gases. Present state of affairs and prospects. Pure Appl Chem,44:335-344(1993) 6. Pedersen-Mjaanes H., Chan L., Mastorakos E., Hydrogen production from rich combustion in porous media. International Journal of Hydrogen Energy,30(6):579 92(2005) 7. Pedersen-Mjaanes H., Mastorakos E., Optimization of hydrogen production from rich combustion of methane in porous media, Proceedings of the European Combustion Meeting, 2005. 8. Hsu P.F., Hoewll J.R., Mettews R.D., A Numerical Investigation of Premixed Combustion Within Porous Inert Media, ASME J. of Heat Transfer, 115(3):744-750(1993) Fig. 11. Hydrogen mole fraction variation with diameter of granule, D=30mm & L= 30cm 9. Howell,J.R., Hall, M.J. and Ellzey, J.L., Combustion within Porous Inert Medium, ASME HTD, Heat Transfer in Porous Media and Two-Phase Flow, 302:1-21(1995.) 10. Miguel A.A. Mendes*, Jose M.C. Pereira, Jose C.F. Pereira, Numerical study of methane TPOX within a small scale Inert Porous Media based reformer, International journal of hydrogen energy, 39(9): 4311-4321(2014) 11. Han S., Park J., Song S., Min Chun K., Experimental and numerical study of detailed reaction mechanism optimization for syngas (H2 + CO) production by non-catalytic partial oxidation of methane in a flow reactor, International journal of hydrogen energy, 35(16):8762-8771(2010) Fig. 12. CO mole fraction variation with diameter of granule, D=30mm & L= 30cm media is applied in this study. The various parameters that may effect on the performance of the rector for hydrogen production were investigated. Regarding to the different concept of running the test set, the maximum difference between the models and experimental results are about 10% that shows a good match between the model and experimental tests. REFERENCES 1. Al-Hamamrea, Z., Voßb, S., Trimis, D, Hydrogen production by thermal partial oxidation of hydrocarbon fuels in porous media based reformer, International Journal of hydrogen energy, 34:827-835(2009) 2. Trimis D., Durst F., Combustion in a porous medium advances and applications, Combust Sci Technol, 68:121-153(1996) 12. Pan H.L., Pickenäcker O., Pickenäcker K., Trimis D., Weber T., Experimental Determination of Effective Heat Conductivities of Highly Porous Media, 5th European Conference on Industrial Furnaces and Boilers, Porto, 11-14(2000) 13. Malico I., Pereira J.C.F., Numerical Study on the Influence of Radiative Properties in Porous Media Combustion, ASME Journal of Heat Transfer, 123(5):951-957(2001) 14. Al-Hamamre Z., Voss S., Trimis D., Hydrogen production by thermal partial oxidation of hydrocarbon fuels in porous media based reformer. International Journal of Hydrogen Energy, 34(2):827 832(2009) 15. Zhdanok S.A., Porous media combustion based hydrogen production, European Combustion Meeting, 2003. 16. Babkin V.S., Propagation of premixed gaseous explosion flame in porous media. Combust Flame, 87(2):182 190(1991) 17. Takeno T., Sato K., An excess enthalpy theory. Combust Sci Technology,20(1-2):73 84(1979) 3. Abdul Mujeebu, M., Abdullah, M.Z., Abu Bakar, M.Z., Mohamad A.A., Muhad, R.M.N., Abdullah, M.K., Combustion in porous media and its applications A comprehensive survey, Journal of Environmental Management, 90: 2287 2312(2009) 4. Al-Hamamre Z., Trimis D., Wawrzinek K., Hydrogen Production by Thermal partial oxidation of methane in porous burners, International Journal of Hydrogen Energy, 35(11): 5367 5377(2010)