Centro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal
|
|
- Daniela Rodgers
- 5 years ago
- Views:
Transcription
1 DEVELOPING CLOSURES FOR TURBULENT FLOW OF VISCOELASTIC FENE-P FLUIDS F. T. Pinho Centro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal C. F. Li Dep. Energy, Environmental and Chemical Engineering, Washington University of St. Louis, St Louis, MO, USA B. A. Younis Dep. Civil and Environmental Engineering, University of California, Davis, USA R. Sureshkumar Dep. Energy, Environmental and Chemical Engineering, Washington University of St. Louis, St Louis, MO, USA Acknowledgments:FCT (POCI/EQU/56342/24), Gulbenkian Foundation Conference on Whither rheology? 2 nd -4 th April 27 Lake Vyrnwy Hotel, Mid Wales, UK
2 Summary Relevance: some facts about drag reduction Governing equations for FENE-P in RANS form Reynolds decomposition of FENE-P equations Closure needs and analysis of DNS case (LDR) A simplified closure with a priori testing of DNS data Some concerns regarding limiting cases Some preliminary results Conclusions 2
3 Relevance: drag reduction in turbulent pipe flow 1-1 f Blasius eq. 6 u % PAA Re=429.25% CMC Re= 166.3% CMC Re=43.9/.9% CMC/XG Re= % PAA.2% PAA.25% CMC.2% XG.9%/.9% CMC/XG Virks MDRA asymptote Re Reduction of shear Reynolds stress (DR) - Increase of normal streamwise Reynolds stress - Dampening of normal radial and tangential Reynolds stress Deficit of Reynolds stress y + Escudier et al JNNFM (1999) 3
4 Time-average governing equations: turbulent flow & FENE-P Continuity:!U i!x i = Momentum balance:! "U i "t ^- instantaneous Overbar or capital letter- time-average or small letter- fluctuations "U +!U i k = # "p " +$ 2 U i s # " (!u i u "x k "x i "x k "x k "x k )+ "% ik, p k "x k Rheological constitutive equation: FENE-P! ij = 2" s S ij +! ij, p ˆ! ij, p = " p & # f ( Ĉ kk )Ĉij $ f ( L )% ( ij ) $ f ( Ĉ kk )Ĉij +! "Ĉij "Ĉij +Ûk # "t "x Ĉ "Ûi "Û # jk k "x Ĉik j & ) = f (L) ij % k "x k ( # % $!Ĉij!t +Ûk!Ĉij!x k " Ĉ jk!ûi!û "!x Ĉik j k!x k & ) ( = Ĉij = " ˆ ij, p + p 4
5 DNS case: LDR DNS, DR=18% (LDR) Fully-developed channel flow u 2,y 1,x 2h We! = 25, Re! = 395 " =.9, L 2 = 9 We! = "u 2! Re #! = hu! " 5
6 Reynolds decomposition of conformation tensor ˆB = B + b where b = Time average polymeric stress Function: f ( C kk ) = L2! 3 L 2! C kk! ij, p = " p & # f ( C kk )C ij $ f ( L)% ij ( ) + " p # f ( C + c kk kk )c ij Time average conformation tensor equation f c ij Neglected (see slide) " + #c %! C ij+! ij #u uk $ c i #u - kj + c j ik #x,- k & #x k #x k (. )/ = $ +, f ( C kk )C ij $ f ( L)1 ij + f ( C kk + c kk )c ij. /! "c ij $ "u C ij + uk # c i "u j kj + c ik "x & k % "x k "x k ) ( = # ij, p + p CT ij NLT ij 6
7 Simplifying assumptions: justification from DNS times smaller c 12 /[L 2 -(C kk +c kk )] C 12 /[L 2 -(C kk )] y + f ( C kk )C 12 >> f c 12 = f ( C kk + c kk )c 12 = ( L 2! 3) c 12 ( ) L 2! C kk + c kk 7
8 Polymer stress 15 1 (a) xx component CT NLT M D S= -! p xx xy 1 5 (b) xy component CT NLT M D S = -! p (c) yy component y + CT NLT M D S= -! p y + 2 yy CT(J) aver y + Housiadas et al (25) Li et al (26) JNNFM 8
9 Modeling requirements! $ ij, p %c & = #C ij # ij %u uk + c i %u kj + c j ik " p %x ( k %x k %x k ) + M ij NLT ij : originates in Oldroyd derivative not negligible, must be modeled What about: CT ij 1) Reynolds stresses? 2) Turbulent kinetic energy? : originates in advective term, it is negligible no need for modeling DNS: Housiadas et al (25), Li et al (26) JNNFM 3) Dissipation of turbulent kinetic energy or of Reynolds stresses? 9
10 Transport equation for the Reynolds stresses and k! "u i u j "t "u +!U i u j k = P ij +Q ij +Q V ij + D ij,n + # ij $!% N V ij $!% ij "x k Q V ij =! ( u i " jk, p + u j " ) ik, p!x k! V ij = 1 % " # jk, p & $u i + # ik, p $x k $u j $x k Viscoelastic turbulent transport due to fluctuations polymeric stresses ( ) Viscoelastic work of polymer chains: dissipation of energy plus stored free energy (< ou >)! Dk Dt +!u "U i "k iu k = #!u i # "pu i " 2 k "u +$ s #$ i "u i s + "% ik, p # % ik, p "x k dx i "x i "x i "x i "x k "x k "x k u i "u i "x k ε Ν Q V! V 1
11 Transport equation of ε Ν "u 2! i " $ # Du i s & "x m "x m % Dt )+ 2! s ( "u i "x m " $ & % "x m #uk "U i "x k "u )+ 2! i " $ # "u u i k s & ( "x m "x m % "x k " #"p& "ui 2 " #" & 2 " #"+ "ui ui "ui +2! s % () 2! s % () 2! ik, p % 2 s "x m "x m $ "x i "x m "x m $ "x k "x m "x m $ "x k ) ( & ( = New term As for Newtonian fluids, the whole equation will be approximated 11
12 Budget of Reynolds stress (a) uu Budget Pxx Q N xx Q V xx D N xx!xx " N xx " V xx (c) vv Budget Pyy Q N yy Q V yy D N yy!yy " N yy " V yy y y+ - Need to model viscoelastic stress work - Need to model pressure strain (effect of elasticity)- Advanced mod. - Viscoelastic turbulent transport is not so important 12
13 Viscoelastic work! V " 1 # $ ik, p %u i = & ( p %x k # C f C + c ik mm mm ) ( ) %u i ( ) %u i %x k + c ik f C mm + c mm %x k + -, 35 Assumption: (same reasons) C ik f ( C mm + c mm )!u i < c ik f ( C mm + c mm )!u i!x k!x k 3 C kk 2 c kk C kk > 2 c kk! V " # p $% f C &u ( mm )c i ik = C &x! V k # p $% f ( C mm ) NLT nn y + C! V = 1.76 (from DNS- next slide) Needs model 13
14 Performance of the viscoelastic work model 1 5 (e)! V > Negative sign: different definitions -5-1! ( V + Re " ) ! V (DNS) -1.7f(C kk )NLT kk /2 versus f ( C kk ) NLT ij y + 1 3
15 Modeling NLT ij 1 Key ideas: 1) Write down exact equation- complex 4 lines long 2) Get inspiration from it: physical insight, trial-and-error, luck 3) Make assumptions 4) Do a priori testing of each term 5) Select appropriate combination and dimensional homogeneity 6) Try in code - under investigation f ( C mm ) NLT ij!!c kj!c kj!c kj u i u m + u i u m " Coef # u i u m!x m!x m!x m $ = function S ij,w ij,c ij," N ij, #u iu j, #C ij, #NLT ij &, M ij,u i u j % #x k #x k #x n ) ( f ( C mm ) NLT ij! #u = f µ1 C k u n #C ij "1 + C u 2 E3 $ C 2 kk u i u j + C & ) 14 u i u k W kn C nj + u j u k W kn C ni + u k u i W jn C nk ( #x n #x k % % ( ), + 15
16 Model for NLT ij 2 Blue model Red model C E3 =.4;C! 1 = 3;C "14 =.15 C E3 =.35;C! 1 = ;C "14 = f(ckk)nltxy f(ckk)nltii Pred f(ckk)nltxy Pred f(ckk)nltii f(ckk)nltyy Pred f(ckk)nltyy 25 2 f(ckk)nltxy f(ckk)nltii Pred f(ckk)nltxy Pred f(ckk)nltii f(ckk)nltyy Pred f(ckk)nltyy y y + 1 f µ1 = ( 1! exp (! y ) ) 2 16
17 Viscoelastic turbulent transport Q V! "# ik, pu i = $ p "x k % " & C ik f ( C mm + c mm )u i + c ik f ( C mm + c mm )u i "x k ( ) C kk > f 2 c kk ( Ĉ ) mm = L 2! 3 ( ) L 2! C mm + c mm Weak coupling between c kk and c ij, u i C ik f ( C mm + c mm )u i < c ik f ( C mm + c mm )u i Neglect of this term is irrelevant because non-neglected term is modeled Q V =! p " # #x $% f ( C mm )CU iik & k Needs model General case ( CU ) ijk 17
18 Model for CU ijk - Same modelling approach as with NLT ij 4 f ( C mm )CU ijk! + % $C = f µ2 "C #1 u i u kj $C m + u j u ik ( m " C # - 7 & $x m $x m )! f C mm, ( ) + ± u 2 C ± u 2 C j ik i jk,-.. / / 3 2 QFii2 Pred QFii2 C!1 =1.3;C!7 =.37 " % f µ2 = 1! exp $! y+ # 26.5 & y + 18
19 Final equations: low Re k-ε type model for channel flow Momentum: d % dy! s & du dy + " # $uv ( p,xy # dp ) dx =! xy, p = " p # f ( C kk )C xy du f ( C kk )C xy =!C yy dy +!NLT xy f ( C kk )C yy =!NLT yy +1 du f ( C kk )C xx = 2!C xy dy +!NLT xx +1 f ( C kk )C zz =!NLT zz +1 f ( C kk ) = L 2! 3 ( ) L 2! C xx + C yy + C zz Reynolds stress: du!"uv = "# T with dy! = C f T µ µ k 2!" N 19
20 k and ε transport equations = d + % dy! + "# ( T s dk. d P k 1 "!2 N 1 "D +! p -,& $ k ) dy/ dy, f ( C mm ) 3 CU nny 2 " d k %! N =! N + D N D N = 2! s $ # dy & ( ). / 1! f C mm p 3 2 NLT nn 2 = d,& dy! + "# T.( s - $ % ) + d!% N dy /!% N 1+ " f 1 C %1 k P k " 2 " f C % N 2 2 % 2 k + "E + E 3 p E =! s " # T 1$ f µ % ( ) d 2 U & dy 2 ( ) 2 f 1 = 1 f 2 = 1!.3exp ( 2!R ) T ( " % + f µ = 1! exp $!y+ - ) # 26.5 &, 2 based on Newtonian model of Nagano & Hishida (1984) 2
21 Some results 1: Re τ = 395; We τ = 25; β=.9, L 2 =9 3 u Virk s asymptote X DNS Black: Newtonian! =! s +! p! =! wall! =! wall ;same!q;re " = Newtonian log law FENE-P simulations C!1 = 1.3;C! 7 =.37;C "14 = 1.5 # 1 $4 C! V = 1.76 With! p C E 3 = 1.93 " 1 # y + Without! p C E 3 = 2.86 " 1 #4 21
22 Some results 2: Re τ = 395; We τ = 25; β=.9, L 2 =9 5 k + k + uv + 1 uv y y 1 22
23 Some results 3: Re τ = 395; We τ = 25; β=.9, L 2 =9 25 NLT ii 2 NLT ii 35 NLT xy 3 NLT xy y y + 23
24 Some results 4: Re τ = 395; We τ = 25; β=.9, L 2 =9 NLT xx NLT yy 15 5 NLT xx NLT yy y y + 24
25 Some results 5: Re τ = 395; We τ = 25; β=.9, L 2 =9 C xy C yy 25 C xy C yy y y + 25
26 Some results 6: Re τ = 395; We τ = 25; β=.9, L 2 =9 5 C xx 1 Stresses + C xx 4 Stresses +.8 -uv + +! p,xy +! s,xy SUMSTRESS du f ( C kk )C xx = 2!C xy dy +!NLT +1 xx y y 26
27 Some results 7: Re τ = 395; We τ = 25; β=.9, L 2 =9 8 + CU iiy CU ii y + 27
28 Conclusions - Developed simplified k-ε model: code and closures are working - Viscoelastic stress power well modeled by NLT ij - Viscoelastic turbulent transport (CU ijk ) is not that relevant at 18% - NLT ij is also required for C ij -Closure for NLT ij has deficiencies and needs significant improvement - Excessive dissipation of turbulence - Need to model viscoelastic turbulence production close to wall - Isotropic turbulence does not allow a good model - Need to consider anisotropic turbulence: anisotr. k-ε and RSM - Closure for CU ijk is fair but also needs improvement: small impact 28
F. T. Pinho. Centro de Estudos de Fenómenos de Transporte, Universidade do Porto & Universidade do. Minho, Portugal
F. T. Pinho Centro de Estudos de Fenómenos de Transporte, Universidade do Porto Universidade do Minho, Portugal Colaboradores P. R. Resende Centro de Estudos de Fenómenos de Transporte, Universidade do
More informationA k-! LOW REYNOLDS NUMBER TURBULENCE
A - LOW REYOLDS UMBER TURBULECE MODEL FOR TURBULET CHAEL FLOW OF FEE-P FLUIDS P. R. Resende Centro de Estudos de Fenómenos de Transporte, Universidade do Porto, Portugal F. T. Pinho Centro de Estudos de
More informationDevelopment of a Low-Reynolds-number k-ω Model for FENE-P Fluids
Flow Turbulence Combust 2013) 90:69 94 DOI 10.1007/s10494-012-9424-x Development of a Low-Reynolds-number k-ω Model for FENE-P Fluids P. R. Resende F. T. Pinho B. A. Younis K. Kim R. Sureshkumar Received:
More informationSOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG
SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG L. Thais Université de Lille Nord de France, USTL F9 Lille, France Laboratoire de Mécanique de Lille CNRS, UMR 817
More informationPERFORMANCE OF THE k-ε AND REYNOLDS STRESS MODELS IN TURBULENT FLOWS WITH VISCOELASTIC FLUIDS.
Proceedings of the th Brazilian Congress of Thermal Sciences and Engineering -- ENCIT 6 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Curitiba, Brazil,- Dec. -8, 6 Paper CIT6-8 PERFORMANCE
More informationDEVELOPMENT OF A REYNOLDS STRESS MODEL TO PREDICT TURBULENT FLOWS WITH VISCOELASTIC FLUIDS
CONFERÊNCIA NACIONAL DE MÉTODOS NUMÉRICOS EM MECÂNICA DOS FLUIDOS E TERMODINÂMICA 6 FCT-UNL, Monte de Caparica, 8 e 9 de Junho, 6 APMTAC, Portugal, 6 DEVELOPMENT OF A REYNOLDS STRESS MODEL TO PREDICT TURBULENT
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationTURBULENCE MODELING FOR POLYMER SOLUTIONS
TURBULENCE MODELING FOR POLYMER SOLUTIONS Fernando T. Pinho 1 Mohammadali Masoudian 1 Carlos B. da Silva 2 1 Transport Phenomena Research Center (CEFT), Mechanical Engineering Department, Faculdade de
More informationOutline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow.
Victor S. L vov: Israeli-Italian March 2007 Meeting Phenomenology of Wall Bounded Turbulence: Minimal Models In collaboration with : Anna Pomyalov, Itamar Procaccia, Oleksii Rudenko (WIS), Said Elgobashi
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationHEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS
HEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS Paulo J. Oliveira Departamento de Engenharia Electromecânica, Universidade da Beira Interior Rua Marquês D'Ávila e Bolama, 600
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationModel-based analysis of polymer drag reduction in a turbulent channel flow
2013 American Control Conference ACC Washington, DC, USA, June 17-19, 2013 Model-based analysis of polymer drag reduction in a turbulent channel flow Binh K. Lieu and Mihailo R. Jovanović Abstract We develop
More informationA Simple Turbulence Closure Model
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationPolymer stress tensor in turbulent shear flows
Polymer stress tensor in turbulent shear flows Victor S. L vov, Anna Pomyalov, Itamar Procaccia, and Vasil Tiberkevich Department of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100,
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationEvents of High Polymer Activity in Drag Reducing Flows
Flow Turbulence Combust (27) 79:23 32 DOI.7/s494-7-975-5 Events of High Polymer Activity in Drag Reducing Flows C. Tesauro B. J. Boersma M. A. Hulsen P. K. Ptasinski F. T. M. Nieuwstadt Received: May 26
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationXIIIth International Workshop on Numerical Methods for non-newtonian Flows. Hôtel de la Paix, Lausanne, 4-7 June 2003 PROGRAMME
XIIIth International Workshop on Numerical Methods for non-newtonian Flows Hôtel de la Paix, Lausanne, 4-7 June 23 PROGRAMME Wednesday Afternoon, 4 June 23... 2 Novel Numerical Methods I (Chairman : Roland
More informationA Simple Turbulence Closure Model. Atmospheric Sciences 6150
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U
More informationDept. Engineering, Mechanical Engineering, University of Liverpool Liverpool L69 3GH, UK,
R. J. Poole pt. Engineering, Mechanical Engineering, University of Liverpool Liverpool L69 3GH, UK, robpoole@liv.ac.uk M. A. Alves partamento de Engenharia uímica, CEFT, Faculdade de Engenharia da Universidade
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationTurbulence Solutions
School of Mechanical, Aerospace & Civil Engineering 3rd Year/MSc Fluids Turbulence Solutions Question 1. Decomposing into mean and fluctuating parts, we write M = M + m and Ũ i = U i + u i a. The transport
More informationSOME THOUGHTS ON DIFFERENTIAL VISCOELASTIC MODELS AND THEIR NUMERICAL SOLUTION
CONFERENCE ON COMPLEX FLOWS OF COMPLEX FLUIDS University of Liverpool, UK, March 17-19, 2008 SOME THOUGHTS ON DIFFERENTIAL VISCOELASTIC MODELS AND THEIR NUMERICAL SOLUTION Paulo J. Oliveira Universidade
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationInterpretation of the mechanism associated with turbulent drag reduction in terms of anisotropy invariants
J. Fluid Mech. (27), vol. 577, pp. 457 466. c 27 Cambridge University Press doi:1.117/s221127583 Printed in the United Kingdom 457 Interpretation of the mechanism associated with turbulent drag reduction
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationAdvances in the Analysis and Prediction of Turbulent Viscoelastic Flows
Old Dominion University ODU Digital Commons CCPO Publications Center for Coastal Physical Oceanography 2014 Advances in the Analysis and Prediction of Turbulent Viscoelastic Flows T. B. Gatski Old Dominion
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationOn the role of pressure in elasto-inertial turbulence
Journal of Turbulence ISSN: (Print) 1468-5248 (Online) Journal homepage: http://www.tandfonline.com/loi/tjot2 On the role of pressure in elasto-inertial turbulence Vincent E. Terrapon, Yves Dubief & Julio
More informationEuclidean invariance and weak-equilibrium condition for the algebraic Reynolds stress model
J. Fluid Mech. 2006), vol. 569, pp. 399 408. c 2006 Cambridge University Press doi:10.1017/s0022112006003041 Printed in the United Kingdom 399 Euclidean invariance and weak-equilibrium condition for the
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationMaximum drag reduction asymptotes and the cross-over to the Newtonian plug
J. Fluid Mech. (26), vol. 551, pp. 185 195. c 26 Cambridge University Press doi:1.117/s221125795 Printed in the United Kingdom 185 Maximum drag reduction asymptotes and the cross-over to the Newtonian
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationHYBRID LES-RANS: Inlet Boundary Conditions for Flows With Recirculation
1 HYBRID LES-RANS: Inlet Boundary Conditions for Flows With Recirculation Lars Davidson Div. of Fluid Dynamics Dept. of Applied Mechanics Chalmers University of Technology, Göteborg, Sweden E-mail lada@chalmers.se
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationReynolds stress budgets in Couette and boundary layer flows
Reynolds stress budgets in Couette and boundary layer flows By Jukka Komminaho and Martin Skote Dept. of Mechanics, KTH, SE-1 44, Stockholm, Sweden Reynolds stress budgets for both Couette and boundary
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
More informationVictor S. L vov, Baku, January 21, Turbulent Drag Reduction by Dilute Solution of Polymers
Victor S. L vov, Baku, January 21, 2009 Turbulent Drag Reduction by Dilute Solution of Polymers In collaboration with T.S. Lo (TL), A. Pomyalov (AP), I. Procaccia (IP), V. Tiberkevich (VT). Weizmann Inst.,
More informationGoverning Equations for Turbulent Flow
Governing Equations for Turbulent Flow (i) Boundary Layer on a Flat Plate ρu x Re x = = Reynolds Number µ Re Re x =5(10) 5 Re x =10 6 x =0 u/ U = 0.99 層流區域 過渡區域 紊流區域 Thickness of boundary layer The Origin
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationPolymer maximum drag reduction: A unique transitional state
Center for Turbulence Research Annual Research Briefs 47 Polymer maximum drag reduction: A unique transitional state By Y. Dubief, C. M. White, E. S. G. Shaqfeh AND V. E. Terrapon. Motivation and objectives
More information5 The Oldroyd-B fluid
5 The Oldroyd-B fluid Last time we started from a microscopic dumbbell with a linear entropic spring, and derived the Oldroyd-B equations: A u = u ρ + u u = σ 2 pi + η u + u 3 + u A A u u A = τ Note that
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
More informationAtmospheric Boundary Layer Studies with Unified RANS-LES and Dynamic LES Methods
5st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 7 - January 23, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 23-747 Atmospheric Boundary Layer Studies with
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationOutline. Motivation Governing equations and numerical methods Results: Discussion:
Bifurcation phenomena in strong extensional flows (in a cross-slot geometry) F. A. Cruz 1,*, R. J. Poole 2, F. T. Pinho 3, P.J. Oliveira 4, M. A. Alves 1 1 Departamento de Engenharia Química, CEFT, Faculdade
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationarxiv: v1 [physics.flu-dyn] 16 Nov 2018
Turbulence collapses at a threshold particle loading in a dilute particle-gas suspension. V. Kumaran, 1 P. Muramalla, 2 A. Tyagi, 1 and P. S. Goswami 2 arxiv:1811.06694v1 [physics.flu-dyn] 16 Nov 2018
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationNumerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model
Numerical Simulation of Turbulent Buoyant Helium Plume by Algebraic Turbulent Mass Flux Model Hitoshi Sugiyama 1),Naoto Kato 1), Masahiro Ouchi 1) 1) Graduate School of Engineering,Utsunomiya University
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More informationChristel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011
Notes by: Andy Thaler Christel Hohenegger A simple model for ketchup-like liquid, its numerical challenges and limitations April 7, 2011 Many complex fluids are shear-thinning. Such a fluid has a shear
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationUniformity of the Universe
Outline Universe is homogenous and isotropic Spacetime metrics Friedmann-Walker-Robertson metric Number of numbers needed to specify a physical quantity. Energy-momentum tensor Energy-momentum tensor of
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
More informationElliptic relaxation for near wall turbulence models
Elliptic relaxation for near wall turbulence models J.C. Uribe University of Manchester School of Mechanical, Aerospace & Civil Engineering Elliptic relaxation for near wall turbulence models p. 1/22 Outline
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationGeneration of initial fields for channel flow investigation
Generation of initial fields for channel flow investigation Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (Mai 2) In the framework of the DFG-funded research
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationApplication of a Volume Averaged k-ε Model to Particle-Laden Turbulent Channel Flow
Application of a Volume Averaged k-ε Model to Particle-Laden Turbulent Channel Flow J. D. Schwarzkopf C. T. Crowe P. Dutta e-mail: dutta@mme.wsu.edu School of Mechanical and Materials Engineering, Washington
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More informationLast update: 15 Oct 2007
Unsteady flow of an axisymmetric annular film under gravity HOUSIADAS K, TSAMOPOULOS J -- PHYSICS OF FLUIDS --10 (10)2500-2516 1998 Yang XT, Xu ZL, Wei YM Two-dimensional simulation of hollow fiber membrane
More informationUNIVERSITY OF WALES INSTITUTE OF NON-NEWTONIAN FLUID MECHANICS CONFERENCE ON INDUSTRIAL RHEOLOGY. Hoole Hall Hotel, Chester 5-7 April 2004 PROGRAMME
UNIVERSITY OF WALES INSTITUTE OF NON-NEWTONIAN FLUID MECHANICS CONFERENCE ON INDUSTRIAL RHEOLOGY Hoole Hall Hotel, Chester 5-7 April 24 PROGRAMME 27. Mar. 24 PROGRAMME Monday 5 April - Afternoon 2:3 Lunch
More informationTurbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media
Copyright 011 Tech Science Press CMES, vol.71, no., pp.149-155, 011 Turbulentlike Quantitative Analysis on Energy Dissipation in Vibrated Granular Media Zhi Yuan Cui 1, Jiu Hui Wu 1 and Di Chen Li 1 Abstract:
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationStrategy in modelling irregular shaped particle behaviour in confined turbulent flows
Title Strategy in modelling irregular shaped particle behaviour in confined turbulent flows M. Sommerfeld F L Mechanische Verfahrenstechnik Zentrum Ingenieurwissenschaften 699 Halle (Saale), Germany www-mvt.iw.uni-halle.de
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationDYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1551-1555 1551 DYNAMIC STABILITY OF NON-DILUTE FIBER SHEAR SUSPENSIONS by Zhan-Hong WAN a*, Zhen-Jiang YOU b, and Chang-Bin WANG c a Department of Ocean
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationRheology of Soft Materials. Rheology
Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ
More informationSkin-Friction Drag Reduction by Dilute Polymer Solutions in Turbulent Channel Flow
Skin-Friction Drag Reduction by Dilute Polymer Solutions in Turbulent Channel Flow by Dong-Hyun Lee A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationMathematical Theory of Non-Newtonian Fluid
Mathematical Theory of Non-Newtonian Fluid 1. Derivation of the Incompressible Fluid Dynamics 2. Existence of Non-Newtonian Flow and its Dynamics 3. Existence in the Domain with Boundary Hyeong Ohk Bae
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More information