• Energy Research

Dynamic stall on airfoils negatively impacts their aerodynamic performance and can lead to structural damage. Accurate prediction and modelling of the dynamic stall loads are crucial for a more robust design of wings and blades that operate under unsteady conditions susceptible to dynamic stall and for widening the range of operation of these lifting surfaces. Many dynamic stall models rely on empirical parameters that need to be obtained from experimental or numerical data which limits their generalisability. Here, we introduce physically derived time scales to replace the empirical parameters in the Goman–Khrabrov dynamic stall model. The physics-based time constants correspond to the dynamic stall delay and the decay of post-stall load fluctuations. The dynamic stall delay is largely independent of the type of motion, the Reynolds number and the airfoil geometry, and is described as a function of a normalised instantaneous pitch rate. The post-stall decay is independent of the motion kinematics and is related to the Strouhal number of the post-stall vortex shedding. The general validity of our physics-based time constants is demonstrated using three sets of experimental dynamic stall data covering various airfoil profiles, Reynolds numbers varying from 75 000 to 1 000 000, and sinusoidal and ramp-up pitching motions. The use of physics-based time constants generalises the Goman–Khrabrov dynamic stall model and extends its range of application.

Dynamic stall on airfoils negatively impacts their aerodynamic performance and can lead to structural damage. Accurate prediction and modelling of the dynamic stall loads are crucial for a more robust design of wings and blades that operate under unsteady conditions susceptible to dynamic stall and for widening the range of operation of these lifting surfaces. Many dynamic stall models rely on empirical parameters that need to be obtained from experimental or numerical data which limits their generalisability. Here, we introduce physically derived time scales to replace the empirical parameters in the Goman–Khrabrov dynamic stall model. The physics-based time constants correspond to the dynamic stall delay and the decay of post-stall load fluctuations. The dynamic stall delay is largely independent of the type of motion, the Reynolds number and the airfoil geometry, and is described as a function of a normalised instantaneous pitch rate. The post-stall decay is independent of the motion kinematics and is related to the Strouhal number of the post-stall vortex shedding. The general validity of our physics-based time constants is demonstrated using three sets of experimental dynamic stall data covering various airfoil profiles, Reynolds numbers varying from 75 000 to 1 000 000, and sinusoidal and ramp-up pitching motions. The use of physics-based time constants generalises the Goman–Khrabrov dynamic stall model and extends its range of application.

Airfoil stall plays a central role in the design of safe and efficient lifting surfaces. We typically distinguish between static and dynamic stall based on the unsteady rate of change of an airfoil's angle of attack. Despite the somewhat misleading denotation, the force and flow development of an airfoil undergoing static stall are highly unsteady and the boundary with dynamic stall is not clearly defined. We experimentally investigate the forces acting on a two-dimensional airfoil that is subjected to two manoeuvres leading to static stall: a slow continuous increase in angle of attack with a reduced pitch rate of 1.3e-4 and a step-wise increase in angle of attack from 14.2�� to 14.8�� within 0.04 convective times. We systematically quantify the stall reaction delay for many repetitions of these two manoeuvres. The onset of flow stall is marked by the distinct drop in the lift coefficient. The reaction delay for the slow continuous ramp-up manoeuvre is not influenced by the blade kinematics and its occurrence histogram is normally distributed around 32 convective times. The static reaction delay is compared with dynamic stall delays for dynamic ramp-up motions with reduced pitch rates ranging from 9e-4 to 0.14 and for dynamic sinusoidal pitching motions of different airfoils at higher Reynolds numbers up to 1e6. The stall delays for all conditions follows the same power law decrease from 32 convective times for the most steady case down to an asymptotic value of 3 for reduced pitch rates above 0.04. Static stall is not phenomenologically different than dynamic stall and is merely a typical case of stall for low pitch rates. Based on our results, we suggest that conventional measurements of the static stall angle and the static load curves should be conducted using a continuous and uniform ramp-up motion at a reduced frequency around 1e-4.

