Model complexity reduction of chemical reaction networks using mixed-integer quadratic programming. (English) Zbl 1339.92103
Summary: The model complexity reduction problem of large chemical reaction networks under isobaric and isothermal conditions is considered. With a given detailed kinetic mechanism and measured data of the key species over a finite time horizon, the complexity reduction is formulated in the form of a mixed-integer quadratic optimization problem where the objective function is derived from the parametric sensitivity matrix. The proposed method sequentially eliminates reactions from the mechanism and simultaneously tunes the remaining parameters until the pre-specified tolerance limit in the species concentration space is reached. The computational efficiency and numerical stability of the optimization are improved by a pre-reduction step followed by suitable scaling and initial conditioning of the Hessian involved. The proposed complexity reduction method is illustrated using three well-known case studies taken from reaction kinetics literature.
MSC:
| 92E20 | Classical flows, reactions, etc. in chemistry |
| 80A30 | Chemical kinetics in thermodynamics and heat transfer |
Keywords:
chemical reaction networks; model reduction; mixed-integer quadratic programming; sensitivity analysisReferences:
| [1] | Theodoropoulos, C., Optimisation and linear control of large scale nonlinear systems: a review and a suite of model reduction-based techniques, (Gorban, A. N.; Roose, D., Coping with Complexity: Model Reduction and Data Analysis. Coping with Complexity: Model Reduction and Data Analysis, Lecture Notes in Computational Science and Engineering, vol. 75 (2011), Springer: Springer Heidelberg, Dordrecht, New York), 37-61 |
| [2] | Xie, W.; Bonis, I.; Theodoropoulos, C., Off-line model reduction for on-line linear MPC of nonlinear large-scale distributed systems, Computers and Chemical Engineering, 35, 750-757 (2011) |
| [3] | Hangos, K. M., Engineering model reduction and entropy-based Lyapunov functions in chemical reaction kinetics, Entropy, 12, 772-797 (2010) · Zbl 1229.80021 |
| [4] | Toth, J.; Li, G.; Rabitz, H.; Tomlin, A. S., The effect of lumping and expanding on kinetic differential equations, SIAM Journal on Applied Mathematics, 57, 6, 1531-1556 (1997) · Zbl 0891.34030 |
| [5] | Lall, S.; Marsden, J. E.; Glavaski, S., A subspace approach to balanced truncation for model reduction of nonlinear control systems, International Journal of Robust and Nonlinear Control, 12, 6, 519-535 (2002) · Zbl 1006.93010 |
| [6] | Anderson, J.; Chang, Y.-C.; Papachristodoulou, A., Model decomposition and reduction tools for large-scale networks in systems biology, Automatica, 47, 1165-1174 (2011) · Zbl 1235.93018 |
| [7] | Gorban, A. N.; Karlin, I. V.; Zinovyev, A. Y., Invariant grids for reaction kinetics, Physica A, 33, 106-154 (2004) |
| [8] | Androulakis, I. P., Kinetic mechanism reduction based on an integer programming approach, AIChE Journal, 46, 361-371 (2000) |
| [9] | Csendes, T., Nonlinear parameter estimation by global optimization—efficiency and reliability, Acta Cybernetica, 8, 361-370 (1988), (Szeged) · Zbl 0654.90074 |
| [10] | Edwards, K.; Edgar, T. F.; Manousiouthakis, V. I., Reaction mechanism simplification using mixed-integer nonlinear programming, Computers and Chemical Engineering, 24, 67-79 (2000) |
| [11] | Maurya, M. R.; Katare, S.; Patkar, P. R.; Rundell, A. E.; Venkatasubramanian, V., A systematic framework for the design of reduced-order models or signal transduction pathways from a control theoretic perspective, Computers and Chemical Engineering, 30, 437-452 (2006) |
| [12] | Radulescu, O.; Gorban, A. N.; Zinovyev, A.; Lilienbaum, A., Robust simplifications of multiscale biochemical networks, BMC Systems Biology, 2, 1, 86 (2008) |
| [13] | Horn, F.; Jackson, R., General mass action kinetics, Archive for Rational Mechanics and Analysis, 47, 81-116 (1972) |
| [14] | Maurya, M. R.; Bornheimer, S. J.; Venkatasubramanian, V.; Subramaniam, S., Mixed-integer nonlinear optimisation approach to coarse-graining biochemical networks, IET Systems Biology, 3, 1, 24-39 (2009) |
| [15] | Mitsos, A.; Oxberry, G. M.; Barton, P. I.; Green, W. H., Optimal automatic reaction and species elimination in kinetic mechanisms, Combustion and Flame, 155, 118-132 (2008) |
| [16] | Chellaboina, V.; Bhat, S. P.; Haddad, W. M.; Bernstein, D. S., Modeling and analysis of mass-action kinetics—nonnegativity, realizability, reducibility, and semistability, IEEE Control Systems Magazine, 29, 60-78 (2009) · Zbl 1395.93098 |
| [17] | Érdi, P.; Tóth, J., Mathematical Models of Chemical Reactions. Theory and Applications of Deterministic and Stochastic Models (1989), Manchester University Press, Princeton University Press: Manchester University Press, Princeton University Press Manchester, Princeton · Zbl 0696.92027 |
| [18] | Feinberg, M., On chemical kinetics of a certain class, Archive for Rational Mechanics and Analysis, 46, 1-41 (1972) |
| [19] | Llaneras, F.; Picó, J., Stoichiometric modelling of the cell metabolism, Journal of Bioscience and Bioengineering, 105, 1-11 (2008) |
| [20] | Otero-Muras, I.; Franco-Uría, A.; Alonso, A. A.; Balsa-Canto, E., Dynamic multi compartmental modelling of metal bioaccumulation in fish: identifiability implications, Environmental Modeling & Software, 25, 344-353 (2010) |
| [21] | Angeli, D., A tutorial on chemical network dynamics, European Journal of Control, 15, 398-406 (2009) · Zbl 1298.92126 |
| [22] | Lin, Y.-C.; Fan, L. T.; Shafiea, S.; Bertók, B.; Friedler, F., Generation of light hydrocarbons through Fischer-Tropsch synthesis: identification of potentially dominant catalytic pathways via the graph-theoretic method and energetic analysis, Computers and Chemical Engineering, 33, 1182-1186 (2009) |
| [23] | Vajda, S.; Valko, P.; Turányi, T., Principal component analysis of kinetic models, International Journal of Chemical Kinetics, 17, 55-81 (1985) |
| [24] | Turányi, T., Applications of sensitivity analysis to combustion chemistry, Reliability Engineering and System Safety, 57, 41-49 (1997) |
| [25] | Nagy, T.; Turányi, T., Reduction of very large reaction mechanisms using methods based on simulation error minimization, Combustion and Flame, 156, 417-428 (2009) |
| [26] | Huang, H.; Fairweather, M.; Griffiths, J. F.; Tomlin, A. S.; Brad, R. B., A systematic lumping approach for the reduction of comprehensive kinetic models, Proceedings of the Combustion Institute, 30, 1309-1316 (2005) |
| [27] | Apri, M.; de Gee, M.; Molenaar, J., Complexity reduction preserving dynamical behavior of biochemical networks, Journal of Theoretical Biology, 304, 16-26 (2012) · Zbl 1397.92248 |
| [28] | Feinberg, M., Chemical reaction network structure and the stability of complex isothermal reactors—I. The deficiency zero and deficiency one theorems, Chemical Engineering Science, 42, 2229-2268 (1987) |
| [29] | Snow, R. H., A chemical kinetics computer program for homogeneous and free-radical systems of reactions, The Journal of Physical Chemistry, 70, 9, 2780-2786 (1966) |
| [30] | Turányi, T., Reduction of large reaction mechanism, New Journal of Chemistry, 14, 794-803 (1990) |
| [31] | Hannemann, R.; Marquardt, W.; Gendler, B.; Naumann, U., Discrete first- and second-order adjoints and automatic differentiation for the sensitivity analysis of dynamic models, Procedia Computer Science, 1, 297-305 (2010) |
| [32] | Brendel, M.; Bonvin, D.; Marquardt, W., Incremental identification of kinetic models for homogeneous reaction systems, Chemical Engineering Science, 61, 16, 5404-5420 (2006) |
| [33] | Anderson, E.; Bai, Z.; Bischof, C.; Blackford, S.; Demmel, J.; Dongarra, J.; Du Croz, J.; Greenbaum, A.; Hammarling, S.; McKenney, A.; Sorensen, D., LAPACK Users’ Guide (1999), Society for Industrial and Applied Mathematics: Society for Industrial and Applied Mathematics Philadelphia, PA · Zbl 0934.65030 |
| [34] | Hannemann-Tamás, R.; Tillack, J.; Schmitz, M.; Förster, M.; Wyes, J.; Nöh, K.; von Lieres, E.; Naumann, U.; Wiechert, W.; Marquardt, W., First- and second-order parameter sensitivities of a metabolically and isotopically non-stationary biochemical network model, (Otter, M.; Zimmer, D., Proceedings of the 9th International Modelica Conference (2012), Modelica Association), 641-648 |
| [35] | Naumann, U., (The Art of Differentiating Computer Programs—An Introduction to Algorithmic Differentiation. The Art of Differentiating Computer Programs—An Introduction to Algorithmic Differentiation, Software, Environments, Tools, vol. 24 (2012), SIAM) · Zbl 1275.65015 |
| [36] | Fourer, R.; Gay, D. M.; Kernighan, B. W., AMPL: A Mathematical Programming Language (1989), AT&T Bell Laboratories |
| [37] | ILOG Inc., ILOG AMPL CPLEX System—User’s Guide. Version 11.0 (2008) |
| [38] | Vardanyan, I. A.; Sachyan, G. A.; Nalbandyan, A. B., The rate constant of the reaction \(HO_2\)+CO⇄\(CO_2\)+OH, International Journal of Chemical Kinetics, 7, 23-31 (1975) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
