α 2 )(k x ũ(t, η) + k z W (t, η)
|
|
- Raymond Francis
- 5 years ago
- Views:
Transcription
1 1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is infinite in the x and z directions. It is bounded at y = by a wall and also at y = 1. In our scenario, this fluid is incompressible and has a constant density. What I m going to do is show you the steps to go from the full uncontrolled flow governed by the Navier Stokes equations to a controlled flow that looks like its governed by heat equations. The controllers we re going to use to do this are the following: Note they are in Fourier Space Ṽc t = 1 Re Ṽ yy 1 Ṽyy α 2 Ṽ c β Ṽ t, η2η 1 coshα1 ηdη, 1 = 1 2πi k x ũy t, ũ y t, 1 + k z Wy t, Re W y t, 1 α 2 Ṽ c β Ṽ t, η2η 1 coshα1 ηdη, 2 Ũ c = 4π2 α 2 k x Ỹ t, 1 + k z ωt, 1 4π 2 = α 2 k x K1, η + k z Γ1, η W c = 4π2 α 2 k z Ỹ y, 1 k x ωt, 1 = 4π 2 α 2 k z K1, η k x Γ1, η and now for the steps to derive these controllers. k x ũt, η + k z W t, η k x ũt, η + k z W t, η 3 dη 4 5 dη. 6 2 Notation We re going to start with notation just to make sure everyone is on the same page. index notation, vector calculus grad, div, laplacian, material derivative, 1
2 2-D Fourier Transform and diff property. = [ x z y ] = 7 2 = = xx + zz + yy = 2 x j x j 9 D Dt Here are some useful examples: 8 1 = t + u. = t + u i 11 p = p = [p x p z p y ] 12.u = u i = u x + w z + v y 13 2 u = 2 u i = u xx + u zz + u yy w xx + w zz + w yy 14 x j x j v xx + v zz + v yy u. u = u j u i x j = uu x + wu z + vu y uw x + ww z + vw y uv x + wv z + vv y 15 We will be using the 2-D Fourier Transform, the forward transform looks like this: F t, k x, k z, y = i nftyf t, x, z, ye 2πikxx+kzz dxdz 16 Recall the useful property of the FT under differentiation F x = 2πik x F 17 F z = 2πik z F 18 Remember this only works for differentiation in the x and z directions as there are the directions in which we are taking the FT. 2
3 3 Navier Stokes and our equilibrium So we start with the Navier Stokes equations Du Dt = 1 Re 2 u p 19 u t = 1 Re 2 u u. u p 2.u = 21 uy = = 22 uy = 1 = 1 23 Remember 19 is the same as 2, only we have written out the material derivative operator. Note that these equations are already non dimensionalized. We assumed density is constant, and that gravity and other forces are neglected. We are including viscosity and then we have a pressure gradient to force the flow. The Reynolds number Re is a ratio of the inertial forces to the viscous forces. It is the characteristic velocity times the characteristic length divided by the kinematic viscosity. Re = UL ν The second equation you see is the incompressible condition. This means the density of a parcel of fluid is constant. It does not mean the density of the flow is constant. We need the additional assumption that density is constant to drop it altogether. Our boundary conditions mean two things. One is that the flow cannot penetrate the walls at the top and bottom, so the normal velocity V must be zero. The other is that there is no slip at the walls, ie there is viscosity, so the tangential velocities U and W must be zero at the boundaries. We can use 21 to derive an equation for the pressure. We take the divergence of the velocity equations.u t, plug in 21, and we arrive at the following equation for the pressure. We get the boundary conditions by looking at the V t equation and plugging in and 1..u t = 1 Re. 2 u.u. u. p 24.u t = 1 Re 2.u u j u i x j 2 p 25 = u j u i + u j x j 2 u i p 26 x j 2 p x j x j = u j u i x j 27 p = U 2 x + W 2 z + V 2 y + 2U z W x + 2U y V y + 2W y V z 28 p y = 1 Re V yy 29 p y 1 = 1 Re V yy1 3 3
4 We said we were looking at the Poiseuille flow. What is that? well, it means we have a flow which is parallel to the walls and is an equilibrium soln to the NS eqns. In our case, this flow looks like this U e = 4y1 y 31 V = W = 32 P e = P 8 Re x 33 When the system remains at this equilibrium, the flow is smooth - laminar. However, when the Reynolds number is high, the flow is unstable and we do not remain at this equilibrium. First we see delaminarization into a transition regime, and for higher Re the flow goes fully turbulent. We want to stabilize the system around the Poiseuille flow so that turbulence is reduced. 4 linearization around the equilibrium Now we shift the system so that we are looking at just the perturbations and then we linearize this resulting system; Our perturbation variables are the following u = U U e, p = P P e, 34 4
5 And now we plug them in. First we shift the system. u + U e t = 1 Re u + U e xx + u + U e zz + u + U e yy u + U e u + U e x V u + U e y W u + U e z p + P e x, 35 u t = 1 Re u xx + u zz + u yy U e yy u + U e u x V u y + U e y W u z p x P 3 x, 36 u t = 1 Re u xx + u zz + u yy u + U e u x V u y + U e y W u z p x, 37 W t = 1 Re W xx + W zz + W yy u + U e W x V W y W W z p + P e z, 38 V t = 1 Re V xx + V zz + V yy u + U e V x V V y W V z p + P e y, 39 u + U e x + W z + V y = 4 u x + W z + V y = 41 p + P e xx + p + P e zz + p + P e yy = 2V 2 y 2u + U e x W z 2u + U e y V x 2u + U e z W x 2V z W y,42 p xx + p zz + p yy = 2V 2 y 2u x W z 2u y + U e y V x 2u z W x 2V z W y, 43 And now we linearize it u t = 1 Re u xx + u zz + u yy U e u x U e y V p x, 44 W t = 1 Re W xx + W zz + W yy U e W x p z, 45 V t = 1 Re V xx + V zz + V yy U e V x p y, 46 p xx + p zz + p yy = 2U e y V x 47 u x + W z + V y = 48 5 stability of linear system Now this linearized system is still unstable. To see this, we look at two coupled equations. To find the first equation, we take the Laplacian of the V t equation 5
6 46 and use the pressure equation 47 to make it autonomous. The second equation consists of the vorticity in the y direction u z W x. The resulting equations look like this: 2 V t = 1 Re 4 V U 2 V x + U yy V x 49 ω t = 1 Re 2 ω U e ω x U e y V z 5 These are the equations from which the Orr-Sommerfeld 197,198 and Squire1933 equations are derived. For those of you who have studied flow stability - take the normal modes of the form V = ṽye iαx+βz σt 51 ω = ωye iαx+βz σt 52 and you arrive at the OS and Squire equations. h iσ + iαud 2 α 2 + β 2 iαu 1 `D2 i α 2 + β 2 2 ṽ = 53 Re h iσ iαu 1 `D2 α 2 + β i 2 ω = iβu ṽ 54 Re Notice that the Squire equation 5 54 is forced by the OS equation Without that forcing, the Squire equation is stable. This leads to the theorum, that for any unstable 3-D flow, there exists a 2-D flow at a *lower* Reynolds number which is also unstable. Now we are going to use 5 to find our controllers, but instead of using 49 we are going to do something similar to what we did in the beam example, as we dont like those higher order derivatives. What we shall do instead, is use V y and p similar to how we used α and u in the beam example. 6 Fourier Transform We are also not going to take normal modes to find out controllers, but we are going to do something similar. Instead of making assumptions in both space and time, we will only make assumptions in space. We shall take the fourier transform of the system in the x and z directions, as our system is infinite in these directions. This gives us wavenumbers in these directions, k x and k z. Why do we do this? Well, we started with 1 system defined by 4 3-D nonlinear pdes. Now, with linearization and the FT, we have an infinite number of systems, each paramaterized by the wave numbers k x and k z, and defined by 4 1-D linear pdes. Yeah! - we only have to look at 1-D linear systems. This approach is not just used for flow control, we can use it in other applications, for example, if one had an infinite or periodic beam - one could use the Fourier Transform. 6
7 After taking the Fourier Transform of the previous linear system 44 48, the resulting equations are α 2 = 4π 2 k 2 x + k 2 z, β = 16πk x i, 55 u t = 1 Re α2 u + u yy + β 2 yy 1u + 42y 1V 2πik xp, 56 W t = 1 Re α2 W + W yy + β 2 yy 1W 2πik zp, 57 V t = 1 Re α2 V + V yy + β 2 yy 1V p y, 58 uy = = V y = = W y = = 59 uy = 1 = U c 6 W y = 1 = W c 61 V y = 1 = V c 62 α 2 p + p yy = β2y 1V, 63 p y = 1 Re V yy 64 p y 1 = 1 Re α2 V c + V yy 1 V c t 65 k x u + k z W = V y 2πi 66 7 Solution to p and V c So we want to find our 2 equations one for V y and p and then other, 5. However, first, lets solve the p equation. The p equation, is nondynamic - similar to the α equation in the beam exercise. We will now solve this equation and substitute it back into the other 3 equations. Looking at our equation for p we see α 2 p + p yy = β2y 1V, 67 p y = 1 Re V yy 68 p y 1 = 1 Re α2 V c + V yy 1 V c t 69 7
8 We have an ODE in p forced by V. The solution to this ODE and p is { y p = β V t, η2η 1 sinhαy ηdη β coshαy V t, η2η 1 coshα1 ηdη sinhα coshα1 y p y t, x, z, sinhα + coshαy } sinhα p yt, x, z, 1, 7 y p = β V t, η2η 1 sinhαy ηdη { + coshαy β V t, η2η 1 coshα1 ηdη sinhα } + 1 Re α2 V c + V yy 1 V c t coshα1 y p y t, x, z, 71 sinhα This next step is again similar to the beam case, we design one controller to put the system into a strict feedback form so that we can use backstepping, ie we dont like the V 2η 1 coshα1 ηdη term, so lets get rid of it. How do we do this? Notice the V c and V c t terms in the equation. We will use these to get rid of our unwanted term by setting V c in the following way. V yy 1 V yy α 2 V c V c t = 1 Re β V t, η2η 1 coshα1 ηdη, 72 Notice we added a V yy term which we have to account for in our resulting p equation. p = β y V t, η2η 1 sinhαy ηdη + coshαy 1 sinhα Re V yy coshα1 y 1 sinhα Re V yy 73 8 Y and ω and backstepping 8
Stability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More information4 Shear flow instabilities
4 Shear flow instabilities Here we consider the linear stability of a uni-directional base flow in a channel u B (y) u B = 0 in the region y [y,y 2. (35) 0 y = y 2 y z x u (y) B y = y In Sec. 4. we derive
More informationSG2221 Wave Motion and Hydrodinamic Stability. MATLAB Project on 2D Poiseuille Flow Alessandro Ceci
SG2221 Wave Motion and Hydrodinamic Stability MATLAB Project on 2D Poiseuille Flow Alessandro Ceci Base Flow 2D steady Incompressible Flow Flow driven by a constant pressure gradient Fully developed flow
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationStability Analysis of a Plane Poisseuille Flow. Course Project: Water Waves and Hydrodynamic Stability 2016 Matthias Steinhausen
Stability Analysis of a Plane Poisseuille Flow Course Project: Water Waves and Hydrodynamic Stability 2016 Matthias Steinhausen y Plane Poiseuille Flow (PPF): Characteristics y 1 0.5 Velocity profile for
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationPerturbation dynamics in laminar and turbulent flows. Initial value problem analysis
Perturbation dynamics in laminar and turbulent flows. Initial value problem analysis Francesca De Santi 1 1 Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Italy 15th of April
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationFeedback on vertical velocity. Rotation, convection, self-sustaining process.
Feedback on vertical velocity. Rotation, convection, self-sustaining process. Fabian Waleffe//notes by Andrew Crosby, Keiji Kimura, Adele Morrison Revised by FW WHOI GFD Lecture 5 June 24, 2011 We have
More informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationLinear stability of plane Poiseuille flow over a generalized Stokes layer
Linear stability of plane Poiseuille flow over a generalized Stokes layer Maurizio Quadrio 1,, Fulvio Martinelli and Peter J Schmid 1 Dip. Ing. Aerospaziale, Politecnico di Milano, Campus Bovisa, I-0156
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More information9. Boundary layers. Flow around an arbitrarily-shaped bluff body. Inner flow (strong viscous effects produce vorticity) BL separates
9. Boundary layers Flow around an arbitrarily-shaped bluff body Inner flow (strong viscous effects produce vorticity) BL separates Wake region (vorticity, small viscosity) Boundary layer (BL) Outer flow
More informationInitial-value problem for shear flows
Initial-value problem for shear flows 1 Mathematical Framework The early transient and long asymptotic behaviour is studied using the initialvalue problem formulation for two typical shear flows, the plane
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationROUTES TO TRANSITION IN SHEAR FLOWS
ROUTES TO TRANSITION IN SHEAR FLOWS Alessandro Bottaro with contributions from: D. Biau, B. Galletti, I. Gavarini and F.T.M. Nieuwstadt ROUTES TO TRANSITION IN SHEAR FLOWS Osborne Reynolds, 1842-1912 An
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationStability of Shear Flow: part 2
Stability of Shear Flow: part 2 Fabian Waleffe notes by Matthew Chantry and Lindsey Ritchie Revised by FW WHOI GFD Lecture 4 June 23, 2011 We derive necessary conditions for linear instability of shear
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationLaminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations
Laminar Boundary Layers Answers to problem sheet 1: Navier-Stokes equations The Navier Stokes equations for d, incompressible flow are + v ρ t + u + v v ρ t + u v + v v = 1 = p + µ u + u = p ρg + µ v +
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationThe Orchestra of Partial Differential Equations. Adam Larios
The Orchestra of Partial Differential Equations Adam Larios 19 January 2017 Landscape Seminar Outline 1 Fourier Series 2 Some Easy Differential Equations 3 Some Not-So-Easy Differential Equations Outline
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationHydrodynamic stability: Stability of shear flows
Hydrodynamic stability: Stability of shear flows November 13, 2009 2 Contents 1 Introduction 5 1.1 Hystorical overview......................................... 5 1.2 Governing equations........................................
More informationLecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev
Lecture Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev Introduction In many cases in nature, like in the Earth s atmosphere, in the interior of stars and planets, one sees the appearance
More informationFLOW-NORDITA Spring School on Turbulent Boundary Layers1
Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationÉCOLE DE PHYSIQUE DES HOUCHES. Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at
ÉCOLE DE PHYSIQUE DES HOUCHES STABILITY OF FLUID FLOWS Carlo COSSU Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at Département de Mécanique, École Polytechinique http://www.imft.fr/carlo-cossu/
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationInfo. No lecture on Thursday in a week (March 17) PSet back tonight
Lecture 0 8.086 Info No lecture on Thursday in a week (March 7) PSet back tonight Nonlinear transport & conservation laws What if transport becomes nonlinear? Remember: Nonlinear transport A first attempt
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationLinear stability of plane Poiseuille flow over a generalized Stokes layer
Linear stability of plane Poiseuille flow over a generalized Stokes layer M. Quadrio, Fulvio Martinelli, Peter Schmid To cite this version: M. Quadrio, Fulvio Martinelli, Peter Schmid. Linear stability
More informationTopics in Fluid Dynamics: Classical physics and recent mathematics
Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationThe Generalized Boundary Layer Equations
The Generalized Boundary Layer Equations Gareth H. McKinley (MIT-HML) November 004 We have seen that in general high Reynolds number flow past a slender body such as an airfoil can be considered as an
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More information0 = p. 2 x + 2 w. z +ν w
Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationarxiv: v4 [physics.comp-ph] 21 Jan 2019
A spectral/hp element MHD solver Alexander V. Proskurin,Anatoly M. Sagalakov 2 Altai State Technical University, 65638, Russian Federation, Barnaul, Lenin prospect,46, k2@list.ru 2 Altai State University,
More informationThe Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations
The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationPAPER 331 HYDRODYNAMIC STABILITY
MATHEMATICAL TRIPOS Part III Thursday, 6 May, 016 1:30 pm to 4:30 pm PAPER 331 HYDRODYNAMIC STABILITY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal
More informationBäcklund Transformations: a Link Between Diffusion Models and Hydrodynamic Equations
Copyright 2014 Tech Science Press CMES, vol.103, no.4, pp.215-227, 2014 Bäcklund Transformations: a Link Between Diffusion Models and Hydrodynamic Equations J.R. Zabadal 1, B. Bodmann 1, V. G. Ribeiro
More informationApplied Mathematics 505b January 22, Today Denitions, survey of applications. Denition A PDE is an equation of the form F x 1 ;x 2 ::::;x n ;~u
Applied Mathematics 505b January 22, 1998 1 Applied Mathematics 505b Partial Dierential Equations January 22, 1998 Text: Sobolev, Partial Dierentail Equations of Mathematical Physics available at bookstore
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More information7 EQUATIONS OF MOTION FOR AN INVISCID FLUID
7 EQUATIONS OF MOTION FOR AN INISCID FLUID iscosity is a measure of the thickness of a fluid, and its resistance to shearing motions. Honey is difficult to stir because of its high viscosity, whereas water
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationChapter 1: Basic Concepts
What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Stokes' law and Reynold number Dr. Zifei Liu The motion of a particle in a fluid environment, such as air or water m dv =F(t) - F dt d - 1 4 2 3 πr3
More informationCriteria for locally fully developed viscous flow
1 Criteria for locally fully developed viscous flow Ain A. Sonin, MIT October 00 Contents 1. Locally fully developed flow.. Criteria for locally fully developed flow. 3 3. Criteria for constant pressure
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationVISCOUS FLOW DUE TO A SHRINKING SHEET
QUARTERLY OF APPLIED MATHEMATICS VOLUME, NUMBER 0 XXXX XXXX, PAGES 000 000 S 0000-0000(XX)0000-0 VISCOUS FLOW DUE TO A SHRINKING SHEET By M. MIKLAVČIČ (Department of Mathematics, Michigan State University,
More informationThin airfoil theory. Chapter Compressible potential flow The full potential equation
hapter 4 Thin airfoil theory 4. ompressible potential flow 4.. The full potential equation In compressible flow, both the lift and drag of a thin airfoil can be determined to a reasonable level of accuracy
More informationLecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s
Lecture 19: Solving linear ODEs + separable techniques for nonlinear ODE s Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationLINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP
MAGNETOHYDRODYNAMICS Vol. 48 (2012), No. 1, pp. 147 155 LINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP N. Vetcha, S. Smolentsev, M. Abdou Fusion Science and Technology Center, UCLA,
More informationRecent progress on the vanishing viscosity limit in the presence of a boundary
Recent progress on the vanishing viscosity limit in the presence of a boundary JIM KELLIHER 1 1 University of California Riverside The 7th Pacific RIM Conference on Mathematics Seoul, South Korea 1 July
More informationd v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =
Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationBorel Summability in PDE initial value problems
Borel Summability in PDE initial value problems Saleh Tanveer (Ohio State University) Collaborator Ovidiu Costin & Guo Luo Research supported in part by Institute for Math Sciences (IC), EPSRC & NSF. Main
More informationShear Turbulence. Fabian Waleffe. Depts. of Mathematics and Engineering Physics. Wisconsin
Shear Turbulence Fabian Waleffe Depts. of Mathematics and Engineering Physics, Madison University of Wisconsin Mini-Symposium on Subcritical Flow Instability for training of early stage researchers Mini-Symposium
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationCandidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More information1. Fluid properties and concept of continuum
1. Fluid properties and concept of continuum Fluid: a substance that deforms continuously under an applied shear stress (e.g., air, water, upper mantle...) Not practical/possible to treat fluid mechanics
More informationexamples of equations: what and why intrinsic view, physical origin, probability, geometry
Lecture 1 Introduction examples of equations: what and why intrinsic view, physical origin, probability, geometry Intrinsic/abstract F ( x, Du, D u, D 3 u, = 0 Recall algebraic equations such as linear
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More information