SOME THOUGHTS ON DIFFERENTIAL VISCOELASTIC MODELS AND THEIR NUMERICAL SOLUTION
|
|
- Cynthia Horton
- 5 years ago
- Views:
Transcription
1 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 da Beira Interior Dep. Eng. Electromecânica UBI
2 TOPICS Derivation of a generalized Oldroyd-B model The traceless stress tensor formulation Law-of-the-wall for viscoelastic flows
3 DERIVATION OF GENERALISED MODEL First TOPIC
4 GOVERNING EQUATIONS iu= 0 Du ρ Dt τ η 2 = p+ i τtot = p+ s u+ iτ
5 OLDROYD-B EQUATIONS Differential formulation Original: τ λ τ η D λ D τ + ( τ ) = 2 + tot tot 0 r Alternative: = + = 2 D+ τ τ τ η tot s p s τ + λ ( τ ) = 2η D p τ With: η0 = ηs+ ηp ηs λr β = = η λ 0
6 OLDROYD-B: molecular formulation ( ) A= A I η p λ τ = ( A I) λ (1) (2) ( ) λ η A p η I τ = A I = p + 2D λ τ λτ = ηp + 2D η λ τ + τ = 2η pd p ( λ and ηp are constants) with: I = 2D
7 Equations with conservative property Idea based on conservativeness property (at continuum level) Convected variables should appear inside derivative: D λ τ ( ) ( ) Dt λ τ = + t div( λuτ ) so that balances of fluxes (out-in): ( ) n+ 1 ( ) n V P λτ λτ P P + λu iaτ λu iaτ = t flux e flux w e W w P e w E... x
8 Example of Conservative stress models FENE-MCR: not conservative FENE-CR: conservative PTT: not conservative (apparently) λ τ + ( τ ) = 2η p D f λ τ + τ = 2η p D f f τ + λ τ += η D ( ) 2 p FENE-P: original equation also apparently not conservative D lnz λ( τ ) + Z τ λ + ( 1 εb) nktλ = 2( 1 εb) nktλ Dt τ I D
9 Example of Conservative stress models: FENE-P (cont) with function: Z 3 tr( τ ) = 1+ ( 1 εb) + 2 b 3nkT ε = bb ( + 2) But re-arrangement shows conservativeness: λ ( η / ) p aη D a f p τ + τ = 2 D I f f Dt with: f Z a b+ 5 b = ηp = nktλ b + 2 b + 5
10 CONCLUSION So, in general all previous equations can be cast in terms of variable relaxation times and viscosity coefficients: as: η f p λ λ ηp τ + τ = 2 D+... f f f with different functions for the various models:
11 TYPICAL STRESS FUNCTIONS: FENE-CR: f [ τ] = L 2 + ( λ / η ) Tr( τ ) 0 p0 L 2 3 FENE-P: f 3 λ = Tr( ) b [ τ] 0 τ + 2 3aηp 0 PTT: f λε = 1 + Tr( τ ) [ τ] 0 η p0
12 GENERALISATION: Now, both viscosity and relaxation time are variable, functions of invariants of τ : η = η ( τ ) p p λ = λ( τ ) We shall use for the constant zero-shear rate values: In general: ( fa) η λ p0 0 Df = A+ f A Dt
13 DERIVATION OF GENERALISED MODEL: Assume as still valid: ( ) A= A I η p λ τ = ( A I) λ (1) Equilibrium between rotation/stretching and relaxation (2) Kramers expression Upper convected derivative of Eq. (2) gives: Dη Dη λτ ηpa ηpi A+ ηpa I+ ηpd Dt Dt ( ) ( ) ( ) p p = = 2
14 DERIVATION OF GENERALISED EQUATION Substituting Eq. (1): Dη ( ) p λτ = 2η ( ) pd+ A I ηp Dt Eq. (2) again: λ Dη ( ) p λτ + τ = 2ηpD+ τ η Dt Finally: Dη λ + p = η Dt p ( λ τ ) τ 1 2 p η p A I ( ) D λ
15 GENERALISED EQUATION SPECIAL CASE Recall: η p0 λ η = λ= 0 p f f assume same function f τ for η and λ. Dη λ [ ] + p = η Dt p ( λ τ ) τ 1 2 Gives: p η Dλ λ Dη λ + + = p D 0 p τ τ τ 1 2ηpD Dt ηp0 Dt
16 MANIPULATION FOR SPECIAL CASE: ( λ / ) D( η / ) p f λ D f λ η τ τ τ D f Dt ηp0 Dt f p0 τ+ τ + τ 1 = 2 τ fτ 2η 0 p 0 λ + = D the PTT equation!!
17 COMPARISON WITH FENE-P: FENE-P under conservative form: ( λ τ ) τ 2 p D η τ τ η D p I + = + = Dt or in a more compact and conservative form: ( λ τ ) τ ( η Ι ) p with: λ = λ 0 /f η = aη f p p0 / Present model: + = ( λ τ ) τ 2 η p D λτ + η p D η Dt p Advantage: permits to control functions separately
18 COMPARISON WITH FENE-P: material functions Stress growth upon inception of shear flow: viscosity 1 FENE-P and FENE-M, L 2 =100 2 L b = + 3 η + (t,γ)/η p We=1.0 We=10 We= t/λ FENE-P (solid lines) goes beyond the linear viscoelastic limit envelope; cf. DeAguiar (1983)
19 COMPARISON WITH FENE-P: material functions Stress growth upon inception of shear flow: viscosity 1 FENE-P, L 2 =10 1 FENE-M, L 2 = η + (t,γ)/η p We=1.0 We=10 We=100 η + (t,γ)/η p We=1.0 We=10 We= t/λ t/λ Same as before, for lower extensibility, L2=10
20 COMPARISON WITH FENE-P: material functions Similar for start up of first normal stress coefficient 10 FENE-P, L 2 = FENE-M, L 2 = Ψ Ψ 1 + (t,γ) We=0.1 We=1.0 We=10 We= We=0.1 We=1.0 We=10 We= t/λ t/λ
21 COMPARISON WITH FENE-P: material functions Start-up of uniaxial elongational flow 200 FENE-P and FENE-M, We=5 η E L2=25 L2=50 L2= t Modified model allows a less steep rise of extensional viscosity (more realistic?? cf. contrast FENE FENE-P: Keunings 97; van Heel et al 98; Herrchen &Ottinger 97)
22 COMPARISON WITH FENE-P: material functions Modified model in uniaxial elongational flow 50 FENE-M, We=5, L 2 = η E L2vis=25 L2vis=10 L2vis= t Variation of L2 in viscosity funtion permits some control
23 TRACELESS STRESS TENSOR FORMULATION Second TOPIC
24 TRACELESS STRESS TENSOR FORMULATION P-correction solver Standard equations of motion: u x j j = ρu ρuu p i + j i = + t x x x 0 τ j i j ij Coupling term (1 β )De Traceless tensor: = 1τ δ ij ij 3 kk ij τ τ so that: ( τ ) = = 0 Tr τ kk
25 TRACELESS MOMENTUM EQ. Substitution in original eqs. gives: ( 1 ρ u ρuu p τ ) τ i j i 3 kk + = + t x x x j i j ρ u ρuu i j i p τ ij + = + t x x x j i j ij with: p = p 1τ 3 kk (modified pressure)
26 TRACELESS STRESS: OLDROYD-B Oldroyd-B equation, indicial notation: Dτ u u τ + λ = η + + λg ij i j ij p ij Dt xj x i with generation tensor term: G u = τ + τ j ij ik jk xk Using the TST we get: u x Dτ ij u u i j τ ij+ λ = η p + + λ ij 3 δij Dt xj x i i k ( 2 G G )
27 TRACELESS STRESS TENSOR: increased viscosity with: and: uj ui G ij = τ ik + τ jk xk xk ul G ll = 2τ lk 2G x l (in turbulence modelling, G= τ is the generation rate) lk xk u k An increased viscosity arises: = + 1λτ p p 3 kk η η very high, where it matters! (link with AVSS )
28 Evolution of the trace Relations for the evolution of the trace of the stress: Contracting indices: G= G + Dτkk uk τkk + λ = 2η p + 2λG Dt x k = τ kk u xj = 0 enabling calculation of the trace j
29 TRACELESS STRESS TENSOR for OLDROYD-B Resumé in tensor notation: τ + λ ( τ ) = 2η D p 1 τ = τ Tr( τ ) Ι 3 D = D 1Tr( D) Ι 3 η = η + 1λ τ and the extended definition: Dτ ( τ ) = τ + τ 3 τ : Dt Tr p p 3 ( ) { T 2 i u u i ( u) I}
30 TRACELESS STRESS TENSOR: RESULTS (1) UCM in channel flow, 20x20, L=10, H=1 Total number of iterations to solve pressure equation: We Outer iterations Total number of inner iterations for pressure Standard method Traceless method Ratio
31 TRACELESS STRESS TENSOR: RESULTS (1a) UCM: 2D channel flow, 20x20 pressure iterations Standard Traceless 0 total iterations We We
32 TRACELESS STRESS TENSOR: RESULTS (2) UCM, 4:1 plane contraction flow (mesh 2960 CV) Total number of iterations to solve pressure equation: We 0.0 Outer iterations 798 Total number of inner iterations for pressure Standard method 8641 Traceless method -- Ratio Ratio proportional to relative CPU times
33 TRACELESS STRESS TENSOR: RESULTS (2) UCM, 4:1 plane contraction flow (mesh 2960 CV) pressure iterations UCM: 4:1 plane contraction Standard Traceless We (ratio proportional to relative CPU times)
34 TRACELESS STRESS TENSOR: CONCLUSIONS Advantages: TST offers better coupling and much less CPU Number of pressure iterations approximately constant Problems: Creation of strong normal stress normal to a wall Difficulty to implement such BCs (oscillations) Standard: TST: e.g. UCM, channel aligned with x: τ xx xx ληγ τ = 2 = 2 yy 0 τ = ληγ τ 2 2 yy = 3ληγ xy τxy = ηγ τ = ηγ Problems are particular to our FVM (collocated) (?) Would like to see other attempts with FEM
35 LAW-OF-WALL FOR VISCOELASTIC FLOWS Third TOPIC
36 LAW-OF-WALL FOR VISCOELASTIC FLOWS At present, 2D meshes are so fine that any attempt to extend to 3D is out of question (ex: cylinder) NC=44280; even a coarser mesh NC=11040x100=
37 LAW-OF-WALL FOR TURBULENT, FLOWS Flow depends on: Dimensional analysis gives: and a wall layer: u u τ y τw ρ µ τ u w τ ρ y + ρyu µ u u + = f1( y ) f ( ) uv = τ y = f 2 3( y ) y y u τ τ = 0 τ = τl + τt const. y Small acceleration implies: (constant-stress layer) Experiments show, for large y+: f3 1 uτ = uv f2 const. = 1/ Κ Κ= 0.41 τ u u = τ y Κy u u 1 = ln( Ey ) Κ + + u τ (the log-law)
38 LAW-OF-WALL FOR VISCOELASTIC FLOWS Carry same ideas for viscoleastic flows near walls List of dependent variables includes relaxation time: y p x ρ η λ Giving the characteristic scales: tc λ y = η = c c 2 λp λ x u = η p Instead of a logarithmic variation, power law attempts were followed in the next fittings of shear flows (following Zagarola et al., Phy Fluids 1997) x
39 LAW-OF-WALL FOR VISCOELASTIC FLOWS Pipe flow: Linear SPTT
40 LAW-OF-WALL FOR VISCOELASTIC FLOWS Pipe flow: Exponential SPTT
41 LAW-OF-WALL FOR VISCOELASTIC FLOWS Channel flow: Leonov fluid
42 LAW-OF-WALL FOR VISCOELASTIC FLOWS Lot to be done. In future attempts, a wall law should be sought for the stresses (and not velocity), based on the notion of local equilibrium
43 CONCLUSIONS Theoretical (empirical): a naive derivation led to a modified FENE-P equation with some advantages for time-dependent flows. It is the PTT when functions are the same. Numerics: the traceless approach allows tighter coupling of equations in decoupled methods. Unsolved problems with normal stresses at wall. Pratical applications (3D): proposal for a loglaw -like treatment to bridge the near wall BLs.
44 ACKNOWLEDGMENTS: Prof. Manuel A. Alves (Univ. of Porto, FEUP) (law-of-wall) Project POCI/EME/58657/2004
Dept. 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 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 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 informationON THE EFFECT OF CORNER CURVATURE IN LAMINAR FLOW THROUGH A CROSS-SLOT CHANNEL
XXIX CILAMCE Iberian Latin American Congress on Computational Methods in Engineering November 4-7, 2008, Maceió-Alagoas, Brazil ON THE EFFECT OF CORNER CURVATURE IN LAMINAR FLOW THROUGH A CROSS-SLOT CHANNEL
More informationVISCOELASTIC FLOW IN AXISYMMETRIC CONTRACTIONS: THE EFFECT OF CONTRACTION RATIO
Proceedings of the 0th Brazilian Congress of Thermal Sciences and Engineering -- ENCIT 004 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, io de Janeiro, Brazil, Nov. 9 -- Dec. 03, 004 VISCOELASTIC
More informationElastic-driven instabilities in microfluidic flow-focusing devices
Elastic-driven instabilities in microfluidic flow-focusing devices M. S. N. Oliveira 1, F. T. Pinho 2, M. A. Alves 1 1: Departamento de Engenharia Química, CEFT Faculdade de Engenharia da Universidade
More informationThe FENE-L and FENE-LS closure approximations to the kinetic theory of finitely extensible dumbbells
J. Non-Newtonian Fluid Mech. 87 (1999) 179 196 The FENE-L and FENE-LS closure approximations to the kinetic theory of finitely extensible dumbbells G. Lielens, R. Keunings,V.Legat CESAME, Division of Applied
More informationCentro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal
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
More informationNumerical study of the flow of FENE-CR viscoelastic fluids in a planar cross-slot geometry
Numerical study of the flow of FENE-CR viscoelastic fluids in a planar cross-slot geometry G.N. Rocha 1, R.J. Poole 2, M.A. Alves 3 and P.J. Oliveira 1 1 Universidade da Beira Interior, Departamento de
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 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 informationA numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 6 (2016, 152 161 Research Article A numerical approximation with IP/SUPG algorithm for P-T-T viscoelastic flows Lei Hou a, Yunqing Feng a,, Lin
More informationPaper CIT UNSTEADY VISCOELASTIC FLOW PAST A CYLINDER IN A CHANNEL
Proceedings of the 10o Brazilian Congress of Thermal Sciences and Engineering -- ENCIT 004 Braz. Soc. of Mechanical Sciences and Engineering -- ABCM, Rio de Janeiro, Brazil, Nov. 9 -- ec. 03, 004 Paper
More informationBoundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity
Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with
More informationPhan-Thien-Tanner Modeling of a Viscoelastic Fluid in the Stick-Slip Scenario. Andrew Griffith
Phan-Thien-Tanner Modeling of a Viscoelastic Fluid in the Stick-Slip Scenario Andrew Griffith Supervisor: Bruce A. Finlayson Department of Chemical Engineering University of Washington March 1, 7 Introduction
More informationStudy of steady pipe and channel flows of a single-mode Phan-Thien Tanner fluid
J. Non-Newtonian Fluid Mech. 101 (2001) 55 76 Study of steady pipe and channel flows of a single-mode Phan-Thien Tanner fluid Manuel A. Alves a, Fernando T. Pinho b,, Paulo J. Oliveira c a Departamento
More informationOldroyd Viscoelastic Model Lecture Notes
Oldroyd Viscoelastic Model Lecture Notes Drew Wollman Portland State University Maseeh College of Engineering and Computer Science Department of Mechanical and Materials Engineering ME 510: Non-Newtonian
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationOUTLINE: Manuel A. Alves. Paulo J. Oliveira. Fernando T. Pinho
OUTLINE: Paulo J. Oliveira Universidade da Beira Interior Dep. Eng. Electromecânica UBI Fernando T. Pinho Universidade do Minho Escola de Engenharia Manuel A. Alves Distinctive aspects of Computational
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 informationCESAME, Louvain, Belgium
Thursday Morning Thursday :2 Johann-Peter-Hebel FM3 PGD-based parametric solutions of double-scale Fokker-Planck equations: application to suspensions involving aggregates of rods Francisco Chinesta, Emmanuelle
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
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 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 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 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 informationNumerical and Experimental Results
(G) Numerical and Experimental Results To begin, we shall show some data of R. D. Moser, J. Kim & N. N. Mansour, Direct numerical simulation of turbulent flow up to Re τ = 59, Phys. Fluids 11 943-945 (1999)
More informationA Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface
A Multigrid LCR-FEM solver for viscoelastic fluids with application to problems with free surface Damanik, H., Mierka, O., Ouazzi, A., Turek, S. (Lausanne, August 2013) Page 1 Motivation Polymer melts:
More informationLecture Administration. 7.2 Continuity equation. 7.3 Boussinesq approximation
Lecture 7 7.1 Administration Hand back Q3, PS3. No class next Tuesday (October 7th). Class a week from Thursday (October 9th) will be a guest lecturer. Last question in PS4: Only take body force to τ stage.
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 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 informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ
Chapter 2 Mathematical Models The following sections present the equations which are used in the numerical simulations documented in this thesis. For clarity, equations have been presented in Cartesian
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationMethod for time-dependent simulations of viscoelastic flows: vortex shedding behind cylinder
J. Non-Newtonian Fluid Mech. 101 (2001) 113 137 Method for time-dependent simulations of viscoelastic flows: vortex shedding behind cylinder Paulo J. Oliveira Departamento de Engenharia Electromecânica,
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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationTransport processes. 7. Semester Chemical Engineering Civil Engineering
Transport processes 7. Semester Chemical Engineering Civil Engineering 1. Elementary Fluid Dynamics 2. Fluid Kinematics 3. Finite Control Volume Analysis 4. Differential Analysis of Fluid Flow 5. Viscous
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
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 informationMicro and Macro in the Dynamics of Dilute Polymer Solutions
Micro and Macro in the Dynamics of Dilute Polymer Solutions Ravi Prakash Jagadeeshan Complex Fluid Mechanics Closed form equations Constitutive Equations Stress Calculator Simulations (BDS etc) Homogeneous
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 informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More information1 Conduction Heat Transfer
Eng6901 - Formula Sheet 3 (December 1, 2015) 1 1 Conduction Heat Transfer 1.1 Cartesian Co-ordinates q x = q xa x = ka x dt dx R th = L ka 2 T x 2 + 2 T y 2 + 2 T z 2 + q k = 1 T α t T (x) plane wall of
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
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 informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
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 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 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 informationViscosity and Polymer Melt Flow. Rheology-Processing / Chapter 2 1
Viscosity and Polymer Melt Flow Rheology-Processing / Chapter 2 1 Viscosity: a fluid property resistance to flow (a more technical definition resistance to shearing) Remember that: τ μ du dy shear stress
More informationSemana da Escola de Engenharia October 20-26, 2011 ON THE EFFECT OF THE WALL SLIP BOUNDARY CONDITON
October 2-26, 2 ON THE EFFECT OF THE WALL SLIP BOUNDARY CONDITON Luis L. Ferrás, J. Miguel Nóbrega, Fernando T. Pinho 2 and Olga S. Carneiro IPC/I3N Institute for Polymers and Composites, University of
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 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 informationA new numerical framework to simulate viscoelastic free-surface flows with the finitevolume
Journal of Physics: Conference Series PAPER OPEN ACCESS A new numerical framework to simulate viscoelastic free-surface flows with the finitevolume method Related content - Gravitational collapse and topology
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 informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationLecture 2: Constitutive Relations
Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt
More informationNumerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
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 informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More informationChapter 1. Viscosity and the stress (momentum flux) tensor
Chapter 1. Viscosity and the stress (momentum flux) tensor Viscosity and the Mechanisms of Momentum Transport 1.1 Newton s law of viscosity ( molecular momentum transport) 1.2 Generalization of Newton
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2011 Outline Outline Part I First Lecture Connection between time and ensemble average Ergodicity1 Ergodicity
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
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 informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationAuthor's personal copy
J. Non-Newtonian Fluid Mech. 139 (2006) 1 20 Numerical investigation of the velocity overshoots in the flow of viscoelastic fluids inside a smooth contraction A. Afonso a, F.T. Pinho a,b, a Centro de Estudos
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav
More informationMacromolecular Hydrodynamics Quiz Solutions. (i) To start, we recognize the following relationships on the stress and strain
Question 1 i To start, we recognize the following relationships on the stress and strain γ = γ k + γ 2 1 τ = G k γ k + μ k γ k = μ 2 γ 2 Therefore, the following relationships are also true γ = γ k + γ
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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationL8. Basic concepts of stress and equilibrium
L8. Basic concepts of stress and equilibrium Duggafrågor 1) Show that the stress (considered as a second order tensor) can be represented in terms of the eigenbases m i n i n i. Make the geometrical representation
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 informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
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 informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationMicro-Scale Gas Transport Modeling
Micro-Scale Gas Transport Modeling Continuum & Slip Flow Regimes: Navier-Stokes Equations Slip Boundary Conditions U g U w = ( σ ) σ U Kn n Slip, Transitional & Free Molecular: Direct Simulation Monte
More informationMultiscale Hydrodynamic Phenomena
M2, Fluid mechanics 2014/2015 Friday, December 5th, 2014 Multiscale Hydrodynamic Phenomena Part I. : 90 minutes, NO documents 1. Quick Questions In few words : 1.1 What is dominant balance? 1.2 What is
More informationChapter 7. Kinematics. 7.1 Tensor fields
Chapter 7 Kinematics 7.1 Tensor fields In fluid mechanics, the fluid flow is described in terms of vector fields or tensor fields such as velocity, stress, pressure, etc. It is important, at the outset,
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
More informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
More informationScale invariant fluid dynamics for the dilute Fermi gas at unitarity
Scale invariant fluid dynamics for the dilute Fermi gas at unitarity Thomas Schaefer North Carolina State University Fluids: Gases, Liquids, Plasmas,... Hydrodynamics: Long-wavelength, low-frequency dynamics
More informationZPG TBLs - Results PhD defense - Kapil Chauhan Preliminaries ZPG TBLs PG TBLs
ZPG TBLs - Results PhD defense - Kapil Chauhan Preliminaries ZPG TBLs PG TBLs ZPG TBLs - Results PhD defense - Kapil Chauhan Preliminaries ZPG TBLs PG TBLs Logarithmic composite profile Compared with past
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
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 informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More information0 = p. 2 x + 2 w. z +ν w
Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p
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 informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More information