The balancing of a wind turbine rotor is a crucial step affecting the machine’s performance, reliability, and safety, as it directly impacts the dynamic loads on the entire structure. A rotor can be balanced either statically or dynamically. A method of rotor balancing was developed that achieves both the simplicity of static balancing and the accuracy of dynamic balancing. This method is best suited, but not limited, to hollow composite blades of any size. The method starts by quantifying the mass and center of gravity of each blade. A dynamic calculation is performed to determine the theoretical shaking force on the rotor shaft at the design operating speed. This force is converted to a net counterbalance mass required for each blade. Despite the most careful methodology, there may still be large errors associated with these measurements and calculations. Therefore, this new method includes a physical verification of each blade’s individual balance against all other blades on the rotor, with the ability to quantify the discrepancy between blades, and make all balance adjustments in situ. The balance weights are aluminum plugs of varying lengths inserted into the root of each blade with a threaded steel rod running through the middle. The balance adjustment is thus not visible from outside. The weight of the plug and rod represent the coarse counterbalance of each blade, based on the dynamic calculations. The threaded steel rod acts as a fine adjustment on the blades’ mass moment when traveled along the plug. A dedicated blade-balance apparatus, designed and constructed in-house, is used to verify and fine-tune each individual blade and compare it to all other blades on the rotor. The resulting blade assembly is verified on a full rotor static balancing apparatus. The full rotor apparatus measures the steady state tilt of the rotor when balanced on a point. Next, the rotors' tilt is related to its overall level of imbalance with quantifiable error. Most error comes from the fact that the hub, comparable in mass to the blades, creates a false righting moment of the assembly not present in operation. The fully assembled rotor is tested, pre and post balance, in operation on the turbine at a series of predetermined speeds. This is accomplished with a 3-axis accelerometer mounted on the main turbine shaft bearing and a control system which regulates and records turbine speed at 100 Hz

