---

Blade is a sensitive and important part of turbines, and a few companies can produce it. Airfoil of blade has three-dimensional surfaces; therefore, it is necessary to have specific equipment for dimensional control of it. The purpose of this project is to design and manufacture a mechanical system for dimensional control of the airfoil. The foregoing device can produce two-dimensional contours of airfoil on the screen of the profile projector using fine pins. In the mentioned system, the blade is located on the table of device and two sets of pins approach it. In this situation pins are moved forward along their axis until their tips touch the surface of the blade, therefore the tips of pins, shape the contour of the airfoil. Then two sets of pins get away from each other and are moved to the focal area of lens of profile projector via a precision linear system. Then two sets of pins approach each other and reshape the previous contour again. In this situation, tips of pins are projected and contour of airfoil is made on the screen of the profile projector with a predefined scale. Produced contours can be compared with reference ones that have been printed on transparent sheets.

---

Blade is a sensitive and important part of turbines, and a few companies can produce it. Airfoil of blade has three-dimensional surfaces; therefore, it is necessary to have specific equipment for dimensional control of it. The purpose of this project is to design and manufacture a mechanical system for dimensional control of the airfoil. The foregoing device can produce two-dimensional contours of airfoil on the screen of the profile projector using fine pins. In the mentioned system, the blade is located on the table of device and two sets of pins approach it. In this situation pins are moved forward along their axis until their tips touch the surface of the blade, therefore the tips of pins, shape the contour of the airfoil. Then two sets of pins get away from each other and are moved to the focal area of lens of profile projector via a precision linear system. Then two sets of pins approach each other and reshape the previous contour again. In this situation, tips of pins are projected and contour of airfoil is made on the screen of the profile projector with a predefined scale. Produced contours can be compared with reference ones that have been printed on transparent sheets.

---

This paper deals with the effects of duty cycle on improvement of pressure distribution over a NLF0414 airfoil using the plasma actuators. Three Dielectric barrier discharges as the plasma actuators are flush mounted on the airfoil surface in different positions to improve pressure distribution at post-stall angles of attack. The experiments were performed in wind tunnel with pressure tabs measurements at Re_c=7.5×〖10〗^5. The main objective of these experiments is to find the most effective duty cycles for different excitation frequencies and different angles of attack. Results show that the plasma actuators in unsteady excitations are more effective in lower duty cycles on low excitation frequencies but the lower duty cycles lose their effectiveness on higher excitation frequencies.

---

In this study a new approach for investigating the flutter speed of nonlinear aeroelastic systems is proposed. In this approach, the compatibility of nonlinear random vibration analysis based on the statistical properties of response is used and extended to the nonlinear aeroelastic systems to analyze the instability of these systems with using neither time domain analysis nor limit cycle oscillations. To this aim a 2-degree nonlinear airfoil with cubic torsional spring under quasi steady flow is considered as an aeroelastic system. At first, one random Gaussian white noise is added to the aerodynamic lift force then the statistical linearization and the random vibration analysis of the nonlinear systems are used to obtain a nonlinear map of response-variance with flow velocity as the control parameter. This nonlinear map leads to a nonlinear algebraic equation which consists of two parameters as the flow velocity and variance of the response. By solving this nonlinear equation for various flow velocities, the flutter speed is considered as the maximum of response-variance. Finally the jump phenomenon is also investigated where tangent bifurcation occurs.

---

In this study, longitudinal dynamic derivatives of an airfoil of the type NACA 6-series, oscillating in pitching and plunging motions were calculated using variation of pitching moment coefficients with angle of attack in various conditions, based on wind tunnel data. Various parameters of the tests were mean angle of attack, reduced frequency and amplitude of oscillation. To calculate the longitudinal dynamic derivatives in harmonic oscillations, the Taylor's series and integral of Fourier were used. Both the methods had the same results and could be extended to each flight vehicles. The effect of parameters on variation of longitudinal oscillatory derivatives was investigated, in three different regions of oscillation: before, over and post stall conditions. The results showed that variation of the longitudinal oscillatory coefficients with angle of attack is different in the pre-stall and over stall conditions with respect to post-stall region. The effect of reduced frequency on stability of the motion is different for two types of oscillations. Increasing the reduced frequency resulted in reducing the stability of plunging motion, but has a little effect on the stability of pitching motion.

---

