Skin Friction Drag and Shear Stress Distribution on Several Streamlined Bodies of Revolution with Varied Fineness Ratio

Total Page:16

File Type:pdf, Size:1020Kb

Skin Friction Drag and Shear Stress Distribution on Several Streamlined Bodies of Revolution with Varied Fineness Ratio Calhoun: The NPS Institutional Archive Theses and Dissertations Thesis Collection 1964 Skin Friction Drag and Shear Stress Distribution on Several Streamlined Bodies of Revolution With Varied Fineness Ratio. Dunning, Cleveland L. Monterey, California. Naval Postgraduate School http://hdl.handle.net/10945/26252 NPS ARCHIVE 1964 DUNNING, C. RKIN FRICTION DRAG AND SMEAR STRESS DISTRIBUTION ON SEVERAL STREAMLINED BODIES OF REVOLUTION WITH VARIED FINENESS RATIO E v'". .i VI LAND L PUNNING ROBERT F. OVERMYER COLBEN K. SIME, JR. is s A ffWmmM' LIBRARY U.S. NAVAL POSTGRADUATE SCHOOL MONTEREY, CALIFORNIA SKIN FRICTION DRAG AND SHEAR STRESS DISTRIBUTION Oil SEVERAL STREAMLINED BODIES OF REVOLUTION WITH VARIED FINENESS RATIO 3)C >fC >f< 3fC 3fi Cleveland L. Dunning Robert F. Overmyer and Colben K. Sime, Jr. f t mm* # LIBRARY OS. NAVAL POSTGRADUATE SCHOOL MONTEREY, CALIFORNIA SKIN FRICTION DRAG AND SHEAR STRESS DISTRIBUTION ON SEVERAL STREAMLINED BODIES OF REVOLUTION WITH VARIED FINENESS RATIO by Cleveland L. Dunning Captain, United States Marine Corps Robert F. Overmyer Captain, United States Marine Corps and Colben K. Sime, Jr. Captain, United States Marine Corps Submitted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE with major in Aeronautics United States Naval Postgraduate School Monterey, California 1964 SKIN FRICTION DRAG AND SHEAR STRESS DISTRIBUTION ON SEVERAL STREAMLINED BODIES OF REVOLUTION WITH VARIED FINENESS RATIO by Cleveland L. Dunning Robert F. Overmyer and Colben K. Sime, Jr. This thesis is accepted in partial fulfillment for the degree of MASTER OF SCIENCE with major in Aeronautics from the United States Naval Postgraduate School ACKNOWLEDGEMENT The authors wish to express their sincere gratitude to Professor George J. Higgins for his great interest, help and encouragement in the completion of this work. Professor Higgins., in addition to his full schedule as a faculty member and on short notice, proposed this project and has graciously provided time and guidance toward the fulfillment of this objective. In addition, many stimulating discussions by Professor Higgins in previous courses in fluid mechanics and aerodynamics are acknowledged. 11 ABSTRACT Experimental techniques were employed to determine the coefficient of skin friction drag and the distribution of local frictional intensity for several low drag bodies of revolution with mixed laminar and turbulent flow The data were obtained as a function of fineness ratio and Reynolds number from experiments conducted in the U. S. Naval Postgraduate School subsonic wind tunnel under conditions equivalent to incomDressible flow. Skin friction drag coefficients were determined first by the indirect method of subtracting pressure drag from total drag. A second more direct method was the deter- mination of local shear distributions and integration over the surface area to obtain total frictional drag* The possibility of inaccuracies introduced by assumptions and agreement between results of the two methods are discussed 111 TABLE OF CONTENTS Section Title Page 1. Introduction 1 2. Equipment and Procedures 2 a. Total Drag Measurements 3 b. Support Drag Lf c. Local Pressure Distribution Measurements 5 d. Boundary Layer Probe 6 3. Presentation of Results 3 h. Skin Friction Determined From Total Drag Minus Pressure Drag 8 a. Total Drag 3 b. Pressure Drag 13 c. Skin Friction Drag 18 5. Direct Determination of Local Skin Friction Intensity 21 a. Analytical and Theoretical Methods 26 b. Experimental Estimate of Shear Stresses 30 6. Comparison of Results 38 7. Conclusion *+0 References *+3 Tables ^5 Illustrations 9^ IV LIST OF TABLES Table Page I Diameter-Length Coordinates, Streamlined Bodies, 66(215) - OXX h5 II Model Characteristics, 66(215) - CXX 1+6 III Pressure Orifice Locations, Streamlined Bodies, 66(215) - OXX h7 IV Coefficient of Drag Streamlined Body of Revolution, L/D 1 hS V Coefficient of Drag Streamlined Body of Revolution, L/D = 2 h9 VI Coefficient of Drag Streamlined Body of Revolution, L/D = 3 50 VII Coefficient of Drag Streamlined Body of Revolution, L/D = h 51 VIII Coefficient of Drag Streamlined Body of Revolution, L/D = 5 52 IX Coefficient of Drag Streamlined Body of Revolution, L/D = 6 53 X Distribution of Normal Pressure Forces and Local Velocities, L/D = 1 5h XI Distribution of Normal Pressure Forces and Local Velocities, L/D = 2 59 XII Distribution of Normal Pressure Forces and Local Velocities, L/D = 3 6k XIII Distribution of Normal Pressure Forces and Local Velocities, L/D = h 69 XIV Distribution of Normal Pressure Forces and Local Velocities, L/D = 5 7h XV Distribution of Normal Pressure Forces and Local Velocities, L/D = 6 79 XVI Lengthwise Velocity Distribution In The Boundary Layer At 0.009 Inches From Surface Table Page XVII Experimentally Estimated Skin Friction, Body of Revolution, L/D =3 36 XVIII Experimentally Estimated Skin Friction, Body of Revolution, L/D = h 88 XIX Experimentally Estimated Skin Friction, Body of Revolution, L/D =5 90 XX Experimentally Estimated Skin Friction, Body of Revolution, L/D =6 92 vi LIST OF ILLUSTRATIONS Figure Page 1. The Streamlined Bodies Tested 9^ 2. Schematic Diagram of the U. So Naval Post- graduate School Closed Circuit Wind Tunnel 95 3. Aerolab Three Component Eeam Balance 96 ha. Equipment Set Up For Total Drag Measurement 97 hl>. Side View of Model And Test Section 98 kc. Front View of Model and Test Section 98 5. Model With Dummy SupDort Strut 99 6. Equipment Set Up For Pressure Measurement 99 7. Drag Coefficients For Bodies of Revolution, CDwet Versus RL , 100 8. Orientation of Mutually Perpendicular Pressure and Frictional Forces 101 9. Distribution of Normal Pressure And Inte- gration of Pressure Drag On A Body of Revolution 102 10. Pressure Drag Coefficients For Bodies Of Revolution, Cd versus R^ 108 11. Skin Friction Drag Coefficients For Bodies Of Revolution, C - 109 Df = CDwet Cop 12. Percentage of Friction Drag to Total Drag Versus The Ratio of Diameter Over Length 110 13. Comparison of Laminar Boundary Layer Profile Slopes on Bodies of Revolution with a Flat Plate 111 Ik, Pressure Distribution Over Bodies of Revolution = L/D 1, 2, 3, *+i 5, and 6 112 15. Lengthwise Velocity Distribution in The Boundary Layer of Streamlined Bodies of Revolution, L/D = 3, h, 118 vi 1 Figure Page 16. Lengthwise Velocity Distribution in The Boundary Layer of Streamlined Bodies of Revolution, L/D = 5, 6. 119 17. Estimated Lengthwise Variation of Non- dimensional Shear Stress on a Body of Revolution, L/D = 3. 120 18. Estimated Lengthwise Variation of Non- dimensional Shear Stress on a Body of Revolution, L/D = k. 121 19. Estimated Lengthwise Variation of Non- dimensional Shear Stress on a Body of Revolution, L/D = 5. 122 20. Estimated Lengthwise Variation of Non- dimensional Shear Stress on a Body of Revolution, L/D = 6. 123 21. Skin Friction Drag Coefficients From the Integration of Local Frictional Shear Stress* 12m- 22. Comparison of Estimated Friction Drag Coef- ficients with Differences Between Total and Pressure Drag Coefficients, L/D = 3, 5- i2 5 23. Comparison of Estimated Friction Drag Coef- ficients With Differences Between Total and Pressure Drag Coefficients, L/D - M-, 6 C 126 viii TABLE OF SYKBCLS Symbol Definition Units Skin friction drag coefficient based on wetted surface area CDwet Total drag coefficient based on vetted surface area CDp Pressure drag coefficient based on wetted surface area CDt Coefficient of total drag - includes support and buoyancy drag ^Dsupp Coefficient of support drag Cf Flat plate drag coefficient Cflam Laminar flat plate drag coefficient cfturb Turbulent fla>t plate drag coefficient cftrans Flat plate drag coefficient at transition Coefficient of buoyancy drag ^CDs+B Coefficient of buoyancy plus supnort drag Df Friction drag •*' lbs Dt Total drag lbs„ D Diameter fto L Axial length fto 2 P Local static pressure Ibo/fto Po Stream static pressure lb /ft AP Local static minus stream static pressure Ibc/fto 2 Qo Free stream dynamic pressure Ibo/fto r Radius . ft c RL Reynolds number based on axial length Local Reynolds number based on developed length ix Symbol Definition Units s Variable developed length ft. S Developed length fto S Frontal area ft. 2 Sw Wetted surface area fto V Local velocity f to/sec. V Free stream velocity f to/seco X Variable axial length fto (r£) gradient surface f e Velocity at to/seCo fto m Slone of non-dimensional velocity profile 1/ft1/fto 2 f> Density of air lb s =<sec a fteA Angle between local tangent to surface and body axis dego ^ Absolute viscosity lbo°sec *° Local shear stress at the surface lbo/fto 2 1. Introduction, The aerodynamic forces on a body immersed in an incom- pressible fluid flow consist of contributions from normal pressures and tangential shear stresses The problem of de= termining the coefficient of skin friction drag and distri- bution of local frictional intensity for several bodies of revolution as a function of fineness ratio and Reynolds number was undertaken at the United States Naval Postgraduate School, Monterey, California . Data for this analysis were obtained from subsonic wind tunnel tests on six models with fineness ratios from one to six at zero angle of attack and through a Reynolds number range of Ool to 5°0 x 10 . The model shapes were based upon the NACA 66(215) -0XX airfoilo The skin friction drag based on wetted surface area was to be initially determined from the difference of the total drag and pressure drag. A second method involving the inte- gration of the local shear stresses in the boundary layer over the developed length of the body was later utilized to deter- mine the skin friction drag. The expected agreement in results from the two methods was dependent upon an accurate determination of the drag com- ponents of the mutually perpendicular pressure and shearing forces acting normal and tangential to the model surface c It was hoped that these results could be correlated with and com- pared to other analytically and experimentally determined skin friction drag coefficients of a flat plate at zero incidence and other axially symmetric bodies* 2.
Recommended publications
  • Prediction of Laminar/Turbulent Transition in Airfoil Flows
    PREDICTION OF LAMINAR/TURBULENT TRANSITION IN AIRFOIL FLOWS by Jeppe Johansen and Jens Norkaer Sorensen DTU FLUID MECHANICS H ENERGY ENGINEERING TECHNICAL UNIVERSITY OF DENMARK / DANISH CENTER FOR APPLIED MATHEMATICS AND MECHANICS Scientific Council Poul Andersen Dept, of Naval Architecture and Offshore Engineering Martin P. Bends0e Dept, of Mathematics Ove Ditlevsen Dept, of Structural Engineering and Materials Ivar G. Jonsson Dept, of Hydrodynamics and Water Resources Wolfhard Kliem Dept, of Mathematics Steen Krenk Dept, of Structural Engineering and Materials P. Scheel Larsen Dept, of Energy Engineering Frithiof I. Niordson Dept, of Solid Mechanics Pauli Pedersen Dept, of Solid Mechanics P. Temdrup Pedersen Dept, of NavalArchitecture and Offshore Engineering- Dan Rosbjerg Dept. of Hydrodynamics and Water Resources Jens Nprkser Sprensen Dept, of Energy Engineering P. Grove Thomsen Dept, of Mathematical Modelling Hans True Dept, of Mathematical Modelling Viggo Tvergaard Dept, of Solid Mechanics Secretary Pauli Pedersen, Docent, Dr.techn. Department of Solid Mechanics, Building 404 Technical University of Denmark DK-2800 Lyngby, Denmark DISCLAIMER Portions of this document may be illegible in electronic image products. Images are produced from the best available original document. Prediction of Laminar/Turbulent Transition in Airfoil Flows Jeppe Johansen* and Jens N. Sprensen* * Ris0 National Laboratory, Denmark and * Department of Energy Engineering, Technical University of Denmark Presented as Paper 98-0702 at the AIAA 36th Aerospace Sciences Meeting Sc Exhibit, Reno, NV, January 12*15, 1998 tphD student, Wind Energy and Atmospheric Physics Department, Rise National Laboratory, DK-4000 RoskUde, Denmark * Associate Professor, Department of Energy Engineering, Technical University of Denmark, DK- 2800 Lyngby, Denmark 1 Abstract The prediction of the location of transition is important for low Reynolds number airfoil flows.
    [Show full text]
  • Wingtip Vortices and Free Shear Layer Interaction in The
    WINGTIP VORTICES AND FREE SHEAR LAYER INTERACTION IN THE VICINITY OF MAXIMUM LIFT TO DRAG RATIO LIFT CONDITION Dissertation Submitted to The School of Engineering of the UNIVERSITY OF DAYTON In Partial Fulfillment of the Requirements for The Degree of Doctor of Philosophy in Engineering By Muhammad Omar Memon, M.S. UNIVERSITY OF DAYTON Dayton, Ohio May, 2017 WINGTIP VORTICES AND FREE SHEAR LAYER INTERACTION IN THE VICINITY OF MAXIMUM LIFT TO DRAG RATIO LIFT CONDITION Name: Memon, Muhammad Omar APPROVED BY: _______________________ _______________________ Aaron Altman Markus Rumpfkeil Advisory Committee Chairman Committee Member Professor; Director, Graduate Aerospace Program Associate Professor Mechanical and Aerospace Engineering Mechanical and Aerospace Engineering _______________________ _______________________ Jose Camberos Wiebke S. Diestelkamp Committee Member Committee Member Adjunct Professor Professor & Chair Mechanical and Aerospace Engineering Department of Mathematics _______________________ _______________________ Robert J. Wilkens, PhD., P.E. Eddy M. Rojas, PhD., M.A., P.E. Associate Dean for Research and Innovation Dean, School of Engineering Professor School of Engineering ii © Copyright by Muhammad Omar Memon All rights reserved 2017 iii ABSTRACT WINGTIP VORTICES AND FREE SHEAR LAYER INTERACTION IN THE VICINITY OF MAXIMUM LIFT TO DRAG RATIO LIFT CONDITION Name: Memon, Muhammad Omar University of Dayton Advisor: Dr. Aaron Altman Cost-effective air-travel is something everyone wishes for when it comes to booking flights. The continued and projected increase in commercial air travel advocates for energy efficient airplanes, reduced carbon footprint, and a strong need to accommodate more airplanes into airports. All of these needs are directly affected by the magnitudes of drag these aircraft experience and the nature of their wingtip vortex.
    [Show full text]
  • ANGLE of ATTACK (A of A)
    ANGLE OF ATTACK (A of A) The angle at which the airfoil meets the relative airflow is called the Angle of Attack. CENTER OF PRESSURE (C of P) If we consider all the distributed pressures to be equivalent to a single force, this force will act through a straight line. The point where this line cuts the chord of an airfoil is called the Center of Pressure. As the angle of attack is increased lift and drag increase and the Center of Pressure moves forward. This situation continues until the stall is reached. Beyond this point, it will move back. The movement of the Center of Pressure causes an airplane to be unstable. THE BOUNDARY LAYER The boundary layer is a very thin sheet of air lying over the surface of the wing and all other surfaces of the airplane. Because air has viscosity, this layer tends to stick to the wing. As the wing moves forward through the air the boundary layer at first flows smoothly over streamlined shape of the airfoil. Here the flow is called the Laminar Layer. As the boundary layer approaches the center of the wing it begins to lose speed due to skin friction and it becomes thicker and turbulent (turbulent layer). The point at which the boundary layer changes from laminar to turbulent is called the Transition Point. Where the boundary layer becomes turbulent, drag due to skin friction is relatively high. AIRFOIL DESIGNS There are many different airfoil designs. The type of operation for which an airplane is intended has a direct influence on the design and shape of the wing for that airplane.
    [Show full text]
  • Modeling Boundary-Layer Transition in DNS and LES Using Parabolized Stability Equations
    Center for Turbulence Research 29 Annual Research Briefs 2016 Modeling boundary-layer transition in DNS and LES using parabolized stability equations By A. Lozano-Dur´an, M. J. P. Hack, G. I. Park AND P. Moin 1. Motivation and objectives The modeling of the laminar-turbulent transition remains one of the key challenges in the numerical simulation of boundary layers, especially at coarse grid resolutions. The issue is particularly relevant in wall-modeled large-eddy simulations (LES), which require 10 to 100 times more grid points in the thin laminar region than in the turbulent regime to properly capture the instabilities preceding the transition (Slotnick et al. 2014). Our study examines the potential of the nonlinear parabolized stability equations (PSE) to provide an accurate yet computationally efficient treatment of the disturbances in the pre-transitional zone. Direct numerical simulation (DNS) of the Navier-Stokes equations has been frequently used as a tool to investigate transitional flows. However, its computational cost is un- affordable in most practical settings. In order to explore more computationally efficient approaches, Sayadi & Moin (2012) conducted LES of K- and H-type transitional bound- ary layers. They found that constant coefficient models for the subgrid-scale stress tensor could not predict the rapid rise in skin friction at the onset of transition. The reasons were traced back to the non-negligible turbulent viscosity in the laminar region, which dampens the amplification of instabilities. Dynamic models sufficiently reduced the tur- bulent viscosity in the laminar flow and allowed the growth of disturbances. When the grid was fine enough, LES reproduced the skin-friction overshoot observed in DNS.
    [Show full text]
  • Two-Dimensional Aerodynamics
    Flightlab Ground School 2. Two-Dimensional Aerodynamics Copyright Flight Emergency & Advanced Maneuvers Training, Inc. dba Flightlab, 2009. All rights reserved. For Training Purposes Only Our Plan for Stall Demonstrations Remember Mass Flow? We’ll do a stall series at the beginning of our Remember the illustration of the venturi from first flight, and tuft the trainer’s wing with yarn your student pilot days (like Figure 1)? The before we go. The tufts show complex airflow major idea is that the flow in the venturi and are truly fun to watch. You’ll first see the increases in velocity as it passes through the tufts near the root trailing edge begin to wiggle narrows. and then actually reverse direction as the adverse pressure gradient grows and the boundary layer The Law of Conservation of Mass operates here: separates from the wing. The disturbance will The mass you send into the venturi over a given work its way up the chord. You’ll also see the unit of time has to equal the mass that comes out movement of the tufts spread toward the over the same time (mass can’t be destroyed). wingtips. This spanwise movement can be This can only happen if the velocity increases modified in a number of ways, but depends when the cross section decreases. The velocity is primarily on planform (wing shape as seen from in fact inversely proportional to the cross section above). Spanwise characteristics have important area. So if you reduce the cross section area of implications for lateral control at high angles of the narrowest part of the venturi to half that of attack, and thus for recovery from unusual the opening, for example, the velocity must attitudes entered from stalls.
    [Show full text]
  • Infinite Swept-Wing Reynolds-Averaged Navier-Stokes Computations with Full En Transition Criterion
    27TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES INFINITE SWEPT-WING REYNOLDS-AVERAGED NAVIER-STOKES COMPUTATIONS WITH FULL EN TRANSITION CRITERION Zhang Kun, Song Wenping National Key Laboratory of Science and Technology on Aerodynamic Design and Research, School of Aeronautics, Northwestern Polytechnical University, Xi’an, China [email protected] [email protected] Keywords: full en method, transition, three- dimensional RANS, swept wing, linear stability Abstract accuracy of boundary layer heat conductivity n will reduce at least 25% [1]. A full e transition prediction method is coupled On the other hand, the method which can to the three-dimensional Reynolds-Averaged accurately predict the transition location is one Navier-Stokes (RANS) solver to predict the of the key technologies for designing the natural transition point automatically during the laminar flow wing. In order to improve the simulation of the flow around the infinite swept performance of the aircraft and to reduce air wings. The three-dimensional linear stability pollution during the cruise, the cruise drag equations are solved using the Cebeci- needs to be reduced. In general, for a typical Stewartson eigenvalue formulation. The swept-winged transport aircraft at cruise, the locations of the calculated transition points are frictional drag accounts for about 35% of the validated by the experimental results. With the total drag [2], so among the various drag reliable transition information, the accuracy of reduction technologies, the laminar flows drag the infinite swept wing’s aerodynamic reduction is one of the most promising performance calculated by the RANS solver has technologies. However, the design of natural been improved.
    [Show full text]
  • Demonstrating the Potential of Transitional CFD for Sailplane Design
    The Pennsylvania State University The Graduate School Department of Aerospace Engineering DEMONSTRATING THE POTENTIAL OF TRANSITIONAL CFD FOR SAILPLANE DESIGN A Thesis in Aerospace Engineering by Christopher J. Axten Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Science December 2019 ii The thesis of Christopher J. Axten was reviewed and approved* by the following: Mark Maughmer Professor of Aerospace Engineering Thesis Advisor Sven Schmitz Associate Professor of Aerospace Engineering Amy Pritchett Professor of Aerospace Engineering Head of the Department of Aerospace Engineering *Signatures are on file in the Graduate School iii ABSTRACT Traditional computational fluid dynamics solvers either model the flow as laminar or with assuming the presence of turbulence. If the flow is modeled with turbulence the initial influence of turbulence is minimal, and the flow can be considered laminar-like, but as the flow develops the amount of turbulence grows until it acts as a fully turbulent boundary layer. Neither approach properly models flow dynamics for the flight regime of a sailplane. To demonstrate the potential of using computational fluid dynamics for sailplane design a racing sailplane is analyzed with computational fluid dynamics using a recently developed transition model to accurately model viscous effects. The results of the analysis are validated against a conventional sailplane analysis program and are found to agree well. Regions with complex flows, such as the wing-fuselage juncture and the empennage juncture, are examined to highlight the potential for utilizing computational fluid dynamics to refine junctures in ways not possible with conventional design methods. Practical uses for computational fluid dynamics in sailplane analysis, such as investigating the stall characteristics and evaluating the tailwheel and pushrod fairing drags, are also discussed along with notable gains in aircraft performance.
    [Show full text]
  • TRANSITION in BOUNDARY LAYER FLOWS by IAIN D GARDINER
    TRANSITION IN BOUNDARY LAYER FLOWS by IAIN D GARDINER, BSc (Hons) Department of Mechanical and Industrial Engineering Thesis submitted to the C.N.A.A. for the Degree of Doctor of Philosophy March 1987 o 4' Dundee College of Technology Boundary Layer Wind Tunnel Facility DECLARATION I hereby declare that the following work has been composed by myself and that this dissertation has not been presented for any previous award of the C.N.A.A. or any other University. Iain D Gardiner 11 Transition in Boundary Layer Flows by Iain D Gardiner Abstract An experimental investigation of transition in boundary layer flows under the influence of various freestream conditions is described. Velocity profiles are obtained automatically by means of a stepper- motor driven traverse mechanism which carries a hot wire probe connected to a constant temperature anemometer and associated instrumentation. This was achieved by use of a data acquisition and control facility centred around a microcomputer with a Eurocard rack mounted extension. The automatic boundary layer traverse is software controlled and the data obtained is stored in a disc file for subsequent analysis and graphical display. As an integral part of this facility a successful method of obtaining reliable intermittency values from a hot wire signal was developed. The influence of freestream turbulence and pressure gradient upon transition within a boundary layer developing on a flat plate is elucidated by a series of controlled experiments. From the data accumulated, the concept of statistical similarity in transition regions is extended to include moderate non-zero pressure gradients, with the streamwise mean intermittency distribution described by the normal distribution function.
    [Show full text]
  • Aerodynamic Characteristics of Different Airfoils Under Varied
    applied sciences Article Aerodynamic Characteristics of Different Airfoils under Varied Turbulence Intensities at Low Reynolds Numbers Yang Zhang , Zhou Zhou *, Kelei Wang and Xu Li College of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China; [email protected] (Y.Z.); [email protected] (K.W.); [email protected] (X.L.) * Correspondence: [email protected] Received: 17 December 2019; Accepted: 14 February 2020; Published: 2 March 2020 Featured Application: The study of airfoil affected by the unsteady jet flow of turboelectric distributed propulsion aircraft was conducted, which lays a good foundation for future design of turboelectric distributed propulsion (TeDP) unmanned aerial vehicles (UAVs). Abstract: A numerical study was conducted on the influence of turbulence intensity and Reynolds number on the mean topology and transition characteristics of flow separation to provide better understanding of the unsteady jet flow of turboelectric distributed propulsion (TeDP) aircraft. By solving unsteady Reynolds averaged Navier-Stokes (URANS) equation based on C-type structural mesh and γ Ref transition model, the aerodynamic characteristics of the NACA0012 airfoil at − θt different turbulence intensities was calculated and compared with the experimental results, which verifies the reliability of the numerical method. Then, the effects of varied low Reynolds numbers and turbulence intensities on the aerodynamic performance of NACA0012 and SD7037 were investigated. The results show that higher turbulence intensity or Reynolds number leads to more stable airfoil aerodynamic performance, larger stalling angle, and earlier transition with a different mechanism. The generation and evolution of the laminar separation bubble (LSB) are closely related to Reynolds number, and it would change the effective shape of the airfoil, having a big influence on the airfoil’s aerodynamic characteristics.
    [Show full text]
  • Transition Prediction in Incompressible Boundary Layer with Finite-Amplitude Streaks
    energies Article Transition Prediction in Incompressible Boundary Layer with Finite-Amplitude Streaks Juan Ángel Martín 1,* and Pedro Paredes 2 1 E.T.S.I. Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Plaza Cardenal Cisneros 3, 28040 Madrid, Spain 2 National Institute of Aerospace (NIA), Hampton, VA 23681, USA; [email protected] * Correspondence: [email protected] Abstract: Modulating the boundary layer velocity profile is a very promising strategy for achieving transition delay and reducing the friction of the plate. By perturbing the flow with counter-rotating vortices that undergo transient, non-modal growth, streamwise-aligned streaks are formed inside the boundary layer, which have been proved (theoretical and experimentally) to be very robust flow structures. In this paper, we employ efficient numerical methods to perform a parametric stability investigation of the three-dimensional incompressible flat-plate boundary layer with finite- amplitude streaks. For this purpose, the Boundary Region Equations (BREs) are applied to solve the nonlinear downstream evolution of finite amplitude streaks. Regarding the stability analysis, the linear three-dimensional plane-marching Parabolized Stability Equations (PSEs) concept constitutes the best candidate for this task. Therefore, a thorough parametric study is presented, analyzing the instability characteristics with respect to critical conditions of the modified incompressible zero-pressure-gradient flat-plate boundary layer, by means of finite-amplitude linearly optimal and suboptimal disturbances or streaks. The parameter space is extended from low- to high- amplitude Citation: Martín, J.A.; Paredes, P. streaks, accurately documenting the transition delay for low-amplitude streaks and the amplitude Transition Prediction in threshold for streak shear layer instability or bypass transition, which drastically displaces the Incompressible Boundary Layer with transition front upstream.
    [Show full text]
  • Surface Structure and Its Effect on Reducing Drag
    Dissertations and Theses 5-2015 Surface Structure and Its Effect on Reducing Drag Ragini Ramachandran Follow this and additional works at: https://commons.erau.edu/edt Part of the Aerospace Engineering Commons Scholarly Commons Citation Ramachandran, Ragini, "Surface Structure and Its Effect on Reducing Drag" (2015). Dissertations and Theses. 277. https://commons.erau.edu/edt/277 This Thesis - Open Access is brought to you for free and open access by Scholarly Commons. It has been accepted for inclusion in Dissertations and Theses by an authorized administrator of Scholarly Commons. For more information, please contact [email protected]. SURFACE STRUCTURE AND ITS EFFECT ON REDUCING DRAG by Ragini Ramachandran A Thesis Submitted to the College of Engineering Department of Aerospace Engineering in Partial Fulfillment of the Requirements for the Degree of Master of Science in Aerospace Engineering Embry-Riddle Aeronautical University Daytona Beach, Florida May 2015 Acknowledgements I would like to express my gratitude towards my advisor, Dr. David J. Sypeck, for accepting this research as a thesis topic. He has constantly guided me and has been instrumental in designing and fabricating the miniature wind tunnel, which would not have been possible without him. He has spent many long hours carefully reviewing and purchasing the fabrication parts for the wind tunnel and test models, preparing the wing models with various fabrics, conducting the experiments, and also going through the thesis document. I would like to acknowledge prior support from the National Science Foundation (8802 Servohydraulic Test System), the ERAU College of Engineering (load cell), and the ERAU Department of Aerospace Engineering (small equipment and supplies).
    [Show full text]
  • Space-Time Structure Around the Transition Point in Channel Flow Revealed by the Stochastic Determinism
    Sixth International Symposium on Turbulence and Shear Flow Phenomena Seoul, Korea, 22-24 June 2009 SPACE-TIME STRUCTURE AROUND THE TRANSITION POINT IN CHANNEL FLOW REVEALED BY THE STOCHASTIC DETERMINISM Hiromu Shimiya Faculty of Science and Engineering, Waseda University 3-4-1 Ookubo, Shinjuku-ku, Tokyo, 169-8555, Japan Ken Naitoh Faculty of Science and Engineering, Waseda University 3-4-1 Ookubo, Shinjuku-ku, Tokyo, 169-8555, Japan [email protected] ABSTRACT difference method extended based on the multi-level The transition to turbulence (Reynolds, 1883) has formulation (Naitoh and Kuwahara, 1992), which can attracted people. Large eddy simulation (LES) and direct calculate spatial derivatives of physical quantities and numerical simulation (DNS) of the transition to turbulence integrated quantities accurately, leads us to the new stage of in straight channels employed the spatial cyclic boundary computational fluid dynamics. conditions between the inlet and outlet of the channel. Starting point of the methodology is related to the equation (Moin and Kim, 1982; Kawamura and Kuwahara, 1985) of the divergence of velocity in the multi-level formulation, Thus, these previous researches capture only the transition which is transformed from the compressible Navier-Stokes in time, although the spatial transition point where the equation while maintaining the Gibbs formula. (Naitoh and laminar flow changes to turbulence could not be computed. Kuwahara, 1992.) The formulation shows that the second Recently, some approaches tried to compute the transition derivative of velocity, the divergence of velocity, controls point in straight channel for the flows having large the physical quantities such as pressure, velocity, and disturbances at the inlet.
    [Show full text]