<p> Control structure design is an essential step in control system synthesis and has big impact on achievable closed-loop performance. This thesis develops a systematic approach of selecting optimal control structures based on closed-loop dynamic performance and other criteria, such as integrity.</p> <p> The main contribution of this thesis is a rigorous mathematical formulation for control structure design problem that includes full closed-loop transient analysis with additional integrity requirement. The multi-objective framework is extendable so that different control performance objectives can be easily added. Unique process requirements and engineer inputs can be taken into account as additional constraints. The proposed formulation is a Mixed Integer Nonlinear Programming (MINLP) with complementarity constraints. The research scope is limited to linear process models and linear controller algorithms.</p> <p> The tailored solving strategy that makes this challenging problem computationally tractable is introduced in this thesis. The modified Branch and Bound algorithm takes advantage of the special problem structure by using control knowledge to generate valid lower bound efficiently. Prior knowledge can be cooperated as heuristic tuning parameters to guide the solving process so that a reasonably good solution can be found early in the solving process. The complexity study shows the solving strategy can attack design problem size up to 8x8. Considering the percentage of good structures needing evaluation will decrease with problem size even larger problems will be tractable.</p> <p> The common control structures in process industries, such as square and nonsquare Single-Input-Single-Output (SISO) loop pairing using PID controller and block-centralized structure using Model Predictive Controller (MPC), are addressed in this thesis. The usefulness of this research has been demonstrated by several case studies, include Tennessee Eastman problem. The proposed methodology finds a physically sound pairing with good performance for Tennessee Eastman problem in less than one hour, while several off-the-shelf NLP, MINLP and global solvers cannot find a solution in five days.</p> / Thesis / Doctor of Philosophy (PhD)

L’objectif de ces travaux est de développer une méthode d’aide à la conception des réseaux de chaleur urbains (RCU). Cette méthode utilise un modèle de type MINLP (Mixed Integer Non Linear Programming) pour l’optimisation simultanée de la configuration et du dimensionnement d’un RCU. Aux variables continues pour l’aide au dimensionnement (température, vitesse, diamètre, aire des échangeurs), s’ajoutent des variables binaires aidant à définir la configuration du réseau (maillage et choix des technologies). La fonction objectif à minimiser est le coût total (capex et opex), qui est soumise à un ensemble de contraintes non linéaires (p. ex. pertes thermiques et de charge, bilans). La méthode développée dans ce manuscrit offre la possibilité de connecter en cascade des consommateurs n’ayant pas les mêmes besoins en température, et de réaliser des réseaux bouclés (une canalisation par tranchée). Elle permet aussi de choisir : les consommateurs à connecter au RCU, le ou les sites de production ainsi que le type de technologie utilisée. Enfin la bonne prise en compte de la physique permet de choisir le meilleur compromis entre pertes thermiques et pertes de charge, sur une large gamme de température. Cette formulation permet donc d’optimiser des réseaux de 4éme génération et de démontrer la rentabilité de l’intégration d’EnR&R sur le long terme (30 ans). Un premier travail est réalisé afin de proposer une méthodologie de résolution en plusieurs étapes permettant l’obtention de l’optimum global. Différents cas d’études académiques sont utilisés pour présenter les intérêts multiples de cette formulation. Enfin la comparaison avec un réseau existant a permis de démontrer la cohérence des résultats du modèle et a servi de base pour l’optimisation d’un cas d’étude de grande dimension. Plusieurs études de sensibilité post-optimale sont réalisées afin de démontrer l’intérêt de cet outil pour l’aide à la conception initiale ou l’extension de RCU existants. / The aim of this thesis is to develop a method that provides design assistance for District Heating Network (DHN). This tool allows simultaneously the optimization of the configuration and its sizing, thanks to an MINLP formulation (Mixed Integer Non-Linear Programming). Binary variables help to choose the optimal configuration (network layout and technologies of production), whereas continuous variables help DHN sizing (temperature, diameter, velocity, heat exchanger area, thermal generating capacity …). The objective function to minimize is the total cost (capex and opex), subjected to numerous nonlinear constraints (e.g. thermal losses, pressure drop, energy balance).This method enables to design temperature cascade between consumers, when consumer temperature requirements are different, and also looped network (only one pipe in one trench). It helps also the decision to connect (or not) consumers to the main network and also the location(s) and type(s) of the heating plant. Moreover, the arbitrage between heat losses and pressure drops is taken into account thanks to physical considerations (non-linear equations). Eventually, it is possible to design 4th generation DHN and prove their financial profitability over the long terms (30 years). First a multi-step resolution strategy is proposed to ensure finding global optimum of the complex MINLP problem. Then academic study cases are analyzed to underline the numerous assets of the formulation. Finally, the optimal design compared to an existing DHN ensures the consistency of the method and allows to build a study case at a wider scale, which can be solved thanks to the comprehensive strategy developed. The design assistance method is available for initial design as well as for extension of existing DHN.

