CALCULATING LIMITS TO PRODUCTIVITY IN REACTOR-SEPARATOR
SYSTEMS OF ARBITRARY DESIGN
DISSERTATION
Presented in Partial Fulfillment of the Requirements for
The Degree Doctor of Philosophy in the Graduate School of
The Ohio State University
By
Yangzhong Tang, B.S.
* * * * * *
The Ohio State University
2005
Dissertation Committee: Approved by
Dr. Martin R. Feinberg, Adviser
Dr. James F. Rathman Adviser Graduate Program in Dr. Bhavik R. Bakshi Chemical Engineering
ABSTRACT
This thesis aims to solve the following major problem in the process synthesis
area of chemical engineering: Given a network of chemical reactions with its associated
kinetics, given certain available resources such as the catalyst mass, given specific feed rates, given permissible temperature and pressure ranges within reactor units, and given
certain environmental constraints, what is the maximum production rate that can be
achieved for a certain desired species over all possible steady-state reactor-separator systems consistent with those specifications?
The most difficult part in solving this problem is the huge varieties of the reactor- separator configurations that can be applied in the process. Besides the classical types of reactors that are widely known (i.e. continuous flow stirred tank reactors or plug flow reactors), there could be reactors that are beyond our imagination. So on top of the complexity already present in consideration of the great variety of already-existent reactor-separator systems, one must also allow for consideration of the as-yet-unknown types of reactors.
ii While it seems that it is almost impossible to solve the problem, an idea called the
CFSTR Equivalence Principle for Reactor-Separator Systems was provided in 2001 by
Feinberg and Ellison as a very important tool. This principle says that to sort out all attainable effluents, it suffices to look at an exceptional design in which the only reactors are continuous flow stirred tank reactors (CFSTRs) whose number is determined by the network of chemical reactions. So the search for the maximum production rate of a certain desired species can be carried out among different variations of this CFSTR-only design rather than among the vast spectrum of different combinations of many types of reactors and separators.
However, the CFSTR Equivalence Principle only serves as a theoretical tool for the problem. What this thesis tries to do is to provide concrete implementation methods, based on the principle, for searching the kinetic bounds on productivity for stipulated resources and specified process constraints. Five linear methods and one nonlinear method are explored for this purpose and a software package called Productivity Limit
Calculator (PLC) is developed with all these methods imbedded.
Case studies are carried out by using the software package both for academic test examples and for more practical industrial processes. It can be seen that in general the nonlinear method is better than the linear methods in performance. The nonlinear method uses a commercial-grade solver for the search of the maximum production rate. It gives good results for both large and small-sized problems. Although linear methods are only
iii good for small-sized problems, they do not depend on any commercial-grade solver, and their underlying algorithms are freely available. For this reason, the linear methods are useful in educational settings.
Besides searching for the kinetic bounds on productivity in reactor-separator
systems of arbitrary design, the software also provides means to examine certain
theoretical questions. In particular, it becomes possible to examine how changes in
available reactor capacity or how tightening of environmental constraints affect the
maximum possible production rate of some desired chemical species. These issues are
examined in the context of two examples.
iv
DEDICATED TO MY PARENTS
v
ACKNOWLEDGMENTS
I would like to express my deepest appreciation to my advisor, Dr. Martin R.
Feinberg, for his advice and guidance throughout my research period. I would also like to thank him for his patient criticism during the writing of this thesis. Working with Dr.
Feinberg has been the most challenging and rewarding experience in my life. He has inspired me in many different ways. Besides being an advisor, Dr. Feinberg is also like a father to me for all these years with his support and encouragement.
My group-mate Thomas Abraham has always been very nice and helpful. It is amazing that he knows so many things besides chemical engineering that I always went to seek answers from him in the first place. Now that he graduated recently and got a job,
I wish him and his families all of the best.
Finally, I am grateful to Chemical and Biomolecular Engineering Department at
The Ohio State University and United States National Science Foundation for their financial supports of my doctoral studies.
vi
VITA
June 25, 1978…………………………… Born – Zhuzhou, Hunan, China
2000…………………………………….. B.S., Control Science and Engineering,
Zhejiang University,
Hangzhou, Zhejiang, China
2000 – 2004…………………………….. Department Fellow,
Chemical Engineering,
The Ohio State University,
Columbus, OH
2004 – 2005…………………………….. Graduate Research Associate
Chemical Engineering,
The Ohio State University,
Columbus, OH
FIELDS OF STUDY
Major Field: Chemical Engineering
vii
TABLE OF CONTENTS
Page
ABSTRACT...... ii
DEDICATION...... v
ACKNOWLEDGMENTS ...... vi
VITA...... vii
LIST OF FIGURES ...... xii
LIST OF TABLES...... xiv
LIST OF ABBREVIATIONS...... xvi
NOMENCLATURE ...... xvii
Chapters:
1. INTRODUCTION ...... 1
1.1 A Big Question in Chemical Engineering ...... 1 1.2 Current Progress in Process Synthesis...... 2 1.2.1 Progress in the Pure Reactor Synthesis Area...... 3 1.2.1.1 Geometric Approach...... 3 1.2.1.2 The Superstructure and Mixed-integer Nonlinear Programming Approach7 1.2.1.3 Target-based Approach...... 10 1.2.1.4 Stochastic Optimization Approach and Other Approaches ...... 11 1.2.2 Progress in the Reactor-Separator Synthesis Area...... 12 viii 1.3 Finding Kinetic Bounds on Productivity in Reactor-Separator Systems of Arbitrary Design ...... 16 1.3.1 Addressing the Problem...... 16 1.3.2 An Important Conceptual Tool The CFSTR Equivalence Principle for Reactor-Separator Systems ...... 21 1.3.3 Beyond the Theory One Step Further ...... 24 1.4 Summary of Chapter 1...... 25
2. THE CFSTR EQUIVALENCE PRINCIPLE FOR REACTOR-SEPARATOR SYSTEMS...... 27
2.1 Some Preliminaries...... 28 2.1.1 A Few Mathematical Preliminaries ...... 28 2.1.2 Stoichiometric and Kinetic Preliminaries ...... 29 2.2 The CFSTR Equivalence Principle for Reactor-Separator Systems...... 32 2.3 Applying the CFSTR Equivalence Principle to Find the Kinetic Bounds on Productivity in Reactor-Separator Systems of Arbitrary Design...... 37 2.3.1 A Few More Mathematical Preliminaries...... 37 2.3.2 How the CFSTR Equivalence Principle Provides Kinetic Bounds...... 39 2.3.3 Constraint Set Ω...... 40 2.3.4 The Role of the CFSTR Equivalence Principle in Finding the Kinetic Bounds ...... 44 2.3.4.1 The CFSTR Equivalence Principle Characterizes the Attainable Effluent Rate Vectors...... 44 2.3.4.2 Another Equivalent Form of Expressing the Attainable Effluents ...... 48 2.3.4.3 Applications of the Two Forms ...... 49 2.3.4.3.1 Application 1: How to Tell If a Molar Effluent Rate Vector is Attainable? ...... 49 2.3.4.3.2 Application 2: Finding the Maximum Productivity without Environmental Constraints...... 50 2.3.4.3.3 Application 3: Finding the Maximum Productivity with Environmental Constraints...... 52 2.4 Summary of Chapter 2...... 55
ix 3. IMPLEMENTATION METHODS FOR SEARCHING FOR KINETIC BOUNDS ON PRODUCTIVITY IN REACTOR-SEPARATOR SYSTEMS OF ARBITRARY DESIGN...... 56
3.1 A Sample Problem...... 57 3.2 The General Problem...... 64 3.2.1 Statement A...... 68 3.2.2 Statement B...... 70 3.3 Linear Methods...... 72 3.3.1 The Qhull Once Algorithm ...... 74 3.3.1.1 Some Preliminaries ...... 74 3.3.1.1.1 Introduction to Qhull...... 74 3.3.1.1.2 The Stoichiometric Subspace...... 79 3.3.1.2 Implementation of the Qhull Once Algorithm...... 85 3.3.1.2.1 First Step: Discretization of the Constraint Set Ω...... 87 3.3.1.2.2 Second Step: Computation of the Convex Hull from Ω* ...... 98 3.3.1.2.3 Third Step: Finding the Maximum Production Rate...... 104 3.3.1.2.4 Summary of the Qhull Once Algorithm...... 110 3.3.2 The Qhull Twice Algorithm...... 112 3.3.3 The Linear Programming (LP) Only Algorithm...... 116 3.3.4 The Qhull Once + Linear Programming Algorithm ...... 120 3.3.5 The Multi-stage Linear Programming Algorithm...... 124 3.3.6 Summary of All Linear Methods ...... 126 3.4 A Nonlinear Method ...... 128
4. A SOFTWARE PACKAGE: THE PRODUCTIVITY LIMIT CALCULATOR...... 132
4.1 Demonstration of the Software Package by the van de Vusse Example ...... 132 4.2 Efficient Computation in the Nonideal Gas Mixture Case...... 145
5. CASE STUDIES...... 153
5.1 The van de Vusse Example Revisited...... 154 5.2 Other Network Examples...... 157 5.3 Summary of the Discussions Based On the Results of the Previous Examples ... 179 x 5.4 Practical Examples...... 181 5.4.1 Cumene Production...... 181 5.4.2 Allyl Chloride Production...... 187 5.4.3 Acrylic Acid Production ...... 193 5.4.4 Phthalic Anhydride Production...... 199 5.5 Summary of All Case Studies...... 205
6. CONCLUSION AND FUTURE RESEARCH WORK ...... 207
6.1 Conclusion ...... 207 6.1.1 An Important Conceptual Tool – The CFSTR Equivalence Principle ...... 209 6.1.2 Different Algorithms in the Software Package Productivity Limit Calculator ...... 211 6.2 Suggested Future Research Work...... 212 6.2.1 Effect of the Environmental Constraint on the Productivity Limit...... 213 6.2.2 Effect of the Reactor Capacity on the Attainable Region...... 221 6.2.3 Sensitivity of the Result to Errors from Different Sources...... 224 6.2.4 True Maximum ...... 224 6.2.5 Improving the Current Design ...... 225
REFERENCES ...... 226
xi
LIST OF FIGURES
Figure Page
1.1 A chemical process with an unspecified design…………………………………. 19
1.2 Varieties of reactor-separator systems…………………………………………… 20
2.1 Illustration of the CFSTR Equivalence Principle………………………………... 34
2.2 Illustration of convex hull………………………………………………………... 38
2.3 The CFSTR-only design…………………………………………………………. 45
3.1 A 4-CFSTR design for the van de Vusse network……………………………….. 59
3.2 A (s+1)-CFSTR design for the general case……………………………………... 66
3.3 A five-point input to Qhull………………………………………………………. 75
3.4 Output of Qhull for the five-point input…………………………………………. 77
3.5 Four points given in R3 live in a 2-d subspace…………………………………... 78
3.6 Illustration of the Qhull Once Algorithm………………………………………… 87
3.7 Valid (c1, c2) combinations in the shadowed area………………………………... 92
3.8 Illustration of discretization of the concentration vector………………………… 94
3.9 A simple discretization of Ω in the van de Vusse example……………………… 98
3.10 A linear map from S to Rs…………………………………………………….… 100
3.11 Comparison of Qhull Once and Qhull Twice Algorithm………………………..114
xii 3.12 Comparison of Qhull Once + LP and Qhull Twice Algorithm…………………. 121
3.13 Illustration of the Multi-stage Linear Programming Algorithm………………... 126
5.1 Cumene production process diagram…………………………………………… 182
5.2 Illustration of the two types of designs in cumene production…………………. 185
5.3 Cumene production rates at different operating temperatures………………….. 186
5.4 Allyl chloride production process flow diagram……………………………….. 188
5.5 Allyl chloride production rates at different operating temperatures……………. 192
5.6 Acrylic acid production process diagram………………………………………. 194
5.7 Phthalic anhydride production process flow diagram…………………………... 200
6.1 Effect of the M3 constraint on the M4 maximum in the Network Example D
in Section 5.2…………………………………………………………………….215
6.2 Effect of the M3 constraint on the M4 maximum in the van de Vusse example
in Section 5.1…………………………………………………………………….217
6.3 Illustration of possible positions of the critical environmental constraint……… 219
6.4 Attainable regions at different reactor volumes for
the van de Vusse example……….……………………………………………… 221
xiii
LIST OF TABLES
Table Page
3.1 Occurrence rate of each reaction in a small reaction network…………………… 80
3.2 Occurrence rate function of each reaction in a small reaction network………….. 81
5.1 Results for the van de Vusse example…………………………………………...155
5.2 Result for the maximum M4 (mol/s) with no limit on M3 in Network A………..160
5.3 Result for the maximum M4 (mol/s) with limit on M3 to be
less than 0.25 mol/s in Network A……………………………………………… 160
5.4 Result for the maximum M4 (mol/s) with no limit on M3 in Network B……….. 165
5.5 Result for the maximum M4 (mol/s) with limit on M3 to be
less than 0.1 mol/s in Network B…………………………………………….…. 165
5.6 Result for the maximum M3 (mol/s) with no other limit in Network C………....168
5.7 Result for the maximum M3 (mol/s) with limit on M2 to be
less than 0.5 mol/s in Network C……………………………………………….. 168
5.8 Result for the maximum M3 (mol/s) with limits on both M2 and M4
to be less than 0.5 mol/s in Network C…………………………………………. 168
5.9 Result for the maximum M4 (mol/s) with no other limit in Network C………... 169
5.10 Result for the maximum M4 (mol/s) with limit on M3 to be
xiv less than 1.0 mol/s in Network C……………………………………………….. 169
5.11 Result for the maximum M4 (mol/s) with limit on M3 to be
less than 0.1 mol/s in Network D……………………………………………….. 173
5.12 Detailed information of the two CFSTRs which gives
the maximum M4 in Network Example D……………………………………… 175
5.13 Manually calculated results of the effluent rates in Network Example D…….... 176
5.14 Process specifications/constraints for cumene production………………………183
5.15 Critical temperature, critical pressure, and Pitzer acentric factor
of all the involved species in the cumene production case……………………... 185
5.16 Process specifications for allyl chloride production……………………………. 190
5.17 Critical temperature, critical pressure, and Pitzer acentric factor
of all the involved species in the allyl chloride production case……………….. 191
5.18 Feed streams information for acrylic acid production………………………….. 193
5.19 Kinetic information for acrylic acid production………………………………... 196
5.20 Process specifications for acrylic acid production……………………………… 196
5.21 Critical temperature, critical pressure, and Pitzer acentric factor
of all the involved species in the acrylic acid production case…………………. 198
5.22 Process specifications for phthalic anhydride production……………………… 202
5.23 Critical temperature, critical pressure, and Pitzer acentric factor
of all the involved species in the phthalic anhydride production case………….. 203
xv
LIST OF ABBREVIATIONS
AR – Attainable region
CFSTR – Continuous flow stirred tank reactor
DLL – Dynamic link library
DSR – Differential side-stream reactor
LP – Linear programming
MINLP – Mixed-integer nonlinear programming
NLP – Nonlinear programming
PFR – Plug flow reactor
PLC – Productivity limit calculator
xvi
NOMENCLATURE
c composition vector
th ci concentration of the i species conv convex hull
M0 molar feed rate vector
M molar effluent rate vector
M the set of all attainable effluent rate vectors
N number of species in a chemical reaction network
P pressure
Pc critical pressure
Pmin minimum allowable pressure
Pmax maximum allowable pressure r(•,•) volumetric species formation rate vector (or function)
RN N-dimensional vector space
N N R + a part of R where every vector has non-negative components s rank of a chemical reaction network
S stoichiometric subspace
T temperature xvii Tc critical temperature
Tmin minimum allowable temperature
Tmax maximum allowable temperature
T a linear map from the stoichiometric subspace S to space Rs
th vi partial molar volume of the i species
V total reactor volume
V maximum total reactor volume
V interval [0, V ]
ω Pitzer acentric factor
Ω constraint set
Ω* discretized constraint set
xviii
CHAPTER 1
INTRODUCTION
1.1 A Big Question in Chemical Engineering
There has been a long-standing question in the process synthesis area of chemical engineering that remained unsolved, and this thesis attempts to answer it in a general way.
Roughly, the question is: Given a system of chemical reactions with its associated kinetics, given certain available resources (e.g., available feed stream, available catalyst mass), and given certain constraints including the temperature-pressure ranges in the reactor units and applicable environmental constraints, what is the maximum production rate for a certain desirable species over all possible steady-state reactor-separator designs consistent with all these specifications?
There are certain ideas that can be used to bound the production rate. One of them is the stoichiometric limit, derived from the mass balance according to the underlying chemical reaction network. However, the stoichiometric limit does not take into account the kinetic information of the underlying chemical reaction network and the available reactor capacity. Therefore, it cannot provide a sharp production bound for a specified
1 reactor size. This thesis tries to provide a much more precise production limit by considering all factors including the kinetics of the underlying chemical reaction network, the available resources (e.g., the available feed streams, the mass of catalyst used in the reactor), and all other constraints including the temperature and pressure ranges in the reactors and applicable environmental constraints.
By saying all possible steady-state reactor-separator configurations, we are considering not only the traditional reactor types such as the continuous flow stirred tank reactor (CFSTR), the plug flow reactor (PFR), and combinations of these, but also many other types including ones that are even unknown to us at this time. Such a huge selection of reactor-separator systems makes the search for the maximum production rate of a certain desired species extremely difficult and seems almost unsolvable. However, this thesis is going to face the challenge and provide ways to answer this big question in chemical engineering.
Before giving a precise statement of the problem that this thesis tries to solve, it is worthwhile to first review the other related work in the process synthesis area.
1.2 Current Progress in Process Synthesis
For many years scholars from different countries have tried different ways and have achieved a lot in solving problems in process design and process synthesis. Some of the processes studied involve only reaction and mixing, and some processes involve reaction, mixing, and separation. For those processes involving only reaction and mixing, we call the problem that is associated with them the pure reactor synthesis problem. And
2 for those processes involving reaction, mixing, and separation, we call the related problem the reactor-separator synthesis problem. In the following sections, achievements and progress in these two areas will be summarized.
1.2.1 Progress in the Pure Reactor Synthesis Area
1.2.1.1 Geometric Approach
A new concept called the attainable region (AR) was proposed by F.J.M. Horn in
1960’s for the purpose of trying to study process design from a geometric point of view
(Horn, 1964). Generally speaking, the attainable region is the set of all possible output states of a process (with a specific feed stream) by using all possible designs that can be subject to specified constraints. If the attainable region could be known, the process designers would be able to tell which output is achievable and which is not. For example, by searching on the boundary of the set of all possible steady-state production rate vectors, the maximum production rate of a certain desirable chemical species over all possible steady-state designs could be obtained. And by comparing that maximum production rate with the rate for a candidate design, designers would be able to tell how far the production rate of the candidate design is from the best that could be achieved.
Based on the attainable region idea, Glasser, Crowe and Hildebrandt (1987) tried to find all possible composition vectors that can be achieved from a specified feed composition vector by all processes involving only reaction and mixing, given a chemical reaction scheme and its kinetics. (Here, assuming that there are N species A1, …, AN
3 involved in the chemical reactions, then the attainable region is the set of all composition
N vectors [c1, …, cN] in R achievable from a given feed composition. Here, ci represents the concentration of species Ai, i = 1, …, N.) Only isothermal systems with no density change during reaction and mixing were considered. Some necessary conditions were derived for the attainable region, one of which being that the attainable region must be convex with non-zero reaction rate vectors on its boundary not pointing outward. Based on this necessary condition, Glasser, Crowe, and Hildebrandt (1987) claimed to construct a region in RN that could not be further extended by any combination of plug flow reactors (PFRs), continuous flow stirred tank reactors (CFSTRs), and recycle reactors.
Later, Hildebrandt, Glasser, and Crowe (1990) claimed to construct another region, extending the previous results for isothermal, fixed-volumed systems (Glasser, Crowe and Hildebrandt, 1987) to include adiabatic and variable-density systems. By only considering certain types of reactor systems in segregated, maximum mixed, and other reactor models, Glasser, Hildebrandt, and Godorr (1994) compared the achieved conversions of these proposed reactor systems with the conversions of a broader selection of reactor systems. In a summary about the synthesis of chemical reactor networks with heat integration (Hildebrandt and Biegler, 1995), two approaches were compared. One approach used a reactor superstructure idea to try to find the optimal reactor structure within a reactor superstructure for a certain performance index, which made the problem appear as an optimization problem. (This superstructure idea will be explained more in
Section 1.2.1.2.) Another approach tried to construct the whole attainable region based on a chosen performance target. The second approach was a geometric approach, and it was
4 severely limited by the dimensionality of the problem. Normally only 2- or 3-dimensional problems could be approached in this way. Nevertheless, this geometric approach gave some insight into the general properties of the optimal reactor structures. The first approach, which was based on the superstructure idea, although not as rigorous as the geometric approach, could extend the problems to higher dimensions. In one application,
Nicol, Hildebrandt, and Glasser (1997) used the attainable region technique to search for the optimal process design for an exothermic reaction system with external cooling consisting of two external cooling utilities. The optimal process structure they claimed to find consisted of a CFSTR followed by a PFR, a cold shot converter, and then the two external cooling utilities. In another application of the attainable region idea, the task choosing an optimal reactor structure as well as optimal operating temperature was assisted by looking at the geometry of the boundary of a candidate attainable region
(Godorr et al., 1999). A practical industrial example based on the attainable region idea was studied in the oxidative dehydrogenation of 1-butene in inert porous membrane reactors (Milne et al., 2004). There, a candidate region was proposed in order to search for the optimal process operating conditions like the partial pressure of oxygen. A systematic approach of trying to construct the attainable region was proposed by Rooney,
Hausberger, Biegler, and Glasser (2000). The idea was to build the attainable region from its lower dimensional attainable region projections. And Kauchali et al. (2004) tried to find the attainable region for a water-gas shift reaction system.
5
Feinberg and Hildebrandt (1997) studied general properties of the attainable region for all possible steady-state designs that involve only reaction and mixing, given a specified feed and the associated reaction network kinetics. It was found that despite the wide varieties of reactors which could be connected in countless ways, it is always true that manifolds of extreme points on the attainable region’s boundary are made up of PFR composition trajectories. Moreover, CFSTRs and differential sidestream reactors (DSRs) provide ways to access those PFR trajectories. Feinberg (2000) [16] pointed out that critical DSRs whose yields were extreme points on the boundary of the attainable region had to obey a certain rule (or a certain design equation), restricting the addition policy of side streams along the reactor. Later, Feinberg (2000) [17] also studied the critical
CFSTRs whose effluent compositions lie on the boundary of the attainable region.
Similar to the fact that critical DSRs need to obey a certain side stream addition rule, critical CFSTRs need to have certain residence times whose values can be determined by the structure of the underlying chemical reaction network and its associated kinetics. To find the boundary of the attainable region, many efforts have been tried. But most of the methods were to extend the existing set of compositions by further reaction and mixing, which can be called a “build from the inside” technique. As a new idea compared with this “build from the inside” method, Abraham and Feinberg (2004) proposed a “bound from the outside” technique to confine the attainable region by a set of bounding hyperplanes. It has been shown for the classical van de Vusse network (Feinberg and
6 Ellison, 2001) that this bounding-hyperplane method provided a region that tightly confined all attainable compositions. This region appeared to be very close to the previously found region by the “build from the inside” method.
1.2.1.2 The Superstructure and Mixed-integer Nonlinear Programming Approach
The superstructure idea has been widely adapted in the reactor synthesis area.
Generally speaking, a superstructure refers to a set of all possible reactor configurations from various combinations of limited numbers of known-types of reactors (e.g., CFSTRs and PFRs), with the same resources (e.g., feed streams, catalyst mass) and the same process constraints (e.g., temperature and pressure ranges in the reactor units and applicable environmental constraints). Once the superstructure is generated, a search can be conducted to find the optimal reactor configuration within this superstructure aiming at a certain objective. Normally, the search for the optimal reactor configuration (within the superstructure) is modeled as a mixed-integer nonlinear optimization problem
(MINLP) because the objective and some constraints are nonlinear functions of the variables that include integer-typed variables for the numbers of certain types of reactors used in the optimal reactor configuration. A big drawback of the superstructure idea, however, is that it can only encompass a limited number of reactor configurations from limited types of reactors, although the generated superstructure could be fairly big. In other words, no matter how large the superstructure is, it is always a subset of all possible reactor configurations that could be used since it only considers a limited number of
7 reactors and limited types of reactors. Thus, the optimal reactor configuration for the chosen superstructure can only get close to, but does not give the true optimum for a certain objective.
Achenie and Biegler (1986, 1990) proposed a superstructure based nonlinear programming (NLP) formulation for optimally generating a reactor network including certain number and types of reactors as well as complex mixing patterns. And by heat addition/removal in the network, this formulation also included nonisothermal reactors.
The solution of this model gave the optimal reactor structure, the types of reactors that needed to be used, and the amount of heat that needed to be added for the chosen superstructure, allowing for serial and parallel connections of all proposed reactor units.
Based on the general properties of the attainable region for processes involving only reaction and mixing (Feinberg and Hildebrandt, 1997; Feinberg, 2000 [16, 17]), a reactor module that consisted of CFSTRs, PFRs, and DSRs was used to synthesize the optimal reactor network within a superstructure with respect to a certain objective function
(Lakshmanan and Biegler, 1996 [41]). The whole problem thus was formulated as a mixed integer nonlinear programming problem (MINLP). As an application to the mixed integer nonlinear programming based reactor network synthesis, Lakshmanan, Rooney, and Biegler (1999) studied a vinyl chloride production process. Here, vinyl chloride was produced from ethylene, chlorine, and oxygen and the objective of the integrated reactor network synthesis was to minimize the waste production. As a result, a candidate flowsheet based on this objective was proposed.
8
Kokossis and Floudas (1989) proposed a superstructure-based method for reactor network synthesis under isothermal condition. This superstructure included many possible configurations of CFSTRs, PFRs (which were approximated as a cascade of
CFSTRs), and recycle reactors (which were used to deal with different feeding, recycling, and bypassing strategies). Based on this superstructure, the reactor network synthesis problem became a mixed-integer nonlinear programming problem and its solution gave the optimal reactor network relative to the chosen superstructure for a specific objective.
For more complex reactor networks, the superstructure-based reactor network synthesis problem then became a large-scale mixed-integer nonlinear programming problem
(Kokossis and Floudas, 1990). Similar to the isothermal case, Kokossis and Floudas
(1994) [37] studied the reactor network synthesis problem under nonisothermal condition.
The superstructure idea was still implemented here with an additional consideration of the temperature control of the intercooled/interheated reactors. The whole problem was still formulated as a mixed integer nonlinear programming problem and the solution provided the optimal temperature control, the optimal temperature profile, and the optimal reactor network with certain types and sizes of reactor units (CFSTRs, PFRs) etc.
Kokossis and Floudas (1995) also studied a synthesis problem for a chemical process where it involved a complex reaction network and two liquid phases. Some discontinuous aspects related with the operational or/and structural alternatives of the chemical process were considered and those aspects were incorporated in the mixed-integer nonlinear programming problem as discrete variables. Another application of the superstructure
9 idea in the reactor network synthesis problem could also be found in Schweiger and
Floudas (1999). Later, Pahor et al. (2001) used the superstructure approach in an industrial example of allyl chloride production.
1.2.1.3 Target-based Approach
In this section, target-based approaches will be introduced. Here, target refers to a bound on a chosen performance index and targeting refers to finding such a bound.
By modifying and extending a two-part model proposed by Ng and Rippin (1965),
Achenie and Biegler (1988) modeled the reactor synthesis problem as the following two- compartment problem:
(a) Find a bound on the chosen performance index (targeting) irrespective of reactor types and configurations;
(b) Find a reactor network (with certain reactor types and sizes) that could achieve the above target.
By studying a variable residence time distribution and different mixing profiles between segregation and maximum mixedness, Achenie and Biegler (1988) claimed to successfully generate targets for many isothermal/nonisothermal reacting systems.
Balakrishna and Biegler (1992) [6] developed a simpler and more efficient formulation for finding a bound on a chosen performance index of an isothermal system by considering special cases in Achenie and Biegler (1988) and by applying some of the attainable region concepts of Glasser et al. (1987). A targeting model was developed that provided a rich representation of possible reactor configurations. This model was
10 formulated as a dynamic optimization problem and different solution strategies were presented. Meanwhile, a nonisothermal reaction system was studied where a similar target model was developed (Balakrishna and Biegler, 1992) [7]. The difference between this work and the previous work (for the isothermal reaction system) was that here both reaction and energy synthesis were considered. The targeting model was still formulated as a dynamic optimization problem where the temperature, the feed distribution function etc. acted as control parameters. And the solution strategy was to successively extend the reactor network based on the performance index until no further extension could be carried out. Based on targeting concepts, a targeting approach was applied for the synthesis of waste-minimizing process (Lakshmanan and Biegler, 1994) due to an increasing concern on the development of clean and efficient processes. Tradeoffs between the waste minimization and the profit maximization were incorporated in this approach. To improve the existing methods which lacked robustness and only gave local or near optimal solutions, a new optimization method was proposed for the synthesis of nonisothermal reactor networks (Mehta and Kokossis, 2000). This new method focused on targets instead of network design details.
1.2.1.4 Stochastic Optimization Approach and Other Approaches
Marcoulaki and Kokossis (1995) proposed to use the stochastic optimization technique as a new optimization approach for the reactor network synthesis problem.
Some examples were studied with isothermal reaction kinetics. Manousiouthakis et al.
(2004) utilized the properties of the attainable region for the reactor network synthesis
11 problem to create a so-called shrink-wrap algorithm in order to approximate the true attainable region to an arbitrarily accurate degree. For the reactor network synthesis problem where reactors include CFSTRs and PFRs under isothermal steady-state condition, a new method called Infinite DimEnsionAl State-space (IDEAS) approach was proposed by Burri et al. (2002). With this new method, the reactor network synthesis problem was formulated as an infinite dimensional convex optimization problem.
Kauchali et al. (2002) used a linear programming (LP) technique in the attainable region analysis to provide a method in constructing candidate regions using fully connected
CFSTRs of arbitrary non-negative volumes.
1.2.2 Progress in the Reactor-Separator Synthesis Area
Kokossis and Floudas (1991) proposed a systematic approach for the synthesis of reactor-separator-recycle systems. This approach was based on a superstructure idea that was previously studied by Kokossis and Floudas (1989). Here, the superstructure included a limited number of classical reactors (e.g., CFSTRs and PFRs) and separators, and it encompassed all possible interconnections of those units. The whole problem was formulated as a mixed-integer nonlinear programming problem where the objective included both integer and continuous variables, and the objective was subject to a set of nonlinear constraints. Different objectives were studied for the synthesis problem, including the minimization of the annual plant cost, the maximization of the profit, and the traditional optimization of the performance of the reactor-separator system, i.e., the productivity or selectivity. Again, a big drawback of this superstructure idea is that there
12 are only limited numbers and types of reactors considered in the superstructure. So the optimum found within this superstructure is not the true optimum that is relative to all possible steady-state reactor-separator systems.
Balakrishna and Biegler (1993) proposed a unified approach for the simultaneous synthesis of reaction, energy, and separation systems. This method was based on the idea of the previously studied target-based approach for reactor networks (Balakrishna and
Biegler, 1992 [6, 7]) and the objective in this study was to optimize energy management.
The synthesis problem was formulated as a mixed-integer optimum control problem where the control profiles included the temperature profile, the separation profile, and the residence time distribution defined for the network. Also, costs of separations were embedded in the model.
McGregor, Glasser, and Hildebrandt (1997) applied the attainable region idea in a binary component system where a more volatile reactant was converted to a less volatile product in a reaction. A number of systems were examined including a reactor-distillation system, a reactor-stripper-recycle system, and two catalytic distillation systems where one was with recycle and the other was at total reflux. Each system was optimized according to a simple linear cost function. It was found that the optimal reactor- distillation-recycle system was the most economical one among all studied systems. As a revisit to this paper where no cost was considered in finding the optimal reflux ratio, the effect of including the cost for the processes used to control the reflux ratio on the operation conditions and the distillation column structure was studied by Kauchali et al.
(2000).
13 Shah and Kokossis (1997) proposed a novel methodology called conceptual programming. Using only some basic information, this method focused on a chosen performance index of chemical processes based on which a few promising candidate processes were screened out to be further tested by more rigorous models. Some complex reactor-separator systems were studied by this new method that proved to be better than some previously studied technologies such as pure mathematical programming.
Another novel method for the synthesis of reaction-separation systems was proposed by Linke et al. (2000). It combined previous achievements in multiphase reactor network synthesis (Mehta and Kokossis, 1996, 1997, 2000) and the superstructure idea (Kokossis and Floudas, 1989) to solve problems in both isothermal and nonisothermal systems. Later, Linke and Kokossis (2002) proposed a systematic decision-making strategy for optimal multiphase reaction and separation system design.
The synthesis scheme utilized a superstructure that consisted of two generic units: the reactor/mass exchange units and the separation task units. This scheme supported the decision making in both early and late process design stages and it could incorporate the design information gained at the early stage into the superstructure formulation at a later stage. Following this technology, Linke and Kokossis (2003) [46] used the stochastic optimization technology in the synthesis of reaction/separation systems. Two stochastic optimization algorithms, the simulated annealing search and the tabu search, were presented.
14 As a powerful tool for the reactor-separator synthesis problem, a theorem called the CFSTR Equivalence Principle for Reactor-Separator Systems was proposed in
Feinberg and Ellison (2001). Using this conceptual tool, Feinberg and Ellison (2001) studied the general kinetic bounds on productivity and selectivity in reactor-separator systems of arbitrary design. A reaction network called the van de Vusse network was presented as an example to see how the CFSTR Equivalence Principle could be used to find the maximum productivity of a certain desirable species given the feed stream information, the maximum capacity of the reactor units, the underlying reaction network with its associated kinetics, and other constraints including temperature and pressure ranges in the reactor units, and an applicable environmental constraint. The van de Vusse example is, however, a very small problem and the result of this example was found in an analytical way. To find the kinetic bounds on productivity for larger and more complicated chemical reaction networks will require more sophisticated methods, which is one of the main objectives of this thesis.
Although a lot of work has been done in both pure reactor synthesis problem and reactor-separator synthesis problem, some mentioned in the previous paragraphs, there are still problems that await us to solve and to implement, one of which is of the primary concern in this dissertation. The following section 1.3 will give a detailed description of the problem.
15
1.3 Finding Kinetic Bounds on Productivity in Reactor-Separator Systems of
Arbitrary Design
1.3.1 Addressing the Problem
The problem that this thesis is trying to solve is different from any of the problems in the previous work. What distinguishes the work in this thesis from all other work is that this thesis is trying to find the kinetic bounds of all possible effluents for a chemical reaction network with given kinetics over all possible steady-state reactor- separator designs. Speaking of all possible designs, we are looking at not only the traditional types of reactors (e.g., CFSTR and PFR), but also other types of reactors that could be even beyond our imagination. None of the previous work has tried this before because they only consider limited types of reactors, ones that we know.
By “all possible” steady-state reactor-separator systems, we need to clarify certain things before we move any further. Just for practical purposes, should we consider extremely high or extremely low temperatures in the reactor units? Should we consider extremely high or extremely low pressures in the reactors? And speaking of the reactor units, is it reasonable to consider infinitely large volumes or infinitely large amount of catalysts used in the reactors? With these practical constraints in mind, we can see that by
“all possible” steady-state reactor-separator systems, we really mean the set of all steady- state reactor-separator systems in which reactors are not operated under unreasonable temperatures or pressures and the total capacity of the reactors is not infinitely large.
16 (Here, we use the total capacity of a reactor to indicate how “big” a reactor is. It could be thought of as the actual physical size of the reactor, or it could be referred to as the total mass of the catalyst used in the reactor.)
Meanwhile, in almost all chemical processes, environmental impact is one of the concerns besides the productivity and the economics of the process. While trying to increase the production rate of a certain desirable chemical species, the production rate of some environmentally unfavorable species needs to be controlled within a certain limit.
This adds another constraint in addition to the previous constraints including the temperature-pressure ranges in the reactor units and the available reactor capacity.
With all these constraints born in our mind, the question that this thesis tries to answer is stated as follows: Given a system of chemical reactions (or a chemical reaction network) with its associated kinetics, given a specified feed stream, what is the highest production rate of a certain desirable chemical species over all possible steady-state reactor-separator systems of arbitrary design that are consistent with all constraints including the temperature-pressure ranges in the reactor units, the available reactor capacity, and applicable environmental concerns? In other words, what are the kinetic bounds on the productivity in reactor-separator systems of arbitrary design?
The importance of addressing this problem is that the answer to this question provides us an important guideline in designing chemical processes. The reason is that once we know the maximum productivity of a certain desired chemical species for a chemical reaction network over all possible steady-state reactor-separator systems of arbitrary design, we can compare this maximum with the productivity of a candidate
17 design. If the candidate design is economical and at the same time yields a production rate close to that maximum, it means that the candidate design is doing very well. On the other hand, if the candidate design is not very economical and gives a production rate far from the maximum, we might want to take a look at the candidate design and see if we can come up with something better that can yield more of the desired species and cost less at the same time. In other words, knowledge of the maximum production rate of a desired chemical species over all possible steady-state reactor-separator designs provides an important benchmark against which existing designs can be measured.
One thing to keep in mind is that it is not the purpose of this thesis to come up with an economical design for a chemical process. The purpose of this thesis is to tell the best possible yield of a desired chemical species over all possible steady-state reactor- separator systems given the underlying chemical reaction network with its associated kinetics, given the feed stream information, and given the constraints like the temperature and pressure ranges in the reactor units, the available reactor capacity, and applicable environmental constraints. To achieve that best possible yield may require a very costly process. However, knowing the upper limit of how much the desired species can be produced gives designers a picture of where their own designs stand and therefore helps them to make decisions about whether or not they need to consider improving their own designs.
To get a better picture of the problem, let us suppose that we have a chemical process operated at steady state with unspecified design. We can envision this process as a “black box” into which certain feed streams enter and from which certain effluent
18 steams exit, as shown in Figure 1.1. There are N species A1, A2, …, AN involved in the
00 0 reactions. At the “entrance” of the chemical process, we use symbols M12,MM ,..., N to
0 represent the total molar feed rates of species A1, A2, …, AN. (That is, M i ≥ 0 represents the total rate at which moles of species Ai are carried into the process by means of all feed streams.) For convenience, we can define a molar feed rate vector M0 such that
0000 M = [MM12 , ,..., MN ] . Similarly, at the “exit” of the chemical process, we use M1,
M2, …, MN to represent the total molar effluent rates of species A1, A2, …, AN and define a molar effluent rate vector M such that M = [M1, M2, …, MN]. (That is, Mi ≥ 0 represents the total rate at which moles of species Ai are carried out from the process by means of all effluent streams.) Suppose that we know the reactions and the associated kinetics, and we want to know what the highest production rate is for a certain desired chemical species, say, A1, over all possible steady-state reactor-separator systems of arbitrary design, but consistent with broadly specified constraints.
Figure 1.1: A chemical process with an unspecified design.
19 Looking at Figure 1.1, we can see that even if we constrain our thinking to reactor-separator systems under certain temperature, pressure and reactor capacity constraints, there is still an enormously wide range of possible reactor-separator configurations which can be employed in the chemical process. Figure 1.2 shows a small sample of feasible reactor-separator systems.
Figure 1.2: Varieties of reactor-separator systems.
20
Figure 1.2 tells us two things. First, we can consider all kinds of combinations of reactors that we already know. That is, by combining different reactors in different ways we can have hundreds of designs. For example, we can have one single CFSTR, or one single PFR. We can have two CFSTR in series, or we can put them in parallel. Second, besides all the designs we can imagine, there could be something that is even beyond our imagination. And for those unknown reactors, we can continue our game of combining reactors and separators in endless ways. So from the figure, it can be seen that to find out the maximum productivity by examining all possible reactor-separator configurations seems extremely difficult and almost impossible, because the possible combinations of reactors and separators (including knowns and unknowns) are just countless.
1.3.2 An Important Conceptual Tool The CFSTR Equivalence Principle for
Reactor-Separator Systems
Facing such a difficulty described as above, does that mean that there is no way of knowing the maximum productivity (of a certain desired species) over all possible steady-state reactor-separator configurations? The answer is no. A very important and powerful conceptual tool for leading us out of this dilemma is a theory that creates a
“bypass” to find that maximum. This theory is called the CFSTR Equivalence Principle for Reactor-Separator Systems (Feinberg and Ellison, 2001). By this theory, we no longer need to look at all reactor-separator configurations while at the same time we are
21 able to get the maximum productivity of a certain desired species for any chemical reaction network. We will give details of the theory and explain why this theory can help us to find the maximum productivity in Chapter 2.
Before we get into the details of the theory, there is something to think about.
What, really, gives the breadth of reactor-separator configurations? Let us imagine the situation where no reactor is involved in the process except that only separators are allowed. We can use different types of separators and combine all these separators in whatever way we want. But no matter how complicated the separation system will be, the molar effluent rate vector M will always be the same as the molar feed rate vector M0 because no reaction is involved in the process, so there should be no change of the mole numbers of all species. From this, we can clearly see that it is the variety of reactors that gives the breadth of possible molar effluent rate vectors. And it is mainly the variety of reactors that makes the problem of finding the maximum productivity of the desired species so difficult. Since reactors play a somewhat more important role than separators, we thereby make an important assumption of arbitrary separation ability of the separation system in the CFSTR Equivalence Principle. To make such an assumption is not to make the problem appear simpler. Rather, doing this is to ensure that the bounds of the attainable effluents are as broad as possible, because what can be attained with the existence of arbitrary separations also includes what can be attained from limited separation ability. Therefore, the maximum productivity (of a certain desired species) obtained from the arbitrary separation assumption is the upper limit of all feasible
22
productivity that can be found. (We can also regard the kinetic bounds of what is attainable from the assumption of arbitrary separation ability the theoretical limit of the attainable effluents.)
There is a small issue to address about the arbitrary separation ability. As we have seen in the literature review of some previous work, one branch of process synthesis problem that is the pure reactor synthesis problem studies processes involving only reaction and mixing. However, what we are interested in here is the reactor-separator synthesis problem where separators are not only allowed, they are also presumed to have arbitrary separation ability. It is obvious to see that the difference between these two problems is that one does not use separation at all, and the other uses separation to an arbitrary degree. It might be wondered why we should study the pure reactor synthesis problem when separation is widely used in chemical processes. The reason is that sometimes separation could be very expensive so it is worth knowing the best that can be achieved if no separation is involved. So while separation is commonly used, it does not mean that the pure reactor synthesis problem loses its importance. However, in this thesis we only study the reactor-separator synthesis problem for a couple of reasons. First, since separation is a common practice, why not allow it in our study? Second, what we can get allowing arbitrary separation ability will also include what we can get if no separation is used. In other words, the bound of what is attainable for the reactor-separator synthesis problem is also a bound of what is attainable for the pure reactor synthesis problem. It
23 may not be a sharp bound for the pure reactor synthesis problem. Nevertheless, it gives a coarse estimation and provides some idea at the initial stage of the pure reactor synthesis problem.
1.3.3 Beyond the Theory One Step Further
We have seen in Section 1.3.1 that the problem that the thesis is trying to solve is to find kinetic bounds on productivity in reactor-separator systems of arbitrary design.
And Section 1.3.2 tells us that the CFSTR Equivalence Principle for Reactor-Separator
Systems is a very important tool to solve this problem. Although the theory provides an idea for solving the problem in general, how to get the actual answers (that is, values of the maximum production rates of the desired species) for individual cases hasn’t been answered. So the most important task of the thesis is to use the idea of the CFSTR
Equivalence Principle to come up with a brand-new practical computational tool for calculating the kinetic bounds on productivity in reactor-separator systems of arbitrary design. It is a big step from having a conceptual tool to having a computational tool that everybody can use (even if the user doesn’t understand the background theory), either for educational purposes or for application to real chemical processes. The importance of having such a computational tool is not only that it can calculate the maximum production rate (of certain desired species) for each individual case, but also that the results from the computational tool can sometimes lead to other important questions and
24 thinking which may not be seen otherwise. In other words, the computational tool is not just for computational purposes, but also has a deeper use in exploring more interesting questions and in making further progresses in the process synthesis area.
1.4 Summary of Chapter 1
A big question in the process synthesis area in chemical engineering has been introduced. That is, given a system of chemical reactions with its associated kinetics, the feed stream information, the available reactor capacity, the allowable temperature- pressure ranges within the reactor units, and applicable environmental constraints, what is the highest production rate of a certain desirable species over all possible steady-state reactor-separator systems of arbitrary design consistent with all the constraints?
To answer this question, there is an available conceptual tool (or theory) that is called the CFSTR Equivalence Principle for Reactor-Separator Systems. However, by the theory itself, the answer to this question can only stay on a conceptual level. Trying to push that answer from the conceptual level to a practical level is what the thesis is mainly about. And the ultimate goal of the thesis is to utilize the idea in the CFSTR Equivalence
Principle to come up with a convenient computational tool that can be used to calculate the kinetic bounds on productivity in reactor-separator systems of arbitrary design.
Having stated the main objective of this thesis, it is time to briefly look at how the contents of this thesis are organized. Chapter 2 will give a detailed explanation of the
CFSTR Equivalence Principle for Reactor-Separator Systems, which is the primary
25 conceptual tool in this thesis for calculating the kinetic bounds on productivity in reactor- separator systems of arbitrary design. From this theory, characterization of attainable effluents will be given in a proposition. After that, several practical implementation methods of finding the kinetic bounds of attainable effluents will be given in Chapter 3.
Based on these methods, a software package is developed to serve as a computational tool to get the productivity limit in the reactor-separator synthesis problem; this will be demonstrated in Chapter 4. By using this software package, we will look at some concrete examples in Chapter 5. At the end of this thesis, conclusions will be given and some future research work will be proposed in Chapter 6.
26
CHAPTER 2
THE CFSTR EQUIVALENCE PRINCIPLE FOR REACTOR-
SEPARATOR SYSTEMS
We already mentioned in Chapter 1 that to find the kinetic bounds on productivity of a certain desired species for a chemical reaction network over all possible steady-state reactor-separator systems is difficult and almost impossible if the way to do this is to enumerate and examine all possible reactor-separator configurations. And we also mentioned that there is an available conceptual tool that can be used to solve this problem.
The conceptual tool (or theory) is called the CFSTR Equivalence Principle for Reactor-
Separator Systems. Basically, what the theory says is that it is not necessary to look at all specific designs to find out the kinetic bounds of the attainable effluent rates given a chemical reaction network with its associated kinetics, the feed stream information, and other information including the temperature-pressure ranges in the reactor units, the available reactor capacity, and the environmental constraints, if there are any. Rather, the theory proves that, for the sole purpose of getting that kinetic bound, it suffices to look at
27 a surprisingly small class of reactor configurations. What that small class of reactor configurations is, why and how it can give the kinetic bounds on the productivity will be explained later in this chapter in detail.
Before explaining the theory, there are some preliminaries that need to be introduced.
2.1 Some Preliminaries
2.1.1 A Few Mathematical Preliminaries
Here we introduce some basic mathematical concepts in preparation of the introduction of the CFSTR Equivalence Principle.
Linear Independence and Linear Subspace:
The usual vector space of N-tuples of real numbers is denoted as RN. And the set
N N of vectors in R with entirely nonnegative components is denoted as R+ . The standard
N basis for R is denoted as e1, e2, …, eN. That is,
e1 = [1, 0, 0, …, 0], e2 = [0, 1, 0, …, 0], …, eN = [0, 0, 0, …, 1].
N Suppose we have vectors v1, v2, …, vk in R . We say that v1, v2, …, vk are linearly dependent if there is a set of real numbers {α1, α2, …, αk}, not all zero, such that
α1 v1 + α2 v2 +… + αk vk = 0.
Otherwise, v1, v2, …, vk are linearly independent.
28 A set Ω of vectors in RN is said to have rank r (r is an integer) if Ω contains r linearly independent vectors but no set of r + 1 linearly independent vectors.
A nonempty set U of vectors in RN is a linear subspace if, for all vectors u and u’ in U and all real numbers α and α’, the vector αu + α’u’ is also a member of U. When U is a linear subspace, the rank of U is called the dimension of U.
2.1.2 Stoichiometric and Kinetic Preliminaries
Reaction Vectors
For a network of chemical reactions having N species, we associate with each reaction a vector in RN as indicated in the following example:
Example: Suppose we have a chemical reaction network which has five species
A1, A2, …, A5. The network is shown as follows:
AA12+→ A 3
2AA34→ AA14+ A 5
5 For reaction AA12+→ A 3 we associate with it a reaction vector in R which is
e3 – e1 – e2 = [-1, -1, 1, 0, 0]; for reaction 2A34→ A we associate with it a reaction
vector e4 – 2e3 = [0, 0, -2, 1, 0]; for reaction A14+AA→ 5 we associate with it a
reaction vector e5 – e1 – e4 = [-1, 0, 0, -1, 1]; and for reaction AAA514→+ we associate with it a reaction vector e1 + e4 – e5 = [1, 0, 0, 1, -1]. 29
Rank of a Chemical Reaction Network
The rank of a chemical reaction network is the rank of the network’s associated reaction vectors.
For the reaction network in the above example
AA12+→ A 3
2AA34→ AA14+ A 5 the set of associated reaction vectors in R5 is:
{[-1, -1, 1, 0, 0],
[0, 0, -2, 1, 0],
[-1, 0, 0, -1, 1],
[1, 0, 0, 1, -1]}.
The first three vectors [-1, -1, 1, 0, 0], [0, 0, -2, 1, 0], [-1, 0, 0, -1, 1] are linearly independent. But all the four vectors are linearly dependent. (The last vector [1, 0, 0, 1, -1] can be obtained by multiplying the third vector [-1, 0, 0, -1, 1] by -1.) So the rank of the set of these four vectors is 3, which means that the rank of the chemical reaction network is 3.
In another chemical reaction network shown as follows:
AAA→ 123
AA12+→2 A 3 there are four reaction vectors in R3 listed as follows:
30 {[-1, 1, 0],
[0, -1, 1],
[0, 1, -1],
[-1, -1, 2]}.
The first two vectors [-1, 1, 0], [0, -1, 1] are linearly independent and the other two vectors are linear combinations of the first two vectors. So the rank of this chemical reaction network is 2.
Species Formation Rate Vector
For a chemical reaction network having N species, we define a composition vector c in RN to be:
c = [c1, c2, …, cN]
th where ci denotes the (molar) concentration of the i species, i = 1, …, N.
For a local composition vector c and a local temperature T, we define a
(volumetric) species formation rate vector r(c, T) in RN to be:
r(c, T) = [r1(c, T), r2(c, T), …, rN(c, T)]
th where ri(c, T) denotes the molar species formation rate of the i species (i = 1, …, N) at composition c and temperature T.
For example, for the following reaction network with mass action kinetics
k1()T AA12+→ A 3