Research Article  Open Access
Weipeng Yue, Yu Xue, Yan Liu, "High Humidity Aerodynamic Effects Study on Offshore Wind Turbine Airfoil/Blade Performance through CFD Analysis", International Journal of Rotating Machinery, vol. 2017, Article ID 7570519, 15 pages, 2017. https://doi.org/10.1155/2017/7570519
High Humidity Aerodynamic Effects Study on Offshore Wind Turbine Airfoil/Blade Performance through CFD Analysis
Abstract
Damp air with high humidity combined with foggy, rainy weather, and icing in winter weather often is found to cause turbine performance degradation, and it is more concerned with offshore wind farm development. To address and understand the high humidity effects on wind turbine performance, our study has been conducted with spread sheet analysis on damp air properties investigation for air density and viscosity; then CFD modeling study using Fluent was carried out on airfoil and blade aerodynamic performance effects due to water vapor partial pressure of mixing flow and water condensation around leading edge and trailing edge of airfoil. It is found that the high humidity effects with water vapor mixing flow and water condensation thin film around airfoil may have insignificant effect directly on airfoil/blade performance; however, the indirect effects such as blade contamination and icing due to the water condensation may have significant effects on turbine performance degradation. Also it is that found the foggy weather with microwater droplet (including rainy weather) may cause higher drag that lead to turbine performance degradation. It is found that, at high temperature, the high humidity effect on air density cannot be ignored for annual energy production calculation. The blade contamination and icing phenomenon need to be further investigated in the next study.
1. Introduction
In recent years, there has been an increase in large offshore wind installations to harness the huge potential that wind energy offers off the coasts of the China. While offshore wind resources are abundant, stronger, and blow more consistently than landbased wind resources, offshore wind turbines are often operated in high humidity environment; however, our wind turbine power curve IEC standard measurements are measured onshore and often much low to mid humidity environment; also the wind shear data from offshore wind may be different from onshore; what do these factors mean to wind turbine performance needs to be addressed as the offshore wind energy projects are evaluated for annual energy production (AEP) calculation. When air contains high water vapor for the offshore high humidity condition, the condensation may happen at rotating blade surface, which may cause the blade performance degradation due to contamination or icing during winter weather, as show in Figure 1. Also in the morning, the foggy weather (or rainy day) may often happen around offshore wind farm, in which the air is mixed with water droplets in micrometer scale; when they impinge with rotating blade, this may also cause higher drag, thus having an impact on turbine performance. With greater potential for corrosion from exposure to seawater and high humidity, offshore wind turbines must be designed more robustly (i.e., requiring less maintenance) than landbased turbines due to the high costs of transporting maintenance crews and replacement components to and from offshore wind plant sites, which are the other factors we have to address.
In this paper, however, we will focus on analyzing the humidity aerodynamic effects in terms of air density and wind turbine airfoil/blade/rotor performance; we will further investigate water condensation effects on wind turbine blade using CFD analysis, as the water condensation and droplet thin firm formation often cause the result of blade contamination and icing formation in the winter, that cause turbine performance degradation. The wind shear data of offshore wind and corrosion issue are not included in this paper.
2. Spreadsheet Aerodynamic Density Effect of Damp Air
Wind is a form of solar energy and is a result of the uneven heating of the atmosphere by the sun, the irregularities of earth’s surface, and the rotation of the earth. The terms wind energy and wind power describe the process by which the wind is used to generate mechanical power or electricity. Wind turbines convert the kinetic energy in the wind into mechanical power.
This mechanical power can be used for specific tasks or a generator can convert this mechanical power into electricity. The most popular modern wind turbines are horizontalaxis wind turbines which typically have either two blades or three blades. These threebladed wind turbines are operated “upwind,” with the blades facing into the wind. Larger wind turbines are more cost effective and are grouped together into wind farms, which provide bulk power to the electrical grid. The power of wind turbine generated can be calculated as follows: where the wind power is directly proportional to the density of air and the aerodynamic coefficient of wind turbine rotor. Table 1 was computed from Barani Air Density Calculator where it shows that humidity (quantified as relative humidity (RH) in percentage) effect on air density can be ranging from 0.57 to 1.92% depending on temperature and thus the same effect on power. Relative humidity is the ratio of the partial pressure of water vapor to the equilibrium vapor pressure of water at a given temperature.
 