In this paper an improved immersed boundary method is used for simulating sinusoidal pitching oscillations of a symmetric airfoil. Immersed boundary methods because of using a fixed Cartsian grid are well suited for such moving boundary problems. Two test cases are used to validate the proposed method and the effects of oscillation frequency and amplitude on the flow field are investigated. Flow field vorticity and kinetic energy contours are reported in this paper. It is found that the deflected wake start to be appeared for Strouhal number more than 0.4 at a fixed pitching amplitude 0.71. A chaotic flow can be observed at oscillation amplitude 2.80, for a fixed Strouhal number, 0.22. Kintic energy contour shows that for Strouhal number 0.1, the airfoil performs work and transfers momentum to flow but the fluid energy loss due to the enlargement of flow separation zone decreases the momentum and kinetic energy behind the airfoil. Deficit momentum and kinetic energy behind the airfoil results in drag force increasing. By increasing the oscillation frequency and amplitude more momentum transfers to flow filed behind the airfoil which results in drag force decreasing.

---

In this research, the plunging motion of an airfoil by a numerical method based on finite volume in different Reynolds numbers is simulated and the thickness effect, amplitude and reduced frequency on the aerodynamic coefficients are investigated. In this process, SIMPLEC algorithm, implicit solver, high order scheme and dynamic mesh method is used in unsteady simulation and the flow is supposed viscous, incompressible and laminar. Simulations are in three Reynolds 1000, 11000 and 50000, respectively, in accordance with the flight of the insects, small birds and pigeons are done in two amplitudes and three reduced frequencies. The simulation results are compared with published data to confirm the validity of research. This comparison shows comprisable agreement. Pressure distribution and Vortex shedding around airfoils show that the thickness of the airfoils delays vortex shedding and changes time-averaged thrust coefficient. Reduced frequency and amplitude of oscillation are two important parameters in this simulation, but the reduce frequency is more effective than amplitude. The response surface methodology (RSM) was used to optimize the plunging airfoil. Optimization shows that airfoil with 0.29% thickness, 3.08 reduced frequency and 0.5 dimensionless oscillating amplitude produce maximum trust coefficient.

---

In this paper, Effects of runback ice accretion on NACA 23012 airfoil have been studied experimentally and numerically. For this purpose, experiments were applied on runback ice within Reynolds number of 0.6×〖10〗^6 over the angle of attack from 0 degree to 20 degree and then results were compared with the results of clean airfoil. Generally, Having examined behavior of the flow pattern and aerodynamic coefficients of the iced airfoil the results of which were compared to that of the clean airfoil, it can be concluded that icing phenomenon affects aerodynamic performance of the airfoil in two ways; in the first way that occurs at low angles of attack prior to stalling of the airfoil the effect is local .In this case ice accretion on the airfoil contributes to formation of a flow separation bubble behind the ice ridge on the upper surface of the airfoil. After numerical simulation of flow field, flow separation bubble behind the ice ridge was observed. The main effect of icing which is related to the second way occurs at angles of attack close to stall and post-stall. In this case flow pattern around the airfoil as well as aerodynamic coefficients undergo a fundamental change. In addition, it was made clear that runback ice causes stall angle decreases 2 degree and maximum lift reduces about 8 percent.

---

The application of wing and stabilizer in aerospace vehicle is most important to stability and flight motion. Nonlinear 2D wing is estimated. Nonlinear damping and stiffness with freeplay in plunging and pitching motion is assumed. 2nd order Damping nonlinearity and 3rd order stiffness nonlinearity in pitching and plunging motion is assumed. Fully nonlinear structure with nonlinear 3rd order piston theory aerodynamic is assumed for the first time and result evaluated with different references. The equations are defined with Hamilton principle with the use of kinetic and potential energy and virtual work. They are solved in the state space via the ruge-kuta numerical method to determine chaotic and limit cycle oscillation motion of supersonic airfoil. The result show that as the speed increases, the behavior of 2D wing is softening type with the use of nonlinear rotational stiffness. But, It shows hardening type with the use of transversal nonlinear stiffness. The effect of transversal and rotational freeplay is more complicated than other parameters and increases instability in low speed. In other hand the stability increases with freeplay in high speed. As shown, increase velocity decrease damping effect in post flutter behavior.

---

This article investigates experimental study of the flow field on a blunt airfoil. For this purpose, PIV technique based on instantaneous flow structures is used in order to view and two dimensional investigation of flow field around unmodified and blunt airfoil and at different times. This study is performed on flows at very low Reynolds number(Reynolds number lower than 4500). This flow regime is very similar to dominant condition on micro air vehicles (MAVs). In order to validate the method used in this study, flow field around cylinder is considered and in continue, instantaneous and mean velocities fields, streamlines and mean vortices field around unmodified and blunt airfoils are obtained. The results show that there are prominent differences on the structure of wake around airfoils and sizes of separation region for blunt and simple airfoils. Meanwhile separation of the flow for both blunt and simple airfoils at this very low Reynolds number, is occurred at angle of attack 5 (at low angle of attack). Also generation of vortex at wake region and their position and circulation at different times, are discussed.

---

