Analytical Vortex Solutions to the Navier-Stokes Equation
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Energy, Vorticity and Enstrophy Conserving Mimetic Spectral Method for the Euler Equation
Master of Science Thesis Energy, vorticity and enstrophy conserving mimetic spectral method for the Euler equation D.J.D. de Ruijter, BSc 5 September 2013 Faculty of Aerospace Engineering · Delft University of Technology Energy, vorticity and enstrophy conserving mimetic spectral method for the Euler equation Master of Science Thesis For obtaining the degree of Master of Science in Aerospace Engineering at Delft University of Technology D.J.D. de Ruijter, BSc 5 September 2013 Faculty of Aerospace Engineering · Delft University of Technology Copyright ⃝c D.J.D. de Ruijter, BSc All rights reserved. Delft University Of Technology Department Of Aerodynamics, Wind Energy, Flight Performance & Propulsion The undersigned hereby certify that they have read and recommend to the Faculty of Aerospace Engineering for acceptance a thesis entitled \Energy, vorticity and enstro- phy conserving mimetic spectral method for the Euler equation" by D.J.D. de Ruijter, BSc in partial fulfillment of the requirements for the degree of Master of Science. Dated: 5 September 2013 Head of department: prof. dr. F. Scarano Supervisor: dr. ir. M.I. Gerritsma Reader: dr. ir. A.H. van Zuijlen Reader: P.J. Pinto Rebelo, MSc Summary The behaviour of an inviscid, constant density fluid on which no body forces act, may be modelled by the two-dimensional incompressible Euler equations, a non-linear system of partial differential equations. If a fluid whose behaviour is described by these equations, is confined to a space where no fluid flows in or out, the kinetic energy, vorticity integral and enstrophy integral within that space remain constant in time. Solving the Euler equations accompanied by appropriate boundary and initial conditions may be done analytically, but more often than not, no analytical solution is available. -
Kármán Vortex Street Energy Harvester for Picoscale Applications
Kármán Vortex Street Energy Harvester for Picoscale Applications 22 March 2018 Team Members: James Doty Christopher Mayforth Nicholas Pratt Advisor: Professor Brian Savilonis A Major Qualifying Project submitted to the Faculty of WORCESTER POLYTECHNIC INSTITUTE in partial fulfilment of the requirements for the degree of Bachelor of Science This report represents work of WPI undergraduate students submitted to the faculty as evidence of a degree requirement. WPI routinely publishes these reports on its web site without editorial or peer review. For more information about the projects program at WPI, see http://www.wpi.edu/Academics/Projects. Cover Picture Credit: [1] Abstract The Kármán Vortex Street, a phenomenon produced by fluid flow over a bluff body, has the potential to serve as a low-impact, economically viable alternative power source for remote water-based electrical applications. This project focused on creating a self-contained device utilizing thin-film piezoelectric transducers to generate hydropower on a pico-scale level. A system capable of generating specific-frequency vortex streets at certain water velocities was developed with SOLIDWORKS modelling and Flow Simulation software. The final prototype nozzle’s velocity profile was verified through testing to produce a velocity increase from the free stream velocity. Piezoelectric testing resulted in a wide range of measured dominant frequencies, with corresponding average power outputs of up to 100 nanowatts. The output frequencies were inconsistent with predicted values, likely due to an unreliable testing environment and the complexity of the underlying theory. A more stable testing environment, better verification of the nozzle velocity profile, and fine-tuning the piezoelectric circuit would allow for a higher, more consistent power output. -
Claus Kahlert and Otto E. Rössler Institute for Physical and Theoretical Chemistry, University of Tübingen Z. Naturforsch
Chaos as a Limit in a Boundary Value Problem Claus Kahlert and Otto E. Rössler Institute for Physical and Theoretical Chemistry, University of Tübingen Z. Naturforsch. 39a, 1200- 1203 (1984); received November 8, 1984 A piecewise-linear. 3-variable autonomous O.D.E. of C° type, known to describe constant- shape travelling waves in one-dimensional reaction-diffusion media of Rinzel-Keller type, is numerically shown to possess a chaotic attractor in state space. An analytical method proving the possibility of chaos is outlined and a set of parameters yielding Shil'nikov chaos indicated. A symbolic dynamics technique can be used to show how the limiting chaos dominates the behavior even of the finite boundary value problem. Reaction-diffusion equations occur in many dis boundary conditions are assumed. This result was ciplines [1, 2], Piecewise-linear systems are especially obtained by an analytical matching method. The amenable to analysis. A variant to the Rinzel-Keller connection to chaos theory (Smale [6] basic sets) equation of nerve conduction [3] can be written as was not evident at the time. In the following, even the possibility of manifest chaos will be demon 9 ö2 — u u + p [- u + r - ß + e (ii - Ö)], strated. o t oa~ In Fig. 1, a chaotic attractor of (2) is presented. Ö The flow can be classified as an example of "screw- — r = - e u + r , (1) 0/ type-chaos" (cf. [7]). A second example of a chaotic attractor is shown in Fig. 2. A 1-D projection of a where 0(a) =1 if a > 0 and zero otherwise; d is the 2-D cross section through the chaotic flow in a threshold parameter. -
Lecture 18 Ocean General Circulation Modeling
Lecture 18 Ocean General Circulation Modeling 9.1 The equations of motion: Navier-Stokes The governing equations for a real fluid are the Navier-Stokes equations (con servation of linear momentum and mass mass) along with conservation of salt, conservation of heat (the first law of thermodynamics) and an equation of state. However, these equations support fast acoustic modes and involve nonlinearities in many terms that makes solving them both difficult and ex pensive and particularly ill suited for long time scale calculations. Instead we make a series of approximations to simplify the Navier-Stokes equations to yield the “primitive equations” which are the basis of most general circu lations models. In a rotating frame of reference and in the absence of sources and sinks of mass or salt the Navier-Stokes equations are @ �~v + �~v~v + 2�~ �~v + g�kˆ + p = ~ρ (9.1) t r · ^ r r · @ � + �~v = 0 (9.2) t r · @ �S + �S~v = 0 (9.3) t r · 1 @t �ζ + �ζ~v = ω (9.4) r · cpS r · F � = �(ζ; S; p) (9.5) Where � is the fluid density, ~v is the velocity, p is the pressure, S is the salinity and ζ is the potential temperature which add up to seven dependent variables. 115 12.950 Atmospheric and Oceanic Modeling, Spring '04 116 The constants are �~ the rotation vector of the sphere, g the gravitational acceleration and cp the specific heat capacity at constant pressure. ~ρ is the stress tensor and ω are non-advective heat fluxes (such as heat exchange across the sea-surface).F 9.2 Acoustic modes Notice that there is no prognostic equation for pressure, p, but there are two equations for density, �; one prognostic and one diagnostic. -
The Vorticity Equation in a Rotating Stratified Fluid
Chapter 7 The Vorticity Equation in a Rotating Stratified Fluid The vorticity equation for a rotating, stratified, viscous fluid » The vorticity equation in one form or another and its interpretation provide a key to understanding a wide range of atmospheric and oceanic flows. » The full Navier-Stokes' equation in a rotating frame is Du 1 +∧fu =−∇pg − k +ν∇2 u Dt ρ T where p is the total pressure and f = fk. » We allow for a spatial variation of f for applications to flow on a beta plane. Du 1 +∧fu =−∇pg − k +ν∇2 u Dt ρ T 1 2 Now uu=(u⋅∇ ∇2 ) +ω ∧ u ∂u 1 2 1 2 +∇2 ufuku +bω +g ∧=- ∇pgT − +ν∇ ∂ρt di take the curl D 1 2 afafafωω+=ffufu +⋅∇−+∇⋅+∇ρ∧∇+ ωpT ν∇ ω Dt ρ2 Dω or =−uf ⋅∇ +.... Dt Note that ∧ [ ω + f] ∧ u] = u ⋅ (ω + f) + (ω + f) ⋅ u - (ω + f) ⋅ u, and ⋅ [ω + f] ≡ 0. Terminology ωa = ω + f is called the absolute vorticity - it is the vorticity derived in an a inertial frame ω is called the relative vorticity, and f is called the planetary-, or background vorticity Recall that solid body rotation corresponds with a vorticity 2Ω. Interpretation D 1 2 afafafωω+=ffufu +⋅∇−+∇⋅+∇ρ∧∇+ ωpT ν∇ ω Dt ρ2 Dω is the rate-of-change of the relative vorticity Dt −⋅∇uf: If f varies spatially (i.e., with latitude) there will be a change in ω as fluid parcels are advected to regions of different f. Note that it is really ω + f whose total rate-of-change is determined. -
A Concept of the Vortex Lift of Sharp-Edge Delta Wings Based on a Leading-Edge-Suction Analogy Tech Library Kafb, Nm
I A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY TECH LIBRARY KAFB, NM OL3042b NASA TN D-3767 A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY By Edward C. Polhamus Langley Research Center Langley Station, Hampton, Va. NATIONAL AERONAUTICS AND SPACE ADMINISTRATION For sale by the Clearinghouse for Federal Scientific and Technical Information Springfield, Virginia 22151 - Price $1.00 A CONCEPT OF THE VORTEX LIFT OF SHARP-EDGE DELTA WINGS BASED ON A LEADING-EDGE-SUCTION ANALOGY By Edward C. Polhamus Langley Research Center SUMMARY A concept for the calculation of the vortex lift of sharp-edge delta wings is pre sented and compared with experimental data. The concept is based on an analogy between the vortex lift and the leading-edge suction associated with the potential flow about the leading edge. This concept, when combined with potential-flow theory modified to include the nonlinearities associated with the exact boundary condition and the loss of the lift component of the leading-edge suction, provides excellent prediction of the total lift for a wide range of delta wings up to angles of attack of 20° or greater. INTRODUCTION The aerodynamic characteristics of thin sharp-edge delta wings are of interest for supersonic aircraft and have been the subject of theoretical and experimental studies for many years in both the subsonic and supersonic speed ranges. Of particular interest at subsonic speeds has been the formation and influence of the leading-edge separation vor tex that occurs on wings having sharp, highly swept leading edges. -
On the Effect of Pipe Boundary Layer Growth on the Formation of a Laminar Vortex Ring Generated by a Piston/Cylinder Arrangement
Theoret. Comput. Fluid Dynamics (2002) 15: 303–316 Theoretical and Computational Fluid Dynamics Springer-Verlag 2002 On the Effect of Pipe Boundary Layer Growth on the Formation of a Laminar Vortex Ring Generated by a Piston/Cylinder Arrangement Michael Shusser and Morteza Gharib Graduate Aeronautical Laboratories, California Institute of Technology, Pasadena, CA 91125, U.S.A. Moshe Rosenfeld Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel Kamran Mohseni Department of Aerospace Engineering Sciences, University of Colorado, Boulder, CO 80309-0429, U.S.A. Communicated by R.D. Moser Received 19 January 2000 and accepted 17 August 2001 Abstract. The growth of a boundary layer at the nozzle wall during laminar vortex ring formation by a nozzle flow generator (piston/cylinder arrangement) is analysed theoretically and numerically and used for modelling the formation of real vortex rings. The predictions of the model are in good agreement with previous experimental and numerical results. 1. Introduction Vortex rings are usually generated in the laboratory by the motion of a piston pushing a column of fluid of length L through an orifice or nozzle of diameter D. This results in a separation of the boundary layer at the edge of the orifice or nozzle and its subsequent spiral roll-up. The piston/cylinder arrangement has been extensively used to address the problem of vortex ring forma- tion (Shariff and Leonard, 1992; Lim and Nickels, 1995). Recently Gharib et al. (1998) in their experimental study of vortex ring formation addressed the question of the largest circulation that a vortex ring can attain, by increasing L/D while keeping the average piston velocity fixed. -
Three Types of Horizontal Vortices Observed in Wildland Mass And
1624 JOURNAL OF CLIMATE AND APPLIED METEOROLOGY VOLUME26 Three Types of Horizontal Vortices Observed in Wildland Mas~ and Crown Fires DoNALD A. HAINES U.S. Department ofAgriculture, Forest Service, North Central Forest Experiment Station, East Lansing, Ml 48823 MAHLON C. SMITH Department ofMechanical Engineering, Michigan State University, East Lansing, Ml 48824 (Manuscript received 25 October 1986, in final form 4 May 1987) ABSTRACT Observation shows that three types of horizontal vortices may form during intense wildland fires. Two of these vortices are longitudinal relative to the ambient wind and the third is transverse. One of the longitudinal types, a vortex pair, occurs with extreme heat and low to moderate wind speeds. It may be a somewhat common structure on the flanks of intense crown fires when burning is concentrated along the fire's perimeter. The second longitudinal type, a single vortex, occurs with high winds and can dominate the entire fire. The third type, the transverse vortex, occurs on the upstream side of the convection column during intense burning and relatively low winds. These vortices are important because they contribute to fire spread and are a threat to fire fighter safety. This paper documents field observations of the vortices and supplies supportive meteorological and fuel data. The discussion includes applicable laboratory and conceptual studies in fluid flow and heat transfer that may apply to vortex formation. 1. Introduction experiments showed that when air flowed parallel to a heated metal ribbon that simulated the flank of a crown The occurrence of vertical vortices in wildland fires fire, a thin, buoyant plume capped with a vortex pair has been well documented as well as mathematically developed above the ribbon along its length. -
D'arcy Wentworth Thompson
D’ARCY WENTWORTH THOMPSON Mathematically trained maverick zoologist D’Arcy Wentworth Thompson (May 2, 1860 –June 21, 1948) was among the first to cross the frontier between mathematics and the biological world and as such became the first true biomathematician. A polymath with unbounded energy, he saw mathematical patterns in everything – the mysterious spiral forms that appear in the curve of a seashell, the swirl of water boiling in a pan, the sweep of faraway nebulae, the thickness of stripes along a zebra’s flanks, the floret of a flower, etc. His premise was that “everything is the way it is because it got that way… the form of an object is a ‘diagram of forces’, in this sense, at least, that from it we can judge of or deduce the forces that are acting or have acted upon it.” He asserted that one must not merely study finished forms but also the forces that mold them. He sought to describe the mathematical origins of shapes and structures in the natural world, writing: “Cell and tissue, shell and bone, leaf and flower, are so many portions of matter and it is in obedience to the laws of physics that their particles have been moved, molded and conformed. There are no exceptions to the rule that God always geometrizes.” Thompson was born in Edinburgh, Scotland, the son of a Professor of Greek. At ten he entered Edinburgh Academy, winning prizes for Classics, Greek Testament, Mathematics and Modern Languages. At seventeen he went to the University of Edinburgh to study medicine, but two years later he won a scholarship to Trinity College, Cambridge, where he concentrated on zoology and natural science. -
An Overview of Impellers, Velocity Profile and Reactor Design
An Overview of Impellers, Velocity Profile and Reactor Design Praveen Patel1, Pranay Vaidya1, Gurmeet Singh2 1Indian Institute of Technology Bombay, India 1Indian Oil Corporation Limited, R&D Centre Faridabad Abstract: This paper presents a simulation operation and hence in estimating the approach to develop a model for understanding efficiency of the operating system. The the mixing phenomenon in a stirred vessel. The involved processes can be analysed to optimize mixing in the vessel is important for effective the products of a technology, through defining chemical reaction, heat transfer, mass transfer the model with adequate parameters. and phase homogeneity. In some cases, it is Reactors are always important parts of any very difficult to obtain experimental working industry, and for a detailed study using information and it takes a long time to collect several parameters many readings are to be the necessary data. Such problems can be taken and the system has to be observed for all solved using computational fluid dynamics adversities of environment to correct for any (CFD) model, which is less time consuming, inexpensive and has the capability to visualize skewness of values. As such huge amount of required data, it takes a long time to collect the real system in three dimensions. enough to build a model. In some cases, also it As reactor constructions and impeller is difficult to obtain experimental information configurations were identified as the potent also. Pilot plant experiments can be considered variables that could affect the macromixing as an option, but this conventional method has phenomenon and hydrodynamics, these been left far behind by the advent of variables were modelled. -
Sensitivity Analysis of Non-Linear Steep Waves Using VOF Method
Tenth International Conference on ICCFD10-268 Computational Fluid Dynamics (ICCFD10), Barcelona, Spain, July 9-13, 2018 Sensitivity Analysis of Non-linear Steep Waves using VOF Method A. Khaware*, V. Gupta*, K. Srikanth *, and P. Sharkey ** Corresponding author: [email protected] * ANSYS Software Pvt Ltd, Pune, India. ** ANSYS UK Ltd, Milton Park, UK Abstract: The analysis and prediction of non-linear waves is a crucial part of ocean hydrodynamics. Sea waves are typically non-linear in nature, and whilst models exist to predict their behavior, limits exist in their applicability. In practice, as the waves become increasingly steeper, they approach a point beyond which the wave integrity cannot be maintained, and they 'break'. Understanding the limits of available models as waves approach these break conditions can significantly help to improve the accuracy of their potential impact in the field. Moreover, inaccurate modeling of wave kinematics can result in erroneous hydrodynamic forces being predicted. This paper investigates the sensitivity of non-linear wave modeling from both an analytical and a numerical perspective. Using a Volume of Fluid (VOF) method, coupled with the Open Channel Flow module in ANSYS Fluent, sensitivity studies are performed for a variety of non-linear wave scenarios with high steepness and high relative height. These scenarios are intended to mimic the near-break conditions of the wave. 5th order solitary wave models are applied to shallow wave scenarios with high relative heights, and 5th order Stokes wave models are applied to short gravity waves with high wave steepness. Stokes waves are further applied in the shallow regime at high wave steepness to examine the wave sensitivity under extreme conditions. -
A Method of Constructing Phyllotaxically Arranged Modular Models by Partitioning the Interior of a Cylinder Or a Cone
A method of constructing phyllotaxically arranged modular models by partitioning the interior of a cylinder or a cone Cezary St¸epie´n Institute of Computer Science, Warsaw University of Technology, Poland [email protected] Abstract. The paper describes a method of partitioning a cylinder space into three-dimensional sub- spaces, congruent to each other, as well as partitioning a cone space into subspaces similar to each other. The way of partitioning is of such a nature that the intersection of any two subspaces is the empty set. Subspaces are arranged with regard to phyllotaxis. Phyllotaxis lets us distinguish privileged directions and observe parastichies trending these directions. The subspaces are created by sweeping a changing cross-section along a given path, which enables us to obtain not only simple shapes but also complicated ones. Having created these subspaces, we can put modules inside them, which do not need to be obligatorily congruent or similar. The method ensures that any module does not intersect another one. An example of plant model is given, consisting of modules phyllotaxically arranged inside a cylinder or a cone. Key words: computer graphics; modeling; modular model; phyllotaxis; cylinder partitioning; cone partitioning; genetic helix; parastichy. 1. Introduction Phyllotaxis is the manner of how leaves are arranged on a plant stem. The regularity of leaves arrangement, known for a long time, still absorbs the attention of researchers in the fields of botany, mathematics and computer graphics. Various methods have been used to describe phyllotaxis. A historical review of problems referring to phyllotaxis is given in [7]. Its connections with number sequences, e.g.