Reference is Barani Air Density Calculator as atmosphere pressure assumed to be the same. 
We take one offshore wind project as an example from southern east China. As the firstorder effect, the average temperature is 20 degrees C and average humidity is 80% RH; the effect on power will be 0.41% ( from Table 1), less power than standard power evaluated at 25% RH. If the wind farm power is 200 MW powered by 4 MW turbine, the power lost due to high humidity is about 4,500,000 kWh per year from theoretical calculation, which is considerable.
In the next sessions, we will study the high humidity effects on wind turbine aerodynamic coefficient, through analyzing the aerodynamic effects on airfoil and blade using CFD method.
3. Airfoil CFD Modeling Study of High Humidity Effects
For the airfoil modeling study, we utilized Fluent CFD analysis following 2D mesh as shown in Figure 2 of a typical airfoil NACA63418 to start with species transport model and film condensation model and then 3D airfoil modeling with Eulerian Wall Film (EWF) model. We use Fluent species transport model without reactions to solve water vapor and air mixing flow around airfoil and then film condensation Fluent UserDefined Function (UDF) model to study the water film formation around airfoil. Last, we use multiphase Eulerian model with film accretion to solve water droplet and air mixing flow around 3D airfoil.
3.1. Species Transport Model
The mixing and transport of the air water vapor mixture were modeled by solving the conservation equations for every species individually, within a phase [1]. The number of phases depends on the way volumetric condensation is treated, which is described in the following equations. A mixture template of air and water vapor was defined with respective material properties. An incompressible ideal gas approximation was used to calculate mixture density. The diffusion energy source for the species equations was considered as shown below, as diffusion might be significant in the low speed nearwall region.
The governing equations are the equations of conservation of mass, momentum, and energy and are written as in [1–4]:
Overall mass is
Momentum is
Energy is
Species mass iswhere is the mixture density, is the pressure, is the shear stress tensor, and is the body force vector. In (4), represents the net source due to viscous dissipation and other work and heat interactions and denotes the net heat flux due to molecular conduction, radiation, and interspecies diffusion.
In (5), is the local mass fraction of each species. is mass diffusion in turbulent flows in the following form:where is the turbulent Schmidt number, is the turbulent viscosity, and is the mass diffusion coefficient for species in the mixture.
The species transport model is best utilized to analyze the aerodynamic effects from mixing air with water vapor without condensation and droplet formation happen.
3.2. Film Condensation Model
When water condensation happens with thin film formation around airfoil, the model for wall condensation is required with the work adopted from the study of Bell (2003) [2] for condensation of superheated water on cold flat plate. It is conceptualized that condensation occurs as water vapor diffuses through the mass transfer boundary layer at the surface and, on contact with the cold surface, releases the latent heat to form a liquid film of water. The liquid film attains the wall temperature and the vapor diffusion and condensation continue at the liquid filmgas phase (airvapor mixture) interface. This process is shown schematically in Figure 3.
A UserDefined Function (UDF) in Fluent has been developed to predict the condensation rate for a mixture of air and water vapor flowing past a specified temperature surface. The effect of the condensation process on the flow and species distribution in the vapor phase is incorporated into the flow calculations through a customized source term applied in walladjacent cells [3]. The model on which the UDF is based contains several assumptions as follows:(1)Vapor phase contains a binary ideal gas mixture of air and water vapor.(2)The liquid phase consists of only water.(3)Local thermodynamic equilibrium exists at the liquidvapor interface.(4)Film condensation only is considered. Dropwise condensation does not take place.(5)The thermal resistance of the liquid film is neglected.
As a consequence of the assumptions on which the model is based, the condensation rate is governed by the rate of diffusion of water vapor towards the cold surface. This allows the condensation rate at the liquid vapor interface to be predicted as follows. Species mass fluxes for air and water vapor at the liquidvapor interface can be written as follows (Bird et al., 1960 [4]): where is condensation mass flux (kg/m^{2}·s). In the UDF, condensation is assumed to take place if the wall temperature is less than or equal to the saturation temperature corresponding to the partial pressure of water vapor at the surface. If the temperature is above the saturation temperature, the water vapor mass fraction is set to a value equal to the value in the wal adjacent cell. In order to satisfy Assumption (), if the temperature is less than or equal to the saturation temperature, the water vapor mass fraction at the wall used to represent the liquidvapor interface is assigned a value such that the partial pressure of the water vapor equals the saturation pressure for water at the local wall temperature as follows:where is volumetric mass source (kg/m^{3}·s). Equations (8) are used as a source terms within the source code of the UDF (Table 2).

