5.3 Inviscid instability mechanism of parallel flows

Size: px
Start display at page:

Download "5.3 Inviscid instability mechanism of parallel flows"

Transcription

1 1 Lecture Notes on Flui Dynamics 1.63J/2.21J) by Chiang C. Mei, Invisci instability mechanism of parallel flows We now turn to an oler problem of the instability of parallel flow without stratification an gravity, such as channel flows, jets, wakes an bounary layers Rayleigh s equation In a shear flow a necessary coniton for instablity is that there must be a point of inflection in the velocity profile U), i.e, 2 U 2 = 5.3.1) somewhere in the flow. Let us begin with an invisci uniform flow. U = U i, P = 5.3.2) an a an infinitesimal isturbance so that the total velocity an pressure fiels are u, w, p) =U + u,w,p ) 5.3.3) where primes enote isturbances which are functions of x,, t). The governing equations are u x + w = 5.3.4) ρu t +U + u )u x + w U + u )= p x 5.3.5) ρw t +U + u )v x + w w = p 5.3.6) The bounray conitions are w =, =, 5.3.7) where can be unboune. Lineariing the momentum equations by omitting terms quaratic in isiturbances, we get ρ u t + Uu x + U ) w = p x 5.3.8) ρw t + Uw x )= p 5.3.9) Eliminating p by cross-ifferntiation, ρ u t + Uu x + U ) w = p x = ρw t + Uw x ) x

2 2 Introucing the stream function by u = ψ,w = ψ x,weget ψ t + Uψ x U ) ψ x = ψ xt + Uψ xx ) x Let us consier a wave like isturbance ψ = f)e ikx Ct), where C = ω k 5.3.1) then ikcf + ikuf ik U ) f = ik k 2 Cf + k 2 Uf ) or U) C) f k 2 f ) f = ) which is known as Rayleigh s equation in hyroynamic instability. The bounfay oncitoins are: f =, =, ) Rayleigh s necessary conition for instability Let us rewrite ) as f k 2 f an multiply by the complex conjugate f f ) k 2 ff which can be rewritten by partial integration U C f = ) U C ff = f f ) f 2 k 2 f 2 U C f 2 = Integrating from =to = an applying the bounary conitions we get simply U C f 2 = f 2 + k 2 f 2) ) Now we let C = C r + ic i to be complex, then on the left-han sie, f 2 U C = U C ) f 2 = U C r + ic i ) f 2 U C 2 U C 2 Equation ) becomes U C r + ic i ) f 2 = f U C k 2 f 2) )

3 3 We get from the imaginary part of ), f 2 U C 2 = ) if C i. Thus for instability it is necessary that is positve for some y an negative elsewhere, i.e., must vanish somewhere in the flow, i.e., the velocity profile must have an inflection point insie the flow. This is Rayleigh s necessary conition for instability Physical explanation Lin, 1945) Figure 5.3.1: a) Velocity profile with an inflection point. b) Vorticity profile. c) Effect of exchanging flui parcels on vorticity In Figure we sketch the velocity profile with a inflection point at =anthe corresponing profile of vorticity ζ = U which has a maximum at =. Note the in the right-hane system the y axis points into the paper, therefore the flow of a positive vortex is clockwise. Consier three layers in the region where U) increases with monotonically, i.e., below the inflection point, with ζ 1 ) <ζ 2 ) <ζ 3 ). If a flui parcel escens from level 3 to level 2, it brings with its vorticity without change, accoring the vorticity transport law in an invisci flui. An excess vorticity is create at level 2 which tens to replace the flui on the right in layer 2 by flui with higher vorticity, an the flui on the left in layer 2 by flui with lower voriticity. The net consequence is to force the original excess vortex to return to layer 3. Similarly if a flui parcel ascens from level 1 to level 2, it brings with its

4 4 vorticity, hence creates an excess vorticity efect with tens to replace the flui on the left in layer 2 by flui with higher vorticity, an the flui on the right in layer 2 by flui with lower vorticity. The net consequence is to force the original efect votex to return to layer 3. Thus an acciental iplacement of a flui parcel tens to it original level; the flow is stable. By a similar reasoning, if U) ecreases monotonically in the flow is also stable. However, if there is a level of vorticity extremum, say level at =, then a flui element arriving at this layer is not force back to its origin. A flui element on one sie of level is equally at home on the opposite sie. The flow is unstable Fjortoft s stronger conition Further information can be obtaine from the real part of ), ) U U C r f 2 = f U C k 2 f 2) ) Now let U I be the velocity at the point of inflection, we use ) so that U I C r ) f 2 U C 2 = ) The ifference of ) an ) is U U U I ) f 2 = f U C k 2 f 2) < ) For this to hol it is necessary that th prouct U U I ) < 5.3.2) for all in the flow. For impoications of Fjortoft s theorem, see Figure Jets an wakes satisfy Rayleigth- Fjortoft criterion, but Poiseuille flows in a pipe an bounary layers on a flat plate o not. Why then is Poiseuille flow known to be unstable beyon Re = 21 Reynols)?