Dirigé par la société Fonroche Géothermie, un consortium de dix partenaires participe au projet FONGEOSEC qui s’inscrit dans le cadre des Investissements d’Avenir de l’ADEME. Ce projet a pour but de concevoir et de réaliser un démonstrateur innovant de centrale géothermique haute enthalpie. L’énergie, ainsi récupérée en profondeur, servira à la cogénération d’électricité et de chaleur. L’une des étapes du projet correspond à l’objectif de cette thèse : développer une méthodologie pour la conception optimale des activités de surface de la centrale géothermique. Il s’agit donc de formuler le problème d’optimisation, de proposer une stratégie de résolution robuste et enfin, de mettre en oeuvre cette stratégie grâce à un outil logiciel.Dans l’outil ainsi développé, la répartition entre la production d’électricité et de chaleur s’effectue en parallèle. Le fluide géothermal est séparé en deux courants, l’un alimentant un Cycle Organique de Rankine (ORC : Organic Rankine Cycle) pour la production d’électricité, et l’autre étant relié à un Réseau de Chaleur Urbain (RCU) pour la distribution de la chaleur. Chaque constituant de l’ORC est dimensionné et la topologie du RCU est déterminée. Cet outil permet alors de déterminer simultanément :quelle est la meilleure répartition entre production d’électricité et de chaleur,quelles sont les meilleures dimensions pour les composants de l’ORC,et quelle est la meilleure topologie du RCU.Concernant l’ORC, l’outil permettra de savoir si l’utilisation d’un éventuel récupérateur de chaleur interne (IHE : Internal Heat Exchanger) est avantageuse ou non. Du point de vue du RCU, tous les consommateurs (sous-stations) envisagés ne sont pas obligatoires. L’outil permettra de choisir quels consommateurs relier au réseau et dans quelle disposition. L’utilisation de variables discrètes est alors nécessaire et le problème d’optimisation ainsi résolu est un problème de type MINLP (Mixed Integer Non Linear Programming).Une méthodologie de résolution permettant l’obtention d’une solution de « confiance » (probablement, mais non certainement, l’optimum global) est proposée. Cette stratégie de résolution est testée pour différents cas d’étude proches des conditions du projet FONGEOSEC. La stabilité et la robustesse de cette stratégie sont alors mises en avant. Une analyse économique et une analyse énergétique sont réalisées. La résolution multi-objectif est alors effectuée dans le but de fournir le meilleur compromis entre bénéfices annuels nets et destruction d’exergie. Pour finir, la diversité des résultats montre qu’il n’est pas satisfaisant de dissocier les études des deux systèmes (ORC et RCU) et démontre l’intérêt de l’outil développé. / A consortium of ten partners, led by “FONROCHE Géothermie”, works on the FONGEOSEC project, an “Investissement d’Avenir” organized by the French Agency for Environment and Energy (ADEME). The aim of this project is to design and create an innovative demonstrator of a high-energy geothermal power plant. The geothermal energy will be used to produce electricity and heat. Among other tasks, this project aims to develop a support tool for the optimal design of the structure for the surface activities in the geothermal plant.Within the developed tool, the repartition between electricity and heat production is in parallel. The geothermal fluid is split in two streams, one is used for an Organic Rankine Cycle (ORC) for electricity production, and the other is connected to a District Heating Network (DHN) for the heat distribution. This tool enables to determine simultaneously:which is the best repartition between electricity and heat,which is the best sizing for ORC components,which is the best configuration for the DHN.About the ORC, the tool will enable to decide if the use of an Internal Heat Exchanger (IHE) is interesting or not. For the DHN point of view, all the consumers envisaged are not mandatory. The tool will enable to choose which consumers it is better to connect to the network and in which disposition. The use of discrete variables is necessary and the optimization problem to be solved is a MINLP (Mixed Integer Non Linear Programming) problem.A solution strategy is implemented in order to obtain a confident solution with a determinist algorithm. This strategy is tested for different study cases close to FONGEOSEC conditions. Stability and Robustness of this strategy are then highlighted. An economic and an exergetic analysis are carried out. In order to find a good compromise between the two objectives, a multi-objective solution is performed. Finally, the diversity of results obtained shows it is not suitable to dissociate ORC and DHN studies and shows the interest of the developed tool.