3.3. Eulerian Wall Film
The Fluent Eulerian Wall Film (EWF) model can be used to predict the creation and flow of thin liquid films on the surface of walls [5].
In EWF model, the film is not resolved using mesh normal to the surface; instead this model uses a virtual film on the surface so the film does not affect the core flow momentum field. In many cases, you would be interested in modeling a thin film flow which can separate, strip, and evaporate but does not affect the core flow field. This kind of thin film modeling is computationally expensive in general multiphase framework where you will need very fine mesh to model the film and also need to calculate interphase flux accurately. EWF model can only be used with 3D solver and EWF model assumes that film always flows parallel to the surface so normal component of film velocity is zero. The film is assumed to have a parabolic velocity profile and bilinear temperature profile across its depth.
Conservation of mass for a twodimensional film in a threedimensional domain iswhere is the liquid density, the film height, the surface gradient operator, the mean film velocity, and the mass source per unit wall area due to water droplet collection.
Conservation of film momentum is given as follows:where
In an Eulerian multiphase flow, the secondary phase to be collected on a solid surface must have the same material as defined in the Eulerian Wall Film model. When the secondary phase is captured by the wall surface, its mass and momentum are removed from the multiphase flow and added as source terms to the continuity and momentum equations, respectively, of the wall film.
The mass source term is given by , where is the secondary phase volume fraction, is the secondary phase density, is the phase velocity normal to the wall surface, and is the wall surface area. The momentum source is given by , where is the secondary phase velocity vector.
The collection efficiency () is calculated as the ratio of the amount or mass of water collected on the surface to the total mass of droplets being injected at the inlet. Eulerian methods can be used to determine the collection efficiency.
Collection efficiency is defined aswhere is the impingement mass flux, that is, the flow rate of the particle stream impinging on the face and is given by . is the farfield liquid water content (kg/m^{3}) and , the phase velocity, is the farfield flow velocity (m/s). In the calculation of collection efficiency, these values are used as reference values only and do not affect the flow field in any way. At least one iteration is required to do collection efficiency calculation. Easily specifying dropletfree stream condition and smooth impingement region can be achieved. Each droplet size as a separated phase with one set of transport equation can be solved.
It is always good to first achieve a converged multiphase flow solution and then solve for Eulerian Wall Film. It is always good to start with a smaller film timestep size. So, for steadystate calculations, we select the Adaptive Time Stepping option to set the Max. Courant Number and the Initial Time Step (the Eulerian Wall Film model is always transient).
3.4. Solution Method and Results
A variety of spatial and temporal discretization schemes and turbulent models are also available in Fluent. For the present simulations, the pressurebased solver was chosen, as the flow can be considered to be in the incompressible fluid. The full threedimensional NavierStokes equations were solved using a steady, implicit approach. The SemiImplicit Pressure Linked EquationsConsistent (SIMPLEC) algorithm 2 in [6] was used to treat pressurevelocity coupling for stability. The thirdorder Monotone UpstreamCentered Schemes for Conservation Laws (MUSCL) were used to derive the face values of different variables for the spatial discretization, which were used to compute the convective fluxes. The RhieChow scheme 3 in [6] was selected to interpolate pressure at the control volume face and to satisfy the mass conservation without pressure oscillations. For simulating turbulence, the shear stress transport (SST) model was used.
The aerodynamic results of lift and drag coefficient are summarized in Table 3. Through the comparison of two different models of moisture and humidity at RH 25% and RH 90% of the lift coefficient, drag coefficient, and lift drag ratio, it was found that, in the species transport model and VF model, 25% air humidity of the lift to drag ratio is higher than 90% humidity, and lift coefficient and drag coefficient curve are not exactly the same; difference comparison between the two found that numerical difference is very small (usually less than 1%), considering that the numerical calculation error can be due to moisture changes on the lift and drag coefficient of little effect.

