Turbulent Flows. Computational Modelling of Turbulent Flows. Overview. Turbulent Eddies and Scales
|
|
- Brianne Shields
- 5 years ago
- Views:
Transcription
1 School of Mechanical Aerospace and Civil Engineering Turbulen Flows As noed above, using he mehods described in earlier lecures, he Navier-Sokes equaions can be discreized and solved numerically on complex domains. Conens: Compuaional Modelling of Turbulen Flows T. J. Craf George Begg Building, C43 In principle, by refining he grid, and using high-order convecion schemes, numerical errors can generally be reduced o negligible levels. For laminar flows, i is usually feasible o obain grid-independen soluions in he above manner. However, many flows of engineering ineres are urbulen, and hese inroduce a number of furher complicaions. Inroducion o CFD Numerical soluion of equaions Reading: Finie difference mehods J. Ferziger, M. Peric, Compuaional Mehods for Fluid Finie volume mehods Dynamics Pressure-velociy coupling H.K. Verseeg, W. Malalasekara, An Inroducion o Compuaional Fluid Dynamics: The Finie Volume Mehod Solving ses of linear equaions nseady problems S.V. Paankar, Numerical Hea Transfer and Fluid Flow Turbulence and oher physical modellingnoes: hp://cfd.mace.mancheser.ac.uk/mcfd Body-fied coordinae sysems - People - T. Craf - Online Teaching Maerial Turbulen flows are sill described by he Navier-Sokes equaions. However, we need o consider some of he differences in behaviour one finds in laminar and urbulen flows. - p. 3 Overview Turbulen Eddies and Scales In earlier lecures we have examined schemes for discreizing, and subsequenly solving, generic fluid flow equaions. We have noed ha he same mehods can be applied o oher ranspor equaions when considering emperaure or mass fracion ranspor, for example. Prof. Osborne Reynolds demonsraed he ransiion from laminar o urbulen pipe flow a Mancheser in 883. When examining finie volume schemes we saw how he full Navier-Sokes equaions could be discreized, and how a pressure-velociy coupling scheme could be applied o obain he pressure field. In his lecure we consider how we migh use he mehods already me o model urbulen flows which are ofen of ineres in engineering problems. One quesion o be addressed is wheher we could simply apply he mehods already me, or wheher furher modelling or approximaion needs o be inroduced. As migh be expeced, we will examine a number of alernaive reamens, wih he aim of oulining when, and in wha circumsances, hey migh be used. He idenified he parameer D/ν as a dimensionless group ha could be used o characerize wheher he flow would be laminar or urbulen. Flow visualizaion of a urbulen boundary layer shows a range of eddy sizes, mixing he fluid. - p. 2 - p. 4
2 In a finie volume scheme (or any oher numerical mehod), he grid mus be fine enough o resolve all he deails of he flow. In a urbulen flow, where here is a large range of urbulen eddy sizes, he compuaional mesh would have o be significanly finer han he smalles of hese eddies. In a flow round an airfoil secion, for insance, he oal domain size migh be of he order of several meres. However, a a high Reynolds number, he smalles eddies may be less han a millimere in size. Since he urbulen eddies are hree-dimensional, a direc soluion of he Navier-Sokes equaions for he above flow would require a full 3-dimensional, ime-dependen, soluion on a grid wih a leas several housand nodes in each of he hree coordinae direcions. These requiremens are well beyond he capabiliies of any curren compuers (and will, mos likely, say so for a long ime). Hence, a number of differen approaches for compuing urbulen flows have been developed. Large Eddy Simulaion (LES) This sraegy also employs a 3-D, ime-dependen, soluion of he Navier-Sokes equaions, bu wih a relaively coarse grid and ime sep, so he fine-scale eddies are no resolved. A model hus has o be inroduced o accoun for he effec of he small-scale eddies. In mos cases his consiss of effecively increasing he viscosiy by an amoun ermed he "sub-grid-scale viscosiy" (which may vary in space and ime). This sub-grid-scale model accouns for he effecs of he small, energy-dissipaing eddies, which are oo small for he compuaional grid o resolve. Mos of he sandard finie volume soluion echniques examined earlier can usually be adoped for LES, alhough a few adapaions are someimes necessary. LES schemes are now beginning o be used a research level for a variey of complex engineering flow. The compuaional ime is sill unaccepably long, however, for mos rouine indusrial applicaions. - p. 5 - p. 7 Direc Numerical Simulaion (DNS) For very low Reynolds number (and simple geomerical) flows, i is possible o simply solve he Navier-Sokes equaions numerically, fully resolving all he eddy sizes and imescales. DNS can provide highly deailed flow daa, bu requires supercompuing class machines. Reynolds Averaging In mos engineering siuaions i is he average velociy ad pressure ha are of ineres, and he fine deails of all he urbulen eddies are no required. As originally proposed by Reynolds, he velociy and pressure fields can be spli ino a mean (or average) and a flucuaing par: e i = i + u i ep = P + p () () = + u() u Since he range of urbulen eddy sizes grows rapidly as he Reynolds number increases, compuaional requiremens soon become prohibiive, as noed above. High order numerical schemes are also normally required. The underlying second-order accuracy of mos finie volume schemes is ofen no sufficien, and specral mehods are frequenly employed. These generally resric one o using simple geomeries. Whils useful for fundamenal flow sudies, DNS is no herefore a ool for rouine engineering calculaions. - p In many cases he flow field may be seady on average, in which case he decomposiion of variables can be easily done by defining 0 Z T i = lim e i d (2) T T 0 - p. 8
3 In flows where here is a large-scale unseadiness a slighly differen definiion is needed. For example, in periodic flows a phase-averaging may be inroduced: NX i () = lim e i ( + NT) (3) N N In ohers, an ensemble-averaging may be used: NX i () = lim e i () (4) N N In such cases he mean quaniies are sill funcions of ime, bu hey only represen he larges unseady moions, and none of he fine-scale urbulence is included in hem. Since only he mean flow is resolved in he calculaions, symmeries (and/or 2-dimensionaliy or seadiness) can be aken advanage of o reduce compuaional requiremens. This approach is herefore he mos widely used in engineering applicaions. T - p. 9 Eddy Viscosiy Models The Reynolds averaged Navier-Sokes equaions can be wrien as i + i j = P + ρ x i ν «i u i u j From hese, i appears ha he u i u j erms simply ac as addiional sresses on op of he viscous sresses ν i /. A simple way o model he Reynolds sresses is hus by inroducing a urbulen or eddy viscosiy, ν, and modelling u i u j in erms of ν and he mean velociy gradiens: i u i u j = ν + «j + 2/3kδ ij (8) x i where δ ij is a special ensor, known as he Kronecker dela, which is defined as δ ij = ( if i = j 0 if i j The quaniy k is he urbulen kineic energy, defined as k = /2(u 2 + u2 2 + u2 3 ) (0) (7) (9) - p. The Reynolds Averaged Navier-Sokes Equaions Having decomposed velociies, pressure, ec. ino mean and flucuaing pars, he nex sep is o subsiue hese ino he Navier-Sokes equaions and average he resuling expressions in order o arrive a a se of equaions governing he behaviour of he mean flow field. Deails of his derivaion are no given here, bu he resul is a se of equaions raher similar o he original Navier-Sokes sysem: Coninuiy: Momenum: i = 0 x i (5) i + ( i j ) = P + ν «i u i u j ρ x i (6) However, he Reynolds sresses, u i u j, are no known (and canno be calculaed, since we are no now resolving he urbulen eddies). Models mus herefore be inroduced o approximae he Reynolds sresses. The subjec of urbulence modelling as a whole is beyond his course. Here we can only noe a few aspecs of some of he more common approaches found in many CFD codes. Wih his model for u i u j, he momenum equaion can simply be wrien as i + i j = P + ρ x i wih he effecive viscosiy ν eff = ν + ν. i ν eff + ««j x i Since he above momenum equaion is essenially he same form as ha already considered, he same numerical procedures already examined can be employed o obain a se of discreized equaions for he mean velociy a each node. The main difference is ha in he diffusive fluxes he viscosiy ha has o be compued a he cell faces is now he effecive viscosiy, ν eff. If he above eddy-viscosiy approach is adoped, he remaining modelling problem is how o prescribe he eddy-viscosiy, ν. Noe ha whils ν is a fluid propery, ν is no i is simply a funcion of he local flow behaviour. () A is simples level, he modelling of he urbulen sresses is done by relaing u i u j o he mean velociy gradiens and some prescribed lenghscale, derived from he geomery of he flow. More complex models solve addiional ranspor equaions for quaniies such as he urbulence energy. - p. 0 - p. 2
4 Zero-Equaion Models In he simple mixing-lengh, or one-equaion model, a mixing-lengh l m is prescribed algebraically, and ν is obained from or, in a more general formulaion: ν = l 2 m ν = lm 2 y (2) s i + j x i l m «2 (3) Two-Equaion Models One problem wih he above models is ha he urbulence lenghscale has o be prescribed, which is ofen difficul o do in a reliable manner for complex flows and geomeries. In a wo-equaion model, a second variable (ofen ε) is solved for: Dk D = P k ε +» (ν + ν /σ k ) k Dε D = c ε ε and he urbulen viscosiy is modelled as k P k c ε2 ε 2 k +» (ν + ν /σ ε ) ε (6) (7) ν = c µk 2 /ε (8) The lenghscale l m is aken o vary linearly wih disance from he wall, bu wih a damping erm o accoun for viscous effecs in he near-wall viscous sublayer. y The model consans c ε, c ε2, c µ, and σ ε are uned for simple equilibrium shear flows. A fairly sandard k-ε model akes c ε =.44, c ε2 =.92, c µ = 0.09, and σ ε =.3. In he near-wall region, where viscous effecs are imporan, a damping funcion is inroduced ino c µ, and addiional source erms are someimes included in he ε equaion. Two-equaion models hus avoid he need o prescribe a lenghscale, and are he mos widely used class of model in ypical indusrial CFD simulaions. - p. 3 - p. 5 One-Equaion Models In a one-equaion model, a ranspor equaion is solved for he urbulen kineic energy k: Dk D = P k ε +» (ν + ν /σ k ) k where P k = u i u j i / and he dissipaion rae ε is modelled as ε = k 3/2 /l ε. The urbulen viscosiy is hen aken as (4) ν = c µk /2 l µ (5) More Advanced Models There are a number of more advanced modelling schemes. For example, full sress ranspor models solve modelled ranspor equaions for each individual sress componen: Du i u j D = P ij + φ ij ε ij + diff ij (9) For a 3-dimensional flow his requires 6 ranspor equaions o be solved for he u i u j ensor (and usually a furher one for he dissipaion rae, ε). Such schemes are herefore more compuaionally expensive han wo-equaion models, bu are poenially more reliable. Boh lenghscales l µ and l ε are ypically prescribed as increasing linearly wih disance from he wall alhough hey may have differen viscous damping erms associaed wih hem o accoun for very near-wall viscous effecs. l y Noe ha he addiional ranspor equaions for urbulence quaniies are of he same generic form as hose already considered when discreizing oher equaions. Broadly he same mehods already covered can hus be applied o discreize and solve hem. In he ypical segregaed finie volume solver srucure examined in earlier lecures, one hus simply inroduces he addiional sep(s) wihin he main ieraion loop of discreizing each urbulence variable ranspor equaion and solving he resuling se of equaions o updae each variable. - p. 4 - p. 6
5 Near-Wall Treamens Close o a wall viscous effecs become imporan, and damping funcions have o be inroduced ino some of he model coefficiens and funcions. There are also usually large gradiens of velociy, urbulence energy (and ε), and a very fine grid mus herefore be employed in order o resolve hese. In he fully urbulen region of a simple zero pressure gradien boundary layer, he velociy profile is known o saisfy he log-law: + = κ log(ey+ ) (2) where he non-dimensional velociy, +, and wall disance, y +, are defined as + = /(τ w/ρ) /2 and y + = y(τ w/ρ) /2 /ν (22) The high aspec raio cells can also affec he numerical soluion sabiliy. In large 3-D calculaions he need o use very fine near-wall grids can make compuaions exremely expensive. If his is assumed o hold a node P, we obain p (τ w/ρ) /2 = κ log(ey p(τ w /ρ) /2 /ν) (23) In an ieraive soluion procedure, equaion (23) can be solved o obain an esimae for he wall shear sress, τ w, which is hen used in updaing he velociy, as described above. Such an approach removes he need o fully resolve he near-wall viscous layer of he flow. Wall funcions are a mehod of reaing he near-wall region of he flow in an approximae manner, wihou fully resolving he hin viscosiy-affeced flow region. - p. 7 However, i does assume ha all near-wall flows saisfy he log-law, which is no rue! There are more advanced wall funcion approaches, which are more widely applicable, bu hese will no be covered here. - p. 9 Wall Funcions The idea of a wall funcion approach is o place he firs compuaional node ouside he viscous sublayer, and make suiable assumpions abou how he near-wall velociy profile behaves, in order o obain he wall shear sress. y n As seen earlier, he momenum equaion in he near-wall cell can be discreized ino an equaion of he form a p p = a e e + a w w + a n n + a s s + S u (20) where a s arises from he shear sress a he souh face of he cell (ie. he wall shear sress). x P y v However, since we are no resolving he near-wall sublayer, esimaing he velociy gradien from a linear variaion of, as we do a oher posiions, will no give an accurae approximaion of he gradien, and hence he shear sress, a he wall. If, insead, we could obain he wall shear sress from an assumed velociy profile, we could se a s o zero and hen add he wall shear sress conribuion, τ w( x), direcly ino he source erm. - p. 8
Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationUnsteady Flow Problems
School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationDiffusion & Viscosity: Navier-Stokes Equation
4/5/018 Diffusion & Viscosiy: Navier-Sokes Equaion 1 4/5/018 Diffusion Equaion Imagine a quaniy C(x,) represening a local propery in a fluid, eg. - hermal energy densiy - concenraion of a polluan - densiy
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationSummary of shear rate kinematics (part 1)
InroToMaFuncions.pdf 4 CM465 To proceed o beer-designed consiuive equaions, we need o know more abou maerial behavior, i.e. we need more maerial funcions o predic, and we need measuremens of hese maerial
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationTHE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES
Kragujevac J. Sci. 3 () 7-4. UDC 53.5:536. 4 THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Hazem A. Aia Dep. of Mahemaics, College of Science,King Saud Universiy
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationFlow-Induced Vibration Analysis of Supported Pipes with a Crack
Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationComputation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM
Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationFailure of the work-hamiltonian connection for free energy calculations. Abstract
Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationNew effective moduli of isotropic viscoelastic composites. Part I. Theoretical justification
IOP Conference Series: Maerials Science and Engineering PAPE OPEN ACCESS New effecive moduli of isoropic viscoelasic composies. Par I. Theoreical jusificaion To cie his aricle: A A Sveashkov and A A akurov
More informationThis document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC
his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More informationFinite Element Analysis of Structures
KAIT OE5 Finie Elemen Analysis of rucures Mid-erm Exam, Fall 9 (p) m. As shown in Fig., we model a russ srucure of uniform area (lengh, Area Am ) subjeced o a uniform body force ( f B e x N / m ) using
More informationLAB # 2 - Equilibrium (static)
AB # - Equilibrium (saic) Inroducion Isaac Newon's conribuion o physics was o recognize ha despie he seeming compleiy of he Unierse, he moion of is pars is guided by surprisingly simple aws. Newon's inspiraion
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationAdvanced FDTD Algorithms
EE 5303 Elecromagneic Analsis Using Finie Difference Time Domain Lecure #5 Advanced FDTD Algorihms Lecure 5 These noes ma conain coprighed maerial obained under fair use rules. Disribuion of hese maerials
More informationLinear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation
Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More information1 Review of Zero-Sum Games
COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationIB Physics Kinematics Worksheet
IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?
More informationNumerical investigation of Ranque-Hilsch energy separation effect A.S. Noskov 1,a, V.N. Alekhin 1,b, A.V. Khait 1,a
Applied Mechanics and Maerials Online: 2013-01-11 ISSN: 1662-7482, Vol. 281, pp 355-358 doi:10.4028/www.scienific.ne/amm.281.355 2013 Trans Tech Publicaions, Swizerland Numerical invesigaion of Ranque-Hilsch
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationMaintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011
Mainenance Models Prof Rober C Leachman IEOR 3, Mehods of Manufacuring Improvemen Spring, Inroducion The mainenance of complex equipmen ofen accouns for a large porion of he coss associaed wih ha equipmen
More informationExplaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015
Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become
More informationEE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?
EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of
More informationRaNS.1 RaNS + PDE Turbulence Closure Models
RaNS.1 RaNS + PDE Turbulence Closure Models Time-averaged Navier-Sokes PDE closure models one-equaion: k wih l specified υ, Re ranspor wo-equaion: k, ε ranspor k,ω ranspor Turbulen kineic energy ranspor
More informationLecture 4 Notes (Little s Theorem)
Lecure 4 Noes (Lile s Theorem) This lecure concerns one of he mos imporan (and simples) heorems in Queuing Theory, Lile s Theorem. More informaion can be found in he course book, Bersekas & Gallagher,
More informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationNumerical Simulation of the Overall Flow Field for Underwater Vehicle with Pump Jet Thruster
Available online a www.sciencedirec.com Procedia Engineering 31 (2012) 769 774 Inernaional Conference on Advances in Compuaional Modeling and Simulaion Numerical Simulaion of he Overall Flow Field for
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationReading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.
PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationTechnical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.
Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear
More informationMechanical Fatigue and Load-Induced Aging of Loudspeaker Suspension. Wolfgang Klippel,
Mechanical Faigue and Load-Induced Aging of Loudspeaker Suspension Wolfgang Klippel, Insiue of Acousics and Speech Communicaion Dresden Universiy of Technology presened a he ALMA Symposium 2012, Las Vegas
More information1 Evaluating Chromatograms
3 1 Evaluaing Chromaograms Hans-Joachim Kuss and Daniel Sauffer Chromaography is, in principle, a diluion process. In HPLC analysis, on dissolving he subsances o be analyzed in an eluen and hen injecing
More informationMorning Time: 1 hour 30 minutes Additional materials (enclosed):
ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationTurbulence in Fluids. Plumes and Thermals. Benoit Cushman-Roisin Thayer School of Engineering Dartmouth College
Turbulence in Fluids Plumes and Thermals enoi Cushman-Roisin Thayer School of Engineering Darmouh College Why do hese srucures behave he way hey do? How much mixing do hey accomplish? 1 Plumes Plumes are
More informationUnsteady Mass- Transfer Models
See T&K Chaper 9 Unseady Mass- Transfer Models ChEn 6603 Wednesday, April 4, Ouline Conex for he discussion Soluion for ransien binary diffusion wih consan c, N. Soluion for mulicomponen diffusion wih
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationFrom Particles to Rigid Bodies
Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and
More informationSterilization D Values
Seriliaion D Values Seriliaion by seam consis of he simple observaion ha baceria die over ime during exposure o hea. They do no all live for a finie period of hea exposure and hen suddenly die a once,
More informationA Shooting Method for A Node Generation Algorithm
A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan
More informationSection 7.4 Modeling Changing Amplitude and Midline
488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More informationnon -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.
LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationBU Macro BU Macro Fall 2008, Lecture 4
Dynamic Programming BU Macro 2008 Lecure 4 1 Ouline 1. Cerainy opimizaion problem used o illusrae: a. Resricions on exogenous variables b. Value funcion c. Policy funcion d. The Bellman equaion and an
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationLES at high Re numbers: a near-wall eddy-viscosity formulation
a high Re numbers: a near-wall eddy-viscosiy formulaion Georgi Kalizin, Jeremy A. Templeon & Gorazd Medic Mechanical ngineering Deparmen, Sanford Universiy, Sanford, CA 943, USA Absrac A near-wall eddy-viscosiy
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationCSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering
CSE 3802 / ECE 3431 Numerical Mehods in Scienific Compuaion Jinbo Bi Deparmen of Compuer Science & Engineering hp://www.engr.uconn.edu/~jinbo 1 Ph.D in Mahemaics The Insrucor Previous professional experience:
More informationA First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18
A Firs ourse on Kineics and Reacion Engineering lass 19 on Uni 18 Par I - hemical Reacions Par II - hemical Reacion Kineics Where We re Going Par III - hemical Reacion Engineering A. Ideal Reacors B. Perfecly
More informationNUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD FACING STEP WITH A FOUR PARAMETER LOGARITHMIC TURBULENCE MODEL
h European Conference on Compuaional Mechanics (ECCM ) 7h European Conference on Compuaional Fluid Dynamics (ECFD 7) 11 June, Glasgow, UK NUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD
More informationVerification of a CFD benchmark solution of transient low Mach number flows with Richardson extrapolation procedure 1
Verificaion of a CFD benchmark soluion of ransien low Mach number flows wih Richardson exrapolaion procedure S. Beneboula, S. Gounand, A. Beccanini and E. Suder DEN/DANS/DMS/STMF Commissaria à l Energie
More informationAnd the solution to the PDE problem must be of the form Π 1
5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release
More informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationE β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.
Noes, M. Krause.. Problem Se 9: Exercise on FTPL Same model as in paper and lecure, only ha one-period govenmen bonds are replaced by consols, which are bonds ha pay one dollar forever. I has curren marke
More informationApplications of the Basic Equations Chapter 3. Paul A. Ullrich
Applicaions of he Basic Equaions Chaper 3 Paul A. Ullrich paullrich@ucdavis.edu Par 1: Naural Coordinaes Naural Coordinaes Quesion: Why do we need anoher coordinae sysem? Our goal is o simplify he equaions
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More informationModule 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II
Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2
More informationDemodulation of Digitally Modulated Signals
Addiional maerial for TSKS1 Digial Communicaion and TSKS2 Telecommunicaion Demodulaion of Digially Modulaed Signals Mikael Olofsson Insiuionen för sysemeknik Linköpings universie, 581 83 Linköping November
More informationEE650R: Reliability Physics of Nanoelectronic Devices Lecture 9:
EE65R: Reliabiliy Physics of anoelecronic Devices Lecure 9: Feaures of Time-Dependen BTI Degradaion Dae: Sep. 9, 6 Classnoe Lufe Siddique Review Animesh Daa 9. Background/Review: BTI is observed when he
More information