Modeling, identification, estimation and adaptation for the control of power-generating kites

Large-scale kites, flying high-force crosswind trajectories, have been proposed for wind power generation. A two phase operational cycle generates net positive power using a ground-based motor/generator. In the traction phase the kite flies a high-force trajectory while reeling out the generator-connected tethers. A low-force retraction phase reels in the tethers and returns the kite to the start of the cycle. Highly variable conditions and significant uncertainty in the dynamics pose challenges to autonomous, well-controlled flight. The control task is divided into trajectory generation and tracking components and the most uncertain parameters in the model are identified online. The control structure uses these parameters in a robust framework resulting in an experimentally verified adaptive control scheme

A generic template for FSSIM

Agricultural and Food Policy, Environmental Economics and Policy, Land Economics/Use,

A generic template for FSSIM

The aim of this report is to describe the generic template of the Farm System Simulator FSSIM) designed to be applied for any farming systems across Europe for simulating farmer behaviour and assessing policy impacts. FSSIM is a bio-economic farm model which integrates biophysical processes, farm decision making and resource endowment. It consists of an agricultural management module (FSSIM-AM) and a mathematical programming model (FSSIM-MP). FSSIM-AM aims to identify current and alternative activities and to quantify their input output coefficients (both yields and environmental effects). FSSIM-MP seeks to describe the farmer’s behaviour given a set of biophysical, socio-economic and policy constraints and to predict his/her reactions under news technologies, policy and market changes

Produced water management - A mathematical model to trade-off economic cost and environmental impact for infrastructure utilisation

A substantial amount of wastewater, known as produced water (PW), is generated during oil and gas extraction. Given that PW can have a detrimental effect on the environment, it must be appropriately managed and treated before reuse. Globally, PW management is one of the greatest challenges in the oil and gas industry due to the costly treatment methods and large amounts involved, and there is a lack of expertise in the knowledge of PW management. Kuwait is a leading oil producer, and PW management poses a severe threat to the sustainability of Kuwait’s oil fields in terms of cost control and environmental safety. Here, life cycle and economic assessments are used to develop a mathematical framework for analysing trade-offs between the financial costs and environmental impacts of PW management operations. Specifically, a multi-objective mixed-integer linear programming framework is formulated for Kuwait Oil Company’s (KOC) PW supply chain management with different operational and regulatory constraints. A model solution for sustainable operations over the short, mid and long term that aligns with KOC’s strategic policy on PW management is presented. A global sensitivity analysis (GSA) also performed to further assess the economic and operational factors that influence KOC’s PW management. Finally, risk assessment is conducted to identify and evaluate risks associated with PW utilisation. Results indicate that treatment operations account for half of total system costs, and electricity consumption accounts for most of environmental impact, affecting the sustainability of the PW supply chain system most significantly. Moreover, KOC's PW supply chain system is impacted by a number of factors, including discount rates, electricity costs, and water treatment costs. In the case of PW utilisation, several types of risks may be posed that may negatively affect health, technology, the environment, and the economy. The findings of this study can be used to assess and guide PW supply chain management at KOC

Lifetime prediction for power converters

Renewable energy is developing rapidly and gaining more and more commercial viability. High reliability of the generation system is essential to maximize the output power. The power inverter is an important unit in this system and is believed to be one of the most unreliable parts. In the case of wind power generation, especially in off-shore wind, when the system reliability requirement is high, a technique to predict the inverter lifetime is invaluable as it would help the inverter designer optimize his design for minimal maintenance. Previous researchers studying inverter lifetime prediction, focus either at device level such as device fatigue damage models, or at system level which require experimental data for their selected device. This work presents a new method to estimate the inverter lifetime from a given mission profile within a reasonable simulation time. Such model can be used as a converter design tool or an on-line lifetime estimation tool after being configured to a real converter system. The key contribution of this work is to link the physics of the power devices to a large scale system simulation within a reasonable framework of time. With this technique, the system down time can be reduced and therefore more power can be generated. Also, the failure damage to the system is avoided which reduces the maintenance cost. A power cycling test is designed to gather the lifetime data of a selected IGBT module. Die-attach solder fatigue is found out to be the dominant failure mode of this IGBT module. The accuracy of widely accepted Miner’s rule, which accumulates damage linearly, is discussed and a nonlinear accumulation method is promoted to predict the lifetime of power inverters

