Optimization of Reduced Kinetic Models for Reactive Flow Simulations

Author(s):  
P. Gokulakrishnan ◽  
R. Joklik ◽  
D. Viehe ◽  
A. Trettel ◽  
E. Gonzalez-Juez ◽  
...  

A robust optimization scheme, known as rkmGen, for reaction rate parameter estimation has been developed for the generation of reduced kinetics models of practical interest for reactive flow simulations. It employs a stochastic optimization algorithm known as simulated annealing (SA), and is implemented in C++ and coupled with Cantera, a chemical kinetics software package, to automate the reduced kinetic mechanism generation process. Reaction rate parameters in reduced order models can be estimated by optimizing against target data generated from a detailed model or by experiment. Target data may be of several different kinds: ignition delay time, blow-out time, laminar flame speed, species time-history profiles, and species reactivity profiles. The software allows for simultaneous optimization against multiple target data sets over a wide range of temperatures, pressures, and equivalence ratios. In this paper, a detailed description of the optimization strategy used for the reaction parameter estimation is provided. To illustrate the performance of the software for reduced kinetic mechanism development, a number of test cases for various fuels were used: one-step, three-step, and four-step global reduced kinetic models for ethylene, Jet-A and methane, respectively, and a 50 step semiglobal reduced kinetic model for methane. The 50 step semiglobal reduced kinetic model was implemented in the Star*CCM+ commercial CFD code to simulate Sandia Flame D using laminar flamelet libraries and compared with the experimental data. Simulations were also performed with the GRI3.0 mechanism for comparisons.

Author(s):  
P. Gokulakrishnan ◽  
R. Joklik ◽  
D. Viehe ◽  
A. Trettel ◽  
E. Gonzalez-Juez ◽  
...  

A robust optimization scheme, known as rkmGen, for reaction rate parameter estimation has been developed for the generation of reduced kinetics models of practical interest for reactive flow simulations. It employs a stochastic optimization algorithm known as Simulated Annealing, and is implemented in C++ and coupled with Cantera, a chemical kinetics software package, to automate the reduced kinetic mechanism generation process. Reaction rate parameters in reduced order models can be estimated by optimizing against target data generated from a detailed model or by experiment. Target data may be of several different kinds: ignition delay time, blow-out time, laminar flame speed, species time-history profiles and species reactivity profiles. The software allows for simultaneous optimization against multiple target data sets over a wide range of temperatures, pressures and equivalence ratios. In this paper, a detailed description of the optimization strategy used for the reaction parameter estimation is provided. To illustrate the performance of the software for reduced kinetic development, a number of test cases for various fuels were used: one-step, three-step and four-step global reduced kinetic models for ethylene, Jet-A and methane, respectively, and a fifty-step semi-global reduced kinetic model for methane. The fifty-step semi-global reduced kinetic model was implemented in the Star*CCM+ commercial CFD code to simulate Sandia Flame D using laminar flamelet libraries and compared with the experimental data. Simulations were also performed with the GRI3.0 mechanism for comparisons.


Author(s):  
P. Gokulakrishnan ◽  
S. Kwon ◽  
A. J. Hamer ◽  
M. S. Klassen ◽  
R. J. Roby

The reduced kinetic mechanism for syngas/methane developed in the present work consists of a global reaction step for fuel decomposition in which the fuel molecule breaks down into CH2O and H2. A detailed CH2O/H2/O2 elementary reaction sub-set is included as the formation of intermediate combustion radicals such as OH, H, O, HO2, and H2O2 is essential for accurate predictions of non-equilibrium phenomena such as ignition and extinction. Since the chemical kinetics of H2 and CH2O are the fundamental building blocks of any hydrocarbon oxidation, the inclusion of detailed kinetic mechanisms for CH2O and H2 oxidation enables the reduced mechanism to predict over a wide range of operating conditions provided the reaction rate parameters of fuel-decomposition reaction is optimized over those conditions. Therefore, the rate coefficients for the fuel-decomposition step are estimated and optimized for the ignition delay time measurements of CH4, H2, CH4/H2, CH4/CO and CO/H2 mixtures available in the literature over a wide range of pressures, temperatures and equivalence ratios that are relevant to gas turbine operating conditions. The optimized reduced mechanism, consisting of 15 species and around 40 reactions, is able to predict the ignition delay time and laminar flame speed measurements of CH4, H2, CH4/H2, CH4/CO and CO/H2 mixtures fairly well over a wide range conditions. The model predictions are also compared with that of GRI3.0 mechanism. The reduced kinetic mechanism predicts the ignition delay time of CH4 and CH4/H2 mixtures far better than GRI mechanism at higher pressures. To demonstrate the predictive capability of the model in reactive flow systems, the reduced mechanism was implemented in Star-CD/KINetics commercial code using a RANS turbulence model to simulate CH4/air premixed combustion in a backward facing step. The CFD model predictions of the stable species in the exhaust gas agree well with the GRI mechanism predictions in a chemical reactor network modeling by approximating the backward facing step with a series of perfectly-stirred reactor and plug-flow reactor.


