Constrained-Equilibrium Modeling of Methane Oxidation in Air

2013 ◽  
Author(s):  
Ghassan Nicolas ◽  
Mohammad Janbozorgi ◽  
Hameed Metghalchi

The Rate-Controlled Constrained-Equilibrium (RCCE) has been further developed and applied to model methane/air combustion process. The RCCE method is based on local maximization of entropy or minimization of a relevant free energy at any time during the non-equilibrium evolution of the system subject to a set of constraints. The constraints are imposed by slow rate-limiting reactions. Direct integration of the rate equations for the constraint potentials has been employed. Once the values of the potentials are obtained, the concentration of all species can be calculated. A set of constraints has been developed for methane/air mixtures in the method of Rate-Controlled Constrained-Equilibrium (RCCE). The model predicts the ignition delay times, which have been compared to those predicted by detailed kinetic model (DKM) and with shock tube experimental measurements. The DKM includes 60 H/O/C1–2/N species and 352 reactions. The RCCE model using 16 constraints has been applied for combustion modeling in a wide range of initial temperatures (900–1200 K), pressures (1–50 atmospheres) and fuel-air equivalence ratio (0.6–1.2). The predicted results of using RCCE are within 5% of those of DKM model and are in excellent agreement with experimental measurements in shock tubes.

2014 ◽  
Vol 136 (3) ◽  
Author(s):  
Ghassan Nicolas ◽  
Mohammad Janbozorgi ◽  
Hameed Metghalchi

Rate-controlled constrained-equilibrium method has been further developed to model methane/air combustion. A set of constraints has been identified to predict the nonequilibrium evolution of the combustion process. The set predicts the ignition delay times of the corresponding detailed kinetic model to within 10% of accuracy over a wide range of initial temperatures (900 K–1200 K), initial pressures (1 atm–50 atm) and equivalence ratios (0.6–1.2). It also predicts the experimental shock tube ignition delay times favorably well. Direct integration of the rate equations for the constraint potentials has been employed. Once the values of the potentials are obtained, the concentration of all species can be calculated. The underlying detailed kinetic model involves 352 reactions among 60 H/O/N/C1-2 species, hence 60 rate equations, while the RCCE calculations involve 16 total constraints, thus 16 total rate equations. Nonetheless, the constrained-equilibrium concentrations of all 60 species are calculated at any time step subject to the 16 constraints.


2015 ◽  
Vol 137 (6) ◽  
Author(s):  
Ghassan Nicolas ◽  
Hameed Metghalchi

The rate-controlled constrained-equilibrium (RCCE) method is a reduction technique based on local maximization of entropy or minimization of a relevant free energy at any time during the nonequilibrium evolution of the system subject to a set of kinetic constraints. In this paper, RCCE has been used to predict ignition delay times of low temperatures methane/air mixtures in shock tube. A new thermodynamic model along with RCCE kinetics has been developed to model thermodynamic states of the mixture in the shock tube. Results are in excellent agreement with experimental measurements.


Author(s):  
Sergio Ugarte ◽  
Mohamad Metghalchi ◽  
James C. Keck

Modeling of a non-equilibrium combustion process involves the solution of large systems of differential equations with as many equations as species present during the process. The process of chemical reaction and combustion is complicated since it may be governed by hundreds, sometimes thousands of microscopic rate processes. Integration of these equations simultaneously becomes more difficult with the complexity of the combustible. In order to reduce the size of these systems of equations, the Rate-Controlled Constrained-Equilibrium method (RCCE) has been proposed to model non-equilibrium combustion processes. This method is based on the Second Law of Thermodynamics, assuming that the evolution of a complex system can be described by a small number of rate-controlling reactions which impose slowly changing constraints on all allowed states of the system, therefore a non-equilibrium system will relax to its final equilibrium state through a sequence of rate controlled constrained equilibrium states. Oxidation induction times and concentration of species during a combustion process are found in a less complicated way with this method, as equations for constraints rather than for species determine the composition and evolution of the system. The time evolution of the system can be reduced since the number of constraints is much smaller than the number of species presents, so the number of equations to solve. The RCCE method has been applied to the stoichiometric combustion of mono-carbon fuels using 29 chemical species and 139 chemical reactions at different sets of pressure and temperature, ranging from 1 atm to 100 atm, and from 900 K to 1600 K respectively. Results of using 8, 9, 10 and 11 constraints compared very well to those of the detailed calculations at all conditions for the cases of formaldehyde (H2CO), methanol (CH3OH) and methane (CH4). For these systems, ignition delay times and major species concentrations were within 5% of the values given by detailed calculations, and computational saving times up to 50% have been met.


