A Reformulation of Degree of Disequilibrium Analysis for Automatic Selection of Kinetic Constraints in the Rate-Controlled Constrained-Equilibrium Method

2021 ◽  
pp. 1-25
Author(s):  
Fatemeh Hadi ◽  
Shrabanti Roy ◽  
Omid Askari ◽  
Gian Paolo Beretta

Abstract 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 introduce an alternative ASVDADD algorithm. Unlike the original ASVDADD algorithm that require the direct computation of the DKM-derived DoDs and the identification of the set of linearly independent reactions, in the alternative algorithm, the components of the overall degree of disequilibrium vector can be computed directly by casting the DKM as an RCCE simulation considering a set of linearly independent constraints equalling the number of chemical species in size. The effectiveness and robustness of the derived constraints from the alternative procedure is examined in hydrogen/oxygen and methane/oxygen ignition delay simulations and the results are compared with those obtained from DKM.

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.


2019 ◽  
Author(s):  
Shrabanti Roy ◽  
Fatemeh Hadi ◽  
Omid Askari

Abstract In this study, the fully automatable Approximate Singular Value Decomposition of the Actual Degrees of Disequilibrium (ASVDADD) method is used for combustion modeling of ethanol. Due to the importance of ethanol as one of the most common type of biofuels, modeling its reaction kinetics and chemical composition evolution during combustion is necessary. The detailed kinetic mechanism (DKM) considered here is generated by authors using reaction mechanism generator (RMG) technique and it consists of 66 species and 1031 reactions. Tracking this number of species and chemical reactions in computational fluid dynamic (CFD) analysis of engineering problems is prohibitive. To alleviate this issue, Rate-Controlled Constrained-Equilibrium (RCCE) model reduction scheme for chemical kinetics is employed. It describes the evolution of a complex chemical system with acceptable accuracy with a number of rates controlling constraints on the associated constrained-equilibrium states of the system, much lower than the number of species in the underlying DKM. Successful approximation of the constrained equilibrium states requires accurate identification of the constraints. One promising procedure is the 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, ASVDADD is used to derive the RCCE constraints and ethanol combustion is modeled using both DKM and RCCE. Comparison of RCCE results with those of DKM shows the effectiveness of the ASVDADD derived constraints which demonstrates the potential of the RCCE method for combustion modeling of heavy and complex fuels.


Author(s):  
Luca Rivadossi ◽  
Gian Paolo Beretta

The Rate-Controlled Constrained-Equilibrium (RCCE) model reduction scheme for chemical kinetics provides acceptable accuracies in predicting hydrocarbon ignition delays by solving a smaller number of differential equations than the number of species in the underlying Detailed Kinetic Model (DKM). To yield good approximations, the method requires accurate identification of the rate controlling constraints. Until recently, a drawback of the RCCE scheme has been the absence of a fully automatable and systematic procedure capable of identifying the best constraints for a given range of thermodynamic conditions and a required level of approximation. A recent paper [1] has proposed a new methodology for such identification based on a simple algebraic analysis of the results of a preliminary simulation of the underlying DKM, focused on the behaviour of the degrees of disequilibrium (DoD) of the individual chemical reactions. The new methodology is based on computing an Approximate Singular Value Decomposition of the Actual Degrees of Disequilibrium (ASVDADD) obtained from the DKM simulation. The effectiveness and robustness of the method has been demonstrated in [1] for some cases of methane/oxygen ignition by considering a C1/H/O (29 species/133 reactions) sub-mechanism of the GRI-Mech 3.0 scheme and comparing the results of a DKM simulation with those of RCCE simulations based on increasing numbers of ASVDADD constraints. The RCCE results are in excellent agreement with DKM predictions for relatively small numbers of RCCE constraints. Here we provide a demonstration of the new method for some cases of shock-tube ignition of a natural gas/air mixture, with higher hydrocarbons approximately represented by propane according to the full (53 species/325 reactions) GRI-Mech 3.0 scheme.


2017 ◽  
Vol 140 (5) ◽  
Author(s):  
Luca Rivadossi ◽  
Gian Paolo Beretta

The rate-controlled constrained-equilibrium (RCCE) model reduction scheme for chemical kinetics provides acceptable accuracies in predicting hydrocarbon ignition delays by solving a smaller number of differential equations than the number of species in the underlying detailed kinetic model (DKM). To yield good approximations, the method requires accurate identification of the rate controlling constraints. Until recently, a drawback of the RCCE scheme has been the absence of a systematic procedure capable of identifying optimal constraints for a given range of thermodynamic conditions and a required level of approximation. A recent methodology has proposed for such identification an algorithm based on a simple algebraic analysis of the results of a preliminary simulation of the underlying DKM, focused on the degrees of disequilibrium (DoD) of the individual chemical reactions. It is based on computing an approximate singular value decomposition of the actual degrees of disequilibrium (ASVDADD) obtained from the DKM simulation. The effectiveness and robustness of the method have been demonstrated for methane/oxygen ignition by considering a C1/H/O (29 species/133 reactions) submechanism of the GRI-Mech 3.0 scheme and comparing the results of a DKM simulation with those of RCCE simulations based on increasing numbers of ASVDADD constraints. Here, we demonstrate the new method for shock-tube ignition of a natural gas/air mixture, with higher hydrocarbons approximately represented by propane according to the full (53 species/325 reactions) GRI-Mech 3.0 scheme including NOx formation.


