Jc Propeller Design Program V. 3.11

PROPULSION

JC-10: Terms, Definitions, and Symbols (17) Apply JC-10: Terms, Definitions, and Symbols filter JC-11: Mechanical Standardization (600) Apply JC-11: Mechanical Standardization filter.

Sections :

Blade Element Theory for Propellers

The differentiation in the operator V is made with respect to the (jc, y, z). DESIGN: UniCopter Pusher. The drawing currently shows the ground clearance below the prop at 11.5'. The chord length, c, is the straight- line width of the propeller at a given distance along the radius. Objective was met by creating a computer program to design a propeller that meets user- defined aircraft performance requirements within the limitations of the motor, material, and manufacturing methods. CAESES provides comprehensive functionality for propeller and fan designers so that it can be used as an expert blade design software. Basically, any kind of propeller blade (e.g. Boat propeller, aircraft propeller, blowers, fans etc.) for any application can be created with it.

A relatively simple method of predicting the performance of a propeller (as well as fans or windmills) is the use of Blade Element Theory.

In this method the propeller is divided into a number of independent sections along the length. At each section a force balance is applied involving 2D section lift and drag with the thrust and torque produced by the section. At the same time a balance of axial and angular momentum is applied. This produces a set of non-linear equations that can be solved by iteration for each blade section. The resulting values of section thrust and torque can be summed to predict the overall performance of the propeller.

The theory does not include secondary effects such as 3-D flow velocities induced on the propeller by the shed tip vortex or radial components of flow induced by angular acceleration due to the rotation of the propeller. In comparison with real propeller results this theory will over-predict thrust and under-predict torque with a resulting increase in theoretical efficiency of 5% to 10% over measured performance. Some of the flow assumptions made also breakdown for extreme conditions when the flow on the blade becomes stalled or there is a significant proportion of the propeller blade in windmilling configuration while other parts are still thrust producing.

The theory has been found very useful for comparative studies such as optimising blade pitch setting for a given cruise speed or in determining the optimum blade solidity for a propeller. Given the above limitations it is still the best tool available for getting good first order predictions of thrust, torque and efficiency for propellers under a large range of operating conditions.

Blade Element Subdivision

A propeller blade can be subdivided as shown into a discrete number of sections.For each section the flow can be analysed independently if the assumption is made that for each there are only axial and angular velocity components and that the induced flow input from other sections is negligible. Thus at section AA (radius = r) shown above, the flow on the blade would consist of the following components.V₀ -- axial flow at propeller disk, V₂ -- Angular flow velocity vector

V₁ -- section local flow velocity vector, summation of vectors V₀ and V₂

Since the propeller blade will be set at a given geometric pitch angle (θ) the local velocity vector will create a flow angle of attack on the section. Lift and drag of the section can be calculated using standard 2-D aerofoil properties. (Note: propellers use a changed reference line : zero lift line not section chord line). The lift and drag components normal to and parallel to the propeller disk can be calculated so that the contribution to thrust and torque of the compete propeller from this single element can be found.

The difference in angle between thrust and lift directions is defined as

$$φ=θ-α$$

The elemental thrust and torque of this blade element can thus be written as

$$ΔT=ΔLcos(φ)-ΔDsin(φ)$$
$${ΔQ}/r=ΔDcos(φ)+ΔLsin(φ)$$

Substituting section data (C_L and C_D for the given $α$) leads to the following equations.

$$ΔL=C_L 1/2 ρV_1^2 c.dr$$
$$ΔD=C_D 1/2 ρV_1^2c.dr$$
per blade

where $ρ$ is the air density, c is the blade chord so that the lift producing area of the blade element is c.dr.

If the number of propeller blades is (B) then,

$$ΔT=1/2 ρV_1^2c(C_Lcos(φ)-C_Dsin(φ))B.drtext' (1)'$$
$$ΔQ=1/2 ρ V_1^2c(C_Dcos(φ)+C_Lsin(φ))B.drtext' (2)'$$

Inflow Factors

Jc Propeller Design Program V. 3.11 Download

A major complexity in applying this theory arises when trying to determine the magnitude of the two flow components V₀ and V₂. V₀ is roughly equal to the aircraft's forward velocity (V_∞) but is increased by the propeller's own induced axial flow into a slipstream. V₂ is roughly equal to the blade section's angular speed ($Ωr$) but is reduced slightly due to the swirling nature of the flow induced by the propeller. To calculate V₀ and V₂ accurately both axial and angular momentum balances must be applied to predict the induced flow effects on a given blade element. As shown in the following diagram, the induced components can be defined as factors increasing or decreasing the major flow components.

A typical streamtube of flow passing through section AA would have velocities

$$V_{θ}=bΩrtext' and 'V_x=V_{∞}(1+a)$$

So for the velocities V₀ and V₂ as shown in the previous section flow diagram,

$$V_0=V_{∞}+aV_{∞}$$

where a is the axial inflow factor and

$$V_2 = Ωr-bΩr$$

where b is the angular inflow factor (swirl factor)