Entropy ◽  
2021 ◽  
Vol 23 (4) ◽  
pp. 387
Author(s):  
Yiting Liang ◽  
Yuanhua Zhang ◽  
Yonggang Li

A mechanistic kinetic model of cobalt–hydrogen electrochemical competition for the cobalt removal process in zinc hydrometallurgical was proposed. In addition, to overcome the parameter estimation difficulties arising from the model nonlinearities and the lack of information on the possible value ranges of parameters to be estimated, a constrained guided parameter estimation scheme was derived based on model equations and experimental data. The proposed model and the parameter estimation scheme have two advantages: (i) The model reflected for the first time the mechanism of the electrochemical competition between cobalt and hydrogen ions in the process of cobalt removal in zinc hydrometallurgy; (ii) The proposed constrained parameter estimation scheme did not depend on the information of the possible value ranges of parameters to be estimated; (iii) the constraint conditions provided in that scheme directly linked the experimental phenomenon metrics to the model parameters thereby providing deeper insights into the model parameters for model users. Numerical experiments showed that the proposed constrained parameter estimation algorithm significantly improved the estimation efficiency. Meanwhile, the proposed cobalt–hydrogen electrochemical competition model allowed for accurate simulation of the impact of hydrogen ions on cobalt removal rate as well as simulation of the trend of hydrogen ion concentration, which would be helpful for the actual cobalt removal process in zinc hydrometallurgy.


1989 ◽  
Vol 67 (5) ◽  
pp. 857-861 ◽  
Author(s):  
Shin-Ichi Miyamoto ◽  
Tetsuo Sakka ◽  
Matae Iwasaki

The reaction rate of hydrogen isotope exchange between D2 and H2O catalyzed by platinum plate is studied. The exchange reaction is described with the kinetic model which is the modification of that for the exchange reaction catalyzed by alumina-supported platinum catalyst. For the comparison of experimental results with this model relative amount of the number of sites for hydrogen adsorption was estimated from the initial rate of hydrogen isotope exchange between H2 and D2 on the same surface. The results show that the kinetic model is applicable for the plate catalyst if the number of the sites for hydrogen absorption, which is very sensitive to the surface state of the catalyst, was estimated not from the macroscopic surface area but from our scheme. Keywords: hydrogen isotope exchange reaction, platinum plate as catalyst.


The Analyst ◽  
2015 ◽  
Vol 140 (9) ◽  
pp. 3121-3135
Author(s):  
Fereshteh Emami ◽  
Marcel Maeder ◽  
Hamid Abdollahi

Schematic of intertwined equilibrium-kinetic model at time = 0,1,2…T when both equilibrium and kinetic models are solved explicitly.


2018 ◽  
Vol 32 (1) ◽  
pp. 867-874 ◽  
Author(s):  
Mario Díaz-González ◽  
Cesar Treviño ◽  
Juan C. Prince

Author(s):  
P. Gokulakrishnan ◽  
M. S. Klassen ◽  
R. J. Roby

Ignition delay times of a “real” synthetic jet fuel (S8) were measured using an atmospheric pressure flow reactor facility. Experiments were performed between 900 K and 1200 K at equivalence ratios from 0.5 to 1.5. Ignition delay time measurements were also performed with JP8 fuel for comparison. Liquid fuel was prevaporized to gaseous form in a preheated nitrogen environment before mixing with air in the premixing section, located at the entrance to the test section of the flow reactor. The experimental data show shorter ignition delay times for S8 fuel than for JP8 due to the absence of aromatic components in S8 fuel. However, the ignition delay time measurements indicate higher overall activation energy for S8 fuel than for JP8. A detailed surrogate kinetic model for S8 was developed by validating against the ignition delay times obtained in the present work. The chemical composition of S8 used in the experiments consisted of 99.7 vol% paraffins of which approximately 80 vol% was iso-paraffins and 20% n-paraffins. The detailed kinetic mechanism developed in the current work included n-decane and iso-octane as the surrogate components to model ignition characteristics of synthetic jet fuels. The detailed surrogate kinetic model has approximately 700 species and 2000 reactions. This kinetic mechanism represents a five-component surrogate mixture to model generic kerosene-type jets fuels, namely, n-decane (for n-paraffins), iso-octane (for iso-paraffins), n-propylcyclohexane (for naphthenes), n-propylbenzene (for aromatics) and decene (for olefins). The sensitivity of iso-paraffins on jet fuel ignition delay times was investigated using the detailed kinetic model. The amount of iso-paraffins present in the jet fuel has little effect on the ignition delay times in the high temperature oxidation regime. However, the presence of iso-paraffins in synthetic jet fuels can increase the ignition delay times by two orders of magnitude in the negative temperature (NTC) region between 700 K and 900 K, typical gas turbine conditions. This feature can have a favorable impact on preventing flashback caused by the premature autoignition of liquid fuels in lean premixed prevaporized (LPP) combustion systems.


Sign in / Sign up

Export Citation Format

Share Document