Onshore wind turbine foundations are mainly subjected to large overturning moments. The wind action imposes cyclic and dynamic loading conditions which occur in extreme and service scenarios. Deep foundations, when used, transfer this large overturning moment through a pile group which combines the axial and lateral resistance of all piles. Several authors noticed that cyclic axially loaded piles could have their resistance reduced due to shaft friction degradation. Considerable efforts have been made to understand this degradation phenomenon. However, the design and performance of cyclic axial loaded piles require greater advances. Therefore, this research sought to assess the shaft friction degradation of axially loaded piles in wind turbine deep foundations, especially for those located in Brazil. Several issues related to the main objective of this study had to be addressed because onshore wind turbine foundations are an area of recent research, especially in Brazil. The first stage of this research explored the key aspects of onshore wind turbine foundations in Brazil and compared them with the worldwide status. The main reason to explore this subject is that several authors consider onshore wind turbine foundations a well-understood topic; however, limited data from actual situations have been published, especially in developing countries where wind energy projects have recently started. Thus, a survey on Brazilian energy companies and foundation designers was conducted, and the first Brazilian database of wind turbine foundations was created. This database contains data from more than three thousand Brazilian wind turbine foundations. The key aspects, types and dimensions of these foundations were summarized. Worldwide, concrete gravity foundations are the most commonly used foundation type for onshore wind turbines. In Brazil, 43.3% of the wind turbines had shallow foundations, essentially concrete gravity, and 56.7% had deep foundations, mostly continuous flight auger piles. The foundation type was chosen according to the local foundation expertise and geotechnical conditions, which included soil type, water table level, soil layer resistance, the extent of porous soil layers and bedrock depth. This first stage of the research identified that Brazilian wind turbine foundations are significantly different from other countries. Approximately 70% of Brazilian wind turbine deep foundations used continuous flight auger piles, most of them embedded in sandy soils. Therefore, experimental investigations of the sand-concrete interface response under monotonic and cyclic loading are essential. The second and third stages of this research sought to investigate the sand-concrete interface response based on this recent insight. The second stage assessed the sand-concrete interface response through monotonic interface direct shear tests under different confinement conditions. The role of surface structural characteristics, confinement condition, sand mean diameter, particle morphology, sand gradation and relative density were evaluated. A nonlinear conceptual model of the interfacial-to-internal friction angle ratio was proposed according to normalized roughness and normalized waviness. Additionally, multiple regression was used to estimate the sand-concrete interface strength by the effect of constant normal stiffness. The results were essential to understand and to predict the sand-concrete interface response of concrete piles under static axial loading. The third stage explored the shaft friction degradation of cyclic axial loaded piles through sand-concrete and sand-steel cyclic interface direct shear tests. In geotechnical engineering practice, field and experimental tests are usually performed to evaluate the number of cycles until failure occurs under constant cyclic amplitude. According to this approach, cyclic failure can either occur quickly or not at all. From a practical viewpoint, the cyclic test times are unpredictable, which makes these tests difficult to plan and to perform. Therefore, a new approach based on increasing cyclic amplitude is proposed to overcome the conventional method. A simple cumulative damage model established a relationship between the cyclic loaded tests under constant and increasing cyclic amplitudes. The new approach provides additional insights into the cyclic interfacial response, such as the effect of previous cycles, the cyclic amplitude at failure and the displacement development throughout cycling. This new approach can be effortlessly extended to other experimental and field investigations. From a practical viewpoint, this new approach can reduce the cost and duration of projects. The author believes that this dissertation brought breakthroughs to the wind energy companies and to the geotechnical engineering community. However, further studies on onshore wind turbine foundations are still required. / As fundações de torres eólicas onshore são submetidas a grandes momentos de tombamento. A ação do vento impõe carregamentos cíclicos e dinâmicos que ocorrem em condições operacionais e extremas. Fundações profundas, quando usadas, transferem esse momento de tombamento através de um grupo de estacas que, por sua vez, são submetidas a esforços axiais e laterais. Muitos estudos observaram que estacas submetidas ao carregamento cíclico axial estão sujeitas à degradação do atrito lateral. A fim de compreender essa degradação, diversas pesquisas foram realizadas. No entanto, o dimensionamento e desempenho de estacas submetidas ao carregamento cíclico axial requer maiores avanços. Dessa forma, esta pesquisa buscou avaliar a degradação do atrito lateral em estacas de torres eólicas, principalmente aquelas localizadas no Brasil. Da mesma forma, foram abordados assuntos relacionados ao objetivo principal deste estudo uma vez que pesquisas sobre fundações de torres eólicas onshore ainda são incipientes, especialmente no Brasil. A primeira etapa da pesquisa explorou as principais características das fundações de torres eólicas onshore no Brasil e comparou-as com a prática internacional. Investigou-se este assunto uma vez que diversos autores consideram as fundações de torres eólicas onshore um tópico já compreendido; entretanto, a quantidade de informações e dados publicados são ínfimos, especialmente em países em desenvolvimento onde os projetos eólicos iniciaram recentemente. Por essa razão, realizou-se uma pesquisa com as empresas e projetistas do setor a fim de criar o primeiro banco de dados brasileiro sobre fundações de torres eólicas onshore. Esse banco de dados possui mais de três mil fundações cadastradas. Os principais aspectos, tipos e dimensões dessas fundações foram apresentados. Internacionalmente, o tipo de fundação mais empregado para torres eólicas onshore são fundações superficiais de gravidade. No Brasil, 43.3% das torres eólicas apresentam fundações superficiais, essencialmente fundações de gravidade, e 56.7% das torres eólicas apresentam fundações profundas, principalmente por grupo de estacas hélice continua. O tipo da fundação foi determinado com base na expertise local e condições geotécnicas, incluindo o tipo de solo, nível de água, resistência do solo, espessura de camadas porosas e profundidade do topo rochoso. Identificou-se, através da primeira etapa, que as fundações de torres eólicas onshore no Brasil são significantemente diferentes de outros países. Aproximadamente 70% das fundações profundas são por grupo de estacas hélice contínua, sendo na maior parte em contato com solos arenosos. Desta forma, há necessidade de realizar investigações experimentais da interface areia-concreto. A segunda e terceira etapas desta pesquisa concentraram-se no comportamento estático e cíclico da interface areia-concreto. A segunda etapa avaliou o comportamento estático da interface areia-concreto através de ensaios de cisalhamento direto em diferentes condições de confinamento. Foram avaliadas as influências das características da superfície sólida, diâmetro médio dos grãos, morfologia dos grãos, distribuição granulométrica e densidade relativa. Um modelo não linear do ângulo de atrito na interface de acordo com a rugosidade e ondulação normalizadas foi proposto. Além disso, empregou-se uma regressão múltipla para estimar a resistência da interface areia-concreto de acordo com a constante de rigidez. Os resultados foram essenciais na compreensão e previsão do comportamento estático da interface de estacas de concreto. A terceira etapa explorou a degradação do atrito lateral em estacas através de ensaios de cisalhamento direto cíclicos na interface areia-concreto e areia-aço. Usualmente, campanhas experimentais são conduzidas para avaliar o número de ciclos até ruptura com uma amplitude cíclica constante. A ruptura geotécnica pode ocorrer rapidamente ou nunca ocorrer; e, desta maneira, a duração dos ensaios é imprevisível, dificultando planejamento e execução das obras. Por essa razão, uma nova abordagem foi proposta baseada em amplitudes cíclicas crescentes. Um modelo de dano acumulado estabeleceu a relação entre os ensaios com amplitude cíclica constante e cíclica. A nova abordagem fornece informações adicionais do comportamento cíclico da interface, como o efeito de ciclos anteriores, amplitude cíclica na ruptura e o acúmulo de deslocamento permanente. Essa nova abordagem pode ser facilmente aplicada em outras investigações experimentais e provas de cargas cíclicas. Do ponto de vista prática, essa abordagem pode reduzir o custo e tempo de projetos, além de melhorar a previsão do desempenho dessas fundações. Espera-se que essa pesquisa tenha trazido avanços para as empresas do setor e comunidade geotécnica. Contudo, novas pesquisas sobre fundações de torres eólicas onshore ainda são necessárias.

