Numerical simulations of the Navier-Stokes equations are performed to achieve a better understanding of the behaviour of wakes generated by wind turbines. The simulations are performed by combining the in-house developed computer code EllipSys3D with the actuator line and disc methodologies. In the actuator line and disc methods the blades are represented by a line or a disc on which body forces representing the loading are introduced. The body forces are determined by computing local angles of attack and using tabulated aerofoil coefficients. The advantage of using the actuator disc technique is that it is not necessary to resolve blade boundary layers. Instead the computational resources are devoted to simulating the dynamics of the flow structures. In the present study both the actuator line and disc methods are used. Between approximately six to fourteen million mesh points are used to resolve the wake structure in a range from a single turbine wake to wake interaction in a farm containing 80 turbines. These 80 turbines are however represented by 20 actuator discs due to periodicity because of numerical limitations. In step one of this project the objective was to find a numerical method suitable to study both the flow structures in the wake behind a single wind turbine and to simulate complicated interaction between a number of turbines. The study resulted in an increased comprehension of basic flow features in the wake, but more importantly in the use of a numerical method very suitable for the upcoming purpose. The second objective of the project was to study the basic mechanisms controlling the length of the wake to obtain better understanding of the stability properties of wakes generated by wind turbine rotors. The numerical model was based on large eddy simulations of the Navier-Stokes equations using the actuator line method to generate the wake and the tip vortices. To determine critical frequencies the flow is disturbed by inserting a harmonic perturbation. The results showed that instability is dispersive and that growth occurs only for specific frequencies and mode types. The study also provides evidence of a relationship between the turbulence intensity and the length of the wake. The relationship however needs to be calibrated with measurements. In the last project objective, full wake interaction in large wind turbine farms was studied and verified to measurements. Large eddy simulations of the Navier-Stokes equations are performed to simulate the Horns Rev off-shore wind farm 15 km outside the Danish west coast. The aim is to achieve a better understanding of the wake interaction inside the farm. The simulations are performed by using the actuator disc methodology. Approximately 13.6 million mesh points are used to resolve the wake structure in the park containing 80 turbines. Since it is not possible to simulate all turbines, the 2 central columns of turbines have been simulated with periodic boundary conditions. This corresponds to an infinitely wide farm with 10 turbines in downstream direction. Simulations were performed within plus/minus 15 degrees of the turbine alignment. The infinitely wide farm approximation is thus reasonable. The results from the CFD simulations are evaluated and the downstream evolution of the velocity field is depicted. Special interest is given to what extent production is dependent on the inflow angle and turbulence level. The study shows that the applied method captures the main production variation within the wind farm. The result further demonstrates that levels of production correlate well with measurements. However, in some cases the variation of the measurement data is caused by the different measurement conditions during different inflow angles. / QC 20100720

No description available.

Des projets d'installation de parcs de plusieurs hydroliennes rapprochées en rivières ou dans les océans ont été récemment démarrés, afin de développer cette source d'énergie renouvelable. Dans ces parcs, les interactions de sillage entre les hydroliennes doivent être calculées puisqu'elles peuvent affecter leur puissance produite. Un modèle CFD stationnaire de type disque d'action couplé aux équations RANS est développé dans ce travail pour calculer la puissance produite et l'écoulement au sein d'un parc d'hydroliennes Darrieus. Ce modèle utilise des répartitions détaillées de force dont l'intensité dépend de la position sur la turbine. Elles sont obtenues par des calculs préliminaires URANS de l'écoulement sur la géométrie de la turbine en rotation. De nouvelles lois sont obtenues pour les coefficients de puissance et de force en utilisant la vitesse locale (vitesse au niveau de la turbine) au lieu de la vitesse amont dans leur définition. Ces coefficients deviennent alors indépendants du confinement de la turbine. Ces lois servent à construire un modèle qui calcule les distributions de force représentant chaque turbine du parc en fonction de la vitesse locale du fluide, pour simuler chaque turbine fonctionnant proche de son point de maximum d'efficacité. Une validation du modèle est réalisée par comparaison à de nouvelles expériences d'une turbine Darrieus à échelle réduite. Différentes configurations de parcs sont ensuite simulées par le modèle 3D, ainsi que par une version 2D du modèle. Les distances entre turbines qui permettent d'obtenir une puissance produite par le parc maximale sont notamment recherchées. / New projects have recently been launched to build farms of several tidal or river turbines, which are part of the renewable energy systems. The turbine wake interactions in the farm must be considered, as they can affect the power production of the turbines. A steady-state Actuator force model using the RANS equations is developed in the present work to calculate the power production and the flow through arrays of tidal or river Darrieus turbines. It uses detailed three dimensional force distributions depending on the position on the turbine, obtained beforehand by a set of blade-resolved URANS simulations of the turbine. New power coefficient and force coefficient laws depending on the local velocity (flow velocity at the machine position) instead of the upstream velocity are established and appear to be independent from the local turbine blockage in an array. Those laws are used to construct a model that adapt the Actuator force distributions to the local velocity of the flow reaching each turbine, in order to simulate each turbine functioning close to its maximum efficiency point. The model is validated against experimental measurements on a reduced-scale Darrieus turbine. Different farm configurations are simulated and compared to results of the same model adapted in two dimensions. The distances between turbines that are optimizing the farm power production are especially investigated.