2020 ◽  
Vol 142 (7) ◽  
Author(s):  
Shrabanti Roy ◽  
Omid Askari

Abstract Reduction of the detail chemical kinetic mechanism is important in solving complex combustion simulation. In this work, a model reduction scheme rate-controlled constrained-equilibrium (RCCE) is considered in predicting the oxidation of ethanol. A detail kinetic mechanism by Merinov from Lawrence Livermore National Laboratory (LLNL) is used in modeling this reduction technique. The RCCE method considers constrained equilibrium states which subjected to a lower number of constraints compared to the number of species. It then has to solve a smaller number of differential equations compared to the number of equations required in solving the detailed kinetic model (DKM). The accuracy of this solution depends on the selection of the constraint. A systematic procedure which will help in identifying the constraint at an optimal level of accuracy is an essential for RCCE modeling. A fully automated Approximate Singular Value Decomposition of the Actual Degrees of Disequilibrium (ASVDADD) method is used in this study to derive the constraint for RCCE simulation. ASVDADD uses an algorithm which follows the simple algebraic analysis on results of underlying DKM to find the degree of disequilibrium (DoD) of the individual chemical reactions. The number of constraints which will be used in RCCE simulation can be selected to reduce the number of equations required to solve. In the current work, this ASVDADD method is applied on ethanol oxidation to select the constraint for RCCE simulation. Both DKM and RCCE calculations on ethanol fuel are demonstrated to compare the result of temperature distribution and an ignition delay time for validating the method.


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.


2020 ◽  
Vol 45 (1) ◽  
pp. 59-79 ◽  
Author(s):  
Guangying Yu ◽  
Fatemeh Hadi ◽  
Ziyu Wang ◽  
Hameed Metghalchi

AbstractDeveloping an effective model for non-equilibrium states is of great importance for a variety of problems related to chemical synthesis and combustion. Rate-Controlled Constrained-Equilibrium (RCCE), a model order reduction method that originated from the second law of thermodynamics, assumes that the non-equilibrium states of a system can be described by a sequence of constrained-equilibrium states kinetically controlled by a relatively small number of constraints within acceptable accuracy. The full chemical composition at each constrained-equilibrium state is obtained by maximizing (or minimizing) the appropriate thermodynamic quantities, e. g., entropy (or Gibbs functions), subject to the instantaneous values of RCCE constraints. Regardless of the nature of the kinetic constraints, RCCE always guarantees a correct final equilibrium state. This paper reviews the fundamentals of the RCCE method, its constraints, constraint potential formulations, and major constraint selection techniques, as well as the application of the RCCE method to combustion of different fuels using a variety of combustion models. The RCCE method has been proven to be accurate and to reduce computational time in these simulations.


Author(s):  
John Ross ◽  
Igor Schreiber ◽  
Marcel O. Vlad

In a chemical system with many chemical species several questions can be asked: what species react with other species: in what temporal order: and with what results? These questions have been asked for over one hundred years about simple and complex chemical systems, and the answers constitute the macroscopic reaction mechanism. In Determination of Complex Reaction Mechanisms authors John Ross, Igor Schreiber, and Marcel Vlad present several systematic approaches for obtaining information on the causal connectivity of chemical species, on correlations of chemical species, on the reaction pathway, and on the reaction mechanism. Basic pulse theory is demonstrated and tested in an experiment on glycolysis. In a second approach, measurements on time series of concentrations are used to construct correlation functions and a theory is developed which shows that from these functions information may be inferred on the reaction pathway, the reaction mechanism, and the centers of control in that mechanism. A third approach is based on application of genetic algorithm methods to the study of the evolutionary development of a reaction mechanism, to the attainment given goals in a mechanism, and to the determination of a reaction mechanism and rate coefficients by comparison with experiment. Responses of non-linear systems to pulses or other perturbations are analyzed, and mechanisms of oscillatory reactions are presented in detail. The concluding chapters give an introduction to bioinformatics and statistical methods for determining reaction mechanisms.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Meku Maruyama ◽  
Riku Matsuura ◽  
Ryo Ohmura

AbstractHydrate-based gas separation technology is applicable to the CO2 capture and storage from synthesis gas mixture generated through gasification of fuel sources including biomass. This paper reports visual observations of crystal growth dynamics and crystal morphology of hydrate formed in the H2 + CO2 + tetrahydropyran (THP) + water system with a target for developing the hydrate-based CO2 separation process design. Experiments were conducted at a temperature range of 279.5–284.9 K under the pressure of 4.9–5.3 MPa. To simulate the synthesis gas, gas composition in the gas phase was maintained around H2:CO2 = 0.6:0.4 in mole fraction. Hydrate crystals were formed and extended along the THP/water interface. After the complete coverage of the interface to shape a polycrystalline shell, hydrate crystals continued to grow further into the bulk of liquid water. The individual crystals were identified as hexagonal, tetragonal and other polygonal-shaped formations. The crystal growth rate and the crystal size varied depending on thermodynamic conditions. Implications from the obtained results for the arrangement of operating conditions at the hydrate formation-, transportation-, and dissociation processes are discussed.


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.


Sign in / Sign up

Export Citation Format

Share Document