This thesis is on the space--time variational multiscale (ST-VMS) computation of wind-turbine rotor and tower aerodynamics. The rotor geometry is that of the NREL 5MW offshore baseline wind turbine. We compute with a given wind speed and a specified rotor speed. The computation is challenging because of the large Reynolds numbers and rotating turbulent flows, and computing the correct torque requires an accurate and meticulous numerical approach. The presence of the tower increases the computational challenge because of the fast, rotational relative motion between the rotor and tower. The ST-VMS method is the residual-based VMS version of the Deforming-Spatial-Domain/Stabilized ST (DSD/SST) method, and is also called ``DSD/SST-VMST'' method (i.e., the version with the VMS turbulence model). In calculating the stabilization parameters embedded in the method, we are using a new element length definition for the diffusion-dominated limit. The DSD/SST method, which was introduced as a general-purpose moving-mesh method for computation of flows with moving interfaces, requires a mesh update method. Mesh update typically consists of moving the mesh for as long as possible and remeshing as needed. In the computations reported here, NURBS basis functions are used for the temporal representation of the rotor motion, enabling us to represent the circular paths associated with that motion exactly and specify a constant angular velocity corresponding to the invariant speeds along those paths. In addition, temporal NURBS basis functions are used in representation of the motion and deformation of the volume meshes computed and also in remeshing. We name this ``ST/NURBS Mesh Update Method (STNMUM).'' The STNMUM increases computational efficiency in terms of computer time and storage, and computational flexibility in terms of being able to change the time-step size of the computation. We use layers of thin elements near the blade surfaces, which undergo rigid-body motion with the rotor. We compare the results from computations with and without tower, and we also compare using NURBS and linear finite element basis functions in temporal representation of the mesh motion.

Fresh water resources in many of the world's oil producing regions, such as western Texas, are scarce, while produced water from oil wells is plentiful, though unfit for most applications due to high salinity and other contamination. Disposing of this water is a great expense to oil producers. This research seeks to advance a technology developed to treat produced water by reverse osmosis and other means to render it suitable for agricultural or industrial use, while simultaneously reducing disposal costs. Pilot testing of the process thus far has demonstrated the technology's capability to produce good-quality water, but process optimization and control were yet to be fully addressed and are focuses of this work. Also, the use of renewable resources (wind and solar) are analyzed as potential power sources for the process, and an overview of reverse osmosis membrane fouling is presented. A computer model of the process was created using a dynamic simulator, Aspen Dynamics, to determine energy consumption of various process design alternatives, and to test control strategies. By preserving the mechanical energy of the concentrate stream of the reverse osmosis membrane, process energy requirements can be reduced several fold from that of the current configuration. Process control schemes utilizing basic feedback control methods with proportional-integral (PI) controllers are proposed, with the feasibility of the strategy for the most complex process design verified by successful dynamic simulation. A macro-driven spreadsheet was created to allow for quick and easy cost comparisons of renewable energy sources in a variety of locations. Using this tool, wind and solar costs were compared for cities in regions throughout Texas. The renewable energy resource showing the greatest potential was wind power, with the analysis showing that in windy regions such as the Texas Panhandle, wind-generated power costs are approximately equal to those generated with diesel fuel.

This thesis presents an aerodynamic shape optimization methodology for rotor blades of horizontal axis wind turbines. Genetic Algorithm and Blade Element Momentum Theory are implemented in order to find maximum power production at a specific wind speed, rotor speed and rotor diameter. The potential flow solver, XFOIL, provides viscous aerodynamic data of the airfoils. Optimization variables are selected as the sectional chord length, the sectional twist and the blade profiles at root, mid and tip regions of the blade. The blade sections are defined by the NACA four digit airfoil series or arbitrary airfoil profiles defined by a Bezier curve. Firstly, validation studies are performed with the airfoils and the wind turbines having experimental data. Then, optimization studies are performed on the existing wind turbines. Finally, design optimization applications are carried out for a 1 MWwind turbine.

