Arxiv:2104.13904V2 [Physics.Flu-Dyn] 25 Jun 2021

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:2104.13904V2 [Physics.Flu-Dyn] 25 Jun 2021 Are Superfluids Lifting? A Novel Variational Theory of Lift Cody Gonzalez∗ and Haithem E. Tahay University of California, Irvine The long-standing controversy about the ability of superfluids to generate lift has persisted mainly because of the lack of a lift theory for superfluids. Here, we revive the historical, yet less-often utilized, Hertz' principle of least curvature, exploiting it to develop a new variational analogue of Euler's equations for the dynamics of an ideal fluid. Using this new variational formulation, we generalize the century-old problem of the flow over a two-dimensional body, to find that lift is a direct consequence of curvature, thus reconciling the seemingly contradicting results in the literature of superfluids[1, 2]. The developed variational principle reduces to the classical Kutta-Zhukovsky condition in the special case of a sharp-edged airfoil, which challenges the accepted wisdom about the Kutta condition being a manifestation of viscous effects. Rather, we found that it represents conservation of momentum. Moreover, the developed variational principle provides, for the first time, a theoretical model for lift over smooth shapes without sharp edges where the Kutta condition is not applicable. This theory resolves the classical debate about the ability of a superfluid to generate lift as it provides conditions on the geometry of the body to be lifting in a dissipation-free medium. I. INTRODUCTION to support this hypothesis. He studied the flow of a su- perfluid Helium II over a flat plate and an ellipse and The classical problem of the flow over a two- found a vanishingly small lift force over these objects at dimensional body (e.g., flat plate, airfoil) is a century- non-zero angles of attack. Hence, he concluded that su- old problem that is ubiquitously found in almost ev- perfluids are non-lifting and that \the classical viscosity ery textbook in aerodynamics and fluid mechanics [3{8]. boundary condition at the trailing edge (the Kutta con- The problem is analytically solvable thanks to three ele- dition) does not apply". ments. First, the potential-flow around a circular cylin- Note that it is already known that the smaller the vis- der is readily known since the 1877 seminal paper of cosity (i.e., the larger the Reynolds number), the larger Lord Rayleigh [9]. Second, the Riemann mapping the- the lift. So, as viscosity is decreased, lift either increases orem (dating back to Riemann's PhD Dissertation 1851) or remains constant. That is, the lift value in the limit ensures that any simply connected domain can be (bi- to a vanishingly small viscosity is finite: the Kutta- holomorphically) mapped to the open disc. So, the flow Zhukovsky lift. Therefore, the above question about around any two-dimensional shape can be easily con- the non-lifting nature of superfluids is interesting from a structed from the cylinder flow via conformal mapping physics perspective because if the hypothesis were true, it between the cylinder and the shape of interest. How- would represent one example in Nature where the physics ever, this solution is not unique. One can always add is not continuous in the limit; the limit of lift as viscosity a circulation of arbitrary strength at the center of the approaches zero is different from the lift value at strictly cylinder, which does not affect the no-penetration bound- zero viscosity (i.e., inviscid). ary condition at all. Interestingly, this circulation is In contrast to the above hypothesis, Musser et al. [2] of paramount importance for lift production; in fact, recently performed quantum simulations of the Gross- it solely dictates the amount of lift generated. There- Pitaevskii equation governing the flow dynamics over an fore, the potential-flow theory alone cannot predict the airfoil in a superfluid. Their simulations show a quan- generated lift force; a closure condition must be pro- tized version of the Kutta-Zhukovsky lift despite the lack vided to fix the dynamically-correct amount of circula- of viscosity in their simulations. However, the continuum tion. Third, the Kutta-Zhukovsky condition has tradi- hypothesis may not be applicable in their ultra-small- tionally provided such a closure. scale simulations. arXiv:2104.13904v2 [physics.flu-dyn] 25 Jun 2021 The question of whether a superfluid (of zero viscos- From the above discussion, there seems to be contra- ity/friction) can generate lift when it flows over a nicely diction between Craig's experiment [1] and the computa- designed body has been lingering for decades. The ac- tions of Musser et al. [2]. The question is: How can one cepted norm asserts that superfluids are non-lifting, and obtain the flow field of a purely inviscid (e.g., a super- the experiment by Craig [1] at Caltech in the 1950s seems fliud) flow over a two-dimensional shape (e.g., flat plate or airfoil)? Ironically, insofar this question seems simple with a well known answer, it is actually not. First of all, why can we not simply perform computational sim- ∗ [email protected] ulations of Euler's equations of an ideal fluid? In fact, y [email protected] it is well known that 2-D inviscid Euler's computations 2 achieve a Kutta-Zhukovsky lift without viscosity [10, 11]. field is determined from purely kinematic analysis with- However, Euler's computations may not be decisive; the out any consideration for dynamical aspects. Therefore, artificial viscosity or numerical dissipation needed for sta- it is fair to expect that such a pure kinematic analysis is ble computations of Euler's equations may imply that the not sufficient to uniquely determine the flow field; the fix obtained computational results are not for a purely in- must come from a dynamical consideration. viscid fluid, but rather for a viscous flow in the limit of a vanishingly small viscosity (see Hirsch [12], Sec. 19.4). Second, the well known potential-flow solution is not B. Variational Formulation is the Solution complete; it requires an auxiliary condition, which is typi- cally provided by the Kutta condition as discussed above. Based on the above discussion, a proper closure condi- But, the accepted wisdom implies that the Kutta condi- tion (i.e., a condition that provides circulation dynamics) tion is a manifestation of viscous effects; it implicitly ac- in potential flow must come from dynamical considera- counts for viscous effects in a potential-flow framework. tions. The challenge is: Can we project Euler's dynam- So, if this common belief is accepted, the Kutta condition ical equations on a one-dimensional manifold to extract should not be used as a closure condition when studying the dynamics of circulation alone? Dynamical equations a purely inviscid fluid. The question then is: What are of motion can be determined either from a Newtonian the closure conditions for the flow of a superfluid? Un- mechanics perspective or an analytical mechanics one. fortunately, there are none! The former stipulates isolating fluid particles and writ- In this paper, we develop a new variational theory of ing the equations of motion for each individual parti- lift that generalizes the century-old classical theory of lift cle even if the free variables in the system are signif- by Kutta and Zhukovksy; the new theory dispenses with icantly fewer than the total degrees of freedom of all the Kutta condition. We revive Hertz' principle of least individual particles due to kinematic or geometric con- curvature from the history of analytical mechanics to de- straints. However, the analytical (Lagrangian or vari- velop, for the first time, a closure condition for the flow ational) mechanics approach allows accepting the kine- of a superfluid over a two-dimensional body. In contrast matical constraints, ignoring the unknown forces that to the Kutta condition, the developed closure condition maintain them, and hence focusing on the relevant equa- is derived from first principles. The new variational con- tions of motion; it provides directly the relevant equa- dition reduces to the Kutta condition in the special case tions of motion for the free variables. of a sharp-edged airfoil, which challenges the accepted Projecting this discussion on the potential-flow case, wisdom about the Kutta condition being a manifestation one finds that the kinematical constraints of potential of viscous effects. Rather, it is found that the Kutta con- flow allows one to construct the entire flow field from the dition is a momentum conservation mechanism for the circulation free variable only. That is, while there are inviscid flow over a sharp-edged airfoil. Morever, the de- infinite degrees of freedom for the infinite fluid particles, veloped theory, unlike the classical theory, allows treat- there is only one free variable (the circulation) which, via ment of smooth shapes not necessarily with sharp trail- the potential-flow kinematical constraints, can be used ing edges. Using this new theory, we resolve the seeming to recover the motion of these infinite degrees of free- contradiction between various results in literature about dom [13]. Hence, the analytical/variational mechanics the lifting capability of superfluids. We provide condi- appears to be specially well-suited for this problem; it tions on the geometry of the body to generate lift in a will provide a single equation for the unknown circula- superfluid. tion without paying attention to the irrelevant degrees of freedom of the fluid particles or the unknown forces that maintain kinematical constraints. Simply, the first vari- II. ILL-POSEDNESS OF POTENTIAL-FLOW ation of the \objective function" with respect to circula- PROBLEMS tion must vanish|and this necessary condition provides a single dynamical equation in the unknown circulation. A. The Potential-Flow Theory Lacks Dynamical Based on the above discussion, two important conclu- Features sions are drawn: (i) A true closure/auxiliary condition for potential flow must come from dynamical considerations; To solve for the flow field of an incompressible and (ii) Variational principles would be particularly use- fluid, both the continuity (kinematics) and momentum ful to derive such dynamics.
Recommended publications
  • Clearing Certain Misconception in the Common Explanations of the Aerodynamic Lift
    Clearing certain misconception in the common explanations of the aerodynamic lift Clearing certain misconception in the common explanations of the aerodynamic lift Navinder Singh,∗ K. Sasikumar Raja, P. Janardhan Physical Research Laboratory, Ahmedabad, India. [email protected] October 30, 2018 Abstract Air travel has become one of the most common means of transportation. The most common question which is generally asked is: How does an airplane gain lift? And the most common answer is via the Bernoulli principle. It turns out that it is wrongly applied in common explanations, and there are certain misconceptions. In an alternative explanation the push of air from below the wing is argued to be the lift generating force via Newton’s law. There are problems with this explanation too. In this paper we try to clear these misconceptions, and the correct explanation, using the Lancaster-Prandtl circulation theory, is discussed. We argue that even the Lancaster-Prandtl theory at the zero angle of attack needs further insights. To this end, we put forward a theory which is applicable at zero angle of attack. A new length scale perpendicular to the lower surface of the wing is introduced and it turns out that the ratio of this length scale to the cord length of a wing is roughly 0.4930 ± 0.09498 for typical NACA airfoils that we analyzed. This invariance points to something fundamental. The idea of our theory is simple. The "squeezing" effect of the flow above the wing due to camber leads to an effective Venturi tube formation and leads to higher velocity over the upper surface of the wing and thereby reducing pressure according the Bernoulli theorem and generating lift.
    [Show full text]
  • How Do Wings Generate Lift? 2
    GENERAL ARTICLE How Do Wings Generate Lift? 2. Myths, Approximate Theories and Why They All Work M D Deshpande and M Sivapragasam A cambered surface that is moving forward in a fluid gen- erates lift. To explain this interesting fact in terms of sim- pler models, some preparatory concepts were discussed in the first part of this article. We also agreed on what is an accept- able explanation. Then some popular models were discussed. Some quantitative theories will be discussed in this conclud- ing part. Finally we will tie up all these ideas together and connect them to the rigorous momentum theorem. M D Deshpande is a Professor in the Faculty of Engineering The triumphant vindication of bold theories – are these not the and Technology at M S Ramaiah University of pride and justification of our life’s work? Applied Sciences. Sherlock Holmes, The Valley of Fear Sir Arthur Conan Doyle 1. Introduction M Sivapragasam is an We start here with two quantitative theories before going to the Assistant Professor in the momentum conservation principle in the next section. Faculty of Engineering and Technology at M S Ramaiah University of Applied 1.1 The Thin Airfoil Theory Sciences. This elegant approximate theory takes us further than what was discussed in the last two sections of Part 111, and quantifies the Resonance, Vol.22, No.1, ideas to get expressions for lift and moment that are remarkably pp.61–77, 2017. accurate. The pressure distribution on the airfoil, apart from cre- ating a lift force, leads to a nose-up or nose-down moment also.
    [Show full text]
  • Lecture 4: Panel Methods G. Dimitriadis Introduction
    Aerodynamics Lecture 4: Panel methods G. Dimitriadis Introduction • Until now we’ve seen two methods for modelling wing sections in ideal flow: –Conformal mapping: Can model only a few classes of wing sections –Thin airfoil theory: Can model any wing section but ignores the thickness. • Both methods are not general. • A more general approach will be presented here. 2D Panel methods • 2D Panel methods refers to numerical methods for calculating the flow around any wing section. • They are based on the replacement of the wing section’s geometry by singularity panels, such as source panels, doublet panels and vortex panels. • The usual boundary conditions are imposed: – Impermeability – Kutta condition Panel placement • Eight panels of length Sj • For each panel, sj=0-Sj. • Eight control (collocation) points (xcj,ycj) located in the middle of each panel. • Nine boundary points (xbj,ybj) • Normal and tangential unit vector on each panel, nj, tj • Vorticity (or source strength) on each panel j(s) (or j(s)) Problem statement • Use linear panels • Use constant singularity strength on each panel, j(s)=j (j(s)=j). • Add free stream U at angle . • Apply boundary conditions: – Far field: automatically satisfied if using source or vortex panels – Impermeability: Choose Neumann or Dirichlet. • Apply Kutta condition. • Find vortex and/or source strength distribution that will satisfy Boundary Conditions and Kutta condition. Panel choice • It is best to chose small panels near the leading and trailing edge and large panels in the middle: x 1 = ()cos +1 , = 0 2 c 2 Notice that now x/c 17 29 begins at 1, passes 16 30 through 0 and then goes 2 1 back to 1.
    [Show full text]
  • Explicit Kutta Condition Correction for Rotary Wing Flows A. D'alascio, A
    Transactions on Modelling and Simulation vol 18, © 1997 WIT Press, www.witpress.com, ISSN 1743-355X Explicit Kutta Condition Correction for Rotary Wing Flows A. D'Alascio, A. Visingardi, and P. Renzoni CIRA, Centra Italiano Ricerche Aerospaziali EMail: [email protected], [email protected], p.renzoni@cira. it Abstract An explicit correction to the Kutta condition applied to rotary wing potential flows is described. Indeed, the simple application of the classical steady Kutta condition, which ensures the uniqueness of the solution of the potential flow problems [1,2], is unable to guarantee a zero pressure jump at the blade trailing edge. Several investigations performed on this problem have indicated the presence of considerable non-zero pressure jumps for highly 3D flows (blade root and tip) [3] and for unsteady flows characterized by high reduced frequencies [4]. The application of a Kutta condition correction in CIRA's BEM code RAMSYS [3] has shown the ability to ensure a continuous pressure across the blade trailing edge. 1 Introduction The flow field around helicopter rotors is unsteady and highly three-dimensio- nal. A low frequency unsteadiness is induced in the flow by the relative motion between rotating and non-rotating components of the helicopter. The interaction between each main rotor blade with the wake of preceding blades is responsible for mid frequency unsteadiness whereas high frequency unsteadiness is generated by blade/tip vortex encounters. The rotor blades during their motion are subject to cross flows which make the flow field three-dimensional, expecially in the vicinity of the ro- Transactions on Modelling and Simulation vol 18, © 1997 WIT Press, www.witpress.com, ISSN 1743-355X 780 Boundary Elements tor blade root and tip.
    [Show full text]
  • Explicit Role of Viscosity in Generating Lift
    AIAA JOURNAL Vol. 55, No. 11, November 2017 Technical Notes Explicit Role of Viscosity I. Introduction in Generating Lift HIS Note tries to elucidate the subtle but critical role of the fluid T viscosity in generating the aerodynamic lift. The viscous origin of the lift has not been widely understood or sufficiently highlighted Tianshu Liu∗ in the literature. In classical aerodynamics textbooks [1–4], the lift Western Michigan University, Kalamazoo, Michigan 49008 theory of an airfoil is usually developed based on the classical potential-flow theory, in which the Kutta–Joukowski (K-J) theorem and † ‡ is used as a key element in calculating the lift. However, the airfoil Shizhao Wang and Guowei He circulation in the K-J theorem cannot be automatically determined in Chinese Academy of Sciences, 100190 Beijing, the potential-flow framework alone because the vorticity cannot be People’s Republic of China physically generated in an inviscid flow. This fundamental difficulty is intrinsically related to D’Alembert’s paradox stating that the DOI: 10.2514/1.J055907 integrated pressure force of a body is zero in a steady inviscid irrotational incompressible flow [5,6]. To circumvent this problem, the Kutta condition is imposed at the airfoil trailing edge to determine Nomenclature the circulation. To a great extent, the success of the classical Cd = drag coefficient circulation theory in predicting the lift depends on the clever Cl = lift coefficient application of the Kutta condition and its generalized forms as the Cp = pressure coefficient phenomenological models of the viscous-flow effect on the lift c = chord length, m generation [7].
    [Show full text]
  • Joukowski Theorem for Multi-Vortex and Multi-Airfoil Flow (A
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Elsevier - Publisher Connector Chinese Journal of Aeronautics, (2014),27(1): 34–39 Chinese Society of Aeronautics and Astronautics & Beihang University Chinese Journal of Aeronautics [email protected] www.sciencedirect.com Generalized Kutta–Joukowski theorem for multi-vortex and multi-airfoil flow (a lumped vortex model) Bai Chenyuan, Wu Ziniu * School of Aerospace, Tsinghua University, Beijing 100084, China Received 5 January 2013; revised 20 February 2013; accepted 25 February 2013 Available online 31 July 2013 KEYWORDS Abstract For purpose of easy identification of the role of free vortices on the lift and drag and for Incompressible flow; purpose of fast or engineering evaluation of forces for each individual body, we will extend in this Induced drag; paper the Kutta–Joukowski (KJ) theorem to the case of inviscid flow with multiple free vortices and Induced lift; multiple airfoils. The major simplification used in this paper is that each airfoil is represented by a Multi-airfoils; lumped vortex, which may hold true when the distances between vortices and bodies are large Vortex enough. It is found that the Kutta–Joukowski theorem still holds provided that the local freestream velocity and the circulation of the bound vortex are modified by the induced velocity due to the out- side vortices and airfoils. We will demonstrate how to use the present result to identify the role of vortices on the forces according to their position, strength and rotation direction. Moreover, we will apply the present results to a two-cylinder example of Crowdy and the Wagner example to demon- strate how to perform fast force approximation for multi-body and multi-vortex problems.
    [Show full text]
  • Kutta Condition
    SNS COLLEGE OF TECHNOLOGY (An Autonomous Institution) COIMBATORE–35 DEPARTMENT OF AERONAUTICAL ENGINEERING KUTTA CONDITION: The Kutta condition is a principle in steady-flow fluid dynamics, especially aerodynamics, that is applicable to solid bodies with sharp corners, such as the trailing edges of airfoils. A body with a sharp trailing edge which is moving through a fluid will create about itself a circulation of sufficient strength to hold the rear stagnation point at the trailing edge. In fluid flow around a body with a sharp corner, the Kutta condition refers to the flow pattern in which fluid approaches the corner from both directions, meets at the corner, and then flows away from the body. None of the fluid flows around the sharp corner. Prepared by Ms.X.Bernadette Evangeline AP/AERO 16AE208 AERODYNAMICS-I SNS COLLEGE OF TECHNOLOGY (An Autonomous Institution) COIMBATORE–35 DEPARTMENT OF AERONAUTICAL ENGINEERING The Kutta condition is significant when using the Kutta–Joukowski theorem to calculate the lift created by an airfoil with a sharp trailing edge. The value of circulation of the flow around the airfoil must be that value which would cause the Kutta condition to exist. In 2-D potential flow, if an airfoil with a sharp trailing edge begins to move with an angle of attack through air, the two stagnation points are initially located on the underside near the leading edge and on the topside near the trailing edge, just as with the cylinder. As the air passing the underside of the airfoil reaches the trailing edge it must flow around the trailing edge and along the topside of the airfoil toward the stagnation point on the topside of the airfoil.
    [Show full text]
  • A Kutta Condition Enforcing BEM Technology for Airfoil Aerodynamics
    Transactions on Modelling and Simulation vol 15, © 1997 WIT Press, www.witpress.com, ISSN 1743-355X A Kutta Condition Enforcing BEM Technology for Airfoil Aerodynamics G. lannelli*, C. Grille*, L. Tulumello* Department of Mechanical and Aerospace Engineering and Eng. Science, University of Tennessee, Knoxville, USA * Department of Mechanical and Aeronautical Engineering, University of Palermo, Italy Abstract An analytical perturbation stream function procedure is established for po- tential flows about lifting airfoils, which exactly enforces the fundamental Kutta condition. Using Green's identity, this formulation is transformed into an equivalent boundary integral equation. A boundary finite element discretization of this equation then yields a linear system for the efficient determination of the airfoil surface tangential velocity and associated lift 1 Introduction For many years boundary integral equation procedures, or panel methods, have been widely used for determining the potential flow field about com- plex aerodynamics configurations*. In the potential function procedure, the potential values are related to the boundary potentials and associated normal derivatives, whereas in the singularity procedure the local potential function value is expressed in terms of boundary sources, doublets and vor- tices of appropriate intensity. In either case, the reliability of the computed results critically depends on accurate boundary distributions of singularities or potentials. Furthermore, the Kutta condition, requiring the aft stagna- tion point to coincide with the trailing edge, must be adequately enforced. Several implementations of the potential function procedure indirectly en- force the Kutta condition. Conversely, the singularity method employs a trailing edge discretization node where the vanishing velocity condition can be enforced. Nevertheless, it has been reported that the substitution of the trailing edge integral equation with the vanishing velocity condition may lead to solution non-uniqueness.
    [Show full text]
  • AA200 Ch 11 Two-Dimensional Airfoil Theory Cantwell.Pdf
    CHAPTER 11 TWO-DIMENSIONAL AIRFOIL THEORY ___________________________________________________________________________________________________________ 11.1 THE CREATION OF CIRCULATION OVER AN AIRFOIL In Chapter 10 we worked out the force that acts on a solid body moving in an inviscid fluid. In two dimensions the force is F d = Φndlˆ − U (t) × Γ(t)kˆ (11.1) ρ oneunitofspan dt ∫ ∞ ( ) Aw where Γ(t) = !∫ UicˆdC is the circulation about the body integrated counter-clockwise Cw about a path that encloses the body. It should be stated at the outset that the creation of circulation about a body is fundamentally a viscous process. The figure below from Van Dyke’s Album of Fluid Motion depicts what happens when a viscous fluid is moved impulsively past a wedge. Figure 11.1 Starting vortex formation about a corner in an impulsively started incompressible fluid. 1 A purely potential flow solution would negotiate the corner giving rise to an infinite velocity at the corner. In a viscous fluid the no slip condition prevails at the wall and viscous dissipation of kinetic energy prevents the singularity from occurring. Instead two layers of fluid with oppositely signed vorticity separate from the two faces of the corner to form a starting vortex. The vorticity from the upstream facing side is considerably stronger than that from the downstream side and this determines the net sense of rotation of the rolled up vortex sheet. During the vortex formation process, there is a large difference in flow velocity across the sheet as indicated by the arrows in Figure 11.1. Vortex formation is complete when the velocity difference becomes small and the starting vortex drifts away carried by the free stream.
    [Show full text]
  • On the Circulation and the Position of the Forward Stagnation Point on Airfoils
    On the circulation and the position of the forward stagnation point on airfoils J. Meseguer (corresponding author), G. Alonso, A. Sanz-Andres and I. Perez-Grande IDR/UPM. EJ.S.I. Aeronauticos, Universidad Politecnica de Madrid, E-28040 Madrid, Spain E-mail: jmeseguer@idr. upm. es Abstract Lift and velocity circulation around airfoils are two aspects of the same phenomenon when airfoils are not stalled and the Kutta-Joukowski theorem applies. This theorem establishes a linear dependence between lift and circulation, which breaks when stalling occurs. As the angle of attack increases beyond this point, the circulation vanishes. Since the circulation determines to a great extent the position of the forward stagnation point on an airfoil, the measurement of this position is an easy and simple way to determine the circulation, which is of help in understanding the role of the latter in the generation of aerodynamic forces on airfoils. Keywords Kutta-Joukowski theorem; airfoil; stagnation point Introduction Wing profiles are aerodynamic devices which are designed to provide high values of the ratio lift/aerodynamic drag, lid. Airfoils are characterised by the camber line, the thickness distribution and the angle of attack. These streamlined two- dimensional bodies generally have a rounded leading edge and a wedged trailing edge, the latter being of paramount importance in generating lift. As is well known, due to the airfoil's shape, in normal operation the flow is accel­ erated on the airfoil's upper side and decelerated on its lower side. Hence, the pres­ sure decreases on the upper side and increases on the lower side, and thus a net force appears.
    [Show full text]
  • Chapter 1 Introduction to Fluid Flow
    CHAPTER 1 INTRODUCTION TO FLUID FLOW 1.1 INTRODUCTION Fluid flows play a crucial role in a vast variety of natural phenomena and man- made systems. The life-cycles of stars, the creation of atmospheres, the sounds we hear, the vehicles we ride, the systems we build for flight, energy generation and propulsion all depend in an important way on the mechanics and thermody- namics of fluid flow. The purpose of this course is to introduce students in Aeronautics and Astronautics to the fundamental principles of fluid mechanics with emphasis on the development of the equations of motion as well as some of the analytical tools from calculus needed to solve practically important problems involving flows in channels along walls and over lifting bodies. 1.2 CONSERVATION OF MASS Mass is neither created nor destroyed. This basic principle of classical physics is one of the fundamental laws governing fluid motion and is a good departure point for our introductory discussion. Figure 1.1 below shows an infinitesimally small stationary, rectangular control volume 6x6y6z through which a fluid is assumed to be moving. A control vol- ume of this type with its surface fixed in space is called an Eulerian control volume. The fluid velocity vector has components UUVW= (),, in the xxyz= (),, directions and the fluid density is l . In a general, unsteady, com- pressible flow, all four flow variables may depend on position and time. The law of conservation of mass over this control volume is stated as bjc 1.1 3/26/13 Conservation of mass ¨¬Rate of mass ¨¬Rate of mass ¨¬Rate of mass «««««« ««accumulation ««flow «« flow ©­= ©­– ©­(1.1) ««inside the control ««into the control ««out of the control «««««« ª®volume ª®volume ª®volume z lW z + 6z ()x +y6x,,+z6y + 6z W(x,y,z,t) lV y + 6y V(x,y,z,t) U U(x,y,z,t) lU x + 6x l()xyzt,,, lU x 6z y ()xyz,, lV lW y 6y z 6x x Figure 1.1 Fixed control volume in a moving fluid.
    [Show full text]
  • Unsteady Kutta Condition and Vortex-Sheet Generation
    XXIV ICTAM, 21-26 August 2016, Montreal, Canada UNSTEADY KUTTA CONDITION AND VORTEX-SHEET GENERATION Xi Xia1 and Kamran Mohseni ∗1;2 1Department of Mechanical and Aerospace Engineering, University of Florida, Gainesville, FL, USA 2Department of Electrical and Computer Engineering, University of Florida, Gainesville, FL, USA Summary Vortex method has been widely used for low-Reynolds-number flights as it reduces the computation domain from the entire flow field to only finite vortical structures. An essential challenge of the vortex method lies in the prediction of the vortex sheets separated from the solid body. To tackle this problem, this study extends the classical Kutta condition to unsteady situations based on the physical sense that flow cannot turn around a sharp edge. This unsteady Kutta condition can be readily applied to an airfoil with cusped trailing edge, where the forming vortex sheet is known to be tangential to the trailing edge. For a finite-angle trailing edge, previous studies indicate that the direction of the forming vortex sheet is ambiguous. Therefore, this study proposes a novel analytical formulation to determine the angle of the trailing-edge vortex sheet based on momentum conservation in the direction normal to the forming vortex sheet. INTRODUCTION Natural flies are observed to be superior than man-made aerial vehicles in terms of aerodynamic performance, especially the high lift-generation mechanisms. Early experimental investigations suggested that the existence of an attached leading edge vortex (LEV) contributes to enhanced lift generation for flapping wings [2]. To understand the fundamental mechanism of the LEV, recent analytical studies have been extensively focused on using vortex method (diagram shown in Fig.
    [Show full text]