Author(s):  
Fatemeh Hadi ◽  
Mohammad Janbozorgi ◽  
M. Reza H. Sheikhi ◽  
Hameed Metghalchi

AbstractThe rate-controlled constrained-equilibrium (RCCE) method is employed to study the interactions between mixing and chemical reaction. Considering that mixing can influence the RCCE state, the key objective is to assess the accuracy and numerical performance of the method in simulations involving both reaction and mixing. The RCCE formulation includes rate equations for constraint potentials, density and temperature, which allows taking account of mixing alongside chemical reaction without splitting. The RCCE is a dimension reduction method for chemical kinetics based on thermodynamics laws. It describes the time evolution of reacting systems using a series of constrained-equilibrium states determined by RCCE constraints. The full chemical composition at each state is obtained by maximizing the entropy subject to the instantaneous values of the constraints. The RCCE is applied to a spatially homogeneous constant pressure partially stirred reactor (PaSR) involving methane combustion in oxygen. Simulations are carried out over a wide range of initial temperatures and equivalence ratios. The chemical kinetics, comprised of 29 species and 133 reaction steps, is represented by 12 RCCE constraints. The RCCE predictions are compared with those obtained by direct integration of the same kinetics, termed detailed kinetics model (DKM). The RCCE shows accurate prediction of combustion in PaSR with different mixing intensities. The method also demonstrates reduced numerical stiffness and overall computational cost compared to DKM.


Author(s):  
Sergio Ugarte ◽  
Mohamad Metghalchi ◽  
James C. Keck

Methanol oxidation has been modeled using the Rate-Controlled Constrained-Equilibrium method (RCCE). In this method, composition of the system is determined by constraints rather than by species. Since the number of constraints can be much smaller than the number of species present, the number of rate equations required to describe the time evolution of the system can be considerably reduced. In the present paper, C1 chemistry with 29 species and 140 reactions has been used to investigate the oxidation of stoichiometric methanol/oxygen mixture at constant energy and volume. Three fixed elemental constraints: elemental carbon, elemental oxygen and elemental hydrogen and from one to nine variable constraints: moles of fuel, total number of moles, moles of free oxygen, moles of free oxygen, moles of free valence, moles of fuel radical, moles of formaldehyde H2CO, moles of HCO, moles of CO and moles of CH3O were used. The four to twelve rate equations for the constraint potentials (LaGrange multipliers conjugate to the constraints) were integrated for a wide range of initial temperatures and pressures. As expected, the RCCE calculations gave correct equilibrium values in all cases. Only 8 constraints were required to give reasonable agreement with detailed calculations. Results of using 9 constraints showed compared very well to those of the detailed calculations at all conditions. For this system, ignition delay times and major species concentrations were within 0.5% to 5% of the values given by detailed calculations. Adding up to 12 constraints improved the accuracy of the minor species mole fractions at early times, but only had a little effect on the ignition delay times. RCCE calculations reduced the time required for input and output of data in 25% and 10% when using 8 and 9 constraints respectively. In addition, RCCE calculations gave valuable insight into the important reaction paths and rate-limiting reactions involved in methanol oxidation.


2018 ◽  
Vol 141 (2) ◽  
Author(s):  
Guangying Yu ◽  
Hameed Metghalchi ◽  
Omid Askari ◽  
Ziyu Wang

The rate-controlled constrained-equilibrium (RCCE), a model order reduction method, has been further developed to simulate the combustion of propane/oxygen mixture diluted with nitrogen or argon. The RCCE method assumes that the nonequilibrium states of a system can be described by a sequence of constrained-equilibrium states subject to a small number of constraints. The developed new RCCE approach is applied to the oxidation of propane in a constant volume, constant internal energy system over a wide range of initial temperatures and pressures. The USC-Mech II (109 species and 781 reactions, without nitrogen chemistry) is chosen as chemical kinetic mechanism for propane oxidation for both detailed kinetic model (DKM) and RCCE method. The derivation for constraints of propane/oxygen mixture starts from the eight universal constraints for carbon-fuel oxidation. The universal constraints are the elements (C, H, O), number of moles, free valence, free oxygen, fuel, and fuel radicals. The full set of constraints contains eight universal constraints and seven additional constraints. The results of RCCE method are compared with the results of DKM to verify the effectiveness of constraints and the efficiency of RCCE. The RCCE results show good agreement with DKM results under different initial temperature and pressures, and RCCE also reduces at least 60% CPU time. Further validation is made by comparing the experimental data; RCCE shows good agreement with shock tube experimental data.