Airfoil stall plays a central role in the design of safe and efficient lifting surfaces. We typically distinguish between static and dynamic stall based on the unsteady rate of change of an airfoil's angle of attack. Despite the somewhat misleading denotation, the force and flow development of an airfoil undergoing static stall are highly unsteady and the boundary with dynamic stall is not clearly defined. We experimentally investigate the forces acting on a two-dimensional airfoil that is subjected to two manoeuvres leading to static stall: a slow continuous increase in angle of attack with a reduced pitch rate of 1.3e-4 and a step-wise increase in angle of attack from 14.2�� to 14.8�� within 0.04 convective times. We systematically quantify the stall reaction delay for many repetitions of these two manoeuvres. The onset of flow stall is marked by the distinct drop in the lift coefficient. The reaction delay for the slow continuous ramp-up manoeuvre is not influenced by the blade kinematics and its occurrence histogram is normally distributed around 32 convective times. The static reaction delay is compared with dynamic stall delays for dynamic ramp-up motions with reduced pitch rates ranging from 9e-4 to 0.14 and for dynamic sinusoidal pitching motions of different airfoils at higher Reynolds numbers up to 1e6. The stall delays for all conditions follows the same power law decrease from 32 convective times for the most steady case down to an asymptotic value of 3 for reduced pitch rates above 0.04. Static stall is not phenomenologically different than dynamic stall and is merely a typical case of stall for low pitch rates. Based on our results, we suggest that conventional measurements of the static stall angle and the static load curves should be conducted using a continuous and uniform ramp-up motion at a reduced frequency around 1e-4.

We propose here a method to experimentally quantify unsteady leading-edge flow separation on aerofoils with finite thickness. The methodology relies on the computation of a leading-edge suction parameter based on measured values of the partial circulation around the leading edge and the stagnation point location. We validate the computation of the leading-edge suction parameter for both numerical and experimental data under steady and unsteady flow conditions. The leading-order approximation of the definition of the leading-edge suction parameter is proven to be sufficiently accurate for the application to thin aerofoils such as the NACA0009 without a priori knowledge of the stagnation point location. The higher-order terms including the stagnation point location are required to reliably compute the leading-edge suction parameter on thicker aerofoils such as the NACA0015. The computation of the leading-edge suction parameter from inviscid flow theory does not assume the Kutta condition to be valid at the trailing edge which allows us to compute its value for separated flows. The relation between the leading-edge suction parameter and the evolution of the shear layer height is studied in two different unsteady flow conditions, a fixed aerofoil in a fluctuating free-stream velocity and a pitching aerofoil in a steady free stream. We demonstrate here that the instantaneous value of the leading-edge suction parameter based on the partial circulation around the leading edge is unambiguously defined for a given flow field and can serve as a directly quantitative measure of the degree of unsteady flow separation at the leading edge.

We propose here a method to experimentally quantify unsteady leading-edge flow separation on aerofoils with finite thickness. The methodology relies on the computation of a leading-edge suction parameter based on measured values of the partial circulation around the leading edge and the stagnation point location. We validate the computation of the leading-edge suction parameter for both numerical and experimental data under steady and unsteady flow conditions. The leading-order approximation of the definition of the leading-edge suction parameter is proven to be sufficiently accurate for the application to thin aerofoils such as the NACA0009 without a priori knowledge of the stagnation point location. The higher-order terms including the stagnation point location are required to reliably compute the leading-edge suction parameter on thicker aerofoils such as the NACA0015. The computation of the leading-edge suction parameter from inviscid flow theory does not assume the Kutta condition to be valid at the trailing edge which allows us to compute its value for separated flows. The relation between the leading-edge suction parameter and the evolution of the shear layer height is studied in two different unsteady flow conditions, a fixed aerofoil in a fluctuating free-stream velocity and a pitching aerofoil in a steady free stream. We demonstrate here that the instantaneous value of the leading-edge suction parameter based on the partial circulation around the leading edge is unambiguously defined for a given flow field and can serve as a directly quantitative measure of the degree of unsteady flow separation at the leading edge.

Abstract Vertical-axis wind turbines are great candidates to enable wind power extraction in urban and off-shore applications. Currently, concerns around turbine efficiency and structural integrity limit their industrial deployment. Flow control can mitigate these concerns. Here, we experimentally demonstrate the potential of individual blade pitching as a control strategy and explain the flow physics that yields the performance enhancement. We perform automated experiments using a scaled-down turbine model coupled to a genetic algorithm optimiser to identify optimal pitching kinematics at on- and off-design operating conditions. Optimal blade kinematics yields a three-fold power coefficient increase at both operating conditions compared to the non-actuated turbine and a 70% reduction in structure-threatening load fluctuations at off-design conditions. Based on flow field measurements, we uncover how blade pitching manipulates the flow structures to enhance performance. Our results can aid vertical-axis wind turbines to increase their much-needed contribution to our energy needs.