L'éolienne Darrieus connaît un intérêt accru ces dernières années parce qu'elle représente une solution alternative potentielle de production d'électricité dans les milieux urbains. En particulier,une éolienne de forte solidité peut être choisie car certaines de ses propriétés peuvent être avantageuses pour son implantation proche de zones habitées. A l'inverse, certaines difficultés aérodynamiques émergent. Ce type d'éolienne fonctionne à de faibles vitesses réduites pour lesquelles le décrochage dynamique a un rôle très significatif. L'objectif de ce travail de thèse consiste à compléter la connaissance du phénomène de décrochage dynamique sur une éolienne à axe vertical afin d'améliorer les modèles numériques de prédiction existants. Cette étude s'appuie sur une analyse combinée de résultats numériques et expérimentaux. Les simulations numériques sont produites avec une méthode des panneaux bidimensionnelle instationnaire. Les effets de la viscosité sont introduits par des corrections utilisant notamment un modèle semi-empirique de décrochage dynamique. Le travail expérimental s'est concentrée sur la dynamique tourbillonnaire à proximité immédiate du rotor résultante du décrochage dynamique. Le montage se compose d'une éolienne à pale droite placée dans une soufflerie. Des mesures instationnaires de la répartition de pression pariétale le long de la corde et des mesures de champ de vitesse par vélocimétrie par images de particules ont été accomplies. Les résultats révèlent la manière dont les caractéristiques du décrochage dynamique sont conditionnées par la vitesse réduite. Le retard au décrochage, l'intensité de l'effet du tourbillon de décrochage dynamique et sa convection ont été quantifiés. Enfin, un examen critique de l'applicabilité du modèle de Leishman-Beddoes pour simuler efficacement les effets du décrochage dynamique a été réalisé. / The Darrieus wind turbine has entered a period of renewed interest over the last years because it may stand for an alternative solution to produce electricity in urban areas. In particular, high solidity wind turbine can be chosen to take benefit from some of its key properties for use near populated city areas. Conversely, some aerodynamic problems arise. This type of wind turbine operates at low tip speed ratio for which dynamic stall has a very significant role. The goal of this work is to provide valuable data to complement the knowledge of the dynamic stall phenomenon that occurs on a vertical axis wind turbine in order to improve existing numerical models. This study relies on a combined analysis of numerical and experimental results. The numerical simulations are based on a bidimensional unsteady vortex panel method. Effects of viscosity are introduced by adding corrections computed with a semi-empirical dynamic stall model. The experimental work focuses on the dynamics of the shed vortices existing in the vicinity of the rotor as a result of dynamic stall. The set-up consists of a straight-bladed wind turbine tested in a wind tunnel. Unsteady pressure distribution measurements along the chord and velocity fields measurements by particle image velocimetry were carried out. Results indicate how the characteristics of dynamic stall are conditioned by the tip speed ratio. Stall inception delay, magnitude of the dynamic stall vortex effects and its convection velocity were evaluated. Blade/Vortex interaction was analyzed through the observation of the vortical system downstream of the rotor. In addition, a critical review of the suitability of the Leishman-Beddoes model to effectively simulate the effects of dynamic stall was accomplished.

A renewed commitment in the United States and abroad to electricity from renewable resources, such as wind, along with the recent deployment of very large turbines that rise to new heights, makes obtaining the most efficient and safe designs of the structures that support them ever more important. Towards this goal, the present research seeks to understand how optimization concepts and Microsoft Excel's optimization capabilities can be used in the design of wind turbine towers and foundations. Additionally, this research expands on the work of previous researchers to study how considering the tower and foundation as an integral system, where tower support conditions are not perfectly rigid, affects the optimal design. Specifically, optimization problems are formulated and solved with and without taking into account the effect of deflections, resulting from the foundation's rotational and horizontal stiffness, on natural frequency calculations. The general methodology used to transcribe the design of wind turbine towers and foundations into an optimization problem includes: 1) collecting information on design requirements and parameter values 2) deciding how to analyze the structure 3) formulating the optimization problem 4) implementation using Microsoft Excel. Key assumptions include: 1) use of an equivalent lumped mass method for estimating natural frequency 2) International Electrotechnical Commission (IEC) 61400-1 extreme loading condition controls design (i.e. fatigue loading condition is not considered) 3) extreme loads are obtained from manufacturer provided structural load document that satisfies loading cases outlined in IEC 61400-1 4) wind forces on the tower are calculated in accordance with IEC 61400-1 5) optimization variables are continuous. The sum of the tower material and fabrication cost and the total foundation cost is taken as the objective function. Important conclusions from this work include: 1) optimization concepts and Microsoft Excel's optimization capabilities can be used to obtain reasonable conceptual level designs and cost estimates 2) detailed designs and cost estimates could be achieved using a solver capable of handling discrete optimization problems 3) considering the tower and foundation as an integral system results in a more expensive, but safer, design 4) for the assumed parameter values, the constraint on the tower's natural frequency was found to control the tower design and the bearing capacity constraint was found to control the foundation design 5) relaxing or tightening the limit on the natural frequency will result in the greatest benefit or penalty, respectively, on the optimum solution.

Forced convective drying using a wind turbine mechanically connected to a ventilation fan was hypothesized for low cost and rapid grain drying in developing countries. The idea was tested using an expandable wind turbine blade system with variable pitch, at low wind speeds in a wind tunnel. The design was based on empirical and theoretical models embedded in a graphical user interface (GUI) created to estimate airflow-power requirements for drying ear corn. Output airflow (0.0016 - 0.0052 m3kg-1s-1) increased within the study wind speed range (2.0 - 5.5 m/s). System efficiency peak (8.6%) was observed at 3.5 m/s wind speed. Flow resistance was overcome up to 1m fill depth in 0.5 m x 0.5 m wide drying bin. Drying study at different airflow rates (no forced convection, 0.002 m3kg-1s-1 and 0.008 m3kg-1s-1) were conducted in a controlled environment at 35oC and 45% relative humidity with mean drying time; 40.3, 37.9 and 22.9 h respectively, that reduced with increasing airflow while drying the ear corn from 22% to 15% moisture content. The overall result supports the hypothesis that the wind convection system increased grain drying rates and should be further developed.

In this thesis, the effect of pitch based optimisation was explored for a 80 turbine wind farm. Using a modified Jensen wake model and the Particle Swarm Optimisation (PSO) model, a pitch optimisation strategy was created for the dominant turbulence and atmospheric condition for the wind farm. As the wake model was based on the FLORIS model developed by P.M.O Gebraad et. al., the wake and power model was compared with the FLORIS model and a -0.090% difference was found. To determine the dynamic predictive capability of the wake model, measurement values across a 10 minute period for a 19 wind turbine array were used and the wake model under predicted the power production by 17.55%. Despite its poor dynamic predictive capability, the wake model was shown to accurately match the AEP production of the wind farm when compared to a CFD simulation done in FarmFlow and only gave a 3.10% over-prediction. When the optimisation model was applied with 150 iterations and particles, the AEP production of the wind farm increased by 0.1052%, proving that the pitch optimisation method works for the examined wind farm. When the iterations and particles used for the optimisation was increased to 250, the power improvement between optimised results improved by 0.1144% at a 222.5% increase in computational time, suggesting that the solution has yet to fully converge. While the solutions did not fully converge, they converged sufficiently and an increase in iterations gave diminishing results. From the results, the pitch optimisation model was found to give a significant increase in power production, especially in wake intensive wind directions. However, the dynamic predictive capabilities will have be improved upon before the control strategy can be applied to an operational wind farm.

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.