Blade Element Momentum Theory

Similar documents
Aerodynamic Performance 1. Figure 1: Flowfield of a Wind Turbine and Actuator disc. Table 1: Properties of the actuator disk.

Rotor reference axis

CHAPTER 3 ANALYSIS OF NACA 4 SERIES AIRFOILS

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH

Calculation of Wind Turbine Geometrical Angles Using Unsteady Blade Element Momentum (BEM)

Multidisciplinary Design Optimization of a Heavy-Lift Helicopter Rotor in Hover

Multidisciplinary Design Optimization Of A Helicopter Rotor Blade

Analysis of the high Reynolds number 2D tests on a wind turbine airfoil performed at two different wind tunnels

Validation of Chaviaro Poulos and Hansen Stall Delay Model in the Case of Horizontal Axis Wind Turbine Operating in Yaw Conditions

Definitions. Temperature: Property of the atmosphere (τ). Function of altitude. Pressure: Property of the atmosphere (p). Function of altitude.

AERODYNAMIC ANALYSIS OF THE HELICOPTER ROTOR USING THE TIME-DOMAIN PANEL METHOD

GyroRotor program : user manual

Verification of a Modelica Helicopter Rotor Model Using Blade Element Theory

Design of Propeller Blades For High Altitude

VORTEX METHOD APPLICATION FOR AERODYNAMIC LOADS ON ROTOR BLADES

Wings and Bodies in Compressible Flows

Propeller theories. Blade element theory

DESIGN AND ANALYSIS METHODS FOR UAV ROTOR BLADES

Development of a Simulation and Optimization Framework for Improved Aerodynamic Performance of R/C Helicopter Rotor Blades

Small-Scale Propellers Operating in the Vortex Ring State

A Compact, Closed Form Solution for the Optimum, Ideal Wind Turbines

Chapter three. Two-dimensional Cascades. Laith Batarseh

Generally, there exists an optimum tip-speed-ratio, λ that maximized C p. The exact λ depends on the individual wind turbine design

Aerodynamic Rotor Model for Unsteady Flow and Wake Impact

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

INCLUSION OF A SIMPLE DYNAMIC INFLOW MODEL IN THE BLADE ELEMENT MOMENTUM THEORY FOR WIND TURBINE APPLICATION

Random Problems. Problem 1 (30 pts)

Research on Propeller Characteristics of Tip Induced Loss

Actuator Surface Model for Wind Turbine Flow Computations

Ideal Optimum Performance of Propellers, Lifting Rotors and Wind Turbines

A case study on optimum tip speed ratio and pitch angle laws for wind turbine rotors operating in yawed conditions

Small-scale Aeroelastic Rotor Simulation, Design and Fabrication

Performance Investigation of Ducted Aerodynamic Propulsors

Air Loads. Airfoil Geometry. Upper surface. Lower surface

Performance Evaluation of a Flexible Rotor in Extreme Ground Effect

Effect of Geometric Uncertainties on the Aerodynamic Characteristic of Offshore Wind Turbine Blades

Aerodynamic Performance Assessment of Wind Turbine Composite Blades Using Corrected Blade Element Momentum Method

Mechanical Engineering for Renewable Energy Systems. Dr. Digby Symons. Wind Turbine Blade Design

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

Assessment of low-order theories for analysis and design of shrouded wind turbines using CFD

Notation. a 0, a 1. c fm c D c L c M. α κ η λ ( c M / η)/(c fm / η)

Numerical Study on Performance of Curved Wind Turbine Blade for Loads Reduction

Assessment of Comprehensive Analysis Calculation of Airloads on Helicopter Rotors

This version was downloaded from Northumbria Research Link:

PERFORMANCE ANALYSIS OF A COAXIAL ROTOR SYSTEM IN HOVER: THREE POINTS OF VIEW

Research on Dynamic Stall and Aerodynamic Characteristics of Wind Turbine 3D Rotational Blade

Propellers and Ducted Fans

The Pennsylvania State University The Graduate School College of Engineering AERODYNAMIC ANALYSIS OF HELICOPTER ROTORS USING

Control Volume Analysis For Wind Turbines

A Numerical Blade Element Approach to Estimating Propeller Flowfields

Wind Turbine Blade Analysis using the Blade Element Momentum Method. Version 1.0

A simplified model for a small propeller with different airfoils along the blade

Sensitivity of Key Parameters in Aerodynamic Wind Turbine Rotor Design on Power and Energy Performance

A LOW ORDER METHOD FOR CO-AXIAL PROPELLER AND ROTOR PERFORMANCE PREDICTION

