A criterion based on computational singular perturbation for the identification of quasi steady state species: A reduced mechanism for methane oxidation with NO chemistry

2008 ◽  
Vol 154 (4) ◽  
pp. 761-774 ◽  
Author(s):  
Tianfeng Lu ◽  
Chung K. Law
2013 ◽  
Vol 78 (8) ◽  
pp. 1177-1188
Author(s):  
Zuozhu Wu ◽  
Xinqi Qiao ◽  
Zhen Huang

A criterion based on the computational singular perturbation (CSP) method is proposed in order to determine the number of quasi-steady state (QSS) species. This criterion is employed for the reduction of a detailed chemical kinetics mechanism for the oxidation of dimethyl ether (DME), involving 55 species and 290 reactions, leading to a 20 steps reduced mechanism which involves 26 species. A software package, named I-CSP, was developed to make the reduction process algorithmic. Input to the I-CSP includes (i) the detailed mechanism, (ii) the numerical solution of the problem for a specific set of operating conditions, (iii) the number of quasi steady state (QSS) species. The resulting reduced mechanism was validated both in homogenous reactor, including auto-ignition and PSR, over a wide range of pressures and equivalence ratios, and in a one-dimensional, unstretched, premixed, laminar steady DME/Air flame. Comparison of the results calculated with the detailed and the reduced mechanisms shows excellent agreement in the case of homogenous reactor, but discrepancies can be observed in the case of the premixed laminar flame.


2019 ◽  
Author(s):  
Δημήτριος Πατσατζής

