Airfoil Tools
Search 1638 airfoils Google+
You have 0 airfoils loaded.
Your Reynold number range is 50,000 to 1,000,000. (set)

Lift & drag polars

Many of the airfoils have polar diagrams which can be viewed in the details and comparison section sections of the site. These show the change in lift coefficient (Cl), drag coefficient (Cd) and pitching moment (Cm) with angle of attack (alpha). There is also a graph of lift coefficient (Cl) against drag coefficient (Cd) which gives the theoretical glide angle of the airfoil.


All the polar diagrams currently available have been produced using Xfoil, an application created by Mark Drela and Harold Youngren for the design and analysis of subsonic airfoils. There is more information and details on the limitations of the analysis on the Xfoil page.

I wanted to create polar diagrams for all the airfoils in the database for a range of Reynolds numbers and ncrit values that would typically be applicable for model or light aircraft or small wind turbines. The polars are keyed by three values.

  • Reynolds number - These range from 50,000 to 1,000,000 in approximatly logarithmic steps. The Reynolds number is a dimensionless value that depends on the velocity, wing chord and fluid. There is more information on the Reynolds number calculator page.
  • Ncrit value - This is used to model of the turbulence of the fluid or roughness of the airfoil. The Xfoil documentation has a section explaining the method and also the table (reproduced below) of suggested values:-
    sailplane12 to 14
    motorglider11 to 13
    clean wind tunnel10 to 12
    average wind tunnel9
    dirty wind tunnel4 to 8
  • Mach number - This has been left at the default value of zero.

The polars were initially generated for an angle of attack (alpha) range -20° to +20° at 0.25° intervals. Some points never converge so this range may not be a continuous. The Xfoil lift and drag predictions are considered valid just beyond maximum lift (Cl max). Depending on the shape of the airfoil and Reynolds number this range may be some way past Cl max. After generating the polar files they were scanned to detect the maximum and minimum Cl values and truncated just past these values. Please be aware that most of this process was scripted so the cut off points may not be ideal. If too few points remain after this process then the polar was not included.

Xfoil has a command line user interface and this was automated using an expect script running under Cygwin on a Windows PC. The airfoil dat files are run through the parser to convert them to a Selig format and remove any extra text so they could be read by Xfoil. The commands to load the dat file and set up the session was fairly straight forward but it was more difficult to handle the user interaction when the iterative calculations failed to converge and restart Xfoil when it froze.

The commands for a typical scripted Xfoil session to generate a polar are below. The xfoil ALFA command was used for each individual angle of attack rather than the ASEQ as it was easier to automate, determine the succcess of the last calculation calculation and "step over" any failures. The polars were generated stating at ALFA = 0 in first a positive and then a negative direction. If a solution was found the next ALFA value was immediatley calculated. If no solution was found an INIT command was run before moving on to the next ALFA command.

LOAD dat/e1211-il.datLoad the dat file
MDESGo to the MDES menu
FILTSmooth any variations in the dat file data
EXECExecute the smoothing
Back to main menu
PANESet the number and location of the airfoil points for analysis
OPERGo to the OPER menu
ITER 70Max number of iterations set to 70 for convergence
RE 50000Set Reynolds number (required?)
VISC 50000Set viscous calculation with Reynolds number
PACCStart polar output file
polar/e1211-il_50000.txtThe output polar file name
No dump file
ALFA 0Calculate lift and drag at 0° angle of attack
ALFA 0.25... 0.25°
ALFA 0.5... 0.5° ...
......more alpha calculations here ...
ALFA 3.5At 3.5° no convergence
ALFA 3.5... try again ...
ALFA 3.5... and again
INIT Run INIT to reinitialise
ALFA 3.75Skip to 3.75° of alpha calculations here ...
PACCClose polar file
VISCReinitialise viscous calculation (required?)
Down to main menu
QUITExit Xfoil