The lift coefficient calculated by the Eulerian Wall Film model for foggy weather condition is smaller than that calculated by other models, and the drag coefficient is significantly larger than that of the other models, which leads to a significant reduction in the lift drag ratio. It can be seen that, with the increase of humidity, the lift coefficient decreases, the drag coefficient increases, and the lift drag ratio decreases. The similar results were shown with heavy rain fall study from [7–10]. Because the drag coefficient is larger than other models, there may be some uncertainties; so this paper is mainly to study the condensation distribution on the surface of the airfoil by the Euler wall model. The condensation of water vapor and the impact of droplets on the surface of an airfoil or a blade can lead to other consequences of aerodynamic degradation, such as blade contamination and icing as shown in Figures 4 and 5.
When the condensation model is used, it is found that when the temperature difference between the damp air and airfoil surfaces becomes larger, the water mass flow rate of the airfoil surface is larger, and when the humidity is higher, the condition is more obvious. But around airfoil, Figure 6 shows the higher temperature difference and the more condensation mass flow rate curve have a similar shape as pressure curve occurs. With different humidities, airfoil surface pressure coefficient curves remained basically unchanged. Airfoil surface condensation cloud mass flow rate can be found as shown in Figure 7 that the condensate mass flow rate larger region at the leading edge and the trailing edge region, mainly airfoil leading edge surface and the upper surface of the trailing edge, so, in the actual wind turbine operation, it would be the region due to the presence of surface water condensation, is more likely to cause contamination on the surface of (insects, such as particles) adhesion as shown in Figure 8.
Figure 9 shows the approximate separation point air velocity vector diagram of the different height and various working conditions of the Eulerian Wall Film model. A crosssectional comparison of different wind speed, temperature, humidity, and velocity vector found the same temperature and wind speed and humidity 90% compared with 25% humidity premature separation; when humidity is higher, separation is stronger and position is closer to the front, which more easily leads to separation.
Figure 10 shows the flow pattern and water vapor condensation area at leading edge and some at trailing edge of airfoil. Figure 11 shows that RH 90% condition has much higher water film thickness in comparison with RH 25% condition and high angle of attack (AOA) has more water condensation in the range of 70 μm water film thickness. Figure 12 shows that the EWF modeling is transient and the water condensation film thickness changes with time.
Figures 13 and 14 are corresponding for water vapor condensation in terms of water film thickness around airfoil for both AOA = 6 and 10 degrees, both showing the airfoil surface water content analysis; it is found that the water condense most in front of the airfoil near the leading edge but much less near the trailing edge. And it is found that the negative pressure side has more water condensation than pressure side, while the AOA = 10 degrees is much significant than AOA = 6 degrees.
Table 4 summarized the RH effects at different AOA for Cl, Cd, water volume fraction, sum of water condensation film, and water film thickness.

In general, through on airfoil surface wet air flow analysis, it is found that the moisture of airfoil dynamic performance impact is small; the moisture effects on airfoil are mainly moisture in the airfoil surface condensation or convergence, that is easy to cause pollution of adhesion. By analyzing the water vapor condensation on the airfoil surface and the distribution of the surface water film, it provides reference for the understanding airfoil degradation because water condensation can lead to blade icing and contamination.
4. Rotating Blade CFD Modeling Study of High Humidity Effect
4.1. Result Analysis Based on Film Condensation Model
Like airfoil study, we also utilized Fluent CFD analysis; to study the effect of humidity on blade performance, used two models, species transport and film condensation UDF as stated above.
Computation using a typical 1.5 MW blade model, for threeblade rotor, considers the 120degree rotation cycle; numerical calculation is only performed on the leaf. The spherical region of 15 times the impeller radius is chosen as the computational domain. In order to improve the calculation accuracy and computation time saving, we used full hexahedral grid. Figure 15 shows the overall outline of the grid and local amplification. Fine mesh is used in the near wake region as blade root and tip. The mesh size of the tip wake area is less than 0.05 m, near the rotor area is less than 2.5 m, and outside of compute domain is less than 20 m. The total number of grids nodes is about 6.2 million. omega SST model is the most commonly used in the field of wind turbine. Table 4 shows the component transport model and film condensation model calculation of the aerodynamic parameters of the blade; it can be seen that, in the two models, with the humidity increasing, the power is reduced, and 90% relative humidity compared to the 25% relative humidity power decreases by about 1%.
Table 5 summarized the result that higher humidity lower the power it generated primarily due to density effect. Figure 16 shows that the condensation mass flow rate distribution on the surface of the blade is calculated by the condensation model. Negative on behalf of condensation water, you can see that most of the area of the surface of the blade has the condensation of water, along the blade airfoil section direction, and the front position has maximum mass flow rate. Along the direction of the blade, the blade root to tip, the whole condensation mass flow rate increases.