Many Departments of Transportation (DOTs), including Texas Department of Transportation (TxDOT), have been challenged by inadequate funding from traditional federal and state fuel taxes, increasing construction cost, aging highway system, traffic congestions, and recent natural disasters, compromising their primary mission to provide safe vehicle transportation routes with adequate capacity. Furthermore, environmental awareness and sustainability concept have strengthened and sparked debates in Congress, culminating with several regulatory policies that affect, inclusively, transportation projects. This scenario has prompted DOTs to pursue innovative ways to reduce maintenance cost (at minimum) and generate revenue (at maximum) exploiting their assets, and to meet the new regulations. Likewise, the Center of Transportation Research at The University of Texas at Austin undertook a comprehensive research study to identify and determine when, where, and under what circumstances TxDOT should pursue the implementation of which Value Extraction Application (VEA), and how to effectively recognize and involve key stakeholders. As a result, 11 VEAs were identified. In addition, a methodological framework – embedding a multi-attribute criteria analysis matrix as the decision making method - was devised to guide TxDOT throughout the process of identifying, evaluating, comparing, and selecting the most appropriate VEA while a list of stakeholders associated with each VEA and a stakeholder analysis framework was provided to help TxDOT to identify and reach out key stakeholders. / text

Αντικείμενο της παρούσας διδακτορικής διατριβής είναι η πειραματική και υπολογιστική διερεύνηση αεροδυναμικής συμπεριφοράς πτερύγων σε διφασική ροή αέρα–νερού και η εφαρμογή σε πτερύγια ανεμοκινητήρων. Αρχικά, γίνεται πειραματική και υπολογιστική μελέτη μονοφασικής ροής αέρα γύρω από αεροτομές, πτέρυγες και πτερύγιο ανεμοκινητήρα και στη συνέχεια μελέτη διφασικής ροής αέρα-νερού γύρω από τα ίδια σώματα. Η σύγκριση μεταξύ των αποτελεσμάτων της μονοφασικής ροής με τα αντίστοιχα της διφασικής ροής αέρα-νερού είναι αναγκαία ώστε να μελετηθούν οι επιπτώσεις της διφασικής ροής αέρα–νερού στην αεροδυναμική απόδοση. Η πειραματική ανάλυση αφορά τη διεξαγωγή πειραμάτων για τη μελέτη της αεροδυναμικής συμπεριφοράς αεροτομών και πτερύγων σε συνθήκες μονοφασικής και διφασικής ροής. Για την προσομοίωση συνθηκών διφασικής ροής αέρα-νερού τροποποιείται η αεροσήραγγα που διαθέτει ήδη το Εργαστήριο με την προσαρμογή ειδικών ακροφυσίων ψεκασμού νερού (συνθήκες βροχής). Για τις ανάγκες των πειραμάτων χρησιμοποιούνται τα μοντέλα αεροτομών και πτερύγων NACA 0012 που συνοδεύουν την αεροσήραγγα και κατασκευάζονται αεροτομή και πτέρυγες S809. Τα πειράματα μονοφασικής και διφασικής ροής γίνονται για την ίδια ταχύτητα αέρα. Για τη διφασική ροή αέρα-νερού εξετάστηκαν τέσσερις διαφορετικές πυκνότητες περιεχόμενης βροχής. Η υπολογιστική ανάλυση γίνεται με το υπολογιστικό πακέτο ANSYS CFD-Fluent. Αρχικά, γίνονται προσομοιώσεις για μονοφασική ροή αέρα γύρω από την αεροτομή NACA 0012, για την οποία υπάρχει πλήθος δημοσιευμένων αποτελεσμάτων, με τρία διαφορετικά μοντέλα τύρβης ώστε να βρεθεί το καταλληλότερο. Ο συντελεστής άνωσης υπολογίζεται με μεγάλη ακρίβεια, σε αντίθεση με το συντελεστή αντίστασης. Το πρόβλημα αυτό οφείλεται στην αδυναμία του Fluent να υπολογίσει το σημείο μετάβασης του οριακού στρώματος από στρωτό σε τυρβώδες. Κρίνεται επομένως αναγκαίο να γίνει σύγκριση του συντελεστή αντίστασης με πειραματικά δεδομένα για πλήρως τυρβώδες οριακό στρώμα. Για ακόμα πιο ακριβή αποτελέσματα αναπτύσσεται αλγόριθμος για τον υπολογισμό του σημείου μετάβασης από στρωτό σε τυρβώδες οριακό στρώμα και γίνονται προσομοιώσεις ορίζοντας την περιοχή αριστερά από το σημείο μετάβασης ως στρωτή και δεξιά από αυτό ως τυρβώδη. Υπολογίζονται οι κατανομές πίεσης και ταχύτητας γύρω από την αεροτομή, καθώς επίσης και τα σημεία ανακοπής, μέγιστης ταχύτητας, αποκόλλησης και επανακόλλησης του οριακού στρώματος. Παρουσιάζονται επίσης οι ροϊκές γραμμές και τα διανύσματα της ταχύτητας γύρω από την αεροτομή. Αντίστοιχες προσομοιώσεις γίνονται και για την αεροτομή S809. Για τη μελέτη του τρισδιάστατου χαρακτήρα της ροής, γίνονται προσομοιώσεις γύρω από πτέρυγα S809. Υπολογίζονται οι συντελεστές άνωσης και αντίστασης, τα σημεία ανακοπής, μέγιστης ταχύτητας, αποκόλλησης και επανακόλλησης του οριακού στρώματος. Επίσης παρουσιάζονται κατανομές της έντασης της τύρβης στην άνω επιφάνεια της πτέρυγας και της συνισταμένης ταχύτητας, της ταχύτητας στη z-διεύθυνση, της έντασης της τύρβης και της επιτάχυνσης της ροής πίσω από την πτέρυγα. Για τη μελέτη της ροής γύρω από περιστρεφόμενο πτερύγιο γίνονται προσομοιώσεις γύρω από το πτερύγιο Phase IV της NREL. Γίνεται μελέτη της κατανομής της αξονικής ταχύτητας πίσω από το δρομέα, της κατανομής της στατικής πίεσης και της έντασης της τύρβης πάνω στην επιφάνεια του πτερυγίου και της κατανομής της στατικής πίεσης σε διάφορα σημεία πάνω στο πτερύγιο. Η υπολογιστική μελέτη της διφασικής ροής αέρα-νερού γίνεται αρχικά για την αεροτομή NACA 0012 με πυκνότητα περιεχόμενης βροχής LWC=30 g/m³, επειδή υπάρχουν αντίστοιχα έγκυρα πειραματικά αποτελέσματα ώστε να γίνει σύγκριση για την εγκυρότητα της διαδικασίας της προσομοίωσης. Στη συνέχεια γίνονται προσομοιώσεις για διφασική ροή αέρα-νερού γύρω από την αεροτομή S809, την πτέρυγα S809 και το περιστρεφόμενο πτερύγιο Phase IV της NREL. Προσομοιώσεις γίνονται επίσης για διαφορετικές πυκνότητες περιεχόμενης βροχής για τη ροή γύρω από τις αεροτομές σε χαμηλό αριθμό Reynolds. Τα αποτελέσματα της διφασικής ροής αέρα-νερού συγκρίνονται με τα αντίστοιχα της μονοφασικής ροής ώστε να προκύψουν συμπεράσματα για τις επιπτώσεις της βροχής στην αεροδυναμική απόδοση. Γίνεται επίσης υπολογισμός του συντελεστή ισχύος του ανεμοκινητήρα σε συνθήκες μονοφασικής ροής αέρα και διφασικής ροής αέρα-νερού. Σε συνθήκες διφασικής ροής αέρα-νερού παρατηρείται υποβάθμιση της αεροδυναμικής απόδοσης, συγκεκριμένα μείωση της άνωσης με παράλληλη αύξηση της αντίστασης. Δυο είναι οι βασικοί μηχανισμοί που επικρατούν και έχουν ως αποτέλεσμα την υποβάθμιση αυτή. Στην επιφάνεια της αεροτομής δημιουργείται ανομοιόμορφο φιλμ νερού που αυξάνει την τραχύτητα και το πάχος της αεροτομής. Τα σταγονίδια καθώς προσκρούουν πάνω στο φιλμ νερού δημιουργούν «κρατήρες» αυξάνοντας την τραχύτητα της αεροτομής. Επίσης, τα σωματίδια νερού διασπώνται κατά την πρόσκρουσή τους πάνω στην αεροτομή σε άλλα σταγονίδια μικρότερης διαμέτρου και μειωμένης ταχύτητας. Αυτό έχει ως αποτέλεσμα τα σταγονίδια αυτά, επαναεπιταχυνόμενα από τη ροή του αέρα να αποσπούν ποσό ενέργειας από το οριακό στρώμα καθιστώντας το πιο ευάλωτο σε αποκόλληση. Στόχος της μελέτης της αεροδυναμικής συμπεριφοράς των πτερυγίων σε διφασική ροή αέρα-νερού είναι η κατασκευή ανεμοκινητήρων υψηλού βαθμού απόδοσης και η παραγωγή φθηνής ενέργειας από την όσο το δυνατόν καλύτερη αξιοποίηση της αιολικής ενέργειας. / The aim of the present doctoral thesis is the experimental and computational study of the aerodynamic behavior of wings in two-phase flow and the application on wind turbine blades. First of all, experimental and computational study of one-phase flow over airfoils, wings and wind turbine blade and afterwards study of two-phase flow over the same bodies is conducted. The comparison of the results between dry and wet conditions is necessary in order to show the effects of two-phase flow at the aerodynamic performance. Wind tunnel tests were conducted to show the aerodynamic behavior of airfoils and wings in one-phase and two-phase flows. To simulate two-phase flow, the wind tunnel of the Fluid Mechanics Laboratory has to be configured with adding commercial rain simulated nozzles. For the experiments NACA 0012 airfoils and wings which come along the wind turbine are utilized and airfoil and wings S809 are constructed. The experiments of one-phase flow and two-phase flow are conducted for the same air velocity. For the two-phase flow four different Liquid Water Contents are examined. For the computational analysis the commercial CFD code ANSYS Fluent is used. In first place, simulations of one-phase flow over the NACA 0012 airfoil are done with three different turbulence models. The NACA0012 airfoil is chosen because it has been studied in depth and has a precise data base to compare the results of the simulation with. The lift coefficients are computed with accuracy in contrast to the drag coefficient. The overprediction of drag is expected since the actual airfoil has laminar flow over the forward half. The turbulence models cannot calculate the transition point from laminar to turbulent and consider that the boundary layer is turbulent throughout its length. Therefore, it is necessary to compare the computational results with experimental data of a fully turbulent boundary layer. In order to get more accurate results, the computational domain could be split into two different domains to run mixed laminar and turbulent flow. The contours of pressure and velocity over the airfoil are presented, as well as stagnation, maximum velocity, detachment and reattachment points of the boundary layer are computed. Streamlines and velocity vectors over the airfoil are also presented. Similar simulations are conducted for the S809 airfoil. In order to study the tree-dimensional effects of the flow, simulations over the S809 wing are made. Lift and drag coefficients, stagnation, maximum velocity, detachment and reattachment points of the boundary layer are computed. Moreover, contours of turbulent intensity on the upper surface of the wing and velocity, z-velocity, turbulence intensity and helicity behind the wing are presented. Simulations over the Phase IV blade of NREL are also conducted. The axial velocity behind the rotor, the static pressure and the turbulence intensity contribution on the blade’s surface and the static pressure contours at several blade cross-sections are studied. First of all, the computational study of the two-phase flow over a NACA 0012 airfoil and Liquid Water Content LWC=30 g/m3 is conducted, because there are published experimental data for comparison, in order to validate the CFD developed model. After that, simulations of two-phase flow over the S809 airfoil, S809 wing and Phase IV blade are made. In addition, computational study of the effects of different Liquid Water Content on the aerodynamic performance of NACA 0012 and S809 airfoil at low Reynolds number is made. The results from two-phase flow are compared with the corresponding results from one-phase flow in order to show the effects of two-phase flow at the aerodynamic performance. The influence of two-phase flow on the power coefficient of a wind turbine is also investigated. The results show that the aerodynamic performance degrades when encountering rain, especially lift is degreased and drag is increased. The aerodynamic degradation is caused by the water film formation on the airfoil’s surface and the cratering effects from the raindrops impact. The presence of uneven water film on the airfoil surface roughens the airfoil surface and increases the airfoil thickness. The cratering effects from the water droplets impact on the water film layer increase also the airfoil thickness. Moreover, the droplets splash-back when they impact the airfoil and as a result droplets with smaller diameter and velocity are formed. The acceleration of the splashed-back droplets by the air flowfield acts as a momentum sink, deenergizing the boundary layer and leaving it more susceptible to separation. The aim of the study of the aerodynamic behavior of blades in two-phase flow is the construction of wind turbines with greater efficiency and the production of energy from wind with low cost.

