Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.

Similar documents
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Turbulent Flow. Turbulent Flow

Turbulent Transport in Single-Phase Flow. Peter Bernard, University of Maryland

THE COUPLED LES - SUBGRID STOCHASTIC ACCELERATION MODEL (LES-SSAM) OF A HIGH REYNOLDS NUMBER FLOWS

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

PHYS 705: Classical Mechanics. Newtonian Mechanics

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

In this section is given an overview of the common elasticity models.

Lecture 20: Noether s Theorem

Problem 1: To prove that under the assumptions at hand, the group velocity of an EM wave is less than c, I am going to show that

Analysis of the two-point velocity correlations. in turbulent boundary layer ows. By M. Oberlack

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models

2 Finite difference basics

Introduction to Turbulence Modeling

Mathematical Preparations

CHAPTER 14 GENERAL PERTURBATION THEORY

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Rigid body simulation

1. Why turbulence occur? Hydrodynamic Instability. Hydrodynamic Instability. Centrifugal Instability: Rayleigh-Benard Instability:

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6

Turbulent Flows and their Numerical Modeling

FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS

Thermodynamics General

PHYS 705: Classical Mechanics. Canonical Transformation II

6. Stochastic processes (2)

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

6. Stochastic processes (2)

Integrals and Invariants of Euler-Lagrange Equations

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Shock Acceleration at an Interplanetary Shock: A Focused Transport Approach

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

Turbulence and its Modelling

Lecture 10 Support Vector Machines II

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Mechanics Physics 151

Lecture 12. Modeling of Turbulent Combustion

8.6 The Complex Number System

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Advanced Quantum Mechanics

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Implicit Integration Henyey Method

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

The Feynman path integral

Turbulence Modeling in Computational Fluid Dynamics (CFD) Jun Shao, Shanti Bhushan, Tao Xing and Fred Stern

Polynomial Regression Models

Mechanics Physics 151

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

NUMERICAL DIFFERENTIATION

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Lecture 5.8 Flux Vector Splitting

Formal solvers of the RT equation

Flow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved

Canonical transformations

EXPERIMENTAL STUDY OF NEAR WALL TURBULENCE USING PIV

Mechanics Physics 151

Introduction to Turbulence. by Håkan Gustavsson Division of Fluid Mechanics Luleå University of Technology

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Lecture 14: Forces and Stresses

2) For a two-dimensional steady turbulent flow in Cartesian coordinates (x,y), with mean velocity components (U,V), write

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)

Numerical Heat and Mass Transfer

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

Macroscopic Momentum Balances

Waveguides and resonant cavities

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

One Dimension Again. Chapter Fourteen

The Geometry of Logit and Probit

PHY2049 Exam 2 solutions Fall 2016 Solution:

Lecture 19. Endogenous Regressors and Instrumental Variables

Lagrangian Field Theory

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Differentiating Gaussian Processes

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Difference Equations

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

