Model-Based Design of Experiments in Pyomo and Its Application to Adsorptive CO2 Capture Systems
- Type: Conference Presentation
- Conference Type: AIChE Annual Meeting
- Presentation Date: November 11, 2021
- Duration: 15 minutes
- Skill Level: Intermediate
- PDHs: 0.50
In this work, we present a general Python package for MBDoE using Pyomo models. Aiming at reducing the uncertainty in the system, MBDoE is conducted by minimizing the variance of the Fisher Information Matrix(FIM), which can be interpreted as minimizing the volume of the covariance ellipsoid. Classical alphabetic criteria such as D-optimality (determinant), A-optimality (trace), E-optimality (minimal eigenvalue), and a modified E-optimality (condition number) are considered to measure the size of the FIM . Nonlinear programming sensitivity analysis, available through sIPOPT  or k_aug , is used to reduce the computational cost of assembling the FIM by over an order of magnitude.
The new capability is demonstrated in two case studies. Case study 1 considers a nonlinear reaction kinetic model with four highly-correlated kinetic parameters, which shows how MBDoE formalism quickly identifies model identifiability challenges. This case study acts as a tutorial for new users. Case study 2 considers CO2 adsorption in a fixed-bed breakthrough experiment to characterize novel Metal-Organic Frameworks (MOFs) materials for CO2 capture. The partial differential-algebraic equation (PDE) model couples mass and momentum transport phenomena with adsorption equilibria (isotherms) and kinetics; discretization in space (method of lines) and time (backward finite difference or collocation) results in over 20,000 sparse resulting in over 20,000 sparse algebraic constraints.
The goal of MBDoE is to infer the heat transfer coefficient (non-adiabatic operation) and a lumped kinetics transport parameter. Through MBDoE, we address the following questions: (1) Is the model identifiable with the current experimental configuration? (2) What is the value of modifying the experimental system (e.g., adding temperature sensors)? (3) How to conduct fixed bed experiments so that we can improve the accuracy of parameter estimation by DoE? The limited literature that considers DoE for a fixed bed adsorption system focuses on âblack-box methodsâ such as factorial methods and response surfaces methods which do not offer insights into the fundamental physical process, in contrast to our proposed MBDoE approach.
In the future, validation mathematical models will be used for optimization and techno-economic analysis of novel CO2 capture processes as part of the Carbon Capture Simulation for Industrial Impact (CCSI2) project.
Franceschini, G., & Macchietto, S. (2008). Model-based design of experiments for parameter precision: State of the art. Chemical Engineering Science, 63(19), 4846-4872.
Waldron, C., Pankajakshan, A., Quaglio, M., Cao, E., Galvanin, F. and Gavriilidis, A., 2020. Model-based design of transient flow experiments for the identification of kinetic parameters. Reaction Chemistry & Engineering, 5(1), pp.112-123.
Balsa-Canto, E., Rodriguez-Fernandez, M. and Banga, J.R., 2007. Optimal design of dynamic experiments for improved estimation of kinetic parameters of thermal degradation. Journal of Food Engineering, 82(2), pp.178-188.
Chakrabarty, A., Buzzard, G.T. and Rundell, A.E., 2013. Modelâbased design of experiments for cellular processes. Wiley Interdisciplinary Reviews: Systems Biology and Medicine, 5(2), pp.181-203.
Box, G.E. and Lucas, H.L., 1959. Design of experiments in non-linear situations. Biometrika, 46(1/2), pp.77-90.
Pirnay, H., LÃ³pez-Negrete, R. and Biegler, L.T., 2012. Optimal sensitivity based on IPOPT. Mathematical Programming Computation, 4(4), pp.307-331.
Thierry, D. and Biegler, L.T., 2019. Dynamic realâtime optimization for a CO2 capture process. AIChE Journal, 65(7), p.e16511.
Nocedal, J., WÃ¤chter, A. and Waltz, R.A., 2009. Adaptive barrier update strategies for nonlinear interior methods. SIAM Journal on Optimization, 19(4), pp.1674-1693.
WÃ¤chter, A., 2002. An interior point algorithm for large-scale nonlinear optimization with applications in process engineering (Doctoral dissertation, PhD thesis, Carnegie Mellon University).
HSL, A., 2007. collection of Fortran codes for large-scale scientific computation. See http://www. hsl. rl. ac. Uk.
 Dowling, A.W., Vetukuri, S.R. and Biegler, L.T., 2012. Largeâscale optimization strategies for pressure swing adsorption cycle synthesis. AIChE journal, 58(12), pp.3777-3791.
|AIChE Member Credits||0.5|
|AIChE Graduate Student Members||Free|
|AIChE Undergraduate Student Members||Free|