Στην παρούσα διπλωματική εργασία πραγματοποιείται μελέτη της επίδρασης ενός συστήματος διασύνδεσης ανεμογεννήτριας με το δίκτυο μέσης τάσης, όσον αφορά την έγχυση αρμονικών από την ανεμογεννήτρια προς το δίκτυο. Η μελέτη περιλαμβάνει τη μοντελοποίηση κάθε τμήματος ενός δικτύου μέσης τάσης με πραγματικά στοιχεία που πάρθηκαν από τη ΔΕΗ. Επιπλέον, μοντελοποιείται μια ανεμογεννήτρια τύπου σύγχρονης μόνιμων μαγνητών, καθώς και όλοι οι ηλεκτρονικοί μετατροπείς ισχύος που απαιτούνται για τη διασύνδεση αυτής με το δίκτυο μέσης τάσης. Όλα τα επιμέρους μοντέλα διασυνδέονται μεταξύ τους ώστε να αποτελέσουν το συνολικό μοντέλο που περιλαμβάνει το δίκτυο μέσης τάσης με τις γραμμές μεταφοράς, τους μετασχηματιστές, τα φορτία και τους πυκνωτές αντιστάθμισης που διαθέτει και την ανεμογεννήτρια μαζί με τους μετατροπείς συνδεδεμένη σε ένα ζυγό του δικτύου. Από το μοντέλο αυτό εξάγεται το αρμονικό περιεχόμενο τόσο στο ζυγό της ανεμογεννήτριας όσο και σε διπλανούς ζυγούς και συγκρίνεται κάθε φορά με το πρότυπο IEC 61000-3-6 που θέτει τα όρια για την έγχυση αρμονικών στο δίκτυο μέσης τάσης. Αρχικά, γίνεται μια αναφορά στα πλεονεκτήματα της αιολικής ενέργειας έναντι των άλλων συμβατικών μορφών ενέργειας. Επιπλέον, αναφέρονται τα προβλήματα που δημιουργούνται από την παρουσία ανώτερων αρμονικών στο δίκτυο. Έπειτα, πραγματοποιείται θεωρητική ανάλυση κάθε τμήματος του δικτύου μέσης τάσης. Για κάθε τμήμα δημιουργείται, επίσης, ένα μοντέλο στο Matlab/Simulink που βασίζεται σε στοιχεία από τη ΔΕΗ και εξισώσεις που αφορούν το εκάστοτε τμήμα. Στη συνέχεια, αναλύεται και μοντελοποιείται στο Matlab/Simulink και κάθε τμήμα του συστήματος διασύνδεσης της ανεμογεννήτριας με το δίκτυο μέσης τάσης που περιλαμβάνει την ανεμογεννήτρια τύπου σύγχρονης μόνιμων μαγνητών, ανορθωτική διάταξη, μετατροπέα συνεχούς τάσης σε συνεχή, τριφασικό αντιστροφέα, φίλτρο και μετασχηματιστή. Για κάθε μετατροπέα αναλύεται και μοντελοποιείται και έλεγχος κλειστού βρόχου. Τέλος, τα επιμέρους μοντέλα συγκροτούνται σε ένα συνολικό μοντέλο. Μεταβάλλοντας το συντελεστή φόρτισης των μετασχηματιστών του δικτύου, το συντελεστή ισχύος και το πλήθος των ανεμογεννητριών που συνδέονται στον ίδιο ζυγό, εξάγεται το αρμονικό περιεχόμενο για διάφορους ζυγούς και ελέγχεται αν αυτό συμφωνεί με το πρότυπο IEC 61000-3-6. Επιπλέον, εξετάζεται το αρμονικό περιεχόμενο της τάσης, όταν ο ζυγός διασύνδεσης αλλάζει θέση μέσα στο υπό μελέτη δίκτυο, καθώς και όταν αλλάζει η ταχύτητα του ανέμου. / In the present diploma thesis a study of the effect of a wind turbine system interconnection with medium voltage grid is conducted, in terms of harmonic injection of the turbine to the grid. The study includes the modeling of each part of a medium voltage grid with actual data taken from the Public Power Corporation of Greece. Furthermore, a synchronous wind turbine with permanent magnets is modeled along with all the electronic power converters that are required for the interconnection with the medium voltage grid. All the individual models are connected together to form the overall model, including the medium voltage transmission lines, transformers, loads and their compensation capacitors, along with the wind turbine and the power converters. The wind turbine is connected in a bus of the medium voltage grid. Then, the harmonic content in both the bus of the wind turbine and in other buses of the grid is extracted and compared with the standard IEC 61000-3-6 that sets the limits for harmonic injection in medium voltage grid. Initially, the advantages of use of wind energy are mentioned compared to other conventional forms of energy. Furthermore, the problems caused by the presence of harmonics in the grid are discussed. A theoretical analysis of each section of the medium voltage grid is conducted. For each section, a model in Matlab / Simulink is also created, based on data from the Public Power Corporation of Greece and on equations related to each section. Then, every part of the system connecting the wind turbine to the medium voltage grid is analyzed and modeled in Matlab/Simulink. This system consists of a synchronous wind turbine with permanent magnets, DC-to-DC converter, three-phase inverter, filter and transformer. A closed loop control for each converter is also analyzed and modeled. Finally, the aforementioned models are built up into one single model. By varying the load factor of the grid transformers, the power factor and the number of wind turbines connected to the same bus, the harmonic content of various buses is extracted and checked if it complies with the standard IEC 61000-3-6. Moreover, the harmonic content of the voltage is calculated, when the system of wind turbine changes position in the present medium voltage grid, as well as when the wind speed is changed.