Utilization of hydro- power as renewable energy source is developed in the world now significantly. Using very low head axial turbines in rivers is one of ways for obtain this energy. In this research, design and optimization of an axial hydraulic turbine with very low head(2.9m) was done. The first step in the optimization of turbine is generation a suitable initial geometry. For this purpose one dimensional design approach based on Euler law was used. Development of computation algorithms is very efficient and suitable in hydraulic turbine design and performance investigation. In this research mesh was generated with mesh-ANSYS software and then the default domain was simulated by solving the 3-D Navier Stokes equations through the runner passages in the CFX software. Optimization geometry is obtained by optimization of Drag to Lift coefficient ratio for different blade midspans. For parameterization of airfoils, the “CST” method and for extraction of flow characteristics of airfoils XFOIL software were utilized. Then airfoil coefficients by minimization of Drag to Lift ratio with fminsearch algorithm in MATLAB software were corrected. The results show that the efficiency in design point is increased about 2.4%.

---

In this study, the effects of Co-Flow jet and injection temperature on the enhancement of airfoil performance in the compressible flow are investigated numerically. Co-Flow jet is a method of increasing lift to drag ratio and varying the Stall Degree which works via injecting the air from the edge of airfoil and suction from the tail. The much number of studied flow changes from 0.4 to 0.6. Clark-Y airfoil has been chosen for this study because of its application in compressible flow, it is the base airfoil for development of new airfoils. A validation is performed for Clark-Y airfoil by comparing the present numerical result and available experimental data in the literature. Results indicate that the enhancement induced by the Co-Flow jet on the compressible flow is less than one in the incompressible flow. The drag and lift coefficients reduces and increases by increasing the jet momentum coefficient, respectively. Using the Co-Flow Jet increase the stall degree. The maximum of lift decrement and drag increment occurs around the stall degree. Increasing the temperature increases lift coefficient slightly where it seems to be better choice in comparison with increment of Jet momentum coefficient due to ease of operation.

---

When a Horizontal axis wind turbine works under yaw condition, each blade element can be considered as an oscillating pitch airfoil while the free stream velocity oscillates horizontally. The unsteady free stream velocity, which is usually ignored, oscillates with the same frequency as the airfoil oscillations and has a great impact on the periodic forces produced by the airfoil oscillation. In order to study the effects of unsteady free stream
velocity on the aerodynamic loads, a 2D NACA0012 oscillating airfoil at Reynolds number of 135000 has been simulated. In this simulation, reduced frequency, reduced amplitude and the phase difference between the free stream velocity oscillation and the airfoil angle of attack oscillation are 0.1≤k≤0.25 ، 0.2≤λ≤0.8 و ϕ=0 ,π, respectively. Results show that free stream oscillations affect the aerodynamic loads, vortex strengths
and dynamic stall characteristics. The lift force can be increased by more than 7 times than that of static case and 3 times compared to the load from steady free stream velocity. Depending on 𝜙 value, the dynamic stall angle of attack can be advanced 1 degree or delayed by more than 7 degrees by increase of reduced amplitude. Also, increase of k always causes delay in leading edge vortex formation and consequently delay in dynamic stall occurrence.

---

Dynamic stal behavour of a NACA0012 airfoil undergoing pitching motion has been studied by a numerical approach. The turbulence intensity, oscillation frequency and amplitude and the Reynolds number were found to be the major contributors in dynamic stall. The flowfield structure and the associated vortices for this airfoil as well as the impact of the oscillation frequency on aerodynamic efficiency were also studied. The simulations were two dimensinal and the k-ω SST turbulence model were utilized for the present analysis. The results show that increasing the oscillation frequency and amplitude and the turbulence intensity, postpones the dynamic stall to higher angles of attack. Furthermore, as increasing the Reynolds number, both the lift coefficient and the width of the associated hysteresis loop decrease. The airfoil aerodynamic efficiency variation with oscillation frequncy has been shown to have a maximum point for all angles of attack considered. The flowfield structure revealed that the main cause of the dynamic stall is a series of low pressure vortices formed at the leading edge which shed into downstream and separate from the surface. A secondary vortex will then appear and increases the lift coefficient dramatically. The present simulation results are in a good agreement whith those found in the literature.

---

