Surrogate Model Generation Using the Concepts of Self-Optimizing Control
|
|
- Bryan Stewart Lambert
- 5 years ago
- Views:
Transcription
1 Surrogate Model Generation Using the Concepts of Self-Optimizing Control Julian Straus, Sigurd Skogestad Department of Chemical Engineering Norwegian University of Science and Technology Trondheim, Norway {julian.straus, Abstract This paper presents the application of selfoptimizing control in the concept of surrogate model generation. Surrogate model generation generally has problems with a large number of independent variables resulting in a large sampling case. Utilizing the concept of self-optimizing control and keeping measurements constant, it is on the one hand possible to reduce the number of independent variables resulting in a reduced sample space. On the other hand, this allows to map a close-tooptimal response surface. As the surrogate model is subsequently used for optimization, this is as well advantageous; regions in which the system is suboptimal are not mapped this way. The proposed method is studied using an ammonia reactor which for disturbances experiences limit-cycle behaviour and/or reactor extinction. Using SOC, it is possible to reduce the number of manipulated variables by three and map the optimal response surface. I. INTRODUCTION Surrogate models are an emerging field and may be utilized for the optimization of integrated flowsheets [1]. Using a sequential-modular software like Aspen Plus, Aspen Hysys, SimSci PRO/II, or UniSim Design Suite, problems may arise in the convergence, whereas equation-orientated solver may be difficult for complicated models. It is possible to simplify the optimization through the incorporation of recycles within the surrogate models. The calculation of surrogate models struggles however if the number of independent variables (n u ) is high. This is generally the case, if the surrogate model i is designed for the latter combination with other surrogate models j. In this case, the connection variables x i, j may increase n u. This so-called curse of dimensionality says that the number of sampling points n p grows exponentially with the number of independent variables n u. This may result in an excessive number of points one has to sample. One alternative is the application of PLS regression to define new latent variables u with dim(u ) < dim(u) [2]. As a second alternative, we proposed to use the concepts of self-optimizing control [3] to reduce n u [1]. This allows the mapping of the region we are actually interested in and will hence be further investigated in this paper. This paper is structured as follows; Section II introduces the problem in the context of surrogate model generation for a submodels. Section III explains how self-optimizing control is applied in this context and the properties of self-optimizing control. Section IV first introduces the utilized case study and then shows results from the application of self-optimizing control in surrogate model generation. Section V discusses the applicability of the proposed procedure and addresses limitations and problems if self-optimizing control is utilized. II. PROBLEM STATEMENT Consider a large scale process to be optimized. This is generally difficult and the process is hence split into several submodels. Each of these submodels should be reformulated as a surrogate model to simplify the optimization. An optimization problem as shown below can then be defined for some of the submodels i min J i (x i,u i,d i ) x i,u i s.t. 0 = g i (x i,u i,d i ) 0 h i (x i,u i,d i ) where J i : R n x R n u R n d R describes the economic objective of the system, g i : R n x R n u R n d R n x the plant model and h i : R n x R n u R n d R n h i the operational inequality constraints. The inlet connection variables, x i,in, are hereby considered to be the disturbances d i R n d. The manipulated variables of a submodel i, u i R n u,mv, are the real degrees of freedom for the submodel. However, we consider the total number of independent variables, n u = n u,mv + n d, for surrogate model generation. The sampling space is then given by bounds as (1) d min d d max (2) u min u u max (3) and the sampling is performed using e.g. Latin hypercube sampling or regular grid sampling. III. PROCEDURE A direct possibility to reduce the number of independent variables by n u,mv is given by fitting a surrogate model for the optimal response surface. If n u,mv is relatively large, this would reduce the number of points which one has to sample. This approach would require the solution of n p nonlinear problems. A further advantage of this approach is that the surrogate model is sampled in the direction interesting for the subsequent optimization. An analogy to this approach is the creation of road maps. Theoretically, it is possible to include
2 in a map every single small road and even hiking paths. This would however make the map complicated. Hence, maps only feature the roads that are useful for the desired application. Alternatively, self-optimizing control can be applied to the submodels and additional n c equality constraints given by 0 = c i H i y i (4) with c i R n c being the SOC variables and y i R n y being chosen measurements of the system are added to the nonlinear system given by g i in Eq. (1) [4]. The matrix H i is the optimal selection matrix and can be calculated in different ways [5]. The advantage of using self-optimizing control in this context can be seen in the reduced number of optimization problems one has to solve compared to mapping the optimal response surface. These correspond to n d + 1 nonlinear problems with additional n u systems of nonlinear equations in the definition of the SOC selection matrix according to Eq. (4). The sampling then consists of solving n p nonlinear systems of equations. A. Calculation of the selection matrix H The optimal selection matrix H can be calculated using the nullspace method [6] or the exact-local method [7]. It is given by the solution to the following optimization problem min H HY F s.t. HG y = J 1/2 uu with G y R n y n u,mv representing the measurement gain matrix with respect to the input u. Y is given by Y = [ ] FW d W n y (6) The optimal sensitivity matrix F = y opt d F = ( G y J 1 uu J ud G y d (5) can be calculated as ) (7) using the disturbance gain G y d Rn y n d, the hessian of the cost function J uu R n u,mv n u,mv, and J uu R n u,mv n d the second order derivative of J with respect u and d. W d and W n are the disturbance and measurement scaling matrices given by d = W d d ; n y = W n yn y (8) in which the vectors d and n y are assumed to satisfy [ y ] d 2 1 (9) n This scaling is necessary as this allows us to define selection matrices which minimizes the loss if disturbances occur. It is important to note for the scaling matrices that the 2- norm is used in the scaling in Eq. (9). This implies, that all disturbances and measurements may not be at their upper or lower limit simultaneously. In the case of control, this seems reasonable and a detailed discussion for using the 2 -norm is given by Halvorsen et al. [4]. However, this is not the case, if we want to use the self-optimizing control variables for variable reduction in surrogate model definition. We actually want to sample these so-called corner points to avoid extrapolation. This could be circumvented by multiplying the W d by nd and would make a difference, if measurement noise is considered. We propose to set the disturbance scaling matrix W d in the case of surrogate model generation to W d = diag(max(d d min,d max d)) (10) and the measurement noise scaling matrix to W n y = 1 k exp diag(1) (11) with 1 being a vector of ones with length n y. The parameter k exp can be chosen arbitrarily as the measurements will not have any noise in the case of surrogate model sampling. However, two necessities arise for the parameter k exp k exp is large enough so that YY T is nonsingular; k exp should be small compared to the entries of W d to reduce the effect of measurement noise in the calculation of the selection matrix H; The solution to this problem in its simplified version is given by [8]: H T = ( YY T) 1 G y (12) One could argue that the nullspace method [6] can be used as an alternative as measurement noise is not existing and the number of measurement is not important. Hence, the dynamic properties of this specific control structure do not play a role and it is possible to increase the number of measurements This results in a reduced loss as shown in [8]. Furthermore, it is possible to extent the measurements to states, which are generally not considered as they are hard to measure, e.g. concentrations. The application of the nullspace method would require that n y n d + n u,mv (13) which may not be possible to satisfy for submodels with a lot of manipulated variables and connection variables. Additionally, problems may arose in the solving of the flowsheet, if a sequential-modular flowsheet solver is used. In this situation, the setpoints to the original manipulated variables u is obtained in an iterative manner which may slow down the solution of each sampling point. This may lead to problems especially in the case of interacting manipulated variables u. Hence, it is generally advisable to choose only a limited number of measurements. B. Measurement Selection In order to select the an optimal subset of measurements n y,sel, Yelchuru and Skogestad [8] developed a mixed integer quadratic programming approach. It requires the reformulation of the problem given in Eq. (5) in vectorized form: min h δ σ δ h T δ F δ h δ s.t. G y δ T h δ = j δ n y j=1 σ δ, j = n y,sel (14)
3 where σ j {0,1} with j = 1 n y are binary variables to indicate, whether measurements are used in the selection matrix and the optimization matrix. The quadratic cost term is given by F δ = Y δ Y T δ (15) in T in,3 u 3 T in,2 u 2 T in,1 u 1 and is block diagonal. The same holds true for G y T δ whereas h δ and j δ are a vectorized form of H and J uu respectively. Further constraints have to be imposed on h δ to guarantee that h jk = 0 for σ k = 0 and input u j and measurement y k. In this problem, the big-m approach is chosen. This results in bounds for the entries in the selection matrix H given by m h 1k m m. σ h 2k k. m. σ k σ k, k 1,2,...,n y (16) m h nu k m Bed 3 T 30 w 30 T 20 out Bed 2 Bed 1 T 10 For a detailed description and derivation, the reader is referred to [8]. Unfortunately, this approach does not handle structural zeros in the selection matrix H. That is, certain measurements should be candidates for certain manipulated variables and would lead to e.g. block diagonal selection matrices. This may then lead to complicated adjustments in the flowsheet solver. Hence, we propose to calculated optimal selection matrices H i R 1 n y k for input u k using close measurements and not the overall measurement set. This does guarantee the optimal measurement combinations and may lead to cases, where problems may arise. However, it is not possible to generalize when problems may occur and when not. Unfortunately, there are no simple methods for solving this problem in a convincing way [5]. IV. CASE STUDY - AMMONIA REACTOR The case study for the investigation of this procedure is given by a heat-integrated ammonia reactor. The investigated reactor submodel is shown in Figure 1. This reactor was previously used in stability analysis [9] as well as the feasibility of the application of economic NMPC [10] for the given system. A detailed model description can be found in [10]. The subscript i is dropped for this submodel in the following to improve the understanding. At the optimal operation point, small disturbances lead to limit-cycle behaviour and/or reactor extinction [10]. Hence, varying the manipulated variables u individually results on the one hand in creating a response surface including regions, in which reactor extinction is present. In regions where limit cycle behaviour occurs, it is on the other hand not possible to define a steady state for the system. Hence, the response surface is complicated and it is necessary to sample a lot of points to achieve a decent surrogate model. A. Model Adaptation In order to increase the applicability of the resulting surrogate model, the hydrogen to nitrogen molar ratio is not Fig. 1. Heat-integrated 3 bed reactor system of the ammonia synthesis gas loop. considered to be fixed anymore. Instead, the ratio of hydrogen to nitrogen given by R H2 /N 2, j = ṅh 2, j ṅ N2, j (17) in each reaction section j is introduced. This results in 3n additional algebraic constraints given by 0 = R H2 /N 2, j ṅh 2, j 1 + r H2, jm cat, j /M H2 ṅ N2, j 1 + r N2, jm cat, j /M N2 (18) in which M i is the respective molar mass and r i, j the reaction rate in [kg i/kg cat h]. Alternatively, it is possible to use the respective mass balance of hydrogen or nitrogen. It seems however useful, to use a fifth inlet variable, whose influence is directly known to process operators and which offers a direct relation to the processes before the reactor. The cost function for the optimization problem shown in Eq. (1) is given by the (mass) extent of reaction ξ ξ = ṁ in (w NH3,30 w NH3,in) (19) whereas the equality constraints are given by the ammonia mass balance and the temperature balance described in [10] with addition of Eq. (18) for each CSTR j in the CSTR cascade. For this system, the independent variables are given by the three split ratios as well as the inlet to the system u = [ u 1 u 2 u 3 ] T (20) d = [ ṁ in p in T in w NH3,in R H2 /N 2,in] T (21) This implies that n u for this submodel corresponds to 8 of which 3 are real manipulated variables. The other 5 variables are input connection variables. Their respective bounds are
4 TABLE I BOUNDS AND UNITS FOR THE CONNECTION VARIABLES. ṁ in p in T in w NH3,in R H2/N2,in [kg/s] [bar] [ C] [wt.%] [-] Lower Bound Nominal Point Upper Bound given in Table I. Two output variables have to be fitted due to the assumption of constant pressure and a single outlet stream. These correspond to the outlet temperature T out and the (mass) extent of reaction ξ. The outlet ratio R H2 /N 2,out can hereby be calculated through the respective outlet molar flows ṅ i,out which in turn are calculated using ξ. This furthermore guarantees mass conservation in the resulting surrogate model. The system was modelled using CasADi [11] and optimized using IPOPT [12]. The number of CSTRs in each bed is given by n = 10. B. Application of SOC As n u,mv = 3, three SOC variables have to be calculated. According to Skogestad [3], a self-optimizing control variable has to satisfy among other requirements 1) the optimal value of the controlled variable should be insensitive to disturbances, which is equivalent to the entries in the scaled sensitivity matrix for the respective measurement is small; 2) the gain from the input i to the controlled variable c i should be large. This corresponds to a flat optimum with respect to c i ; 3) the controlled variables c should not be closely related; The chosen measurements are give by a local analysis for each reactor bed to reduce the correlation in-between the measurements. This does not imply that the measurements and their combination are as well optimal considering the overall process. The MIQP approach as proposed by Yelchuru and Skogestad [8] is applied hence individually for each bed resulting in a block diagonal matrix given by H = H H 2 0 (22) 0 0 H 2 In order to have a small number of measurements, n y,sel = 1 and n y,sel = 2 are used and compared to the intuitive control structure, where the input and output temperatures of the respective beds are used. The investigated measurements for the MIQP approach are given by y 1 = [ T in,1 T T 1:10] T for Bed 1, u 1 y 2 = [ T in,2 T11:20] T T for Bed 2, u 2 (23) y 3 = [ ] T in,3 T T T 21:30 for Bed 3, u 3 The scaling matrix W d according to Eq. (10) and Table I is given by W d = diag ([ ]) (24) Fig. 2. Infinity norm for the row k of the scaled sensitivity matrix FW d and gain G y,k j as a function of the measurement y comb,k in the system and input u j. whereas the parameter k exp in the calculation of W n y is given by k exp = 3. Based on the requirements for SOC variables, it may be interesting to look at the sensitivity of single measurements with respect to disturbances as well as the gain. The infinity norm of the rows in the scaled sensitivity matrix FW d as well as the gain G y, jk of each temperature k with respect to the corresponding input u j are given in Figure 2. The temperature coordinate k is given by entry k in the vector y comb = [ ] T. y T 1 y T 2 y T 3 The application of the infinity norm instead of the l 1 or l 2 norm can be reasoned by the fact, that we are interested in the worst case scenario. From this Figure, it would make sense to use as a single measurement a temperature close to the end of each bed. This would correspond to a maximized gain and minimized scaled sensitivity matrix and is indeed the solution of the optimization problem (14) with n y,sel = 1 H QP1,1 = T 9 H QP1,2 = T 18 H QP1,3 = T 25 (25) with m = 100 in the big-m approach of Eq. (16). The solution to the optimization problem (14) with n y,sel = 2 is given in Table II and unfortunately cannot be explained by Figure 2. Similar to the results reported by Yelchuru and Skogestad [8], the chosen measurements change depending on the chosen number of measurements n y,sel. C. Fitting of the Surrogate Model The surrogate models are cubic B-splines fitted through the application of the SPLINTER library [13] which requires a regular grid. The regular grid is given by four points for each of the varied variable d resulting in n p = 1024 sampling points. The advantage of using B-splines of order two or higher is given by the continuity of the first derivative of the surrogate model. This gives advantages for the subsequent optimization in which the surrogate model should be used. If self-optimizing control is not used for variable reduction,
5 TABLE II OPTIMAL MEASUREMENT SUBSET FOR EACH INPUT AND THE CORRESPONDING OPTIMAL SELECTION MATRIX H i WITH n y,sel = 2 (MIQP) AS WELL AS THE OPTIMAL SELECTION MATRIX FOR A FIXED SELECTION (IN-OUT). MIQP In-Out Chosen Variables Selection Matrix [ ] H QP2,1 T 4, T [ ] H QP2,2 T in,2, T [ ] H QP2,3 T 28, T H IO,1 T in,1, T 10 [ ] H IO,2 T in,2, T 20 [ ] H IO,3 T in,3, T 30 [ ] this would correspond to 4 8 = sampling points. Alternatively, other surrogate model structures like Kriging [14] or the Alamo approach [15] could be used. D. Results The resulting surrogate models for the outlet temperature T out and ξ were evaluated using 5000 randomly sampled validation points. These validation points are the optimal response surface for this model. This implies, that the surrogate model may theoretically give perfect fit for self-optimizing control model structure and it would not be seen in the results. However, this is only of minor interest as the aim of the surrogate model is to utilize it in further optimization. The relative error ε for the dependent variable is subsequently calculated with the optimal response surface. In order to compare the different methods, the maximum absolute relative error max ε and the mean absolute relative error ε are calculated. The results of the three different combination matrices can be found in Table III. Based on the presented results, we can see that arbitrarily chosen measurements do not necessarily result in a good surrogate model fit. Through application of the MIQP approach in selecting the optimal measurement variables, it is possible to reduce the maximum and mean error with respect to the optimal response surface by more than a factor of two for both the outlet temperature and the extent of reaction. Increasing the number of measurements in the MIQP approach to 2 results in a decrease of one order of magnitude with respect to 1 measurement. Varying the inlet variables and the manipulated variables as an alternative to utilizing self-optimizing control is for this case study not advisable. This leads in the case of all inlet variables being at their lower bound (see Table I) and the manipulated variables at their nominal optimum to reactor extinction. Contrary, if all inlet and manipulated variables are at their nominal value, the reactor is at its optimum. Hence, TABLE III RESULTS FOR THE THREE SOC VARIABLES USING DIFFERENT SELECTION MATRICES H. H definition Extent of Reaction ξ Outlet Temperature T out max ε ε max ε ε In-Out % % % % MIQP % % % % MIQP % % % % x Fig. 3. Outlet temperature of Bed 3 with a pressure drop of p in = 15 bar at t = 10 min with a constant input u at the optimal point. including the split ratios in the independent variable space would require the mapping of regions in which the reactor is extinct as well as crossing the limit-cycle region in which no upper steady-state solution exists. This region is exemplified in Figure 3 where the inlet pressure is at its lower bound and the other disturbances at their nominal value. We can directly see, that the system is in limit-cycle behaviour and it is not possible to define a steady-state value for this operating point. These regions are not important for the subsequent optimization, and hence, should not be sampled. V. DISCUSSION The proposed utilization self-optimizing control to map the optimal response surface can be a promising new method in the generation of surrogate models. The main advantages are given by 1) a reduced number of sampling points compared to sampling the total number of independent variables; 2) a response surface which is close to the optimal response surface with a reduced number of optimization problems to be solved; The reduced number of sampling points is directly visible through the definition of setpoints of controlled variables. The latter allows as well to sample only regions which we are interested in and neglect regions, which will not be approached in the practical application. In the proposed case study, it is in fact not possible to use the splitratios as independent variables and a variable transformation would be required independently of the application of self-optimizing control. This is similar to the variable transformation utilizing the existing control structure as proposed by Straus and Skogestad [1]. An alternative can be seen in sampling the optimal response surface which is limited by the computational demand of this approach; it would require the solution to n p nonlinear problems whereas in the application of the proposes method, only n d +1 have to be solved in the calculation of the optimal sensitivity matrix F.
6 If it is necessary to have surrogate models for other variables, it is possible to calculate them as well, e.g. for the actual split ratios u or for additional potential measurements. However, certain limitations can be identified and need to be addressed as well: Number of Manipulated Variables The number of manipulated variables is important in considering the applicability of the approach. If n u,mv = 1, the reduction in sampling points can be negligible in cases where other sampling methods than regular grid sampling are used. Hence, there may not be a clear advantage in applying selfoptimizing control for variable reduction. However, the simpler response surface through application of self-optimizing control may still hold. If n u,mv > 1, measurement selection may result in problems associated with the non-convexity of the MIQP approach if structural zeros are defined in the selection matrix. The approach of treating each input independently can work and give good results as it is the case for the presented ammonia reactor. It is also advantageous to use a structured selection matrix H if one is using a sequentialmodular flowsheet software. Here, the manipulated variables are adjusted iteratively, and hence, it is beneficial if there is a reduced coupling between the SOC variables c. However, it cannot be generalized for all problems, if the application of a structured selection matrix improves the performance. Furthermore, it cannot be said how the solution to a nonconvex optimization problem (14) for a structured selection matrix would differ to the one calculated through treating the reactor beds independently. Application in Flowsheeting Software One may argue that the application in flowsheeting software is difficult, mostly due to the iterative procedure in which variables are adjusted in sequential-modular software. This can be especially problematic in the investigated case study, where each change in one of the manipulated variables affects the value of all controlled variables directly. As a result, the calculation time for each sampled point is higher than in the case of not applying the proposed method. This can be negated, if the measurement variables are localized close to the respective manipulated variable. As a results, a structured selection matrix H will be calculated. The problem is however less pronounced in equation-orientated flowsheeting solver. There, the application of surrogate model-based optimization results in smaller models, and hence, a simpler initialization of the models. Then, the only change is an increase in the number of nonlinear equality constraints and potentially a slightly more complex system of equations. VI. CONCLUSION Combining principles from control and surrogate modelling, a new method was derived to reduce the number of independent variables of surrogate models. In addition, the corresponding response surface may be simpler. This is caused by omitting regions in which the submodel is for example in reactor extinction. Furthermore, the method is independent of the structure of the surrogate model. Hence, it is possible to combine it with other approaches in the literature for the calculation of surrogate models. VII. ACKNOWLEDGMENT The authors gratefully acknowledge the financial support provided by Yara International ASA. REFERENCES [1] J. Straus and S. Skogestad, Minimizing the complexity of surrogate models for optimization, in 26th European Symposium on Computer Aided Process Engineering, ser. Computer Aided Chemical Engineering, Z. Kravanja and M. Bogataj, Eds. Elsevier, 2016, vol. 38, pp [2] J. Straus and S. Skogestad, Variable reduction for surrogate modelling, in Proceedings of Foundations of Computer-Aided Process Operations 2017, Tucson, AZ, USA, Jan [3] S. Skogestad, Plantwide control: the search for the selfoptimizing control structure, Journal of Process Control, vol. 10, no. 5, pp , [4] I. J. Halvorsen, S. Skogestad, J. C. Morud, and V. Alstad, Optimal selection of controlled variables, Industrial & Engineering Chemistry Research, vol. 42, no. 14, pp , [5] J. Jäschke, Y. Cao, and V. Kariwala, Self-optimizing control a survey, Annual Reviews in Control, vol. 43, no. Supplement C, pp , [6] V. Alstad and S. Skogestad, Null space method for selecting optimal measurement combinations as controlled variables, Industrial & engineering chemistry research, vol. 46, no. 3, pp , [7] V. Alstad, S. Skogestad, and E. S. Hori, Optimal measurement combinations as controlled variables, Journal of Process Control, vol. 19, no. 1, pp , [8] R. Yelchuru and S. Skogestad, Convex formulations for optimal selection of controlled variables and measurements using mixed integer quadratic programming, Journal of Process Control, vol. 22, no. 6, pp , [9] J. C. Morud and S. Skogestad, Analysis of instability in an industrial ammonia reactor, AIChE Journal, vol. 44, no. 4, pp , [10] J. Straus and S. Skogestad, Economic NMPC for heatintegrated chemical reactors, in st International Conference on Process Control (PC), June 2017, pp [11] J. Andersson, A General-Purpose Software Framework for Dynamic Optimization, PhD thesis, Arenberg Doctoral School, KU Leuven, Department of Electrical Engineering (ESAT/SCD) and Optimization in Engineering Center, Kasteelpark Arenberg 10, 3001-Heverlee, Belgium, October [12] A. Wächter and L. T. Biegler, On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming, Mathematical Programming, vol. 106, no. 1, pp , [13] B. Grimstad et al., SPLINTER: a library for multivariate function approximation with splines, , accessed: [14] J. A. Caballero and I. E. Grossmann, An algorithm for the use of surrogate models in modular flowsheet optimization, AIChE Journal, vol. 54, no. 10, pp , [15] A. Cozad, N. V. Sahinidis, and D. C. Miller, Learning surrogate models for simulation-based optimization, AIChE Journal, vol. 60, no. 6, pp , 2014.
CONTROLLED VARIABLE AND MEASUREMENT SELECTION
CONTROLLED VARIABLE AND MEASUREMENT SELECTION Sigurd Skogestad Department of Chemical Engineering Norwegian University of Science and Technology (NTNU) Trondheim, Norway Bratislava, Nov. 2010 1 NTNU, Trondheim
Incorporating Feedforward Action into Self-optimizing Control Policies
Incorporating Feedforward Action into Self-optimizing Control Policies Lia Maisarah Umar, Yi Cao and Vinay Kariwala School of Chemical and Biomedical Engineering, Nanyang Technological University, Singapore
COMBINATIONS OF MEASUREMENTS AS CONTROLLED VARIABLES: APPLICATION TO A PETLYUK DISTILLATION COLUMN.
COMBINATIONS OF MEASUREMENTS AS CONTROLLED VARIABLES: APPLICATION TO A PETLYUK DISTILLATION COLUMN. V.Alstad, S. Skogestad 1 Department of Chemical Engineering, Norwegian University of Science and Technology,
Self-Optimizing Control of a Two-Stage Refrigeration Cycle
Preprint, 11th IFAC Symposium on Dynamics and Control of Process Systems, including Biosystems June 6-8, 16. NTNU, Trondheim, Norway Self-Optimizing Control of a Two-Stage Refrigeration Cycle Adriaen Verheyleweghen,
Local Self-Optimizing Control with Average Loss Minimization
Local Self-Optimizing Control with Average Loss Minimization Vinay Kariwala, Yi Cao and S. Janardhanan Division of Chemical & Biomolecular Engineering, Nanyang Technological University, Singapore 637459
Model Predictive Control for the Self-optimized Operation in Wastewater Treatment Plants
Krist V. Gernae, Jakob K. Huusom and Rafiqul Gani (Eds.), 12th International Smposium on Process Sstems Engineering and 25th European Smposium on Computer Aided Process Engineering. 31 Ma 4 June 2015,
Optimal operation of simple refrigeration cycles Part II: Selection of controlled variables
Computers and Chemical Engineering 31 (2007) 1590 1601 Optimal operation of simple refrigeration cycles Part II: Selection of controlled variables Jørgen Bauck Jensen, Sigurd Skogestad Department of Chemical
Time scale separation and the link between open-loop and closed-loop dynamics
Time scale separation and the link between open-loop and closed-loop dynamics Antonio Araújo a Michael Baldea b Sigurd Skogestad a,1 Prodromos Daoutidis b a Department of Chemical Engineering Norwegian
Studies on Selection of Controlled Variables
Vidar Alstad Studies on Selection of Controlled Variables Doctoral thesis for the degree of doktor ingeniør Trondheim, May 2005 Norwegian University of Science and Technology Faculty of Natural Sciences
Implementation issues for real-time optimization of a crude unit heat exchanger network
1 Implementation issues for real-time optimization of a crude unit heat exchanger network Tore Lid a, Sigurd Skogestad b a Statoil Mongstad, N-5954 Mongstad b Department of Chemical Engineering, NTNU,
Lecture 1: Basics Concepts
National Technical University of Athens School of Chemical Engineering Department II: Process Analysis and Plant Design Lecture 1: Basics Concepts Instructor: Α. Kokossis Laboratory teaching staff: Α.
I would like to thank the following people for their contributions to this project:
1 Abstract This report presents optimization and control structure suggestions for a model of petroleum production wells. The model presented in the report is based the semi-realistic model of the Troll
Optimizing Control of Petlyuk Distillation: Understanding the Steady-State Behavior
Optimizing Control of Petlyuk Distillation: Understanding the Steady-State Behavior Ivar J. Halvorsen and Sigurd Skogestad Abstract. Norwegian University of Science and Technology, Department of Chemical
THE DOS AND DON TS OF DISTILLATION COLUMN CONTROL
THE DOS AND DON TS OF DISTILLATION COLUMN CONTROL Sigurd Skogestad Department of Chemical Engineering, Norwegian University of Science and Technology, N-7491 Trondheim, Norway The paper discusses distillation
ECONOMIC PLANTWIDE CONTROL: Control structure design for complete processing plants
ECONOMIC PLANTWIDE CONTROL: Control structure design for complete processing plants Sigurd Skogestad Department of Chemical Engineering Norwegian University of Science and Tecnology (NTNU) Trondheim, Norway
Feedback: Still the simplest and best solution
Feedback: Still the simplest and best solution Sigurd Skogestad Department of Chemical Engineering Norwegian Univ. of Science and Tech. (NTNU) Trondheim, Norway skoge@ntnu.no Abstract Most engineers are
Linear and nonlinear optimization with CasADi
Numerical optimization course computer exercise inear and nonlinear optimization with CasADi Prof. Dr. Moritz Diehl 1, Prof. Dr. Angelika Altmann-Dieses 2, Adrian Bürger 1,2 1 Systems Control and Optimization
The dos and don ts of distillation column control
The dos and don ts of distillation column control Sigurd Skogestad * Department of Chemical Engineering Norwegian University of Science and Technology N-7491 Trondheim, Norway Abstract The paper discusses
A Definition for Plantwide Controllability. Process Flexibility
A Definition for Plantwide Controllability Surya Kiran Chodavarapu and Alex Zheng Department of Chemical Engineering University of Massachusetts Amherst, MA 01003 Abstract Chemical process synthesis typically
Direct and Indirect Gradient Control for Static Optimisation
International Journal of Automation and Computing 1 (2005) 60-66 Direct and Indirect Gradient Control for Static Optimisation Yi Cao School of Engineering, Cranfield University, UK Abstract: Static self-optimising
Application of Decomposition Methodology to Solve Integrated Process Design and Controller Design Problems for Reactor-Separator-Recycle Systems
Proceedings of the 9th International Symposium on Dynamics and Control of Process Systems (DYCOPS 2010), Leuven, Belgium, July 5-7, 2010 Mayuresh Kothare, Moses Tade, Alain Vande Wouwer, Ilse Smets (Eds.)
This article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
Machine learning, ALAMO, and constrained regression
Machine learning, ALAMO, and constrained regression Nick Sahinidis Acknowledgments: Alison Cozad, David Miller, Zach Wilson MACHINE LEARNING PROBLEM Build a model of output variables as a function of input
CHAPTER 2: QUADRATIC PROGRAMMING
CHAPTER 2: QUADRATIC PROGRAMMING Overview Quadratic programming (QP) problems are characterized by objective functions that are quadratic in the design variables, and linear constraints. In this sense,
Self-consistent inventory control
Elvira Marie B. Aske,, and Sigurd Skogestad, Department of Chemical Engineering, Norwegian University of Science and Technology, Trondheim, Norway E-mail: skoge@ntnu.no Abstract Inventory or material balance
Real-Time Feasibility of Nonlinear Predictive Control for Semi-batch Reactors
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Real-Time Feasibility of Nonlinear Predictive Control
The dos and don ts of distillation column control
The dos and don ts of distillation column control Sigurd Skogestad * Department of Chemical Engineering Norwegian University of Science and Technology N-7491 Trondheim, Norway Abstract The paper discusses
Numerical Optimal Control Mixing Collocation with Single Shooting: A Case Study
Preprint, 11th IFAC Symposium on Dynamics and Control of Process Systems, including Biosystems Numerical Optimal Control Mixing Collocation with Single Shooting: A Case Study Anders Albert Lars Imsland
THIS paper deals with robust control in the setup associated
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 50, NO 10, OCTOBER 2005 1501 Control-Oriented Model Validation and Errors Quantification in the `1 Setup V F Sokolov Abstract A priori information required for
Consistent inventory control
Consistent inventory control Elvira Marie B. Aske,, and Sigurd Skogestad, Department of Chemical Engineering, Norwegian University of Science and Technology, Trondheim, Norway Abstract Inventory or material
Sensitivity of Optimal Operation of an Activated Sludge Process Model
UKACC International Conference on Control 2012 Cardiff, UK, 3-5 September 2012 Sensitivity of Optimal Operation of an Activated Sludge Process Model Antonio Araujo, Simone Gallani, Michela Mulas and Sigurd
15-780: LinearProgramming
15-780: LinearProgramming J. Zico Kolter February 1-3, 2016 1 Outline Introduction Some linear algebra review Linear programming Simplex algorithm Duality and dual simplex 2 Outline Introduction Some linear
PRACTICAL CONTROL OF DIVIDING-WALL COLUMNS
Distillation Absorption 2010. Strandberg, S. Skogestad and I.. Halvorsen All rights reserved by authors as per DA2010 copyright notice PRACTICAL CONTROL OF DIVIDING-WALL COLUMNS ens Strandberg 1, Sigurd
Academic Editor: Dominique Bonvin Received: 26 November 2016; Accepted: 13 February 2017; Published: 27 February 2017
Article Sensitivity-Based Economic NMPC with a Path-Following Approach Eka Suwartadi 1, Vyacheslav Kungurtsev 2 and Johannes Jäschke 1, * 1 Department of Chemical Engineering, Norwegian University of Science
Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati
Process Design Decisions and Project Economics Prof. Dr. V. S. Moholkar Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 2 Flowsheet Synthesis (Conceptual Design of
Tightening a Discrete Formulation of the Quadratic Assignment Problem
A publication of 1309 CHEMICAL ENGINEERING TRANSACTIONS VOL. 32, 2013 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 2013, AIDIC Servizi S.r.l., ISBN 978-88-95608-23-5; ISSN 1974-9791 The Italian
Robust control for a multi-stage evaporation plant in the presence of uncertainties
Preprint 11th IFAC Symposium on Dynamics and Control of Process Systems including Biosystems June 6-8 16. NTNU Trondheim Norway Robust control for a multi-stage evaporation plant in the presence of uncertainties
Optimal Controlled Variables for Parallel Process Units
Parallel Process Johannes Jäschke, Sigurd Skogestad Department of Chemical Engineering NTNU Trondheim, Norway NPCW 2012 NPCW 2012 1 / 20 Main Idea & Outline of this talk Main idea Avoid online re-optimization
Improved PI control for a surge tank satisfying level constraints
Integral-Derivative Control, Ghent, Belgium, May 9-, 28 FrBT4.6 Improved PI control for a surge tank satisfying level constraints Adriana Reyes-Lúa Christoph Josef Backi Sigurd Skogestad, Department of
Economic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples
1 Economic Model Predictive Control Historical Perspective and Recent Developments and Industrial Examples Public Trial Lecture Candidate: Vinicius de Oliveira Department of Chemical Engineering, Faculty
Facing New Products Demand through Simultaneous Structural and Operational Decisions in the Design of the Control Recipe
Facing New Products Demand through Simultaneous Structural and Operational Decisions in the Design of the Control Recipe Marta Moreno-Benito, Antonio Espuña* Chemical Engineering Department, Univesitat
Real-Time Implementation of Nonlinear Predictive Control
Real-Time Implementation of Nonlinear Predictive Control Michael A. Henson Department of Chemical Engineering University of Massachusetts Amherst, MA WebCAST October 2, 2008 1 Outline Limitations of linear
TKP Kjemisk prosessteknologi, fordypningsprosjekt
TKP4580 1 Kjemisk prosessteknologi, fordypningsprosjekt Kandidat 10023 Oppgaver Oppgavetype Vurdering Status 1 Handing in the specialization report. Filopplasting Manuell poengsum Levert TKP4580 1 Kjemisk
Dynamic Real-Time Optimization: Linking Off-line Planning with On-line Optimization
Dynamic Real-Time Optimization: Linking Off-line Planning with On-line Optimization L. T. Biegler and V. Zavala Chemical Engineering Department Carnegie Mellon University Pittsburgh, PA 15213 April 12,
Synergy between Data Reconciliation and Principal Component Analysis.
Plant Monitoring and Fault Detection Synergy between Data Reconciliation and Principal Component Analysis. Th. Amand a, G. Heyen a, B. Kalitventzeff b Thierry.Amand@ulg.ac.be, G.Heyen@ulg.ac.be, B.Kalitventzeff@ulg.ac.be
Mathematical Modeling of Chemical Processes. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University
Mathematical Modeling of Chemical Processes Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University Chapter Objectives End of this chapter, you should be
Invariants for Optimal Operation of Process Systems
Johannes Jäschke Invariants for Optimal Operation of Process Systems Doctoral thesis for the degree of philosophiae doctor Trondheim, June 2011 Norwegian University of Science and Technology The Faculty
Model-based Conceptual Design and Tool Support for the Development of Continuous Chemical Processes
Ian David Lockhart Bogle and Michael Fairweather (Editors), Proceedings of the 22nd European Symposium on Computer Aided Process Engineering, 17-20 June 2012, London. 2012 Elsevier B.V. All rights reserved
Static Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
Journal of Process Control
Journal of Process Control 3 (03) 404 44 Contents lists available at SciVerse ScienceDirect Journal of Process Control j ourna l ho me pag e: www.elsevier.com/locate/jprocont Algorithms for improved fixed-time
A Generalization of Principal Component Analysis to the Exponential Family
A Generalization of Principal Component Analysis to the Exponential Family Michael Collins Sanjoy Dasgupta Robert E. Schapire AT&T Labs Research 8 Park Avenue, Florham Park, NJ 7932 mcollins, dasgupta,
Efficient Handling of the Implicit Constraints Problem for the ASPEN MINLP Synthesizer
2006 Ind. Eng. Chem. Res. 1993,32, 2006-2011 Efficient Handling of the Implicit Constraints Problem for the ASPEN MINLP Synthesizer U. M. Diwekar' and E. S. Rubin Environmental Institute h Department of
Paris'09 ECCI Eduardo F. Camacho MPC Constraints 2. Paris'09 ECCI Eduardo F. Camacho MPC Constraints 4
Outline Constrained MPC Eduardo F. Camacho Univ. of Seville. Constraints in Process Control. Constraints and MPC 3. Formulation of Constrained MPC 4. Illustrative Examples 5. Feasibility. Constraint Management
Economic Plant-wide Control. Sigurd Skogestad
JWST-c JWST-Rangaiah November, 0 : Printer: Yet to come 0 0 0 Economic Plant-wide Control. Introduction Sigurd Skogestad Department of Chemical Engineering, Norwegian University of Science and Technology
ROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE. Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt
ROBUST STABLE NONLINEAR CONTROL AND DESIGN OF A CSTR IN A LARGE OPERATING RANGE Johannes Gerhard, Martin Mönnigmann, Wolfgang Marquardt Lehrstuhl für Prozesstechnik, RWTH Aachen Turmstr. 46, D-5264 Aachen,
SUPPORT VECTOR MACHINE FOR THE SIMULTANEOUS APPROXIMATION OF A FUNCTION AND ITS DERIVATIVE
SUPPORT VECTOR MACHINE FOR THE SIMULTANEOUS APPROXIMATION OF A FUNCTION AND ITS DERIVATIVE M. Lázaro 1, I. Santamaría 2, F. Pérez-Cruz 1, A. Artés-Rodríguez 1 1 Departamento de Teoría de la Señal y Comunicaciones
MPC Infeasibility Handling
MPC Handling Thomas Wiese, TU Munich, KU Leuven supervised by H.J. Ferreau, Prof. M. Diehl (both KUL) and Dr. H. Gräb (TUM) October 9, 2008 1 / 42 MPC General MPC Strategies 2 / 42 Linear Discrete-Time
Semi-empirical Process Modelling with Integration of Commercial Modelling Tools
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Semi-empirical Process Modelling with Integration
Convexification of Mixed-Integer Quadratically Constrained Quadratic Programs
Convexification of Mixed-Integer Quadratically Constrained Quadratic Programs Laura Galli 1 Adam N. Letchford 2 Lancaster, April 2011 1 DEIS, University of Bologna, Italy 2 Department of Management Science,
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title Bidirectional branch and bound for controlled variable selection Part II : exact local method for self-optimizing
SPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
Theoretical Models of Chemical Processes
Theoretical Models of Chemical Processes Dr. M. A. A. Shoukat Choudhury 1 Rationale for Dynamic Models 1. Improve understanding of the process 2. Train Plant operating personnel 3. Develop control strategy
Active Fault Diagnosis for Uncertain Systems
Active Fault Diagnosis for Uncertain Systems Davide M. Raimondo 1 Joseph K. Scott 2, Richard D. Braatz 2, Roberto Marseglia 1, Lalo Magni 1, Rolf Findeisen 3 1 Identification and Control of Dynamic Systems
Example 8: CSTR with Multiple Solutions
Example 8: CSTR with Multiple Solutions This model studies the multiple steady-states of exothermic reactions. The example is from Parulekar (27) which in turn was modified from one by Fogler (1999). The
Notes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
Written Examination
Division of Scientific Computing Department of Information Technology Uppsala University Optimization Written Examination 202-2-20 Time: 4:00-9:00 Allowed Tools: Pocket Calculator, one A4 paper with notes
Modeling and Optimization of Multistream Heat Exchangers with Phase Transition
Department of Chemical Engineering Specialisation Project Modeling and Optimization of Multistream Heat Exchangers with Phase Transition Written by: Axel Lødemel Holene Supervisors: Sigurd Skogestad Johannes
Course on Model Predictive Control Part II Linear MPC design
Course on Model Predictive Control Part II Linear MPC design Gabriele Pannocchia Department of Chemical Engineering, University of Pisa, Italy Email: g.pannocchia@diccism.unipi.it Facoltà di Ingegneria,
MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands
MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands Axel Nyberg Åbo Aademi University Ignacio E. Grossmann Dept. of Chemical Engineering, Carnegie Mellon University,
6. Iterative Methods for Linear Systems. The stepwise approach to the solution...
6 Iterative Methods for Linear Systems The stepwise approach to the solution Miriam Mehl: 6 Iterative Methods for Linear Systems The stepwise approach to the solution, January 18, 2013 1 61 Large Sparse
Elvira Marie B. Aske
Elvira Marie B. Aske Design of plantwide control systems with focus on maximizing throughput Thesis for the degree of philosophiae doctor Trondheim, March 2009 Norwegian University of Science and Technology
Nonlinear Behaviour of a Low-Density Polyethylene Tubular Reactor-Separator-Recycle System
European Symposium on Computer Arded Aided Process Engineering 15 L. Puigjaner and A. Espuña (Editors) 2005 Elsevier Science B.V. All rights reserved. Nonlinear Behaviour of a Low-Density Polyethylene
Chapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
An Optimization-based Approach to Decentralized Assignability
2016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 2016 Boston, MA, USA An Optimization-based Approach to Decentralized Assignability Alborz Alavian and Michael Rotkowitz Abstract
The Growth of Functions. A Practical Introduction with as Little Theory as possible
The Growth of Functions A Practical Introduction with as Little Theory as possible Complexity of Algorithms (1) Before we talk about the growth of functions and the concept of order, let s discuss why
Simultaneous Synthesis of Heat Exchanger Networks with Operability Considerations: Flexibility and Controllability
Simultaneous Synthesis of Heat Exchanger Networks with Operability Considerations: Flexibility and Controllability Marcelo Escobar *, Jorge O. Trierweiler Chemical Engineering Department, Federal University
Adaptive Dual Control
Adaptive Dual Control Björn Wittenmark Department of Automatic Control, Lund Institute of Technology Box 118, S-221 00 Lund, Sweden email: bjorn@control.lth.se Keywords: Dual control, stochastic control,
A Survey for the Selection of Control Structure for Distillation Columns Based on Steady State Controllability Indexes
Iranian Journal of Chemical Engineering Vol. 6, No. 2 (Spring), 2009, IAChE A Survey for the Selection of Control Structure for Distillation Columns Based on Steady State Controllability Indexes K. Razzaghi,
Software for Integer and Nonlinear Optimization
Software for Integer and Nonlinear Optimization Sven Leyffer, leyffer@mcs.anl.gov Mathematics & Computer Science Division Argonne National Laboratory Roger Fletcher & Jeff Linderoth Advanced Methods and
Dividing wall columns for heterogeneous azeotropic distillation
Dividing wall columns for heterogeneous azeotropic distillation Quang-Khoa Le 1, Ivar J. Halvorsen 2, Oleg Pajalic 3, Sigurd Skogestad 1* 1 Norwegian University of Science and Technology (NTNU), Trondheim,
Majorization Minimization - the Technique of Surrogate
Majorization Minimization - the Technique of Surrogate Andersen Ang Mathématique et de Recherche opérationnelle Faculté polytechnique de Mons UMONS Mons, Belgium email: manshun.ang@umons.ac.be homepage:
NONLINEAR PREDICTIVE CONTROLLER: INFLUENCE OF ITERATIONS NUMBER ON REGULATION QUALITY
NONLINEAR PREDICIVE CONROLLER: INFLUENCE OF IERAIONS NUMBER ON REGULAION QUALIY P. Škrabánek Faculty of Electrical engineering and Informatics Department of process control University of Pardubice Abstract
Model Predictive Control
Model Predictive Control Davide Manca Lecture 6 of Dynamics and Control of Chemical Processes Master Degree in Chemical Engineering Davide Manca Dynamics and Control of Chemical Processes Master Degree
Interior-Point Methods for Linear Optimization
Interior-Point Methods for Linear Optimization Robert M. Freund and Jorge Vera March, 204 c 204 Robert M. Freund and Jorge Vera. All rights reserved. Linear Optimization with a Logarithmic Barrier Function
Elvira Marie B. Aske
Elvira Marie B. Aske Design of plantwide control systems with focus on maximizing throughput Doctoral thesis for the degree of philosophiae doctor Trondheim, September 2008 Norwegian University of Science
inputs. The velocity form is used in the digital implementation to avoid wind-up [7]. The unified LQR scheme has been developed due to several reasons
A LQR Scheme for SCR Process in Combined-Cycle Thermal Power Plants Santo Wijaya 1 Keiko Shimizu 1 and Masashi Nakamoto 2 Abstract The paper presents a feedback control of Linear Quadratic Regulator (LQR)
An Investigation into the Nonlinearity of Polyethylene Reactor Operation. Department of Chemical Engineering Queen s University
An Investigation into the Nonlinearity of Polyethylene Reactor Operation M.-A. Benda,, K. McAuley and P.J. McLellan Department of Chemical Engineering Queen s University Outline Motivation Background:
Data reconciliation and optimal operation of a catalytic naphtha reformer
Data reconciliation and optimal operation of a catalytic naphtha reformer Tore Lid Statoil Mongstad 5954 Mongstad Sigurd Skogestad Department of Chemical Engineering Norwegian Univ. of Science and Technology
3.1 Overview 3.2 Process and control-loop interactions
3. Multivariable 3.1 Overview 3.2 and control-loop interactions 3.2.1 Interaction analysis 3.2.2 Closed-loop stability 3.3 Decoupling control 3.3.1 Basic design principle 3.3.2 Complete decoupling 3.3.3
Constrained Nonlinear Optimization Algorithms
Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu Institute for Mathematics and its Applications University of Minnesota August 4, 2016
Nonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control
Nonlinear Programming (Hillier, Lieberman Chapter 13) CHEM-E7155 Production Planning and Control 19/4/2012 Lecture content Problem formulation and sample examples (ch 13.1) Theoretical background Graphical
Matrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
2. Review on Material Balances
2. Review on Material Balances Objectives After completing this chapter, students should be able to recall the principle of the conservation of mass recall the principle of the stoichiometry of chemical
Finding normalized and modularity cuts by spectral clustering. Ljubjana 2010, October
Finding normalized and modularity cuts by spectral clustering Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu Ljubjana 2010, October Outline Find
Support Vector Regression Machines
Support Vector Regression Machines Harris Drucker* Chris J.C. Burges** Linda Kaufman** Alex Smola** Vladimir Vapnik + *Bell Labs and Monmouth University Department of Electronic Engineering West Long Branch,
Lecture Notes: Geometric Considerations in Unconstrained Optimization
Lecture Notes: Geometric Considerations in Unconstrained Optimization James T. Allison February 15, 2006 The primary objectives of this lecture on unconstrained optimization are to: Establish connections
Indicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
Integrated Knowledge Based System for Process Synthesis
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Integrated Knowledge Based System for Process Synthesis
Global Optimization of Ordinary Differential Equations Models
Global Optimization of Ordinary Differential Equations Models Angelo Lucia*, Meghan L. Bellows and Leah M. Octavio Department of Chemical Engineering, University of Rhode Island, Kingston, RI 02881 Abstract
A three-level MILP model for generation and transmission expansion planning
A three-level MILP model for generation and transmission expansion planning David Pozo Cámara (UCLM) Enzo E. Sauma Santís (PUC) Javier Contreras Sanz (UCLM) Contents 1. Introduction 2. Aims and contributions