5 Figure 5.3.2: Fjortoft Theorm 5

5.3 Inviscid instability mechanism of parallel flows

5.3 Inviscid instability mechanism of parallel flows 1 Lecture Notes on Flui Dynamics 1.63J/2.21J) by Chiang C. Mei, 27 5.3 Invisci instability mechanism of parallel flows We now turn to an oler problem of the instability of parallel flow without stratification

More information

CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY

CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY 1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 00 CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY References: Drazin: Introduction to Hydrodynamic Stability Chandrasekar: Hydrodynamic

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Shear instabilities. Chapter Energetics of shear instabilities

Shear instabilities. Chapter Energetics of shear instabilities Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever

More information

Linear Hydrodynamic Stability Analysis Summary and Review. Example of Stability Problems. Start with Base Velocity Profile

Linear Hydrodynamic Stability Analysis Summary and Review. Example of Stability Problems. Start with Base Velocity Profile Linear Hydrodynamic Stability Analysis Summary and Review OCEN 678 Fall 2007 Scott A. Socolofsky Example of Stability Problems Shear Flows: flows with strong velocity gradients Classic examples Wakes Jets

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Effect of Rotation on Thermosolutal Convection. in a Rivlin-Ericksen Fluid Permeated with. Suspended Particles in Porous Medium

Effect of Rotation on Thermosolutal Convection. in a Rivlin-Ericksen Fluid Permeated with. Suspended Particles in Porous Medium Av. Theor. Appl. Mech., Vol. 3,, no. 4, 77-88 Effect of Rotation on Thermosolutal Convection in a Rivlin-Ericksen Flui Permeate with Suspene Particles in Porous Meium A. K. Aggarwal Department of Mathematics

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

S10.G.1. Fluid Flow Around the Brownian Particle

S10.G.1. Fluid Flow Around the Brownian Particle Rea Reichl s introuction. Tables & proofs for vector calculus formulas can be foun in the stanar textbooks G.Arfken s Mathematical Methos for Physicists an J.D.Jackson s Classical Electroynamics. S0.G..

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Completely passive natural convection

Completely passive natural convection Early View publication on wileyonlinelibrary.com (issue an page numbers not yet assigne; citable using Digital Object Ientifier DOI) ZAMM Z. Angew. Math. Mech., 1 6 (2011) / DOI 10.1002/zamm.201000030

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Stability of Stratified Couple-Stress Dusty Fluid in the Presence of Magnetic Field through Porous Medium

Stability of Stratified Couple-Stress Dusty Fluid in the Presence of Magnetic Field through Porous Medium vailable at http://pvamu.eu/aam ppl. ppl. Math. ISSN: 93-9466 Vol. 6, Issue (December ), pp. 5 5 pplications pplie Mathematics: n International Journal (M) Stability of Stratifie Couple-Stress Dusty Flui

More information

V q.. REASONING The potential V created by a point charge q at a spot that is located at a

V q.. REASONING The potential V created by a point charge q at a spot that is located at a 8. REASONING The electric potential at a istance r from a point charge q is given by Equation 9.6 as kq / r. The total electric potential at location P ue to the four point charges is the algebraic sum

More information

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10

DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3

1.2 - Stress Tensor Marine Hydrodynamics Lecture 3 13.021 Marine Hyroynamics, Fall 2004 Lecture 3 Copyright c 2004 MIT - Department of Ocean Engineering, All rights reserve. 1.2 - Stress Tensor 13.021 Marine Hyroynamics Lecture 3 Stress Tensor τ ij:. The

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Thermal Modulation of Rayleigh-Benard Convection

Thermal Modulation of Rayleigh-Benard Convection Thermal Moulation of Rayleigh-Benar Convection B. S. Bhaauria Department of Mathematics an Statistics, Jai Narain Vyas University, Johpur, Inia-3400 Reprint requests to Dr. B. S.; E-mail: bsbhaauria@reiffmail.com

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

Problem set 2: Solutions Math 207B, Winter 2016

Problem set 2: Solutions Math 207B, Winter 2016 Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1

More information

UNIVERSITY OF EAST ANGLIA

UNIVERSITY 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 information

The effect of dissipation on solutions of the complex KdV equation

The effect of dissipation on solutions of the complex KdV equation Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,

More information

