;
Oh+'0HP
$TU Delft Repository search results0TU Delft Repository search results (max. 1000)TU Delft LibraryTU Delft Library@?~5@?~5՜.+,0HPX`hp
x
WorksheetFeuilles de calcul
B=%r8X"1Calibri1Calibri1Calibri1
Calibri 83ffff̙̙3f3fff3f3f33333f33333.^ETU Delft Repositoryg b uuidrepository linktitleauthorcontributorpublication yearabstract
subject topiclanguagepublication type publisherisbnissnpatent
patent statusbibliographic noteaccess restrictionembargo datefaculty
departmentresearch group programmeprojectcoordinates)uuid:56a2f8d7fb2640f0bb9200271f9492dfDhttp://resolver.tudelft.nl/uuid:56a2f8d7fb2640f0bb9200271f9492df7Optimizing the gradient ring of a lowfield MRI scanneroMeijer, Angeline (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Applied Sciences)tvan Gijzen, Martin (mentor); van Dongen, Koen (mentor); Delft University of Technology (degree granting institution)This report describes the study on the optimization of the magnetic field generated by a lowfield magnetic resonance imaging modality. In the optimization the Genetic Algorithm and a quasiNewton Algorithm are used.;Optimization; MRI; lowfield; Magnetic field; Gradient ringenbachelor thesis%Applied Mathematics  Applied Physics
LowField MRI)uuid:6818b21f038d4700afcde6dac79e7e82Dhttp://resolver.tudelft.nl/uuid:6818b21f038d4700afcde6dac79e7e82`On The Cap Set Problem: upper bounds on maximal cardinalities of caps in dimensions seven to tenqVersluis, Nina (Electrical Engineering, Mathematics and Computer Science; Delft Institute of Applied Mathematics)Gijswijt, Dion (mentor); van Elderen, Emiel (graduation committee); de Groot, Joost (graduation committee); Delft University of Technology (degree granting institution)'This thesis concerns the cap set problem in affine geometry. The problem is illustrated by the card game SET and its geometrical interpretation in ternary affine space. The maximal cardinality of a cap is known for the dimension one to six. For the four lowest dimensions, a maximal cap is constructed and the optimality of its size proven. From there, two recursive methods are described and applied to obtain upper bounds for the maximal size of caps in dimensions seven to ten. The best found upper bounds are 291, 771, 2070 and 5619, respectively.*Affine caps; Affine geometry; Optimization)uuid:5ecc206d60bb43ef9874ee72c6263c6aDhttp://resolver.tudelft.nl/uuid:5ecc206d60bb43ef9874ee72c6263c6a2Three Way Duels: Infinite Games on the Unit SquareLHa, Deon (TU Delft Electrical Engineering, Mathematics and Computer Science)WFokkink, Robbert (mentor); Delft University of Technology (degree granting institution)With the growing wealth and economy of a country, there are an increasing amount of small and big businesses. Every company has its own marketing strategy that it uses in order to lure customers away from their competition and increase their sales. Choosing the perfect time to advertise or discount several products is of essence for a company to gain more money than their competition. These type of marketing games are all slight variations of duels. The purpose of this report is to research how this duel is played most optimal when there are two or more participants. Several types of twoplayer duels shall be analysed first in order to understand and analyse a threeplayer duel.>Game Theory; Optimization; Differential Equations; ProbabilityApplied Mathematics)uuid:0fc3b7460d3f459d8fc2c9e6aea821dfDhttp://resolver.tudelft.nl/uuid:0fc3b7460d3f459d8fc2c9e6aea821dfTreechild Network ContainmentUHuijsman, Robbert (TU Delft Electrical Engineering, Mathematics and Computer Science)van Iersel, Leo (mentor); Murakami, Yuki (mentor); Janssen, Remie (mentor); van den Dries, Bart (graduation committee); Delft University of Technology (degree granting institution)Interest in phylogenetic trees for histories of species and DNA has spawned many problems, one of which is TreeContainment; a problem that asks whether a tree is contained within a network. The TreeContainment problem is proven to be NPhard for general trees and networks, however it is solvable in polynomial time for networks that meet the treechild restriction. An algorithm to solve TreeContainment for binary treechild networks has been created previously with quadr< atic running time (van Iersel, Semple, Steel, 2010). Janssen and Murakami have recently created a new algorithm that solves a larger problem NetworkContainment, for semibinary treechild networks (Janssen, Murakami, 2019). This new algorithm uses treechild sequences introduced by Linz and Semple, but there has not been an implementation of it until now. In this paper I show an implementation (using Python) of this algorithm, in which I have made a modification that increases its speed on networks with large indegrees. Furthermore I have proven in this paper that the output of this algorithm remains correct under this modification, and that the running time of the modified algorithm is now linear without requiring a constant maximum indegree at all.=Phylogenetic Network; Algorithm; Implementation; OptimizationApplied Physics)uuid:b7ccbf2752cf4d50ba649282aec248e4Dhttp://resolver.tudelft.nl/uuid:b7ccbf2752cf4d50ba649282aec248e4RSimulation Verification and Optimization of a Vertical Axis Wind Turbine using CFD
Kortleven, M.Kenjeres, S. (mentor); Heemink, A.W. (mentor)NIn this research computational fluid dynamics (CFD) is used to model a vertical axis wind turbine. Well known turbulence models  i.e. the kepsilon and komega SST model  are used and compared in their performance against each other and a set of experimental data. The formal research question is "What models are available for the simulation of vertical axis wind turbines and what considerations should be taken into account to most effectively obtain the power coefficient of those turbines?" The results show that the 2D kepsilon simulations approximate the experimental data the closest, but that the difference between the simulations and the experiment in the komega model can be explained by the simulations being 2D and the experiment being 3D. To improve the results and verify that a 3D simulation indeed does produce better results, another study should be conducted in 3D when more computation power is readily available. As to what considerations should be taken into account: Meshing determines largely how long a simulation will take and to what accuracy a result can be calculated. The sliding mesh method should be used to decrease calculation time, so that the mesh does not need to be recalculated every time step. The orders of the different variables should be set to second order accuracy for the pressure p, the momentum, the turbulent kinetic energy k and the time derivative. The second order calculation of these variables results in a significant increase in accuracy of the simulations relative to simulations conducted in first order accuracy. For further improvement a research could be set up that involves both the experimental and the CFD part. That way the practical limitations to building a wind turbine could be modelled in the simulations, possibly resulting in even better synergy between the measurements and the simulations.9CFD; Simulation; Vertical Axis Wind Turbine; OptimizationApplied SciencesChemE/Chemical EngineeringTransport Phenomena)uuid:437ca7513fd642c3933f72771f025eeeDhttp://resolver.tudelft.nl/uuid:437ca7513fd642c3933f72771f025eee]Diglimit Optimization: Mixed Integer Linear Programming for grade control in openpit mining>Buist, Jan Willem (TU Delft Civil Engineering and Geosciences)uWambeke, Tom (mentor); van Duijvenbode, Jeroen (mentor); Delft University of Technology (degree granting institution).In openpit mining on a bench level, the diglimit optimization problem is deciding whether to classify a Selective Mining Unit(SMU) as waste or whether to classify it as ore. SMU s with high ore grades are protable for the mining operation, but due to equipment size a whole frame of SMU s has to be classied as either ore or waste. In this thesis the diglimit optimization problem will be solved using Mixed Integer Linear Programming (MILP). The model proposed in this thesis takes the frame constraint and the average grade constraint into account. The main results of the MILP model are in table 1. The results show that M< ILP can successfully be used to determine the optimal value of a mining operation, however the determination of the diglimit for larger grids leads to exponential computation time.<brYopenpit mining; Mixed integer linear programming; grade control; diglimit; Optimization)uuid:6b8900808a60418989dcb12f23a8677aDhttp://resolver.tudelft.nl/uuid:6b8900808a60418989dcb12f23a8677aOrienting phylogenetic networksXVerzijlbergen, Esmay (TU Delft Electrical Engineering, Mathematics and Computer Science)van Iersel, Leo (mentor); Jones, Mark (graduation committee); Keijzer, Marleen (graduation committee); Delft University of Technology (degree granting institution)1 This thesis provides you with basic information on graph theory as well as phylogenetic networks, it studies the relationship between undirected (unrooted) and directed (rooted) phylogenetic networks, based on the manuscript 'Rooting for phylogenetic networks' . Undirected phylogenetic networks can be oriented to become a directed network. In this manuscript the authors come up with multiple algorithms for orienting phylogenetic networks meeting different characteristics. A network is built up of reticulation vertices (where lineages merge) and tree vertices (where lineages separate). A network can be binary, meaning every node in the body of the network has a degree of three or a network can be nonbinary, meaning there are no restrictions to the amount of edges a vertex can have. <br/>When a network is binary, an algorithm is described that given the location of the root as well as a set of reticulation points, is used to find an orientation (Algorithm 1). If a network is nonbinary, an algorithm is described that given a location of the root as well as the indegree of each vertex, is used to find an orientation (Algorithm 2 ). Once an orientation is found for a certain undirected network this orientation can be checked to see whether it meets the characteristics of a certain network class. Three network classes are considered. First of all an orientation can be treechild, meaning that every nonleaf vertex has a child that is not a reticulation. Secondly an orientation can be stackfree, meaning that no reticulation has a child which is a reticulation. And last, an orientation can be valid, meaning that it is stackfree and deleting a single reticulation edge and suppressing its endpoints does not give parallel arcs. When we either did not find an orientation in the class we wanted or want to know all the orientations for a certain network in a certain class, we use Algorithm 3 or 4. The location of the root and the set of reticulation points were necessary input for Algorithm 1; whereas Algorithms 3 and 4 do not require this information. Algorithm 3 returns you the first found orientation in the class you wanted and Algorithm 4 returns you a collection of all orientations in the class. All orientations found by the program are returned in an output format that can be read by a network visualisation website.5Optimization; Phylogenetic network; Orientation model)uuid:0126e698df974313a4d178f4c3f8c6d0Dhttp://resolver.tudelft.nl/uuid:0126e698df974313a4d178f4c3f8c6d03Dynamic Response Optimization of an Acoustic Guitarvan Boven, M.M.A.Langelaar, M. (mentor)m Currently the production process of acoustic guitar factories does not result in instruments with equal dynamical behaviour, due to the varying material properties of wood from tree to tree. Experienced individual guitar builders are able to handle those material variations, but this is a labour intensive process because there is no structured methodology available about how to obtain a certain desired frequency response. In this study it is investigated if the perceived sound produced by a guitar is controllable by using dynamic response optimization methods in order to adjust geometric design parameters during the production process. The goal is to find a method which is both able to improve the consistency between massproduced instruments and to decrease the required building time of handbuilt guitars. Unt< il now no research is available which uses optimization methods in order to shape the dynamic response of a guitar over a range of frequencies. Previous research presented an analytical model which is able to give a quantitative description of a guitar s dynamic response. The response is characterized by a measurement of the sound pressure level at 1 meter distance over a range of excitation frequencies. It was shown that the sound of a guitar is mainly determined by the dynamic behaviour of its top plate. This study uses a finite element model in order to estimate the effect of geometric design variations. The finite element model is used in combination with an optimization algorithm which shapes the dynamic response to a certain desired response. Due to this optimization approach the finite element model needs to be evaluated for many combinations of design variables. A computationally cheap finite element model is needed to approximate the physical behaviour of the guitar s top plate. Based on the results of the optimization algorithm it is concluded that it is indeed possible to control the perceived sound produced by a guitar, by using dynamic response optimization methods in order to adjust geometric design parameters during the production process. Although some aspects need to be worked out further before the approach presented in this research can be implemented in a production process, the theoretical background will definitely be useful in order to improve consistency between instruments by using frequency response measurements during the building process.Optimization; Guitar
master thesis.Mechanical, Maritime and Materials Engineering&Precision and Microsystems Engineering)uuid:3e54884398e74a18a50e017309ae95d8Dhttp://resolver.tudelft.nl/uuid:3e54884398e74a18a50e017309ae95d8RThe optimization of IceWorlds' wooden boarding system for recreational ice rinks.:van Erven, Jeroen (TU Delft Industrial Design Engineering)van de Geer, Stefan (mentor); Kroon, Caroline (graduation committee); Delft University of Technology (degree granting institution)Ice skating is (still) very popular in the Netherlands, and it is gaining popularity across the globe. Temporary ice rinks on squares and in parks is an increasing solution in order to attract visitors in the winter and holiday season. IceWorld International, the world leader in temporary ice rinks since 1992, has experienced a significant increase in demand for temporary ice rink rental in the Netherlands. As a premium supplier, IceWorld wants to keep offering its customers high quality and innovative solutions. Therefore, IceWorld asked to explore opportunities for updating the boarding systems with a new modular boarding system that fits every rink dimension and could replace current systems. This graduation project has resulted in a redesign of the wooden boarding system. A boarding system is used as a barrier around an ice rink to prevent ice skaters from sliding or falling off the ice. It is a modular system, as every ice rink has different dimensions. Currently, Ice World carries two different types of boarding systems: a transparent boarding and a wooden boarding. The wooden boarding is preferred by European customers, while international customers prefer transparent boarding. Existing problems with the boarding especially relate to the wooden system. The focus of this project was set to the wooden boarding because of increasing market demand (within the Benelux), and it is the core business of IceWorld.fIce rink boarding; Product development; Sustainability; Prototyping; IceWorld; Optimization; RedesignIntegrated Product design)uuid:3f3e955b36bc4a09a9ccb7065a59d39bDhttp://resolver.tudelft.nl/uuid:3f3e955b36bc4a09a9ccb7065a59d39bIOptimal model distribution in multi model adaptive control with switching Ghosh, D.Baldi, S. (mentor) Multi model adaptive control is an emerging field that has proven to be successful in mitigating limitations of classical adaptive control by providing more quicker and faster adaptation. In this method instead of one si< ngle parameter varying controller, multiple fixed parameter controllers pertaining to different operation regimes are utilized. A variant of this technique known as the multi model unfalsfied adaptive control is a recently proposed method. This methodology although promising needs further systematic analysis regarding the determination of the minimum number of controllers required and their corresponding locations (in the uncertain parameter space) to address stability and performance issues. This research work presents the basics of Adaptive Control with Multiple Models, introduces the problem of optimal problem distribution of the multiple models, and finally demonstrates a novel optimization based methodology to determine the number of controllers and their optimal locations for a given uncertainty.6adaptive control; multi model; switching; Optimization$Delft Center for Systems and Control,Mechanical Engineering  Control Engineering)uuid:f9a658fe9d0b4a1d808a908698c8e8c2Dhttp://resolver.tudelft.nl/uuid:f9a658fe9d0b4a1d808a908698c8e8c2@Analysis of an over the wing based distributed propulsion system3Khajehzadeh, Arash (TU Delft Aerospace Engineering)oVeldhuis, Leo (mentor); Hulshoff, Steven (mentor); Delft University of Technology (degree granting institution)This project aims to investigate propeller wing interaction of an over the wing (OTW) based distributed propulsion (DP) system with the addition of a secondary wing. This research explored the opportunities that OTW configuration provides to DP concept, such as improving the lifting performance of system by inducing the flow over the upper surface of wing. Secondary wing is oriented in a biplane configuration above propeller and intends to improve propeller s performance by decreasing propeller s inflow velocity. This research aimed to investigate the influence of shape and position modification on the system s aerodynamic performance and overall propulsive efficiency. The influence of spacing between the wings, position of propeller, secondary wing s angle of attack and secondary wing s initial lift coefficient on aerodynamic performance and overall propulsive efficiency of DP system was investigated in this project.<br/>The analysis of this particular configuration was performed with the help of Euler calculations. Eliminating viscous effects from the analysis, reduced the demanded computational cost of this study and could help this research to perform a broader investigation. The influence of isothermal flow condition on Euler calculation was examined to reduce the computational cost of optimization study; This study showed that the drag coefficient of system is sensitive to this assumption and that only by simulating the initial and final points of the optimization study with adiabatic flow condition, the use of isothermal flow condition could be justified. The secondary wing s shape was optimized to study the influence of its shape variation on the system s aerodynamic performance. The first optimization study aimed to decrease the drag coefficient of system, and as a result, overall propulsive efficiency of was improved. The second optimization study aimed to improve the lifting performance of system. As a result, the drag coefficient of system significantly increased, however, the overall propulsive efficiency was again improved. The flow decelerated between main wing and secondary wing since the lifting performance of secondary wing improved, which increased the propeller s thrust according to blade element method (BEM) calculation.<brpropeller wing interaction; numerical analysis; over the wing; Optimization; CFD; cfx; distributed propulsion; actuator disc)uuid:87d1d976a0ba40d28fa4437bd65c5b91Dhttp://resolver.tudelft.nl/uuid:87d1d976a0ba40d28fa4437bd65c5b91.Transverse 3D Electrical Properties TomographyjAn, Xin (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Circuits and Systems)QRemis, Rob (mentor); Delft University of Technology (degree granting institution)MRbased Electrical Properties Tomography (EPT) is a me< dical imaging modality to reconstruct the electrical conductivity and permittivity inside the human body based on the magnetic field data obtained from an MRI system. Currently, 2D CSIEPT implementation is able to generate reliable reconstructions under special assumption of electromagnetic fields, while 3D implementation of CSIEPT is able to achieve accurate reconstruction. However, much effort is needed to reduce computational cost for practical application in a MRI setting. Several experiments show that the radiofrequency(RF) fields on the midplane of an birdcage RF coil is Epolarized. In this study, we assume that RF field is Epolarised in the midplane of the RF birdcage coil and no other field assumptions are imposed. The electromagnetic fields relations governed by Maxwell s equations in the birdcage RF coil are simplified by the negligible<br/>transverse electric fields into two equations, referred to as midplane equations. In order to find the parameters that satisfy the midplane equations, an objective function is created by evaluating the difference between the measured and modeled data. A novel alternatingdirection iterative method is proposed for minimizing the multivariate objective function where the electrical properties of biological tissues are retrieved in an iterative fixoneoptimizefortheother way. Within each iteration, a twostep procedure is conducted by updating electric field parameter and dielectric parameter in turn. Three different optimization algorithms are implemented for updating the electric field parameter. The magnetic fields data for numerical simulations are generated by placing a male head model or a longitudinal uniform head model inside<br/>a simulated magnetic resonance scanner. The proposed methods are demonstrated by reconstructions from the simulated data, which are in good agreement with true electrical properties values.<brJElectrical Property Tomography; Optimization; Alternating Direction Method
20240827)uuid:7fe5a51181c944e590baa2b479011077Dhttp://resolver.tudelft.nl/uuid:7fe5a51181c944e590baa2b479011077fRobust optimal sensor placement for optimal estimation: A framework and application to adaptive opticsyde Groot, Sjoerd (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Delft Center for Systems and Control)Verhaegen, Michel (mentor); Doelman, Reinier (graduation committee); Ferrari, Riccardo (graduation committee); Delft University of Technology (degree granting institution)rThe goal of this thesis is to present a framework which can be used to optimize the observer and controller performance by placement of sensors and actuators on a continuous domain. The framework provides choices of model structure, cost functions and optimization methods to perform the placement. Systems which would be suitable for the use of the framework are for example adaptive optics systems, active dampening of windmill blades or climate control systems in buildings. To demonstrate the use of the framework, two applications are presented. Both applications are sensor placement for an adaptive optics system. The first implementation uses a Linear Time Invariant (LTI) state space model and uses a Vector Autoregressive (VAR) approach for identification. The second implementation uses a Linear Parameter Varying (LPV) state space model with VAR identification approach. The second implementation is meant to demonstrate how the placement can be made more robust to changing conditions. Finally we show how parametrization of the sensor locations can be used to gain insight into what properties an optimal placement exhibits.&Control; Optimization; Adaptive OpticsSystems and Control)uuid:78ae5715d0794a7e983cd999ba761c16Dhttp://resolver.tudelft.nl/uuid:78ae5715d0794a7e983cd999ba761c16Towards the Automation of Preliminary Bridge Designs with Parametric Design Software: Assessment and optimization of multiple steelconcrete bridge typologies<Fransen, Edward (TU Delft Civil Engineering and Geosciences)Nijsse, Rob (graduation committee); Eigenraam, Peter (graduation com< mittee); Smits, Joris (graduation committee); Freriks, W.G.M. (mentor); Delft University of Technology (degree granting institution)This thesis concerns research into the parametric modelling and optimization of steelconcrete bridge typologies typically found in the Netherlands, in order to further the automation of preliminary bridge designs. The research question posed is: "How can a parametric bridge design model be composed to fit into the <br/>overall road design process in order to get a comparative overview of structurally sound bridge designs of different typologies?". Within the literature review some background information is provided about the design of steelconcrete bridges, parametric design, finite element modelling and evolutionary solvers for optimization. With this information seven different bridges are designed according to boundary conditions and requirements posed in a fictitious case. These bridge designs cover the beam, truss, arch and cable stayed<br/>bridge typologies. The designs are then parametrically modelled and converted to a finite element model in a single parametric design environment, after which they are structurally analysed and optimized for minimum costs. The result is that, for the fictitious case considered, the Warren truss without verticals is the<br/>solution with the lowest costs. This bridge design is able to carry the loads, according to load model 1 from the Euro Codes,with a lower amount of structural steel compared to the girder beam bridges, and the absence of expensive cables in the structure makes it less expensive than the arch and cable stayed bridges. Due to the scope defined at the start of the research, there are some limitations. The costs calculated for the bridges are purely based on the bridge structure itself, foundations and approach structures are not considered and this has to be kept in mind when comparing the bridges. Finally, the research into how to implement this design and optimization tool into the overall design process has not been performed, however, the model is composed in such a way that it will be compatible with minor changes to the start and end of the process.EParametric design; Steel; Bridges; Concrete; Optimization; Automation
20190221Building Engineering)uuid:056d28defe7e48ae8a9d0af5bd9aa2c0Dhttp://resolver.tudelft.nl/uuid:056d28defe7e48ae8a9d0af5bd9aa2c0UOptimization of User Equilibrium container transportation problems using toll pricingnBruijns, Lianne (TU Delft Applied Sciences; TU Delft Electrical Engineering, Mathematics and Computer Science)Aardal, Karen (mentor); Dubbeldam, Johan (graduation committee); Phillipson, F. (graduation committee); Sangers, A (graduation committee); Delft University of Technology (degree granting institution)VThe goal of this thesis is to provide a method to obtain a User Equilibrium (UE) in a synchromodal transportation network in which we transport containers for multiple customers. <br/>We use the Space Time Network (STN), in which the travel time of modalities is fixed and independent of the occupancy of the network.<br/>We define the User Equilibrium solution as the solution wherein each customer can travel via their cheapest paths possible, and no customer is harmed by the route choice of other customers. We define the System Optimal (SO) solution as the solution in which the total costs of the network are minimized.<br/>We expand the goal of finding a UE on the STN to finding a solution where the solution is both SO as well as UE.<br/>The solving method we developed in this thesis, consists of finding SO solutions by solving the Minimum Cost MultiCommodity Flow (MCMCF) problem on the STN to find an SO solution, and create tolls schemes to create a UE solution, while maintaining the SO solution.<br/>We investigate several types of tolls in order to obtain sufficient both SO and UE solutions.,Optimization; User Equilibrium; Toll Pricing)uuid:f3c8e60a65d245e1a838001a5c21ccebDhttp://resolver.tudelft.nl/uuid:f3c8e60a65d245e1a838001a5c21ccebContesting Conservatism with Complexity: < Exploring Linear Relaxations of the Positivstellensatz for Robust Stability Analysis of Nonlinear SystemsRFischerGundlach, Moritz (TU Delft Mechanical, Maritime and Materials Engineering)pWolfslag, Wouter (mentor); Wisse, Martijn (mentor); Delft University of Technology (degree granting institution)iLyapunov's 2nd method can be formulated as a convex optimization problem by means of SumofSquares (SOS) optimization. This gives a powerful tool for robust stability analysis of nonlinear systems. The use of SOS optimization is limited to small dynamical systems, because of the demanding numerical complexity underlying SOS optimization. For systems where scaling of SOS is a concern, diagonally dominant SOS (DSOS) is a possible alternative. DSOS approximates the SOS problem from within and results in a less complex optimization problem. The shortcoming of DSOS, however, is conservatism within the robust stability analysis, entailed by the approximation.<br/><br/>In this thesis, we consider the problem of conservatism in robust stability analysis due to the application of DSOS. We propose new DSOSbased formulations which aim to reduce conservatism in robust stability analysis. Furthermore, we show how to formulate optimization problems which don't require an initial Lyapunov function candidate, and stabilize the numerical solutions of the DSOS problems. The main tool we leverage in this study is the Positivstellensatz in combination with DSOS. <br/><br/>We demonstrate the proposed formulations on three numerical experiments. We estimate the region of attraction of the van der Pol oscillator, investigate stability of an uncertain nonlinear system and proof stability of a nonlinear system, which is documented to suffer from numerical issues when approached with SOS methods. The results show a reduction in conservatism and improvement in numerical robustness when compared to the DSOS formulation found in literature.%SOS; Stability Analysis; Optimization
20180825"Biomechanical Design / Biorobotics)uuid:232d91048cd34ace89b3042e50490871Dhttp://resolver.tudelft.nl/uuid:232d91048cd34ace89b3042e50490871Aircraft Trajectory Optimization for European Flights: Assessment of the trajectory variability considering a tradeoff between operating costs and climate impact.Blom, Martijn (TU Delft Aerospace Engineering)TGrewe, Volker (mentor); Delft University of Technology (degree granting institution)In order to reduce the climate impact induced by the aircraft industry, the variability between cost optimized and climate optimized flight trajectories was evaluated. This was done to establish how the climate impact reduction affects the operating costs and vice versa. To evaluate the climate impact of air traffic, a climatechemistry model was used to simulate the earth its atmosphere and submodels to study aircraft trajectories. This study considered 85 flights spread out over the European airspace, analyzed for six months, including realistic weather systems throughout the entire analysis spectrum. Cost and climate optimized trajectories were compared for flight altitude, horizontal flight paths, emissions and operating costs. The presence of a seasonal effect was examined which marks the difference between a winter and a summer season, and four different flight directions were considered that showed the effect of the wind speed and direction on the flight properties.BTrajectory; Optimization; Climate; Impact; Operating; Costs; ATM4E)uuid:841a9a8dc7d44bdba1a47a4e7960070bDhttp://resolver.tudelft.nl/uuid:841a9a8dc7d44bdba1a47a4e7960070bGContinuoustime System Identification: A Bilinear Optimization Approach{Krijgsman, Wiegert (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Delft Center for Systems and Control)sVerhaegen, Michel (mentor); Doelman, Reinier (mentor); Delft University of Technology (degree granting institution)NIn the noiseless case, the identification of a greybox model can be posed as a feasibility problem, i.e. determining if existent  and if so  finding a parameter vector such that the para< metric model equals the actual model (or its associated inputoutput data). Being that in this thesis we are interested in continuoustime greybox identification, we shall be dealing with models that allow for forming a direct relationship with physical meaningful quantities. Such models include the state space representation and the matrix differential equation. In general, identifying such greybox models turns out to be a nonconvex problem. In this thesis, we initially review a framework which allows us to solve feasibility problems which have bilinear constraints. It turns out that most of the aforementioned nonconvexities can be captured into a single bilinear matrix equation. However, the resulting feasibility problem, including the bilinear matrix equations, makes the overall search for the actual parameter vector NPhard. In order to come up with numerical tractable algorithms, we use a heuristic known as Sequential Convex Relaxation to relax the bilinear equality constraints. This iterative scheme is flexible enough to allow for additional (in)equality constraints, possibly resembling any other physical constraints. We explore two different approaches to identify both the state space model and the matrix differential equation; one, by directly identifying the model from the given frequency response function; two, by first identifying a blackbox model before performing a small scale optimization problem, transforming the blackbox model such that it fits the greybox parameterization. In addition, we present a novel method which uses the Power Spectral Density to estimate a 2nd order model. All methods are numerically validated.<br/><br/><brOptimization; GreyBox; Bilinear; statespace models; matrix differential equation; Frequency Response Function; parametric approach; physical modelling; Identification; power spectrum
20220101)uuid:00eb0f82521346578fa42630c3594463Dhttp://resolver.tudelft.nl/uuid:00eb0f82521346578fa42630c35944639Decentralized coordination for area surveillance purposesMDemetriou, Kyriakos (TU Delft Mechanical, Maritime and Materials Engineering)Sijs, Joris (mentor); Fransman, Jeroen (mentor); De Schutter, Bart (graduation committee); Ferrari, Riccardo (graduation committee); Delft University of Technology (degree granting institution)The area surveillance problem is the problem of surveying a known or an unknown area with the main purpose of detecting objects. This thesis will tackle the problem of how to employ a group of mobilesensors for surveying an unobstructed area in an optimal manner. The mobilesensors should make use of their onboard computers and they should iteratively compute their waypoints until they successfully surveyed the area. To solve the area surveillance problem in an optimal manner, the mobilesensors should be able to coordinate their actions. The Distributed Constraint Optimization Problem (DCOP) framework will be employed to define the area surveillance problem. In the DCOP, the mobilesensors can define their optimal position for the next discrete time step by communicating with the other mobilesensors that participate in the problem.<br/>For solving this DCOP, the Distributed Pseudotree Optimization Procedure (DPOP) will be utilized. The DPOP is a complete solver that can find the optimal solution of a DCOP in a decentralized manner. However, the main limitation of DPOP is that the size of the largest message that the mobilesensors will have to exchange is spaceexponential in the induced width of the pseudotree. Considering that the mobilesensors have to use their onboard computers for solving the problem, exchanging huge size of messages is not desired. To overcome this limitation of DPOP, a new extension, known as the MSTDPOP will be presented in this Thesis. The MSTDPOP makes use of the Maximum Spanning Tree (MST) algorithm to reduce the size of the largest message that the mobilesensors have to exchange. Employing MST along with DPOP can bound the size of the largest message and the required computations for constructing the utility messages. Ho< wever, the new extension cannot guarantee that its solution will be the optimal. The experimental results shows that the MSTDPOP is able to define a solution with an average error less than $2\%$. Moreover, the MSTDPOP requires on average around $1$ discrete time step more than the DPOP to solve the area surveillance problem. Consequently, given that the MSTDPOP overcomes the high memory requirements of DPOP, it is preferred for solving the area surveillance problem.4Optimization; DCOP; DPOP; Decentralized coordination,Mechanical Engineering  Systems and Control)uuid:62636a6807ba424480362f4f9bc133a1Dhttp://resolver.tudelft.nl/uuid:62636a6807ba424480362f4f9bc133a1WDesign Optimization of VariableStiffness Composite Panels under Thermomechanical LoadsGvan Loo, Rick (Aerospace Engineering; Aerospace Structures & Materials)Fatemi, Javad (mentor); Abdalla, Mostafa (mentor); Peeters, Danil (mentor); Bergsma, Otto (graduation committee); Turteltaub, Sergio (graduation committee); Delft University of Technology (degree granting institution)Striving to improve structural efficiency the aerospace industry shows increasing interest in variablestiffness composite laminates. Advanced fiber placement is a hybrid manufacturing technique that offers the flexibility of both filament winding and automated tape laying. With the development of this novel system curved tows can be placed and a spatially variablestiffness laminate can be designed with continuous changing stiffness from point to point. <br/><br/>The increased design freedom to tailor a structure by inplane stiffness variation leads to a challenging design optimization problem. A multistep framework is developed by the aerospace structures and materials department to optimize variablestiffness laminates. Variablestiffness laminate design allows for sophisticated designs. Based on this premise it is investigated how such design could improve the structural performance of an engine thrust frame, a structural application that transfers the thrust loads from the rocket engine to the rest of the launch system. The engine thrust frame is subject to cryogenic thermal loads, something not incorporated in the available optimization framework. The goal of this work is to add thermal loads to the laminate analysis routine and to adjust the optimization routine to incorporate thermal influences.<br/><br/>With the thermomechanical optimization framework in place the engine thrust frame is modeled. Conceptual design optimization of the engine thrust frame under thermomechanical loads is performed to increase buckling resistance. A mismatch in the coefficient of thermal expansion is used by the optimal variablestiffness design. The stiffened areas contract less than the interstiffener bay regions. Consequently a stabilizing tensile stress is induced in the prone to buckling bay regions, whereas compressive stresses are distributed to the stiffened areas. Based on the stabilizing thermal stresses and load distribution significant gains in performance are found. <br^Thermomechanical; VariableStiffness; Composite; Laminate Analysis; Optimization; Thermal Load
20220621)uuid:0ad5f3cca5e945c0a37e75e139713409Dhttp://resolver.tudelft.nl/uuid:0ad5f3cca5e945c0a37e75e1397134099Thermal Morphology: A Geometrically Optimized Trombe Wall?Farrugia, Eve (TU Delft Architecture and the Built Environment)rTenpierik, Martin (mentor); Turrin, Michela (mentor); Delft University of Technology (degree granting institution)Recently, there has been a growing interest in passive climate design aimed at reducing the large energy demand arising from existing buildings. Sustainable consumption of quantities of energy by buildings has necessitated research into innovative designs for comfortable and energyefficient buildings. This graduation project branches out from this research field and focuses on experimentation with geometry and its effects on thermal mass in the context of architecture.<br/><br/>It is known that geometry has an impact on thermal performance, however information regarding < alternatives for material distribution to improve thermal performance is somewhat limited. The main objective of this research is to identify how geometry could influence and enhance the thermal performance to reduce the overall energy consumption in a specific context. This graduation project focuses on thermal mass in the form of a trombe wall in office building typologies. It also targets the use of translucent Phase Change Materials (PCM) as a thermal energy storage material for lightweight building components in the context of two different climates; a hot and dry climate and a temperate climate.<br/><br/>A combination of simulation and modelling software is used to determine the effect of geometry on the thermal performance of the trombe wall with PCM in two different climates. A stepped methodology is used to morph a flat wall following a series of design strategies and two different geometry scales. By analysing the differences in air and PCM temperatures it is determined whether geometrical morphology is beneficial in terms of a comfortable air temperature in an office.<br/><br/>General trends are identified which highlight the potential of material distribution in thermal performance in the context of PCM and trombe walls. These trends are then transformed into design guidelines and two different designs for a trombe wall with translucent PCM in the two climates.mMorphology; Thermal Performance; PCM; Trombe Wall; Optimization; Energy Efficiency; Thermal comfort; Geometry3Building Technology  Sustainable Graduation Studio)uuid:5137791874be4ec4bad2115276e35cd8Dhttp://resolver.tudelft.nl/uuid:5137791874be4ec4bad2115276e35cd8MTensegrity Full Body Control: An Inequality Constrained Optimization ApproachrGudmundsson, Kjartan (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Biomechanical Engineering)van der Weijde, Joost (mentor); Vallery, Heike (graduation committee); Peternel, Luka (graduation committee); Delft University of Technology (degree granting institution)Tensegrity is a structural form that is defined as a set of rigid elements suspended in a net of continuous tension. This structure shows potential for compliance, impact tolerance and mechanical robustness. However, its nonlinear coupled dynamics and often complex geometry require advanced control strategies. An actuator reference planning strategy to bring tensegrity robots closer to controlled full body movements was proposed by Guido Tournois \cite{GuidoThesis} in 2017. This strategy, called the Full Body Reference Planner (FBRP), finds a sequence of equilibrium configurations for a tensegrity structure, predominantly to follow a trajectory in space. However, the method is incapable of incorporating inequality constraints while obtaining said equilibrium configurations. This is a problem when dealing with certain restrictions, e.g., actuator limitations and stability of the structure. <br/><br/>In this thesis we implemented a robust way to account for inequality constraints while utilizing the FBRP. That was done by means of optimization, i.e., an implementation of a Sequential Quadratic Programming method to ensure inequality constraints were respected for each configuration. The approach was validated in scenarios related to practical applications where inequality constraints were enforced. The results showed advancements towards practical feasibility. Furthermore, the robustness, efficiency and accuracy of the method were validated. The extended implementation depicted robustness to parameter variations and good results in terms of accuracy. However, given the iterative nature of the method, it was more computationally expensive than its precursor."Tensegrity; Robotics; OptimizationMechanical Engineering)uuid:c6924adf9f7b48508099056eb06bc11aDhttp://resolver.tudelft.nl/uuid:c6924adf9f7b48508099056eb06bc11aBOptimization Strategies of CO2 Injection for Sequestration and EOR7Chen, Yuan (TU Delft Civil Engineering and Geosciences)TVoskov, Denis (mentor); Delft University of Technology (degree granting institution)The curr< ent situation with green gas emission requires the development of low carbon energy solutions. However, a significant part of the modern energy industry still relies on fossil fuels. To combine these two contradictory targets, we investigate a strategy based on a combination of CO2 sequestration with Enhanced Oil Recovery (EOR) in the hydrocarbon reservoirs. In such technology, the development of miscibility is the most attractive strategy from both technological and economic aspects. <br/><br/>Modeling of this process involves solving complex nonlinear problem describing compositional flow and transport in highly heterogeneous porous media. An accurate capture of the miscibility development usually requires an extensive number of components to be present in the compositional problem which makes simulation runtime prohibitive for optimization. Here, we apply a multiscale reconstructing of compositional transport to the optimization of CO2 injection. In this approach, a restriction operator, based on the parametrization of injection and production tielines, is constructed following the fractional flow theory. This operator is tabulated as a function of pressure and pseudocomposition which then is used in the OperatorBased Linearization (OBL) framework for simulation. As a result, a pseudo twocomponent solution of the multidimensional problem will match the position of trailing and leading shocks of the original problem which helps to accurately predict phase distribution. Applying a simple prolongation operator, based on interpolation between injection and production compositions, a compositional solution with the correct reconstruction of twophase region can be obtained.<br/><br/>The reconstructed multicomponent solution can be used then as an effective proxymodel mimicking the behavior of the original multicomponent system. In this study, we use this proxymodel in the optimization procedure which helps to improve the performance of the simulation with the increasing intricacy of the reservoir model. In this work, the nonlinear constrained optimization function is applied to find optimal pressure values for CO2 injection process to acquire a highest Net Present Values (NPV). Starting with the limited number of the controls, the global extrema of the objective function can be determined from both full physics model and proxy model. The result shows that both models can converge to the similar extrema given various initial guesses. It is noting that for multicomponent system, the proxy model can also have salient accuracy to predict the NPV extrema compared with the original full model. <br/><br/>As a result, the reconstructed proxy model can be significantly cheaper than a full conventional compositional model. An additional benefit of the proposed methodology is based on the fact that important technological features of CO2 injection process can be captured with lower degrees of freedom which accelerates the optimization process.fCO2 injection; Optimization; Compositional modeling; Enhanced Oil Recovery; Carbon capture and storageApplied Earth Sciences)uuid:538fd43d5ea841ca89927f40f8abc450Dhttp://resolver.tudelft.nl/uuid:538fd43d5ea841ca89927f40f8abc450<Optimum design of freeformenabled space optical instruments;Menduia Fernndez, Alvaro (TU Delft Aerospace Engineering)SKuiper, Hans (mentor); Delft University of Technology (degree granting institution)?Nowadays space missions rely heavily on optical payloads to carry out a wide range of tasks including earth observation, weather monitoring, astrophysics research and communications. Until recently, the design of such systems was done according to design principles such as rotational symmetry as it simplifies the theory, reduces technological risks, manufacturing costs and cuts down assembly, integration and testing efforts.<br/><br/>Nevertheless, the evergrowing need to develop more compact and lightweight payloads with enhanced performance often forces designers to adopt tilted offaxis configurations which break the rotational symmetry of the system. In such c< ircumstances, rotationally symmetric optical surfaces cannot compensate the<br/>tiltinduced optical aberrations and thus no longer provide optimum performance over the complete field of view. The solution to this problem is to abandon the conventional approach and adopt a new design paradigmbased on the use of nonrotationally symmetric surfaces, also known as freeform optics.<br/><br/>Most of the stateoftheart methods used for the design of conventional optical payloads are not entirely suitable for freeform optics because: first of all, they were developed on the basis of rotational symmetry (which is no longer applicable) and secondly, because they do not cope well with the large amount of additional degrees of freedom that freeform optics usually entails. Therefore, this Master Thesis was devoted to the development of a novel methodology for the design and optimization of payloads based on freeform optics.:Freeform; optical instruments; Space Systems; OptimizationSpace Exploration)uuid:b67eae7d4d5a47d2a320cf0e14944852Dhttp://resolver.tudelft.nl/uuid:b67eae7d4d5a47d2a320cf0e14944852yAutomatic Selection of an Optimal Power Plant Configuration: Using Client Preferences and time based Operational ProfilesIvan Dijk, Jitte (TU Delft Mechanical, Maritime and Materials Engineering)de Vos, Peter (mentor); Boogaart, R (mentor); Visser, Klaas (graduation committee); Atasoy, Bilge (graduation committee); Delft University of Technology (degree granting institution)TVessels are becoming more and more optimized towards their mission and an increasing amount of systems are considered during this optimization. Most of these considerations occur at an early design stage, where data is scarce and their impact is large. <br/>This research aims to aid designers in considering this growing forest of possible system using concept exploration. During this exploration both the mission of the vessel and the wishes of the owner of said vessel (or client) are considered. <br/>To do so a methodology is developed, which compares a large set of possible power plant configurations based on fuel consumption, emissions, system mass and volume.<br/>To obtain data for these four characteristics the design of each power plant is required. To obtain this the systems inside a power plant have to be designed individually. During the design of these systems several design choices are encountered. <br/>These choices are solved using the preferences of the client, which are defined as multi criteria weight factors related to the aforementioned characteristics. And as such create individual systems and from that a power plant configuration, which is optimal according to the client. <br/>The entire selection process is demonstrated using three separate cases; a bulk cargo carrier, a harbor tug and a Trailing Hopper Suction Dredger.<brJOptimization; Multi Criteria Analysis; Concept exploration; Design choices
20230615)uuid:58e639f7ea6a4e89baf80bfa64c9a7f4Dhttp://resolver.tudelft.nl/uuid:58e639f7ea6a4e89baf80bfa64c9a7f4UThe feasibility of the energyonly market in a highly renewable European power systemAHollebrandse, Julian (TU Delft Technology, Policy and Management)Weijnen, Margot (mentor); de Vries, Laurens (mentor); Auping, Willem (mentor); van der Veen, W. (mentor); Delft University of Technology (degree granting institution)The Netherlands policy on electricity steers towards a carbon neutral power generation in 2050. Among<br/>the studies that focus on long term power system research, none use the starting point of a 100%<br/>renewable power system in 2050. It is unknown if such a system is technically feasible given the<br/>variability of wind and solar. The generation patterns of intermittent renewables differ continuously,<br/>causing the remaining power sources in the system to mitigate the variable IRES generation.<br/>Additionally, it is also unknown to what extent the energy only market can incentivise investment. By<br/>using the PLEXOS electricity market model, the dispatch of CO2 low generation portfolios and the<br/>market pr< ices and net profits of generators can be determined in various scenarios. These scenarios vary<br/>on demand changes, capacity factors for wind and solar and primary fuel prices. The results show a great<br/>sensitivity of an energy only market under the combination of a demand growth and a meteorological<br/>year of low capacity factors, resulting in blackouts and high wholesale market prices. Storage and<br/>intermittent renewables, however, do not benefit enough from a scenario with a high annual market price.yelectricity market; energy only; Optimization; renewable energy sources; investment; merit order; intermittency; dispatch*Complex Systems Engineering and Management)uuid:a9e6065b1560459f8bde992ddf4ede39Dhttp://resolver.tudelft.nl/uuid:a9e6065b1560459f8bde992ddf4ede390MAXUS: Synergizing water, food and energy policy;Burger, Rogier (TU Delft Civil Engineering and Geosciences)Rutten, Martine (graduation committee); van de Giesen, Nick (graduation committee); Abraham, Edo (graduation committee); van Veelen, Peter (graduation committee); Delft University of Technology (degree granting institution)gPopulation growth, meatfocused diets and emerging industries are increasing stress on water, energy and food (hereafter: WEF) supply around the globe. As stress on the resources rises, the interdependencies between the sectors become more apparent and often lead to unforeseen chain reactions. An example of which is: a drought leads to disappointing hydropower generation that leave groundwater pumps inoperable, which in turn lead to disappointing harvests.<br/><br/>Because of these chain effects, synergizing water, food, and energy policy is no easy task. Unconnected institutional entities further complicate this. Farmers are promoted to grow water intensive rice while at a nearby city faces a drinking water crisis.<br/>Sustainability regulations encourage the replacement of food crops by biofuel crops while elsewhere forests are cut to make place for food production.<br/><br/>The Nexus , that is an acronym for the interrelated WEF system, and the need to obtain integrated policies was globally advocated. Several analytical Nexus models where developed but left a lot to be desired according to critics. The main critique of these models so far has been that many cannot serve as a<br/>decision support tool because they lack the ability to investigate specific governance actions or the implementation of technical interventions. These models generally have intensive data requirements and are not scalable and flexible enough to perform for many Nexus studies within a single model framework.<br/><br/>In this study an optimization model framework is proposed, titled MAXUS , specifically designed to address the shortfalls of current models. It was built to customize a model for a specific Nexus study. To test the methodology of MAXUS it was applied to a case study for Ghana and Burkina Faso. Allocation of water and land resources for the final supply of WEF was optimized over space and time, for the objective and constraints given.<br/>It demonstrates how sectors could respond in harmony to changes occurring in one of the sectors. The model shows nontrivial, multisectoral, spatial and temporal tradeoffs for operational management and infrastructural planning of WEF sectors. For example, changing locations of proposed irrigation capacity because of an increase in electricity demand.<br/><br/>Moreover, the case study demonstrated how it can support in decision making on crossborder cooperation by means of separate optimization for Ghana and Burkina Faso. It showed how sharing resources and open trade would lead to a large reduction in infrastructural requirements. In case of noncooperation, it<br/>would be beneficial for Burkina Faso to expand reservoir storage capacity and irrigation capacity on a much larger basis than in case of cooperation. Mainly because it would not benefit from hydropower production generated in Ghana in case of noncooperation. <br/>It also showed how thermal power production would have to fill up the ga< p in electricity supply and how food imports would need to compensate for the loss in food production.<br/><br/>MAXUS is built to serve in a wide range of nexus studies. Hence, it is scalable in time and space, has an adaptable data structure and allows customization of objectives, balances, constraints dimensions and decision variables.<br/><br/>Using Maxus to customize a model for a specific cases study requires broad expertise. It requires engineers to setup the equations for the model, experts in different sectors that understand the interactions, and local experts that now the<br/>conditions that apply in the case study area. Furthermore, it demands policymakers to think through their ambitions, their preferences, their range of possible interventions and governance actions, and the influences of their decisions on each of the WEF sectors. Exploring strategies for responding to developments in the WEF sectors requires cooperation of all sectors.ONexus; Optimization; Water management; Integration; Policy; Water; Food; Energy7.9527706, 1.0307118)uuid:ae23df6646bd4946ba9c1364996fde4bDhttp://resolver.tudelft.nl/uuid:ae23df6646bd4946ba9c1364996fde4b_Modeling stochastic behavior of state variables within the optimization of departure proceduresSmits, Tommy (TU Delft Aerospace Engineering)pHartjes, Sander (mentor); Mitici, Mihaela (mentor); Delft University of Technology (degree granting institution)4Optimization; Stochastic Analysis; Noise; TrajectoryAerospace Engineering)uuid:3d712eb25d6f45268f2cbecc2eeaa13fDhttp://resolver.tudelft.nl/uuid:3d712eb25d6f45268f2cbecc2eeaa13fGOptimal Design of Wind Turbine Blades with BendTwist Couplings Effects/Roeleven, Niek (TU Delft Aerospace Engineering)Holierhoek, Jessica (mentor); Langelaar, Matthijs (mentor); Stolpe, Mathias (mentor); Alderliesten, Rene (mentor); Lund, E. (mentor); Delft University of Technology (degree granting institution)This thesis concerns the optimal design of wind turbine blades for bend twist coupling effects. Bendtwist coupling effects can prove to be an effective way of reducing fatigue damage and root bending moments for a wind turbine blade. However, how can bendtwist coupling be incorporated into a wind turbine blade considering the many structural requirements? This thesis investigates the optimal design for bendtwist coupling in wind turbine blades. Gradient based<br/>numerical optimization algorithms are applied to obtain composite material layups favorable for bend twist coupling. The thesis work builds a framework that allows for the selection of fiber angles at several sections along the blade, while considering structural constraints. This involves the use of a cross section analysis software and a finite element beam formulation for the evaluation of responses, while an aeroelastic software is used for load mitigation validation.<br/>The structural constraints include tip twists and deflection, eigenfrequencies and material failure criteria of which the analytic sensitivities are calculated, implemented in the framework and validated against central finite difference gradient values. A simple cantilever sandwich beam is considered first as a trial for the framework and four problem formulations are defined to test both the framework and the considered optimization algorithms. A full scale wind turbine blade for a 10MW wind turbine is then subject for the optimization. A converged design for the full scale blade is obtained and subject to aeroelastic simulations in a steady state and turbulent inflow. From this, the load mitigation effects are observed and discussed. With the work done in this thesis project, a step closer to an effective framework for aeroelastic tailoring is achieved.5Optimization; Wind Turbine Blade; Bendtwist coupling)uuid:c23bdb6d72db42f68241c44b475521f1Dhttp://resolver.tudelft.nl/uuid:c23bdb6d72db42f68241c44b475521f1BEstimating the Vulnerability of Ship Distributed System Topologiesvvan Leeuwen, Simon (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Marine and Transport Techno< logy)Duchateau, Etienne (mentor); de Vos, Peter (mentor); Visser, Klaas (graduation committee); Negenborn, Rudy (graduation committee); Delft University of Technology (degree granting institution)
All ships rely on many on board systems to reliably complete their missions (e.g. transporting containers, laying pipes or protecting the seas). However, for many of these ships, the detailed design of these systems is put off until the later stages of the design. For most ships, this will create no problems. However, for those complex vessels for which the system design dictates and drives the size and performance of the vessel, changes in a later design stage can be difficult and costly. For these complex ships (e.g. naval ships, pipelaying vessels, etc), more focus should be spent on system design earlier in the ship design process.<br/><br/>A important aspect of the system design is determining and designing the system topology (e.g. how components within a distributed system are connected). To achieve this, the ongoing PhD research of de Vos [2017] created the "Ship Distributed Systems Automatic Topology Generator" (SDSATG) tool. This tool combines the fundamentals of network theory and marine engineering to automatically generate a multitude of ship distributed system topology concepts. A genetic algorithm is then used to optimize these topologies with respect to system performance and cost. The performance objective can be related to the vulnerability (or robustness) of the found system topologies.<br/><br/>This MSc thesis focusses on the investigation and development of a vulnerability assessment method, based on network theory, that can be used to assess the generated ship distributed system topologies. The assessment method will be used to better define the performance of a system design, based on the vulnerability of system topologies in early stage design.<br/><br/>To create the method proposed within this thesis, it is necessary to evaluate the created SDSATG tool, as well as vulnerability prediction methods and metrics available from literature (and used in other fields). Next, the newly developed vulnerability assessment method is discussed. This discussion is based on the assumptions made to define the method, the steps taken within the method, and the advantages and disadvantages of the method. Finally, the new vulnerability assessment method is verified using hand calculations, and the results of a practical system design testcase are used to show the improvements that can be made with regards to the vulnerability of the system topology design.<br/><br/>To conclude, the developed vulnerability assessment method is seen to greatly increase the capability of the SDSATG tool as an eraly stage design tool.<br`Ship design; System design; Vulnerability; Network Theory; Optimization; Resilience; Reliability&Ship Design, Production and Operations)uuid:38e13f1e84244a00bbad70b9eadf97acDhttp://resolver.tudelft.nl/uuid:38e13f1e84244a00bbad70b9eadf97ac'Orbit Modeling of Galilean Moons Flybys#Nardi, Alex (Aerospace Engineering)RNoomen, Ron (mentor); Delft University of Technology (degree granting institution)dThis research aims at modeling the trajectory followed by the Galileo spacecraft during a variety of flybys about the Galilean moons. The chosen flybys have good Earth s elevation angles and either low or high closestapproach altitudes, so that the comparison of the two can give relevant insight into the accuracy of the corresponding trajectories. By propagating the state of the spacecraft during these flybys, optimizing the nominal initial state of the spacecraft (obtained through the SPICE program) and the spherical harmonics of the moons, and estimating new harmonics, the minimum root mean square error between the resulting trajectory and the ephemerides by Jet Propulsion Laboratory (JPL) is found. The analysis of its components along the Local Orbital axes gives insight into the existing relation between them and the Earth's elevation and azimuth angles. In particular, a lowaltitude flyby implies in general a< higher error, but when two flybys have similar altitudes, then the Earth's elevation plays a relevant role and the flyby with the largest one is more likely to have a larger error too. The root mean square error of the fitted trajectories can vary from 15 cm to 7 m, so always less than the 9 m declared by National Aeronautics and Space Administration (NASA) as the maximum error of the moons ephemerides. Furthermore, the flybys about Ganymede and Callisto show a high error in the alongtrack and crosstrack directions, since the radius of their sphere of influence is quite larger than that of the inner moons' ones, hence there is more time for the perturbations to influence the orbit. A byproduct of this research is the estimate of the Galilean moons' gravity field, in particular the new values for their J2 and C2,2 coefficients led to the conclusion that Io is less hydrostatic, while Europa and Callisto are more hydrostatic than previously thought.Galilean moons; Orbit modeling; Flyby; Hydrostaticity; Gravity field; Spherical harmonics; MATLAB; Integration; Propagation; Optimization; Data fitting)uuid:d034d1e2fd984da28ddcc4f6a71acd99Dhttp://resolver.tudelft.nl/uuid:d034d1e2fd984da28ddcc4f6a71acd992Surrogate Modelling for Airfoil Shape optimization2van Velden, Karel (TU Delft Aerospace Engineering)QElham, Ali (mentor); Delft University of Technology (degree granting institution)Airfoil optimizations require costly highfidelity solvers. By replacing the expensive solver with a cheap to evaluate surrogate model genetic optimization becomes feasible and a global optimum can be obtained. This thesis investigates applicability of surrogate model based optimization for airfoil optimization and provides recommendations regarding parameterization and options for building surrogate models.<br/><br/>Building a surrogate model is done by firstly evaluating an initial sampling with the highfidelity solver. A surrogate model is fitted through the training data. An iterative process follows that adds sequential sampling points and refits the surrogate model until a stopping criterion is met. The quality of the surrogate model is determined through error analysis. <br/><br/>CST, NURBS and BezierParsec parameterization are investigated for surrogate model building. To build an effective, consistent and efficient surrogate models CST or NURBS without weights is recommended. Sequential sampling plans should be created with LOLAVoronoi or generalized probability of feasibility. For the surrogate model Kriging is recommended with a minimum sample size of 10 times the amount of parameters.<brtAirfoil; Shape; Optimization; Metamodel; Surrogate; Model sequential; Sampling; CST; BezierParsec; NURBS; MSES; LHD!Flight Performance and Propulsion)uuid:6f3aec1e9225431cb1a8be635907eb67Dhttp://resolver.tudelft.nl/uuid:6f3aec1e9225431cb1a8be635907eb67eModeling and analysis of high frequency high voltage multiplier circuit for high voltage power supplyxQian, Weijun (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Electrical Sustainable Energy)Mao, Saijun (mentor); Qin, Zian (graduation committee); Bauer, Pavol (graduation committee); Rueda Torres, Jose (graduation committee); Delft University of Technology (degree granting institution)pHigh frequency high voltage power supply has been widely applied in many industrial applications such as the medical Xray machine and eletrostatic precipitators. As a part of the high frequency high voltage power supply, the electrical performances of the voltage multiplier circuit will influence the behaviors of applications like the Xray machine such as the imaging quality. The electrical performances include the output voltage drop and voltage ripple, rise time and decay time of output voltage and power losses. In order to get high imaging quality of the Xray machine and reduce damage to patients, the multiplier circuit is required to be designed with low output voltage drop and voltage ripple as well as fast respond time. <br/><br/>This thesis concentrates on the investi< gation of the electrical performances of the Halfwave series CockcroftWalton(HWCW) voltage multiplier circuit. The operations in startup process and steady state are explained in details and methods to evaluate the electrical performances are introduced. Significant parameters of the multiplier circuit that play a role in determining the electrical performances are investigated. Analysis of impact of the parasitic components on the electrical performances are carried out together with simulations. An analytical power loss model is developed in the thesis by deriviations of currents in the HWCW voltage multiplier circuit. At last, optimization of capacitance distributions are discussed and compared to provide methods when selecting the capacitance values in the circuit. The analyses in the thesis are verified by the simulation results in LTspice.2Voltage multiplier; HWCW; Power loss; Optimization)uuid:cc39b153a7994cceaabb80bf35af44c2Dhttp://resolver.tudelft.nl/uuid:cc39b153a7994cceaabb80bf35af44c26Variational Multiple Shooting: Theory and Applications0Greco, Cristian (TU Delft Aerospace Engineering)Noomen, Ron (mentor); Vasile, Massimiliano (mentor); Schrama, Ernst (graduation committee); Hartjes, Sander (graduation committee); Delft University of Technology (degree granting institution)Electric lowthrust propulsion has nowadays found wide application in space dynamics as it entails considerable savings in spacecraft propellant mass, thanks to the very high specific impulse that this kind of engine is able to generate. However, continuous thrust opens new extensive sets of feasible trajectories, and optimization algorithms are needed to mine the admissible search space and find optimal transfers. Practical methods to solve complex optimal control problems as lowthrust trajectory optimization typically involve high computational times, the major bottleneck of these techniques. The objective of this work is the development of a novel multipleshooting optimization tool, employing a variational approach for quick derivative computation, and the assessment of its performance against a variety of test cases. Indeed, after a first analysis of practical optimization methods, the derivative estimation by finitedifference approximations has been found as the major contributor to the computational burden, and the propagation of the variational dynamics has been selected as an accurate approach to speedup their computation. The theory of variational dynamics for multipleshooting application has been analyzed in detail, and further developed for what concerns the secondorder equations. After practical considerations on the method implementation (scaling procedures, sparsity patterns, et cetera) and its interface with WORHP, the selected nonlinear programming solver, the tool has been applied to a broad range of test cases, spanning from elementary problems to practical applications. For what concerns the latter cases, two complex problems were analyzed and optimized: a CubeSat rendezvous departing from EarthMoon L2 and arriving at asteroid 2000SG334, resulting in a propellant mass convenient trajectory suitable for asteroid reconnaissance well before a proposed NASA manned mission in 2069; The Kessler Run, i.e. the 9th Global Trajectory Optimization Competition (GTOC9), in which the developed tool, employed as last step of the optimization cascade of Strathclyde++ team, managed to make the solution constraintfeasible, valid and further massoptimal.OLowthrust trajectories; Optimization; Multiple shooting; Variational equations)uuid:7879208441584ba591e0c88e6033a216Dhttp://resolver.tudelft.nl/uuid:7879208441584ba591e0c88e6033a216Automation of the conceptual design stage for material handling systems: Development of an automated design model for the conceptual design phase of baggage handling systemsIVijlbrief, Mark (TU Delft Mechanical, Maritime and Materials Engineering)Beelaerts van Blokland, Wouter (mentor); Schott, Dingena (graduation committee); Delft University of Technology (degree granting institution)< oBaggage handling systems (BHS), are complex and integrate systems hidden within an airport. In this thesis, a model is developed which can generate BHS concept design automatically. This means that the model defines the type and amount of equipment needed and is able to develop a 3D model that can be used for volume reservation. The model takes into account the stakeholders desires, operational preferences and the terminal shape and size and can therefore define the optimal solution for the stakeholders. By developing this model, a generalized design framework for material handling systems in general can be created,iBaggage Handling System; Automation of design; Generative Design; Optimization; Material handling systems
20240307#Transport Engineering and Logistics)uuid:bfe097c160604701b0d664c5e9bbb0f3Dhttp://resolver.tudelft.nl/uuid:bfe097c160604701b0d664c5e9bbb0f3LVega Launchers' Trajectory Optimization Using a Pseudospectral TranscriptionCDi Campli Bayard de Volo, Giuseppe (TU Delft Aerospace Engineering)Naeije, Marc (mentor); Roux, Christophe (mentor); Volpi, Matteo (graduation committee); Gabrielli, Alessandro (graduation committee); Delft University of Technology (degree granting institution) Given the current technological limitations, launch missions are extremely costly: rockets are huge, being about 100 times more massive than the payload, and most of the times expendable. The specific cost required to place a spacecraft in LEO is between 10 and 25 k /kg. Consequently, trajectory optimization plays a major role in the design of a launch mission. Universities, companies, agencies and research centers have their own optimization software and continuously try to improve the existing methods. In this context, the main objective of the presented research is to assess advantages and disadvantages of using a pseudospectral approach for launcher trajectory optimization in an industrial context by applying it to the Vega launcher family, namely Vega, Vega C and miniVega E.<br/>The software GPOPS is chosen for the transcription of the continuous optimal control problem and SNOPT for the solution of the associated nonlinear programming problem. The launcher ascent problem is formulated considering all the main constraints: angle of attack, pitch rate, aerodynamic loads and heat flux. In addition, the dynamic model includes thrust profiles and drag curves of the aforementioned vehicles. The control to be optimized is the direction of thrust and it is formulated in two different ways, allowing to select the most suited for the specific application: Cartesian xyz components or pitchyaw angles.<br/>After numerous tests with all kinds of target orbits (Sunsynchronous, polar, equatorial, escape, etc.), it can be concluded that the PS method and the software composed by GPOPS and SNOPT are suited for practical ascent trajectory optimization: first, all the typical constraints are included and fulfilled; second, the problem formulation is general enough to be easily applied to all kinds of configurations (number of stages, types of motors, etc.) and number of phases (single/multiple boosts, coast, etc.); third, problem formulation allows for both payload mass maximization (preliminary studies) and propellant reserve maximization with fixed payload (industrial application); finally, the obtained performances do match perfectly with the performance maps of ELV. In conclusion, this research provides the hosting company with a validated method and a complete tool, ready for a daily utilization."Launcher; Optimization; Trajectory
20220901)uuid:72c525b07b48462f9df15864e4b08de3Dhttp://resolver.tudelft.nl/uuid:72c525b07b48462f9df15864e4b08de3mSmallSignal Stability Modelling and Analysis of Modular Multilevel Converter and its ClosedLoop ControllersRSAKINCI, Ozgur (TU Delft Electrical Engineering, Mathematics and Computer Science)kBauer, Pavol (mentor); Kontos, Minos (mentor); Delft University of Technology (degree granting institution)
As the transformation of the electricity generation progresses, the number of pow< er electronicbased devices connected to the electricity grid is expected to grow. To ensure the stability of the power system in a scenario with high penetration of power electronic converters, the dynamics and stability characteristics of these devices have to be comprehensively studied. This is especially important for bulk converters with capacities in the hundreds of megawatts range such as highvoltage direct current (HVDC) systems.<br/><br/>Smallsignal analysis is extensively used in power system stability studies. The calculation of the system eigenvalues provides insights regarding the dynamic behavior and the nature of the interactions within the system. However, generating a smallsignal representation of power electronic converters is known to be problematic, due to the switched operation of such converters and their complex control structures. Nevertheless, usage of modular multilevel converters (MMC) in HVDC systems is rapidly increasing. This creates the need for extensive studies regarding its dynamic behavior. Due to the MMC unique submodulebased topology, special actions have to be taken if the stability of the converter is to be analyzed.<br/><br/>In this thesis, a smallsignal model of a generic MMC was developed by linearizing all the converter plant and control equations, to analyze the converter stability. Averaging was used to eliminate nonlinearities brought forth by the converter switching actions.<br/><br/>The main contributions of this thesis are the application of classic lineartime invariant theory to a timevarying system by linearizing the full converter model at different points in time, and the application of modal analysis methods to the MMC.<br/><br/>The smallsignal model was verified against a largesignal one, both developed in Matlab/Simulink. Additionally, the system eigenvalues obtained analytically and numerically using Matlab/MuPaD and Simulink Control Design Toolbox, respectively were benchmarked. After obtaining the system statespace matrix, its transfer functions were used to estimate the input impedance of the MMC in the frequency domain. Using the linearized system matrix and an optimization algorithm, namely MVMO, the locations of the converter eigenvalues were optimized by finding a new set of controller gains. Finally, calculating the converter eigenvalues, understanding the effects of these eigenvalues on the system variables and applying advanced tools to tune the controllers of the MMC, and possibly also its main circuit parameters, may aid system designers in controller design and parametrization studies.Smallsignal stability; Modular multilevel converter; Optimization; Smallsignal modelling; Voltagesource converter; Eigenvalue analysis; Modal analysis; Participation factors; MMCElectrical Engineering)uuid:bb7dbde731ae4915b7335742e30b50efDhttp://resolver.tudelft.nl/uuid:bb7dbde731ae4915b7335742e30b50efIntegral Internal Transport Planning and Coordination of Multiple Assets and Resources for turnaround of Aircraft for KLM at Amsterdam Airport SchipholDVos, Peter (TU Delft Mechanical, Maritime and Materials Engineering)Beelaerts van Blokland, Wouter (graduation committee); Negenborn, Rudy (mentor); Delft University of Technology (degree granting institution)xCurrently, KLM's apron processes work in separate silos and there is only minimal integration. The result is that processes tend to work to a local optimum instead of a global optimum. A Transport Logistic Control Tower is a centralized information hub from which integral decisions can be made which could benefit the overall turnaround process at Schiphol airport. There is little knowledge about the effects on the system performance when input parameters are varied. A uniform way of scheduling tasks and making integral decisions is missing. The goal of this research is to determine the integral control scenario that should be used by a Transport Logistic Control Tower in order to optimize the overall performance of KLM's apron processes. The main research question is formulated as: Which integral co< ntrol scenario should be used by a Transport Logistic Control Tower in order to optimize the overall performance of KLM's apron processes? In the current state analysis the Delft Systems Approach is used in combination with swim lane analyses to analyze the current system. Relations between the different processes, how they interact and the decisions that are made to plan their tasks are analyzed to get a better understanding of the apron processes. To asses the performance of the apron processes, the On Time Start (OTS), On Time Performance (OTP) and the On Time Delivery (OTD) are introduced. The number of flight delays and the corresponding costs are used to combine the process KPIs in order to make a final decision on the best control scenario for KLM. A thorough data analysis is performed by combining the multiple data sources of KLM to get insight in the performance of the system. This showed us that the data quality is not sufficient, and data does not correspond well with the real system. It is therefore recommended to use historical data in combination with the current calculation methods to improve the way the task duration is calculated. This also means that the way data is stored in the HHT should be automated to reduce the amount of measurement errors. A simulation model is designed to assess different control scenarios and the effect they have on the KPIs. The one day run shows significant improvements when dynamic flight priorities are used in combination with a departure oriented control scenario that schedules based on departure time. Compared to the same control scenarios using static priorities, an average improvement of 31,1\% or 1,7 million euros annually could be seen in the delay costs. This control scenario is therefore recommended and shows the necessity to use a uniform way of planning and the potential of the proposed dynamic flight priorities. The recommendations are done based on the analyses and the found results, and used to advice KLM on their future control scenario for their apron processes and improvement levers from which they could benefit in the future.Optimization; KLM; Apron
20220327)uuid:66aaf06787634c3991e3ef01fb25948dDhttp://resolver.tudelft.nl/uuid:66aaf06787634c3991e3ef01fb25948d8Optimization of Composite Structures for Crashworthiness(Li, Ang (TU Delft Aerospace Engineering)VBisagni, Chiara (mentor); Delft University of Technology (degree granting institution)Structural optimization for crashworthiness in composite structures has become an important topic of research in aerospace, attributable to its proven benefits to the occupants' safety in an aircraft during a crash event. This thesis provides a comprehensive investigation on LSDYNA modeling and design optimization of composite square tubes for crashworthiness. The objective of this thesis is to construct a optimization framework specialized for the preliminary design of composite structures with a balance of performance and efficiency. Three primary components are involved: firstly, couponlevel simulations in LSDYNA are performed to characterize the material properties for carbon/epoxy composite material IM7/8552; secondly, a square tube is modeled by a singlelayer approach in three mainstream material cards (MAT54, MAT58 and MAT262) with calibrated parameters for crush simulation in LSDYNA. Meanwhile, detailed sensitivity analysis of influential parameters in different material models is also be performed to have a more comprehensive understanding of the complex failure mechanism. The simulation results indicate good correlation to the experiments in terms of energy absorption and maximum peak load, with high computational efficiency and lowcost calibration. Lastly, a twostage singleobjective optimization is performed, which incorporates the fiber orientation for each layer as design variables and design/manufacturing rules as constraints. Two surrogate models are created to formulate the mapping between input design variables and output crashworthiness metrics, including Deep Neural Networks (DNN) and Gradie< nt Boosting Regression Trees (GBRT) ensemble. Followed by MixDiscrete Particle Swarm Optimization (MDPSO) algorithm, the optimal set of design variables are obtained for each surrogate model. The firststage optimization results demonstrate significant improvement in the crashworthiness performance compared with the baseline value, while the secondstage optimization results indicate an excellent transferability of the proposed optimization framework. This applicable, transferable, and datadriven optimization framework can be used in the aerospace industry regarding the crashworthy design and optimization of composite structures.8Optimization; Crashworthiness; Surrogate models; LSDYNA51.9900, 4.3754)uuid:588a8ded1f5942ca9eb594f2e1e996b8Dhttp://resolver.tudelft.nl/uuid:588a8ded1f5942ca9eb594f2e1e996b8Expanding the Abalities of Impact Testing for Rubber Characterization: Using an inverse approach for extracting material propertiesIde Roos, Wouter (TU Delft Mechanical, Maritime and Materials Engineering){de Klerk, Dennis (mentor); Pasma, Eric (graduation committee); Delft University of Technology (degree granting institution) Vibrations are experienced as unwanted dynamic motions and the industry takes a lot of care to reduce these motions to improve comfort and to increase vehicle component durability. The properties of rubber, being generally soft while showing large inherent damping, makes the material a perfect choice to use in vibration isolators. To minimize vibrations the transfer path from vibrational source to the driver almost always includes a highly damped material like rubber. With the current trend of electrification in the automotive industry, higher order vibrations are becoming more pronounced. Although the produced sound is not of a high octave with respect to traditional combustion engines, they are experienced as unpleasant. Considering a large frequency range, the dynamic characterization of a rubber object is proven to be a difficult topic. In the search of improved rubber characterization new methods arise. The goal of this research is to provide a new approach for extracting material properties from the experimental impact method. Literature research is done on the rubber specific behaviour and how this affects the modelling approach. In the research itself both numerical and experimental characterization is used. Through impact testing the receptance of the rubber object is measured and by applying inverse substructuring the dynamic stiffness is obtained. The Finite Element model is used to directly obtain the dynamic stiffness. The numerical and experimental results are made comparable by applying the virtual point transformation method. With the finite element material model being the optimization variable, the numerical dynamic stiffness is optimized to match the experimental results. The rubber properties are obtained from the optimized material model. The presented inverse approach makes it possible to use the impact method for characterizing rubber material properties. It therefore broadens the abilities of impact testing used for rubber characterization. With the new method the drivingpoint properties can be predicted which do not depend on complex decoupling methods. The obtained material properties are material, but not geometry related. The findings can be used for different geometries of the same rubber. The new approach has the ability to significantly reduce experimental effort because the material properties can be extracted from one single loadcase, whereas the experimental should excite all loadcases to fully characterize the dynamic stiffness.6Rubber; Characterization; Impact testing; Optimization)uuid:aae2621eb81148debec139171151aae4Dhttp://resolver.tudelft.nl/uuid:aae2621eb81148debec139171151aae4On the turbulence modeling and conjugate heat transfer study of the cooling system of the Beam Dump Facility target at CERN: Steady and Transient analysisHPandey, Mayank (TU Delft Mechanical, Maritime and Materials Engineering)Boersma, Bendiks Jan (mentor); Nemati, Hassan (g< raduation committee); Sola, Edmundo Lopez (mentor); Calviani, Marco (graduation committee); Delft University of Technology (degree granting institution)The Beam Dump Facility (BDF) is a new fixed target facility proposed to be installed in the North Area of CERN. Currently in its design phase, the BDF is aimed at the Search for Hidden Particles (SHiP) experiment, which purpose is to investigate the origin of dark matter and other lightly interacting particles. The BDF target/dump sits at the core of the installation, and its aim is to fully absorb the high intensity Super Proton Synchrotron (SPS) beam and produce charmed mesons. The average beam power on target is 305 kw due to which very high thermomechanical loads will be generated on the target/dump configuration which can lead to mechanical failure of the target material. This calls for continuous cooling and heat removal from the target material, which requires the design and optimization of a very complex cooling system to ensure the target reliability during operation.<br/>Keeping in view the high velocities that are obtained in the channels of the cooling system of the target assembly, an extensive 3D turbulence modeling of the full scale cooling system is required in order to predict the target water cooling system behavior. A comparative conjugate heat transfer (CHT) study using a simplified 2D geometry of the BDF target for different mesh size was performed for identical y+ values to validate the pressure drop in the cooling channels and the temperature profile in the target blocks with the analytical calculations. Subsequently, a 3D model of the cooling channel was simulated for different Reynolds number and an extensive study was performed to check the behaviour of the flow in the loglaw layer. In addition to that the friction factor was validated with the analytical results and with the available literature for various Reynolds number. Thereafter a full scale 3D steady and transient model was simulated using the information obtained from the previous studies. All these simulations were carried out in Ansys Fluent. The energy deposition in space on the target blocks was obtained via FLUKA MonteCarlo simulations. The variation of HTC in different channels and the fluidsolid interface temperature is found out to be in accordance with the analytical calculations. The pressure drop and the temperature rise from inlet to outlet of the cooling system is also in confirmation with the design parameters. Transient simulations were performed subsequently in order to study the impact of time and space varying energy deposition in the target blocks on the overall flow. With the given inlet velocity, the boiling at the surface of the target blocks is not expected to happen. As a final part of this study, in order to make the numerical model more robust, a mesh sensitivity analysis was done in order to optimize the mesh size, especially in the boundary layer region.Beam Dump Facility; Target; SHiP; SPS; Turbulence Modelling; Conjugate Heat Transfer; FLUKA; Fluent; CERN; Cooling system; Optimization
201812196Mechanical Engineering  Energy and Process Technology4Beam Dump Facility  SHiP  Physics Beyond Particles)uuid:3378e0cf78314432a383c9c5c8dbce2dDhttp://resolver.tudelft.nl/uuid:3378e0cf78314432a383c9c5c8dbce2d;Coorbital motion and its application to JAXA's MMX mission5Bernal Mencia, Pablo (TU Delft Aerospace Engineering)pNoomen, Ron (mentor); Kawakatsu, Yasuhiro (mentor); Delft University of Technology (degree granting institution){In the framework of JAXA's MMX mission to explore the Martian moon of Phobos, an analysis of the stability of threedimensional quasisatellite orbits in the MarsPhobos circular restricted threebody problem was conducted. For this analysis, notions of coorbital motion, interpreted as the slow motion of the guiding center of the trajectory along the disturbing potential of Phobos, were used. After identifying and analyzing different regions of stability for three quasisatellite orbits at 100, 50 and 30 km from the center of Phobos,< several conclusions were drawn regarding the dynamics of the ballistic escape of the spacecraft, interpreted in terms of coorbital motion. By making use of these insights, a novel methodology to find periodic quasisatellite orbits able to reach high latitudes over the surface of Phobos was derived. This methodology consists of two steps: a multiobjective minimization using coorbital parameters as target functions, to isolate regions with potential periodic orbits; followed by a shooting algorithm to arrive at the final periodic orbit. As a result of this new methodology, two periodic orbits were found at 50 and 30 km from the center of Phobos, able to reach latitudes as high as 54 deg and 32 deg respectively. This result represents an important contribution to both the operations and the scientific return of the Phobos proximity phase within MMX.fAstrodynamics; Orbit Design; Optimization; Multiobjective optimization; Periodic Orbits; Mars; Python
20200101)uuid:13e189889b504cc8bf8e392455acfb44Dhttp://resolver.tudelft.nl/uuid:13e189889b504cc8bf8e392455acfb44Optimal Power Management System of EVs Charging from PV System in a Low Voltage Distribution Network: An integration between EVs, PV system and gridShiddiq Sumitro, Muhammad (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Electrical Sustainable Energy)Bauer, Pavol (mentor); Chandra Mouli, Gautham Ram (mentor); Izadkhast, Seyedmahdi (mentor); Smets, Arno (graduation committee); Ramirez Elizondo, Laura (graduation committee); Delft University of Technology (degree granting institution)QAn increase in carbon emission which mostly caused by the transportation sector and electric power generation has been a hot topic nowadays in most countries in the world. To tackle this problem, the share of renewable energy use has been increased by up to 14% in the Netherlands. Moreover, the number of electric vehicles (EVs) on the road also reaches a total of 120,000 EVs in 2018. However, the high penetration of renewable energy sources (RESs) such as solar & wind power and the EVs charging in the distribution network could result in a severe problem. One of the solutions to avoid this problem is that switching the uncontrolled charging of EVs into a controlled charging or called as smart charging. Further, an integration between the EVs, RESs and the distribution grid could potentially lead to technical and economic benefits. The focus of this thesis is to develop an optimal power management system (PMS) between the EVs, PV system, and the distribution network. The goal of the power management system is to obtain the minimum operational cost while also considering the technical grid constraints, which subsequently could avoid the grid violation. The proposed power management system will be modeled in a mixed integer nonlinear programming (MINLP) optimization problem and executed in General Algebraic Modelling System (GAMS) software. To evaluate the performance of the proposed power management system, a comparison between the withgrid and the nogrid constraints case will be performed through several case studies. This study shows that by implementing the proposed power management system of EVs charging from PV system considering the grid constraints, it could decrease the total operational cost remarkably by 18.16%  214.08% when compared to the uncontrolled charging scheme. Besides, the grid problem caused by the uncontrolled charging process such as exceeding the allowable voltage deviation and the transformer rated power could be prevented. However, in comparison to the smart charging without considering the grid constraints, the operational cost is increased by 1.43%  113.20%.;Smart charging; Electric Vehicles; Optimization; PV systemsSustainable Energy Technology)uuid:d8f71f9a6b594b519fe5b0a4358a0d7fDhttp://resolver.tudelft.nl/uuid:d8f71f9a6b594b519fe5b0a4358a0d7f7Parallel Machine Scheduling: with Partition ConstraintsWSwennenhuis, Cline (TU Delft Electrical Engineering, Mathematics and Computer Science)van Iersel, Leo (men< tor); Janssen, Teun (mentor); Aardal, Karen (graduation committee); Delft University of Technology (degree granting institution)0Scheduling; Partition; constraints; Optimization"Applied Mathematics  Optimization)uuid:64606bb7ce5c42ddaa257c1045ef3f80Dhttp://resolver.tudelft.nl/uuid:64606bb7ce5c42ddaa257c1045ef3f80bLong trem transmission expansion planning for Nigerian deregulated power system: A system approach
Egeruoh, C.C.Andres Ramos, G. (mentor)Nigerian suffers from inadequate electricity supply of which her transmission capacity is a major factor. Having deregulated her power system recently , systems planning( generation and transmission) becomes a chicken and egg problem of which its coordination is a complex sociotechnical issue. In an effort to approach this problem ,TEPES ( a long term transmission expansion model developed in IIT, Upcomillas , Spain), was used in this thesis to design a tactical transmission expansion plan that will be indicative for both optimal power wheeling capacity and better generator placement up to year 2030. This approach considered the expected generation mix in the Nigerian energy policy.]Optimization; sociotechnical problem; transmission expansion planning; benders decompositionCampus only!Technology, Policy and ManagementMulti Actor Systems)uuid:d60603b662ba4be9a54fff371228e12aDhttp://resolver.tudelft.nl/uuid:d60603b662ba4be9a54fff371228e12aYOptimization of the maintenance scheduling and production valuation by green certificatesMeriem Bourada, M.Yannick Perez, Y. (mentor)Summary Part I  Optimization of the maintenance scheduling Resource management in the industry is an important feature to make the most of the production, particularly in maintenance field where the availability of people is a very important condition to minimize losses. An optimized maintenance scheduling is of a crucial importance for wind farms. There is great potential to reduce costs by optimizing maintenance for project profitability. We propose here a methodology that will supply us a schedule for resource allocation for the maintenance of wind turbines. Our goal is to create a decisionmaking tool allowing us to allocate resources annually for maintenance of wind farms of the company (or a number determined by the user), while minimizing production losses by taking advantage of production forecasts. Concerning a shortterm maintenance (weekly), interventions have to be led to optimize the use of wind farms and their efficiency. For that purpose, a graphic interface was built for the benefit of the user based on weekly production forecasts offering an alternative to the former planning of interventions established in terms of efficient combination between production and maintenance. Summary Part II  Green Certificates The valuation of production also passes by green certificates. Indeed, the liberalization of energy market has created various instruments compatible with the introduction of competition, such as markets for green certificates that represent a financial instrument for supporting the development of green electricity as adapted to the constitution of liberalized markets. However, in France the process is voluntary and benefits are still low, but an increase is felt from businesses looking to promote their products through that mean in one hand and the purchasers who show interest in renewable energy. The actual efficiency to create a stable framework and incentive for investment is to be measured by current experiments in the EU member state.HOptimization; Maintenance; Wind turbines; Green certificates; SchedulingfrEngineering and Policy AnalysisEPAEMIN)uuid:581a19ab1bac4fbba678bc8e905804d7Dhttp://resolver.tudelft.nl/uuid:581a19ab1bac4fbba678bc8e905804d7qMass optimization of PocketQube structures and deployers: A method demonstrated on the DelfiPQ PocketQube mission.Spronck, Tano (TU Delft Aerospace Engineering)Bouwmeester, Jasper (mentor); Zandbergen, Barry (graduation committee); Sundaramoorthy, Prem (graduation committee); Fatemi, Javad (graduation commi< ttee); Delft University of Technology (degree granting institution)the aim of this thesis is to regard all aspects concerning mass optimization with respect to small satellite structures and deployers, specifically aimed at the DelfiPQ PocketQube mission.PocketQube; mass; Optimization)uuid:50aaee056a0b4000b685572d175d5dd7Dhttp://resolver.tudelft.nl/uuid:50aaee056a0b4000b685572d175d5dd7LOptimizing electrodialysis processes for concentrating ammonium rich streamsTDeckers, Rob (TU Delft Civil Engineering and Geosciences; TU Delft Water Management)Spanjers, Henri (mentor); van Linden, Niels (mentor); van Lier, Jules (graduation committee); Witkamp, GeertJan (graduation committee); Lindeboom, Ralph (graduation committee); Dewil, Raf (graduation committee); Delft University of Technology (degree granting institution)Ammonia in aqueous environments leads to eutrophication and toxicity of the receiving water body (Metcalf & Eddy, 2003). In order to prevent environmental pollution, ammonia needs to be removed from residual water streams, before the water can be discharged to the aqueous environment (Song et al., 2012). Waste water treatment plants (WWTPs) reduce the ammonianitrogen concentration by applying biological treatment technologies. However, the oxidation of ammonia by bacteria requires high amounts of oxygen and consequently has a high energy consumption. The nitrification and denitrification of wastewater requires 15.83 kWh per kgN (Magri et al., 2013). Moreover, the reject water from the digested sludge contains high amounts of nitrogen, which is fed back to the biological treatment, and can contribute 1520% of the nitrogen load (Fux et al., 2002). Concentrating ammonium fluxes lead to the production of a concentrated ammonium stream, used as potential energy source in the N2kWh project, and a sweet stream.[Electrodialysis; Concentrating; Optimization; ammonia; Brine treatment; Ammonium; Fuel Cell
20181107From Pollutant to Power)uuid:446ce1ef5c164ceba31beab72884b4d4Dhttp://resolver.tudelft.nl/uuid:446ce1ef5c164ceba31beab72884b4d4LFast computation of SEP transfers to Mars using analytic curvefit functions0Galletti, Elena (TU Delft Aerospace Engineering)Noomen, Ron (mentor); Visser, Pieter (graduation committee); Cowan, Kevin (graduation committee); Delft University of Technology (degree granting institution)>
The increasing interest in solar electric propulsion techniques, that enable significant propellant mass savings for a wide class of transfers, has brought about a revolution in the approach to trajectory design and optimization, as a result of the complexity and diversity of the problem. To date, solutions of different nature exist, but numerical methods that require significant computational effort and user experience are typically used already in the early stages of mission design, due to the limited availability of reliable mediumtolow fidelity design tools for SEP transfers. This research project proposes a novel method that computes transfer performance parameters for EarthMars massoptimal SEP transfers, by means of empirically derived analytic relations. The method is intended for applications such as concurrent engineering and earlyphase concept development, which require the fast characterization of a broad design space. Besides accommodating a wide range of currently available systems, the method successfully deals with modelling the effect of nonzero infinity velocity at departure and/or arrival.<br/><br/>The methodology that has been applied consists of a first phase of generation and characterization of the transfers, and of the subsequent selection of the model variables, model functions and architecture. Regarding the generation of the transfers, it is assumed that the transfers are coplanar and that the initial and target orbits are circular. Hundreds of transfers are optimized in a semiautomatic way and characterized in terms of thrust profile and transfer performance parameters. In the investigated design space, different regimes are identified, but approximately < 90% of the acceleration range of interest falls into the thrustcoastthrust profile for any combination of departure and arrival infinity velocity. For a proper description of the underlying trends in the transfer parameters, three key variables have been identified, namely the average acceleration, the total infinity velocity and the infinity velocity at arrival (expressed as a function of the total infinity velocity). By means of curvefitting, analytic relations are derived that successfully describe those trends, limited to the thrustcoastthrust class of transfers. <br/><br/>The method that is presented computes nearoptimal transfers in terms of V cost, transfer time, transfer angle and departure date. While the first three parameters are the outputs of the mentioned curvefit model above, the departure date is computed by solving analytically the problem of the phasing with Mars, in a subsequent step. The fit functions that are derived model circletocircle planar transfers with an accuracy in the order of 0.1% with respect to the V , 1.5% to the transfer time and 1.2% to the transfer angle, successfully dealing with the dependence on the departure and arrival infinity velocities and generating instant estimates for all relevant transfer parameters. When the model performance is considered in relation to transfers derived in the full ephemeris model, the errors are within 1% for the V , within 15% for the transfer time and within 12% for the transfer angle, which, together with the demonstrated efficiency and simplicity of implementation, make it suitable both for earlystage assessments and for generation of suitable first guesses.OTrajectory Design; Solar Electric Propulsion; Optimization; Mars; Curve Fitting
20190918Astrodynamics & Space Missions)uuid:5b0b6016b8d242d3a9ba6f18e9c63bcaDhttp://resolver.tudelft.nl/uuid:5b0b6016b8d242d3a9ba6f18e9c63bcaThe Systematic Search With Assignment Problem Solving (SSWAPS) Algorithm: a New Solution Method for the Nurse ReRostering Problemvan den Dool, Vera (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Delft Institute of Applied Mathematics)Zvan Essen, Theresia (mentor); Delft University of Technology (degree granting institution)Hospitals need to be staffed around the clock to ensure satisfactory care for all patients. To this end, a set of shifts is defined which are assigned to nurses in a work schedule. The nurse rerostering problem (NRRP) occurs when due to unforeseen circumstances some shifts become unassigned: one or more nurses are no longer able to work some of the shifts assigned to them. Changes must then be made to the existing schedule to ensure that the shift occupancy demand is once again satisfied. The number of changes should be minimized to disrupt nurses personal plans as little as possible. In this paper, a heuristic solution method (SSWAPS) is proposed for the NRRP. The heuristic systematically searches the solution space and estimates the solution quality using the assignment problem.Optimization
20170930)uuid:611d466d3d264b3299e6e98a805435e8Dhttp://resolver.tudelft.nl/uuid:611d466d3d264b3299e6e98a805435e8Ore Pass Design And Placement<Koivisto, Mikko (TU Delft Civil Engineering and Geosciences)Buxton, Mike (mentor); Rinne, M. (mentor); Lottermoser, B.G. (mentor); Toodu, Elen (mentor); Delft University of Technology (degree granting institution)xA current trend in mining is towards lower grade deposits when easy to access mineral deposits are depleted in the world. Therefore, to make the current deposits accessible, a cost of mining need to be economic. Costs can be decreased in equipment and mine planning to maximize the availability and utilization of equipment for example by decreasing travelling and hauling distances in daily production. Ore passes are one of the possibilities to decrease hauling costs when mine is evolving the greater depths.<br/><br/>Ore pass design and placement is conducted using modified Bieniawski s design strategy for rock structures. Thesis studies in detail all paramete< rs affecting to the ore pass design and placement. With the strategy and data available from mine and literature, a case study for Kittil Mine had been carried out.<br/><br/>Data collection from literature review and data acquired from the mine site was made and applied in the Bieniawski s design strategy. Production related factors were weighted more in the design process in order of mining to assess the ore pass design and placement. The location and design of the ore pass were optimized using a stateoftheart technology and best practices in combination with mathematically optimized location.<br/><br/>Via strategy it was possible to identify the project risks, technical viability and feasibility of the ore pass project. Case study shows that via applying the design strategy it is possible to show that mine can benefit from using ore passes in its operations and that ore passes are financially viable solution if ore passes can be operated the way presented in this study.JOre pass; Mining; Mine Planning; Optimization; Underground Mining; Hauling)uuid:61a8031c9278494d88920fc4c8252b84Dhttp://resolver.tudelft.nl/uuid:61a8031c9278494d88920fc4c8252b84rDeveloping a parametric script for the preliminary design and optimization of concrete balanced cantilever bridges:Scheele, Bram (TU Delft Civil Engineering and Geosciences)Nijsse, Rob (graduation committee); Jonkers, Henk (graduation committee); van der Linden, Lennert (graduation committee); van der Sanden, Janine (mentor); Delft University of Technology (degree granting institution)In this thesis, research is performed into parametric design of balanced cantilever bridges within the total road design. The main research question is: What can be achieved by using a parametric model in the preliminary design phase, with the possibility to optimize and compare concrete bridges on material usage, costs and environmental impact? And to answer this question, the cantilever bridge and all other elements in the road design are investigated and a parametric model is developed. In this model the total vertical alignment of the road design is implemented along with approach bridges, substructure and foundations. Furthermore, the results of the designs are visualized in a 3D environment compared with an overview of the material costs and environmental impact. In the "Discussion", the following is concluded: Although the parametric model for the preliminary design of concrete cantilever bridges took a long time to develop, this time investment can be regained by using the model several times. Because by using the parametric model for the preliminary design of concrete bridges, the design time per procurement is decreased significantly. Furthermore, the design process of the balanced cantilever bridge has become more flexible since changes are automatically implemented for the whole model. Also the chance for (human) errors is decreased, for instance the risk of data loss is minimized due to the model s single environment. Since only a limited amount of calculations is performed by the model, the results of the analysis must be used with caution. Another advantage is the instant results as produced by the model. These results consist of a visualization of the design and an overview of the material amounts, with their costs and environmental impact. Furthermore, the model is also able to optimize these results in a fast way due to a smart evolutionary solver. This all together makes the parametric model a useful tool with high potential, but it should be further<br/>developed and improved.\Parametric design; Concrete; Bridge; cantilever; Materials; Costs; Optimization; Grasshopper
20200219Civil Engineering)uuid:2c79c37252d8448993ae1e2f7c673469Dhttp://resolver.tudelft.nl/uuid:2c79c37252d8448993ae1e2f7c673469^LaserEnhanced Solar Sailing: Modeling and Trajectory Optimization for Interplanetary Missions/Carzana, Livio (TU Delft Aerospace Engineering)Noomen, Ron (mentor); Dachwald, Bernd (graduation committee); Stam, Daphne (graduation committee); Sundaramoorthy, Prem (graduati< on committee); Delft University of Technology (degree granting institution)When sunlight illuminates a body, a tiny pressure is exerted upon its surface due to the photons impacting on it. Such a principle forms the basis of solar sailing, in which the solar radiation pressure is used to accelerate highly reflective lightweight structures called solar sails. Similarly, a laserenhanced solar sail is a solar sail in which also an external laser beam pointed towards the sail is exploited to generate thrust. In this way an additional laser radiation pressure is exerted onto the sail, hence conferring higher propulsive and steering capabilities, and leading to an increased maneuverability of the sailcraft.<br/><br/>The main purpose of this thesis project is to provide a model of the laserenhanced solar sail dynamics and to establish the advantages of laserenhanced solar sailing as compared to "traditional" solar sailing. The analysis has been pursued focusing on interplanetary missions and considering ideal sails, i.e. sails able to perfectly reflect the impinging radiation.<br/>Normally, for lowthrust interplanetary missions the propellant consumption and time of flight required for the transfers to take place play a crucial role. However, since solar sails do not exploit any propellant, the traditional and laserenhanced sailcraft performances have been compared by analyzing their flighttime optimal trajectories, focusing on three different mission scenarios: a Mercury orbit rendezvous, Mars orbit rendezvous and Neptune flyby.<br/>These trajectories have been computed by taking advantage of an evolutionary neurocontrol optimizer, in which newly added functionalities have been implemented with the purpose of optimizing laserenhanced solar sail trajectories.<br/>The trajectory analysis results have shown that, if laserenhanced sailcraft are used instead of traditional sailcraft, flight time gains in the order of 811% can be achieved for the missions to Mercury and Mars orbits, while a smaller 2.5% gain is achieved for the flyby mission to Neptune.hLaserenhanced; Solar Sail; Ideal; Trajectory; Optimization; Machine Learning; Evolutionary Neurocontrol
20180701)uuid:6dfc6eb9676f4a0e99c58462ab70a3cbDhttp://resolver.tudelft.nl/uuid:6dfc6eb9676f4a0e99c58462ab70a3cbOptimal Battery Energy Storage System Sizing and Placement in Low Voltage Distribution Network: Emerging robust electric vehicle charging using photovoltaics\Fakhri Muhammad, Tubagus (TU Delft Electrical Engineering, Mathematics and Computer Science)= Electric Vehicles (EVs) provide a promising sustainable clean mode of transportation with much fewer emissions compared to the traditional Internal Combustion Engine (ICE) vehicles. However, this can only be considered true if the energy sources used for charging EVs are also from clean energy form. Various Renewable Energy Sources (RES) are available to provide the EVs charging demand. Nevertheless, due to intermittency power generation characteristic of RES, there is a high risk of supply and demand power mismatch. A gridconnected topology system is commonly used as a backup for lack or over power supply from RES generation. Even though, this configuration may lead to new problems of under/over voltage issues and maximum transformer/lines loading capacity due to RES fluctuation power production. The storage system is considered necessary to be installed in the system as a buffer of the power mismatch and also reducing the grid stress. Excellent approximation of storage sizing and placement is mandatory to provide satisfactory influence to the overall system technically and economically.<br/><br/>This study proposes optimization strategy for optimal battery sizing and placement with optimal power management system in the EVs charging station using photovoltaics (PV) power in a low voltage distribution network using Mixed Integer Linear Programming (MINLP) algorithm. The lithiumion battery is considered in this project because of its relatively high power and energy density characteristic along with low selfdisch< arging. The radial low voltage CIGRE benchmark is adopted as a model topology to evaluate the performance of the proposed optimization strategy. Dynamic daily electricity price, residential load, and PV power profile of the Netherlands are taken into account in this project. Furthermore, the research project implements several case studies using DICOPT (Discrete and Continuous OPTimizer) solver in GAMS (General Algebraic Modelling System) software with MATLAB interfacing. As the main result, proper EVs charging strategy with optimal battery sizing and placement with the appropriate power management system is successfully maintaining the system away from the grid violation while at the same time reducing the total cost of the system.<br/><br/>Keywords: EVPV charging, optimal battery sizing and placement, MINLP<brElectric Vehicles Charging; Battery Energy Storage System; Storage Sizing and Placement; Photovoltaics Power; Optimization; Mixed Integer Non Linear Programming)uuid:a732d9ba71744b2081aa01f387e5de35Dhttp://resolver.tudelft.nl/uuid:a732d9ba71744b2081aa01f387e5de355Multiinput Embedded Realtime Software Defined RadioJalaludeen, A.EGaydadjiev, G.N. (mentor); Bertels, K. (mentor); Gross, H.G. (mentor)eThis master thesis work is inspired by the practical aspects of the Analytical Constant Modulus Algorithm (ACMA) proposed by AlleJan van der Veen and Arogyaswami Paulraj. The ACMA deals with the beamforming problem associated with constant modulus cochannel signal interference in wireless communication. Cochannel interference occur when multiple signals are transmitted simultaneously at the same frequency from different sources. Beamforming is a technique applied in spatial signal processing to separate out individual signals using an antenna array. ACMA provides an efficient analytical approach to solve beamforming. It is a \textit{blind beamforming} algorithm as it does not require any knowledge of the signals and the channels. The scope of this thesis work is to implement a low cost, low power, embedded Multiinput receiver application. The Multiinput receiver system shall handle partially and fully overlapping constant modulus signals from distinct sources. The signals are modulated using a generic modulation scheme. The Multiinput receiver shall separate out individual signals using the ACMA (blind beamforming) and demodulate each signal. The system shall be software defined such that the beamforming and demodulator are implemented in software on an embedded Digital Signal Processor (DSP) platform. The Multiinput receiver system shall be optimized for the DSP platform in order to achieve realtime performance in terms of speed and stay within the power budget of the system. The receiver shall work efficiently in the presence of interferences and noise such as thermal noise generated by the receiver.8Analytical Constant Modulus Algorithm; DSP; Optimization
201408108Electrical Engineering, Mathematics and Computer Science'Microelectronics & Computer EngineeringEmbedded Systems)uuid:c5ddec13e10349b48b84a8ad013c753cDhttp://resolver.tudelft.nl/uuid:c5ddec13e10349b48b84a8ad013c753cIMultifidelity CoKriging Optimization using Hybrid Injected RANS and LES4Fatou Gomez, Javier (TU Delft Aerospace Engineering)oHickel, Stefan (mentor); Dwight, Richard (mentor); Delft University of Technology (degree granting institution)The computation of complex turbulent flows design optimization processes is currently limited by the lack of accuracy of ReynoldsAveraged NavierStokes (RANS) in massively separated flows and the infeasible cost of multiple LargeEddy Simulation (LES) evaluations. A novel method is presented, injecting data from LES or other highfidelity source such as DNS into the RANS equations, forming a Hybrid Injected RANS (HIRANS) model. The aim is to construct a multifidelity design optimization framework that outperforms singlefidelity RANS and LES variants. Two different formulations, injecting a scaled version of the nondimensional anisotropic part of the Reynolds stress tensor and both is< otropic and anisotropic components, are tested in the periodic hill case. A costeffective LES configuration is assessed, and the agreement of the RANS and HIRANS results with respect to the LES reference is investigated. The original geometry is parametrized using a hill width multiplier, computing several LES evaluations. The injection of LES information from the same geometry into HIRANS and the prediction capabilities when using interpolated data from different geometries are tested. A global design optimization process is computed, using singlefidelity RANS, HIRANS and LES Kriging and multifidelity RANSLES and HIRANSLES CoKriging surrogates. The objective function is based on a combination of turbulent mixing and total pressure losses. <br/><br/>The correction of the mean velocity components required of the injection of both isotropic and anisotropic components for the test case. The LES setup analysis yielded similar results to the reference data with one tenth of grid points and forty percent of its averaging period. The locally corrected HIRANS model successfully reduced the L2 norms of the Reynolds stresses with respect to LES to a third part of the original RANS values in the fiftynine LES samples tested, with a modest improvement in the mean velocity components. The nonlocal corrections yielded irregular results for the mean velocity components, with successful corrections of the Reynolds stresses despite the long distances in the parameter space and different flow features of neighbouring LES cases to interpolate from. In the optimization process, the CoKriging LESHIRANS was not able to outperform the CoKriging LESHIRANS and Kriging LES methods. It improved the initial prediction of the underlying function, but the surrogate yielded artificially low predicted errors far away from the LES samples, leading to an overly exploitative method. An error correction formulation combining two HIRANS fidelity levels was simulated using a modified Kriging believer criterion, outperforming the original formulation and achieving similar results as Kriging LES. The computational efficiency improvements for future research of the CoKriging HIRANS are suggested to be linked to an adequate error estimation integration into the surrogate model.RANS; LES; HIRANS; Turbulence; CFD; Interpolation; Optimization; Kriging; CoKriging; Aerodynamics; Periodic hill; Correction; Prediction
20181214)uuid:efd1c23f4ab741dd88e4e9a1683c4cccDhttp://resolver.tudelft.nl/uuid:efd1c23f4ab741dd88e4e9a1683c4cccComputational Design Method Based on Multidisciplinary Design Optimization and Optioneering Techniques for Energy Efficiency and Cost EffectivenessJRodriguez Garcia, Alvaro (TU Delft Architecture and the Built Environment)vTurrin, Michela (mentor); van der Spoel, Willem (mentor); Delft University of Technology (degree granting institution)One of the cuttingedge interests gaining relevance among design offices, notably in the sphere of architecture firms, is computational design. With its inception embedded in the 60s and later reformed through the CAD revolution that boomed in the 80s, the use of the computer as a design tool has evolved to such lengths that it became indispensable to contemporary practices. Nowadays, under the bright prospect of a challenging future and considering the sustainability hypothesis  sophisticated computational approaches and parametric design can genuinely be considered a most promising alternative for exceeding outcomes in the fields of performance and efficiency. <br/><br/>Thus, under this spectrum, the following research focuses on investigating the evolution of the aforementioned tools, through understanding the current user experience and interface in prevalent architectural offices. Consequently, through the use of existing computational tools  simulations, optimization and design exploration techniques  it will put forward a new method and use pattern that would also project the potential and the future of these technologies, specifically applied to a sustainable, energyefficient and cost< effective vantage.Computational Design; Optimization; Performance Based Simulations; Multidisciplinary Design Optimization; Energy Efficiency; Costs effectiveness; Sports Facilities; Workflow; Optioneering)uuid:1c741d018a3f45b7abe9f32456791dc3Dhttp://resolver.tudelft.nl/uuid:1c741d018a3f45b7abe9f32456791dc3=Optimization of the preliminary design of a radial compressorPietersz, J.W.H.2Olivero, M. (mentor); Van Buijtenen, J.P. (mentor)Optimization is the mathematical process through which the set of conditions is obtained that produces the maximum or minimum value of a specified function. To solve different types of optimization problems several algorithms has been developed. They can be globally distinguished as deterministic or stochastic algorithms. Further on deterministic and stochastic methods can be divided in gradient based and derivativefree methods, simulated annealing and evolutionary algorithms respectively. Optimization methods mentioned above are used in many fields including engineering. This thesis is dedicated to the optimization of the design of a radial compressor using an optimization algorithm. Optimization; radial; compressorProcess and EnergyEnergy Technology)uuid:28671d8aa00a4526a6700c1e80ad50ffDhttp://resolver.tudelft.nl/uuid:28671d8aa00a4526a6700c1e80ad50ff6Optimal Placement and Sizing of Distributed Generation
de Luis, R.H.)Stadler, M. (mentor); Heleno, M. (mentor)
In this thesis work we solve the problem of optimal placement and sizing of distributed generation by using an original Fuzzy Adaptive Particle Swarm Optimization algorithm and a Mixed Integer Linear Programming formulation of the problem. The goal of integrating Fuzzy Logic in Particle Swarm Optimization is to be able to overcome some of the classical disadvantages of the algorithm. Particle Swarm Optimization has been traditionally criticized for the complexity to set the acceleration constants of the algorithm and the low exploration capabilities of the algorithm. In this thesis work, it is proposed a new implementation of Particle Swarm Optimization that avoids the complex setting procedure of the acceleration constants of the algorithm, while aiming at improving the exploration capabilities of the algorithm. In this thesis work it is also analyzed a novel Mixed Integer Linear Programming formulation of the problem of optimal sitting of distributed generation in distribution networks implemented in DER CAM, a tool developed at the Lawrence Berkeley National Laboratory. The results of the two models are analyzed and contrasted against each other for a real case study of an islanded microgrid located in the north of the U.S. Results obtained in the case study depict the differences between the two analyzed approaches to solve the problem. It is found that over and under estimations of voltage magnitudes in high and low loading scenarios of distribution networks have the potential to impact investment decisions in distributed generation capacity for the linear formulation of the problem. Also the models analyzed depict the synergies between renewable energy technologies and thermal generators to increase energy savings while maintaining the operation limits of the grid.Distributed Generation; Investment Planning; Microgrids; Distribution Networks; Optimization; Particle Swarm Optimization; Linear Programming)uuid:f2e84a4aaef445b5a9d0e97afd17b04bDhttp://resolver.tudelft.nl/uuid:f2e84a4aaef445b5a9d0e97afd17b04bGA robust optimization approach to synchromodal container transportationRChiscop, Irina (TU Delft Electrical Engineering, Mathematics and Computer Science)De Klerk, Etienne (mentor); Phillipson, F. (graduation committee); Sangers, A (graduation committee); Nane, Tina (graduation committee); Delft University of Technology (degree granting institution)This thesis addresses synchromodal planning at operational level from the perspective of a logistics service provider. The existing infrastructure and the transportation activities are studied and modeled as an optimization problem with simultaneous vehic< le routing and containertomode assignment. A special characteristic of this problem is the uncertain data. In other words, it is assumed that the release times of the containers belong to an uncertainty interval, and no further statistical information is available. The problem is then classified according to an extensive framework previously developed within the project. An extensive body of literature is reviewed to identify current modeling approaches and their theoretical and practical limitations. This literature study shows that, although discrete time models have been intensively investigated, there are few studies which propose continuous modeling of time. The container routing problem is modeled as a mixed integer program with explicit time variables and lateness penalties. A robust formulation is then proposed to eliminate the uncertain parameters from the objective function and constraints. By solving the new model exactly, with the aid of an optimization solver, robust solutions are obtained corresponding to transportation plans which remain feasible for any realization of the release times within the prespecified uncertainty interval. In order to introduce some flexibility in the transportation plan, the continuous time variables are modeled as affine functions of the uncertain parameters. The resulting twostage decision model is tested for a smallsized instance in both situations, with high and low lateness penalties. The computational results show that the adjustable robust model yields on the one hand, routedependent adjusted solutions for the case of penalized lateness, and on the other hand, a direct improvement of the objective function for the case of tolerated lateness. The results suggest that the adjustable robust optimization framework has sufficient potential to model the synchromodal container routing problem. This thesis concludes with addressing some of the limitations of the proposed model and indicating concrete approaches for countering them."Synchromodal; Robust; Optimization)uuid:6184bd6263d24929ae947f549149f5a5Dhttp://resolver.tudelft.nl/uuid:6184bd6263d24929ae947f549149f5a5{Topology Optimised Pedestrian Bridge: A feasibility study in using Topology Optimisation as a design tool for bridge design;van Baalen, Sander (Architecture and the Built Environment)nSmits, Joris (mentor); Eigenraam, Peter (mentor); Delft University of Technology (degree granting institution)5Research in using Topology optimisation as a tool to improve the design process of bridge design. The research focuses on implementing Topology optimisation in the design process of architects and engineers. As an example topology optimisation is used to propose a design for a pedestrian bridge in Amsterdam.4Topology; Bridge; Design; Optimisation; OptimizationBuilding Technology)uuid:36c9527eabd347bcafb86196e7d5e587Dhttp://resolver.tudelft.nl/uuid:36c9527eabd347bcafb86196e7d5e587Development and Validation of Shortterm Mine Planning Optimization Algorithms for a Sublevel Stoping Operation with Backfillingivan Duijvenbode, Jeroen (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)KMatthus, A. (mentor); Lottermoser, B. (graduation committee); Hennig, A. (graduation committee); Buxton, Mike (graduation committee); Rinne, M. (graduation committee); Delft University of Technology (degree granting institution); RWTHAachen University (degree granting institution); Aalto University (degree granting institution)$Sublevel Stoping; Optimization; MILP
20210731European Mining Course)uuid:888bf27a4dd44b08ab9a4b7c8ef62ef7Dhttp://resolver.tudelft.nl/uuid:888bf27a4dd44b08ab9a4b7c8ef62ef7WThe Aerogravity Assist: Combined Atmospheric and Interplanetary Trajectory OptimizationFodde, Iosto (TU Delft Aerospace Engineering)UDirkx, Dominic (mentor); Delft University of Technology (degree granting institution)Currently one of the major drawbacks of interplanetary travel is the cost of these types of missions. Mission to extraterrestrial bodies require a large amount of propellant to compl< ete, which often takes up a large part of the mass budget of the spacecraft. This requires larger launchers to be used to get the spacecraft into the required orbit, which increases the cost of these missions. An often used method to decrease the propellant used during an interplanetary mission is the gravity assist. The gravity assist uses the gravitational influence of a planetary body to increase its heliocentric velocity. This is done by exchanging some of the planets momentum as it orbits the Sun to the spacecraft. The increase in velocity is limited by several factors: the gravitational influence of the body, the closest distance possible between the body and the spacecraft, and the orientation of the planetary bodies relative to each other. A new method has been proposed by McRonald (1990) called the aerogravity assist, which improves upon the gravity assist by increasing the angle over which the spacecraft can turn using the aerodynamic forces it experiences while travelling through the atmosphere of the planetary body. This increase in the velocity bending angle increases the heliocentric velocity compared to the gravity assist, and also allows for changes in inclination of the interplanetary orbit using less propellant. This report investigates both the atmospheric and interplanetary trajectory of the aerogravity assist, and combines them to understand the benefit of the aerogravity assist for interplanetary trajectories.<br/><br/>A modular simulation environment is designed using the TU Delft Astrodynamics Toolbox to simulate the atmospheric trajectory of the aerogravity assist. This simulator allows the selection of different types of acceleration and environment models. The waverider is selected to be the vehicle used during the aerogravity assist, as it is able to obtain high liftoverdrag ratios which reduce the drag losses during the atmospheric section of the maneuver. The specific aerodynamic characteristics and heating models of the selected vehicle are evaluated and incorporated into the simulator. From this simulation environment, the optimal atmospheric trajectories were investigated using an optimization algorithm. It was found that this specific problem is extremely sensitive to the selected optimization algorithm, the tuning parameters of the optimization algorithm, and the objecive and decision variables that are used. The MOEA/D algorithm was selected to find the optimal atmospheric trajectories for Mars, Earth, and Venus, where the objective variables were chosen to be the atmospheric bending angle and the incoming and outgoing inclination difference. It was found that using the angleofattack and bank angle as control variables, an atmospheric bending angle of 120.5 degrees and an inclination difference of 0.2 degrees could be obtained at Mars. Using the control action instead of the inclination difference as objective variable at Mars showed that the atmospheric bending could be increased to 151.2 degrees while maintaining an inclination difference of below 3 degrees, and smoothed the trajectory in the process. For Earth and Venus the same optimization scheme was employed. They were able to obtain atmospheric bending angles between 70 and 50 degrees, while obtaining inclination differences below 1.5 degrees. For Venus, the heat load was also optimized instead of the inclination difference and it was found that the peak heat flux decreased by around 1000 W/cm$^2$, while decreasing the atmospheric bending angle by around 10 degrees. <br/><br/>A mission planner was designed to be able to investigate the effect of the aerogravity assist on interplanetary trajectories. To increase the accuracy of this mission planner, the numerically found atmospheric trajectories are used to determine the possible interplanetary trajectories. Pareto fronts were created for Mars, Earth, and Venus for a range of incoming velocities that showed the possible velocity bending angles and outgoing velocities that can be achieved using an aerogravity assist. These Pareto fronts were then implemented in the mission planner to deter< mine the influence of the AGA on an interplanetary trajectory. Several different interplanetary trajectories were investigated using both gravity and aerogravity assists. It was found that due to the fact that for Earth and Venus higher velocities are needed to travel through the atmosphere, the aerogravity assist would lose too much velocity and thus not improve upon the gravity assist. However, for Mars an increase was found of 4 km/s compared to the gravity assist for a trajectory to Saturn using flyby's at Mars and Jupiter. Taking into account the extra mass of the thermal protection system, the decrease of mass was found to be around 60 percent, which could be used to lower the cost of the launcher or increase the scientific output of the mission by adding instruments to the payload.<br%reentry; Optimization; gravityassit)uuid:ccc91e95a7904fbcbcafeffedc9dbd4dDhttp://resolver.tudelft.nl/uuid:ccc91e95a7904fbcbcafeffedc9dbd4djTowards the Industrialization of MDAO: Evolving 1st Generation MDAO to an Industry Ready Level of Maturity1de Vries, Danil (TU Delft Aerospace Engineering)La Rocca, Gianfranco (mentor); Binder, Simon (mentor); van Gent, Imco (mentor); Delft University of Technology (degree granting institution)BMultidisciplinary Design Analysis and Optimization (MDAO) promises major advances in aircraft design. However, judging from todays aircraft, MDAO is clearly not fully embraced by the industry yet. This is due to technical difficulties (coupling incompatible tools) and nontechnical barriers (intellectual property).<br/><br/>Recent research, such as the AGILE project, aims at enabling 3rd generation MDO, where collaboration of distributed teams is key. Maturing MDAO is important to eventually make Overall Aircraft Design possible. However, the industry is still in the early days of 1st generation MDO application.<br/><br/>This thesis proposes a new 1st generation MDO 'pipeline' benefiting from the latest tools developed for 3rd generation MDO. This opensource toolchain makes it easy to connect analyses and include gradient information. It is shown that gradient inclusion yields dramatically reduced computational costs. Rapid (re)configuration and inclusion of gradient information using this pipeline is demonstrated by considering the Sellar problem and a wing optimization.bMDO; MDAO; Optimization; AGILE; OpenMDAO; Framework; Multidisciplinary; Aircraft; Design; Analysis)uuid:d25e499e89de4fa5b96e3ada53141754Dhttp://resolver.tudelft.nl/uuid:d25e499e89de4fa5b96e3ada53141754#Optimization of Jetted Radial Wells=Troost, Diederik (TU Delft Civil Engineering and Geosciences)Jansen, Jan Dirk (mentor); Peters, Lies (mentor); Leeuwenburgh, Olwijn (mentor); Delft University of Technology (degree granting institution)Jetted radial wells are a relatively new form of hydraulic stimulation for oil and gas wells, where small diameters holes extend radially from a wellbore. This technique for increasing well connectivity could also be applied to geothermal wells. Radials can be freely placed along a backbone well, making the design of an economically viable well plan time intensive. Numerical optimization has been extensively applied to well design and will be applied for the design of a radial well plan. A simple optimization case is used to compare gradientbased and gradientfree optimization for jetted radial optimization, the gradientfree CMAES algorithm is chosen for the workflow. In addition to CMAES, mixed integer CMAES and uncertaintyhandling CMAES are incorporated in the workflow. Three reservoir cases are evaluated for optimization using these algorithms. The workflow is able to determine well plans with higher NPV in each of these cases though finding a global optimum remains difficult.LPetroleum Engineering; geothermal; Optimization; Jetted radial wells; CMAES
20190109)uuid:e095b2bec4974de4a55b25126e960dbeDhttp://resolver.tudelft.nl/uuid:e095b2bec4974de4a55b25126e960dbe@Robot Placement for Mobile Manipulation in Domestic EnvironmentslChawla, Hemang (TU Delft Mechanical, Mariti< me and Materials Engineering; TU Delft Biomechanical Engineering)Wisse, Martijn (mentor); Kober, Jens (graduation committee); Gaisser, Floris (graduation committee); Delft University of Technology (degree granting institution)
The development of domestic mobile manipulators for unconstrained environments has driven significant research recently. Robot Care Systems has been pioneering in developing a prototype of a mobile manipulator for elderly care. It has a 6 degrees of freedom robotic arm mounted on their flagship robot LEA, a nonholonomic differential drive platform. In order to utilize the navigation and manipulation capabilities of such mobile manipulators, robot placement algorithm that computes a favorable position and orientation of the mobile base is sought, which enables the end effector to reach a desired target. None of the existing approaches perform robot placement while ensuring a high chance of successful planning to target through a short path, while accounting for sensing and actuation errors typical in real world scenarios. This thesis presents a novel robot placement algorithm DeCOWA (Determining Commutation configuration using Optimization and Workspace Analysis) with these characteristics. Since the approach to robot placement is dependent upon the kind of mobile manipulation, a comparative study of sequential and full body methods is performed with respect to criteria important in domestic settings. Sequential mobile manipulation is found to be most suitable, for which a modular mobile manipulation framework encompassing motion planning and robot placement is presented. With sequential mobile manipulation, the ability to successfully reach a target depends upon the kinematic capabilities of the arm. Accordingly, robot placement with DeCOWA determines a favorable location for the arm, and corresponding platform orientation. To find the position of arm s base, an offline manipulator workspace analysis is performed generating the Inverse Reachability and Planability maps. During online use, these maps are combined into an Inverse Fusion Map that ranks different<br/>locations based on the ability of the arm placed there to find a successful and short motion plan to target. This map is filtered to generate a set of feasible locations at the arm s height. Through a ranked iterative search, a suitable collision free arm location is determined followed by minimization of the platform distance from robot s current pose. This approach is evaluated against an unbiased random placement of robot near the target using a sample set of twenty scenes mimicking domestic settings. It is found that DeCOWA is able to generate commutation configurations in fraction of a second, that lead to a high planning success rate, a short path length, and account for goal tolerance of navigation. Also, its modularity allows to use several planability metrics, making it useful for domestic application.Mobile Manipulation; Robot; Robot Placement; Vision; Domestic Environments; Commutation Configuration; Robotics; Base Placement; Reachability; Workspace Analysis; Optimization; Planability; Algorithm
20191110"Biomechanical Design  BioRobotics)uuid:b7911d199b28445ebf1905db3b9f85ddDhttp://resolver.tudelft.nl/uuid:b7911d199b28445ebf1905db3b9f85ddSStructural Design Exploration of Highrise Design Concepts Using Genetic AlgorithmsASchagen, Roland (TU Delft Architecture and the Built Environment)wTurrin, Michela (mentor); Oikonomopoulou, Faidra (mentor); Delft University of Technology (degree granting institution)The adoption of optimization algorithms has been slower in conceptual design compared to the later design phases. The choices made during the conceptual design phase have a large impact on the performance of the result. This research aims to find out why it is so hard to use optimization in conceptual design. Then a method is proposed to solve this issue during the application of optimization on the conceptual design of a case study. As subject of the case study the CCTV building in Beijing is chosen as it is a good example o< f a building were the concept development is nebulous, the shape is difficult to rationalize, and a small change in the concept can have large impact on the overall performance of the result. The case study was used to develop a method for optimizing design concepts. This method involves the parameterization of the CCTV building envelope, the calculation of its structural performance using an abstracted finite elements analysis, and the use of NSGA II. The result is a deeper insight into why optimization is difficult in the conceptual phase, a method that is considerate of these difficulties, and an analysis tool for this method.Genetic Algorithm; Illdefined problems; Optimization; conceptual design; NSGA II; FEM analysis; CCTV tower; structural analysis39.91347, 116.45805)uuid:2b6974777d594083b51d70c0d3137261Dhttp://resolver.tudelft.nl/uuid:2b6974777d594083b51d70c0d3137261)Time Optimization for Laser Sailing RacesQPardo Arteaga, Gabriela (TU Delft Mechanical, Maritime and Materials Engineering)VHeemink, Arnold (mentor); Delft University of Technology (degree granting institution)Sailing against thewind follows a zigzag trajectory. To find the optimal time and pathtrajectory to move from one target to another during a race, this study analyzes how the wind determines the optimal timepath of Lasers, one of the smallest sailboats that compete during the summer Olympic Games.<br/>To answer this question, this study usesfour wind models in an algorithm developed in MATLAB to find the optimal timepath. The wind models were forecasts and wind measurements for the area of race R1 from the World Cup Series 2018 at Hyres, France. One of the wind models used were a Weather Research and Forecasting model (WRF) with a grid resolution of 1km, a time step of 10 minutes. On the other hand, the most basic model was a constant and uniform wind field. The race R1 hasthree lines, limited by two buoys, andone point(buoy), one line and one point define a leg. R1 hasfivelegs and two of them areagainst the wind. The results, time and pathtrajectories, of each of the wind models,were comparedwith the results of the top 10 winners of the race. They showed that the legs sailed against the windare also characterizedby the location of the sailboats on the start line. The times of these legs using the WRFwind model and the racetime had an error of less than 5%.Forthe prediction of the startlocation, it was the same as the winner of the race. However, the direction of the pathswas not predictedaccurately for these legs. Using the constant and uniform windscenario, the percentage error of the racetime respect to the winner is about 7%.However, the direction of leg 2 is not even similar to the winner. To review the effects of the heights of the waves this study proposes to model the sailboat in 3dimensions including the Xcoordinate of the sailman position. In addition, a 3D modelallowstheanalysisofhow the center of effort (CE) on the sailis affectedbythe current and waves.(Sailing; time; Optimization; Laser; WindMechanical Engineering  BioMechanical Design)uuid:04f50f3b85724c68b75be6990c639be0Dhttp://resolver.tudelft.nl/uuid:04f50f3b85724c68b75be6990c639be0WDesign and Implementation of a Development Platform for Indoor Quadrotor Flight ControlNavia Vela, Libardo (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Delft Center for Systems and Control)nKeviczky, Tamas (mentor); Zheng, Zixuan (mentor); Delft University of Technology (degree granting institution)Research and development of flight control of quadrotor Unmanned Aerial Vehicles (UAVs) have gained popularity during the past few years due to their deployment flexibility and wide range of applications: agriculture, entertainment, cinematography, package delivery, search and rescue, etc. This thesis project follows up this interest aiming to build a realtime development platform for experimentation and test of indoor flight controllers in the Networked Embedded Robotics in Delft laboratory. For this goal, we employ a Matrice < 100 quadrotor from DJI alongside its ROS C++ Software Development Kit (SDK), develop a system identification routine to model its autopilot, vertical and horizontal displacement, and based on the attained results, we design and implement a control topology for realtime position control of the UAV inside the laboratory by means of a Model Predictive Controller (MPC) and a LinearQuadratic Regulator. Finally, we enunciate and discuss potential applications in which this development platform can be used.0UAV; MPC; Optimization; ROS; C++; LQR; Quadrotor)uuid:2694dcb5292f41e2a9f7e79b0dccc8c0Dhttp://resolver.tudelft.nl/uuid:2694dcb5292f41e2a9f7e79b0dccc8c0tInventory and Modelling of Different Objective Functions and their Impact on Optimal Design of an Offshore Wind FarmRMurali, Vishal (TU Delft Electrical Engineering, Mathematics and Computer Science)Watson, Simon (graduation committee); Zaayer, Michiel (mentor); Langelaar, Matthijs (graduation committee); Delft University of Technology (degree granting institution)The design of an offshore wind farm (OWF) is multidisciplinary in nature as it involves the design of many disciplines such as the wake effects, support structure, electrical cables etc. For the optimal design of an OWF, an optimization procedure is required where all the disciplines are optimized simultaneously. The objective function plays a significant role in optimization as it expresses the main aim of the model which is to be either minimized or maximized. So far, cost of energy (COE) and annual energy production (AEP) are one of the commonly used objective functions for OWF optimization as far as the author is aware. However, there might be other objective functions that may influence the optimal design of an OWF as well. This may include maximizing the profit, minimizing the environmental impact, reducing their carbon emissions etc. Hence, this thesis investigates the overview of different objective functions and understand its impact on the optimal design of an OWF.<br/><br/>An inventory of different objective functions is prepared, and relevant ones are selected for further study. It is observed that even though some objectives are dissimilar, they still depend on the same wind farm parameters and are therefore expected to give similar design results. From the list of objective functions, net present value (NPV) and risk management objectives are chosen for further research.<br/><br/>The selected objective functions are then formulated in a metric for optimization. The price of electricity plays a significant role in determining the NPV. It is learnt that electricity price varies with the power supply depending on the site conditions. The electricity price is low if the supply of power is high in a region where there are many OWF s and vice versa. Moreover, OWF investors value constant power output without any fluctuations. Hence, taking all these aspects into consideration, the electricity price in the NPV function is modelled for a constant value, wind variability and wind power predictability.<br/>The risk management function, on the other hand, aims at minimizing the uncertainty associated with an OWF project. The risk here refers to the uncertainty associated with the profit obtained from the OWF. A set of annual average wind speeds is computed using monte carlo simulations and the AEP and NPV are estimated. The mean( NPV_mean) and standard deviation ( NPV_std ) of NPV are then calculated. NPV_std represents the uncertainty in this scenario and is minimized to reduce the risk.<br/><br/>A suitable method is then identified to deal with multiple objectives. The NPV function is maximized for maximum profit and this objective is evaluated using a single objective optimization technique. The risk management objective involves the calculation of NPV_mean and NPV_std. Both objectives are contrasting in nature as a significant reduction in NPV_std corresponds to an undesirable reduction in NPV_mean. A tradeoff between both these objectives is the best possible solution. Therefore, a multi objective optimizat< ion technique is used, and a list of solutions is obtained by generating a pareto front.<br/><br/>The new approach is then evaluated by implementing different case studies. It is observed that optimum rotor diameter and number of turbines for the single objective optimization technique are influenced by economic indicators such as the real interest rate and lifetime. However, they are not influenced by variation in the electricity price. Nevertheless, the NPV function is sensitive to the economic indicators and variation in the electricity price. <br/><br/>For the multi  objective optimization technique, multi criteria analysis was used to determine the weight to the objective functions while moving along the pareto curve. It was observed that the improvement of one objective led to the deterioration of the other objective. Hence, the pareto front provides opportunities to investors to negotiate and decide on the weight they want to specify for their objectives.<br1Wind Energy; Optimization; wind farm optimization)uuid:8bd18483c1e1460b9f590f5e20726679Dhttp://resolver.tudelft.nl/uuid:8bd18483c1e1460b9f590f5e20726679rThe financial implications of building design: a developer and consumer preferencebased residential design systemIn 't Hout, C.E.2Binnenkamp, R. (mentor); Van Bortel, G.A. (mentor)aAs of yet, there are virtually no practical tools available for incorporating both developer and consumer preferences within a floorplan design. This study aims to deliver a proof of concept for such preferencebased design systems, and their application within the commercial property development process. A research model is presented that is able to generate residential floorplan design alternatives, which are optimized for both developer and consumer preferences. At the foundation of this model lies the concept of Preference Function Modelling. With the use of a case study, the research model is shown to reflect reality to an accurate degree. The results indicate that the research model has generated floorplan design alternatives that are associated with a higher aggregated preference score than the actually realized floorplans for the case study. As this was the intended purpose of the research model, a proof of concept has been delivered. In conclusion, further exploration of the potential utility and applications of similar preferencebased design systems is encouraged by the outcomes of this study.PreferenceBased Design System; Preference Function Modelling; Operations Research; MultiActor MultiCriteria Analysis; Design Alternatives; Optimization; What sBest!&Architecture and The Built Environment#Management in the Built Environment)uuid:dec3c264d6524370864c1f9ea18cbd2aDhttp://resolver.tudelft.nl/uuid:dec3c264d6524370864c1f9ea18cbd2a5Reservoir Characterization using a Geometric ApproachVKhandelwal, Anshul (TU Delft Electrical Engineering, Mathematics and Computer Science)Eisemann, Elmar (mentor); Hildebrandt, Klaus (mentor); van de Giesen, Nick (graduation committee); Liem, Cynthia (graduation committee); Delft University of Technology (degree granting institution)XQuantifying the anthropogenic impacts such as reservoir characterization is a big challenge in the field of water management. In this work, a computer graphics based geometric approach is presented which can predict the underlying topology of largescale reservoirs. The proposed algorithm uses freely available, satellite based landscape data of the surrounding regions to predict reservoir characteristics. The premise of the presented approach is that the slope of the surrounding landscape is an important determinant to understand the underlying landscape of the reservoirs. This method outperforms the existing stateoftheart techniques used to estimate the storage capacities drastically, both in terms of estimated maximum volume stored and estimated volume area curves. Evaluation of the geometric model presented is done on 28 reservoirs using the HydroSHEDs data which was developed using the Shuttle RADAR Topography Mission conducted by NASA. This < HydroSHEDs data was obtained in 2000 which acts as ground truth data for the reservoirs built after 2000. Further, model parameters are introduced to improve the modeling capabilities of the reconstructed reservoirs. This approach further intensifies the case of using computer graphics techniques for raster based analysis and provides a platform for further research in the field of water management.<brgReservoir Characterization; Computer Graphics; Optimization; Water Management; Digital Elevation Models)uuid:c49d975200524db0af2015477951f6fdDhttp://resolver.tudelft.nl/uuid:c49d975200524db0af2015477951f6fd4Optimization of a wind turbine bladeroot connectionVRao, Athul (TU Delft Aerospace Engineering; TU Delft Aerospace Structures & Materials)~Shroff, Sonell (mentor); Bellotto, Valter (graduation committee); Delft University of Technology (degree granting institution)For several years Tbolts have been a popular choice for joints in the field of wind energy, specifically for connecting the blade roots to the hub of the wind turbine. Their use is mandated by the geometry of the joint and they perform very well under pure axial loading. However recent analyses have shown significant bending stresses in the Tbolts preventing their use to their full capacity. These bending stresses are unavoidable due to the presence of a slew bearing between the blade root and the hub. The bending stresses generated while loading the blade root, causes the blade rootbearing joint to gradually open, causing excessive loading of the Tbolt above a certain load. <br/>It is hypothesized that modifying the blade root design to reduce the effects of local bending can open up the possibility of reducing its mass and cost. To test this hypothesis, the blade root is initially studied and the stress ratio is identified as an appropriate joint performance parameter. The performance of the joint is boosted by increasing the pretension of the bolt. After an initial phase of over designing the joint to reduce the constraint stresses, the joint optimization is carried out using the Sequential Quadratic Programming algorithm. The optimization culminates with the mass reducing by roughly 110Kg and the material cost reducing by approximately 13% per blade root. The number of bolts reduces from 88 to 52. Thus, a simpler design is achieved, that promises simpler and cheaper manufacturability, higher reliability and lesser sites for crack nucleation in the laminates. The current design strategy at Suzlon is to employ a greater number of Tbolts with thinner shanks. Curiosity in the field of cost optimization that initiated from within Suzlon has proved that there exists a different design strategy that holds great promise for delivering structurally equivalent if not better designs with improved cost, mass and reliability.7WInd turbine; Blade root; Tbolt; Optimization; Fatigue
2022091452.290023, 6.781693)uuid:fa661ab603244115adc80c8665e2e1d9Dhttp://resolver.tudelft.nl/uuid:fa661ab603244115adc80c8665e2e1d9Technoeconomic assessment of grid connected Power to Heat and Power to Hydrogen technologies: An electricity market scenario approachVChintada, Priyanka (TU Delft Electrical Engineering, Mathematics and Computer Science)Zaayer, Michiel (mentor); MoralesEspana, G (mentor); Watson, Simon (graduation committee); de Vries, Laurens (graduation committee); Delft University of Technology (degree granting institution)= Wind energy is the most promising renewable energy source for the Netherlands in the energy transition period. However, the investment costs of the offshore wind are huge. Until now, the government was supporting the investments in wind farms with subsidy schemes irrespective of the electricity market dynamics. The latest largescale offshore wind projects in Germany and the Netherlands (Hollandse Kust Zuid 1 and 2) have been awarded zerosubsidy. However, in the Netherlands, the transmission costs of the offshore wind electricity are socialized. This would force the unsubsidized windfarms to feed their electricity into the public gri< d exposing them directly to the whole sale electricity market. With large shares of wind penetration in the market, the value of the electricity is declining. At the same time, the intermittent nature of the renewable energy sources generates large demand for flexibility in the electricity system. These pros and cons of renewable energy sources create new opportunities for electrification of other sectors to support achieving Paris climate agreement goals of the country. Heat and hydrogen sectors are identified as potential flexibility options to facilitate the decarbonization of the energy system by 2050. European electricity market simulation model COMPETES is identified as potential tool to assess the integrated systems at national level. However, the flexibility provided by heat and hydrogen demands with thermal and hydrogen storage technologies is not part of the model directly. Therefore, an optimization model is developed to determine the size of power to heat and power to hydrogen systems with least cost. The results of this problem determine the outlook for investments in power to heat and power to hydrogen technologies. The problem also determines the operational strategy of these units in response to the hourly power prices. Consequently, the coupled implications of the emerging electricity load from heat, hydrogen markets on the national power system is studied. Industrial heat, domestic district heat networks currently provided by combined heat and power gas plants in the Netherlands and the developments in hydrogen consumption market are taken as study cases for assessment. This work provided insights about the opportunities, threats and solutions for the involved stakeholders.TPower to heat; Power to hydrogen; Electricity market model; Optimization; Renewables)uuid:621c77c467644d179a0ab3d3dd3025e0Dhttp://resolver.tudelft.nl/uuid:621c77c467644d179a0ab3d3dd3025e0;Guidance and Control for Spinning, GunLaunched Projectilesvan der Geest, Lars (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Delft Center for Systems and Control)van Wingerden, JanWillem (mentor); Bouquet, Franois (mentor); Hellendoorn, Hans (graduation committee); Keijzer, Twan (graduation committee); Delft University of Technology (degree granting institution)
This thesis report is focused on the design of a guidance and control system that is able to minimize the deviation from the desired impact point for a firing range of 1 kilometer, given a 30mm spinning gunlaunched projectile with a novel actuator design. This actuator is fixed on the projectile and offers a single force that can only be switched on or off. <br/><br/>In order to test the effectiveness of the guidance, navigation and control (GNC) system, an accurate model for both the projectile and the proposed actuator are constructed. These models will serve as a substitute for testing with working prototypes. The model for the projectile dynamics is constructed based on existing nonlinear 6DOF rigidbody models that are widelyused and validated. A new, simplified second order model that approximates the dynamics of the actuator, based on available measurements, is constructed.<br/><br/>The overall structure of the GNC loop is defined, with the guidance method of choice being the precalculation of a reference trajectory towards the intended target. This method serves to alleviate underactuation and computation time problems. One of the most important contributions of the thesis is made with the description of a method for transforming the input from a single binary input signal towards two continuous virtual input forces. This method uses a discretization of the rolling motion of the projectile, such that an optimization can be done which results in a transformation from the binary on/off signal to two virtual force inputs. These two virtual force represent the steering forces in the directions perpendicular to the forwards motion that would have the same effect as the binary on/off signal if directly applied on the projectile.<br/><br/>The restructuring of the < input allows for the use of PDcontrollers to track the predefined trajectory. Using a gridsearch method, the controller gains are found that best satisfy the design goal of minimizing the dispersion, which is defined by the sum of the mean error and the standard deviation of the impact points.<br/><br/>Analysis of the integrated solution shows that the proposed solution with the input transformation and the PDcontroller is able to significantly reduce the projectile dispersion from standard deviations of about 50 cm to just a few cm. The exact performance depends on the frequency of the actuator and the spin rate of the projectile. The best performance is reached with high spin rates and actuator frequencies high enough to match these spin rates.7Guidance; Control; Projectile; PD control; Optimization
20240221)uuid:86c056d8f3684239893f07ca3a22e112Dhttp://resolver.tudelft.nl/uuid:86c056d8f3684239893f07ca3a22e1121The Moving Node Approach in Topology OptimizationOvervelde, J.T.B.(Not available because of confidentialityTopology; Optimization
20130418)uuid:3e7cb93b7d3b4b33bcc0d562630f12e7Dhttp://resolver.tudelft.nl/uuid:3e7cb93b7d3b4b33bcc0d562630f12e7)Ready for TakeOff: Scaled Flight TestingVan Gorcum, J.J.Veldhuis, L.L.M. (mentor)Atmospheric free flight scaled flight testing is an affordable way to investigate the dynamic properties of an aircraft, while enabling a wider range of test possibilities than a windtunnel. This research is the first step into the development of a scaled flight testing model and as such will highlighted the difficulties and discrepancies which will be faced when scaling is performed. This study has the objective of investigating the effect of scaling on the aerodynamic properties and of the flow field at the tail location, with the focus on takeoff conditions. The full scale case is a regional aircraft, the ATR72, with an unswept and slightly tapered wing. The scaled case is 14.7% geometric scaled The analysis of the aerodynamic properties is split between the clean wing and the wing with an additional surface in the form of a flap. The aerodynamic properties investigated were the lift, moment and drag curve, and the maximum lift coefficient. The clean wing was analysed using a quasi3D analysis called Q3D, which is a combination of a vortexlattice method (AVL) and a vortex panel method (XFOIL). A modification of XFOIL, called RFOIL, is used to analyse the difference in maximum lift coefficient between the airfoil of the scaled and full scale case. RFOIL is selected because of the better results near stall. The wing with flap is investigated using a vortex lattice method, in this case again AVL, where the airfoil with flap is investigated using the Eulersolver MSES. Finally the maximum lift coefficient for the wing with flap is analysed with the semiempirical Pressure Difference Rule. The Pressure Difference Rule states that there is a ratio between the peak pressure and the trailing edge pressure, on a surface, at which maximum lift occurs. The methods are limited in their incorporation of viscous effects, the 2Danalysis tools (XFOIL, RFOIL and MSES) only include viscous effects in the small boundary layer region. AVL does not include viscous effects directly, but only via a correction factor for the liftcurve slope. However, the results found in literature using these methods are satisfactory. An analysis of the wake field is done using methods developed by the ESDU. The wake properties of both the scaled as well as the full scale are calculated using semiempirical models and use a simplification where to place the vortex sheet. The dimensions and velocity loss of the wake itself can be calculated using either a method by ESDU or a method by Schlichting. Results showed that both these methods provided similar solutions. The scaling of the wing proves to change the aerodynamic properties and the similarity is no longer present for all the investigated aerodynamic properties. For both the wing with and without flap this is the case. These differences are mainly due to< Reynolds number effects, the Mach number effects are only minimal on the lift coefficient. Due to the limits on the selected methods the exact magnitude of the difference found can not be guaranteed, but the trends in the found differences are certain. A basis for this is found in the different boundary layer properties, the scaled wing exhibits a thicker boundary layer, leading to a decambering effect, and laminar separation bubbles are found to occur. Both the clean and wing with flap shows more outboard loading for the scaled wing. The scaled wing exhibits a thicker wake, with a larger velocity loss. The difference between the full and scaled clean wing is thus larger than the difference between the full and scaled wing with flap. The reason for difference between the full and scaled wing can be found in the increased drag (coefficient) of the scaled wing. The reason behind the larger increase in wake for the clean wing is due to the fact that the difference in drag coefficient between the full and scaled case is larger for the clean wing than for the wing with flap. An investigation to minimize the differences between the full and scaled wing is done. It is decided to change the airfoil shape. Improvements are visible for the optimized airfoil, on all the aerodynamic properties. The optimized airfoil shapes tends towards a thinner airfoil and thus an investigation into solely optimizing the thickness of the airfoil is also done. This also proves to give better results than the original airfoil, however, not as good as the shape optimized airfoils.. This thinner profile was also investigated in the configuration with a flap, but here it proved to decrease the performance of the airfoil with flap. The difference between the full and scaled case can be reduced by an optimization. If an optimization is to be done, it must be done on the whole of configurations and flight conditions. As the next step into a development of a scaled flight test model, the exact extend of tests must be determined, only then can it be investigated how scaling affects the results of testing. No direct solutions exist to completely overcome the gap between the full and scaled aerodynamic properties of a flight testing model, however airfoil shape optimization does provide a better similarity.zScaled Flight Testing; Aerodynamic properties; Wake; Flap; Similarity; Scaling; Optimization; MSES; RFOIL; XFOIL; Q3D; AVL)uuid:cc24b8b9fca64387a746f7b50bb32c73Dhttp://resolver.tudelft.nl/uuid:cc24b8b9fca64387a746f7b50bb32c73lQuantifying TakeOver Quality: A Novel Approach Utilizing a ScenarioSpecific Optimized Reference TrajectorygLoosveld, Erik (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Intelligent Vehicles)Happee, Riender (mentor); de Winter, Joost (graduation committee); Doubek, Fabian (graduation committee); Delft University of Technology (degree granting institution)Due to recent advancements in automated driving systems, drivers are able to withdraw themselves from the control loop of the vehicle in certain driving situations. However, system limits occur relatively frequently, which results in the driver having to take back control. It is argued most currently used takeover quality quantification methods, such as identifying the minimum timetocollision or maximum lateral acceleration, neglect inherent criticality of the driving scenario. This can lead to poor overall quality assessment, especially in challenging takeover scenarios. This thesis proposes a novel method of quantifying how well the driver performed in the task of executing an emergency automationtomanual control takeover, by means of a scenariospecific quality reference in the form of an optimized trajectory. Firstly, a humanintheloop driving simulator study was carried out using a high fidelity 6DOF driving simulator (n=25). The simulated environment consisted of a twolane highway setting. During the drive, six different automationtomanual takeover scenarios occurred. Optimized trajectories were also generated for these scenarios, using a Stateof< theArt trajectory planner developed at Volkswagen Group Research Facility, Wolfsburg. The independent variables were the time budget (5s, 7s, and 20s), and the traffic density (5 veh/km, and 10 veh/km). A variety of dependent variables were determined for both the human driver trajectories and the optimized trajectories: minimum timetocollision, minimum distancetolaneboundary, minimum distancetoovertakingvehicle, maximum forced braking of the overtaking vehicle, maximum and standard deviations of accelerations (lateral and longitudinal), maximum and standard deviations of jerk (lateral and longitudinal), and lateral quickness. Based on the findings from this humanintheloop driving simulator study and the generated optimized trajectories, a quantification framework was developed. The framework subdivides takeover quality assessment into two separate parameters. One with regard to safety, and one with regard to comfort. The safety parameter quantifies criticality of the human driver takeover based on the sum of normalized weighted quantification scores. These scores are computed for the metrics minimum timetocollision, minimum distancetolane boundary, minimum distancetoovertaking vehicle, and maximum forced braking of the overtaking vehicle, and are based on the difference to the metric values corresponding to the optimized<br/>trajectory. The comfort parameter quantifies comfort in a similar manner based on the metrics lateral quickness, maximum lateral and longitudinal jerk, and standard deviation of lateral and longitudinal jerk. The final quantification scores are comprised of a value between 0 and 1 for both the safety and comfort parameters. A perfect score of 0 would indicate all metrics were identical, or better than the optimized reference trajectory. A value of 1 indicates a crash for the safety parameter, and highly uncomfortable driving for the comfort parameter. The resulting quantification framework is scenariospecific as the quantification is performed relative to an optimized reference trajectory.`Automated driving systems; takeover requests; Quantification; Optimization; trajectory planning
20240601)uuid:044eec5fc3f94021a7d81f3303d3492fDhttp://resolver.tudelft.nl/uuid:044eec5fc3f94021a7d81f3303d3492fqThermoeconomic Optimization of Organic Rankine Cycle (ORC) for Low Grade Heat Recovery at Temperatures below 100!Arumugam, Senthil Kumar (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Technology, Policy and Management; TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Process and Energy)Hemmes, Kas (mentor); Woudstra, Theo (graduation committee); Blok, Kornelis (graduation committee); Purushothaman Vellayani, Aravind (graduation committee); Stougie, Lydia (graduation committee); Delft University of Technology (degree granting institution)Lowgrade heat sources are abundant on earth but are majorly untapped due to lower thermodynamic efficiency at low temperatures and cost considerations. A costeffective technology is needed to convert this energy resource into useful forms of energy. This work aims at optimizing Organic Rankine Cycle (ORC) based heat engine and a cogeneration system developed to produce electricity and refrigeration from a heat source below 100!, from both thermodynamic and economic point of view. Exergoeconomics, an algebraic thermoeconomic method, was used to analyze and optimize the systems for costeffectiveness and exergetic efficiency. Also, the prototype of the cogeneration system was experimentally tested. The results exergoeconomic optimization show that the costeffectiveness of the cogeneration system can be significantly improved by design parameter changes. The experimental results obtained were comparable with the results obtained from theoretical simulations.Uexergoeconomics; thermoeconomics; lowgrade heat; organic rankine cycle; Optimization
20201025)uuid:72d06941c9754d518d4ba976a7c0b9cfDhttp://resolver.tudelft.nl/uuid:72d06941c9754d518d4ba976a7c0b9cfOA practical maintenance task packaging mode< l applicable to aircraft maintenance.Witteman, Max (TU Delft Aerospace Engineering)uDeng, Qichen (mentor); Lopes Dos Santos, Bruno (mentor); Delft University of Technology (degree granting institution)3This research concerns the problem of scheduling aircraft maintenance tasks, that must be carried in multiple maintenance checks to keep the aircraft airworthy. During this research the allocation of maintenance tasks to their maintenance opportunities is referred to as the task allocation problem. It is a complex combinatorial problem that is solved daily by aircraft operators.We propose a novel twostage framework capable of solving the task allocation problem, for an entire fleet. A MixedInteger Linear Programming (MILP) formulation was developed for both stages of the framework. In the first stage an exact method is used to predetermine the workforce to allocate per day to each aircraft under maintenance. In the second stage, the task allocation problem is solved in dependently for each aircraft, allocating task at the work shift level and per workforce skill. Moreover, an approximation algorithm is proposed, using bin packaging problem solution techniques, for both stages of the framework.Both stages were tested and validated using data of 45 aircraft from a European airline. The compu tational performance of the approximation algorithm is benchmarked using a MILP formulation, solved with a commercial solver. Results indicate that the solution quality of the proposed algorithm remains excellent and runs up to 26 times faster than the exact method.This research is the first to present an optimization model to solve the maintenance task allocation problem, being solved at the work shift level. Furthermore, it presents a binpackaging approximation algorithm that can be solved within minutes for the 45 aircraft fleet and a planning horizon of 4 years. The framework was validated in practice and the results give insights on the potential benefit of using such framework to solve this daily recurrent problem.>Maintenance; Optimization; Task allocation; Practical; Airline Aerospace Transport & Operations)uuid:2cc94bdd857d48fe91405b7e8cfa001fDhttp://resolver.tudelft.nl/uuid:2cc94bdd857d48fe91405b7e8cfa001f+Optimal Wind Farm Lifetime Power Production+Cheng, Ken (TU Delft Aerospace Engineering)Odgaard, Peter Fogh (mentor); Bierbooms, Wim (mentor); Kim, Taeseong (mentor); Jrgensen, John Bagterp (mentor); Delft University of Technology (degree granting institution); Denmark Technical University (degree granting institution)This project is a collaboration between Goldwind, Technical University of Denmark (DTU) and Delft University of Technology (TUDelft), of which focusing in developing an optimization model that could determinate the optimal operation schedule for wind turbine under the uncertainty of the electricity price and limitation in lifetime equivalent fatigue loads, for the purpose of improving and maximizing the revenue or cash flow for wind farm project.<br/><br/>The project comprises of two main parts, including the design and validation of up and derated operation modes, and the model development and optimization of the revenue. The design of the different operation modes provide varies rated power without sacrificing the lifetime structural integrity, as such offer extra flexibility in wind turbine operation schedule according to the fluctuation in electricity price. Lifetime equivalent loads would be calculated for the different operation modes design and a weighting factor, that indicate the rate of fatigue consumption relatively to the normal operation, would be computed and used in the optimization process. An optimization model is developed based on the turbine structural constraints as well as the financial profile. The optimization model processes the electricity price data and other wind turbine data, for example the fatigue consumption, to reproduce the optimal operation schedule and operation price thresholds that would yield utmost financial benefits.<br/><br/>The optimizer has a proven performa< nce improvement in revenue generation of 1.5% and 5% reduction in the project related interest payment under one of the case study under electricity price uncertainty. The model is effective in scheduling the turbine operation upon the varying electricity price. This optimization model is providing game changing insights and operation strategies that could be used throughout the turbine lifetime under the ever changing electricity price of the energy market. And further drive down the cost of energy in the highly competitive market. By improvising continuous input to the optimizer including bending moment sensors inputs, the optimization process could reproduce a more accurate revenue improvement.;Wind Energy; Fatigue Analysis; Optimization; Cost of Energy
20230816#European Wind Energy Masters (EWEM))uuid:ab92c1fba5b14302b2b4fcd36b84f0dcDhttp://resolver.tudelft.nl/uuid:ab92c1fba5b14302b2b4fcd36b84f0dcHStructural topology optimization of an active motion compensated gangwaytVergeer, Michael (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Marine and Transport Technology)Miedema, Sape (graduation committee); Langelaar, Matthijs (graduation committee); Helmons, Rudy (graduation committee); Kromwijk, PJ (mentor); Delft University of Technology (degree granting institution)uIn the offshore industry, there is a growing demand for designing efficient, sustainable and competitive products. In order to fulfil the component requirements, a method named topology optimization can be applied. This is a mathematical design method which can be used in the early phases of the design process.At IHC there is an interest of applying topology optimization for their equipment development process. Therefore the possibilities and limitations of the method should be investigated thoroughly.<br/><br/>The research covers the optimization process of a motion compensated gangway. A motion compensated gangway is a walkway which can be used to provide access from the transport vessel to the offshore structure. Its function is to transport people and cargo safely from the ship to the offshore structure or vice<br/>versa. The goal of this research was to determine to what extent topology optimization can be used in the design of a motion compensated gangway. Finding an optimized result in terms of weight and stiffness by using this mathematical method which satisfies all the requirements. The structural optimization is carried out with several commercial software packages which are compared by using a multicriteria analysis.<br/><br/>During the optimization process it has been found that there are essentially two stages in the optimization process. In the first stage, the topology or beam orientation of the structure is defined by the topology optimization process. In this part the concept of the design is generated. Variation of the optimization parameters was used in order to develop an efficient structure. The objective for the optimizer was to minimize the compliance of the structure for a certain volume fraction.<br/><br/>In the second stage, the dimensions of all the beams and elements are defined by performing a size optimization. A line model is generated which represents the orientation of the members in the structure. During the size optimization the shape and the dimensions of the members are defined in order to fulfill the<br/>objective. The objective is to minimize the mass of the structure while constraints are defined for the maximum allowable stresses in the members and the maximum vertical deflection of the structure. This postprocessing step is required in order to obtain a feasible design. The structural stability of the gangway was improved by performing a linear buckling analysis and by adapting the structure in order to reduce the buckling behaviour.<br/><br/>In the final step of the optimization process, a CAD drawing is generated. This model is analysed by performing a finite element analysis. This showed that the new optimized design satisfies all the requirements which are stated by the DNV for designing a< motion compensated gangway. The combination of the topology and size optimization resulted in a new design which yielded a weight reduction of 36,4% compared to the current design. The weight was reduced from 13,08 ton to 8,31 ton, while still satisfying all the constraints. Therefore it can be concluded that this methodology can be used in the optimization of a motion compensated gangway structureOptimization; Topology optimization; Shape optimization; Motion compensated gangway; Structural Optimization; Optistruct; Hyperworks; Size optimization; Design methodology
20201009!Offshore and Dredging Engineering)uuid:98fd5f3841d744b885e8c54e5aeb9c40Dhttp://resolver.tudelft.nl/uuid:98fd5f3841d744b885e8c54e5aeb9c40aExplore Pole Vaulting Strategies by control optimization: A biomechanical model based explorationPvan den Bosch, Maarten (TU Delft Mechanical, Maritime and Materials Engineering)van der Helm, Frans (mentor); Bergsma, Otto (mentor); Schwab, Arend (mentor); Delft University of Technology (degree granting institution)The development of the pole vaulting record is only improved by 1 centimeter over the last 25 years. Multiple studies show signs of different control strategies in elite pole vaulters, without pointing out what the cause of the differences is The aim of this research is to explore pole vaulting strategies by a control torque optimization. The optimization is performed for a simplified biomechanical model. The vaulter contains three segments: arms, a trunk including the head and legs. The joints are represented by torque actuators, located at the hands, shoulder joint and hip joint. The torque profile of each actuator is optimized. The optimization starts just after the pole is planted. It ends when the vaulter releases the pole. The pole length, pole stiffness, the vaulter s length, the vaulter s mass and the takeoff angle are varied to discover the influence of these parameters on the control strategy. Three control strategies are found, as well as two power management strategies. Global trends are found. An increased pole stiffness, a decrease in length and mass of the athlete and a decreased takeoff angle improved the performance. A possible optimal pole length is found.5Pole Vault; Optimization; MATLAB; Biomechanical ModelBMechanical Engineering  BioMechanical Design  Sports Engineering)uuid:2a15115f6d704057abdaeabf1adc72a7Dhttp://resolver.tudelft.nl/uuid:2a15115f6d704057abdaeabf1adc72a7TThe Improvement of Simulation Correlation for Formula 1: With a focus on kerb ridingbKoek, Bart (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Systems and Control)Hellendoorn, Hans (mentor); Alirezaei, Mohsen (mentor); Shyrokau, Barys (graduation committee); Delft University of Technology (degree granting institution)The highly competitive nature of Formula one makes every team eager to find areas of improvement on their cars. In the case of Scuderia Toro Rosso, it was found that there is a gap to improve performance while driving on kerb stones. Analysis of competitors showed that they can take wider driving lines with more speed. In order to improve this performance, the need for a high fidelity dynamic car model arose. This model can be used to understand the car dynamics, which can be used to choose a better suspension setup.<br/><br/>The baseline model used for this study was developed in MSC.Adams. After a discussion with field experts, two areas of model improvement were identified. In the baseline model, the tire is modeled using a MFTire model, where the vertical dynamics was modeled as a single spring and damper. In order to simulate the tire dynamics on kerb stones, a more advanced model is necessary. For this research, FTire was chosen as the best option. The other area of improvement is modeling the compliance of the suspension system. In the baseline model,<br/>all suspension members are modeled as rigid links, without any compliance. To model this compliance, an equivalent stiffness model was proposed, where all stiffnesses are combined into flexible joints.<br< /><br/>To use the FTire model, a set of parameters needs to be found for the Formula one tires. To find these parameters, measurements were done to capture the dynamics of the tire. A method is developed to filter out any rig disturbance, and make the data useful for parameter estimation. The estimation was done using Matlab, with MSC.Adams and FTire intheloop, and made use of a Genetic Algorithm (GA). Comparing the simulated behavior with the measurement showed good correlation for various conditions.<br/><br/>Since every part in the car contains some compliance, modeling and measuring every single part would be an extensive task. In order to reduce the complexity, all compliances are combined into a simplified model. In this simplified model, all compliances are modeled in the joints. A PatternSearch algorithm in Matlab is used in combination with MSC.Adams to find the joint stiffnesses in such way the model matches measurements of the overall car compliance. <br/><br/>To measure the overall improvement, both models are combined into a full car model. The baseline and proposed model are validated against measured track data. From this compare, it can be concluded that the proposed model shows a significant better correlation on kerb simulations. This is mainly due to the fact that the baseline tire model is not capable of simulating the harsh road input. The proposed model provides a more powerful tool for engineers to better understand the car dynamics, and to find the best suspension setup for<br/>kerb riding.LFormula one; Automotive; Motorsport; Compliance; Tire Dynamics; Optimization
20230216)uuid:49fa597132ae4252845dfcdd14412e53Dhttp://resolver.tudelft.nl/uuid:49fa597132ae4252845dfcdd14412e53Improving Rock Fragmentation at the Ahafo Gold Mine Site: Enhancing Improvements in Production Efficiency and Reducing Operational Costs Through Drill & Blast Pattern Optimizationcvan der Waal, Kaj (TU Delft Civil Engineering and Geosciences; TU Delft Geoscience and Engineering)UBuxton, Mike (mentor); Keersemaker, Marco (graduation committee); Rinne, M. (graduation committee); Lottermoser, B.G. (graduation committee); Coenen, J.H. (graduation committee); Delft University of Technology (degree granting institution); Aalto University (degree granting institution); RWTH Aachen University (degree granting institution)An efficient production is essential for an economic continuation of most gold mining companies, due to a moderate gold price and a competitive market. Recently, the Ahafo South Mining operations owned by the Newmont Mining Corporation and located on a Proterozoic gold belt in the south west of Ghana, indicated that the cost of production must be reduced in order to viably continue the mining operations. For this research, the Ahafo South Mine Site was used as a case study. A site visit of 90 days provided insight in where the efficiency of productivity could possibly be improved. During this visit, historical production and financial data [20142016] were gathered from a diverse group of engineers on site. Several audits were conducted to form a solid hypothesis with regards to the optimization of rock fragmentation. The hypothesis of this research states that improved rock fragmentation will result in a reduction of operational costs and delays along with an increase in crusher and mill throughput.<br/><br/>This case study was used to investigate the effect of a varying fragmentation on shovel loading time and crusher and mill performance. The financial impact of a potential blast pattern change was also investigated. After a statistical analysis of the available data on postblast fragmentation, it was observed that the desired run of mine (ROM) P80 target of 350 mm is rarely met. It was concluded that the "KuzRam" model used to predict the fragment size distribution is incorrectly applied, which resulted in a planned blast design that produces a coarse fragmentation with an excessive presence of boulders. This is referred to as "bad fragmentation" and proved to contribute significantly to redundant wea< randtear of mining equipment, delays at the shovel loading operations and the primary crushing station. A correlation between fragmentation and loading time was found after a statistical analysis on the daily average loading time of three Liebherr 9400 shovels. It was concluded that the shovels are only able to produce at the desired rate of 2.4 minutes per truck load when the fragmentation of the ROM meets the P80 target. The Six Sigma principle Define, Measure, Analyze, Improve, Control (DMAIC) was applied to define an accepted interval range of shovel loading time variation (2.4 m. 20 s.). This so called target interval was generally not achieved by the shovels. All three shovels show an increase in average loading time during period where the ROM fragmentation was coarser.<br/><br/>Data generated by a VisioRock"! camera showed, in combination with the throughput of the mills, that an appropriate fragmentation is crucial for an efficient milling process. A detailed analysis of the available data, categorized in "periods of high interest", demonstrated that the power consumption of the mills significantly increased when the feed consisted of relatively coarse fragments. It was observed that the postcrushed fragmentation was more coarse when a high fraction of the feed originated directly from the pit. This observation proved that the primary crusher is not able to eliminate all boulder related issues prior to the crushing operations as a result of inadequately planned blasts. Due to inaccurate or unreliable data on mill feed, it was concluded that a full mine to mill survey is required to obtain data needed for an accurate optimization of mill performance.<br/><br/>Newmont Ghana Gold Limited (NGGL) stated that at least 2.3M US$ has been saved on a yearly basis since 2015, as a result of the reduction of explosives used through a blast pattern widening. The cost impact analysis showed that it is uneconomical to save costs on explosives by widening the blast pattern, due to the potential downstream effects when the fragment size distribution is affected. The financial data was inconsistent and it is therefore recommended to conduct a detailed analysis to provide insight in the actual overall (financial) benefits of a blast pattern optimization. It is recommended to investigate both "doublebenching" and the purchase of a continuous surveying camera system. Another recommendation is to narrow the blast pattern down to 3.5 m. by 4 m. so that fragmentation related issues are avoided which are currently to be expected for the continuing production in the Subika and Awonsu Pit.<brRock Fragmentation; Gold Mine; Blast Pattern Selection; KuzRam; Shovel Loading Time; Mill Throughput; Powder Factor; ROM; Mining; Processing; Six Sigma; Optimization
20221222KEuropean Mining, Minerals and Environmental Program: European Mining Course)uuid:ee6dbfd2f5d2436686f6d4dcea609d7bDhttp://resolver.tudelft.nl/uuid:ee6dbfd2f5d2436686f6d4dcea609d7bUOptimization of the capacity of a rose sorting system using discrete event simulationwVreugdenhil, Wouter (TU Delft Mechanical, Maritime and Materials Engineering; TU Delft Marine and Transport Technology)iPang, Yusong (mentor); Bac, Wouter (mentor); Delft University of Technology (degree granting institution)In Rose cultivation companies in the Netherlands, there is a demand for a higher sorting capacity on the existing sorting systems. The objective for this research is to advice which part of the sorting system needs to be adjusted to gain a higher sorting capacity. For the current sorting systems, five bottlenecks are defined, the bottlenecks limit the sorting capacity. To be able to forecast the effect of machine adjustments, a discrete event simulation model has been constructed, using Simulink and Matlab. This simulation model is verified, matched and validated using data of on an existing rose sorting system. Results of a single day validation showed that the time to process the roses can be simulated with a 97% accuracy. Subsequently, five different simulations are executed. In each simu< lation, one of the five bottlenecks is removed or reduced. With the results of these simulations the capacity limitation due to each bottleneck is quantified. However entirely removing a bottleneck is not feasible in reality for all bottlenecks. A last situation is simulated where all feasible bottleneck reductions are combined. This showed that the time to sort all roses is reduced by 35%.=Discrete Event Simulation; Rose; Sorting system; Optimization
20200905)uuid:2355ab39550c44478c5661ffabd53353Dhttp://resolver.tudelft.nl/uuid:2355ab39550c44478c5661ffabd53353XMission Planner for HeatingOptimal ReEntry Trajectories with Extended Range CapabilityPapp, Z.Mooij, E. (mentor)
Atmospheric entry is a process defined by extremes, from the great velocity at the point of entry to the high deceleration loads and aerodynamic heating in the lower layers of the atmosphere. In addition to the vehicle maintaining structural integrity, an equally important aspect of the reentry mission is safety. This safety can only be guaranteed if the vehicle is controllable throughout the reentry process. Mission design plays a role in adhering to operational and safety requirements in the form of trajectory planning: the development of a reentry trajectory wherein the vehicle remains within its operational constraints and is controllable throughout. The focus of this thesis lies on the designtime trajectory development part of mission design, i.e., the development and analysis of feasible trajectories under certain requirements to be flown for a specific mission. The research goal of this thesis is formulated as follows: to what extent can optimal reentry trajectories be developed in the designtime phase of mission development for a winged entry vehicle that provide a maximumrange capability under the objective of minimizing heat loads and adhering to operational constraints? To answer this question, a mission design tool is developed in four successive stages: 1) Development of the reentry simulator, 2) Design of the guidance algorithm, 3) Development of the mission planner, and 4) Mission planner testing. The capacity to quickly and reliably simulate reentry trajectories is paramount to a mission planner. For this purpose, a simulator was developed with the specific goal of later integration with the mission planner. Steering is achieved by modulating the vehicle s attitude in terms of its angle of attack and bank angle. The guidance profile is based on the specification of attitude commands at specific points in the trajectory related to the instantaneous energy of the vehicle. The course of an entire trajectory can be specified by its guidance profile. The purpose of the mission planner is to develop guidance profiles for trajectories that keep the vehicle within its operational constraints, minimize the heat load, and provide the largest possible range under these conditions. The objectives of minimum heat load and maximum range are conflicting. Trajectories with minimum heat load requirements are generally short in duration with smooth heat flux profiles where the heat flux is maintained close to its constraint value. The total heat load is minimized by keeping the duration of the reentry as short as possible. This, however, is in direct opposition to the objective of maximizing range, where keeping the vehicle aloft for as long as possible is beneficial. The mission planner develops trajectories by specifying the individual parameters in this guidance profile by performing multiobjective optimization to determine the combinations of parameters that result in optimal trajectories in terms of the mission objectives. The mission planner consistently provides a set of optimal trajectories over a diverse range of objective values and under the provision that operational constraints are met. Optimal trajectories are determined based on these conflicting objectives over a range of objective values, wherein the relative priorities of the objectives are varied. In all cases is the mission planner able to provide trajectories with an exte< nded range capability, even when minimizing the heat load is considered the main priority.@Reentry; Guidance; Optimization; Mission planning; Trajectories Astrodynamics and Space Missions)uuid:40d4af9516b14b3f9182ac96929e2525Dhttp://resolver.tudelft.nl/uuid:40d4af9516b14b3f9182ac96929e2525A change of state: A thermodynamic and costeffective optimized Trombe wall based on latent heat storage (LHS) for year round applicationDHendriks, Kees Jan (TU Delft Architecture and the Built Environment)Tenpierik, Martin (mentor); Klein, Tillmann (graduation committee); Turrin, Michela (graduation committee); Delft University of Technology (degree granting institution) The following graduation thesis "A change of state: A thermodynamic and costeffective optimized Trombe wall based on latent heat storage (LHS) for year round application aims to investigate and optimize the thermal energy performance of a Phase change material (PCM) Trombe wall to create an economically feasible product. This passive system reduces the total cost of ownership of the energy system, inside an office building located in The Netherlands, by reducing the energy demand and the maximum peakloads from the mechanical system. In the European Union (EU) buildings are responsible for almost 40% of the total energy consumption, the energy efficiency of these buildings needs to increase by at least 32,5% by 2030 according to the Sustainable development goals (SDGs) from the EU. Research has shown that the information on the economic feasibility and optimization of PCM within the built environment is somewhat limited, this research will give an indepth insight in the actual performance on thermodynamics and the costeffectiveness. The PCM Trombe wall will be optimized by means of the Research through design  method, using different design strategies based on the knowledge of the thermodynamic principles of the PCM. An initial simulation model is used for the energy performance calculations, this model is extended and developed according to the parameters defined for the optimization. A combination of MATLAB and Simulink, a simulation environment based on textual and graphical programming, together with modeFRONTIER is employed. The results from these simulations are first validated with DesignBuilder Software Ltd to verify the legitimacy of the results from the added components and the changed location input data. The actual performance of the PCM on the reduction of the heating and cooling energy demand, the maximum powerload and the investment cost for the product are determined according to dynamic setpoint calculations. A detailed study together with an yearly performance analysis is used, the results from all the different simulations and optimizations are summarized in a design guideline. This guideline gives a clear indication on the differences between the input parameters and the results from the performance of the system. In the end, these results are all brought together in an adaptive and integrated design solution for the application in an office building in the Netherlands.Phase change materials; Optimization; Costeffectiveness; Thermodynamics; Total Cost of Ownership; Economic feasibility; PCM,Architecture, Urbanism and Building Sciences)uuid:43fec58c7ec34995be9e1224ded49995Dhttp://resolver.tudelft.nl/uuid:43fec58c7ec34995be9e1224ded49995pDevelopment of three dimensional PV structures as shading devices for a Decentralized Facade Unit of the Future.Sampatakos, D.CBergsma, A. (mentor); Martin, C.L. (mentor); Van Manen, W. (mentor)The first part of the Thesis consists of a research on topics related to Solar Energy, PV technologies and Materials used as part of active solar structures (from the active materials to substrates) as well as their integration in the form of custom shapes. The analysis of these topics reveals the potential of new and emerging PV technologies, which combine very lightweight and flexible characteristics, with a very good performance in partial shading and diffuse light conditions, but most importantly a v< ery high price reduction potential due to specific fabrication methods and materials. In a next step a (four phase) Methodology is developed in order to observe how the process of shaping a three dimensional PV structure in order to increase its surface area per facade space can lead to performance advantages compared to conventional flat and rigid PV structures. This part consists of a first phase, where possible light control/shading systems are compared according to specific requirements related to shading efficiency but more importantly to PV integration potential and PV surface increase potential. The second phase is a shape research in terms of Surface Area Potential and a comparison between 3d shapes using parametric tools. Thereby the ability of different shapes with various folding or curvature patterns to increase their surface for predefined total dimensions is evaluated. The results lead to structures which can provide much more active PV area for a given facade area. The third phase is a solar analysis of 56 of the three dimensional shapes tested and compared in phase 2. Thereby the real average solar radiation reaching the surface of various prismatic or corrugated 3d patterns is calculated and graphically presented. Incident Solar Energy values and partial shading percentages within every shape clearly illustrate which 3d shapes have the highest potential for making maximum use of solar radiation. The relation between surface area increase through denser folding or curving patterns and total received solar energy constitutes the most important finding of this phase. In a fourth phase Laboratory tests are performed at the PVLAB of the Faculty of Electrical Engineering TUDelft. Different 3d shape versions are tested using flexible Thin Film solar panels. The fact that solar cells performance (output) does not decrease linearly proportional to partial shading and total incident solar energy changes makes real Lab Tests an import ant verification of the solar analysis results. Finally the best performing shapes are used as the basis to create design proposals in the form of a light control/shading device. The PV optimized 3d surface is thereby developed as a facade structure with possible ways of operation and integration into the NEXT FACADE Concept by Alcoa Architectuursystemen. Additionally ways of using the produced energy and creating a complete efficient PV based system are proposed.2Facade; Solar Cells; PV; 3d; Shading; Optimization
20140930Internation Facade Master)uuid:6ce5a7a8ccba422d99ffe5159773ba5fDhttp://resolver.tudelft.nl/uuid:6ce5a7a8ccba422d99ffe5159773ba5fJFlexibility in Future Power Systems with High Renewable Energy Penetration\Hernandez Serna, Ricardo (TU Delft Electrical Engineering, Mathematics and Computer Science)Bauer, Pavol (mentor); van der Blij, Nils (mentor); Ramirez Elizondo, Laura (graduation committee); Delft University of Technology (degree granting institution)Future energy systems will present high penetrations of variable renewable energy (VRES). Due to the uncertainty that is related to the intermittency of wind and solar energy, and loads such as electric vehicles, the need for flexibility will be higher. This research identifies the possible sources of flexibility in a power system. Firstly, using statistical errors and probability theory, the required reserves of flexibility that a system needs to ensure a proper operation are quantified. Secondly, with the increase of VRES, the amount of synchronous generators is expected to decrease. This will mean that the inertial response used to keep a stable frequency will be reduced. On this paper, a method to quantify inertia and the concept of minimum required reaction time of flexibility are explained. Finally, an optimization problem is proposed to allocate the different types of flexibility taking into consideration their technical capabilities. This problem minimizes the cost of reserving different flexibility sources in a distribution grid.:Flexibility; Power Systems; Optimization; Renewable Energy$Master sustainable< energy technology)uuid:6b16e0f9b0384d53bc3cb07d7fbfe957Dhttp://resolver.tudelft.nl/uuid:6b16e0f9b0384d53bc3cb07d7fbfe957dDesign of an Optimized Acoustic Array for Aeroacoustic Research in an Openjet Anechoic Wind Tunnel7Luesutthiviboon, Salil (TU Delft Aerospace Engineering)Snellen, Mirjam (mentor); Sijtsma, Pieter (mentor); Malgoezar, Anwar (mentor); Delft University of Technology (degree granting institution)This work aims to determine the optimal microphone placement on an acoustic array of TU Delft s Vtunnel which is used for beamforming in aeroacoustic studies. The beamforming performance is driven by two parameters; the Maximum Side lobe Level (MSL) and the Main Lobe Width (MLW). The array design should give a good tradeoff between these parameters. The proposed optimization method has two optimization loops. First, the main loop consists of design variables used to collectively describe the distribution of microphones. Then the nested loop generates arrays which satisfy the geometry descriptions from the main loop. Finally, the main loop searches for the optimal design variables. The optimized array is able to achieve the MSL below 15 dB up to the distance approximately four times the MLW around the main lobe. Experimental validation was also carried out to compare the optimized array s performance with a benchmark array and an array from the beginning of the optimization.$Acoustic array; Optimization; Design)uuid:24ce69c5999442ef89ce26114555e970Dhttp://resolver.tudelft.nl/uuid:24ce69c5999442ef89ce26114555e970The Hybrid Power Plant: Optimal design and operation of battery energy storage as an addition to onshore wind farms in the NetherlandsJSijtsma, Lennard (TU Delft Mechanical, Maritime and Materials Engineering)de Jong, Wiebren (mentor); Heijnen, Petra (graduation committee); Ramaekers, Lou (graduation committee); Delft University of Technology (degree granting institution),
In the Netherlands, wind power has the highest potential for future development compared to alternatives such as solar PV and biomass. The intermittent nature of the wind power resource calls for novel approaches to supply and demand management, as well as power quality assurance. Energy storage technologies allow for the separation between power generation and power supply to the grid. The energy and power densities, efficiency and response times of secondary batteries make them highly suitable for utility scale application to wind farms. However, a knowledge gap exists on exact configurations of battery energy storage systems for wind farms within the existing power markets. This study proposes a hybrid power plant approach, combining an onshore wind farm with a battery energy storage system. Literature research identifies balance between power obligation and production as a main design criterium, along with the profitability of a configuration and operational strategy. A technical analysis establishes power fluctuation on multiple time scales as relevant wind power characteristics, and maximum depth of discharge, discharge rates and efficiency as important battery parameters. A detailed analysis of the Dutch power market structure identifies three accessible markets for trade: the Day Ahead Market, the primary reserve (FCR) and the secondary reserve (aFRR). For modelling and optimization purposes, FCR is not considered due to excessive uncertainties regarding bid acceptance and bid activation. The design phase of the study is charged with finding the optimal battery capacity for a given wind farm production, while assuring the obligation to the Day Ahead Market is met. A linear optimization model is developed with the objective to maximize the obtained revenue. It is shown that the optimal battery capacity scales with the power capacity of the wind farm. For a 60 MW wind farm, the optimal battery capacity lies within the range of 60 80 MWh. An operational model is designed based on the model of the plant established in the design phase, elaborated with a Model Predictive Control approach. Based on wind power generati< on forecasts, the algorithm is able to adjust the charge/discharge and imbalance settlement strategy continuously based on expected market price series. However, the low level of accuracy of the applied market price approximations lead to unprofitable results in all simulated cases. Market price approximations with increased accuracy, as well as reduction of battery installation costs may lead to profitable operation of a hybrid power plant.:Wind Energy; Batteries; Optimization; Modelling; Wind Farm)uuid:278180cf809d4ab2864bd1841982c35fDhttp://resolver.tudelft.nl/uuid:278180cf809d4ab2864bd1841982c35fParameter estimation for lateral dynamic tire models using empirical inertial data: A comparative case study of gradientbased algorithms vs genetic algorithmsLden Hartog, Daniel (TU Delft Mechanical, Maritime and Materials Engineering)pHappee, Riender (mentor); Shyrokau, Barys (mentor); Delft University of Technology (degree granting institution)5Selfdriving vehicles are the future of automotive engineering. Systems that take over control from the driver are developed to be able to interact with the conditions of the road and other obstacles. To develop these systems, developers use vehicle models to simulate the behaviour of the moving vehicle. The systems developed using these models are transferred the vehicle and are expected to perform accordingly. Therefore it is important that the vehicle models reflect the behaviour of the actual vehicle.<br/><br/>This thesis investigates the methods used to estimate the model parameters which influence the resulting vehicle simulation. Experimental tests were carried out to acquire vehicle body inertial measurement data. This data was used to minimize lateral acceleration and yawrate error produced by the output of a single track vehicle model using a simplified Magic formula developed by H. Pacejka.<br/>Deterministic gradientbased algorithms, such as multiobjective Sequential Quadratic Programming (SQP), are frequently used as a way to optimize an initial approximation of the model parameter set. It is shown that this initial parameter strongly influences results of the algorithm due to local minima.<br/><br/>As an alternative, genetic algorithms were investigated to minimize the influence of the initial parameter set. An adaptive component, Temporal Difference Q Learning (TDQL), was also added to further reduce the input of the user and to increase the performance of the algorithm.<br/><br/>Four algorithms, of various complexity, were implemented in Matlab and executed. The performance of the resulting parameter sets, as well as the performance of the algorithms themselves, were analysed and compared. It is concluded in this thesis that the adaptive genetic algorithm performs slightly better than a simple gradientbased algorithm with respect to the objectives and the duration of the algorithm. However, a genetic algorithm without TDQL is recommended for its good performance and the simulation flexibility of the results for this simple vehicle and tire model.oVehicle dynamics; Machine Learning; Parameter Estimation; Optimization; Tire Dynamics; Magic Formula; Algorithm)uuid:34ce9384935241bc99ce2a54bd1f3361Dhttp://resolver.tudelft.nl/uuid:34ce9384935241bc99ce2a54bd1f3361;Informative Path Planning for Search and Rescue using a UAVKAnil Meera, Ajith (TU Delft Mechanical, Maritime and Materials Engineering)QSiegwart, Roland (mentor); Wisse, Martijn (mentor); Popovi, Marija (mentor); Millane, Alexander (mentor); Pan, Wei (graduation committee); Alonso Mora, Javier (graduation committee); Mohajerin Esfahani, Peyman (graduation committee); Delft University of Technology (degree granting institution); ETH Zrich (degree granting institution)Target search in an obstacle filled environment is a practically relevant challenge in robotics that has a huge impact in the society. The wide range of applications include searching for victims in a search and rescue operation, detecting weeds in precision agriculture, patrolling borders for military and navy, automated census of endan< gered species in a forest etc. An efficient target search algorithm provides a data acquisition platform with least human intervention, thus improving the quality of life of humans. This thesis aims at introducing a general path planning algorithm for UAVs flying at different heights in an obstacle filled environment, searching for targets in the ground field. An adaptive informative path planning (IPP) algorithm is introduced that simultaneously trade off between area coverage, field of view, height dependent sensor performance and obstacle avoidance. It plans under uncertainties in the sensor measurements at varying heights, and is robust against wrong target detections. It generates an optimal fixed horizon plan in the form of a 3D minimumsnap trajectory that maximizes the information gain in minimum flight time by providing maximum area coverage, without any collision with the obstacles. The resulting planner is modular in terms of the mapping strategy, environment complexity, different target, changes in the sensor model and optimizer used. The planner is tested against varying environmental complexities, demonstrating its capability in handling a wide range of possible environments. The planner outperforms other planners like nonadaptive IPP planner, coverage planner and random sampling planner, by demonstrating the fastest decrease in map error while flying for a fixed time budget. A proof of concept for the algorithm is provided through real experiments by running the algorithm on a UAV flying inside a lab environment, searching for targets lying on the ground. All the targets were successfully found and mapped by the algorithm, demonstrating its applicability in a reallife target search problem.OPath Planning; Robotics; Unmanned Aerial Vehicle; Computer Vision; Optimization$Mechanical Engineering  Biorobotics)uuid:4e54cfb33024420a820548c7eb767ce9Dhttp://resolver.tudelft.nl/uuid:4e54cfb33024420a820548c7eb767ce9XParametric Modeling and Optimization of Advanced Propellers for NextGeneration Aircraft^Klein, Peter (TU Delft Aerospace Engineering; TU Delft Aerodynamics, Wind Energy & Propulsion)qSinnige, Tomas (mentor); Vitale, Salvatore (mentor); Delft University of Technology (degree granting institution)As traditional fossil fuels become scarcer and more attention is given to environmental impact of the combustion of fossil fuels, for nextgeneration aircraft, the focus of development will mainly be on reducing fuel consumption. Open rotor engines have the advantage over conventional turbofans that they are expected to perform 25% to 30% better in terms of fuel consumption. The topic of this thesis focuses on one of the current challenges related to open rotor configurations; the unwanted unsteady effects imposed on the aerodynamic characteristics of a propeller, imposed by an upstream pylon. There is a potential for reducing undesired propellerpylon installation effects by taking these effects into account in aerodynamic design process of the propeller. A propeller design optimization routine is implemented that includes a parametric modeling tool, a low and highfidelity performance analysis method and an existing installation effects performance model.gPropeller; CFD; Installed pusher configuration; Optimization; Parametric design tool; Noise; Pylon wake)uuid:71ec065fcec04b67a8cd90882c3a86a7Dhttp://resolver.tudelft.nl/uuid:71ec065fcec04b67a8cd90882c3a86a7Evaluating cooperation policies for rail utilization in the port to hinterland freight transport system: A combined method approachBKarampelas, Dimitrios (TU Delft Civil Engineering and Geosciences)Tavasszy, Lri (mentor); Maknoon, Yousef (mentor); Duinkerken, Mark (mentor); Kourounioti, Ioanna (mentor); Delft University of Technology (degree granting institution)BAs the margins for improvements in the current freight transport system become limited, researchers address more and more the importance of collaboration between the actors which is crucial for the implementation of new, more efficient transport concepts as synchromodality. In additio< n, rail is concerned as a sustainable mode of transport that can also achieve economies of scale due to its ability to haul large quantities of goods. This study investigates cooperation policies that affect actors behavior to better utilize the rail use and lead to a more efficient system. We propose an innovative approach that combines gaming, simulation and optimization as a mixed method to test and evaluate these policies. The port to hinterland freight transportation system in the range of Port of Rotterdam is used as a case study. First, gaming sessions are organized in order to observe actors behavior and collect data. The game that is used was initiated by Port of Rotterdam, especially to identify the problems in this system. Subsequently, by assessing the observed data, a simulation model is developed and different policy scenarios are simulated to quantify their performance. In addition, the optimization model is developed, which sets the upper bound for performance and used as a solid base for comparison between the policy alternatives. Finally, the explanation of the difference between the policies and the optimized performance can give an insight on what are the root causes of the inefficiency, what is the best allocation of the resources and where the solutions should be focused.kcooperation; Policies; Collaboration; Freight Transport; Port; hinterland; gaming; simulation; Optimization'Transport, Infrastructure and Logistics)uuid:053099df8008436ea51ccca6e5598861Dhttp://resolver.tudelft.nl/uuid:053099df8008436ea51ccca6e5598861GPerformance optimization of wind turbine rotors for wind farm operation Koole, R.Van Bussel, G.J.W. (mentor) Due to fierce competition on the electricity market, many wind energy related research projects are currently aimed at cost reduction. In order to establish cost reduction, more and more wind turbines are being clustered in (offshore) wind farms. Major drawback of wind farms however, is the energy production loss due to aerodynamic wake interaction between wind turbines. Optimizing wind turbine rotor designs on wind farm scale instead of focussing on wind turbine level, can increase the energy production of a wind farm. In the current project, optimization problems have been investigated that manipulate wind turbine thrust and power curves. The thrust and power curves were manipulated by Matlab's optimization toolbox, using builtin optimization algorithms. The optimization algorithm aims to reduce wake losses, which results in the maximum annual energy production of a wind farm. ECN's wind farmwake simulation tool FarmFlowfast was used to calculate the annual energy production of a wind farm, for each evaluated set of optimization parameter values. The current project involved two main objectives. First, the FarmFlowfast model was validated for optimization purposes. This was done by comparing simulation data with power measurement data from the Horns Rev and Lillgrund wind farms. Secondly, the potential of wind turbine performance curve optimization was investigated. This was done by investigating multiple scenarios, in which the increase in annual energy production of a wind farm was calculated as function of the wind turbine thrust and power curves. It is found that FarmFlowfast is capable to be used for optimization purposes. The most suitable optimization algorithm, concerning time consumption and objective function value decrement, was the pattern search algorithm. Unique optimum performance curves were found, depending on the wind farm layout, accounted wind speed and direction ranges and the parametrization and optimization method. Although the potential increase in annual energy production was rather small, the thrust reduction on wind turbine rotors around rated wind speed might as well be an important outcome of the optimization runs. The thrust reduction can result in cheaper wind turbine designs, which might even be a more promising result in terms of cost of energy reduction than the increase in annual energy production. The results of this project contri< bute to choosing new design parameters for future wind turbine rotor designs, in order to minimize the costs of energy.1Windenergy; Wind farm; Wind turbine; OptimizationAerodynamics & Wind EnergyWind Energy52.78028, 4.67763)uuid:6e27409599204d209e11d5b76363e709Dhttp://resolver.tudelft.nl/uuid:6e27409599204d209e11d5b76363e709WFuel Cell and Battery Hybrid System Optimization: Towards Increased Range and Endurance2Hoogendoorn, Joey (TU Delft Aerospace Engineering)Voskuijl, Mark (mentor); de Boer, Alte (mentor); van Benthem, Roel (mentor); Veldhuis, Leo (graduation committee); Roling, Paul (graduation committee); Delft University of Technology (degree granting institution)yAircraft manufacturers are increasingly exploring emissionfree flight or emission reduction for larger passenger aircraft. The low energy density of stateoftheart battery technology limits the application to small, electric, fixed wing aircraft up to a flight time of approximately one hour. To overcome these limits, a combination of fuel cells and batteries to exploit the benefits of battery power density and hydrogen energy density was studied. Current LithiumIon battery cells reach approximately 1.6 kW/kg of maximum power density, much higher than fuel cell systems. On the other hand, the energy storage capacity of suitable hydrogen storage methods is much larger than battery cells, the latter have an energy density of 240 Wh/kg.<br/><br/>Because most demonstrated applications are for fixed wing aircraft, the unmanned GeoCopter GC201 helicopter was used for performance requirements, weight and volume analysis. The study focuses on the preliminary sizing of the powertrain and the optimization of fuel cell and mission profile variables for this vehicle. Helicopter performance modelling, fuel cell static behavior as well as a battery discharge simulation are combined with lower fidelity models for other components. <br/><br/>The study results in a comparison of batteryonly and fuel cellbattery configurations through payloadrange diagrams, allowing for a quick evaluation of application areas. These mainly show that batteries excel at high payload, low range applications whereas a fuel cellbattery combination shows clear advantages at low payload, longer range applications. Liquid hydrogen will be shown to be comparable to the current micro gas turbine powered rotorcraft, with 400 and 500 km range capabilities respectively. Range capabilities for 300 bar and 700 bar compressed gas tank storage options show 140 and 180 km, with batteryonly reaching a maximum range of 80 km.Fuel Cell; Batteries; PEMFC; Helicopter; UAV; LithiumIon; Preliminary Sizing; Optimization; Endurance; Range; Hybrid; Unmanned; Rotorcraft)uuid:068ff836099b4abe8d9ecf96706169dfDhttp://resolver.tudelft.nl/uuid:068ff836099b4abe8d9ecf96706169df_Classifying Candida species using Mixed Integer Optimization based optimal classification treesRvan Dijk, Mick (TU Delft Electrical Engineering, Mathematics and Computer Science)van Iersel, Leo (mentor); Stougie, Prof. dr. L. (mentor); Kelk, Steven (mentor); Boekhout, Prof. dr. T. (graduation committee); Aardal, Karen (graduation committee); Delft University of Technology (degree granting institution)Global medical use of azole antifungals and echinocandins has led to an enormous increase in resistant <i>Candida</i> species, that are most commonly associated with fungal infections. A possible mechanism causing resistance are single or simultaneous point mutations in the genes responsible for encoding antifungal target enzymes. The aim of this thesis is to apply and compare several classification algorithms, in particular decision tree algorithms, on <i>Candida</i> data sets received from the Westerdijk Fungal Biodiversity Institute. Bertsimas and Dunn recently introduced a novel formulation based on Mixed Integer Optimization to generate optimal classification trees. We have implemented this method and applied it on <i>C. albicans </i>and <i>C. glabrata</i> data sets to construct univariate and multivariate classification tree< s. We were able to correctly classify 6872% of the <i>C. albicans</i> isolates and 76.582.5% of <i>C. glabrata</i> isolates. Moreover, by changing the objective function and adding constraints to the original MIO formulation, we constructed trees that take into consideration false negative errors, decreasing this type of error by 6480% for <i>C. albicans</i> and 5666% for <i>C. glabrata</i>. To deal with ambiguous nucleotides in the <i>C. albicans</i> data set we introduced a novel formulation to construct nonbinary classification trees. It turned out that ternary trees are a good representation of the <i>C. albicans</i> data set, performing strong in terms of outofsample accuracy. Finally, we identified combinations of amino acid substitutions and nucleotide mutations possibly related to resistance in <i>C. albicans</i> and <i>C. glabrata</i>..Optimization; Machine Learning; Bioinformatics)uuid:b980aab7f346403097a3a3cde13a51d6Dhttp://resolver.tudelft.nl/uuid:b980aab7f346403097a3a3cde13a51d6Optimization, Scaling and Economics of Pumping Kite Power Systems: Smallscale Systems for Developing Countries or Utilityscale Units for the European Market? Grete, C.1van Bussel, G.J.W. (mentor); Schmehl, R. (mentor)UOptimization; Scaling; Economics; Developing Countries; Kite Power System; Windenergy4Aerodynamics, Wind Energy & Propulsion / Wind Energy KitePower)uuid:1197caf1d9ad4829bdb88c3867e606e9Dhttp://resolver.tudelft.nl/uuid:1197caf1d9ad4829bdb88c3867e606e9TModelling Approach for Special Properties of Magnetic Materials Applied in ActuatorsNZhang, Jie (TU Delft Electrical Engineering, Mathematics and Computer Science)Polinder, Henk (mentor); Bax, Peter (mentor); Dong, Jianning (graduation committee); Lahaye, Domenico (graduation committee); Delft University of Technology (degree granting institution)9Eddy currents are loops of electrical currents induced within conductors by varying magnetic field in the conductor. Eddy currents flow in the closed loops in the conductor, in planes perpendicular to the magnetic field. The eddy current effect plays an important role in the transformer or an actuator system. <br/>If there is a coil wound around a cylindrical iron core, system, the AC coil can induce eddy currents in the solid core. The eddy current can change the flux in the core, and further influence the impedance of the system. If such a system is applied to an electromagnetic actuator, a proper modeling approach for the eddy current effect can help the researchers to determine the impedance change and predict the output power of the actuator. This thesis deals with the eddy current modeling improvement problem. A concept of the transformer ring topology is proposed to simplify the actuator model. A brass ring and an iron material ring are studied, in the purpose of realizing modeling optimization and test functions respectively. By using analytical calculation and modeling analysis, the eddy current modeling method is improved and optimized. The practical impedance of the model becomes closer to the simulation results. This research can help the researcher to compute the actuator model, and better simulates it.%Eddy current; Modelling; OptimizationElectrical Sustainable Energy)uuid:ac7c2d810fca470bbab7e47e0bf15fa9Dhttp://resolver.tudelft.nl/uuid:ac7c2d810fca470bbab7e47e0bf15fa9Modelbased Individualization of Humanlike Steering Controllers: Parameter Estimation of a Cybernetic Steering Model Using Global Optimization TechniqueXin, G./Abbink, D.A. (mentor); De Schutter, B. (mentor)
Current autonomous steering controllers are designed based on either optimal control strategies or on average human driver steering behavior. Individual drivers may wish to steer differently than those controllers do. Such conflicts can be mitigated by customizing individual steering behavior. This thesis aims to develop a method to individualize humanlike steering controllers, using a parameter estimation technique for a humanlike steering model based on a global optimization algorithm. With this technique, the < accuracy of estimated parameters is investigated in order to help to understand the nature of driver steering behavior. A linear humanlike cybernetic steering model inspired by human physical steering control actions is used in this thesis. In this steering model, seven parameters relate vehicles states to steering activity, of which the parameter values can be interpreted in a physically meaningful way. Preliminary research has demonstrated that this model structure enables parameter estimation of model parameters. But the current implemented parameter estimation technique of the steering model still needs requirements on initial guesses of the parameters and the accuracy of estimated parameters is underinvestigated. In this study, parameter estimation of the steering model is investigated by formulating an optimization problem to minimize the prediction error. To deal with the nonlinearity of the optimization problem, a genetic algorithm is complemented with the LevenbergMarquardt algorithm. This thesis presents the parameter estimation process in two steps. In the first step, the steering model generates a data set by simulating it with known parameters. This data set is used to estimate the steering model parameters and the results are compared with the known parameter set. Accurate parameter estimation results could not be reached in this first step, which is shown to be caused by overparameterization of the steering model. Two parameters in the neuromuscular system of the steering model are then simplified to mitigate the overparameterization problem and the other five parameters are remained. These five parameters can be estimated accurately by examining the metrics of VAF and Euclidean distance after the model simplification. In the second step, a simulation data set from a different validated humanvehicle model is used for parameter estimation to demonstrate the algorithm. It is shown that with the simplified steering model, the proposed parameter estimation technique can estimate consistent parameter values with good VAFs. A crossvalidation test is also implemented to give the corroborative evidence to indicate overparameterization of the steering model."Parameter estimation; OptimizationControl Engineering)uuid:b47ea25951214ffbab5dddee2e81a68bDhttp://resolver.tudelft.nl/uuid:b47ea25951214ffbab5dddee2e81a68b.Optimization of indepth divergence strategies<Sangers, Marvin (TU Delft Civil Engineering and Geosciences)Leeuwenburgh, Olwijn (mentor); Voskov, Denis (mentor); Zitha, Pacelli (graduation committee); Heemink, Arnold (graduation committee); Delft University of Technology (degree granting institution)@This study discusses the potential of the indepth water diversion (IDD) process to enhance oil recovery from subsurface reservoirs and investigates IDD strategies through an ensemble optimization approach. Pilot studies in the North Sea performed in recent years have shown that Sodium Silicate can be applied as blocking agent, diverting water to unswept zones of the reservoir. On the other hand, numerical simulation studies have tried to simulate the IDD behaviour with methods based on weak coupling of a reservoir flow simulator and an external chemical module. This study presents a fully implicit coupled chemicalcompositionalflow implementation to simulate the permeability reduction through silicate gels, since IDD is essentially a coupled flowchemical process. In addition, this study discusses the impact of resolution in space and time on the simulation performance for 2D subsurface petroleum reservoir models. Sensitivity of the oil recovery, and IDD characteristics such as moles of silicate, to design parameters of the IDD process is discussed as well. Since adjoint gradients are not typically available for the parameters describing indepth divergence and the uncertainties are expected to be large, the optimization study uses an ensemblebased methodology to find optimal IDD reservoir management strategies. In addition to discussing IDD strategies in a deterministic setting, the design of optimal< IDD strategies under geological uncertainty is investigated. This study will demonstrate that the indepth divergence process can be used to extend the reservoir field production life time when timing, size of the Sodium Silicate batch and concentration is optimized. Finally this study discusses the issue of computational costs associated with modeling high resolution required for accurate simulation of this coupled process.TOptimization; Indepth divergence; Sodium silicate; Petroleum; ensemble optimization)uuid:c3fac2a7dda443dbbe6e0b0fdee69949Dhttp://resolver.tudelft.nl/uuid:c3fac2a7dda443dbbe6e0b0fdee69949Pheromone networks  a computational approach to improve the social dynamics of the Almere central train station and it's urban network.La Roi, F.M.@Biloria, N.M. (mentor); Bier, H.H. (mentor); Sobota, M. (mentor)The amount of travelers from and to the Dutch city Almere will grow from 40.000 back in 2009 to 100.00 in 2030. This requires intensive infrastructural improvements to the current situation of Almere. One of the main improvements of the network will be the redesign of the central station of the city. The current station is situated around a combination of small dysfunctional public spaces and conflicting flows of traffic. In a series of experiments with the topics of different network speeds, users, behavior and needs is tested to propose a infrastructural landscape and embodiment that is both efficient and pleasant as a new design for the central train station of Almere. The quest is no search for an optimal result, but rather a proposal abstract machine that can generate pretty good solutions while still maintaining intuitive decision making in the design process. The aim of the design is to enhance the connectivity of the dysfunctional part of the Almere city center and to create a multifunctional area with a 18 hour usage pattern. The design method is based on a bottomup research process with the use of different algorithmic principles for generation and performance evaluation like swarm behavior and anttrail and Dijkstra algorithm s. Organizing infrastructure is the main interest in this project, not only to make a more efficient network but also to revitalize the area by influencing the behavior patterns of the users of the area.architecture; parametric; pheromone; ant trail; algorithm; swarm; Infrastructure; Almere; hyperbody; selforganization; Dijkstra; Optimization; CNC; Bridge; Train station; computational
20120704Architecture HyperbodyhHyperbody Graduation Studio: Advanced NonStandard & Interactive Architecture Formations & Embodiments)uuid:18a163a298a1438fafc9a66c500009edDhttp://resolver.tudelft.nl/uuid:18a163a298a1438fafc9a66c500009edImprovement of road traffic sustainability by implementation of priority weights for trucks in predictive signalized intersection controlIHaanstra, Leroy (TU Delft Mechanical, Maritime and Materials Engineering)Jiang, Xiaoli (mentor); Polinder, Henk (graduation committee); Taale, Henk (graduation committee); Delft University of Technology (degree granting institution)o
In the European Union Road freight transport volume is expected to grow 78% between 2000 and 2030, which results in more trucks on the road network. The worldwide estimated trend shows an increase of 150 million freight vehicles and an increase of 240 million passenger vehicles. The growth of both vehicle classes will have a major impact on the road network and the roads will become congested. Especially in dense urban environments with many intersections. Further, trucks have a detrimental impact on traffic flows, especially at intersections, because of their slow dynamics and large size. In addition, a stopping truck results in higher emissions and fuel consumption compared to a car. However, today s traffic controllers are not capable of optimizing traffic flow at intersections based on classification of different vehicles. Therefore, it would be beneficial to all vehicles involved if the number of stops for trucks would be reduced to a minimum, by servicing each vehicle c< lass in a different way.<br/>Throughout this research the focus is to develop a signalized intersection controller which is able to reduce the number of stops for trucks, while maintaining an equal level of service for the other modes of transport. Extensive literature studies provided important insights into the development of signalized intersection controllers. A selection of the techniques found in the literature is used to develop a new truck signal priority controller design. The basic idea is to use a weighted traffic light schedule in combination with priority weights to enable truck signal priority. The design is evaluated in a case study and simulated for multiple configurations. This leads to an overview of the performance in terms of the number of stops and vehicle delay. <br/>From the results several conclusions can be drawn on the optimal weight configuration, which is compared to a stateoftheart model predictive intersection controller. The results showed an reduction of the total number of stops by 751 and total vehicle delay by 304 minutes for trucks over a week (26,6% and 20,4% respectively). While, the total number of stops and total vehicle delay for cars increased, by 155 stops and 201 minutes (0,1% and 0,3% respectively). However, the overall total number of stops and total vehicle delay were reduced by 596 stops and 103 minutes (0,42% and 0,16% respectively). It can be concluded that the proposed truck signal priority controller design can reduce the number of stops for trucks at a signalized intersection, while maintaining the overall traffic flow at least as good as a stateoftheart model predictive intersection controller.SSignalized intersection controller; Truck Signal Priority; Simulation; Optimization
20190531)uuid:00967ce63a9c4f0eaa22974e1f065a8aDhttp://resolver.tudelft.nl/uuid:00967ce63a9c4f0eaa22974e1f065a8aGOptimization of Automotive Structure using Variable Stiffness Laminates'Francis, Sachin (Aerospace Engineering)van Campen, Julien (mentor); Abdalla, Mostafa (mentor); Peeters, Danil (mentor); Delft University of Technology (degree granting institution)oVariable Stiffness Laminates are created by spatially varying the fiber orientation resulting in designs that makes use of curved fibers rather than Unidirectional fibers to tailor the properties to design requirement. The advent of automated manufacturing methods such as Automated Fiber Placement and Tailored Fiber Placement has made variable stiffness laminate more realistic and attractive for use in aerospace and automotive sector. Importance of light weight designs in automotive sector has received new interest due to the stringent emissions rules and entry of designs based on alternative energy. Within this context, this thesis intends to create variable stiffness design for an automotive part to investigate the possible improvements in structural responses compared to a design based on conventional laminates. The design is done based on 2D Finite element analysis coupled with a convex optimization that helps to generate steered fiber designs optimized for strength, stiffness and buckling. One of the important load case studied here is the Inertia Relief. Inertia Relief method is a computationally efficient way of analyzing rigid bodies without doing a dynamic analysis. The method was implemented in the finite element framework and the responses from the analysis were taken for optimization.<br/>Final optimization of the structure showed that significant improvement in the objective responses can be achieved by using variable stiffness laminates over a conventional laminate design (UD). The result also implies that reduction in weight can be achieved if a variable thickness optimization is to be done on the model.<Optimization; VariableStiffness; Composites; Inertia Relief)uuid:ba778a4714c145d0870a24b698096e3aDhttp://resolver.tudelft.nl/uuid:ba778a4714c145d0870a24b698096e3a3FSI Analysis and Optimization of 3D beam structures7Mozo Gutierrez, Alvaro (TU Delft Aerospace Engineering)zTurteltaub, Sergio (mentor); Maut< e, K (graduation committee); Delft University of Technology (degree granting institution) Aiming to enhance the performance of the industrial design process, structural optimization techniques have been proposed as an alternative to traditional design and optimization techniques.<br/>They own the potential of achieving an optimal distribution of the material through the design domain, thus reducing waste loss and weight, and increasing the ability to carry the loads and the overall efficiency of the design.<br/>During this thesis work it has been developed a 3D FluidStructure Interaction model for beamlike structures, such as wind turbine blades. The crosssection geometry of the beam can be optimized using structural optimization techniques, such as size, shape or topology optimization.<br/>The proposed 3D beam is partially based on the formulation of the classical beam element for slender beams (EulerBernoulli), including SaintVenant torsional effects for isotropic materials, and with the addition of the terms related with the coupling between axial and torsion, and bending and torsion contributions, which may arise when using nonlinear materials.<br/>The stiffness information of the beam is interpolated from its crosssection geometries and materials, which can vary along the beam length.<br/>The crosssection geometries are defined on a XFEM mesh. The fluid and solid domains are specified using a Level Set Function. This provides a smooth geometry and crisp representation of the solid/fluid interface without the necessity of remeshing, as in the case of classical FEM. A 2D fluid simulation based on Incompressible Navier Stokes flow at low Reynolds number is carried around each crosssection, in order to obtain the aerodynamic loading over its contour. This aerodynamic loading serves as an input for the beam model, to compute deformation of the beam. This deformation is mapped onto the crosssections, obtaining the updated displacements and rotations of the geometry. With the updated geometry the fluid field is altered and it needs to be updated as well, forming a nonlinear iterative process that loops until a converged structure is obtained.<br/>The 3D FSI model is solved on a monolithic NewtonRaphson solver that treats all the equations involved at once. The Jacobian terms derived for the monolithic solving scheme that has been developed for the forward analysis allow a straightforward computation of the sensitivities using adjoint method. This sensitivity analysis makes possible the optimization of the geometry of the crosssections based on certain criteria and constraints.H3D Beam; FSI; XFEM; Level Set Method; Sensitivity Analysis; Optimization.Aerospace Structures and Materials Engineering)uuid:8b76b21a4f434b52b468d641dcdce02fDhttp://resolver.tudelft.nl/uuid:8b76b21a4f434b52b468d641dcdce02fNOptimization of the fluidization system of the trailing suction hopper dredgerEKamal, Emal (TU Delft Mechanical, Maritime and Materials Engineering)van Rhee, Cees (mentor); Talmon, Arnold (mentor); Alvarez Grima, M (mentor); Delft University of Technology (degree granting institution)kThe Trailing Suction Hopper Dredgers (TSHD) are typically used for the transportation of soil. The soil in the hopper can be removed by hydraulic transport using water jets and the discharge at the bottom level of the hopper. The soil fluidization and transportation is jet assisted. The unloading process can take quite some time; the unloading process can be limited by the soil transportation instead of the pumping system. The scope of this research is the optimization of the jetting system in order to deliver high mixture densities, to increase the unloading production and therefore decrease the unloading time.:Optimization; Fluidization; Trailing suction hopper dredge)uuid:0a67e3e9358a494fa553b65ed9f3c8c6Dhttp://resolver.tudelft.nl/uuid:0a67e3e9358a494fa553b65ed9f3c8c6*Path Optimization of a Pumping Kite SystemgJerez Venegas, Marcos (TU Delft Aerospace Engineering; TU Delft Aerodynamics, Wind Energy & Propulsion)Schmehl, Rola< nd (mentor); Gaunaa, M (mentor); Delft University of Technology (degree granting institution); Technical University of Denmark (degree granting institution)Airborne wind energy is a new research field that aims to harvest winds at high altitudes. The pumping kite concept investigated at the Delft University of Technology is based on leading edge inflatable kites that divide a generation cycle into three phases: traction, retraction, and transition. This research project aims to optimise the trajectory flown by the kite to increase the costcompetitiveness of the system. Thus, the complete system was simulated, and a sensitivity analysis was conducted. The results show that there is an area in the wind window that maximises the power production of the cycle. This area is defined by the wind speed, the shape of the figure, and the reel in/out force. Based on the findings, an algorithm that seeks the optimal path during the flight was developed. In the case study, the algorithm led to a 100% increase in cycle power and a 16% increase in the average traction power.=Airborne; Optimization; Path; Pumping Cycle; Kite; Windenergy4Aerospace Engineering  Aerodynamics and Wind EnergyEuropean Wind Energy Master)uuid:2987341fd34f4477b31f0cac08cd9c5cDhttp://resolver.tudelft.nl/uuid:2987341fd34f4477b31f0cac08cd9c5c:Optimization of the Locations of Virtual Network FunctionsMInan, Alp (TU Delft Electrical Engineering, Mathematics and Computer Science)tKuipers, Fernando (mentor); Klepper, Matthijs (mentor); Delft University of Technology (degree granting institution)In traditional networks, the network functions are implemented as physical appliances which are called middleboxes and these middleboxes need to be deployed in predetermined locations throughout the network. Therefore one of the most important disadvantages of using middleboxes is that each network function deployment requires a great deal of human work. During their operation, middleboxes need maintenance which also demands for massive human work. Therefore using middleboxes for each network function puts a constraint on the flexibility of the network. Network Function Virtualization (NFV) is a proposed architectural framework which introduces the virtualization of network functions. In NFV, each network function is implemented as software and deployed on data center servers which offers greater network flexibility with the ability of runtime network service deployment throughout the network. With this technology, the places of network functions can be determined for satisfying different requirements such as Quality of Service (QoS) for endusers. Determining optimal locations for Virtual Network Functions (VNFs) is an open research area. <br/>In this research, the optimal placement of VNFs is studied. The objective of this work is determined to minimize the endtoend Round Trip Time (RTT) between the end users and the deployed network functions while also minimizing the number of data centers used.<br/>First, the problem is formulated using Integer Linear Programming (ILP) and it is shown that it is NPhard. Therefore an updated problem definition is formulated to propose a heuristic. The proposed heuristic uses the knapsack problem notion to solve the problem in a smaller amount of time with a tradeoff in objective value. Lastly, for the proof of concept, the updated problem definition is also formulated using ILP and the proposed heuristic and the ILP formulation are compared under different scenarios.&NFV; Optimization; Heuristic; VNF; ILP)uuid:1e7dca7a575f49358265822643484344Dhttp://resolver.tudelft.nl/uuid:1e7dca7a575f49358265822643484344bBand gap grading of a SiO<sub>x</sub>:H and its effect on the performance of thin film solar cellsRodriguez Luna, Rogelio (TU Delft Electrical Engineering, Mathematics and Computer Science; TU Delft Electrical Sustainable Energy)Wvan Swaaij, Rene (mentor); Delft University of Technology (degree granting institution)Solar cells are showing significant promise to become the solution for the growing energy needs of our w< orld. However for this to happen, new disruptive technologies with high efficiency and low cost are needed in the market. One possibility comes from multijunction thin film solar cells based on aSi alloys and ncSi. For this purpose a SiO<sub>x</sub>:H is an interesting material since it can have V<sub>OC</sub> values above 1 V and FF above 0.7. However when paired in a tandem structure with a material of an advantageous bandgap, like ncSi:H. The performance of the tandem solar is limited in output current by the a SiO<sub>x</sub>:H layer.<br/><br/> <br/><br/>Research to increase J<sub>SC</sub> by increasing thickness of the absorber layer show it is not practical to increase the thickness of ia SiO<sub>x</sub>:H above 250 nm. Spillover knowledge from other thin film solar cells (GaAs, aSiGe<sub>x</sub> and CIGS), showed bandgap grading was able to increase performance of the electrical parameters. Grading in a solar cell means that in one of the layers a material property is varied continuously in concentration in order to achieve a different performance.<br/><br/>Our aim was to experiment with bandgap grading in the absorber layer of a SiO<sub>x</sub>:H solar cells, to try to achieve a higher J<sub>SC</sub> while still retaining the high V<sub>OC</sub> x FF product. Test layers were deposited at different CO<sub>2</sub>/SiH<sub>4 </sub>ratios to determine the dependence of the bandgap (E<sub>04</sub>) and the deposition rate on the CO<sub>2</sub>/SiH<sub>4</sub> ratio. With the experimental data and fitted polynomial equations a method was devised to vary continuously the bandgap in a step wise manner. Using this grading method, experiments were designed where the intrinsic a SiO<sub>x</sub>:H layer of a total length of 200 nm was subdivided in 3 graded bandgap regions. The first graded region named pi started from the end of the player with a high bandgap (2.1 eV). Decreasing the bandgap over a certain length until reaching a region with no added oxygen with a low bandgap (1.96 eV.) From here the central i region started, maintaining a constant bandgap for a certain length until the bandgap starts increasing again. This marks the 3rd region called in, where the bandgap continues to increase over a certain width until reaching 2.1 eV at the beginning of the nlayer. <br/><br/> <br/><br/>Graded experimental cells results showed that it was beneficial to have a small graded region width in the pi and in region (1030 nm). Since with this grading length J<sub>SC</sub> could be increased significantly (8% increase) from reference cell values without grading performed. However a small compromise in a drop of V<sub>OC</sub> and FF values around (12%) was observed at the same time from reference cell values. The gains in J<sub>SC</sub> are bigger than the loss of V<sub>OC</sub> x FF product and result in a relative efficiency gain of (45%) from reference cell. This technique paired with other cutting edge techniques to increase photocurrent can lead to a new record for an a SiO<sub>x</sub>:H thin film solar cell. Or it can lead to better current matching in tandem or triple junction solar cells.ZSolar cell; amorphous; Thin film; silicon oxide; bandgap; Grading; Profiling; Optimization)uuid:4a40e302bdf1431981dea6aad9376d65Dhttp://resolver.tudelft.nl/uuid:4a40e302bdf1431981dea6aad9376d65dMultipoint aerodynamic shape optimization for airfoils and wings at supersonic and subsonic regimes/Mangano, Marco (TU Delft Aerospace Engineering)la Rocca, Gianfranco (mentor); Martins, Joaquim (mentor); Veldhuis, Leo (graduation committee); Dwight, Richard (graduation committee); Delft University of Technology (degree granting institution)rThe secondgeneration of supersonic civil transport has to match ambitious targets in terms of noise reduction and efficiency to become economically and environmentally viable. Highfidelity numerical optimization offers a powerful approach to address the complex tradeoffs intrinsic to this novel configuration. Past and current research however, despite proving the potential of such design strateg< y, lacks in deeper insight on final layouts and optimization workflow challenges. Stemming from the necessity to quantify and exploit the potential of modern design tools applied to supersonic aircraft design, this work partially fills the gap in previous research by investigating RANSbased aerodynamic<br/>optimization for both supersonic, transonic and subsonic conditions. The investigation is carried out with the stateoftheart, gradientbased MDO framework \textit{MACH}, developed at University of Michigan's MDO Lab  which hosted the author for the 14month research stint. Details of the tool and a brief overview of supersonic aircraft design and modern aerodynamic optimization strategies are reported in the first part of this manuscript.<br/>After circumscribing the research niche, I perform single and multipoint optimization to minimize the drag over an ideal supersonic aircraft flight envelope and assess the influence of physical and numerical parameters on optimization accuracy and reliability. Leading and trailing edge morphing capabilities are introduced to improve the efficiency at transonic and subsonic flight speed by relaxing the tradeoffs on clean shape optimization. Benefits in terms of drag reduction are quantified and benchmarked with fixededges results. It is observed how the optimized airfoils outperform baseline reference shapes from a minimum of 4\% up to 86\% for different design cases and flight<br/>conditions. The study is then extended to the optimization of a planar, lowaspectratio, and lowsweep wing, using the same schematic approach of 2D analysis. I investigate the influence of wing twist alone and twist and shape on cruise performance, obtaining a drag reduction of 6\% and 25\% respectively as the optimizer copes with both viscosity and compressibility effects over the wing. Results for 3D multipoint optimization suggest that the proposed strategy enables a fast and effective design of highlyefficient wings, with drag reduction ranging from a minimum of 24\% up to 74\% for cruise at different speeds and altitudes, including edge deflection. Ultimately, this work provides an extensive and, to the best of author knowledge, unprecedented insight on the optimal design solutions for this specific aircraft configuration and the challenges of the optimization framework. The benefits of RANSbased aerodynamic shape optimization to capture nonintuitive design tradeoffs and offer deeper physical insight are ultimately discussed and quantified. Given the promising results in terms of performance improvements and design efficiency, it is hoped that this work will foster the implementation of this method for more comprehensive fullconfiguration, multidisciplinary supersonic aircraft optimization studies.mOptimization; Aerodynamics; MDO; CFD; Supersonic; wing design; Airfoil; morphing; Gradientbased Optimization)uuid:7d3bb76b295746e49a2685350f8975aeDhttp://resolver.tudelft.nl/uuid:7d3bb76b295746e49a2685350f8975aefFactors affecting container transshipment volumes at ports:: A data driven holistic modelling approachJVeldstra, Tijmen (TU Delft Mechanical, Maritime and Materials Engineering)Frouws, Koos (mentor); Hekkenberg, Robert (mentor); Beelaerts van Blokland, Wouter (mentor); Helbing, Jolke (mentor); Delft University of Technology (degree granting institution) In 20 years, from 1996 to 2016, global cumulative container port throughput had an average annual growth rate of 6.3%. Furthermore, the average global transshipment incidence increased from 11% in 1980 to a stable 27% in 2017. This change in transshipment market dynamics has led to opportunities for potential port investments or disinvestments. In this research, a tool is developed for Royal HaskoningDHV that improves transshipment forecasts. The developed tool calculates the optimal ports for transshipment, transshipment throughput volumes, fleet utilization, and optimal cost routes. This is done by means of minimizing the total transport costs, while meeting container demand requirements and port and vessel constraints. The most sta<tistically significant factors affecting transshipment volumes are incorporated in the tool. By modifying input parameters, the tool facilitates calculations of optimal cost future transshipment port throughput volumes. This could provide key insights in lucrative port investment or disinvestment opportunities.`Container shipping; Hub and Spoke; Transshipment; Mixed Integer Linear Programming; Optimization
20490424)uuid:311d536e542d45b597d9f52b2dd4abf3Dhttp://resolver.tudelft.nl/uuid:311d536e542d45b597d9f52b2dd4abf35Application of biobased FRP on a road traffic bridge]Zarifis, George (TU Delft Civil Engineering and Geosciences; TU Delft Engineering Structures)Veljkovic, Milan (mentor); Pavlovic, Marko (mentor); Hendriks, Max (mentor); Claassen, W. (mentor); Delft University of Technology (degree granting institution)The last decades Fiber Reinforced Polymer (FRP) structures start to find their own place in the field of civil engineering applications and quite lately in bridge engineering. Composites are considered highly versatile materials that offer great tailoring on their design. Synthetic fibers and resins are already quite popular especially in applications where low weight, construction speed and high performance are the objectives. However, environmental concerns over the disposability of nonrecycle composites is driving the engineering world towards the use of natural composites. Certain natural fibers appear to be advantageous not only in terms of mechanical properties but also in terms of sustainability compared to synthetic ones. However, environmental performance of bio materials is difficult to assess and predict due to their large natural variability. Durability assessment through literature review and optimization of the superstructure of an FRP bridge, accounting for failure modes mostly deriving from high concentrated vertical loading (wheel loading), are the main objectives of this report. In order to deal with the large number of variables and constraints for the design according to existing structural standards, advanced optimization algorithms are employed. In addition, due to lack of analytical formulas for local buckling resistance prediction, Finite Element (FE) analysis is used since it offers a better insight into buckling failure of composite laminates. Finally, an optimization workflow is developed that accounts for all variables and constraints and minimizes the objective, which in this case is the weight of the superstructure. The obtained results form recommendations for the optimum design of a flax and glass FRP bridge, regarding the thicknesses and orthotropy of the cross section s laminates with L/300 and L/500 maximum deflection requirement, providing a deck slenderness L/16 and a spacing of the web sw = hSP /3. The assessment and the optimization reveals that flax or BioMid FRP could be promising candidates in replacing synthetic fibers for a bridge application as they can compete in performance and weight the synthetic glass FRP for the case of single span bridge and a span range of 10  30m.$Biobased; FRP; Bridge; Optimization
*+&ffffff?'ffffff?(?)?"dXX333333?333333?U}}}}}}}}}} }
}}}
}}}}}}}}}}}}
@
!
"
#
$
%
&
'@
(
)
*
+
,

