A superstructure based approach to chemical reactor network synthesis

Abstract Frequently the product mix in a reacting system is influenced by both heat and mixing effects in the reactor. This paper presents a nonlinear programming (NLP) formulation for optimally generating reactor networks that would produce the desired effects, given a kinetic mechanism and expressions for the reaction rate. Using a generalized configuration of recycle reactors and heat exchangers (superstructure) as a basis, the NLP is able to extract the optimal subnetwork of units that would maximize an arbitrary objective function defined over a given kinetic scheme. Thus mixing and heat effects that affect reactions occurring in the homogeneous phase are taken into account. The resulting model and adjoint equations of this formulation form a two-point boundary value problem that interfaces with an efficient optimization strategy. Decisions representing network structure, reactor type and the amount of heat addition are made through continuous parameters in the model. The superstructure allows for serial and parallel connections involving reactor units, and is therefore fairly general. Since the method is equation based and since there are numerous ODE solvers it can be applied to large kinetic mechanisms with almost any objective function. Literature examples are presented and solved in order to demonstrate the effectiveness of this approach.

[1]  Philip E. Gill,et al.  Numerically stable methods for quadratic programming , 1978, Math. Program..

[2]  F. Horn,et al.  The use of the adjoint variables in the development of improvement criteria for chemical reactors , 1967 .

[3]  J. D. Paynter,et al.  Determination of optimal reactor type , 1970 .

[4]  Ignacio E. Grossmann,et al.  An outer-approximation algorithm for a class of mixed-integer nonlinear programs , 1986, Math. Program..

[5]  L. Biegler,et al.  Algorithmic synthesis of chemical reactor networks using mathematical programming , 1986 .

[6]  Arthur E. Bryson,et al.  Applied Optimal Control , 1969 .

[7]  Lorenz T. Biegler,et al.  Developing Targets for the Performance Index of a Chemical Reactor Network: Isothermal Systems , 1988 .

[8]  A. L. Ravimohan Optimization of chemical reactor networks with respect to flow configuration , 1971 .

[9]  Shih-Ping Han A globally convergent method for nonlinear programming , 1975 .

[10]  A. M. Geoffrion Generalized Benders decomposition , 1972 .

[11]  S L Ong,et al.  Optimization of CSTRs in series by dynamic programming , 1986, Biotechnology and bioengineering.

[12]  M. J. D. Powell,et al.  A fast algorithm for nonlinearly constrained optimization calculations , 1978 .

[13]  David F. Ollis,et al.  Biochemical Engineering Fundamentals , 1976 .

[14]  Lorenz T. Biegler,et al.  Improved infeasible path optimization for sequential modular simulators—II: the optimization algorithm , 1985 .

[15]  R. Jackson Optimization of chemical reactors with respect to flow configuration , 1968 .