Simulation of different shapes and arrangements of holes over the leading edge of airfoil by blowing to prevent ice accr
- PDF / 5,353,289 Bytes
- 14 Pages / 595.276 x 790.866 pts Page_size
- 6 Downloads / 147 Views
(2020) 42:448
TECHNICAL PAPER
Simulation of different shapes and arrangements of holes over the leading edge of airfoil by blowing to prevent ice accretion Yasser Barzanouni1 · Mofid Gorji‑Bandpy1 · Hassan Basirat Tabrizi2 Received: 17 June 2019 / Accepted: 25 July 2020 © The Brazilian Society of Mechanical Sciences and Engineering 2020
Abstract Effect of blowing out on NACA0012 airfoil surface in order to prevent of ice accretion was studied numerically. Holes with regular pattern were generated over the leading edge of airfoil. Velocity of air injection, diameter, pitch of holes, and symmetrical angle of two rows about of cord line and arrangement of holes, diamond, and rectangle were varied. Shear-stress transport k-ω model was chosen to simulate the turbulence closure model, which makes better prediction for the case of pressure gradient over airfoils. Protective layer of air injections caused to prevent water droplets from the impact and strike to the airfoil surface. Results indicated that the diameter has the most effect on lowering of the ice accretion. In addition, pitch of the holes has the second most important role for reduction of ice weight. Reduction of 89% in the ice accretion in the presence of blowing out was achieved. Keywords Airfoil · Anti-icing · Blowing · CFD · Ice accretion · SST-k-ω List of Symbols AOA Angle of attack Cice Specific heat of ice (J/(kg.K)) Cw Specific heat of water (J/(kg.K)) Cd Drag coefficient CL Lift coefficient d Cylinder diameter or twice the leading-edge radius of airfoil, m Fr Froude number, dimensionless Gk Generation of turbulence kinetic energy due to mean velocity gradients Gw Generation of ω hc Convective heat-transfer coefficient, W/m2 K hf Film thickness (m) kt Turbulent kinetic energy (m2/s2) K Inertia parameter, dimensionless L∞ Characteristic length (m) Lsubl Latent heat of sublimation (J/kg) Technical Editor: André Cavalieri. * Yasser Barzanouni [email protected]; [email protected] 1
Babol Noshirvani University of Technology (NIT), Babol, Mazandaran 47148‑71167, Islamic Republic of Iran
Amirkabir University of Technology (AUT), Tehran, Tehran 15875‑4413, Islamic Republic of Iran
2
Lfusion Latent heat of freezing of water, (J/kg) Levap Latent heat of evaporation of water, (J/kg) LWC Liquid–water content (g/m3) MVD Water droplet median volume diameter (μm) ṁ evap Instantaneous mass of evaporation ṁ ice Instantaneous mass of ice accretion P Static pressure (N/m2) Re Reynolds number, dimensionless T Absolute temperature (K) V Velocity (m/s) V Velocity vector, dimensionless Y Dissipation due to turbulence Y* Dimensionless coordinate t Time α Volume fraction β Collection efficiency, dimensionless Γ Effective diffusivity δ* Dimensionless ice thickness (= δ/cord) ε Solid emissivity or, turbulent dissipation rate (m2/ s 3) θ Angle with respect to cord line μ Viscosity of air (g/m s) ρ Density (kg/m3) σ Boltzman constant (= 5.67 × 10−8 W/(m2 K4)) τ Shear stress tensor (N/m2) ϕ Holes diameter (mm)
Data Loading...