In the present thesis the wake flow behind wind turbines is analyzed numerically using large-eddy simulations. The wind turbine rotors are modeled by using either the actuator disc method or the actuator line method in which the blades are represented by body forces computed with airfoil data. Using these models, the boundary layers of the turbine blades are not resolved and most of the computational power is preserved to simulate the wake flow. The actuator disc method is used for the wake interaction studies of the Lillgrund wind farm. In this study the power production is simulated for two different wind directions and compared to measurements. A grid sensitivity study and a turbulence intensity study are also performed. As a last step the front row turbines are derated in an attempt to increase the total production of the farm. The results show that it is important to impose atmospheric conditions at the inlet in the simulations, otherwise production will be unrealistically low for some turbines in the farm. The agreement between the simulated and measured power is very good. The study based on derating the front row turbines does not show any positive increase on the farm production. The actuator line method is used for near wake analysis of the MEXICO rotor. In this study the near wake is simulated for five different flow cases and compared with particle image velocimetry (PIV) measurements. The analysis is performed by comparing size and circulation of the tip vortices, the radial and streamwise velocity distributions, the spatial expansion of the wake and the axial induction factor. The simulations and measurements generally are in agreement. In some cases, however, the measurements are affected by tunnel effects which are not captured in the simulations. In connection to the actuator disc method a power control strategy for operating conditions below rated power is implemented and tested. The strategy is first validated using an in-house developed blade element momentum code and then is implemented in the actuator disc method used in the EllipSys3D code. The initial tests show that the strategy responds as expected when changing the moment of inertia of the rotor and when varying the inlet conditions. Results from the implementation of the strategy in the actuator disc method in EllipSys3D show that the turbine adapts to the conditions it is operating in by changing its rotational velocity and power output when the inlet conditions are varied. / <p>QC 20130111</p>

<p>In this work, two overarching global crises are addressed with an engineering lens; the COVID-19 pandemic and climate change. Regarding the latter, an investigation is completed into the fluxes of CO2 in the wake of a simple wind farm for identification of potentially beneficial siting of Direct Capture of CO2. In this analysis, large-eddy simulations are used to quantify scalar entrainment in the turbines’ wake for several empirical CO2 profiles. In instances with positive or a combination of CO2 gradients, it was found that the concentration of CO2 increased in wake through downward mixing and entrainment. In a negative CO2 gradient, the opposite was found, with the wind turbine mixing away the increased surface<br> concentration and entraining down lower concentration air from above. These findings were used to make recommendations on scenarios in which wind turbines were beneficial to Direct Capture plants.<br> In addition, as part of the ongoing response to the COVID-19 pandemic, an innovative new technology was designed and constructed; a prototype photoacoustic spectrometer for the rapid detection of viruses. With the vision to become a viral "breathalyzer", the primary stage of development involved the creation of a prototype for proof-of-concept of viral detection using PAS. An extensive literature review was completed to determine optimal<br> design, with several distinct innovations integrated with the end-product in mind; such as a pure silicon resonator cell and a light-emitting diode source for low-cost, portable detection.<br> This was estimated to be of sufficient quality to detect single virions, as found through Finite Element Analysis.<br> Additionally, the validation of a proposed improvement on the medical mask, named Hy-Cu, is shown. Through various tests, Hy-Cu was found to have greater breathability than KN95 or surgical masks, as well as comparable efficiency in filtration of viral droplets.<br> Additionally, the novel inclusion of a diamond-like carbon-coated copper mesh layer resulted in viral inactivation of 99% after a period of 2 hours, allowing Hy-Cu to be safely reused without risk of transmission.<br> </p> <p> </p>

As wind energy technology continues to mature, farm sizes grow and wind farm layout design becomes more difficult, in part due to the number of design variables and constraints. Wind farm layout optimization is typically approached using gradient-free methods because of the highly multi-modal shape of the wind farm layout design space. Gradient-free method performance generally degrades with increasing problem size, making it difficult to find optimal layouts for larger wind farms. However, gradient-based optimization methods can effectively and efficiently solve large-scale problems with many variables and constraints. To pave the way for effective and efficient wind farm layout optimization for large-scale wind farms, we have worked to overcome the primary barriers to applying gradient-based optimization to wind farm layout optimization. To improve model/algorithm compatibility, we adjusted wake and wind farm models, adding more realistic curvature and smoothness to enable optimization algorithms to travel through areas in the design space where they had previously gotten stuck. We reduced the number of function calls required for gradient-based wind farm layout optimization by over three orders of magnitude for large farms by using algorithmic differentiation to compute derivatives. We reduced the multi-modality of the wind farm layout design space using wake expansion continuation (WEC). We developed WEC to work with existing optimization algorithms, enabling them to get out of local optima while remaining fully gradient-based. Across four case studies, WEC found results with lower wake loss, on average, than the other methods we tested. To resolve concerns about optimization algorithms exploiting model inaccuracies, we compared the initial and optimized layouts to large-eddy simulation (LES) results. The simple models predicted an AEP improvement of 7.7% for a low-TI case, and LES predicted 9.3%. For a high-TI case, the simple models predicted a 10.0% improvement in AEP and LES predicted 10.7%. To resolve uncertainty regarding relative solution quality for gradient-based and gradient-free methods, we collaborated with seven organizations to compare eight optimization methods. Each method was managed by researchers experienced with them. All methods found solutions of similar quality, with optimized wake loss between 15.48 % and 15.70 %. WEC with SNOPT was the only purely gradient-based method included and found the third-to-best solution.