FINITE-STATE DYNAMIC WAKE INFLOW MODELLING FOR COAXIAL ROTORS

Lecture 4: Wind energy

UNSTEADY AERODYNAMIC ANALYSIS OF HELICOPTER ROTOR BY USING THE TIME-DOMAIN PANEL METHOD

Study on the Mechanism of the Variable-Speed Rotor Affecting Rotor Aerodynamic Performance

Analysis of Wind Turbine Pressure Distribution and 3D Flows Visualization on Rotating Condition

CFD RANS analysis of the rotational effects on the boundary layer of wind turbine blades

Theory of turbo machinery. Chapter 3

Introduction to Aerospace Engineering

Numerical Study on Performance of Innovative Wind Turbine Blade for Load Reduction

The Pennsylvania State University. The Graduate School. College of Engineering. Inviscid Wind-Turbine Analysis Using Distributed Vorticity Elements

ENGR 4011 Resistance & Propulsion of Ships Assignment 4: 2017

ABSTRACT. Professor Inderjit Chopra and Visiting Professor Omri Rand Department of Aerospace Engineering

Description and validation of the rotordisksource class for propeller performance estimation

Active Pitch Control of a Vertical Axis Wind Turbine Enhancing performance in terms of power and loads including dynamic stall effects

Assessment of aerodynamics models for wind turbines aeroelasticity

Actuator disk modeling of the Mexico rotor with OpenFOAM

FINITE ELEMENT ANALYSIS AND OPTIMIZATION DESIGN OF ALUMINIUM AXIAL FAN BLADE

Investigation on the influence of scaling effects in propeller testing through the use of theoretical prediction codes

OPTIMUM AND DESIGN TRENDS OF COMPOUND HELICOPTERS

A Short History of (Wind Turbine) Aerodynamics

The Truth About Elliptic Spanloads or Optimum Spanloads Incorporating Wing Structural Weight

Introduction to Aeronautics

Dynamic Stall of an Experimental Wind Turbine Blade

A SMOKE STUDY OF THE MOTION OF THE VORTICES SHED FROM THE TIPS OF A TWO-BLADED MODEL HELICOPTER ROTOR IN HOVERING FLIGHT A THESIS.

ANALYSIS OF HORIZONTAL AXIS WIND TURBINES WITH LIFTING LINE THEORY

Theoretical Aerodynamic analysis of six airfoils for use on small wind turbines

Wind turbine power curve prediction with consideration of rotational augmentation effects

STRUCTURAL PITCH FOR A PITCH-TO-VANE CONTROLLED WIND TURBINE ROTOR

ν δ - 1 -

DETERMINATION OF HORIZONTAL AXIS WIND TURBINE PERFORMANCE IN YAW BY USE OF SIMPLIFIED VORTEX THEORY

University of California at Berkeley Department of Mechanical Engineering ME 163 ENGINEERING AERODYNAMICS FINAL EXAM, 13TH DECEMBER 2005

MOST of the coupled ap-lag-torsion aeroelastic stability

Aircraft Structures Design Example

PEMP ACD2505. Finite Wing Theory. M.S. Ramaiah School of Advanced Studies, Bengaluru

Rotor design and matching for horizontal axis wind turbines

Semi-Empirical Modeling of Two-Dimensional and Three-Dimensional Dynamic Stall

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Introduction to Aerospace Engineering

Research in Aeroelasticity EFP-2005

Validation and comparison of aerodynamic modelling approaches for wind turbines

Performance. 5. More Aerodynamic Considerations

Insight into Rotational Effects on a Wind Turbine Blade Using Navier Stokes Computations

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay

Robot Dynamics Rotorcrafts: Dynamic Modeling of Rotorcraft & Control

Helicopter Rotor Unsteady Aerodynamics

Transcription:

Blade Element Theory has a number of assumptions. The biggest (and worst) assumption is that the inflow is uniform. In reality, the inflow is non-uniform. It may be shown that uniform inflow yields the lowest induced power consumption. Slide 1

Model y dy Annulus At a distance y Width dy: area da=2πydy Mass flow rate Slide 2

Thrust and Power The thrust for the annulus is : Or in the non-dimensional form: The induced power coefficient is: Slide 3

Thrust and Power For the hover case λ C =0. The thrust coefficient is then: The power coefficient is: Valid for any form of λ Slide 4

Thrust Coefficient Let s assume that λ(r)= λ tip r n for n 0. We can then relate λ tip with C T Slide 5

Power Coefficient Calculating the power coefficient Substituting in the λ tip equation : Since we know that Slide 6