Evolutionary learning and global search for multi-optimal PID tuning rules

With the advances in microprocessor technology, control systems are widely seen not only in industry but now also in household appliances and consumer electronics. Among all control schemes developed so far, Proportional plus Integral plus Derivative (PID) control is the most widely adopted in practice. Today, more than 90% of industrial controllers have a built-in PID function. Their wide applications have stimulated and sustained the research and development of PID tuning techniques, patents, software packages and hardware modules. Due to parameter interaction and format variation, tuning a PID controller is not as straightforward as one would have anticipated. Therefore, designing speedy tuning rules should greatly reduce the burden on new installation and ‘time-to-market’ and should also enhance the competitive advantages of the PID system under offer. A multi-objective evolutionary algorithm (MOEA) would be an ideal candidate to conduct the learning and search for multi-objective PID tuning rules. A simple to implement MOEA, termed s-MOEA, is devised and compared with MOEAs developed elsewhere. Extensive study and analysis are performed on metrics for evaluating MOEA performance, so as to help with this comparison and development. As a result, a novel visualisation technique, termed “Distance and Distribution” (DD)” chart, is developed to overcome some of the limitations of existing metrics and visualisation techniques. The DD chart allows a user to view the comparison of multiple sets of high order non-dominated solutions in a two-dimensional space. The capability of DD chart is shown in the comparison process and it is shown to be a useful tool for gathering more in-depth information of an MOEA which is not possible in existing empirical studies. Truly multi-objective global PID tuning rules are then evolved as a result of interfacing the s-MOEA with closed-loop simulations under practical constraints. It takes into account multiple, and often conflicting, objectives such as steady-state accuracy and transient responsiveness against stability and overshoots, as well as tracking performance against load disturbance rejection. These evolved rules are compared against other tuning rules both offline on a set of well-recognised PID benchmark test systems and online on three laboratory systems of different dynamics and transport delays. The results show that the rules significantly outperform all existing tuning rules, with multi-criterion optimality. This is made possible as the evolved rules can cover a delay to time constant ratio from zero to infinity based on first-order plus delay plant models. For second-order plus delay plant models, they can also cover all possible dynamics found in practice

An upwind cell centred finite volume method for large strain explicit solid dynamics in OpenFOAM

