International Journal of Engineering Materials and Manufacture (2017) 2(3) 37-48 https://doi.org/10.26776/ijemm.02.03.2017.01 M. H. Hasan Department of Mechanical and Industrial Engineering Ryerson University, Canada E-mail: hasibulhasan@ryerson.ca S. Ahmed Department of Robotics and Mechatronics Engineering Kazi Motahar Hossain Bhaban University of Dhaka, Dhaka 1000, Bangladesh E-mail: shugataahmed.rme@du.ac.bd Reference: Hasan, M. H., and Ahmed, S. (2017). Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating. International Journal of Engineering Materials and Manufacture, 2(3), 37-47. Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating Muhammad Hasibul Hasan and Shugata Ahmed Received: 13 April 2017 Accepted: 08 July 2017 Published: 15 September 2017 Publisher: Deer Hill Publications © 2017 The Author(s) Creative Commons: CC BY 4.0 ABSTRACT Efficiency and durability are critical issues that affect widely-adopted aerofoil-power generator as a sustainable source of electrical power. Even though high wind power density can be achieved; installing wind turbines in desert condition has difficulties including thermal variation, high turbulence and sand storms. Sand blasting on turbine blade surface at high velocities causes erosion resulting turbine efficiency drop. Damage-induced erosion phenomena and aeroelastic performance of the blades needed to be investigated. Suitable coating may prevent erosion to a great extent. A numerical investigation of erosion on NACA 4412 wind turbine blade has been performed using commercial computational fluid dynamics software ANSYS FLUENT 14.5 release. Discrete phase model (DPM) has been used for modelling multi-phase flow of air and sand particles over the turbine blade. Governing equations have been solved by finite volume method (FVM). Conventional 30-70% glass fibre resin and non-conventional jute fibre composite have been used as turbine blade material. Sand particles of 0.3 mm diameter have been injected from 20° , 30°, 45°, 60° and 90° angles at 50℃ temperature. Erosion rate, wall shear stress and strain rate have been calculated for different wind velocities and impingement angles. Simulation results for higher velocities deviate from the results observed at lower wind velocities. In simulation, erosion rate is highest for 30 0 impingement angle at low wind velocities, which has been validated by experiment with a mean absolute error (MAE) of 5.56%. Erosion rate and wall shear stress are higher on jute composite fibre than glass fibre resin. Developed shear stress on wind turbine blade surface is highest for 20 0 impingement angle at all velocities. On the other hand, exerted pressure on turbine blade surface is found highest for 90° angle of attack. Experimental results, with or without Titanium nitride(TiN) nano-coating, also revealed that surface roughness augments with increasing impingement angles. Nano-coating (TiN) by RF sputtering technique reduced the surface roughness significantly as oppose to uncoated samples. Highest roughness has been observed on uncoated blade surface collided with 0.3-0.69 mm diameter brown aluminium oxide particles. Keywords: Turbine blade, Aerofoil, Nano coating, Wear, Wear resistance, Power generation, Turbine efficiency, Fibre composite turbine blade 1. INTRODUCTION Combustion of bio-fuels generates greenhouse gases. Releasing greenhouse gases is one of the major reasons of atmospheric temperature rise. Moreover, amount of underground petroleum storage is limited. Hence, renewable energy, such as wind, solar, geothermal etc., is becoming popular day by day as alternative sources. Wind energy is clean, available and easy to convert to electrical energy. Wind turbines are used to generate electrical power from wind energy. However, selecting turbine blade material is a critical issue. Moreover, in desert condition air flows with sand particles at relatively high velocities. Sand blasting on turbine blade surface results significant erosion. Hence, a sharp decrement of turbine efficiency is observed with time. Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 38 Nomenclature 𝑎 Particle impact angle (°) 𝐴𝑓𝑎𝑐𝑒 Surface area of impacted wall boundary cell (𝑚𝑚 2) 𝑑𝑝 Particle diameter (𝑚𝑚) �⃗� Gravitational acceleration (𝑚𝑠 −2) 𝑘 Turbulence kinetic energy (𝑚2𝑠 −2) �̇�𝑝 Particle mass flow rate (𝑘𝑔𝑠 −1) 𝑛 Material independent index 𝑃 Pressure (𝑃𝑎) 𝑅𝑒 Reynolds number 𝑅𝑒𝑟𝑜𝑠𝑖𝑜𝑛 Erosion rate (𝑘𝑔𝑚 −2𝑠 −1) 𝑇 Temperature (℃) �⃗⃗� Velocity of air (𝑚𝑠 −1) �⃗⃗�𝑝 Particle velocity (𝑚𝑠 −1) 𝑣𝑖 Particle impact velocity (𝑚𝑠 −1) 𝛼 Volume fraction of air 𝜖 Turbulence kinetic energy dissipation rate (𝑚2𝑠 −3) 𝜃 Solid domain temperature (℃) 𝜇 Dynamic viscosity (𝑁. 𝑠𝑚−2) 𝜌 Density of air (𝑘𝑔𝑚−3) 𝜌𝑝 Particle density (𝑘𝑔𝑚 −3) Numerous trials for selecting proper wind turbine blade materials have been remarked in the literature. Mishnaevsky et al. [1] studied the usability of timber as a turbine blade material. Their investigation represents a framework for fabrication, testing and installation of timber blades. Composite materials are also getting attentions for fabricating wind turbine blades due to their low cost and attractive properties. Polymer composites had been considered by researchers to enhance surface resistance to wear such as carbon fibres, short aramids or glass [2-7], polyetheretherketone (PEEK) [8], polyphenylene sulfide (PPS) [9], polyoxymethylene (POM) [10], polyamide (PA) [11] etc. Patnaik et al. [12] numerically and experimentally investigated erosion rate on zinc-aluminum alloy metal (ZA- 27) filled with titania (TiO2) for different particle impact velocity, size, temperature and impingement angle. They observed that highest erosion rate occurs at 60 0 impingement angle. Erosion rate on unfilled base metal was found higher than reinforced alloy. Erosion on composite material was also investigated experimentally by Drensky et al. [13]. It was found that mass loss of the material is highest for 45° impact angle. Wind velocity, particle size and fibre orientation also have great influence on mass loss. Simulation was carried out by Fiore et al. [14] to investigate erosion on wind turbine blade for insects and sand particles impact. Maximum erosion rate was detected on blade’s low pressure side. Wind tunnel testing was conducted by Gaudern et al. [15] to figure out erosion patterns on wind turbine blade surface under varying conditions. Hamed et al. [16] also explored surface damage of turbine blade both numerically and experimentally. Their results show that erosion and surface roughness enhance with increasing impact angle and particle size. Foley et al. [17] studied erosion behaviour of steel AISI 4340 and found that for the same material hardness, the erosion at impact angle of 30 0 was much greater than at 90 0 . On the other hand, Hutchings [18] reported erosion peaks at 90° for AISI 52100 steel. E. Rodríguez et al. [19] found three different wear regions based on impact angles and hardness. Higher amount of wear found for 10° and 20° impact angles at lower hardness value, while, for impact angles of 60°, 75° and 90°, amount of wear is high for higher hardness value. For impact angles of 30° and 40°, wear is almost same for all hardness values. An effective way to prevent erosion on wind turbine blade is to provide coating. Rico et al. [20] accomplished low wear on SAE-42 steel by depositing Al2O3 – 13% TiO2 coating. Liang et al. [21] substantially improved damping properties of wind turbine blade using carbon nanofibre nanocomposite coating. Nano/polymer coating is promising for preventing erosion to a great extent. Zhang et al. [22] tested silica filled transparent acrylate-based coatings and remarked that nanosilica particles improve erosion resistance properties. Figure 1 illustrates coating with different percentage of nanosilica particles. In this paper, a static NACA 4412 single wind turbine blade under desert condition has been simulated by software and laboratory setup. Numerical investigation has been carried out by ANSYS FLUENT 14.5 release. In software, discrete phase model (DPM) has been used to simulate the multiphase flow of air and sand particles. Erosion rate, development of shear stress and strain on wind turbine blade surface have been investigated. Aerodynamic performances of turbine blade have also been explored. Experiment has been carried out using air jet erosion tester. Simulation results show good agreement with experimental data. Hasan and Ahmed (2017): International Journal of Engineering Materials and Manufacture, 2(3), 37-48 39 Figure 1: TEM micrograph of nano/polymer coating with (a) 10wt% and (b) 40wt% nanosilica particles [22]. 2. MATHEMATICAL MODEL Discrete phase model (DPM) has been used to simulate the flow of air and sand particles. In DPM, the fluid phase is considered as continuum and solid particles are treated as discrete phase. Discrete solid particles are tracked in the flow field in a Lagrangian frame of reference, where governing equations are solved for a single particle. On the other hand, transport equations for per unit volume of continuum phase is derived, which is Eulerian approach for modelling multiphase flow. The force balance equation for a single sand particle is following: 𝑑𝑢𝑝 𝑑𝑡 = 𝐹𝐷 (�⃗⃗� − �⃗⃗�𝑝) + �⃗⃗�(𝜌𝑝−𝜌) 𝜌𝑝 + �⃗� (1) As sand particles are very small, it is considered that no gravitational force is acting on them. Hence, gravitational acceleration �⃗� = 0. The term 𝐹𝐷 (�⃗⃗� − �⃗⃗�𝑝) defines the drag force/ unit mass of sand particle where 𝐹𝐷 = 18𝜇𝐶𝐷𝑅𝑒 𝜌𝑝𝑑𝑝 2 24 (2) Relative Reynolds number is defined as 𝑅𝑒 = 𝜌𝑑𝑝|�⃗⃗⃗�𝑝−�⃗⃗⃗�| 𝜇 (3) The governing equations for air is derived by Eulerian method. Continuity equation is the following: 𝜕(𝛼𝜌) 𝜕𝑡 + ∇(𝛼𝜌�⃗⃗�) = 0 (4) Conservation of momentum equation: 𝜕(𝛼𝜌�⃗⃗⃗�) 𝜕𝑡 + 𝛻(𝛼𝜌�⃗⃗��⃗⃗�) = −𝛼𝛻𝑃 + 𝛻𝜏 (5) Standard 𝑘 − 𝜀 turbulence model [12] is used for turbulence flow modelling. Equations related to turbulence: 𝜕(𝜌𝑘) 𝜕𝑡 + ∇(𝜌𝑘𝑢) = ∇(𝛼𝑘 𝜇𝑒𝑓𝑓 ∇𝑘) + 𝐺𝑘 + 𝐺𝑏 − 𝜌𝜀 − 𝑌𝑀 + 𝑆𝑘 (6) 𝜕(𝜌𝜀) 𝜕𝑡 + ∇(𝜌𝜀𝑢) = ∇(𝛼𝜀 𝜇𝑒𝑓𝑓 ∇𝑘) + 𝐶1𝜀 𝜀 𝑘 (𝐺𝑘 + 𝐶3𝜀 𝐺𝑏 ) − 𝐶2𝜀 𝜌 𝜀2 𝑘 − 𝑅𝜀 + 𝑆𝜀 (7) In equation (6) and (7) 𝐺𝑘 = Generation of turbulence kinetic energy due to the mean velocity gradients 𝐺𝑏 = Generation of turbulence kinetic energy due to buoyancy 𝑌𝑀 = Contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate 𝑆𝑘 = Source term for turbulence kinetic energy 𝑆𝜀 = Source term for turbulence kinetic energy dissipation rate The erosion model is included in DPM. Erosion rate is defined by following equation: 𝑅𝑒𝑟𝑜𝑠𝑖𝑜𝑛 = ∑ �̇�𝑝𝐶(𝑑𝑝)𝐹(𝑎)𝑣𝑖 𝑛 𝐴𝑓𝑎𝑐𝑒 𝑁𝑝𝑎𝑟𝑡𝑖𝑐𝑙𝑒𝑠 𝑝=1 (8) Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 40 Here, 𝐶(𝑑𝑝) is particle diameter function and 𝐹(𝑎) is impact angle function. Following boundary conditions are applied in the computational domain to solve the governing equations numerically: Inlet: velocity of air at Y-direction, 𝑢𝑦 = 𝑢𝑖𝑛; other velocity components, 𝑢𝑥 = 𝑢𝑧 = 0; temperature of air and sand particles, 𝑇𝑖𝑛 = 50℃; mass flow rate of sand particles, �̇� = �̇�𝑖𝑛. Outlet: pressure, 𝑝 = 𝑝𝑜𝑢𝑡. Atmospheric pressure is defined at domain outlet. Walls: Fluid velocity adjacent to stationary domain walls, 𝑢 = 0. Outer walls are considered thermally insulated. Hence, temperature gradient normal to the wall, 𝜕𝜃 𝜕𝑛 = 0. 3. NUMERICAL SIMULATION 3.1 Meshing and Grid Independent Test A Hybrid mesh has been generated by ANSYS Workbench. Figure 2 shows meshing of turbine blade surface and the fluid domain with 60° impingement angle. The grid consists of hexahedron elements. Minimum number of elements should be considered in the simulation to minimize computational time. To find out optimum number of elements, grid independent test has been carried out. Erosion rate has been observed for 1.3 × 105 to 1.5 × 105 number of elements and relative error has been calculated, which are shown in Table 1. Finally, 144828 number of cells has been selected for simulation with a relative error of 1.36%. 3.2 Solution methods Governing equations have been solved numerically by finite volume method (FVM). Second order upwind scheme has been used for discretization. Pressure-velocity coupling equation has been solved by Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm, developed by Patankar and Spalding [13]. (a) (b) (c) Figure 2: (a) Wind turbine blade surface meshing; Meshing of fluid domain with 60° impingement angle (b) cross- sectional view, (c) three dimensional view. Table 1: Grid independence test results No. of elements Rerrosion (kg.m -2 .s) Relative error (%) 134425 1.31 × 10−7 13.74 137854 1.38 × 10−7 7.97 144828 1.47 × 10−7 1.36 148710 1.44 × 10−7 3.47 153318 1.49 × 10−7 -- Hasan and Ahmed (2017): International Journal of Engineering Materials and Manufacture, 2(3), 37-48 41 4. EXPERIMENTAL STUDY 4.1 Samples and nanocoating Glass fibre 6 X 2 cm samples were collected from the 3 mm thick turbine blades and jute reinforced epoxy resin composite samples were prepared in the laboratory. Multiple acetone wash removed dusts and gel coat of the turbine blades, if any. RF sputtering technique was adopted for the coating while target material was 99.99% Titanium. Gas flow rate for Argon (Ar) was 50 sccm and for Nitrogen (N) was 5 sccm. Duration of each deposition was 60 minutes. 4.2 Sand blasting and roughness analysis Sand blasting machine with a blasting pressure of 50 psi (20 m/sec) and 75 psi (35 m/sec) were used to erosion rate calculation. The test rig is shown in Figure 3. Target samples were kept 20cm away from the nozzle tip and the duration of sand blasting was 10 seconds for each sample. Brown aluminium oxide of 0.30 to 0.70 mm particles were used as sand particles. Sand blasted samples were analysed by Mitutoyo Surftest (SJ-301) surface roughness analyser. JIS B0601(2001) method was adopted and Gauss filter was used to calculate the roughness (Ra) value. 5. RESULTS AND DISCUSSION 5.1 Simulation results 5.1.1 Erosion rate Figure 4 shows the contour plot of erosion on wind turbine blade surface of glass fibre resin for 60 ms−1 wind velocity and 90° impingement angle. It is observed that highest erosion rate on turbine blade surface is 1.4 × 10−6 kgm−2s−1. In Figure 5, highest erosion rate on wind turbine blade surface of glass fibre resin is plotted for different velocities and impingement angles. It is found that erosion rate elevates with increasing wind velocity for all angles. At low wind velocities, erosion rate is highest for 30° impingement angle, which is followed by 90°, 60°, 45° and 20° angles respectively. However, erosion rate at high velocities for 90° and 30° angles are same. Exceptionally, at high velocities, erosion rate for 45° angle is higher than 60° impingement angle. At 65 ms−1 wind velocity, erosion rate for 45° angle is 8 × 10−7 kgm−2s−1, while for 60° angle it is 6 × 10−7 kgm−2s−1. In Figure 6, erosion rate for 30° impingement angle is plotted for jute fibre composite and 30-70% glass fibre resin. It shows clearly that erosion rate on jute fibre composite blade is higher than glass fibre resin. 5.1.2 Wall shear stress In Figure 7, development of shear stress on wind turbine blade surface is shown by contour plot. It is visible that shear stress is very high at the edge of the blade. Hence, the possible damage from wind flow may occur at the blade edge. Similar to erosion rate, shear stress on the blade surface also elevates with increasing wind velocity for all angles, which is observed in Figure 8. Wall shear stress for 20° impingement angle is much higher than other angles. Exceptionally, wall shear stress for 90° angle at 65 ms−1 wind velocity is higher than 30° and 60° angles. However, at other velocities, shear stresses for 30° and 60° angles are found higher than 90° angle. At 40 ms−1 wind velocity, shear stress for 60° impingement angle is about 50 Pa which is much higher than 30° angle. However, at higher velocities, shear stresses for 30° and 60° angles are found almost same. Lowest shear stress is developed for 45° angle at all velocities. Shear stress is also higher on jute fibre composite blade than glass fibre resin as shown in Figure 9. Figure 3: Test facility for sand blasting on wind turbine blade prototype Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 42 Figure 4: Erosion on wind turbine blade surface of glass fibre-resin for 60 ms−1 wind velocity and 90° impingement angle Figure 5: Erosion rate vs. wind velocity curve for glass fibre resin at different angle of attacks Hasan and Ahmed (2017): International Journal of Engineering Materials and Manufacture, 2(3), 37-48 43 Figure 6: Rate of erosion on two different composite materials for 30° angle of attack Figure 7: Shear stress on wind turbine blade surface of glass fibre-resin for 65 ms−1 wind velocity and sand injection from 90° angle Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 44 Figure 8: Wall shear stress vs. wind velocity curve for glass fibre resin at different angle of attacks Figure 9: Development of shear stress on two different composite materials for 30° angle of attack 5.1.3 Strain rate In Figure 10, alike wall shear stress, strain rate is also high at the edge of the blade. Again, from Figure 11, it is seen that strain rate rises with the increment of wind velocity. Strain rate is found highest for 60° angle of attack at all velocities, which is followed by 30°,20°, 45° and 90° angles respectively. However, at 65 ms−1 wind velocity, same strain rate is observed for 45° and 90° angles. Hasan and Ahmed (2017): International Journal of Engineering Materials and Manufacture, 2(3), 37-48 45 5.1.4 Pressure distribution In Figure 12, pressure distribution on 30-70% glass fibre resin blade surface are shown for various angle of attacks has been plotted. It is perceived that highest pressure developed on blade surface declines with increasing angle. However, for 90° angle of attack, highest pressure developed on the blade suddenly rises to 29.5 KPa, which is larger than pressures developed at all other angles. Figure 10: Strain rate on wind turbine blade surface of glass fibre-resin for 65 ms−1 wind velocity and sand injection from 60° angle Figure 11: Strain rate vs. wind velocity curve for glass fibre resin Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 46 Figure 12: Highest pressure vs. injection angle curve for glass fibre resin, 35 ms−1 wind velocity, 0.3 mm diameter sand particle injection at 0.5 kgs−1 mass flow rate. 5.2 Experimental results 5.2.1 Comparison with simulation results In Figure 13, both experimental and simulation results of erosion rate for 5⁰, 30⁰ and 60⁰ impingement angles have been plotted. Simulation results show good agreement with experimental data for 5⁰ and 60⁰ angles. However, for 30⁰ angle, erosion rate is observed higher in experiment. The mean absolute error (MAE) is defined using following equation: 𝑀𝐴𝐸 = 1 𝑀 ∑ |𝑅𝑒𝑟𝑜𝑠𝑖𝑜𝑛,𝑝𝑟𝑒𝑑−𝑅𝑒𝑟𝑜𝑠𝑖𝑜𝑛,𝑒𝑥𝑝| 𝑅𝑒𝑟𝑜𝑠𝑖𝑜𝑛,𝑒𝑥𝑝 × 100% (9) Here, M represents number of data points. Calculated MAE for 30° impingement angle is 5.56%. Figure 13: Comparison of experimental and simulation results for erosion rate at different angle of attacks Hasan and Ahmed (2017): International Journal of Engineering Materials and Manufacture, 2(3), 37-48 47 5.2.2 Surface roughness Surface roughness increases almost linearly with increasing impingement angle, which can be apprehended from Figure 14. Surface roughness on uncoated blade surface collided with 0.3-0.69 mm diameter particles is highest, which is expected. Surface roughness is found low for finer particle sizes. Figure 14: Experimental results of surface roughness vs. angle of attacks for different coated and uncoated surfaces 6. Conclusions Rate of Erosion, development of shear stress and strain rate on wind turbine blade surface in desert conditions have been explored. Pressure distribution on wind turbine blade for various injection angles has also been observed. Major findings from numerical and experimental investigations are following: 1. In both simulation and experiment, erosion rate from wind turbine blade surface is found highest for 30° impingement angle. 2. Shear stress development on wind turbine blade surface is highest for 20° angle and lowest for 45° angle. Although strain rate has been found lowest for 90° impingement angle, it is also very low for 45° angle. 3. Erosion rate on glass fibre resin has been found lower than jute fibre composite. Hence, it is more suitable as a wind turbine blade material from erosion perspective. 4. Pressure distribution on turbine blade surface is highest for 90⁰ angle of attack. 5. Surface roughness increases linearly with increment of impingement angle. As expected, surface roughness observed on uncoated blade surface is higher than coated surfaces. ACKNOWLEDGEMENT This research was funded by Ministry of Higher Education, Malaysia under Research Grant FRGS 14-131-0372. The authors are grateful to those who contributed directly and indirectly in producing this paper. References 1. Mishnaevsky, L., Freere, P., Sinha, R., Acharya, P., Shrestha, R., & Manandhar, P. (2011). Small wind turbines with timber blades for developing countries: Materials choice, development, installation and experiences. Renewable Energy, 36(8), 2128-2138. 2. Stachowiak, G.W. and Batchelor, A.W. 2001. Engineering Tribology, 2nd, Boston: Butterworth-Heinemann. Wear Resistance Performance of Conventional and Non-Conventional Wind Turbine Blades with TiN Nano-Coating 48 3. Friedrich, K. 1997. Wear Performance of High Temperature Polymers and their Composites. CRC Press, : 221– 246. 4. Briscoe, B.J. 1993. The tribology of composite materials: a preface. Adv. Compos. Tribol, 8: 3–15. 5. Hager, A.M. and Davies, M. 1993. Short-fibre reinforced, high temperature resistant polymers for a wide field of tribological applications. Adv. Compos. Tribol, 8: 107–157. 6. Kong, C., Bang, J., & Sugiyama, Y. (2005). Structural investigation of composite wind turbine blade considering various load cases and fatigue life. Energy, 30(11), 2101-2114. 7. Joosse, P. A., Van Delft, D. R. V., Kensche, C., Hahn, F., Jacobsen, T., & van den Berg, R. (2000, January). Economic Use of Carbon Fibres in Large Wind Turbine Blades. In Proceedings of AIAA/ASME Wind Energy Symposium. Reno, NV. 8. Bahadur, S. and Gong, D. 1992. The role of copper compounds as fillers in the transfer and wear behavior of polyetheretherketone. Wear, 154: 151–165. 9. Schwartz, C.J. and Bahadur, S. 2001. The role of filler deformability, filler polymer bonding, and counterface material on the tribological behavior of polyphenylene sulfide (PPS). Wear, 251: 1532–1540. 10. Yu, L., Yang, S., Wang, H. and Xue, Q. 2000. An investigation of the friction and wear behaviors of micrometer copper particle-filled polyoxymethylene composites. J. Appl. Polymer Sci, 19: 2404–2410. 11. Bahadur, S. and Kapoor, A. 1992. The effect of ZnF2, ZnS and PbS fillers on the tribological behavior of nylon 11. Wear, 155: 49–61. 12. A. Patnaik, T. G. Mamatha, S. Biswas, & P. Kumar, (2012). Damage assessment of titania filled zinc–aluminum alloy metal matrix composites in erosive environment: A comparative study. Materials & Design, 36 (2012), 511- 521. 13. Drensky, G., Hamed, A., Tabakoff, W., & Abot, J. (2011). Experimental investigation of polymer matrix reinforced composite erosion characteristics. Wear, 270(3), 146-151. 14. Fiore, G., & Selig, M. S. A Simulation of Operational Damage for Wind Turbine Blades. 15. N. Gaudern, A practical study of the aerodynamic impact of wind turbine blade leading edge erosion. Journal of Physics: Conference Series, 524(1) (2014). 16. A. A. Hamed, W. Tabakoff, R. B. Rivir, K. Das & P. Arora. Turbine blade surface deterioration by erosion. Journal of Turbomachinery, 127(3) (2005), 445-452. 17. Foley, T., & Levy, A. (1983). The erosion of heat-treated steels. Wear, 91(1), 45-64. 18. Hutchings, I. M. (1992). Tribology: friction and wear of engineering materials. 19. Rodríguez, E., Flores, M., Pérez, A., Mercado-Solis, R. D., González, R., Rodríguez, J., & Valtierra, S. (2009). Erosive wear by silica sand on AISI H13 and 4140 steels. Wear, 267(11), 2109-2115. 20. Rico, A., Poza, P., & Rodríguez, J. (2013). High temperature tribological behavior of nanostructured and conventional plasma sprayed alumina-titania coatings. Vacuum, 88, 149-154. 21. Liang, F., Gou, J., Kapat, J., Gu, H., & Song, G. (2011). Multifunctional nanocomposite coating for wind turbine blades. International Journal of Smart and Nano Materials, 2(3), 120-133. 22. Zhang, H., Tang, L., Zhang, Z., Gu, L., Xu, Y., & Eger, C. (2010). Wear-resistant and transparent acrylate-based coating with highly filled nanosilica particles. Tribology International, 43(1), 83-91. 23. B. E. Launder and D. B. Spalding. Lectures in Mathematical Models of Turbulence. Academic Press, London, England. 1972. 24. S. V. Patankar and D. B. Spalding, A calculation procedure for heat, mass and momentum transfer in three- dimensional parabolic flows. International Journal of Heat and Mass Transfer, 15(10) (1972), 1787-1806.