La thèse présente un modèle d’optimisation mathématique ainsi qu’une méthodologie d’étude pour la conception optimale de réseaux de chaleur et de froid flexibles et intégrant des technologies innovantes (les smart réseaux ou réseaux 4e génération). Dans ce modèle, différentes contraintes énergétiques, économiques et environnementales sont alors considérées. Le modèle considère des réseaux de chaud (ou de froid) à différents niveaux de température, ainsi que différents sites de production et demande d’énergie. Le tracé du réseau ainsi que la configuration de l’ensemble des systèmes de production sont obtenus en minimisant soit l'exergie totale consommée soit le coût total pour l’investissement et l’opération des systèmes. Le modèle mathématique développé est formulé en un problème MINLP multi-période. Les contraintes du modèle sont présentées dans plusieurs modèles imbriqués. Le premier modèle M1 est multi-période et inclut les contraintes d’échanges avec le réseau et un moyen de stockage thermique. Le deuxième modèle, appelé M2, contient l’ensemble du modèle M1 ainsi que des nouvelles contraintes permettant de dimensionner des pompes à chaleur à COP variable. Le troisième modèle M3 inclut le modèle M2 ainsi que de nouvelles contraintes pour le tracé géographique des réseaux et la prise en compte des flux d’électricité. Une méthodologie d’étude est aussi présentée, permettant, à l’aide du modèle M3 développé, de traiter des cas d’étude complexes et réalistes. La méthodologie propose une décomposition du problème selon trois étapes consécutives. Cette méthodologie est alors utilisée et illustrée pour un cas d’étude complexe : l’optimisation d’un nouveau quartier, considérant des PAC, un stockage de chaleur saisonnier et de la production PV. / This thesis presents a mathematical optimization model and a methodology for the optimal design of district heating and cooling networks. The various constraints in energy balances, costs and environmental criteria can be considered in order to propose structures of flexible and innovative networks.The model allows the consideration of heating (or cooling) networks at different temperature levels, as well as different production sites and heat demands. The network’s layout as well as the optimal architecture of the heat production systems are achieved by minimizing either the total exergy consumed or the total cost for the investment and operation of the systems.The developed mathematical model is formulated into a multi-period MINLP problem. The constraints of the model are presented in nested models. The first model M1 considers the multi-period aspect and includes the constraints related to the heat exchanges between production/demand streams with the network and a thermal storage. The second model, called M2, contains the entire model M1 as well as new constraints for sizing heat pumps with variable COP. The model M3 includes the model M2 as well as new constraints for the geographical layout of the networks and the consideration of electricity balance in the problem.A methodology is also presented making it possible, using the model M3, to deal with complex and realistic case studies. The methodology proposes a decomposition of the problem following three consecutive steps. This methodology is then used for the optimization of a new district, considering heat pumps, seasonal heat storage and PV production.