4.2. Result Analysis Based on Eulerian Wall Film Model
The water volume fraction curves of blades with different section positions at 25% and 90% relative humidity are shown in Figure 17. The higher the humidity is, the higher the volume fraction of water is. Humidity is the same, from root to leaf tip, and numerical water volume fraction increases rapidly. In the inner region of the blade (near root), the water is mainly distributed in the position of the droplet impingement region and the trailing edge in suction side of the blade, and the outer side of the blade (near blade tip) is only distributed in the leading edge. Meantime, water distribution can be found clearing in Figure 18; therefore, actual environment, high humidity in cold area, and the location in winter are the most easily freeze position, and arid sandy area is the location with the most polluted areas. These locations require special protection.
(a)
(b)
5. Conclusions
With the background of high humidity environment, the density effect and the airfoil/blade aerodynamic performance of different humidity are simulated by using CFD Fluent with the component transport, UDF thin film condensation model, and the Eulerian Wall Film model. It is found that the high humidity has considerable effect on density at high temperature condition and little effect on airfoil/blade aerodynamic performance, and it was found that the water vapor condensation on the surface of the airfoil or a blade may occur around the leading edge and trailing of airfoil, which then leads to blade contamination in the dirty environment such as dust and insect or icing in winter and low temperature conditions. When that happens, the wind turbine performance may be severely degraded, especially for the airfoil sensitivity to the surface. In a foggy day, the study found that water droplets in the air with high humidity are impacted and collected on the airfoil surface to form a water film, which affects the flow of the boundary layer, resulting in much higher drag and thus losing good performance. The similar behaviors were found by other studies that found the loss of performance in rainy day.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge and thank Dr. Chen Chunmei, Dr. Gao Wei, and Chief Cai Xiaofeng from China Creative Wind Energy Inc. for support and sponsorship of the project. The authors would like to thank Dr. Tang Hongfeng and Dr. Wang Hu from CDT Science and Technology Research Institute for support of this project. The authors would like to thank Professor Wang Xiaodong from China North Electric Power University for many interesting discussions and helpful suggestions.
References
 A. Kumar and S. Mazumder, “Coupled solution of the species conservation equations using unstructured finitevolume method,” International Journal for Numerical Methods in Fluids, vol. 64, no. 4, pp. 409–442, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Bell, Application Brief: Film Condensation of Water Vapor, Fluent, Inc., Lebanon, NH, USA, 2003.
 K. Das, C. Manepally, R. Fedors, and D. Basu, Numerical and Experimental Study of InDrift Heat and Mass Transfer Progesses, Center for Nuclear Waste Regulatory Analyses, San Antonio, Tex, USA, 2011.
 R. B. Bird, W. E. Stewart, and E. N. Lightfoot, Transport Phenomena, John Wiley and Sons, New York, NY, USA, 1960.
 ANSYS Inc, ANSYSFLUENT® Version 12.1 User's Guide, ANSYS, Inc, Canonsburg, Pa, USA, 2009.
 D. K. Kolmogorov, W. Z. Shen, N. N. Sørensen, and J. N. Sørensen, “Fully consistent CFD methods for incompressible flow computations,” Journal of Physics: Conference Series, vol. 524, no. 1, Article ID 012128, 2014. View at: Publisher Site  Google Scholar
 E. C. Douvi and D. P. Margaris, “Aerodynamic performance investigation under the influence of heavy rain of a NACA 0012 airfoil for wind turbine applications,” International Review of Mechanical Engineering, vol. 6, no. 6, pp. 1228–1235, 2012. View at: Google Scholar
 M. Cai, E. Abbasi, and H. Arastoopour, “Analysis of the performance of a windturbine airfoil under heavyrain conditions using a multiphase computational fluid dynamics approach,” Industrial and Engineering Chemistry Research, vol. 52, no. 9, pp. 3266–3275, 2013. View at: Publisher Site  Google Scholar
 A. C. Cohan and H. Arastoopour, “Numerical simulation and analysis of the effect of rain and surface property on windturbine airfoil performance,” International Journal of Multiphase Flow, vol. 81, pp. 46–53, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Wu and Y. Cao, “Numerical simulation of flow over an airfoil in heavy rain via a twoway coupled EulerianLagrangian approach,” International Journal of Multiphase Flow, vol. 69, pp. 81–92, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 “Converter for humidity of air,” http://www.cactus2000.de/uk/unit/masshum.shtml. View at: Google Scholar
 W. C. Reynolds, Thermodynamic Properties in SI, 1979.
 N. Harper, “Detecting Ice on Windturbine Blade,” 2011, http://www.windpowerengineering.com/maintenance/detectingiceonwindturbineblades/. View at: Google Scholar
 M. WadhamGagnon, C. Farley, and C. Arbez, “Wind site considerations: in cold climate conditions,” 2013, http://www.nacleanenergy.com/articles/16824/windsiteconsiderationsincoldclimateconditions. View at: Google Scholar
 L. Rempel, “Rotor blade leading edge erosionreal life experiences,” Wind Systems Magazine, pp. 22–24, 2012. View at: Google Scholar
 C. J. Spruce, “Power performance of active stall wind turbines with blade contamination,” in Proceedings of the European Wind Energy Conference, Athens, Greece, 2006. View at: Google Scholar
Copyright
Copyright © 2017 Weipeng Yue et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.