General introduction to Hydrodynamic Instabilities

General 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 information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

Centrum voor Wiskunde en Informatica

Centrum voor Wiskunde en Informatica Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

ME338A CONTINUUM MECHANICS

ME338A CONTINUUM MECHANICS global vs local balance equations ME338A CONTINUUM MECHANICS lecture notes 11 tuesay, may 06, 2008 The balance equations of continuum mechanics serve as a basic set of equations require to solve an initial

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem

(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem Professor Fred Stern Fall 2014 1 Chapter 6: Viscous Flow in Ducts 6.2 Stability and Transition Stability: can a physical state withstand a disturbance and still return to its original state. In fluid mechanics,

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

4 Shear flow instabilities

4 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 information

Prototype Instabilities

Prototype Instabilities Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall. Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,

More information

SOLUTIONS for Homework #3

SOLUTIONS for Homework #3 SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical

More information

Conservation and Balance Equations

Conservation and Balance Equations Chapter 2 Conservation an Balance Equations In this chapter we consier some applications of Reynol s transport theorem, Theorem 1.1. For a balance equation of the general type u(t,y)y f(t,y)y (2.1) t we

More information

Balance laws on domains with moving interfaces. The enthalpy method for the ice melting problem.

Balance laws on domains with moving interfaces. The enthalpy method for the ice melting problem. Balance laws on omains with moving interfaces. The enthalpy metho for the ice melting problem. Gowin Kakuba 12th March 2008 Outline 1 Balance equations on omains with moving bounaries Introuction Rankine-Hugoniot

More information

The Eady problem of baroclinic instability described in section 19a was shown to

The Eady problem of baroclinic instability described in section 19a was shown to 0. The Charney-Stern Theorem The Eady problem of baroclinic instability described in section 19a was shown to be remarkably similar to the Rayleigh instability of barotropic flow described in Chapter 18.

More information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN

More information

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling.

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling. Chapter 8 Analytic Functions Stuents nee encouragement. So if a stuent gets an answer right, tell them it was a lucky guess. That way, they evelop a goo, lucky feeling. 1 8.1 Complex Derivatives -Jack

More information

Problem Set 6: Workbook on Operators, and Dirac Notation Solution

Problem Set 6: Workbook on Operators, and Dirac Notation Solution Moern Physics: Home work 5 Due ate: 0 March. 014 Problem Set 6: Workbook on Operators, an Dirac Notation Solution 1. nswer 1: a The cat is being escribe by the state, ψ >= ea > If we try to observe it

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

MAE 210A FINAL EXAM SOLUTIONS

MAE 210A FINAL EXAM SOLUTIONS 1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix

More information

6 Two-layer shallow water theory.

6 Two-layer shallow water theory. 6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour

More information

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.

Candidates 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 information

constant ρ constant θ

constant ρ constant θ Answers to Problem Set Number 3 for 8.04. MIT (Fall 999) Contents Roolfo R. Rosales Λ Boris Schlittgen y Zhaohui Zhang ẓ October 9, 999 Problems from the book by Saff an Snier. 2. Problem 04 in section

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

State-Space Model for a Multi-Machine System

State-Space Model for a Multi-Machine System State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x =

cosh x sinh x So writing t = tan(x/2) we have 6.4 Integration using tan(x/2) 2t 1 + t 2 cos x = 1 t2 sin x = 6.4 Integration using tan/ We will revisit the ouble angle ientities: sin = sin/ cos/ = tan/ sec / = tan/ + tan / cos = cos / sin / tan = = tan / sec / tan/ tan /. = tan / + tan / So writing t = tan/ we

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

Advanced Partial Differential Equations with Applications

Advanced Partial Differential Equations with Applications MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.

More information

Curvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor

Curvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor Curvature, Conformal Mapping, an 2D Stationary Flui Flows Michael Taylor 1. Introuction Let Ω be a 2D Riemannian manifol possibly with bounary). Assume Ω is oriente, with J enoting counterclockwise rotation

More information

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2) 1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following

More information

Vorticity in natural coordinates

Vorticity in natural coordinates Vorticity in natural coordinates (see Holton pg 95, section 4.2.) Let s consider the vertical vorticity component only, i.e. ζ kˆ ω, we have ω u dl kˆ ω lim --- lim ----------------- curve is in xy plane

More information

1. Comparison of stability analysis to previous work

1. Comparison of stability analysis to previous work . Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N

1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N Physics 15 Solution Set Winter 018 1. At time t = 0, the wave function of a free particle moving in a one-imension is given by, ψ(x,0) = N where N an k 0 are real positive constants. + e k /k 0 e ikx k,

More information

Jet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities. Mark Fruman