This thesis presents a supervisory hybrid controller for the automatic operation and control of a wind energy conversion and battery storage system. The supervisory hybrid control scheme is based on a radically different approach of modeling and control design, proposed for the subject wind energy conversion and battery storage system. The wind energy conversion unit is composed of a 360kW horizontal axis wind turbine mechanically coupled to an induction generator through a gearbox. The assembly is electrically interfaced to the dc bus through a thyristor-controlled rectifier to enable variable speed operation of the unit. Static capacitor banks have been used to meet reactive power requirements of the unit. A battery storage device is connected to the dc bus through a dc-dc converter to support operation of the wind energy conversion unit during islanded conditions. Islanding is assumed to occur when the tiebreaker to the utility feeder is in open position. The wind energy conversion unit and battery storage system is interfaced to the utility grid at the point of common coupling through a 25km long, 13.8kV feeder using a voltage-sourced converter unit. A bank of static (constant impedance) and dynamic (induction motor) loads is connected to the point of common coupling through a step down transformer. A finite hybrid-automata based model of the wind energy conversion and storage system has been proposed that captures the different operating regimes of the system during grid-connected and in islanded operating modes. The hybrid model of the subject system defines allowable operating states and predefines the transition paths between these operating states. A modular control design approach has been adapted in which the wind energy conversion and storage system has been partitioned along the dc bus into three independent system modules. Traditional control schemes using linear proportional-plus-integral compensators have been used for each system module with suitable modifications where necessary in order to achieve the required steady state and transient performance objectives. A supervisory control layer has been used to combine and configure control schemes of the three system modules to suite the requirements of system operation during any one operating state depicted by the hybrid model of the system. Transition management strategies have been devised and implemented through the supervisory control layer to ensure smooth inter-state transitions and bumpless switching among controllers. It has been concluded based on frequency domain linear analysis and time domain electromagnetic transient simulations that the proposed supervisory hybrid controller is capable of operating the wind energy conversion and storage system in both grid-connected and in islanded modes under changing operating conditions including temporary faults on the utility grid.

