Tools used to create fine meshes used by the Prop Calculator based on propeller test data.
One of the steps in aircraft design process is calculating the performance figures. Max power, max airspeed and max service ceiling are the most prominent examples. In order to calculate them we need to know the engine's Coefficient of Power Cp and the propeller's Coefficient of Efficiency η. Both of these coefficients have their own charts, but we are going to focus on the efficiency. An example η chart for a 2 bladed prop can be seen below:
Where:
- V - Aircraft speed
- n - Propeller speed
- D - Propeller diameter
And - Advance Ratio - a dimensionless velocity describing the flow around the propeller blade
Designing an aircraft powertrain can be a tedious process when done manually. The designer is required to interpolate the data from the chart based on initial conditions in order to calculate the η. When the user chooses one of the angles from the chart the process is straightforward. Otherwise the interpolation becomes very time consuming and the precision of the results may be questioned. The whole process could easily be automated but the data needed to be processed first. That led to the birth of this toolkit.
The application is based on two reports published in 1938 and 1939 respectively by the NASA predecessor - National Advisory Committee for Aeronautics:
- NACA 640 - Tests of propellers having 2, 3 and 4 blades of different airfoils at blade angles up to 45°.
- NACA 658 - Tests of two propellers at blade angles up to 60°.
The higher blade angles are intended for high-speed aircraft. The latter of the reports has data only for propellers with 3 or 4 blades. This is a problem as it limits the configurations the designer can use.
The first step in exploring the data further is making it 3D. A new axis was assigned to the Angle parameter. The result can be seen below with the other blade numbers thrown in for good measure.
Every propeller type has different η characteristics. As you can see the 4 bladed prop is the only one having data up to 60° of Angle. The remaining data was approximated using polynomial regression.
This process can be described in a few simple steps:
- Calculate η ratios for a given pair of propellers (n- and 4-bladed).
- Normalize them and fit a polynomial.
- Use this polynomial to calculate unknown ratios for angles greater than 45°.
- Calculate approximated curves from the 4 bladed prop data and newly obtained ratios.
Normalized η ratios for a 2 bladed propeller. A 9th degree polynomial was chosen as it had the best fit to the underlying data.
As you can see the generated data fits very well with the test curves. The expanded data was then used as a basis for the mesh densing algorithm.
The mesh was created using curve intersections located at the same relative z for every data series. Splines were constructed on these nodes and then used to interpolate the green points. The data series on the back side of the surface have been hidden for better readability.
The mesh step size on the y axis was 0.5°, which was a great compromise between accuracy and performance.
The points were finally recalculated on a regular grid, which made drawing the mesh as a surface possible. This greatly improved accuracy in production usage and enabled clear results visualization.
You can see them in action on the Results tab of Prop Calculator.
- Python
- NumPy - This is the fundamental package for scientific computing with Python
- pandas - This is a fast and easy to use data analysis and manipulation tool
- SciPy - This is a Python library used for scientific and technical computing
- Plotly - This Python graphing library makes interactive, publication-quality graphs