Jet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities. Mark Fruman p. 1/24 Jet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities Mark Fruman Bach Lien Hua, Richard Schopp, Marc d Orgeville, Claire Ménesguen LPO IFREMER, Brest, France IAU

More information

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001

More information

9. 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 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 information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

All s Well That Ends Well: Supplementary Proofs

All s Well That Ends Well: Supplementary Proofs All s Well That Ens Well: Guarantee Resolution of Simultaneous Rigi Boy Impact 1:1 All s Well That Ens Well: Supplementary Proofs This ocument complements the paper All s Well That Ens Well: Guarantee

More information

MAE 210C FLUID MECHANICS III SPRING 2015 Rayleigh s equation for axisymmetric flows Forthestabilityanalysisofparallelaxisymmetricflowswithbasicvelocity profileū = (U x,u r,u θ ) = (U(r),0,0) weintroducemodalperturbationsoftheform(v

More information

2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con

2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con A MIXED FINITE ELEMENT METHOD FOR A THIRD ORDER PARTIAL DIFFERENTIAL EQUATION G. Li 1 an F. A. Milner 2 A mixe finite element metho is escribe for a thir orer partial ifferential equation. The metho can

More information

6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) instability

6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) instability 1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 6.6 Rayleigh-Darcy (or Horton-Rogers-Lapwood) instability in a porous layer 6-6-Lapwood.tex May 11, 3 Nield & Bejan, Chapter 6 Convection

More information

Gravitation & Cosmology. Exercises # µ x = 0 (1)

Gravitation & Cosmology. Exercises # µ x = 0 (1) Gravitation & Cosmology. Exercises # 4.1 - Geoesics a) Show that the Euler-Lagrange equations for the Lagrangian L τ ẋ L µ x = 0 (1) µ L = 1 2 g µνẋ µ ẋ ν (2) are the geoesic equations where, as usual,

More information

Turbulence - Theory and Modelling GROUP-STUDIES:

Turbulence - Theory and Modelling GROUP-STUDIES: Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence

More information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618

More information

1 Lecture 18: The chain rule

1 Lecture 18: The chain rule 1 Lecture 18: The chain rule 1.1 Outline Comparing the graphs of sin(x) an sin(2x). The chain rule. The erivative of a x. Some examples. 1.2 Comparing the graphs of sin(x) an sin(2x) We graph f(x) = sin(x)

More information

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes 1 Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA March 2003

More information

MATHEMATICS BONUS FILES for faculty and students

MATHEMATICS BONUS FILES for faculty and students MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department

More information

Physics 2212 GJ Quiz #4 Solutions Fall 2015

Physics 2212 GJ Quiz #4 Solutions Fall 2015 Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.

More information

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS

COUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND

More information

Turbulence Modeling I!

Turbulence 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 information

Linear stability analysis of coupled parallel flexible plates in an axial flow

Linear stability analysis of coupled parallel flexible plates in an axial flow Linear stability analysis of couple parallel flexible plates in an axial flow Sébastien Michelin,a,b, Stefan G. Llewellyn Smith a a Department of Mechanical an Aerospace Engineering, Jacobs School of Engineering,

More information

Lecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev

Lecture 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 information

Making a Wavefunctional representation of physical states congruent with the false vacuum hypothesis of Sidney Coleman

Making a Wavefunctional representation of physical states congruent with the false vacuum hypothesis of Sidney Coleman Making a Wavefunctional representation of physical states congruent with the false vacuum hypothesis of Siney Coleman A. W. Beckwith Department of Physics an Texas Center for Superconuctivity an Avance

More information

Geometric Algebra Approach to Fluid Dynamics

Geometric Algebra Approach to Fluid Dynamics Geometric Algebra Approach to Flui Dynamics Carsten Cibura an Dietmar Hilenbran Abstract In this work we will use geometric algebra to prove a number of well known theorems central to the fiel of flui

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

[Kumar*, 5(2): February, 2016] ISSN: (I2OR), Publication Impact Factor: 3.785

[Kumar*, 5(2): February, 2016] ISSN: (I2OR), Publication Impact Factor: 3.785 [Kumar*, 5(): February, 6] ISSN: 77-9655 (IOR), Publication Impact Factor: 3.785 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY THERMOSOLUTAL CONVECTION IN A HETEROGENEOUS VISCO-ELASTIC

More information

TEST 2 (PHY 250) Figure Figure P26.21

TEST 2 (PHY 250) Figure Figure P26.21 TEST 2 (PHY 250) 1. a) Write the efinition (in a full sentence) of electric potential. b) What is a capacitor? c) Relate the electric torque, exerte on a molecule in a uniform electric fiel, with the ipole

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information