Η Διδακτορική Διατριβή αυτή πραγματεύεται την ανάλυση δύο βιολογικών συστημάτων πολλαπλών χρονοκλιμάκων με την Υπολογιστική μέθοδο Ιδιόμορφων Διαταραχών (Computational Singular Perturbation, CSP). H CSP είναι μία αλγοριθμική μέθοδος της ασυμπτωτικής ανάλυσης, η οποία βρίσκει εφαρμογή στο ευρύτερο πλαίσιο της Γεωμετρικής Θεωρίας Ιδιόμορφων Διαταραχών (Geometric Singular Perturbation Theory). Η μέθοδος CSP έχει τη δυνατότητα να παρέχει τα κατάλληλα απλοποιημένα μοντέλα, τα κριτήρια για την εγκυρότητα τους και έναν αριθμό διαγνωστικών εργαλείων για την απόκτηση της φυσικής κατανόησης των υπό μελέτη μηχανισμών. Δεδομένου ενός δύσκαμπτου συστήματος (πολλαπλών χρονοκλιμάκων) οποιασδήποτε πολυπλοκότητας, η CSP μπορεί α) να κατασκευάσει έναν αριθμό απλοποιημένων μοντέλων, συμπεριλαμβανομένων και αυτών που παράγονται από τις παραδοσιακές τεχνικές της Προσέγγισης Μερικούς Ισορροπίας (Partial Equilibrium Approximation, PEA) και της Προσέγγισης Οιωνεί ΜόνιμηςΚατάστασης (Quasi Steady-State Approximation, QSSA), β) να καθορίσει τις συνθήκες υπό τις οποίες τα απλοποιημένα αυτά μοντέλα παρέχουν πρώτης τάξης ακρίβεια και ευστάθεια, γ) να αναγνωρίσει τις “γρήγορες” χημικές αντιδράσεις και μεταβλητές, δ) να εντοπίσει τις αντιδράσεις που παράγουν τις “γρήγορες”και “αργές” χρονοκλίμακες, ε) να διακρίνει τις αντιδράσεις που συνεισφέρουν στη δημιουργία του Slow Invariant Manifold που αναπτύσσεται στο χώρο των φάσεων, πάνω στο οποίο εξελίσσεται το σύστημα και στ) να αναγνωρίσει τις αντιδράσεις που επηρεάζουν την αργή εξέλιξη όλων των χημικών ενώσεων: δηλαδή,όλες τις πληροφορίες που απαιτούνται για τη ρύθμιση του μοντέλου. Τα δύο βιολογικά συστήματα πολλαπλών χρονοκλιμάκων που αναλύονται είναι (i) ο μηχανισμός Michaelis-Menten και (ii) ένας μηχανισμός μεταβολισμού γλυκόζης και γαλακτικού οξέος για την παραγωγή ενέργειας στα αστροκύτταρα και τους νευρώνες του ανθρώπινου εγκεφάλου. Ο μηχανισμός πολλαπλών χρονοκλιμάκων Michaelis-Menten (MM) προσομοιώνει τη δράση των ενζύμων σε ενδοκυτταρικό και εξωκυττταρικό επίπεδο και χρησιμοποιείται σε μεγάλη ποικιλία επιστημονικών πεδίων. Το βασικό αποτέλεσμα είναι η εισαγωγή μέσω της CSP μίας νέας έκφρασης για το ρυθμό του ΜΜ μηχανισμού, η οποία είναι έγκυρη σε όλο τον παραμετρικό χώρο που χαρακτηρίζεται από δυναμική πολλαπλών χρονοκλιμάκων. Η προσπάθεια κατασκευής μιας τέτοιας έκφρασης ξεκίνησε με την εισαγωγή του ΜΜ μοντέλου το 1913 και, συνεπώς, υπάρχει εκτενής βιβλιογραφία για τέτοιου είδους εκφράσεις και για τις συνθήκες εγκυρότητάς των. Η πιο διαδεδομένη από αυτές τις εκφράσεις εισήχθηκε το 1925, είναι γνωστή ως ΜΜ εξίσωση (MM equation) και είναι έγκυρη σε περιορισμένη περιοχή του παραμετρικού χώρου. Όλες οι υπάρχουσες εκφράσεις έχουν περιορισμένη εγκυρότητα στον παραμετρικό χώρο και έχουν παραχθεί μέσω των παραδοσιακών PEA/QSSA προσεγγίσεων. Η νέα έκφραση, απλοποιείται σε όλες τις υπάρχουσες εκ φράσεις, στις περιοχές στις οποίες η κάθε μία εξ’ αυτών είναι έγκυρη. Οι απλοποιήσεις αυτές εισάγουν έγκυρα κριτήρια για την ακρίβεια των PEA/QSSA προσεγγίσεων. Αποδείχθηκε ότι η έκφραση που παρείχε η CSP είναι πάντα ευσταθής και η πιο ακριβής από όλες τις γνωστές PEA/QSSA εκφράσεις, γεγονός που το καθιστά ιδανική επιλογή για την ανάλυση των πειραματικών αποτελεσμάτων. Ο μηχανισμός μεταβολισμού της γλυκόζης και του γαλακτικού οξέος στον ανθρώπινο εγκέφαλο, ο οποίος εμπεριέχει ΜΜ εκφράσεις, αναλύθηκε με σκοπό τη φυσική κατανόηση δύο αντικρουόμενων υποθέσεων στη νευροχημεία του εγκεφάλου. Οι αντετιθέμενες υποθέσεις αφορούν το κατά πόσον το γαλακτικό οξύ ρέει από τα αστροκύτταρα στους νευρώνες ή από τους νευρώνες στα αστροκύτταρα. Σημειωτέον ότι υπάρχουν πειραματικά αποτελέσματα, τα οποία ενισχύουν και τις δύο αυτές υποθέσεις. Στη παρούσα ανάλυση, τα αλγοριθμικά εργαλεία της CSP εφαρμόζονται για την ανάλυση της γρήγορης και αργής δυναμικής του συστήματος, με σκοπό την διερεύνηση εγκυρότητας των δύο αυτών υποθέσεων. Αρχικά,αναγνωρίζονται οι ισορροπίες που αναδύονται μέσω της γρήγορης δυναμικής και αποκαλύπτονται τα μονοπάτια μέσω των οποίων μεταφέρεται η γλυκόζη και το γαλακτικό οξύ, για κάθε υπόθεση. Δεδομένων αυτών των πληροφοριών, η αργή δυναμική του συστήματος αναλύεται περαιτέρω, αναδεικνύοντας τις αντιδράσεις και τις μεταβλητές που οδηγούν την αργή δυναμική του συστήματος. Αναδεικνύεται ότι και οι δύο υποθέσεις είναι έγκυρες σε διαφορετικές όμως περιοχές στον παραμετρικό χώρο και καθορίζονται οι βιολογικές διεργασίες που εξασφαλίζουν την εγκυρότητα της κάθε υπόθεσης. Τέλος, αποδεικνύεται ότι οι βιολογικές διεργασίες.