This thesis presents a supervisory hybrid controller for the automatic operation and control of a wind energy conversion and battery storage system. The supervisory hybrid control scheme is based on a radically different approach of modeling and control design, proposed for the subject wind energy conversion and battery storage system. The wind energy conversion unit is composed of a 360kW horizontal axis wind turbine mechanically coupled to an induction generator through a gearbox. The assembly is electrically interfaced to the dc bus through a thyristor-controlled rectifier to enable variable speed operation of the unit. Static capacitor banks have been used to meet reactive power requirements of the unit. A battery storage device is connected to the dc bus through a dc-dc converter to support operation of the wind energy conversion unit during islanded conditions. Islanding is assumed to occur when the tiebreaker to the utility feeder is in open position. The wind energy conversion unit and battery storage system is interfaced to the utility grid at the point of common coupling through a 25km long, 13.8kV feeder using a voltage-sourced converter unit. A bank of static (constant impedance) and dynamic (induction motor) loads is connected to the point of common coupling through a step down transformer. A finite hybrid-automata based model of the wind energy conversion and storage system has been proposed that captures the different operating regimes of the system during grid-connected and in islanded operating modes. The hybrid model of the subject system defines allowable operating states and predefines the transition paths between these operating states. A modular control design approach has been adapted in which the wind energy conversion and storage system has been partitioned along the dc bus into three independent system modules. Traditional control schemes using linear proportional-plus-integral compensators have been used for each system module with suitable modifications where necessary in order to achieve the required steady state and transient performance objectives. A supervisory control layer has been used to combine and configure control schemes of the three system modules to suite the requirements of system operation during any one operating state depicted by the hybrid model of the system. Transition management strategies have been devised and implemented through the supervisory control layer to ensure smooth inter-state transitions and bumpless switching among controllers. It has been concluded based on frequency domain linear analysis and time domain electromagnetic transient simulations that the proposed supervisory hybrid controller is capable of operating the wind energy conversion and storage system in both grid-connected and in islanded modes under changing operating conditions including temporary faults on the utility grid.