2015 ◽  
Vol 138 (2) ◽  
Author(s):  
Ghassan Nicolas ◽  
Hameed Metghalchi

The rate-controlled constrained-equilibrium method (RCCE) has been further developed to model the combustion process of ethanol air mixtures. The RCCE is a reduction technique based on local maximization of entropy or minimization of a relevant free energy at any time during the nonequilibrium evolution of the system subject to a set of kinetic constraints. An important part of RCCE calculation is determination of a set of constraints that can guide the nonequilibrium mixture to the final stable equilibrium state. In this study, 16 constraints have been developed to model the nonequilibrium ethanol combustion process. The method requires solution of 16 differential equations for the corresponding constraint potentials. Ignition delay calculations of ethanol oxidizer mixtures using RCCE have been compared to those of detailed chemical kinetics using 37 species and 235 reactions. Agreement between the two models is very good. In addition, ignition delay of C2H5OH/O2/Ar mixtures using RCCE has been compared with the experimental measurements in the shock tube and excellent agreement has been reached validating the RCCE calculation.


Author(s):  
Mohammad Janbozorgi ◽  
Yue Gao ◽  
Hameed Metghalchi ◽  
James C. Keck

The Rate-Controlled Constrained-Equilibrium method (RCCE) is a powerful technique for simplifying the treatment of chemical reactions in complex systems. The method is based on the assumption that slow chemical reactions impose constraints on the allowed composition of such systems. Since the number of constraints can be very much smaller than the number of species, the number of rate equations to be integrated can be considerably reduced. In the present work, a kinetic scheme with 55 species and 366 reactions has been used to investigate stoichiometric Ethanol-oxygen mixture in a constant energy constant volume chamber. The state of the system was determined by three fixed elemental constraints: elemental carbon, elemental oxygen and elemental hydrogen and six variable constraints: moles of fuel, moles of fuel radicals, total number of moles, moles of free valence, moles of free oxygen, moles of OH+O+H. The 9 rate equations for the constraint potentials (LaGrange multipliers associated with the constraints) were integrated over a range of 1400 K-1700 K for initial temperature and 1atm-10atm for initial pressure. The RCCE calculations were in good agreement with detailed calculations and were faster than detailed calculations, which required integration of 55 species rate equations.


2018 ◽  
Author(s):  
Fatemeh Hadi ◽  
Vreg Yousefian ◽  
Ehsan Sarfaraz ◽  
Gian Paolo Beretta

The Rate-Controlled Constrained-Equilibrium (RCCE) is a model reduction scheme for chemical kinetics. It describes the evolution of a complex chemical system with acceptable accuracy with a number of rate controlling constraints on the associated constrained-equilibrium states of the system, much lower than the number of species in the underlying Detailed Kinetic Model (DKM). Successful approximation of the constrained-equilibrium states requires accurate identification of the constraints. One promising procedure is the fully automatable Approximate Singular Value Decomposition of the Actual Degrees of Disequilibrium (ASVDADD) method that is capable of identifying the best constraints for a given range of thermodynamic conditions and a required level of approximation. ASVDADD is based on simple algebraic analysis of the results of the underlying DKM simulation and is focused on the behavior of the degrees of disequilibrium (DoD) of the individual chemical reactions. In this paper, we propose a method, as part of our work-in-progress efforts, that could expand the applicability of the derived constraints. This method involves running DKM calculations for a wider range of initial conditions, appending the results of all these cases one after the other after normalizing, and finally running the ASVDADD method to get a set of ‘universal’ constraints applicable within that range of conditions. The effectiveness and robustness of the derived constraints is examined in hydrogen/oxygen ignition delay simulations and the results are compared with those obtained from DKM. The proof-of-concept results demonstrate the potential of the method for finding ‘universal’ constraints.


Sign in / Sign up

Export Citation Format

Share Document