Conceptual Design and Performance Optimization of a Tip Device for a Regional Turboprop Aircraft
Abstract
:1. Introduction
1.1. Induced Drag Reduction Techniques
1.2. Non-Planar Devices
1.3. Winglet Studies on Various Aircraft Platforms
- Decrease in fuel burn (or increase in range).
- Increase in flight ceiling.
- Reduced take-off runs.
- Increased time between engine maintenance.
1.4. General Framework for the Use of Winglets in Regional Turboprop Aircraft
2. Methodology
2.1. Geometry Modeling
2.1.1. Reference Model
2.1.2. Wing
2.1.3. Fuselage
2.1.4. Vertical and Horizontal Tailplanes
2.2. Aerodynamic Analysis
2.3. Validation of Reference Aerodynamic Model
2.4. Structural Sizing
2.4.1. Reference Model Structural Validation
2.4.2. Project Model Structural Design
2.5. Aircraft Performance Assessment
2.6. Python Integration Framework
- OpenVSP and VSPAero Runner: Create a pool of design points for the desired range of variables. Then, this is fed to OpenVSP for each design point in an AngelScript++ input file for OpenVSP containing case-specific design parameters. Then, the geometry generated by OpenVSP was read in VSPAero and executed to calculate the aerodynamic data for 1 and 2.5-g conditions.
- EMWET input file parsing: As EMWET requires the geometry, spanload, and the quarter chord pitching moment, the VSPAero output files were parsed and written onto the EMWET initialization and load files. Each individual EMWET case was appended onto a Matlab × m script that could then be run from Matlab. To process the results, two functions were developed in Python to ultimately output a comma-separated values file in ASCII for data visualization and allow oprimization work in the future.
- Parsing of VSPAero aerodynamic coefficients and associated wing weight and joining each case input to its output.
- Writing case input and output in a *.csv file.
3. Results
3.1. Cant Angle Effects
3.2. Span Effects
3.3. Pareto Front
3.4. Optimized Winglet Design
4. Discussion
5. Conclusions
Author Contributions
Funding
Acknowledgments
Conflicts of Interest
Acronyms
AR | Aspect Ratio |
BAe | British Aerospace |
CFD | Computational Fluid Dynamics |
CPFH | Cost Per Flying Hour |
DOC | Direct Operating Cost |
EMWET | Elham Modified Wing Weight Estimation Technique |
ISA | International Standard Atmosphere |
MAC | Mean Aerodynamic Chord |
MTOW | Maximum Take-Off Weight |
MZFW | Maximum Zero-Fuel Weight |
NACA | National Advisory Committee for Aeronautics |
NASA | National Aeronautics and Space Administration |
NFLC | National Flying Laboratory Centre |
OEW | Operating Empty Weight |
RPM | Revolutions Per Minute |
USAF | United States Air Force |
VLM | Vortex Lattice Method |
WRBM | Wing Root Bending Moment |
References
- Morris, J.; Rowbotham, A.; Morrell, P.; Foster, A.; Poll, I.; Owen, B.; Raper, D.; Mann, M.; Ralph, M. UK Aviation: Carbon Reduction Futures. Final Report to the Department for Transport; Cranfield University: Cranfield, UK, 2009. [Google Scholar]
- Farries, P.; Eyers, C. Aviation CO2 Emissions Abatement Potential from Technology Innovation; Office of Climate Change: London, UK, 2008.
- Prandtl, L. Uber TragflUgel des Kleinsten Induzierten Widerstandes. Z. Flugtech. Mot. 1933, 24, 305–306. [Google Scholar]
- Prandtl, L. Application of Modern Hydrodynamics to Aeronautics; NACA Technical Report No. 116; National Advisory Committee for Aeronautics: Washington, DC, USA, 1923.
- Jones, R.T. The Spanwise Distribution of Lift for Minimum Induced Drag of Wings Having a Given Lift and a Given Bending Moment; NACA Technical Note No. 2249; National Advisory Committee for Aeronautics: Washington, DC, USA, 1950.
- Hoerner, S. Aerodynamic Shape of the Wing Tips; United States Air Force Technical Report No. 5752; Wright-Patterson Air Force Base: Dayton, OH, USA, 1952. [Google Scholar]
- Jarrett, P.F.W. Lanchester and the Great Divide. J. Aeronaut. Hist. 2014. Available online: https://www.aerosociety.com/publications/jah-fw-lanchester-and-the-great-divide/ (accessed on 25 September 2019).
- Whitcomb, R.T. A Design Approach and Selected Wind Tunnel Results at High Subsonic Speeds for Wing-Tip Mounted Winglets; Technical Note TN-D-8260; NASA Langley Research Center: Hampton, VA, USA, 1976.
- Heyson, H.H.; Riebe, G.D.; Fulton, C.L. Theoretical Parametric Study of the Relative Advantages of Winglets and Wing-Tip Extensions; NASA Technical Memorandum TM X-74003; NASA Langley Research Center: Hampton, VA, USA, 1977.
- Jones, R.T.; Lasinski, T.A. Effects of Winglets on the Induced Drag of Ideal Wing Shapes; NASA Technical Memorandum TM-81230; NASA Ames Research Center: Moffett Field, CA, USA, 1980.
- Asai, K. Theoretical considerations in the aerodynamic effectiveness of winglets. J. Aircr. 1985, 22, 635–637. [Google Scholar] [CrossRef]
- Kroo, I.; McMasters, J.; Smith, S.C. Highly nonplanar lifting systems. In Proceedings of the Transportation Beyond 2000: Technologies Needed for Engineering Design, NASA Langley Research Center, Hampton, VA, USA, 26–28 September 1995. [Google Scholar]
- Ning, S.A.; Kroo, I. Tip extensions, winglets and C-wings: Conceptual design and optimization. In Proceedings of the 26th AIAA Applied Aerodynamics Conference, Honolulu, HI, USA, 18–21 August 2008. [Google Scholar]
- Elham, A.; van Tooren, M.J.L. Winglet multi-objective shape optimization. Aerosp. Sci. Technol. 2014, 37, 93–109. [Google Scholar] [CrossRef]
- Khosravi, S.; Zingg, D.W. A numerical optimization study on winglets. In Proceedings of the 15th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference, Atlanta, GA, USA, 16–20 June 2014. [Google Scholar]
- Jobe, C.E. Prediction of Aerodynamic Drag; Technical Memorandum AFWAL-TM-84-203; Flight Dynamics Laboratory: Wright-Patterson Air Force Base, OH, USA, 1984.
- National Research Council. Assessment of Wingtip Modifications to Increase the Fuel Efficiency of Air Force Aircraft; The National Academies Press: Washington, DC, USA, 2007. [Google Scholar]
- Lehmkuehler, K.; Wong, K.C. Winglet design for a Fairchild Merlin III using CFD analysis. In Proceedings of the AIAC14 14th Australian International Aerospace Congress, Melbourne, Australia, 28 February–3 March 2011. [Google Scholar]
- Nicolosi, F.; Della Vecchia, P.; Corcione, S. Design and aerodynamic analysis of a twin-engine commuter aircraft. Aerosp. Sci. Technol. 2015, 40, 1–16. [Google Scholar] [CrossRef]
- Della Vecchia, P. Aerodynamic guidelines in the design and optimization of new regional turboprop aircraft. Aerosp. Sci. Technol. 2014, 38, 88–104. [Google Scholar] [CrossRef]
- Conley, N. Winglet toe-out angle optimization for the Gates Learjet longhorn wing. J. Aircr. 1980, 17, 851–855. [Google Scholar] [CrossRef]
- Dees, P.; Stowell, M. 737–800 Winglet Integration; SAE: Warrendale, PA, USA, 2001. [Google Scholar]
- Streit, T.; Himisch, J.; Heinrich, R.; Nagel, B.; Horstmann, K.; Liersch, C. Design of a retrofit winglet for a transport aircraft with assessment of cruise and ultimate structural loads. In New Results in Numerical and Experimental Fluid Mechanics VI. Notes on Numerical Fluid Mechanics and Multidisciplinary Design (NNFM); Tropea, C., Jakirlic, S., Heinemann, H.J., Henke, R., Honlinger, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; Volume 96, pp. 62–70. [Google Scholar]
- Maughmer, M.D. The design of winglets for low-speed aircraft. Tech. Soar. 2006, 30, 6173. [Google Scholar]
- British Aerospace. Jetstream 31 Maintenance Manual; British Aerospace: Farnborough, UK, 1982. [Google Scholar]
- Cooke, A.K. A Simulation Model of the NFLC Jetstream 31; College of Aeronautics Report No. 0402; Cranfield University: Bedford, UK, 2006; Available online: https://dspace.lib.cranfield.ac.uk/handle/1826/7623 (accessed on 25 September 2019).
- Maughmer, M.D.; Swan, T.S.; Willits, S.M. Design and testing of a winglet airfoil for low-speed aircraft. J. Aircr. 2002, 39, 654–661. [Google Scholar] [CrossRef]
- Torenbeek, E. Advanced Aircraft Design, Conceptual Design, Analysis and Optimisation of Subsonic Civil Airplanes; John Wiley and Sons Ltd.: West Sussex, UK, 2013; p. 410. [Google Scholar]
- European Aviation Safety Agency (EASA). Type Certificate Data Sheet Jetstream 3100/3200 Series; EASA: Cologne, Germany, 2015. Available online: https://www.easa.europa.eu/documents/type-certificates/aircraft-cs-25-cs-22-cs-23-cs-vla-cs-lsa/easaa191#main-content (accessed on 25 September 2019).
- Conway, J.T. Exact actuator disc solutions for non-uniform heavy loading and slipstream contraction. J. Fluid Mech. 1998, 365, 235–267. [Google Scholar] [CrossRef]
- Stoll, A.M. Comparison of CFD and experimental results of the LEAPTech distributed electric propulsion blown wing. In Proceedings of the 15th AIAA Aviation Techonlogy, Integration and Operations Conference, AIAA Aviation forum (AIAA 2015-3188), Dallas, TX, USA, 22–26 June 2015. [Google Scholar]
- European Aviation Safety Agency (EASA). Certification Specification-25.331–341 (Large Aeroplanes); EASA: Cologne, Germany, 2007. Available online: https://www.easa.europa.eu/sites/default/files/dfu/CS-25_Amdt%203_19.09.07_Consolidated%20version.pdf (accessed on 25 September 2019).
- Lawson, N.J.; Jacques, H.; Gautrey, J.E.; Cooke, A.K.; Holt, J.C.; Garry, K.P. Jetstream 31 national flying laboratory: Lift and drag measurement and modelling. Aerosp. Sci. Technol. 2017, 60, 84–95. [Google Scholar] [CrossRef] [Green Version]
- Dababneh, O.; Kipouros, T. A review of aircraft wing mass estimation methods. Aerosp. Sci. Technol. 2018, 72, 256–266. [Google Scholar] [CrossRef]
- Elham, A.; La Rocca, G.; van Tooren, M.J.L. Development and implementation of an advanced, design-sensitive method for wing weight estimation. Aerosp. Sci. Technol. 2013, 29, 100–113. [Google Scholar] [CrossRef]
- Mariens, J.; Elham, A.; van Tooren, M.J.L. Influence of weight modelling on the outcome of wing design using multidisciplinary design optimisation techniques. Aeronaut. J. 2013, 117, 871–895. [Google Scholar] [CrossRef]
- Honeywell Official Web Site. Available online: https://aerospace.honeywell.com/en/~/media/aerospace/files/brochures/n61-1491-000-000-tpe331-10turbopropengine-bro.pdf (accessed on 25 September 2019).
- Elham, A. Weight Indexing for Multidisciplinary Design Optimization of Lifting Surfaces. Ph.D. Thesis, TU Delft, Delft, The Netherlands, 2013. [Google Scholar]
- Kundu, A.K. Aircraft Design; Cambridge University Press: Cambridge, MA, USA, 2012; ISBN 9780511844652. [Google Scholar]
- Roskam, J. Airplane Design: Component Weight Estimation (Part V); DAR Corporation: Lawrence, UK, 1999; ISBN 9781884885501. [Google Scholar]
- Raymer, D.P. Aircraft Design: A conceptual Approach, 5th ed.; AIAA Education Series; American Institute of Aeronautics and Astronautics, Inc.: Reston, VA, USA, 2012; p. 495. ISBN 9781600869112. [Google Scholar]
- Forecast International Official Web Site. Jetstream 31 Outlook. Available online: https://www.forecastinternational.com/archive/disp_old_pdf.cfm?ARC_ID=305 (accessed on 25 September 2019).
- Anderson, J. A Survey of Multiobjective Optimization in Engineering Design; Technical Report LiTH-IKP-R-1097; Department of Mechanical Engineering, Linkoping University: Linkoping, Sweden, 2000. [Google Scholar]
- Gautham, N.; Bibin, J. Effect of winglets induced tip vortex structure on the performance of subsonic wings. Aerosp. Sci. Technol. 2016, 58, 328–340. [Google Scholar]
- Gueraiche, D.; Popov, S. Winglet geometry impact on DLF-F4 aerodynamics and an analysis of a hyperbolic winglet concept. Aerospace 2017, 4, 60. [Google Scholar] [CrossRef]
- Guerrero, J.; Sanguineti, M.; Wittkowski, K. CFD study of the impact of variable cant angle winglets on total drag reduction. Aerospace 2018, 5, 126. [Google Scholar] [CrossRef]
- Heller, G.; Kreuzer, P.; Dirmeier, S. Development and integration of a new high performance wingtip device for transonic aircraft. In Proceedings of the ICAS 2002 Congress, Toronto, ON, Canada, 8–13 September 2002. [Google Scholar]
- Liebeck, R.H.; Andrastek, D.A.; Chau, J.; Girvin, R.; Lyon, R.; Rawdon, B.K.; Scott, P.W.; Wright, R.A. Advanced Subsonic Airplane Design and Economic Studies; NASA TECHNICAL Report CR-195443; NASA Lewis Research Center: Cleveland, OH, USA, 1995.
- Lappas, I.; Bozoudis, M. The development of an Ordinary Least Squares parametric model to estimate the Cost Per Flying Hour of ‘unknown’ aircraft types and a comparative application. Aerospace 2018, 5, 104. [Google Scholar] [CrossRef]
Type | Year of First Flight | Maximum Operating Altitude (ft) | Range (Km) | MTOW (Kg) | Aspect Ratio | Winglets Use |
---|---|---|---|---|---|---|
DHC-6 Twin Otter (Series 400) | 2010 | 25,000 | 1480 | 5670 | 10.0 | No |
Harbin Y-12 | 1982 | 23,000 | 1340 | 5300 | 8.7 | No |
Beechcraft 1900D | 1982 | 25,000 | 707 | 7764 | 10.8 | Yes |
Dornier Do-228 | 1981 | 25,000 | 396 | 6575 | 9.0 | No |
BAe Jetstream 31 | 1980 | 25,000 | 1260 | 6950 | 10.0 | No |
Let L-410 Turbolet | 1969 | 27,500 | 510 | 6600 | 11.5 | No |
Fairchild Swearingen Metroliner | 1969 | 25,000 | 2131 | 6577 | 10.5 | No |
Embraer EMB 110 Bandeirante | 1968 | 21,500 | 1964 | 5900 | 8.1 | No |
Variable (Unit) | Values | Justification/Constraints |
---|---|---|
Winglet position | 0.8, 0.98 | As a ratio of wing-tip chord. It was to either start from the leading edge (0.98) or recessed rearwards to minimize the wetted area, but still in front of the front spar position to facilitate integration to the front spar in subsequent detailed design. |
Winglet taper ratio | 0.25, 0.5 | Combination of Lehmkuehler’s [18] and Della Vecchia [20] optimal taper ratios. |
Cant angle (absolute, from XY plane) (deg) | 7, 43.5, 80 | Planar 7 degrees, same as the main wing dihedral. An 80-degree maximum angle, with 43.5 being the middle value of the range. |
Winglet span (or winglet height, as a ratio of wing semispan) | 10%, 15%, 20% | As a ratio of the wing semispan with the upper boundary constrained by the Jetstream 31 certification airport reference code (B-II), stipulating a maximum 24 m of span. The upper boundary was scaled down to a more reasonable maximum 20% of semispan. |
Quarter chord sweep (deg) | 14, 26, 38, 50 | Although mostly suitable as a parameter for aircraft operating in the transonic regime, it is kept to investigate the aerodynamics and the torsional effects on the wing structure. |
Toe-out (deg) | −4, −2.5, −1 | Toe-out angle dictated by the local lift coefficient requirement. |
Twist (deg) | −4, −2.5, −1 | Twist controls the winglet spanload [20]. |
Parameter (Unit) | Value |
---|---|
Projected span (m) | 15.85 |
Area (m2) | 25.902 |
Centerline chord (m) | 2.375 |
Root chord (m) | 2.16 |
Tip chord (m) | 0.79 |
Sweep at 30% of chord (deg) | 0 |
Aspect ratio | 10.0 |
Reference chord: Mean aerodynamic chord (MAC) (m) | 1.87 |
Root setting angle (deg, applied at the 25% of the chord) | 3 |
Parameter (Unit) | Value |
---|---|
Circular maximum diameter (m) | 1.981 |
Length (m) | 13.347 |
Parameter (Unit) | Value |
---|---|
Projected span (m) | 6.60 |
Area (m2) | 7.80 |
Centerline chord (m) | 1.676 |
Tip chord (m) | 0.6855 |
Sweep at 25% of chord (deg) | 7.10 |
Aspect Ratio | 5.60 |
Mean aerodynamic chord (MAC) (m) | 1.181 |
Root aerofoil | NACA 0012 |
Tip aerofoil | NACA 0010 |
Dihedral (deg) | 0 |
Parameter (Unit) | Value |
---|---|
Projected span (m) | 3.32 |
Area (m2) | 6.65 |
Centerline chord (m) | 3.20 |
Tip chord (m) | 0.88 |
Sweep at 25% of chord (deg) | 7.10 |
Aspect Ratio | 5.60 |
Mean aerodynamic chord (MAC) (m) | 2.04 |
Root aerofoil | NACA 0012 |
Tip aerofoil | NACA 0010 |
Parameter (Symbol) (Unit) | Value |
---|---|
Altitude (ft) | 25,000 |
Ground speed (V∞) (km/h) | 425 |
Mach number (M) | 0.38 |
Density (ρ) (International Standard Atmosphere (ISA), kg/m3) | 0.5495 |
Reynolds number (Re) | 7.98 × 106 (MAC) |
Reference chord (m) | 1.87 (MAC) |
Reference area (Sref) (m) | 24.952 |
Reference span (m) | 15.85 |
Parameter | Lawson et al. [32] | VSPAero Prediction for the Developed Model | % Difference |
---|---|---|---|
Lift curve slope | 0.1052 | 0.1101 | 4.66 |
Lift coefficient for zero angle of attack | 0.3305 | 0.3323 | −2.48 |
Lift coefficient at 3-degree angle of attack | 0.6461 | 0.655 | −1.4 |
Zero lift drag coefficient | 0.0422 | 0.0433 | 0.26 |
Drag coefficient at 3-degree angle of attack | 0.0629 | 0.0641 | 1.78 |
L/D at zero angle of attack | 7.472 | 6.985 | −6.52 |
L/D at 3-degree angle of attack | 10.27 | 10.23 | −0.4 |
Maximum L/D | 10.5 | 10.78 | 2.66 |
Fuel Weight Fraction (Mffi) | Turboprop Aircraft |
---|---|
Start and warm up | 0.990 |
Taxi | 0.995 |
Take off | 0.995 |
Climb | 0.985 |
Cruise | Calculated |
Descent | 0.985 |
Landing, taxi, and shutdown | 0.995 |
Weights | Value (Kg) |
---|---|
Reference wing | 449.52 |
MTOW | 7350 |
MZFW | 6350 |
Maximum payload | 1935 |
Maximum fuel | 1491 |
OEW | 4415 |
Parameters | Value |
---|---|
Winglet position | 0.8 |
Cant | 80 (deg.) |
Span | 0.1 |
Sweep | 26 (deg.) |
Toe-out | −1 (deg.) |
Twist | −1 (deg.) |
CD/CDref | −1.186% |
W/Wref | 3.25% |
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Lappas, I.; Ikenaga, A. Conceptual Design and Performance Optimization of a Tip Device for a Regional Turboprop Aircraft. Aerospace 2019, 6, 107. https://doi.org/10.3390/aerospace6100107
Lappas I, Ikenaga A. Conceptual Design and Performance Optimization of a Tip Device for a Regional Turboprop Aircraft. Aerospace. 2019; 6(10):107. https://doi.org/10.3390/aerospace6100107
Chicago/Turabian StyleLappas, Ilias, and Akira Ikenaga. 2019. "Conceptual Design and Performance Optimization of a Tip Device for a Regional Turboprop Aircraft" Aerospace 6, no. 10: 107. https://doi.org/10.3390/aerospace6100107
APA StyleLappas, I., & Ikenaga, A. (2019). Conceptual Design and Performance Optimization of a Tip Device for a Regional Turboprop Aircraft. Aerospace, 6(10), 107. https://doi.org/10.3390/aerospace6100107