An accurate model for RO process has significant importance in the simulation and optimization proposes. A steady state model of RO process is developed based on solution diffusion theory to describe the permeation through membrane and thin film approach is used to describe the concentration polarization. The model is validated against the operation data reported in the literature. For the sake of clear understanding of the interaction of feed temperature and salinity on the design and operation of RO based desalination systems, simultaneous optimization of design and operation of RO network is investigated based on two-stage RO superstructure via MINLP approach. Different cases with several feed concentrations and seasonal variation of seawater temperature are presented. Also, the possibility of flexible scheduling in terms of the number of membrane modules required in operation in high and low temperature seasons is investigated A simultaneous modelling and optimization method for RO system including boron removal is then presented. A superstructure of the RO network is developed based on double pass RO network (two-stage seawater pass and one-stage brackish water pass). The MINLP problem based on the superstructure is used to find out an optimal RO network which will minimize the total annualized cost while fulfilling a given boron content limit. The effect of pH on boron rejection is investigated at deferent seawater temperatures. The optimal operation policy of RO system is then studied in this work considering variations in freshwater demand and with changing seawater temperature throughout the day. A storage tank is added to the RO layout to provide additional operational flexibility and to ensure the availability of freshwater at all times. Two optimization problems are solved incorporating two seawater temperature profiles, representing summer and winter seasons. The possibility of flexible scheduling of cleaning and maintenance of membrane modules is investigated. Then, the optimal design and operation of RO process is studied in the presence of membrane fouling and including several operational variations such as variable seawater temperature. The cleaning schedule of single stage RO process is formulated as MINLP problem using spiral wound modules. NNs based correlation has been developed based on the actual fouling data which can be used for estimating the permeability decline factors. The correlation based on actual data to predict the annual seawater temperature profile is also incorporated in the model. The proposed optimization procedure identified simultaneously the optimal maintenance schedule of RO network including its design parameters and operating policy. The steady state model of RO process is used to study the sensitivity of different operating and design parameters on the plant performance. A non-linear optimization problem is formulated to minimize specific energy consumption at fixed product flow rate and quality while optimizing the design and operating parameters. Then the MINLP formulation is used to find the optimal designs of RO layout for brackish water desalination. A variable fouling profile along the membrane stages is introduced to see how the network design and operation of the RO system are to be adjusted Finally, a preliminary control strategy for RO process is developed based on PID control algorithm and a first order transfer function (presented in the Appendix).

La superposition de topologies virtuelles à la topologie physique d'un réseau est un des principaux mécanismes de l'ingénierie de trafic. Soit un réseau physique d'une certaine topologie et capacité fixées et une matrice de trafic à véhiculer, il s'agit trouver une topologie logique permettant de mapper de manière optimale la matrice de trafic sur le réseau physique. Lors de l'évolution de la matrice de trafic sur des échelles de temps longues, il faudra agir sur le layout. La première contribution concerne la définition de fonctions de coût mieux adaptées à la réalité d'un opérateur, la deuxième contribution concerne la prise en compte du coût de changement du layout. Il s'avère intéressant d'un point de vue opérateur de réduire la complexité du layout, mesurée comme une fonction du nombre de chemins virtuels. Nous avons donc formulé divers problèmes de minimisation de la complexité du layout sous des contraintes de QoS. Il s'agit d'une modélisation réaliste mais qui engendre des modèles difficiles à résoudre. Nous avons développés des heuristiques qui permet de trouver des solutions approchées pour des réseaux de grande taille. Nous avons montré que la complexité des layouts peut être significativement réduite en comparaison avec celle obtenue suite à l'optimisation des fonctions de coût classiques. Le changement du layout implique d'une part un coût d'opération et d'autre part peut engendrer des coupures de service qui affecteront directement le coût d'opération. Nous avons formulé une famille de problèmes prenant en compte le coût de reconfiguration du layout. L'une des heuristiques citées a été adaptée pour analyser ces nouveaux problèmes.

