Utah State University DigitalCommons@USU
All Graduate Theses and Dissertations Graduate Studies
5-2015
Computational Fluid Dynamics Simulations of Oscillating Wings and Comparison to Lifting-Line Theory
Megan Keddington Utah State University
Follow this and additional works at: https://digitalcommons.usu.edu/etd
Part of the Aerospace Engineering Commons
Recommended Citation Keddington, Megan, "Computational Fluid Dynamics Simulations of Oscillating Wings and Comparison to Lifting-Line Theory" (2015). All Graduate Theses and Dissertations. 4473. https://digitalcommons.usu.edu/etd/4473
This Thesis is brought to you for free and open access by the Graduate Studies at DigitalCommons@USU. It has been accepted for inclusion in All Graduate Theses and Dissertations by an authorized administrator of DigitalCommons@USU. For more information, please contact [email protected].
COMPUTATIONAL FLUID DYNAMICS SIMULATIONS OF OSCILLATING
WINGS AND COMPARISON TO LIFTING-LINE THEORY
by
Megan Keddington
A thesis submitted in partial fulfillment of the requirements for the degree
of
MASTER OF SCIENCE
in
Aerospace Engineering
Approved:
______Warren F. Phillips Steven L. Folkman Major Professor Committee Member
______Douglas F. Hunsaker Mark R. McLellan Committee Member Vice President for Research and Dean of the School of Graduate Studies
UTAH STATE UNIVERSITY Logan, Utah
2015
ii
Copyright © Megan Keddington 2015
All Rights Reserved
iii
ABSTRACT
Computational Fluid Dynamics Simulations of Oscillating Wings
and Comparison to Lifting-Line Solution
by
Megan Keddington
Utah State University, 2015
Major Professor: Dr. Warren F. Phillips Department: Mechanical and Aerospace Engineering
Computational fluid dynamics (CFD) analysis was performed in order to compare the solutions of oscillating wings with Prandtl’s lifting-line theory. Quasi-steady and steady-periodic simulations were completed using the CFD software Star-CCM+. The simulations were performed for a number of frequencies in a pure plunging setup. Additional simulations were then completed using a setup of combined pitching and plunging at multiple frequencies. Results from the CFD simulations were compared to the quasi-steady lifting-line solution in the form of the axial-force, normal-force, power, and thrust coefficients, as well as the efficiency obtained for each simulation. The mean values were evaluated for each simulation and compared to the quasi-steady lifting-line solution. It was found that as the frequency of oscillation increased, the quasi-steady lifting-line solution was decreasingly accurate in predicting solutions.
(350 pages)
iv
PUBLIC ABSTRACT
Computational Fluid Dynamics Simulations of Oscillating
Wings and Comparison to Lifting-Line Theory
Megan Keddington
Computational fluid dynamics (CFD) analysis was performed in order to compare the solutions of oscillating wings with Prandtl’s lifting-line theory. Quasi-steady and steady-periodic simulations were completed using the CFD software Star-CCM+. Quasi-steady simulations were completed for both two- dimensional and three-dimensional setups. The steady-periodic simulations were only performed using a three-dimensional setup. The simulations were performed for nine separate frequencies in a pure plunging setup. An additional four simulations were then completed using a setup of combined pitching and plunging at four separate frequencies. Results from the CFD simulations were compared to the quasi-steady lifting-line solution in the form of the axial-force, normal-force, power, and thrust coefficients, as well as the efficiency obtained for each simulation. The mean values were evaluated for each simulation and compared to the quasi-steady lifting-line solution.
It was found that as the frequency of oscillation increased, the quasi-steady lifting-line solution was decreasingly accurate in predicting solutions. It was also observed that the thrust was generated only by plunging, not pitching of the wing in the simulations.
v
ACKNOWLEDGMENTS
I would like to thank Dr. Warren Phillips for imparting to me a small portion of his vast knowledge. He expected the highest, and for that, I am grateful. I would also like to thank Dr. Douglas
Hunsaker for his constant guidance, council, and assistance. I could have never completed this thesis without his support and encouragement. I am indebted to both Dr. Phillips and Dr. Hunsaker for the extensive work they have put into teaching me and helping me along in this process.
I am eternally grateful to my family, and especially my parents, for teaching me that I can do anything. Lastly, I’d like to thank my husband, for always being my rock, my champion, and my sanctuary.
Words cannot express how grateful I am.
Megan Keddington Jensen
vi
LIST OF TABLES
Table Page
1 Correlation coefficients for Eqs. (1.61) and (1.62) ...... 14
ω → ω = α = 2 Mean values of 2D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the exact solution ...... 60
ω → ω = α = 3 Mean values of 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the quasi-steady lifting-line solution ...... 68
ω → ω = α = 4 Mean values of 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution ...... 70
ω → ω = α = − 5 Mean values of 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution ...... 73
ω = ω = α = 6 Mean values of oscillation cycle with ˆ x .0 13414 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 76
ω = ω = α = 7 Mean values of oscillation cycle with ˆ x .0 14905 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 78
ω = ω = α = 8 Mean values of oscillation cycle with ˆ x .0 16162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 81
ω = ω = α = 9 Mean values of oscillation cycle with ˆ x .0 19163 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 83
ω = ω = α = 10 Mean values of oscillation cycle with ˆ x .0 20637 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 86
ω = ω = α = 11 Mean values of oscillation cycle with ˆ x .0 24841 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 88
ω = ω = α = 12 Mean values of oscillation cycle with ˆ x .0 29162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line theory ...... 91
ω = ω = α = 13 Mean values of oscillation cycle with ˆ x .0 33536 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 93
ω = ω = α = 14 Mean values of oscillation cycle with ˆ x .0 38327 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution ...... 96
vii
ω = ω = α = 15 Mean values of oscillation cycle with ˆ x .0 25797 , ˆ y .0 09365, and ˆ y 2.0 compared to the quasi-steady lifting-line solution ...... 99
ω = ω = α = 16 Mean values of oscillation cycle with ˆ x .0 47909 , ˆ y .0 09365 , and ˆ y 2.0 compared to the quasi-steady lifting-line solution ...... 101
ω = ω = α = − 17 Mean values of oscillation cycle with ˆ x .0 17198 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution ...... 104
ω = ω = α = − 18 Mean values of oscillation cycle with ˆ x .0 31939 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution ...... 107
viii
LIST OF FIGURES
Figure Page
1 Parameters used by Theodorsen [9] to describe plunging and pitching motion (from Hunsaker and Phillips [12] with permission) ...... 2
2 Example coarse grid for computational fluid dynamics calculations (from Hunsaker and Phillips [12] with permission) ...... 8
3 Comparison between the CFD solutions and the correlation with Eq. (1.61) (from Hunsaker and Phillips [12] with permission) ...... 14
4 Comparison between the CFD solutions and the correlation with Eq. (1.62) (from Hunsaker and Phillips [12] with permission) ...... 14
5 Comparison of the mean axial-force coefficient predicted from the CFD solutions, Eq. (1.57), and the model first presented by Theodorsen [9] (from Hunsaker and Phillips [12] with permission) ...... 15
6 Comparison of the mean required-power coefficient predicted from the CFD solutions, Eq. (1.56), and the model first presented by Theodorsen [9] (from Hunsaker and Phillips [12] with permission) ...... 15
7 Comparison of the propulsive efficiency predicted from the CFD solutions, Eq. (1.58), and the model first presented by Theodorsen [9] (from Hunsaker and Phillips [12] with permission) ...... 15
8 Constant-i planes for a typical C-O grid (from Phillips, Fugal, and Spall [253] with permission) ...... 32
9 Constant-j planes for a typical C-O grid (from Phillips, Fugal, and Spall [253] with permission) ...... 32
10 Constant-k planes for a typical C-O grid (from Phillips, Fugal, and Spall [253] with permission) ...... 33
11 Example of a constant-i plane at the trailing edge of the wingtip and constant-j plane on the wing surface of a three-dimensional grid (from Phillips, Fugal, and Spall [253] with permission) ...... 35
12 Example of a constant-k plane at the base of the end cap and constant-j plane on the end cap surface of a three-dimensional grid (from Phillips, Fugal, and Spall [253] with permission)...... 35
13 Constant-k plane at the base of the wing cap and constant-j plane on the wing surface near the leading edge of the three-dimensional coarse grid...... 36
14 Coarse grid for two-dimensional CFD calculations ...... 36
15 Overset grid compared to the background grid ...... 37
ix
16 Star-CCM+ farfield background region ...... 38
17 Star-CCM+ symmetry plane and outlet background regions with the airfoil foreground region visible ...... 39
18 Star-CCM+ foreground regions ...... 39
19 Star-CCM+ foreground and background regions ...... 40
20 Comparison of two-dimensional coarse-, medium-, and fine-grids with Richardson Extrapolation to exact lift coefficient ...... 42
21 Coarse-, medium-, and fine-grids compared to the quasi-steady lifting-line solution for the three-dimensional background grid ...... 43
22 Lift coefficient when rotating the overset grid ...... 44
23 Lift coefficient when rotating the freestream velocity ...... 44
24 Lift coefficient when imparting a velocity to the overset grid ...... 45
25 Error for the Richardson Extrapolation for each angle of attack defining technique ...... 45
26 Comparison of the error in the lift coefficient of the closest overset grids as well as the overset grid chosen ...... 46
27 Lift coefficient of the coarse-, medium-, and fine-grid CFD solutions compared to Richardson Extrapolation solution ...... 47
28 Drag coefficient of the coarse-, medium-, and fine-grid CFD solutions compared to Richardson Extrapolation solution ...... 47
29 Error in the axial-force coefficient of the coarse-, medium-, and fine-grid CFD solutions compared to Richardson Extrapolation solution ...... 48
30 Error in the normal-force coefficient of the course-,medium-, and fine-grid CFD solutions compared to Richardson Extrapolation Solution ...... 48
31 Time step resolution shown with the comparison of the power coefficient ...... 53
32 Time step resolution shown with the comparison of the thrust coefficient ...... 54
33 Cycle resolution shown with the comparison of the power coefficient...... 55
34 Cycle resolution shown with the comparison of the thrust coefficient ...... 55
ω → ω = 35 Axial-force coefficient of the 2D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 0.0 compared with the exact solution ...... 58
ω → ω = 36 Normal-force coefficient of the 2D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 0.0 compared with the exact solution ...... 58
x
ω → ω = α = 37 Power coefficient of the 2D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the exact solution ...... 59
ω → ω = α = 38 Thrust coefficient of the 2D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the exact solution ...... 59
ω → ω = 39 Axial-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 0.0 compared with the quasi-steady lifting-line solution, Eq. (3.23) ...... 66
ω → ω = 40 Normal-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 0.0 compared with the quasi-steady lifting-line solution, Eq. (3.28) ...... 66
ω → ω = α = 41 Power coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the quasi-steady lifting-line solution, Eq. (3.33) ...... 67
ω → ω = α = 42 Thrust coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 0.0 compared with the quasi-steady lifting-line solution, Eq. (3.40) ...... 67
ω → ω = 43 Axial-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.23) ...... 68
ω → ω = 44 Normal-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.28) ...... 69
ω → ω = α = 45 Power coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.33) ...... 69
ω → ω = α = 46 Thrust coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.40) ...... 70
ω → ω = 47 Axial-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = − ˆ y .0 02 compared with the quasi-steady lifting-line solution, Eq. (3.23) ...... 71
ω → ω = 48 Normal-force coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and α = − ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.28) ...... 71
ω → ω = α = − 49 Power coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.33) ...... 72
ω → ω = α = − 50 Thrust coefficient of the 3D quasi-steady cycle ˆ x 0 , ˆ y .0 07492 , and ˆ y 2.0 compared with the quasi-steady lifting-line solution, Eq. (3.40) ...... 72
xi
ω = ω = 51 Axial-force coefficient of oscillation cycle with ˆ x .0 13414 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 74
ω = ω = 52 Normal-force coefficient of oscillation cycle with ˆ x .0 13414 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 74
ω = ω = α = 53 Power coefficient of oscillation cycle with ˆ x .0 13414 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 75
ω = ω = α = 54 Thrust coefficient of oscillation cycle with ˆ x .0 13414 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 75
ω = ω = 55 Axial-force coefficient of oscillation cycle with ˆ x .0 14905 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 76
ω = ω = 56 Normal-force coefficient of oscillation cycle with ˆ x .0 14905 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 77
ω = ω = α = 57 Power coefficient of oscillation cycle with ˆ x .0 14905 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 77
ω = ω = α = 58 Thrust coefficient of oscillation cycle with ˆ x .0 14905 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 78
ω = ω = 59 Axial-force coefficient of oscillation cycle with ˆ x .0 16162 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 79
ω = ω = 60 Normal-force coefficient of oscillation cycle with ˆ x .0 16162 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 79
ω = ω = α = 61 Power coefficient of oscillation cycle with ˆ x .0 16162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 80
ω = ω = α = 62 Thrust coefficient of oscillation cycle with ˆ x .0 16162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 80
ω = ω = 63 Axial-force coefficient of oscillation cycle with ˆ x .0 19163 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 81
ω = ω = 64 Normal-force coefficient of oscillation cycle with ˆ x .0 19163 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 82
xii
ω = ω = α = 65 Power coefficient of oscillation cycle with ˆ x .0 19163 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 82
ω = ω = α = 66 Thrust coefficient of oscillation cycle with ˆ x .0 19163 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 83
ω = ω = 67 Axial-force coefficient of oscillation cycle with ˆ x .0 20637 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 84
ω = ω = 68 Normal-force coefficient of oscillation cycle with ˆ x .0 20637 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 84
ω = ω = α = 69 Power coefficient of oscillation cycle with ˆ x .0 20637 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 85
ω = ω = α = 70 Thrust coefficient of oscillation cycle with ˆ x .0 20637 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 85
ω = ω = 71 Axial-force coefficient of oscillation cycle with ˆ x .0 24841 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 86
ω = ω = α = 72 Normal-force coefficient of oscillation cycle with ˆ x .0 24841 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 87
ω = ω = α = 73 Power coefficient of oscillation cycle with ˆ x .0 24841 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 87
ω = ω = α = 74 Thrust coefficient of oscillation cycle with ˆ x .0 24841 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 88
ω = ω = 75 Axial-force coefficient of oscillation cycle with ˆ x .0 29162 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 89
ω = ω = 76 Normal-force coefficient of oscillation cycle with ˆ x .0 29162 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 89
ω = ω = α = 77 Power coefficient of oscillation cycle with ˆ x .0 29162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 90
ω = ω = α = 78 Thrust coefficient of oscillation cycle with ˆ x .0 29162 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 90
xiii
ω = ω = 79 Axial-force coefficient of oscillation cycle with ˆ x .0 33536 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 91
ω = ω = 80 Normal-force coefficient of oscillation cycle with ˆ x .0 33536 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 92
ω = ω = α = 81 Power coefficient of oscillation cycle with ˆ x .0 33536 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 92
ω = ω = α = 82 Thrust coefficient of oscillation cycle with ˆ x .0 33536 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 93
ω = ω = 83 Axial-force coefficient of oscillation cycle with ˆ x .0 38327 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 94
ω = ω = 84 Normal-force coefficient of oscillation cycle with ˆ x .0 38327 , ˆ y .0 07492 , and α = ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 94
ω = ω = α = 85 Power coefficient of oscillation cycle with ˆ x .0 38327 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 95
ω = ω = α = 86 Thrust coefficient of oscillation cycle with ˆ x .0 38327 , ˆ y .0 07492 , and ˆ y 0.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 95
ω = ω = 87 Axial-force coefficient of oscillation cycle with ˆ x .0 25797 , ˆ y .0 09365, and α = ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 97
ω = ω = 88 Normal-force coefficient of oscillation cycle with ˆ x .0 25797 , ˆ y .0 09365, and α = ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 97
ω = ω = α = 89 Power coefficient of oscillation cycle with ˆ x .0 25797 , ˆ y .0 09365, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 98
ω = ω = α = 90 Thrust coefficient of oscillation cycle with ˆ x .0 25797 , ˆ y .0 09365, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 98
ω = ω = 91 Axial-force coefficient of oscillation cycle with ˆ x .0 47909 , ˆ y .0 09365 , and α = ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 99
ω = ω = 92 Normal-force coefficient of oscillation cycle with ˆ x .0 47909 , ˆ y .0 09365 , and α = ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 100
xiv
ω = ω = α = 93 Power coefficient of oscillation cycle with ˆ x .0 47909 , ˆ y .0 09365 , and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 100
ω = ω = α = 94 Thrust coefficient of oscillation cycle with ˆ x .0 47909 , ˆ y .0 09365 , and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 101
ω = ω = 95 Axial-force coefficient of oscillation cycle with ˆ x .0 17198 , ˆ y .0 06243, and α = − ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 102
ω = ω = 96 Normal-force coefficient of oscillation cycle with ˆ x .0 17198 , ˆ y .0 06243, and α = − ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 102
ω = ω = α = − 97 Power coefficient of oscillation cycle with ˆ x .0 17198 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 103
ω = ω = α = − 98 Thrust coefficient of oscillation cycle with ˆ x .0 17198 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 103
ω = ω = 99 Axial-force coefficient of oscillation cycle with ˆ x .0 31939 , ˆ y .0 06243, and α = − ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.23) ...... 105
ω = ω = 100 Normal-force coefficient of oscillation cycle with ˆ x .0 31939 , ˆ y .0 06243, and α = − ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.28) ...... 105
ω = ω = α = − 101 Power coefficient of oscillation cycle with ˆ x .0 31939 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.33) ...... 106
ω = ω = α = − 102 Thrust coefficient of oscillation cycle with ˆ x .0 31939 , ˆ y .0 06243, and ˆ y 2.0 compared to the quasi-steady lifting-line solution, Eq. (3.40) ...... 106
103 Comparison of the CFD mean axial-force coefficients with quasi-steady lifting-line theory ...... 107
104 Comparison of the CFD mean normal-force coefficients with quasi-steady lifting-line theory ...... 108
105 Comparison of the CFD mean power coefficients with quasi-steady lifting-line theory ...... 108
106 Comparison of the CFD mean thrust coefficients with the quasi-steady lifting-line theory ..... 109
107 Comparison of the CFD mean efficiency with the lifting-line theory ...... 109
xv
CONTENTS
ABSTRACT ...... iii
PUBLIC ABSTRACT ...... iv
ACKNOWLEDGMENTS ...... v
LIST OF TABLES ...... vi
LIST OF FIGURES ...... viii
NOMENCLATURE ...... xvi
CHAPTER
I. INTRODUCTION ...... 1
TWO-DIMENSIONAL STUDIES ...... 1 THREE-DIMENSIONAL STUDIES ...... 16 PURPOSE OF THE PRESENT RESEARCH ...... 27
II. COMPUTATIONAL FLUID DYNAMICS METHODOLOGY...... 31
SOFTWARE SELECTION ...... 33 GRID GENERATION ...... 34 QUASI-STEADY SETUP ...... 38 GRID CONVERGENCE AND RESOLUTION...... 40 STEADY-PERIODIC CASE SETUP ...... 49 STEADY-PERIODIC CONVERGENCE ...... 52
III. RESULTS ...... 57
QUASI-STEADY RESULTS ...... 57 STEADY-PERIODIC RESULTS ...... 73 MEAN RESULTS ...... 107
IV. CONCLUSION AND RECOMMENDATION FOR FUTURE WORK ...... 111
CONCLUSION ...... 111 RECOMMENDATION FOR FUTURE WORK ...... 113
REFERENCES ...... 114
APPENDIX A: TWO-DIMENSIONAL SCRIPTS...... 129
APPENDIX B: THREE-DIMENSIONAL SCRIPTS ...... 140
APPENDIX C: REFERENCES, KEYWORDS AND ABSTRACTS ...... 158
APPENDIX D: TABLES OF DATA USED FOR PLOTTING ...... 226
xvi
NOMENCLATURE
A = Fourier coefficient defined in Eq. (1.47)
An = coefficients in the infinite series solution
Ax0− x2 = coefficients defined in Eq. (1.50) Ay1 = coefficient defined in Eq. (1.50)
A12−32 = correlation coefficients used in Eq. (1.61) and defined in Table 1 a = axial location of the axis of rotation nondimensionalized with respect to the airfoil half chord an = planform coefficients in the decomposed Fourier-series solution to the lifting-line equation B = Fourier coefficient defined in Eq. (1.47)
Bx1−x2 = coefficient defined in Eq. (1.50) By1 = coefficient defined in Eq. (1.50)
B12−30 = correlation coefficients used in Eq. (1.62) and defined in Table 1 b = airfoil section half chord bn = washout coefficients in the decomposed Fourier-series solution to the lifting-line equation C = complex Theodorsen function given in Eq. (1.4)
CD = drag coefficient C = induced drag coefficient Di C = mean wing induced-drag coefficient averaged over a complete flapping cycle Di C = wing parasitic-drag coefficient DP CL = instantaneous section lift coefficient
CL = mean section lift coefficient
CL,α = wing lift slope C = lift coefficient of the first, second, and third grid ~L 1, −3 CL,α = airfoil-section lift slope C = section pitching moment coefficient about the airfoil aerodynamic center mac C = section pitching moment coefficient about the airfoil quarter chord mc 4/ C = mean section pitching moment coefficient about the quarter chord mc 4/ CP = instantaneous section power-required coefficient CP = mean section power-required coefficient C = mean available-power coefficient averaged over a complete flapping cycle PA C = instantaneous flapping-power coefficient required to support the flapping rate Pf C = mean flapping-power coefficient averaged over a complete flapping cycle Pf CT = instantaneous thrust coefficient
CT = mean thrust coefficient Cx = instantaneous section axial-force coefficient
Cx = mean section axial-force coefficient
Cx 0 = axial-force coefficient at the zero plunging velocity point in the cycle C y = instantaneous section normal-force coefficient
C y = mean section normal-force coefficient
C0−∞ = coefficients in the Taylor series expansion given in Eq. (1.44) c = airfoil section chord length cn = control coefficients in the decomposed Fourier-series solution to the lifting-line equation dn = plunging coefficients in the decomposed Fourier-series solution to the lifting-line equation en = plunging coefficients in the Fourier-series solution for induced drag F = real part of the complex Theodorsen function given in Eq. (1.5)
xvii
Fx = section axial force Fy = section normal force G = imaginary part of the complex Theodorsen function given in Eq. (1.6) g, g1 3 = cell size for Richardson extrapolation, defined in Eq. (2.3) h, ha = complex plunging displacement at the axis of rotation (positive downward) hc/4 = complex plunging displacement at the quarter chord (positive downward) h3c/4 = complex plunging displacement at the three-quarter chord (positive downward) h&a = complex plunging velocity at the axis of rotation h&& a = complex plunging acceleration at the axis of rotation h&c 4/ = complex plunging velocity at the quarter chord h&&c 4/ = complex plunging acceleration at the quarter chord h&3c 4/ = complex plunging velocity at the three-quarter chord i = square root of −1
J0, J1 = Bessel functions k = reduced frequency used by Theodorsen based on the airfoil half chord, defined in Eq. (1.7) L = instantaneous section lift L = mean section lift
Ma, Mc/4 = complex section pitching moment about the axis of rotation and the quarter chord mac, mc/4 = section pitching moment about the aerodynamic center and the quarter chord mc 4/ = mean section pitching moment about the quarter chord = number of cells P = complex section lift force (positive downward) p = apparent order of convergence, defined in Eq. (2.7)
Pf = instantaneous power required to support the instantaneous flapping rate pˆ = dimensionless angular wing-tip rotation rate about the roll axis pˆ steady = steady dimensionless rigid-body rolling rate, which produces a zero net rolling rate pˆ rms = root mean square of the dimensionless flapping rate averaged over a complete flapping cycle Q = variable used by Theodorsen, defined in Eq. (1.1) q = logarithm of the difference of the grid ratios, defined in Eq. (2.11)
RA = aspect ratio of the wing
RAL = lift amplitude ratio, defined in Eq. (1.26) RAm = quarter-chord pitching-moment amplitude ratio, defined in Eq. (1.27)
RT = wing taper ratio
Rτ = ratio of the current case to the original case
Ry1 = amplitude of the sinusoidal terms Ay1 and By1 r21 = ratio of the second grid to the first grid r32 = ratio of the third grid to the second grid S = wing planform area s = sign of the ratio of the lift coefficients, defined in Eq. (2.12) T = section thrust
TA1−2 = correlation coefficients, used in Eq. (1.61) and defined in Table 1
TB1−2 = correlation coefficients, used in Eq. (1.62) and defined in Table 1 t = time V = freestream airspeed
Vy = y-velocity component of the airfoil aerodynamic center (positive upward)
xviii
Vs = start velocity
Vt = translational velocity x, y, z = streamwise, upward normal, and spanwise coordinates relative to the section quarter chord
Y0, Y1 = Bessel functions y A = aerodynamic-center plunging amplitude y = original aerodynamic-center plunging amplitude A0 yac = aerodynamic-center plunging displacement α = pitching angle of attack α& = angular velocity α&& = angular acceleration α = mean pitching angle
α A = pitching amplitude
α g = geometric angle of attack
α i = induced angle of attack
α L0 = zero-lift angle of attack
αP = plunging angle of attack
αroot = root angle of attack
αˆ y = aerodynamic-center geometric angle-of-attack ratio, defined in Eq. (1.22) Γ = wing section circulation ~ Di = local section induced-drag increments