PHYS 705: Classical Mechanics. Calculus of Variations II

), it produces a response (output function g (x)

Lecture Notes on Linear Regression

Feature Selection: Part 1

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Integrals and Invariants of

Passive Scalar and Scalar Flux in Homogeneous Anisotropic Turbulence

Quantifying Uncertainty

Multicomponent Flows

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Iterative General Dynamic Model for Serial-Link Manipulators

Lecture Note 3. Eshelby s Inclusion II

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March

Extracting turbulent spectral transfer from under-resolved velocity fields

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Tools for large-eddy simulation

Lecture 3: Probability Distributions

Transcription:

Turbulence Lecture 1 Non-lnear Dynamcs Strong non-lnearty s a key feature of turbulence. 1. Unstable, chaotc behavor.. Strongly vortcal (vortex stretchng) 3 s & 4 s Taylor s work on homogeneous turbulence Kolmogorov. 6 s & 7 s Krachnan and followers addressed non-lnear stochastc processes. Deductve theory hasn t been developed. 1. Learned much about non-lnear dynamcs.. Theores are based on ad hoc assumptons. Smplfy the problem and consder homogeneous turbulence. - Contans essental physcs of the non-lnear dynamcs. - Smple as possble. Homogeneous n all spatal drectons. If t s then: a. Energy exchange wth mean flow. b. Dffuson terms are not mportant. Lose some large structures pathologcal (?) flow. (Could be too much of a smplfcaton.) Ideas at least apply to small-scale n hgh Re no. flows. Can generate a homogeneous flow n a lab. 1

Screen; mesh sze, small compared to D. / U Possbly wth constant mean shear to provde energy to perturbatons downstream. Not homogeneous 1. Taylor ntroduced the concepts n the 3 s and developed the experment to study t.. Also can study homogeneous turbulence on the computer. Solve the N.-S. equatons drectly at a hgh enough Re. No. Re 1, λ realzable based on Taylor mcroscale. 3. If the Re No. of the flow s hgh enough, the small scale structures of the flow are homogeneous. A. Homogeneous turbulence can be produced n the presence of background shear. = λx (Corrsn, Champagne) u1 B. Background densty stratfcaton 1 ρ = constant. ρ z C. Magnetc feld. Consder the case where U =, P = Mean momentum equaton says = Choose a coordnate system attached to mean flow gves U =. There s an obvous non-homogenety n x.

L (Integral scale n x) Ut x T s decay tme of turbulence. T L x = u \ So Turbulence level u 1 Then approx. homogeneous n x. U x = Ut Relate lab, spatally decayng flud to computer, tme decayng turb. smulaton. For homogeneous turb. the 1-pt. energy equaton s: u t = µ ss = ε ε - dssppaton rate. It can only decay. The equaton hdes the reason for decay (.e., vortex stretchng) Need -pt statstcs, go to -pt measurements to understand the dynamcs. For hgh Re. No. homogeneous turbulence scalng s the same. ε = E u u tme u 3 u u K approxmately ndependent of R. The mean square vortcty (enstrophy) equaton reduces to ω t ω = ωω s v x / producton dsspaton long tme \ / scale short tme scales 3

Scalng leads to terms on RHS beng n balance. Note: A couple of references for homogeneous turbulence: Batchelor, Monn & Yaglom Vol., Hnze, Leseur. ω can actually grow although t must eventually decay. E( Κ) Κ E / dsspaton rate Κ Startng pont pt spatal correlaton. (, ) = (, ) ( + ) R rt u xt u x rt, or equvalently, the knetc energy spectrum. Homogeneous n x, no mean flow, want an equaton for R ( rt, ) Defne x = x+ r Prme means dealng wth a dfferent pont n space. u s u at x 4

x - equaton: u 1 p = uu + v u t x x at x : k k ρ u 1 p = + t x x uu v u ρ Multply by u and u then add, average. { } 1 p uu = u uu + u u u u u p v u u u + + + u t xk x ρ x x k (Page 61 of Matheu & Scott) Now prmed quanttes can be taken nsde of unprmed dervatve because they do not depend on x. u u = and we can use the notaton: e.g. uk uuu k xk xk ( ) ( uuu k = fk x x = fk r ) where r = x x uuu k r k p Smlarly u u u uu = u and u = u p x r x x u p = up x r u u = uu xk rk u u = uu x r and puttng t together 5

1 R ( rt, ) = { uuu k uuu k} + pu pu + v R rt t rk ρ r r rk I ( r) ( r) P Smlar to / \ eqn. 6.9 two pont pressure velocty 6.95 on pg. 61 cubc moments correlatons n Matheu & Scott (, ) N.B. We have an equaton for the Reynolds stresses (-pt). uu, but now we have ntroduced the closure problem at a hgher level through the trple moment and pressure-velocty correlaton terms that were not present n the R.A.N.S. equatons. Normally a relatonshp between pressure and velocty s found by takng the dvergence of the momentum equatons to form a Posson equaton, e.g. x 1 p = uu ρ x x for nfnte doman, homogeneous flud. 1 1 p( x) = ρ 4π uu x x x x x and the pressure-velocty correlatons. 1 1 x pu = uuku ρ 4πρ x x k x x Pressure velocty correlaton can be wrtten n terms of a trple velocty correlaton. For weak turbulence the trple correlatons are small. R = v R t r Κ 6

R ( rt, ) 11 1 R ( rt, ) 11 Dffuse n tme. 7