Submitted by Daniella Santos (daniella.santos@fgv.br) on 2010-03-23T12:01:36Z No. of bitstreams: 1 Dissertacao_Jefferson_Donizeti_Bertolai.pdf: 274926 bytes, checksum: a60c8343a27883dd5e7f529f517bc2e8 (MD5) / Approved for entry into archive by Andrea Virginio Machado(andrea.machado@fgv.br) on 2010-03-23T12:31:24Z (GMT) No. of bitstreams: 1 Dissertacao_Jefferson_Donizeti_Bertolai.pdf: 274926 bytes, checksum: a60c8343a27883dd5e7f529f517bc2e8 (MD5) / Made available in DSpace on 2010-03-24T12:48:06Z (GMT). No. of bitstreams: 1 Dissertacao_Jefferson_Donizeti_Bertolai.pdf: 274926 bytes, checksum: a60c8343a27883dd5e7f529f517bc2e8 (MD5) Previous issue date: 2009-12-08 / The dificulty in characterizing non-stationary allocations or equilibria is one of the main explanations for the use of concepts and assumptions that trivialize the dynamics of the economy. This di¢ culty is especially critical in Monetary Theory, in which the dimensionality of the problem is high even for very simple models. In this context, this paper reports the computational strategy for implementing the recursive method proposed by Monteiro and Cavalcanti (2006), which allows you to calculate the optimal sequence (possibly non-stationary) of distributions of money in an extension of the model proposed by Kiyotaki and Wright (1989). Three aspects of this calculation are emphasized: (i) the computational implementation of the plannerís problem involves the choice of continuous and discrete variables that maximize a nonlinear function and satisÖes nonlinear constraints; (ii) the objective function of this problem is not concave and constraints are not convex, and (iii) the set of admissible choices is not known a priori. The goal is to document the di¢ culties involved, the proposed solutions and available methods and resources to implement the numerical characterization of e¢ cient monetary dynamics under the assumption of random matching. / A dificuldade em se caracterizar alocações ou equilíbrios não estacionários é uma das principais explicações para a utilização de conceitos e hipóteses que trivializam a dinâmica da economia. Tal dificuldade é especialmente crítica em Teoria Monetária, em que a dimensionalidade do problema é alta mesmo para modelos muito simples. Neste contexto, o presente trabalho relata a estratégia computacional de implementação do método recursivo proposto por Monteiro e Cavalcanti (2006), o qual permite calcular a sequência ótima (possivelmente não estacionária) de distribuições de moeda em uma extensão do modelo proposto por Kiyotaki e Wright (1989). Três aspectos deste cálculo são enfatizados: (i) a implementação computacional do problema do planejador envolve a escolha de variáveis contínuas e discretas que maximizem uma função não linear e satisfaçam restrições não lineares; (ii) a função objetivo deste problema não é côncava e as restrições não são convexas; e (iii) o conjunto de escolhas admissíveis não é conhecido a priori. O objetivo é documentar as dificuldades envolvidas, as soluções propostas e os métodos e recursos disponíveis para a implementação numérica da caracterização da dinâmica monetária eficiente sob a hipótese de encontros aleatórios.

