Turbulence and its Modelling
|
|
- Deirdre Mathews
- 5 years ago
- Views:
Transcription
1 School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows become turbulent. Turbulence and ts Modellng T. J. Craft George Begg Buldng, C41 Contents: Readng: Naver-Stokes equatons F.M. Whte, Flud Mechancs Invscd flows J. Matheu, J. Scott, An Introducton to Turbulent Flow Boundary layers P.A. Lbby, Introducton to Turbulence Transton, Reynolds averagng P. Bernard, J. Wallace, Turbulent Flow: Analyss Measurement & Predcton Mxng-length models of turbulence S.B. Pope, Turbulent Flows Turbulent knetc energy equaton D. Wlcox, Turbulence Modellng for CFD One- and Two-equaton models Notes: Flow management - People - T. Craft - Onlne Teachng Materal We have seen how the Reynolds averagng approach can be used to descrbe turbulent flows. ~ U ~ U(t) = + U u(t) U t 0 0 We saw the Reynolds Averaged Naver-Stokes (RANS) equatons, and the appearance n these of the Reynolds stresses, u u j. Turbulence and ts Modellng 2010/11 2 / 27 t u 0 t Here, we begn to consder how one mght use the Naver-Stokes or the RANS equatons to compute engneerng-related turbulent flows. Ths wll lead to the concluson that we need to ntroduce mathematcal models (approxmatons) for the Reynolds stresses. Before examnng how such models mght be devsed, we wll frst consder some of the physcal features of turbulent flows and the behavour of the Reynolds stresses n certan, farly smple, flows. We can then examne how some of these features are (or ar not) accounted for n some of the more wdely-used, but stll relatvely smple, turbulence models. Turbulence and ts Modellng 2010/11 3 / 27 The Energy Cascade Turbulence s characterzed by a wde range of eddy szes. Flow nstabltes develop across a range of wavenumbers, leadng to large eddes breakng down nto smaller ones, these breakng nto stll smaller ones, etc... The turbulent knetc energy, k = 0.5(u 2 + v 2 + w 2 ), s manly transferred from the mean flow va the larger eddes n the flow. As these eddes break up, knetc energy gets transferred nto smaller eddes and, as these break up, nto even smaller ones. Thus, one gets a cascade of turbulent knetc energy transferred down the scale spectrum, from the largest eddes to the smallest ones. At the smallest scales (whch wll also be examned later) the knetc energy s dsspated nto heat. Turbulence and ts Modellng 2010/11 4 / 27
2 The turbulence energy spectrum shows how much knetc energy s contaned n the dfferent scales of moton. 0.5u u = E(k,t)dk (1) 0 k n the fgure represents the wavenumber of the dfferent scales of moton. Large k corresponds to small eddes, and small k to large eddes. The lowest wavenumbers (largest eddes) are typcally related to the flow geometry. As the Reynolds number ncreases, the sze of the smallest scales present (largest wavenumbers) decreases. Thus, as the Reynolds number ncreases there s a wder range of turbulent eddy szes n the flow. Turbulent Eddy Szes The largest eddes n the flow generally scale wth the flow geometry, eg.: Wth ppe dameter or channel wdth. Wth boundary layer thckness. At the smallest scales the turbulence energy s dsspated to heat by the acton of vscosty. Dmensonal analyss suggests that the sze of these smallest eddes wll be related to (ν 3 /ε) 1/4 where ν s the knematc vscosty of the flud and ε s the dsspaton rate of the turbulent knetc energy. Note that as ν decreases (as the Reynolds number ncreases), the sze of the smallest eddes decreases. Thus, as the Reynolds number ncreases one expects to get a wder range of turbulent eddy szes. The large eddes can be very ansotropc. At the smaller scales the turbulence s generally reckoned to become more sotropc. Turbulence and ts Modellng 2010/11 5 / 27 Turbulence and ts Modellng 2010/11 6 / 27 Computaton of Turbulent Flow The Naver-Stokes equatons do stll descrbe a turbulent flow ncludng all the turbulent eddy detals. In prncple, one could solve these drectly (numercally) and obtan full detals of the turbulent flow feld (averagng the results to obtan mean veloctes, etc., f desred). If we consder a flow wth velocty and length scales U and L, we would need a computatonal doman of sze L. As noted above, the smallest eddy sze s of order (ν 3 /ε) 1/4. If we assume ε can be scaled wth U 3 /L (on dmensonal grounds), then the smallest eddy szes are of order (ν 3 L/U 3 ) 1/4. L U The number of computatonal nodes needed to cover a length L, whlst resolvng the smallest eddes s therefore of order L/(ν 3 L/U 3 ) 1/4 = (UL/ν) 3/4 = Re 3/4. Snce turbulent flows are 3-D, we need ths many ponts n each drecton, so the total number of computatonal nodes scales as Re 9/4. However, a consderaton of the range of scales nvolved shows ths s not actually possble, wth current technology, for most flows of engneerng nterest. The grd requrements thus ncrease rapdly as the Reynolds number of the flows goes up. Turbulence and ts Modellng 2010/11 7 / 27 Turbulence and ts Modellng 2010/11 8 / 27
3 Drect Numercal Smulaton (DNS) For rather low Reynolds numbers the Naver-Stokes equatons can be solved numercally, provdng full resoluton of a turbulent flow. The numercal soluton methods requre more accuracy than the standard fnte volume schemes studed n Modellng & Smulaton or CFD-1, and ths often restrcts ther use to smple geometres. Fully developed turbulent channel flows at Reynolds numbers (based on bulk velocty and channel half-wdth) of around 2500 were performed n the 1980 s. More recent channel flow smulatons at Re have been reported. These were performed on a 2048 processor machne, requred processor hours (the entre machne for around 4 months), and produced around 25TB of raw data. In engneerng applcatons, much hgher Reynolds numbers than ths are often encountered. Large gas ppelnes, for example, can be operatng at Reynolds numbers of several mllon. Turbulence and ts Modellng 2010/11 9 / 27 DNS can gve very useful, detaled, data for fundamental turbulence studes, and can be used to gude turbulence model development. Cross-stream flow & streamwse fluctuatons Vortcty sosurfaces However, t clearly s not a tool for everyday engneerng use. Turbulence and ts Modellng 2010/11 10 / 27 RANS Modellng There are modellng approaches such as Large Eddy Smulaton (LES), where only the large-scale turbulence structures are resolved, and approxmatons ntroduced for the effects of the smaller eddes. These are computatonally cheaper than DNS, but often stll too expensve for routne use n an engneerng envronment. Most turbulent flow smulatons thus employ the RANS equatons, solvng for the mean, or average, flow feld. For ncompressble flow wthout body forces, the governng equatons are then: DU = 1 P + ( ν U ) u Dt ρ x x j x u j (2) j U = 0 (3) x (note the use of tensor summaton conventon). In 3-D we have 3 momentum equatons and contnuty. However, the unknowns are 3 velocty components, pressure, and the 6 Reynolds stresses. To close the system we thus need to ntroduce models (approxmatons) for the Reynolds stresses, uuj. There are many such turbulence models avalable, havng dfferent levels of complexty, and provdng dfferng levels of accuracy n varous flows. To understand whch model to choose for a partcular applcaton t s necessary to have some apprecaton of how they were devsed, and hence some knowledge of ther strengths and/or weaknesses. In order to assess a model s accuracy, we need to know what behavour would be found n a real flud flow. Thus, before consderng how models may be devsed, we examne some features of a few smple turbulent flow stuatons, and the behavour exhbted by the turbulence statstcs. Turbulence and ts Modellng 2010/11 11 / 27 Turbulence and ts Modellng 2010/11 12 / 27
4 Shear Stresses Across a Boundary Layer Wth the usual boundary layer approxmatons the mean U-momentum equaton may be wrtten U 2 x + UV = 1 P ρ x + ( ν U ) uv (4) In a turbulent boundary layer, the turbulent shear stress uv s thus the most mportant component of the Reynolds stress tensor. If we assume that convecton s neglgble, then the U-momentum equaton becomes P x = τ (5) where τ s the total shear stress, gven by the sum of the molecular and turbulent shear stresses: τ = µ U ρuv (6) }{{}}{{} Turbulent Molecular Turbulence and ts Modellng 2010/11 13 / 27 y U o x In a zero pressure gradent boundary layer the total shear stress s thus constant across the layer, τ/ = 0. At the wall, uv = 0 and µ U/ = τw where τw s the wall shear stress. We thus obtan µ U ρuv = Constant = τw (7) Immedately adjacent to the wall the turbulent shear stress s neglgble (snce uv = 0 at the wall), and the molecular stress s the domnant. Outsde the thn vscous sublayer, n the fully turbulent part of the boundary layer, molecular effects are neglgble and the turbulent shear stress s then essentally equal to the total stress. In the fully turbulent near-wall regon of a zero pressure gradent boundary layer the turbulent shear stress s thus equal n magntude to the wall shear stress: uv = τw/ρ Vscous sublayer Turb. shear stress Mol. shear stress Turbulence and ts Modellng 2010/11 14 / 27 τ w/ρ y Wth a favourable pressure gradent we get τ/ < 0 so the total shear stress decreases wth dstance from the wall. Conversely, n an adverse pressure gradent we get τ/ > 0 and the total shear stress ncreases wth dstance from the wall (although as separaton s approached the above approxmatons may not be vald). In practce, the above approxmatons do not hold n the outer boundary layer regons, but are not too far from realty n the near-wall regons. Near-Wall Non-Dmensonalzaton It s often convenent to non-dmensonalze near-wall flow quanttes usng wall parameters, rather than the global velocty/length scales. (τw/ρ) 1/2 s a velocty scale, and ν/(τw/ρ) 1/2 a length scale. The velocty and dstance from the wall are then non-dmensonalzed as U + = U/(τw/ρ) 1/2 and y + = y(τw/ρ) 1/2 /ν In the fully turbulent near-wall regon of a zero pressure gradent boundary layer the velocty profle s found to ft a log-law: U + = 1 κ log(ey+ ) (8) Zero pressure gradent Favourable pressure gradent wth constants κ 0.41 and E 9. Turbulence and ts Modellng 2010/11 15 / 27 Turbulence and ts Modellng 2010/11 16 / 27
5 Plottng the shear stresses aganst y + we can see the vscous sublayer typcally extends to around y (although the exact value s rather flow-dependent). Zero pressure gradent Favourable pressure gradent Turbulent Knetc Energy Transport Equaton In an earler lecture we examned the mean knetc energy transport equaton n a turbulent flow. The turbulent knetc energy s defned by k 0.5u 2 0.5(u 2 + v 2 + w 2 ). To derve ts transport equaton, we can use Dk Dt = D ( ) u 2 Du Dt /2 = u Dt A transport equaton for the fluctuatng velocty, u, can be obtaned by subtractng the Reynolds-averaged momentum equaton from the Naver Stokes equaton: Du Dt = DŨ DU Dt Dt Note that here the mean veloctes are used n the convectve contrbutons to the total dervatve D/Dt. e. we nterpret D Dt t + U j x j (9) Turbulence and ts Modellng 2010/11 17 / 27 Turbulence and ts Modellng 2010/11 18 / 27 The dervaton s left as an exercse, but the result s: u u + U t j = 1 p U u xj ρ j + ν 2 u x xj xj 2 ( ) uu j u u j xj (10) Multplyng ths equaton by the fluctuatng velocty u (note the mpled summaton over the ndex ), and averagng, we can arrve at: k t +U k U j = u x u j ν u u ( u j x j x j x j x j 2 u /2+u p/ρ ν k ) (11) x The second term on the rght hand sde of equaton (11) represents the dsspaton of turbulent knetc energy by vscosty at the smallest scales. It s usually denoted as ε. The last term represents dffuson, as a result of vscous and turbulent mxng. The term u u j U / x j must therefore represent the generaton of turbulent knetc energy. It s often denoted by P k and called the producton, or generaton, rate of k. Turbulence and ts Modellng 2010/11 19 / 27 P k appears n both mean and turbulent knetc energy equatons (wth opposte sgns), and s thus nterpreted as the rate at whch knetc energy s lost from the mean flow and transferred to the turbulent eddes. DK Dt = 1 (PU ρ x j)+ j x j ( ν K ) U x u u j j ( ) U 2 + u U u j ν x j x j }{{} P k Dk Dt = u U u j ν u u ( u x j x j x j x j 2 u /2+u p/ρ ν k ) x } {{ } } {{ } P k ε In most crcumstances P k s postve, representng a transfer of knetc energy from the mean flow to the turbulence. However, there are flow condtons under whch P k can be locally negatve n certan regons. Turbulence and ts Modellng 2010/11 20 / 27
6 Boundary Layer Form of the k Equaton Wth the usual approxmatons for a steady boundary layer (or other thn shear flow) the k equaton reduces to (Uk) + (Vk) = Pk ε + ( u2uu/2 pu2/ρ ν k ) (12) x The generaton rate P k now reduces to P k = uv U Note that, as n the boundary layer U-momentum equaton, t s the shear stress component, uv, that s the most mportant element of the Reynolds stress tensor. We can now explan why the peak levels of turbulent knetc energy generally occur near the wall. y U o x As seen earler, n the fully turbulent regon of the boundary layer we can approxmate the turbulent shear stress as a constant (uv τw/ρ), and the mean velocty fts a logarthmc profle, so dfferentatng gves U/ 1/y. In the fully turbulent regon Pk thus decreases wth dstance from the wall as P k 1/y. At the wall tself, P k = 0, snce the turbulent shear stress vanshes. Hence, Pk must have a maxmum somewhere between the wall and the fully turbulent regon e. wthn the vscous sublayer. In an earler lecture we saw ths peak occurs roughly where the molecular and vscous shear stresses are equal. Turbulence and ts Modellng 2010/11 21 / 27 Turbulence and ts Modellng 2010/11 22 / 27 Plane Channel Flow at Re b 5500 Normal Reynolds Stresses Shear stresses Mean knetc energy budget Turbulent knetc energy Turbulent knetc energy budget In a boundary layer, or other shear flow, the 3 normal stresses (u 2, v 2 and w 2 ) are not equal. In general, the streamwse stress s largest, and the wall-normal the smallest. Transport equatons can be derved for each Reynolds stress component, u u j, smlar to that obtaned for k. The result s of the form Du u j = P Dt j ε j + φ j + d j (13) where the generaton rate P j = u u k U j/ x k u ju k U / x k and ε j, d j and φ j represent dsspaton, dffuson and a redstrbuton term. Turbulence and ts Modellng 2010/11 23 / 27 Turbulence and ts Modellng 2010/11 24 / 27
7 In a smple shear flow the generaton terms are gven by P11 = 2uv U P22 = P33 = 0 As expected, snce k = 0.5(u 2 + v 2 + w 2 ), we see that P11 + P22 + P33 = 2P k. P12 = v 2 U Note, however, that all the turbulent knetc energy generaton actually happens n the u 2 component. Other processes then redstrbute ths knetc energy nto v 2 and w 2. Ths explan why the streamwse normal stress s generally larger than the other two components. Turbulence and ts Modellng 2010/11 25 / 27 Wall-Lmtng Stress Behavour To examne the mmedate near-wall behavour of the stresses, we can express the near-wall veloctes as Taylor seres expansons n powers of y (the wall-normal dstance): u = a1y + b1y 2 + c1y v = a2y + b2y 2 + c2y w = a3y + b3y 2 + c3y where the a s, b s, c s etc. are functons of x, z and t. The above expansons ensure that the veloctes vansh at the wall, but they must also satsfy contnuty ( u/ x + v/+ w/ z = 0). Substtutng the expansons for the veloctes nto the contnuty equaton: y a1/ x + y 2 b1/ x a2 + 2y b2 + 3y 2 c y a3/ z + y 2 b3/ z +... = 0 Consderng the O(1) terms leads to a2 = 0. Turbulence and ts Modellng 2010/11 26 / 27 Hence close to the wall the veloctes behave as u y, w y, but v y 2. The Reynolds stresses thus behave as u 2 y 2, v 2 y 4, w 2 y 2 and uv y 3 for small y. Snce k = 0.5(u 2 + v 2 + w 2 ), then k y 2. Plottng the near-wall stresses on log-log axes shows they do closely match these wall-lmtng forms. Turbulence and ts Modellng 2010/11 27 / 27
Turbulent Flow. Turbulent Flow
http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and
More informationIntroduction to Turbulence Modeling
Introducton to Turbulence Modelng Professor Ismal B. Celk West Vrgna nversty Ismal.Celk@mal.wvu.edu CFD Lab. - West Vrgna nversty I-1 Introducton to Turbulence CFD Lab. - West Vrgna nversty I-2 Introducton
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationPublication 2006/01. Transport Equations in Incompressible. Lars Davidson
Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport
More informationFlow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved
Smulaton of nose generaton and propagaton caused by the turbulent flow around bluff bodes Zamotn Krll e-mal: krart@gmal.com, cq: 958886 Summary Accurate predctons of nose generaton and spread n turbulent
More informationLecture 12. Modeling of Turbulent Combustion
Lecture 12. Modelng of Turbulent Combuston X.S. Ba Modelng of TC Content drect numercal smulaton (DNS) Statstcal approach (RANS) Modelng of turbulent non-premxed flames Modelng of turbulent premxed flames
More informationA Numerical Study of Heat Transfer and Fluid Flow past Single Tube
A Numercal Study of Heat ransfer and Flud Flow past Sngle ube ZEINAB SAYED ABDEL-REHIM Mechancal Engneerng Natonal Research Center El-Bohos Street, Dokk, Gza EGYP abdelrehmz@yahoo.com Abstract: - A numercal
More informationTurbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.
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
More informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationHandout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models
Handout: Large Eddy mulaton I 058:68 Turbulent flows G. Constantnescu Introducton to ubgrd-cale (G) Models G tresses should depend on: Local large-scale feld or Past hstory of local flud (va PDE) Not all
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationTurbulent Transport in Single-Phase Flow. Peter Bernard, University of Maryland
Turbulent Transport n Sngle-Phase Flow Peter Bernard, Unversty of Maryland Assume that our goal s to compute mean flow statstcs such as U and One can ether: 1 u where U Pursue DNS (.e. the "honest" approach)
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationIC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model
IC Engne Flow Smulaton usng KIVA code and A Modfed Reynolds Stress Turbulence Model Satpreet Nanda and S.L. Yang Mechancal Engneerng-Engneerng Mechancs Department Mchgan Technologcal Unversty Houghton,
More informationTools for large-eddy simulation
Center for Turbulence Research Proceedngs of the Summer Program 00 117 Tools for large-eddy smulaton By Davd A. Caughey AND Grdhar Jothprasad A computer code has been developed for solvng the ncompressble
More informationTURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March
TURBULENT FLOW A BEGINNER S APPROACH Tony Saad March 2004 http://tsaad.uts.edu - tsaad@uts.edu CONTENTS Introducton Random processes The energy cascade mechansm The Kolmogorov hypotheses The closure problem
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationIntroduction to Turbulence. by Håkan Gustavsson Division of Fluid Mechanics Luleå University of Technology
Introducton to Turbulence by Håkan Gustavsson Dvson of Flud Mechancs Luleå Unversty of Technology Foreword The ntenton wth these pages s to present the student wth the basc theoretcal concepts of turbulence
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More information2) For a two-dimensional steady turbulent flow in Cartesian coordinates (x,y), with mean velocity components (U,V), write
058:68 Turbulent Flows 004 G. Constantnescu HOMEWORKS: Assgnment I - 01/6/04, Due 0/04/04 1) A cubcal box of volume L 3 s flled wth flud n turbulent moton. No source of energy s present, so that the turbulence
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More information1. Governing Equations
1. Governng Equatons 1a. Governng Equatons for Mean Varables The governng equatons descrbe the varaton n space and tme of the zonal, merdonal and vertcal wnd components, densty, temperature, specfc humdty
More informationSimulation of Turbulent Flow Using FEM
Internatonal Journal of Engneerng and Technology Volume 2 No. 8, August, 2012 Smulaton of Turbulent Flow Usng FEM Sabah Tamm College of Computng, AlGhurar Unversty, Duba, Unted Arab Emrates. ABSTRACT An
More informationTurbulent Flows and their Numerical Modeling
REVIEW Lecture 4: Numercal Flud Mechancs Sprng 015 Lecture 5 Fnte Volume on Complex geometres Computaton of convectve fluxes: x y For md-pont rule: Fe ( vn. ) ds fs e e ( ( vn. ) e Se e m e e e ( Se ue
More informationHomogeneous model: Horizontal pipe and horizontal well. Flow loops can't duplicate field conditions. Daniel D. Joseph. April 2001
Homogeneous model of producton of heavy ol through horzontal ppelnes and wells based on the Naver-Stokes equatons n the ppelne or the well and Darcy's law n the reservor Homogeneous model: Danel D. Joseph
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationPsychology 282 Lecture #24 Outline Regression Diagnostics: Outliers
Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.
More informationMoments 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.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationFEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS
C I T E S 009_ Krasnoyarsk 009 FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS A. F. Kurbatsky Insttute
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationTurbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics
Internatonal Research Journal of Engneerng and Technology (IRJET) e-issn: 2395-56 Volume: 4 Issue: 7 July -217 www.ret.net p-issn: 2395-72 Turbulent Flow n Curved Square Duct: Predcton of Flud flow and
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationSignificance of Dirichlet Series Solution for a Boundary Value Problem
IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More informationCalculation of Aerodynamic Characteristics of NACA 2415, 23012, Airfoils Using Computational Fluid Dynamics (CFD)
Calculaton of Aerodynamc Characterstcs of NACA 2415, 23012, 23015 Arfols Usng Computatonal Flud Dynamcs (CFD) Hmanshu Parashar Abstract A method of solvng the flow over arfols of Natonal Advsory Commttee
More informationThe Discretization Process
FMIA F Moukalled L Mangan M Darwsh An Advanced Introducton wth OpenFOAM and Matlab Ths textbook explores both the theoretcal foundaton of the Fnte Volume Method (FVM) and ts applcatons n Computatonal Flud
More informationNumerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method
Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) Numercal Smulaton of Ld-Drven Cavty Flow Usng the Lattce Boltzmann Method M.A. MUSSA, S. ABDULLAH *, C.S. NOR AZWADI
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationA new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems
Home Search Collectons Journals About Contact us My IOPscence A new ntegrated-rbf-based doman-embeddng scheme for solvng flud-flow problems Ths artcle has been downloaded from IOPscence. Please scroll
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationThermodynamics and statistical mechanics in materials modelling II
Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationResearch & Reviews: Journal of Engineering and Technology
Research & Revews: Journal of Engneerng and Technology Case Study to Smulate Convectve Flows and Heat Transfer n Arcondtoned Spaces Hussen JA 1 *, Mazlan AW 1 and Hasanen MH 2 1 Department of Mechancal
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationBindel, Spring 2014 Applications of Parallel Computers (CS 5220)
Bndel, Sprng 2014 Applcatons of Parallel Computers (CS 5220) 1 Dervng SPH The Naver-Stokes equatons wth gravty are ρa = p + µ 2 v + ρg. The acceleraton s the materal dervatve of velocty, and we usually
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationA NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW
A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW We Zhang and Qngyan Chen Buldng Technology Program Massachusetts Insttute of Technology 77 Mass. Ave., Cambrdge, MA
More informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationConvective Heat Transfer Past. Small Cylindrical Bodies. By Michele S. Titcombe and Michael J. Ward 1. Abstract
Convectve Heat Transfer Past Small Cylndrcal Bodes By Mchele S. Ttcombe and Mchael J. Ward 1 Abstract Convectve heat transfer from an array of small, cylndrcal bodes of arbtrary shape n an unbounded, two-dmensonal
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationEntropy generation in a chemical reaction
Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationSolution of the Navier-Stokes Equations
Numercal Flud Mechancs Fall 2011 Lecture 25 REVIEW Lecture 24: Soluton of the Naver-Stokes Equatons Dscretzaton of the convectve and vscous terms Dscretzaton of the pressure term Conservaton prncples Momentum
More informationBuckingham s pi-theorem
TMA495 Mathematcal modellng 2004 Buckngham s p-theorem Harald Hanche-Olsen hanche@math.ntnu.no Theory Ths note s about physcal quanttes R,...,R n. We lke to measure them n a consstent system of unts, such
More informationFORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO
ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO Mohammad Rahnama*, Seyed-Mad Hasheman*, Mousa Farhad**
More information