Hi All,
since I noticed that there are not only many fans of WWII aircraft, but also many people actually interested in calculation of basic performance indicators, I would like to post some useful equations. I have found the exact approach of aerodynamics engineers quite bothersome to use and tried (successfully) to reduce the equations as much as possible, while retaining sufficient degree of accuracy.
The simplest approach to calculate time to perform sustained turn at 1000 meters height:
t = (W/b) * (Vmax/P)^0,5 * (1/8,80)
where W is weight in kg, b is span in meters, V is max sea level speed in km/h and P is power at the maximum sea level speed (measured in metric horspower, i.e. 735.5 W). The constant is empirically derived, but the rest of the equation follows from following reasoning.
Time to perform sustained turn is inversely proportional to turn rate, which is defined as g*[(n^2)-1)]^0.5/V (g is 9.81 m/s2, n is load factor in g's and speed in m/s in sustained turn).
The speed in sustained turn is
V = [(P*eff*735.5)/(2*ro*S*Cd0)]^(1/3) m/s
and the load in sustained turn is
n = 0.687 * {[(P*eff*735.5)^2*E*ro*S*K]/W}^(1/3) m/s2
P is metric horsepower, eff is propeller efficiency coefficient, ro is air density in kg/m3, S is wing area in m2, Cd0 is zero lift drag coefficient, E is aircraft maximum lift-to-drag ratio, K is induced drag coefficient and W is weight in kg.
K = pi*AR*oswald (oswald means Oswald span efficiency)
AR = b^2/S
E = [(K/Cd0)^0.5]/2
Assuming that propeller efficiency and Oswald is constant for all aircraft (not exactly true) and we are moving at sea level, we can rewrite the above equations without numerical constants (i.e. rewrite just variables, ~ denotes "is in proportion to") as:
V ~ P^1/3 * S^(-1/3) * Cd0^(-1/3)
n ~ P^2/3 * S^(1/3) * Cd0^(-1/6) * K^0.5 * W^(-1)
Cd0 can be estimated numerically, using formula Cd0 = Cd - Cl^2/K, but for the purpose of further simplification it can be assumed as equal (in max speed level flight) to total coefficient of drag. Error from this assumption is negligible (on the order of less than 1% accuracy in final time to turn computation). This is because of the 90 percent share of Cd0 in total drag (as long as you take max speed level flight conditions). Total coefficient of drag is proportional to P * Vmax^(-3) * S^(-1), where the Vmax denotes sea level max speed. Inserting this into equations above (and assuming, that the turn is done at maximum power) we get:
V ~ Vmax
This is not that much surprising. Speed for best load sustained turn (at SL) occurs at constant percentage of max speed for that altitude, usually 64 percent. Variations in actual aircraft seem to be under 1%, when the result was computed exactly with all coefficients. The second equation (after further substitution for K ~ (b^2)/S) yields:
n ~ P^0.5 * Vmax^0.5 * b * W^(-1)
Now because of the subtraction in [(n^2)-1]^0.5, there is not much to simplify here. But since the aircraft I am mostly interested in (WWII pistons, usually fighters) are all making the sustained turns at 2-4 g's, I can simply ignore the subtraction. The ordering of aircraft (i.e. who turns better) should not be affected and the numerical difference between 2.5 g and 3.5 g aircraft is 140 percent (if the subtraction is ignored) and 146,4 percent, when it is not ignored. Sizing the final constant in the middle of aircraft of interest, it is thus possible to have less than 3 percent error in estimates of turn time.
The final equation is therefore:
t ~ P^0.5 * Vmax^(-0.5) * b * W^(-1)
The constant for sustained turn time in 1000 meters height and metric units is 8.80 (see the first equation). For sea level turn it comes as 9.48.
For those using imperial units, 1000 meters constant is 4.69 and sea level 5.05 (P in horses, Vmax in statute miles per hour, b in feet and W in pounds).
Example of use:
Fw-190 A-3, 1560 HP @ SL, 540 km/h @ SL, span 10.5m, 3850 kg.
t = (3850/10.5) * (540/1560)^0.5 * (1/8.80) = 24.5 seconds
Possible errors: Oswald span efficiency and propeller efficiency varies from aircraft to aircraft. Though it is possible to compute them, they vary so slightly, that there should be no significant errors. Above example yields the same result with exact formulas for 0.78 propeller efficiency and 0.868 oswald. If the propeller efficiency was 0.858 (+10 percent), turn time would be 23.2 seconds and with Oswald 0.7812 (-10 percent) turn time rises to 26.1 seconds. Thus +- 10 percent differences produce about 5-6 percent difference in results. With variable pitch propellers of second world war I don't expect much that much effective difference, while the differences in Oswald span efficiency are very low between these high aspect ratio aircraft (for AR between 5 and 7, Raymer computes 0.9 to 0.84 Oswald efficiency). Thus I don't expect these factors to be significant.
Other than this, the power in 1000 meters is different than at sea level. Using sea level power is however convenient, and it introduces usually only small error. E.g. Yak-1 has about 1020 hp at sea level and 1050 at 2000 meters. The error in square root of 1035 versus 1020 is 0.7 percent. In weight terms it is for 2900 kg fighter equal to 20 kg (less than average ammunition load).
The last known problem is data. If you have no sea level speed and power data, you can for first orientation use data at known altitude and correct them for the different air density. Quick and dirty approach would be to replace the speed in above equation with (Valt*exp(-alt/28.0)) and use appropriate horsepower (i.e. power at that altitude); Valt is speed in altitude (km/h), alt is the altitude (m). But remember, garbage in, garbage out. Try not to abuse the god of physics and take care to avoid this approach as much as possible. Some engines are especially unsuitable for this approach, which essentially supposes, that the engine power remains the same with altitude.
Example:
Spitfire Mk VB, 1470 HP @ 6100m, 597 km/h @ 6100m, span 11.23m, 2950 kg.
t = [(exp(-6.1/28.0) * 597)/1470]^0.5 * (2950/11.23) * (1/8.80) = 17.1 seconds
I apologize for such a long post and hope you find the equations useful. If you have any feedback, find any errors or have questions, don't hesitate to ask.
Zdenek,
badger 45
P.S. note that the span loading and speed-to-power coefficients are after all the only things that affect the aircrafts turning rate. Thus it seems that the turn time was really more or less unsubstantial characteristic in WWII design, since every aircraft designer in a quest for performance sought to maximize these and thus inevitably worsened turning characteristics of aircrafts.
since I noticed that there are not only many fans of WWII aircraft, but also many people actually interested in calculation of basic performance indicators, I would like to post some useful equations. I have found the exact approach of aerodynamics engineers quite bothersome to use and tried (successfully) to reduce the equations as much as possible, while retaining sufficient degree of accuracy.
The simplest approach to calculate time to perform sustained turn at 1000 meters height:
t = (W/b) * (Vmax/P)^0,5 * (1/8,80)
where W is weight in kg, b is span in meters, V is max sea level speed in km/h and P is power at the maximum sea level speed (measured in metric horspower, i.e. 735.5 W). The constant is empirically derived, but the rest of the equation follows from following reasoning.
Time to perform sustained turn is inversely proportional to turn rate, which is defined as g*[(n^2)-1)]^0.5/V (g is 9.81 m/s2, n is load factor in g's and speed in m/s in sustained turn).
The speed in sustained turn is
V = [(P*eff*735.5)/(2*ro*S*Cd0)]^(1/3) m/s
and the load in sustained turn is
n = 0.687 * {[(P*eff*735.5)^2*E*ro*S*K]/W}^(1/3) m/s2
P is metric horsepower, eff is propeller efficiency coefficient, ro is air density in kg/m3, S is wing area in m2, Cd0 is zero lift drag coefficient, E is aircraft maximum lift-to-drag ratio, K is induced drag coefficient and W is weight in kg.
K = pi*AR*oswald (oswald means Oswald span efficiency)
AR = b^2/S
E = [(K/Cd0)^0.5]/2
Assuming that propeller efficiency and Oswald is constant for all aircraft (not exactly true) and we are moving at sea level, we can rewrite the above equations without numerical constants (i.e. rewrite just variables, ~ denotes "is in proportion to") as:
V ~ P^1/3 * S^(-1/3) * Cd0^(-1/3)
n ~ P^2/3 * S^(1/3) * Cd0^(-1/6) * K^0.5 * W^(-1)
Cd0 can be estimated numerically, using formula Cd0 = Cd - Cl^2/K, but for the purpose of further simplification it can be assumed as equal (in max speed level flight) to total coefficient of drag. Error from this assumption is negligible (on the order of less than 1% accuracy in final time to turn computation). This is because of the 90 percent share of Cd0 in total drag (as long as you take max speed level flight conditions). Total coefficient of drag is proportional to P * Vmax^(-3) * S^(-1), where the Vmax denotes sea level max speed. Inserting this into equations above (and assuming, that the turn is done at maximum power) we get:
V ~ Vmax
This is not that much surprising. Speed for best load sustained turn (at SL) occurs at constant percentage of max speed for that altitude, usually 64 percent. Variations in actual aircraft seem to be under 1%, when the result was computed exactly with all coefficients. The second equation (after further substitution for K ~ (b^2)/S) yields:
n ~ P^0.5 * Vmax^0.5 * b * W^(-1)
Now because of the subtraction in [(n^2)-1]^0.5, there is not much to simplify here. But since the aircraft I am mostly interested in (WWII pistons, usually fighters) are all making the sustained turns at 2-4 g's, I can simply ignore the subtraction. The ordering of aircraft (i.e. who turns better) should not be affected and the numerical difference between 2.5 g and 3.5 g aircraft is 140 percent (if the subtraction is ignored) and 146,4 percent, when it is not ignored. Sizing the final constant in the middle of aircraft of interest, it is thus possible to have less than 3 percent error in estimates of turn time.
The final equation is therefore:
t ~ P^0.5 * Vmax^(-0.5) * b * W^(-1)
The constant for sustained turn time in 1000 meters height and metric units is 8.80 (see the first equation). For sea level turn it comes as 9.48.
For those using imperial units, 1000 meters constant is 4.69 and sea level 5.05 (P in horses, Vmax in statute miles per hour, b in feet and W in pounds).
Example of use:
Fw-190 A-3, 1560 HP @ SL, 540 km/h @ SL, span 10.5m, 3850 kg.
t = (3850/10.5) * (540/1560)^0.5 * (1/8.80) = 24.5 seconds
Possible errors: Oswald span efficiency and propeller efficiency varies from aircraft to aircraft. Though it is possible to compute them, they vary so slightly, that there should be no significant errors. Above example yields the same result with exact formulas for 0.78 propeller efficiency and 0.868 oswald. If the propeller efficiency was 0.858 (+10 percent), turn time would be 23.2 seconds and with Oswald 0.7812 (-10 percent) turn time rises to 26.1 seconds. Thus +- 10 percent differences produce about 5-6 percent difference in results. With variable pitch propellers of second world war I don't expect much that much effective difference, while the differences in Oswald span efficiency are very low between these high aspect ratio aircraft (for AR between 5 and 7, Raymer computes 0.9 to 0.84 Oswald efficiency). Thus I don't expect these factors to be significant.
Other than this, the power in 1000 meters is different than at sea level. Using sea level power is however convenient, and it introduces usually only small error. E.g. Yak-1 has about 1020 hp at sea level and 1050 at 2000 meters. The error in square root of 1035 versus 1020 is 0.7 percent. In weight terms it is for 2900 kg fighter equal to 20 kg (less than average ammunition load).
The last known problem is data. If you have no sea level speed and power data, you can for first orientation use data at known altitude and correct them for the different air density. Quick and dirty approach would be to replace the speed in above equation with (Valt*exp(-alt/28.0)) and use appropriate horsepower (i.e. power at that altitude); Valt is speed in altitude (km/h), alt is the altitude (m). But remember, garbage in, garbage out. Try not to abuse the god of physics and take care to avoid this approach as much as possible. Some engines are especially unsuitable for this approach, which essentially supposes, that the engine power remains the same with altitude.
Example:
Spitfire Mk VB, 1470 HP @ 6100m, 597 km/h @ 6100m, span 11.23m, 2950 kg.
t = [(exp(-6.1/28.0) * 597)/1470]^0.5 * (2950/11.23) * (1/8.80) = 17.1 seconds
I apologize for such a long post and hope you find the equations useful. If you have any feedback, find any errors or have questions, don't hesitate to ask.
Zdenek,
badger 45
P.S. note that the span loading and speed-to-power coefficients are after all the only things that affect the aircrafts turning rate. Thus it seems that the turn time was really more or less unsubstantial characteristic in WWII design, since every aircraft designer in a quest for performance sought to maximize these and thus inevitably worsened turning characteristics of aircrafts.