The local flow velocity and the angle of attack for the blade section is thus

$$V_1=√{V_0^2+V_2^2}text' (3)'$$
$$α=θ-tan^{-1}(V_0/V_2)text' (4)'$$

Axial and Angular Flow Conservation of Momentum

The governing principle of conservation of flow momentum can be applied for both axial and circumferential directions. For the axial direction, the change in flow momentum along a stream-tube starting upstream, passing through the propeller at section AA and then moving off into the slipstream, must equal the thrust produced by this element of the blade.

To remove the unsteady effects due to the propeller's rotation, the stream-tube used is one covering the complete area of the propeller disk swept out by the blade element and all variables are assumed to be time averaged values.

[ table ΔT, =, text 'Change in Momemtum flow rate through tube at disk'; , =, text 'Mass flow rate in tube x Change in velocity along tube'; , =, ρ.2πr.dr.V_0×(V_s-V_{∞})]

By applying Bernoulli's equation and conservation of momentum, for the three separate components of the tube, from freestream to face of disk, from rear of disk to slipstream far downstream and balancing pressure and area versus thrust, it can be shown that the axial velocity at the disk will be the average of the freestream and slipstream velocities.

$$V_0={V_{∞}+V_s}/2$$

which means that

$$V_s=V_{∞}(1+2a)$$

Therefore

[table ΔT, =,2πrρV_{∞}(1+a)(V_{∞}(1+2a)-V_{∞}).dr; , =,4πrρV_{∞}^2(1+a)a.dr text' (5)']

For angular momentum

[table ΔQ, =,text'Change in angular momentum rate of flow in tube x radius'; , =,text'Mass flow rate in tube x change in circumferential velocity x radius'; , =,ρ.2πr.dr.V_0×(V_{θ(slipstream)}-V_{θ(freestream)})×r ]

By considering conservation of angular momentum in conjunction with the axial velocity change, it can be shown that the angular velocity in the slipstream will be twice the value at the propeller disk.

$$V_{θ(slipstream)}=2bΩrtext' and 'V_{θ(freestream)}=0$$

Thus

[table ΔQ, =, 2πrρV_{∞}(1+a)(2bΩr)r.dr; , =, 4πr^3V_{∞}(1+a)bΩ.drtext' (6)' ]

Because these final forms of the momentum equation balance still contain the variables for element thrust and torque, they cannot be used directly to solve for inflow factors. However there now exists a nonlinear system of equations (1),(2),(3),(4),(5) and (6) containing the four primary unknown variables ΔT, ΔQ, a, b, so an iterative solution to this system is possible.

Jc Propeller Design Program V. 3.11 Full

Iterative Solution procedure for Blade Element Theory.

Jc Propeller Design Program V. 3.11 Pdf

The method of solution for the blade element flow will be to start with some initial guess of inflow factors (a) and (b). Use these to find the flow angle on the blade (equations (3),(4)), then use blade section properties to estimate the element thrust and torque (equations (1),(2)). With these approximate values of thrust and torque equations (5) and (6) can be used to give improved estimates of the inflow factors (a) and (b).

This process can be repeated until values for (a) and (b) have converged to within a specified tolerance.

It should be noted that convergence for this nonlinear system of equations is not guaranteed. It is usually a simple matter of applying some convergence enhancing techniques (ie Crank-Nicholson under-relaxation) to get a result when linear aerofoil section properties are used. When non-linear properties are used, ie including stall effects, then obtaining convergence will be significantly more difficult.

For the final values of inflow factor (a) and (b) an accurate prediction of element thrust and torque will be obtained from equations (1) and (2).

Propeller Thrust and Torque Coefficients and Efficiency.

The overall propeller thrust and torque will be obtained by summing the results of all the radial blade element values.

$$T=ΣΔT text' (for all elements) and 'Q=ΣΔQ text' (for all elements)'$$

The non-dimensional thrust and torque coefficients can then be calculated along with the advance ratio at which they have been calculated.

$$C_T={T}/{ρn^2D^4}$$

and

$$C_Q={Q}/{ρn^2D^5}$$

where n is the rotation speed of propeller in revs per second and D is the propeller diameter.

The efficiency of the propeller under these flight conditions will then be

$$η_{prop}={JC_T}/{2πC_Q}$$

where advance ratio J is $J=V_{∞}/{nD}$.

Software Implementation of Blade Element Theory

Two programming versions of this propeller analysis technique are available. The first is a demonstration program which can be used to calculate thrust and torque coefficients and efficiency for a relatively simple propeller design using standard linearised aerofoil section data. The blade is assumed to have a constant pitch (p) so that the variation of θ with radius is calculated from the standard pitch equation.

$$p=2πrtan(θ)$$

Click here to download the Propeller Analysis Program Propel.exe (MS Windows Executable)

The second option is a MATLAB script file for the implementation of this method. The source code in this script is by default a simple propeller design with linear properties. However with the inclusion of your own propeller geometry and section data a more accurate analysis of the specific propeller design can be obtained.

Click here to download Propeller MATLAB script: Propel.m