no / In this work, a model based MINLP (mixed integer nonlinear programming) optimisation framework is developed for evaluating boron rejection in a reverse osmosis (RO) desalination process. A mathematical model (for the RU process) based on solution diffusion model and thin film theory is incorporated in the optimisation framework. A superstructure of the RU network is developed which includes two passes: (a) seawater pass containing normal two-stage RU system housing seawater membrane modules and (b) the brackish water pass (BW) accommodating brackish water membrane modules. For fixed freshwater demand, the objective of this work is to demonstrate the effectiveness of the MINLP approach for analyzing and optimizing the design and operation of RU network while attaining desired limit on boron concentration in the freshwater produced. The effect of seasonal variation in seawater temperature and pH on boron removal efficiency is also discussed.

An accurate model for RO process has significant importance in the simulation and optimization proposes. A steady state model of RO process is developed based on solution diffusion theory to describe the permeation through membrane and thin film approach is used to describe the concentration polarization. The model is validated against the operation data reported in the literature. For the sake of clear understanding of the interaction of feed temperature and salinity on the design and operation of RO based desalination systems, simultaneous optimization of design and operation of RO network is investigated based on two-stage RO superstructure via MINLP approach. Different cases with several feed concentrations and seasonal variation of seawater temperature are presented. Also, the possibility of flexible scheduling in terms of the number of membrane modules required in operation in high and low temperature seasons is investigated A simultaneous modelling and optimization method for RO system including boron removal is then presented. A superstructure of the RO network is developed based on double pass RO network (two-stage seawater pass and one-stage brackish water pass). The MINLP problem based on the superstructure is used to find out an optimal RO network which will minimize the total annualized cost while fulfilling a given boron content limit. The effect of pH on boron rejection is investigated at deferent seawater temperatures. The optimal operation policy of RO system is then studied in this work considering variations in freshwater demand and with changing seawater temperature throughout the day. A storage tank is added to the RO layout to provide additional operational flexibility and to ensure the availability of freshwater at all times. Two optimization problems are solved incorporating two seawater temperature profiles, representing summer and winter seasons. The possibility of flexible scheduling of cleaning and maintenance of membrane modules is investigated. Then, the optimal design and operation of RO process is studied in the presence of membrane fouling and including several operational variations such as variable seawater temperature. The cleaning schedule of single stage RO process is formulated as MINLP problem using spiral wound modules. NNs based correlation has been developed based on the actual fouling data which can be used for estimating the permeability decline factors. The correlation based on actual data to predict the annual seawater temperature profile is also incorporated in the model. The proposed optimization procedure identified simultaneously the optimal maintenance schedule of RO network including its design parameters and operating policy. The steady state model of RO process is used to study the sensitivity of different operating and design parameters on the plant performance. A non-linear optimization problem is formulated to minimize specific energy consumption at fixed product flow rate and quality while optimizing the design and operating parameters. Then the MINLP formulation is used to find the optimal designs of RO layout for brackish water desalination. A variable fouling profile along the membrane stages is introduced to see how the network design and operation of the RO system are to be adjusted Finally, a preliminary control strategy for RO process is developed based on PID control algorithm and a first order transfer function (presented in the Appendix). / Government grant

This research focuses on large-scale manufacturing systems having a number of stations with multiple tools and product types with different and deterministic processing steps. The objective is to determine the production quantities of multiple products and the tool requirements of each station that maximizes net profit while satisfying strategic constraints such as cycle times, required throughputs, and investment. The formulation of the problem, named OptiProfit, is a mixed-integer nonlinear programming (MINLP) with the stochastic issues addressed by mean-value analysis (MVA) and queuing network models. Observing that OptiProfit is an NP-complete, nonconvex, and nonmonotonic problem, the research develops a heuristic method, Differential Coefficient Based Search (DCBS). It also performs an upper-bound analysis and a performance comparison with six variations of Greedy Ascent Procedure (GAP) heuristics and Modified Simulated Annealing (MSA) in a number of randomized cases. An example problem based on a semiconductor manufacturing minifab is modeled as an OptiProfit problem and numerically analyzed. The proposed methodology provides a very good quality solution for the high-level design and operation of manufacturing facilities.