In this research, pitching and plunging motion of bio inspired and NACA airfoil are simulated numerically and the effects of reduced frequency, pitching and plunging amplitude on aerodynamic coefficients, power-extraction and propulsion efficiency are investigated and compared with each other. The simulation is done at Reynolds number of 1100 which is correspond to insect flight regime, using dynamic mesh capability of OpenFoam and fluid flow is assumed unsteady, viscous and laminar. In order to computation of fluid flow field, control volume method is used and value of variables store at the center of control volume. Reduced frequency, plunging and pitching amplitudes vary between 0.05-0.5, 0.25-1.75, 15-75 respectively and phase difference between pitching and plunging motion is kept constant at 90 degrees. Comparison of result with published data confirms the validation of research. Combination of different motion parameter such as reduced frequency, pitching and plunging amplitudes determine that bio inspired airfoil acting in power-extraction (fluid works on the airfoil), propulsion (airfoil works on the fluid) or feather (no producing power or propulsion) regime, and qualitatively is the same as NACA airfoil. The obtained results shows that with variation of reduced frequency, pitching and plunging amplitudes, whatever close to the feathering regime, bio inspired airfoil shows higher efficiency than NACA airfoil and vice versa.

---

In this paper, the Semi-Empirical and numerical methods that can be used to investigate the effects of dynamic stall are compared with each other, and the capabilities of the methods are studied. The experimental measurements have been used in order to compare the methods. The Semi-Empirical Leishman-Beddoes (L-B), Snel and ONERA methods have been used, and the finite volume method was being used for numerical simulations. The lift coefficient was being calculated by all the methods at various conditions, and the drag coefficient had been computed by the numerical and Leishman-Beddoes methods. The parameters that have been used in order to compare the methods, are the maximum lift coefficient value, the angle of attack of the largest lift coefficient, the error at upstroke phase and the error at down stroke phase. The results show among the semi-empirical models; the L-B method has the highest precision to predict the lift coefficient, and although the numerical method can investigate the flow with more details, but the error percentage at the down stroke phase is higher than expectations. The results from the drag coefficient modeling show that the numerical method can predict this coefficient better than the L-B method. The results also can help other researchers to select the best dynamic stall model in order to investigate the wind-turbine aerodynamics.

---

Aerodynamic and optimal design of a blade of a horizontal axis wind turbine (HAWT) has been performed in order to extract maximum power output with considering the strength of the blade structure resulted from different loads and moments. A design procedure is developed based on the Blade Element Momentum (BEM) theory and suitable correction factors are implemented to include three-dimensionality effects on the turbine performance. The design process has been modified to achieve the maximum power by searching an optimal chord distribution along the blade. Based on the aerodynamic design, the blade loads have been extracted and the blade mechanical strength has been investigated by analyzing the thickness of the blade surface and the blade material. The developed numerical model can be considered as a suitable tool for aerodynamically and mechanically design of a turbine blade. The results for a 500 W turbine show that the turbine performance improves by 5% approximately, by modifying chord radial distribution. Yield stress analysis shows the effect of introduced chord distribution on the blade strength, in different blade thicknesses and different blade materials. In addition, optimum tip speed ratio for having favorable mechanical safety factor is derived. Three different airfoil are examined for this investigation and comparing their mechanical safety factor.

---

---

In this study, we derived the rotating airfoil system of equation considering Loewy aerodynamics. To this end, we define the local coordinate system on airfoil and reference coordinate on the hub. We define the free air velocity vector and the airfoil rotating speed vector according to the reference coordinate. So, the Kinetic and Potential energies are derived based on linear stiffness and linear damping according to the Hamiltonian principle. Wakes behind the rotating blades form into the helix. Therefore, we the equation of motion with Loewy aerodynamic which compensates the wake effects. Stability analysis is performed by the well-known P-K method. Flutter speed and stability boundary are estimated. Comparing the results of stability analysis and the reference validates the applied method. Furthermore, we proposed the PID Control to suppress the flutter speed. the PID controller input and command. The desired time and error tolerance are selected to design PID controller. Unit step response shows that pitch angle response is under-damped. However, step response tracks input well. Besides, disturbance rejection by considering the gain from input to output to remain below the gain value is analyzed. 

---

A wind turbine airfoil was analysed, using computational fluid dynamics (CFD) to study the oscillating effects and slip boundary conditions. The slip boundary condition is due to applying superhydrophobic surface. Fluids on these surfaces are repelled. The superhydrophobic surface can delay the icing on blades. The surfaces is assumed at the leading edge; the icing can occur on this region. The chosen oscillation parameters was enough for modelling dynamic stall. The dynamic stall cause a severe loading on the blade. This phenomenon is depicted by two vortices: leading edge vortex and trailing edge vortex. Three reduced frequencies are considered:  in a range of  slip lengths. In this regard, the Transition-SST model is applied for SD7037 airfoil with. The results showed that applying a superhydrophobic surface with low values of the slip length cannot be appropriate during the oscillating motion; but at the slip lengths larger than 100 microns, the aerodynamic coefficients are significantly changed. At the highest reduced frequency, the lift and drag coefficients are reduced about 12% and 40%, respectively. Increasing the slip length postponed the vortex formation and stall angle.