.@
/
0
1
2
3
4
5
6@
7
8
9
:
;
<
=
>@
?
@
A
B
C
D
E
F
G
H@
I
J
K
L
M
N
O@
P
Q
R
S
T
U
V@
W
X
Y
Z
[
\
]
^
_
` @
a
b
Y
c
d
e
f
g
h
@
i
j
Y
Z
k
l
m
n
o
p
q@
r
s
Y
t
u
v
w
x@
y
z
Y
{

}
~
@
Y
@
Y
@
Y
1
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
@
Y
!
!
!
!
!!@
!
!
!
! Y
"
"!
""
"#
"$"@
"%
"&
"
" Y
#'
#(
#)
#*
#+#@
#,
#
#
# Y
#.
#/
$0
$1
$2
$3
$4$@
$5
$6
$
$ Y
$7
%8
%9
%:
%;
%<%@
%=
%>
%
% Y
%?
&@
&A
&B
&C
&D&@
&E
&F
&
& Y
&G
'H
'I
'J
'K
'L'@
'M
'N
'
' Y
'
'O
(P
(Q
(R
(S
(T(@
(U
(V
(
( Y
)W
)X
)Y
)Z
)[)@
)\
)]
)
) Y
)^
)_
)`
*a
*b
*c
*d
*e*@
*f
*g
*
* Y
*h
*
+i
+j
+k
+l
+m+@
+n
+o
+
+ Y
+p
,q
,r
,s
,t
,u,@
,v
,
, Y
,w
x
y
z
{
p@
}
~

 Y



