Modelling, mathematical analysis, numerical solution and parameter identification in reaction systems

This thesis is divided in three different parts in which an extensive study of reactors is done. It is complemented with three appendices describing some tools and results used . The first part is devoted to the description of the reactors we are interested in. We formulate the models of the main STR reactors. Then, we describe the general convection-diffusion-reaction model. Finally, the FBR model is also described. In second part an extensive study is done for the convection-diffusion-reaction model beginning with the mathematical analysis for the n-dimensional reactor and then numerical solution of the reactor models is designed. The last part deals with the identification of the best kinetic model from a list of proposed functional forms, and also of the values of its corresponding parameters by means of an optimization process. For this purpose, we use a combination of an incremental and an integral method

Concept of Variants and Invariants for Reaction Systems, with Application to Estimation, Control and Optimization

The concept of reaction variants and invariants for lumped reaction systems has been known for several decades. Its applications encompass model identification, data reconciliation, state estimation and control using kinetic models. In this thesis, the concept of variants and invariants is extended to distributed reaction systems and used to develop new applications to estimation, control and optimization. The thesis starts by reviewing the material and heat balances and the concept of variants and invariants for several lumped reaction systems. Different definitions of variants and invariants, in particular the vessel extents, are presented for the case of homogeneous reaction systems, and transformations to variants and invariants are obtained. The extension to systems with heat balance and mass transfer is also reviewed. The concept of extents is generalized to distributed reaction systems, which include many processes involving reactions and described by partial differential equations. The concept of extents and the transformation to extents are detailed for various configurations of tubular reactors and reactive separation columns, as well as for a more generic framework that is independent of the configuration. New developments of the extent-based incremental approach for model identification are presented. The approach, which compares experimental and modeled extents, results in maximum-likelihood parameter estimation if the experimental extents are uncorrelated and the modeled extents are unbiased. Furthermore, the identification problem can be reformulated as a convex optimization problem that is solved efficiently to global optimality. The estimation of unknown rates without the knowledge or the identification of the rate models is described. This method exploits the fact that the variants computed from the available measurements allow isolating the different rates. Upon using a Savitzky-Golay filter for differentiation of variants, one can show that the resulting rate estimator is optimal and obtain the error and variance of the rate estimates. The use of variants and invariants for reactor control is also considered. Firstly, offset-free control via feedback linearization is implemented using kinetic models. Then, it is shown how rate estimation can be used for control via feedback linearization without kinetic models. By designing an outer-loop feedback controller, the expected values of the controlled variables converge exponentially to their setpoints. This thesis presents an approach to speed up steady-state optimization, which takes advantage of rate estimation without rate models to speed up the estimation of steady state for imperfectly known dynamic systems with fast and slow states. Since one can use feedback control to speed up convergence of the fast part, rate estimation allows estimating the steady state of the slow part during transient operation. The application to dynamic optimization is also shown. Adjoint-free optimal control laws are computed for all the types of arcs in the solution. In the case of reactors, the concept of extents allows the symbolic computation of optimal control laws in a systematic way. A parsimonious input parameterization is presented, which approximates the optimal inputs well with few parameters. For each arc sequence, the optimal parameter values are computed via numerical optimization. The theoretical results are illustrated by simulated examples of reaction systems

Efficient hybrid modeling and sorption model discovery for non-linear advection-diffusion-sorption systems: A systematic scientific machine learning approach

This study presents a systematic machine learning approach for creating efficient hybrid models and discovering sorption uptake models in non-linear advection-diffusion-sorption systems. It demonstrates an effective method to train these complex systems using gradientbased optimizers, adjoint sensitivity analysis, and JIT-compiled vector Jacobian products, combined with spatial discretization and adaptive integrators. Sparse and symbolic regression were employed to identify missing functions in the artificial neural network. The robustness of the proposed method was tested on an in-silico data set of noisy breakthrough curve observations of fixed-bed adsorption, resulting in a well-fitted hybrid model. The study successfully reconstructed sorption uptake kinetics using sparse and symbolic regression, and accurately predicted breakthrough curves using identified polynomials, highlighting the potential of the proposed framework for discovering sorption kinetic law structures.Comment: Preprint paper to be submitted soon in Elsevier Journa