Cotutela Universitat Politècnica de Catalunya i Swansea UniversityIn practical engineering applications involving extremely complex geometries, meshing typically constitutes a large portion of the overall design and analysis time. In the computational mechanics community, the ability to perform calculations on tetrahedral meshes has become increasingly important. For these reasons, automated tetrahedral mesh generation by means of Delaunay and advancing front techniques have recently received increasing attention in a number of applications, namely: crash impact simulations, cardiovascular modelling, blast and fracture modelling. Unfortunately, modern industry codes in solid mechanics typically rely on the use of traditional displacement based Finite Element formulations which possess several distinct disadvantages, namely: (1) reduced order of convergence for strains and stresses in comparison with displacements; (2) high frequency noise in the vicinity of shocks; and (3) numerical instabilities associated with shear locking, volumetric locking and pressure checker-boarding. In order to address the above mentioned shortcomings, a new mixed-based set of equations for solid dynamics formulated in a system of first order hyperbolic conservation laws was introduced. Crucially, the new set of conservation laws has a similar structure to that of the well known Euler equations in the context of Computational Fluid Dynamics (CFD). This enables us to borrow some of the available CFD technologies and to adapt the method in the context of solid dynamics. This thesis builds on the work carried out by Lee et al. 2013 by further developing the upwind cell centred finite volume framework for the numerical analysis of large strain explicit solid dynamics and its tailor-made implementation within the open source code OpenFOAM, extensively used in industrial and academic environments. The object oriented nature of OpenFOAM implementation provides a very efficient platform for future development. In this computational framework, the primary unknown variables are linear momentum and deformation gradient tensor of the system. Moreover, the formulation is further extended for an additional set of geometric strain measures comprising of the co-factor of deformation gradient tensor and the Jacobian of deformation, in order to simulate polyconvex constitutive models ensuring material stability. The domain is spatially discretised using a standard Godunov-type cell centred framework where second order accuracy is achieved by employing a linear reconstruction procedure in conjunction with a slope limiter. This leads to discontinuities in variables at the cell interface which motivate the use of a Riemann solver by introducing an upwind bias into the evaluation of numerical contact fluxes. The acoustic Riemann solver presented is further developed by applying preconditioned dissipation to improve its performance in the near incompressibility regime and extending its range to contact applications. Moreover, two evolutionary frameworks are proposed in this study to satisfy the underlying involutions (or compatibility conditions) of the system. Additionally, the spatial discretisation is alternatively represented through a nodal cell centred finite volume framework for comparison purposes. From a temporal discretisation point of view, a two stage Total Variation Diminishing Runge-Kutta time integrator is employed to ensure second order accuracy. Additionally, inclusion of a global posteriori angular momentum projection procedure enables preservation of angular momenta of the system. Finally, benchmark numerical examples are simulated to demonstrate various aspects of the formulation including mesh convergence, momentum preservation and the locking-free nature of the formulation on complex computational domains.En aplicaciones prácticas de ingeniería que implican geometrías extremadamente complejas, el mallado requiere típicamente una gran parte del tiempo total de diseño y análisis. En la comunidad de mecánica computacional, la capacidad de realizar cálculos sobre mallas tetraédricas está siendo cada vez más importante. Por estas razones, la generación automatizada de mallas tetraédricas por medio de técnicas de Delaunay y frente avanzado han recibido cada vez más atención en ciertas aplicaciones, a saber: simulaciones de impacto, modelado cardiovascular, modelado de explosión y fractura. Por desgracia, los códigos en la industria moderna para mecánica de sólidos se basan normalmente en el uso de formulaciones tradicionales de Elementos Finitos formulados en desplazamientos que poseen varias desventajas: (1) menor orden de convergencia para tensiones y deformaciones; (2) ruido de alta frecuencia cerca de las ondas de choque; y (3) inestabilidades numéricas asociadas con el bloqueo a cortante, el bloqueo volumétrico y oscilaciones de presión. Con el fin de abordar estas deficiencias, se introduce un nuevo conjunto de ecuaciones para mecánica del sólido formulada como un sistema de leyes de conservación de primer orden basada en una formulación mixta. Fundamentalmente, el nuevo sistema de leyes de conservación tiene una estructura similar a la de las famosas ecuaciones de Euler en el contexto de la Dinámica de Fluidos Computacional (CFD). Esto nos permite aprovechar algunas de las tecnologías CFD disponibles y adaptar el método en el contexto de la Mecánica de Sólidos. Esta tesis se basa en el trabajo realizado en Lee et al. 2013 mediante el desarrollo de la estructura de volúmenes finitos centrados en celdas upwind para el análisis numérico de dinámica del sólido explícita en grandes deformaciones y su implementación específicamente diseñada dentro del software de código abierto OpenFOAM, ampliamente utilizado ámbito académico e industrial. Además, la naturaleza orientada a objetos de su implementación proporciona una plataforma muy eficiente para su desarrollo posterior. En este marco computacional, las incógnitas básicas de este sistema son el momento lineal y el tensor gradiente de deformación. Asimismo, la formulación se extiende adicionalmente para un conjunto adicional de medidas de deformación que comprenden el cofactor del tensor gradiente de deformación y el jacobiano de deformación, con el fin de simular modelos constitutivos policonvexos que aseguran la estabilidad del material. El dominio se discretiza espacialmente usando un marco centrado en células de tipo Godunov estándar, donde se consigue la precisión de segundo orden empleando un procedimiento de reconstrucción lineal junto con un limitador de pendiente. Esto conduce a discontinuidades en las variables en la interfase de la célula que motivan el uso de un solucionador de Riemann mediante la introducción de un sesgo contra el viento en la evaluación de flujos de contacto numéricos. El presente solucionador acústico de Riemann es posteriormente desarrollado aplicando disipación pre-condicionada para mejorar su rendimiento en el cercano pero incompresibilidad régimen y extender su gama a aplicaciones de contacto. Además, se proponen dos marcos evolutivos en este estudio para satisfacer las involuciones subyacentes (o condiciones de compatibilidad) del sistema. Además, la discretización espacial se representa alternativamente a través de un marco de volumen finito centrado en células nodales para fines de comparación. Desde el punto de vista de la discretización temporal, se emplea un integrador temporal de Runge-Kutta de dos etapas con Disminución de Variación Total para asegurar segundo orden de precision. Finalmente, se simulan ejemplos numéricos de referencia para demostrar varios aspectos de la formulación que incluyen convergencia de malla, conservación de momento y la naturaleza libre de bloqueo de la formulación en dominios computacionales complejos.Postprint (published version

An upwind cell centred finite volume method for large strain explicit solid dynamics in OpenFOAM

In practical engineering applications involving extremely complex geometries, meshing typically constitutes a large portion of the overall design and analysis time. In the computational mechanics community, the ability to perform calculations on tetrahedral meshes has become increasingly important. For these reasons, automated tetrahedral mesh generation by means of Delaunay and advancing front techniques have recently received increasing attention in a number of applications, namely: crash impact simulations, cardiovascular modelling, blast and fracture modelling. Unfortunately, modern industry codes in solid mechanics typically rely on the use of traditional displacement based Finite Element formulations which possess several distinct disadvantages, namely: (1) reduced order of convergence for strains and stresses in comparison with displacements; (2) high frequency noise in the vicinity of shocks; and (3) numerical instabilities associated with shear locking, volumetric locking and pressure checker-boarding. In order to address the above mentioned shortcomings, a new mixed-based set of equations for solid dynamics formulated in a system of first order hyperbolic conservation laws was introduced. Crucially, the new set of conservation laws has a similar structure to that of the well known Euler equations in the context of Computational Fluid Dynamics (CFD). This enables us to borrow some of the available CFD technologies and to adapt the method in the context of solid dynamics. This thesis builds on the work carried out by Lee et al. 2013 by further developing the upwind cell centred finite volume framework for the numerical analysis of large strain explicit solid dynamics and its tailor-made implementation within the open source code OpenFOAM, extensively used in industrial and academic environments. The object oriented nature of OpenFOAM implementation provides a very efficient platform for future development. In this computational framework, the primary unknown variables are linear momentum and deformation gradient tensor of the system. Moreover, the formulation is further extended for an additional set of geometric strain measures comprising of the co-factor of deformation gradient tensor and the Jacobian of deformation, in order to simulate polyconvex constitutive models ensuring material stability. The domain is spatially discretised using a standard Godunov-type cell centred framework where second order accuracy is achieved by employing a linear reconstruction procedure in conjunction with a slope limiter. This leads to discontinuities in variables at the cell interface which motivate the use of a Riemann solver by introducing an upwind bias into the evaluation of numerical contact fluxes. The acoustic Riemann solver presented is further developed by applying preconditioned dissipation to improve its performance in the near incompressibility regime and extending its range to contact applications. Moreover, two evolutionary frameworks are proposed in this study to satisfy the underlying involutions (or compatibility conditions) of the system. Additionally, the spatial discretisation is alternatively represented through a nodal cell centred finite volume framework for comparison purposes. From a temporal discretisation point of view, a two stage Total Variation Diminishing Runge-Kutta time integrator is employed to ensure second order accuracy. Additionally, inclusion of a global posteriori angular momentum projection procedure enables preservation of angular momenta of the system. Finally, benchmark numerical examples are simulated to demonstrate various aspects of the formulation including mesh convergence, momentum preservation and the locking-free nature of the formulation on complex computational domains.En aplicaciones prácticas de ingeniería que implican geometrías extremadamente complejas, el mallado requiere típicamente una gran parte del tiempo total de diseño y análisis. En la comunidad de mecánica computacional, la capacidad de realizar cálculos sobre mallas tetraédricas está siendo cada vez más importante. Por estas razones, la generación automatizada de mallas tetraédricas por medio de técnicas de Delaunay y frente avanzado han recibido cada vez más atención en ciertas aplicaciones, a saber: simulaciones de impacto, modelado cardiovascular, modelado de explosión y fractura. Por desgracia, los códigos en la industria moderna para mecánica de sólidos se basan normalmente en el uso de formulaciones tradicionales de Elementos Finitos formulados en desplazamientos que poseen varias desventajas: (1) menor orden de convergencia para tensiones y deformaciones; (2) ruido de alta frecuencia cerca de las ondas de choque; y (3) inestabilidades numéricas asociadas con el bloqueo a cortante, el bloqueo volumétrico y oscilaciones de presión. Con el fin de abordar estas deficiencias, se introduce un nuevo conjunto de ecuaciones para mecánica del sólido formulada como un sistema de leyes de conservación de primer orden basada en una formulación mixta. Fundamentalmente, el nuevo sistema de leyes de conservación tiene una estructura similar a la de las famosas ecuaciones de Euler en el contexto de la Dinámica de Fluidos Computacional (CFD). Esto nos permite aprovechar algunas de las tecnologías CFD disponibles y adaptar el método en el contexto de la Mecánica de Sólidos. Esta tesis se basa en el trabajo realizado en Lee et al. 2013 mediante el desarrollo de la estructura de volúmenes finitos centrados en celdas upwind para el análisis numérico de dinámica del sólido explícita en grandes deformaciones y su implementación específicamente diseñada dentro del software de código abierto OpenFOAM, ampliamente utilizado ámbito académico e industrial. Además, la naturaleza orientada a objetos de su implementación proporciona una plataforma muy eficiente para su desarrollo posterior. En este marco computacional, las incógnitas básicas de este sistema son el momento lineal y el tensor gradiente de deformación. Asimismo, la formulación se extiende adicionalmente para un conjunto adicional de medidas de deformación que comprenden el cofactor del tensor gradiente de deformación y el jacobiano de deformación, con el fin de simular modelos constitutivos policonvexos que aseguran la estabilidad del material. El dominio se discretiza espacialmente usando un marco centrado en células de tipo Godunov estándar, donde se consigue la precisión de segundo orden empleando un procedimiento de reconstrucción lineal junto con un limitador de pendiente. Esto conduce a discontinuidades en las variables en la interfase de la célula que motivan el uso de un solucionador de Riemann mediante la introducción de un sesgo contra el viento en la evaluación de flujos de contacto numéricos. El presente solucionador acústico de Riemann es posteriormente desarrollado aplicando disipación pre-condicionada para mejorar su rendimiento en el cercano pero incompresibilidad régimen y extender su gama a aplicaciones de contacto. Además, se proponen dos marcos evolutivos en este estudio para satisfacer las involuciones subyacentes (o condiciones de compatibilidad) del sistema. Además, la discretización espacial se representa alternativamente a través de un marco de volumen finito centrado en células nodales para fines de comparación. Desde el punto de vista de la discretización temporal, se emplea un integrador temporal de Runge-Kutta de dos etapas con Disminución de Variación Total para asegurar segundo orden de precision. Finalmente, se simulan ejemplos numéricos de referencia para demostrar varios aspectos de la formulación que incluyen convergencia de malla, conservación de momento y la naturaleza libre de bloqueo de la formulación en dominios computacionales complejos

Advances in foundation design and assessment for strategic renewable energy

In order to meet EU legislation on emissions, significant effort is being invested into the development of cost-effective renewable power generation technologies. The two leading technologies are solar and wind power because of their potential for the lowest levelised cost of energy and for showing a growth in installed capacity and technological development. Various research findings have suggested that significant cost savings in the capital expenditure of renewable energy projects can be made through the optimisation of their support foundations, the understanding of which has formed the main goal of the research. [Continues.