2018 ◽  
Vol 141 (4) ◽  
Author(s):  
David Awakem ◽  
Marcel Obounou ◽  
Hermann Chopkap Noume

This work highlights the ability of the computational singular perturbation (CSP) method to calculate the significant indices of the modes on evolution of species and the degree of participation of reactions. The exploitation of these indices allows us to deduce the reduced models of detailed mechanisms having the same physicochemical properties. The mechanism used is 16 species and 41 reversible reactions. A reduction of these 41 reactions to 22 reactions is made. A constant pressure application of the detailed and reduced mechanism is made in OpenFOAM free and open source code. Following the Reynolds-averaged Navier–Stokes simulation scheme, standard k–ε and partial stirred reactor are used as turbulence and combustion models, respectively. To validate the reduced mechanism, comparison of numerical results (temperature and mass fractions of the species) was done between the detailed mechanism and the simplified model. This was done using the DVODE integrator in perfectly stirred reactor. After simulation in the computational fluid code dynamic (CFD) OpenFOAM, other comparisons were made. These comparisons were between the experimental data of a turbulent nonpremixed diffusion flame of type “DLR-A flame,” the reduced mechanism, and the detailed mechanism. The calculation time using the simplified model is considerably reduced compared to that using the detailed mechanism. An excellent agreement has been observed between these two mechanisms, indicating that the reduced mechanism can reproduce very well the same result as the detailed mechanism. The accordance with experimental results is also good.


2019 ◽  
Vol 142 (6) ◽  
Author(s):  
Cyrille Gnentedem ◽  
David Awakem ◽  
Marcel Obounou ◽  
Henri Paul Ekobena Fouda ◽  
Donatien Njomo

Abstract This work presents an application of a reduced chemical kinetic mechanism using computational singular perturbation (CSP) based on the significant indices of the modes on the evolution of species and the degree of participation of reactions. With this approach, the mechanism of Yang and Pope is reduced to 22 reversible reactions (RR22). In this study, the tabulation of ignition delays is made with Yang and Pope mechanism, GRI 3.0, and the reduced mechanism RR22; the results obtained show a good agreement among the three mechanisms. The “Modèle Intermittent Lagrangien” (MIL) necessary to calculate the chemical source term of the transport equation of the species requires the library of ignition delays determined above and a probability density function (PDF) of the mixture fraction presumed by a beta distribution. The scalar variance, one of the key parameters for the determination of the presumed beta function, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modeled using either an algebraic model or a transport equation. All these models are introduced in the computational fluid dynamics “Code-Saturne” to simulate a turbulent CH4/H2/N2 jet flame (DLR Flame A) performed at the Deutsches Zentrum für Luft-und Raumfahrt (DLR), or German Aerospace Center. A set of comparisons is made and the results of simulations show a good agreement among the three mechanisms as well with the experimental data.


2008 ◽  
Vol 45 ◽  
pp. 161-176 ◽  
Author(s):  
Eduardo D. Sontag

This paper discusses a theoretical method for the “reverse engineering” of networks based solely on steady-state (and quasi-steady-state) data.


Sign in / Sign up

Export Citation Format

Share Document