Induced Inflow In the previous equation when n=0 then κ=1, uniform inflow which has the lowest induced power. For n>0 then κ>1. With increasing n the inflow is more heavily biased towards the blade tip. Therefore what we want to achieved is, through the blade properties (pitch angle and chord), the first situation Slide 7

From BET we have found Then we have Which can be expressed: Slide 8

Or in another way: This quadratic equation has the solution: Slide 9

Ideal Twist In hover (λ C =0) it simplifies to Since the objective is uniform inflow then we need θr=const.= θ tip. The ideal twist is given by: Slide 10

Thrust coefficient in Ideal Twist We can then calculate the thrust coefficient with this ideal twist θr= θ tip Recalling that Slide 11

Blade Pitch angle In hover We can then obtain θ tip Slide 12

Final notes Recall that: Therefore the thrust varies linearly with r. We have therefore a linear (triangular) variation of lift over the blade. Slide 13

Ideal Rotor vs. Optimum Rotor Ideal rotor has a non-linear twist:θ=θ tip /r This rotor will, according to the BEM theory, have a uniform inflow, and the lowest induced power possible. The rotor blade will have very high local pitch angles θ near the root, which may cause the rotor to stall. Ideally Twisted rotor is also hard to manufacture. For these reasons, helicopter designers strive for optimum rotors that minimize total power, and maximize figure of merit. This is done by a combination of twist, and taper, and the use of low drag airfoil sections. Slide 14

Optimum Hovering Rotor Minimum P i requires λ=const. (uniform inflow) Minimum P 0 requires α= α(min C d /C l )= α 1 Then for minimum induce power θ= θ tip /r and each blade element must operate at α 1 With (BEMT) dc T =4λ 2 rdr then: Slide 15

Optimum Hovering Rotor We have seen that the minimum induced power requires a uniform inflow. Therefore the previous equation is constant over the disk. Let s assume that α 1 is the same for all airfoils along the blade and is independent of Re and M From the equation since α 1 =const and we now that λ=const then σr must be constant too. Slide 16

Optimum Hovering Rotor The previous situation is achieved when Slide 17

Optimum Hovering Rotor Let s now study the blade twist for uniform inflow and constant AOA: Expressing now in terms of the blade pitch Slide 18

Optimum Hovering Rotor The total thrust can be calculated The local solidity of the optimum rotor can be written: And the blade twist Slide 19

Optimum Hovering Rotor Slide 20

Optimum Hovering Rotor Slide 21

Power estimates In a real rotor the non uniformity of λ means that we must calculate numerically: We also saw that we could then calculate the induce power factor κ : Slide 22

Power estimates Let s now calculate the rotor profile power. Remember that (BET): We can use for C d the expression: Slide 23

Power estimates Slide 24

Power estimates The profile power is then: Slide 25

Power estimates For a ideally twisted blade λ=const and θr= θ tip : Since Slide 26

Power estimates Slide 27

Figure of Merit Since we now have the expressions for all power coefficients we can calculate the rotor FM: Slide 28

Prandtl s Tip-loss Function The idea of tip losses was already been introduced with the factor B Prandtl provided a solution to the problem of loss of lift near the blade tip. In this solution the curved helical vortex sheets of the rotor wake were replaced by a series of 2D sheets, meaning that the blade tip radius of curvature is large Slide 29

Prandtl s Tip-loss Function Prandtl s final result can be expressed in terms of an induced velocity correction factor F: Where f: Remember that Slide 30

Prandtl s Tip-loss Function Non dimensional radial position r Slide 31

Prandtl s Tip-loss Function Note that when N b (actuator disk) then F 1 The function can be incorporated in the BEMT: Since from BET We can write Slide 32

Prandtl s Tip-loss Function With the solution Since F is a function of λ the solution must be found numerically Slide 33

Compressibility Corrections So far we have considered that all aerodynamic characteristics independent of Mach number. To introduce a correction to take into account the influence of M, let s used Glauert s rule: Slide 34

Compressibility Corrections The local blade M is: Which gives a lift-curve-slope correction of: We had obtained: Slide 35

Compressibility Corrections Assuming the ideal twist (uniform inflow): Calculating the total thrust coefficient: Slide 36

Compressibility Corrections In the previous expression: If M 0 then K 1 and we obtain the incompressible result Slide 37

Weighted Solidity We have seen that for a optimum rotor c varies with r. In these cases where c varies along the blade span the rotor solidity is different from the local blade solidity: The objective of weighted solidity is to help compare performance of different rotors with different blade planforms. Slide 38

Trust Weighted Solidity The thrust coefficient is: Assuming a constant C l : With the equivalent chord Slide 39

Power Weighted Solidity Extending the study to the power coefficient Therefore Slide 40