.
.
.
.
..l@
.
.
.
. Y
.
.
.
.
/
/
/
/
//@
/
/
/
/ Y
0
0
0
0
00@
0
0
0
0 Y
0
0
1
1
1
1
11@
1
1
1
1 Y
1
1
2
2
2
2
22@
2
2
2
2 Y
2
21
3
3
3
3
33@
3
3
3
3 Y
4
4
4
4
44@
4
4
4
4 Y
4
4
5
5
5
5
55@
5
5
5
5 Y
5
6
6
6
6
6m6@
6
6
6
6 Y
6p
7
7
7
7
77l@
7
7
7
7 Y
7
7
7
7
8
8
8
8
88@
8
8
8
8 Y
8
8
9
9
9
9
99@
9
9
9
9 Y
9
:
:
:
:
::h@
:
:
:
: Y
:Z
:
:
;
;
;
;
;;@
;
;
;
; Y
;
;
;
<
<
<
<
<<@
<
<
<
< Y
=
=
=
=
==@
=
=
=
= Y
=
>
>
>
>
>>@
>
>
> Y
>
>
?
?
?
?
??@
?
?
?
? Y
?h
@
@
@
@
@@@
@
@
@
@ Y
A
A
A
A
AA@
A
A
A
A Y
A!
A
B"
B#
B$
B%
B&B@
B'
B(
B
B Y
B)
B*
C+
C,
C
C.
C/C@
C0
C1
C
C Y
C
C2
D3
D4
D5
D6
D7D@
D8
D9
D
D Y
D:
E;
E<
E=
E>
E?E@
E@
EA
E
E Y
E
FB
FC
FD
FE
FFF@
FG
FH
F
F Y
Fp
GI
GJ
GK
GL
GMG@
GN
GO
G
G Y
GP
GQ
HR
HS
HT
HU
HVH@
HW
HX
H
H Y
IY
IZ
I[
I\
I]I@
I^
I_
I
I Y
I`
Ia
Jb
Jc
Jd
Je
JfJ@
Jg
Jh
J
J Y
Ki
Kj
Kk
Kl
KmK@
Kn
Ko
K
K Y
Kp
K
Lq
Lr
Ls
Lt
LVLp@
Lu
Lv
L
L Y
Lw
LZ
L[
Mx
My
Mz
M{
MM@
M}
M~
M
M Y
M
M
N
N
N
N
NN@
N
N
N
N Y
N
O
O
O
O
OO@
O
O
O
O Y
O
P
P
P
P
PP@
P
P
P
P Y
P
Q
Q
Q
Q
QQ@
Q
Q
Q
Q Y
Q
Q
R
R
R
R
RR@
R
R
R
R Y
R
R
S
S
S
S
SS@
S
S
S
S Y
S
T
T
T
T
TT@
T
T
T
T Y
T
U
U
U
U
UU@
U
U
U
U Y
U
U
V
V
V
V
VV@
V
V
V
V Y
V
V/
W
W
W
W
WWx@
W
W
W
W Y
W
W
X
X
X
X
XX@
X
X
X
X Y
X
Y
Y
Y
Y
YYx@
Y
Y
Y
Y Y
Y
YP
Y
Y
Z
Z
Z
Z
ZZ@
Z
Z
Z
Z Y
Z
[
[
[
[
[[@
[
[
[
[ Y[l{I@
\
\
\
\
\\@
\
\
\
\ Y
\_
]
]
]
]
]]@
]
]
]
] Y
^
^
^
^
^^@
^
^
^
^ Y
^
_
_
_
_
__@
_
_
_
_ Y
_
`
`
`
`
``@
`
`
`
` Y
`
a
a
a
a
aa@
a
a
a
a Y
a
a
a
a
b!
b"
b#
b$
b%b@
b&
b'
b
b Y
b
c(
c)
c*
c+
c,c@
c
c.
c
c Y
c1
d/
d0
d1
d2
d3dx@
d4
d
d Y
d
d5
dp
d6
e7
e8
e9
e:
e;e@
e<
e=
e
e Y
e>
f?
f@
fA
fB
fCf@
fD
fE
f
f Y
fZ
f
fF
gG
gH
gI
gJ
gKg@
gL
gM
g
g Y
hN
hO
hP
hQ
hRhp@
hS
hT
h
h Y
hU
hV
hW
hX
iY
iZ
i[
i\
i]i@
i^
i_
i
i Y
i`
i/
ja
jb
jc
jd
jej@
jf
jg
j
j Y
kh
ki
kj
kk
klk@
km
kn
k
k Y
ko
lp
lq
lr
ls
ltl@
lu
lv
l
l Y
l
mw
mx
my
mz
m{m@
m
m}
m
m Y
m~
m
n
n
n
n
nn@
n
n
n
n Y
n?
o
o
o
o
oo@
o
o
o
o Y
op
p
p
p
p
pp@
p
p
p
p Y
p
q
q
q
q
qq@
q
q
q
q Y
q
r
r
r
r
rr@
r
r
r
r Ystuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~>@ddrrryKyKhttp://resolver.tudelft.nl/uuid:56a2f8d7fb2640f0bb9200271f9492dfyKyKhttp://resolver.tudelft.nl/uuid:6818b21f038d4700afcde6dac79e7e82yKyKhttp://resolver.tudelft.nl/uuid:5ecc206d60bb43ef9874ee72c6263c6ayKyKhttp://resolver.tudelft.nl/uuid:0fc3b7460d3f459d8fc2c9e6aea821dfyKyKhttp://resolver.tudelft.nl/uuid:b7ccbf2752cf4d50ba649282aec248e4yKyKhttp://resolver.tudelft.nl/uuid:437ca7513fd642c3933f72771f025eeeyKyKhttp://resolver.tudelft.nl/uuid:6b8900808a60418989dcb12f23a8677ayKyKhttp://resolver.tudelft.nl/uuid:0126e698df974313a4d178f4c3f8c6d0 yKyKhttp://resolver.tudelft.nl/uuid:3e54884398e74a18a50e017309ae95d8
yKyKhttp://resolver.tudelft.nl/uuid:3f3e955b36bc4a09a9ccb7065a59d39byKyKhttp://resolver.tudelft.nl/uuid:f9a658fe9d0b4a1d808a908698c8e8c2yKyKhttp://resolver.tudelft.nl/uuid:87d1d976a0ba40d28fa4437bd65c5b91
yKyKhttp://resolver.tudelft.nl/uuid:7fe5a51181c944e590baa2b479011077yKyKhttp://resolver.tudelft.nl/uuid:78ae5715d0794a7e983cd999ba761c16yKyKhttp://resolver.tudelft.nl/uuid:056d28defe7e48ae8a9d0af5bd9aa2c0yKyKhttp://resolver.tudelft.nl/uuid:f3c8e60a65d245e1a838001a5c21ccebyKyKhttp://resolver.tudelft.nl/uuid:232d91048cd34ace89b3042e50490871yKyKhttp://resolver.tudelft.nl/uuid:841a9a8dc7d44bdba1a47a4e7960070byKyKhttp://resolver.tudelft.nl/uuid:00eb0f82521346578fa42630c3594463yKyKhttp://resolver.tudelft.nl/uuid:62636a6807ba424480362f4f9bc133a1yKyKhttp://resolver.tudelft.nl/uuid:0ad5f3cca5e945c0a37e75e139713409yKyKhttp://resolver.tudelft.nl/uuid:5137791874be4ec4bad2115276e35cd8yKyKhttp://resolver.tudelft.nl/uuid:c6924adf9f7b48508099056eb06bc11ayKyKhttp://resolver.tudelft.nl/uuid:538fd43d5ea841ca89927f40f8abc450yKyKhttp://resolver.tudelft.nl/uuid:b67eae7d4d5a47d2a320cf0e14944852yKyKhttp://resolver.tudelft.nl/uuid:58e639f7ea6a4e89baf80bfa64c9a7f4yKyKhttp://resolver.tudelft.nl/uuid:a9e6065b1560459f8bde992ddf4ede39yKyKhttp://resolver.tudelft.nl/uuid:ae23df6646bd4946ba9c1364996fde4byKyKhttp://resolver.tudelft.nl/uuid:3d712eb25d6f45268f2cbecc2eeaa13fyKyKhttp://resolver.tudelft.nl/uuid:c23bdb6d72db42f68241c44b475521f1yKyKhttp://resolver.tudelft.nl/uuid:38e13f1e84244a00bbad70b9eadf97ac yKyKhttp://resolver.tudelft.nl/uuid:d034d1e2fd984da28ddcc4f6a71acd99!!yKyKhttp://resolver.tudelft.nl/uuid:6f3aec1e9225431cb1a8be635907eb67""yKyKhttp://resolver.tudelft.nl/uuid:cc39b153a7994cceaabb80bf35af44c2##yKyKhttp://resolver.tudelft.nl/uuid:7879208441584ba591e0c88e6033a216$$yKyKhttp://resolver.tudelft.nl/uuid:bfe097c160604701b0d664c5e9bbb0f3%%yKyKhttp://resolver.tudelft.nl/uuid:72c525b07b48462f9df15864e4b08de3&&yKyKhttp://resolver.tudelft.nl/uuid:bb7dbde731ae4915b7335742e30b50ef''yKyKhttp://resolver.tudelft.nl/uuid:66aaf06787634c3991e3ef01fb25948d((yKyKhttp://resolver.tudelft.nl/uuid:588a8ded1f5942ca9eb594f2e1e996b8))yKyKhttp://resolver.tudelft.nl/uuid:aae2621eb81148debec139171151aae4**yKyKhttp://resolver.tudelft.nl/uuid:3378e0cf78314432a383c9c5c8dbce2d++yKyKhttp://resolver.tudelft.nl/uuid:13e189889b504cc8bf8e392455acfb44,,yKyKhttp://resolver.tudelft.nl/uuid:d8f71f9a6b594b519fe5b0a4358a0d7fyKyKhttp://resolver.tudelft.nl/uuid:64606bb7ce5c42ddaa257c1045ef3f80..yKyKhttp://resolver.tudelft.nl/uuid:d60603b662ba4be9a54fff371228e12a//yKyKhttp://resolver.tudelft.nl/uuid:581a19ab1bac4fbba678bc8e905804d700yKyKhttp://resolver.tudelft.nl/uuid:50aaee056a0b4000b685572d175d5dd711yKyKhttp://resolver.tudelft.nl/uuid:446ce1ef5c164ceba31beab72884b4d422yKyKhttp://resolver.tudelft.nl/uuid:5b0b6016b8d242d3a9ba6f18e9c63bca33yKyKhttp://resolver.tudelft.nl/uuid:611d466d3d264b3299e6e98a805435e844yKyKhttp://resolver.tudelft.nl/uuid:61a8031c9278494d88920fc4c8252b8455yKyKhttp://resolver.tudelft.nl/uuid:2c79c37252d8448993ae1e2f7c67346966yKyKhttp://resolver.tudelft.nl/uuid:6dfc6eb9676f4a0e99c58462ab70a3cb77yKyKhttp://resolver.tudelft.nl/uuid:a732d9ba71744b2081aa01f387e5de3588yKyKhttp://resolver.tudelft.nl/uuid:c5ddec13e10349b48b84a8ad013c753c99yKyKhttp://resolver.tudelft.nl/uuid:efd1c23f4ab741dd88e4e9a1683c4ccc::yKyKhttp://resolver.tudelft.nl/uuid:1c741d018a3f45b7abe9f32456791dc3;;yKyKhttp://resolver.tudelft.nl/uuid:28671d8aa00a4526a6700c1e80ad50ff<<yKyKhttp://resolver.tudelft.nl/uuid:f2e84a4aaef445b5a9d0e97afd17b04b==yKyKhttp://resolver.tudelft.nl/uuid:6184bd6263d24929ae947f549149f5a5>>yKyKhttp://resolver.tudelft.nl/uuid:36c9527eabd347bcafb86196e7d5e587??yKyKhttp://resolver.tudelft.nl/uuid:888bf27a4dd44b08ab9a4b7c8ef62ef7@@yKyKhttp://resolver.tudelft.nl/uuid:ccc91e95a7904fbcbcafeffedc9dbd4dAAyKyKhttp://resolver.tudelft.nl/uuid:d25e499e89de4fa5b96e3ada53141754BByKyKhttp://resolver.tudelft.nl/uuid:e095b2bec4974de4a55b25126e960dbeCCyKyKhttp://resolver.tudelft.nl/uuid:b7911d199b28445ebf1905db3b9f85ddDDyKyKhttp://resolver.tudelft.nl/uuid:2b6974777d594083b51d70c0d3137261EEyKyKhttp://resolver.tudelft.nl/uuid:04f50f3b85724c68b75be6990c639be0FFyKyKhttp://resolver.tudelft.nl/uuid:2694dcb5292f41e2a9f7e79b0dccc8c0GGyKyKhttp://resolver.tudelft.nl/uuid:8bd18483c1e1460b9f590f5e20726679HHyKyKhttp://resolver.tudelft.nl/uuid:dec3c264d6524370864c1f9ea18cbd2aIIyKyKhttp://resolver.tudelft.nl/uuid:c49d975200524db0af2015477951f6fdJJyKyKhttp://resolver.tudelft.nl/uuid:fa661ab603244115adc80c8665e2e1d9KKyKyKhttp://resolver.tudelft.nl/uuid:621c77c467644d179a0ab3d3dd3025e0LLyKyKhttp://resolver.tudelft.nl/uuid:86c056d8f3684239893f07ca3a22e112MMyKyKhttp://resolver.tudelft.nl/uuid:3e7cb93b7d3b4b33bcc0d562630f12e7NNyKyKhttp://resolver.tudelft.nl/uuid:cc24b8b9fca64387a746f7b50bb32c73OOyKyKhttp://resolver.tudelft.nl/uuid:044eec5fc3f94021a7d81f3303d3492fPPyKyKhttp://resolver.tudelft.nl/uuid:72d06941c9754d518d4ba976a7c0b9cfQQyKyKhttp://resolver.tudelft.nl/uuid:2cc94bdd857d48fe91405b7e8cfa001fRRyKyKhttp://resolver.tudelft.nl/uuid:ab92c1fba5b14302b2b4fcd36b84f0dcSSyKyKhttp://resolver.tudelft.nl/uuid:98fd5f3841d744b885e8c54e5aeb9c40TTyKyKhttp://resolver.tudelft.nl/uuid:2a15115f6d704057abdaeabf1adc72a7UUyKyKhttp://resolver.tudelft.nl/uuid:49fa597132ae4252845dfcdd14412e53VVyKyKhttp://resolver.tudelft.nl/uuid:ee6dbfd2f5d2436686f6d4dcea609d7bWWyKyKhttp://resolver.tudelft.nl/uuid:2355ab39550c44478c5661ffabd53353XXyKyKhttp://resolver.tudelft.nl/uuid:40d4af9516b14b3f9182ac96929e2525YYyKyKhttp://resolver.tudelft.nl/uuid:43fec58c7ec34995be9e1224ded49995ZZyKyKhttp://resolver.tudelft.nl/uuid:6ce5a7a8ccba422d99ffe5159773ba5f[[yKyKhttp://resolver.tudelft.nl/uuid:6b16e0f9b0384d53bc3cb07d7fbfe957\\yKyKhttp://resolver.tudelft.nl/uuid:24ce69c5999442ef89ce26114555e970]]yKyKhttp://resolver.tudelft.nl/uuid:278180cf809d4ab2864bd1841982c35f^^yKyKhttp://resolver.tudelft.nl/uuid:34ce9384935241bc99ce2a54bd1f3361__yKyKhttp://resolver.tudelft.nl/uuid:4e54cfb33024420a820548c7eb767ce9``yKyKhttp://resolver.tudelft.nl/uuid:71ec065fcec04b67a8cd90882c3a86a7aayKyKhttp://resolver.tudelft.nl/uuid:053099df8008436ea51ccca6e5598861bbyKyKhttp://resolver.tudelft.nl/uuid:6e27409599204d209e11d5b76363e709ccyKyKhttp://resolver.tudelft.nl/uuid:068ff836099b4abe8d9ecf96706169dfddyKyKhttp://resolver.tudelft.nl/uuid:b980aab7f346403097a3a3cde13a51d6eeyKyKhttp://resolver.tudelft.nl/uuid:1197caf1d9ad4829bdb88c3867e606e9ffyKyKhttp://resolver.tudelft.nl/uuid:ac7c2d810fca470bbab7e47e0bf15fa9ggyKyKhttp://resolver.tudelft.nl/uuid:b47ea25951214ffbab5dddee2e81a68bhhyKyKhttp://resolver.tudelft.nl/uuid:c3fac2a7dda443dbbe6e0b0fdee69949iiyKyKhttp://resolver.tudelft.nl/uuid:18a163a298a1438fafc9a66c500009edjjyKyKhttp://resolver.tudelft.nl/uuid:00967ce63a9c4f0eaa22974e1f065a8akkyKyKhttp://resolver.tudelft.nl/uuid:ba778a4714c145d0870a24b698096e3allyKyKhttp://resolver.tudelft.nl/uuid:8b76b21a4f434b52b468d641dcdce02fmmyKyKhttp://resolver.tudelft.nl/uuid:0a67e3e9358a494fa553b65ed9f3c8c6nnyKyKhttp://resolver.tudelft.nl/uuid:2987341fd34f4477b31f0cac08cd9c5cooyKyKhttp://resolver.tudelft.nl/uuid:1e7dca7a575f49358265822643484344ppyKyKhttp://resolver.tudelft.nl/uuid:4a40e302bdf1431981dea6aad9376d65qqyKyKhttp://resolver.tudelft.nl/uuid:7d3bb76b295746e49a2685350f8975aerryKyKhttp://resolver.tudelft.nl/uuid:311d536e542d45b597d9f52b2dd4abf3gg
Root Entry F?~5?~5@SummaryInformation( F<Workbook F2 DocumentSummaryInformation8 F
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~