SPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology
|
|
- Claire Lane
- 5 years ago
- Views:
Transcription
1 SPH for the Modeling of non-newtonian Fluids with Thermal-Dependent Rheology G. Bilotta 1,2 1 Dipartimento di Matematica e Informatica, Università di Catania, Italy 2 Istituto Nazionale Geofisica e Vulcanologia, Sezione di Catania, Italy BAW/RUB workshop on Applied SPH
2 Outline 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
3 Motivations 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
4 Motivations
5 Motivations
6 Motivations Work started as part of LAVA-V3 project (INGV / DPC). Aim: develop a new numerical model for lava flows capable of modeling all aspects of the flow, including: solidification crust/front/levee formation and rupture tunnel formation ephemeral vent opening
7 People Alexis Hérault modeling and GPU implementation Giuseppe Bilotta modeling and GPU implementation Eugenio Rustico multi-gpu implementation Ciro Del Negro project coordiator
8 The problem 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
9 The problem Navier-Stokes equations: Dρ Dt = ρ u, eqn. of state or u = 0, ρdu Dt = P + τ + F, where ρ, u, P, F: density, velocity, pressure, external forces (gravity) τ = µ u: viscous stress
10 The problem Navier-Stokes and heat equations: Dρ Dt = ρ u, eqn. of state or u = 0, where ρdu Dt = P + τ + F, DT ρc p = (κ T) + τ u Dt ρ, u, P, F: density, velocity, pressure, external forces (gravity) τ = µ u: viscous stress c p, κ: specific heat capacity and thermal conductivity
11 Non-Newtonian Fluids with Thermal-Dependent Rheology Thermal depenency: heat equation has a viscous heating term: i,j τ ij u i / x j ; (Costa & Macedonio, 2003) µ = µ(t). Non-Newtonian fluids: τ µ u. (or rather, µ is not independent of u).
12 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
13 Rheology Rheology: study of flow of deformable matter when subject to stress. Also: the law with which deformable matter flows, i.e. the relation between shear stress and rate of strain. Newtonian rheology: linear relation τ = µ u, coefficient µ (viscosity) does not depend on u. (Or conversely: u = φτ, with φ = 1/µ fluidity, independent from τ.) Most fluids are actually non-newtonian: stress/strain relation depends on the stress/strain rate, but we can always talk about the apparent viscosity (ratio of stress to strain rate).
14 Non-Newtonian rheologies Categories: shear-thinning or pseudo-plastic: (apparent) viscosity decreases with stress/shear rate (e.g. paint); shear-thickening or dilatant: (apparent) viscosity increases with stress/shear rate (e.g. cornstarch in water); Can be time-dependent (shear thinning: thixotropic; shear thickening: rheopectic) or time-independent: Bingham plastic, power-law rheologies, Herschel-Bulkley fluids, etc.
15 Beyond Newton One-dimensional formulations: Bingham plastic : τ = τ 0 + Kγ, no flow if τ < τ 0, Newton-like behavior afterwards with τ stress, γ shear rate, τ 0 limit stress, K consistency index, n power index.
16 Beyond Newton One-dimensional formulations: Bingham plastic : τ = τ 0 + Kγ, no flow if τ < τ 0, Newton-like behavior afterwards Power law : τ = Kγ n ; Newton for n = 1, pseudo-plastic for n > 1, dilatant for n < 1; with τ stress, γ shear rate, τ 0 limit stress, K consistency index, n power index.
17 Beyond Newton One-dimensional formulations: Bingham plastic : τ = τ 0 + Kγ, no flow if τ < τ 0, Newton-like behavior afterwards Power law : τ = Kγ n ; Newton for n = 1, pseudo-plastic for n > 1, dilatant for n < 1; Herschel-Bulkley : τ = τ 0 + Kγ n, general case that includes Newtonian, Bingham and power-law behavior; with τ stress, γ shear rate, τ 0 limit stress, K consistency index, n power index.
18 Rheology plots/i shear stress more viscous less viscous strain rate
19 Rheology plots/ii strain rate less viscous more viscous shear stress
20 Rheology plots/iii dilatant strain rate newtonian pseudoplastic Bingham plastic Herschel-Bulkley pseudoplastic shear stress
21 Herschel-Bulkley General formulation, with temperature dependency: { µ µ(t) = 0 (T) if ɛ ɛ 0 (T), K(T)ɛ n(t) + τ 0(T) ɛ if ɛ > ɛ 0 (T) (ɛ: second invariant of rate of strain tensor) infinite viscosity for null strain
22 Herschel-Bulkley General formulation, with temperature dependency, regularized (Zhu, following Papanastasiou): { µ µ(t) = 0 (T) if ɛ ɛ 0 (T), K(T)ɛ n(t) + τ 0 (T) 1 exp( αɛ) ɛ if ɛ > ɛ 0 (T) (ɛ: second invariant of rate of strain tensor) (α some large regularizing constant)
23 SPH discretization 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
24 Smoothed Particle Hydrodynamics Smoothed Particle Hydrodynamics. Lagrangian meshless method with smoothed field (position, velocity, temperature) interpolation: f (x) = f (y)δ(x y)dy f (y)w(x y, h)dy Ω Ω j f (x j )W(x x j, h)v j where W is a smoothing kernel that approximates the Dirac s δ, V j is the volume of particle j, kernel with compact support (radius h) = summation extended to neighbors only.
25 SPH kernel properties E.g. density: ρ(x) j ρ(x j )W(x x j, h)v j = j m j W(x x j, h) where m j = V j ρ(x j ) is the mass of particle j. Kernel normalization (Shepard): W(x x i, h) = W(x x i ) j W(x x j)v j allows exact reconstruction of constant fields: normalized kernels ensure zero-order accuracy, symmetric kernels ensure first-order accuracy, positive kernels ensure isolated particles have nonnegative density.
26 SPH gradients/i For gradients, we use Green s formula to get f ( x) = f ( y)δ( x y)d y f ( y)w( x y, h)d y Ω Ω f ( y) y W( x y, h)d y + Ω + f ( y)w( x y, h) nds j Ω m j ρ j f j j W( x x j, h) +??? = = j m j ρ j f j i W( x x j, h) +???. If supp W i Ω, then??? = 0. Otherwise, correction terms have to be added.
27 SPH gradients/ii Other possible formulas: f (x i ) m j (f ρ j ± f i ) i W(x i x j, h), j j f (x i ) 1 m ρ j (f i f j ) i W(x i x j, h), i f (x i ) ρ i j j m j ( fj ρ 2 j + f i ρ 2 i ) i W(x i x j, h). Chosen for different fields for different stability properties. Common property: explicit computation, no need to invert linear systems.
28 SPH discretization/i Discretized continuity equations (W ij = W(x i x j, h), f ij = f i f j ): ( ) dρ = m dt j u ij i W ij or i or (summation density): ( ) dρ m j = ρ dt i u i ρ ij i W ij, j ρ i = j m j W ij. (but summation density can only be used reliably in absence of a free surface).
29 SPH discretization/ii Discretized momentum and thermal equations (W ij = W(x i x j, h), f ij = f i f j ): ( ) du = ( Pj m dt j i ρ 2 j + m j µ i + µ j ρ i ρ j ( ) dt c pi = 4m j dt i ρ i ρ j 2m j ρ i ρ j + P i ρ 2 i ) i W ij + r ij i W ij u ij + g, r 2 ij κ i κ j r ij i W ij T κ i + κ ij j rij 2 + µ i µ j µ i + µ j (r ij u ij ) 2 r ij i W ij r 4 ij
30 Weakly compressible SPH To preserve direct computation, slightly compressible flows with equation of state P = ρ (( ) 0c s ρ γ 1) γ (ρ 0 density at rest, c s (artificial) speed of sound, taken an order of magnitude larger than the maximum expected speed of the phenomenon). Speed of sound is chosen ficticious because it limits time-stepping in explicit methods. (For lava we don t really worry about that, because the viscosity is actually what limits the timestep.) ρ 0
31 Viscosity Non-Newtonian fuid: µ i is the apparent viscosity of the particle, computed separately before each force computation. Strain rate tensor: u i = j u ij i W ij V j Second invariant: ɛ i = 1 ( (tr u) 2 2 tr( u 2 ) ) Then compute µ i from (regularized) Herschel-Bulkley
32 Boundary conditions: ground/i Lava ows on complex topographies, we assume them de ned through a Digital Elevation Model (DEM): 2D grid of elevations, cell side p. Interpolate to nd z0 height at projection (x0, y0) of particle on DEM, again to nd height at (x0 + p, y0), (x0, y0 + p), nd plane by three points, interact with plane.!y P(x,y,z)!!x " z0!y P(x,y,z)!!x " z0 z2 z1
33 Boundary conditions: ground/ii! n h! u r! u T Friction:!! u! f = µs 1 u T f = µ T! n! " ( h2 # r 2 ) with S 1 ground particle surface. Ground particle surface can be estimated by S 1 = S g /N g or set as problem parameter.
34 Boundary conditions: ground/iii Conduction to ground: particle close to the DEM lose T = κs 1 c p m T T ground t r due to heat exchange with ground. Ground temperature T ground is assumed constant.
35 Boundary conditions: free surface/i Need to detect free surface for radiative losses. Standard approach: check for neighbors in cone pointing towards gradient of characteristic function. (a) (b) Trick: lava flows are essentially two-dimensional, check for presence of particles in upwards cone. (a is on the surface, b is not).
36 Boundary conditions: free surface/i Radiative losses: T = K BεS 2 c p m (T 4 T 4 ext ) T. with K B Stefan-Boltzmann constant, ε emissivity. External temperature T ext is assumed constant. Free-surface particle surface S 2 is either a problem parameter, or can be computed as S 2 = S f /N f.
37 Phase change q: fraction of latent heat for phase change; q = 0: liquid, q = 1: solid; T : new temperature from integration of heat equation; if T < T s and q < 1, then: T = T s increment q by c p (T s T )/L if q 1, then: T = T s L(q 1)/c p q = 1 otherwise T = T.
38 Implementation/I Weakly compressible SPH is embarrasingly parallel (each particle can be computed concurrently and independently from the others, by just reading their data).
39 Implementation/II SPH: excellent candidate for GPU implementation, > 150 faster than (serial) CPU, 50 faster than parallel CPU implementation. GPUs also give us free interpolation (in hardware), good for topography.
40 Implementation/III Based on GPUSPH, open source implementation of SPH on GPU using CUDA, available from implements SPH for Newtonian fluid-dynamics multiple SPH formulations variety of smoothing kernel choices variety of viscosity models variety of boundary conditions single GPU, multi-gpu, multi-node
41 Lava flow simulations 1 Motivations 2 The problem 3 Non-Newtonian Rheologies 4 SPH discretization 5 Lava flow simulations 6 Conclusions
42 Parameters parameter symbol values unit density ρ 2400 kg m 3 heat capacity C p 1150 J kg 1 K 1 latent heat L J kg 1 conductivity κ 2.2 W m 1 K 1 emissivity ε 1 vent temperature T e 1350 K solidus temperature T s 1163 K ground temperature T ground 293 K air temperature T ext 293 K
43 Temperature-dependent parameters Consistency factor from Giordano & Dingwell (2003): log K(T) = [H 2O] T ln([h 2 O]) with T in Kelvin and [H 2 O] is the water content in weight percentage. In our tests [H 2 O] = 0.005%. Yield strength from Miyamoto (1997): log τ 0 (T) = (T T e ) Regularizing coefficient α = 1000, and for Herschel-Bulkley: n(t) = 0.1 (shear-thinning).
44 Apparent viscosities K 1263 K K 1263 K Bingham Herschel-Bulkley
45 Newtonian Constant viscosity Temperature-dependent viscosity
46 Bingham Constant viscosity Temperature-dependent viscosity
47 Herschel-Bulkley Constant viscosity Temperature-dependent viscosity
48 Conclusions we can simulate both Newtonian and non-newtonian fluids with SPH; coupled thermal/dynamic problems are possible, including phase transition and thermal dependent rheology; GPU implementation gives us good performance at low cost, and some things for free (e.g. DEM interpolation); our model can be used to explore non-newtonian fluids, and help finding what is the most correct rheology for lava.
Computational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationMECHANICAL PROPERTIES
MECHANICAL PROPERTIES Rheology S.C. BAYNE, 1 J.Y. Thompson 2 1 University of Michigan School of Dentistry, Ann Arbor, MI 48109-1078 sbayne@umich.edu 2 Nova Southeastern College of Dental Medicine, Ft.
More informationRheology and Constitutive Equations. Rheology = Greek verb to flow. Rheology is the study of the flow and deformation of materials.
Rheology and Constitutive Equations Rheology = Greek verb to flow Rheology is the study of the flow and deformation of materials. The focus of rheology is primarily on the study of fundamental, or constitutive,
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationPreliminary validation of lava benchmark tests on the GPUSPH particle engine
ACCEPTED ON ANNALS OF GEOPHYSICS, 61, 2018; Doi: 10.4401/ag-7870 Preliminary validation of lava benchmark tests on the GPUSPH particle engine V. Zago a,c, G. Bilotta a, A. Cappello a, R. A. Dalrymple b,
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More informationPlease remember all the unit that you use in your calculation. There are no marks for correct answer without unit.
CHAPTER 1 : PROPERTIES OF FLUIDS What is fluid? A fluid is defined as a substance that deforms continuously when acted on by a shearing stress at any magnitude. In a fluid at rest, normal stress is called
More informationChapter 3 Non-Newtonian fluid
Chapter 3 Non-Newtonian fluid 3-1. Introduction: The study of the deformation of flowing fluids is called rheology; the rheological behavior of various fluids is sketchen Figure 3-1. Newtonian fluids,
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 information1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts)
1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts) Introduction, general information Course overview Fluids as a continuum Density Compressibility Viscosity Exercises: A1 Fluid mechanics Fluid
More informationThis chapter is a study of the shear stress as a function of the shear rate for Newtonian and non-newtonian biological materials.
Lecture 10 and Flow (Ch. 6) This chapter is a study of the shear stress as a function of the shear rate for Newtonian and non-newtonian biological materials. 1 Lecture 10 and Flow (Ch. 6) When a fluid
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries Section 4: Non-Newtonian fluids and rheometry (PART 1)
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 4: Non-Newtonian fluids and rheometry (PART 1) Globex Julmester 2017 Lecture #3 05 July 2017 Agenda Lecture #3 Section
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationLecture 3. Properties of Fluids 11/01/2017. There are thermodynamic properties of fluids like:
11/01/2017 Lecture 3 Properties of Fluids There are thermodynamic properties of fluids like: Pressure, p (N/m 2 ) or [ML -1 T -2 ], Density, ρ (kg/m 3 ) or [ML -3 ], Specific weight, γ = ρg (N/m 3 ) or
More informationViscous dissipation and temperature dependent shear thinning rheology
Viscous dissipation and temperature dependent shear thinning rheology M. Filippucci 1, A. Tallarico 1, M. Dragoni 2 1 Dipartimento di Scienze della Terra e Geoambientali, Università degli Studi di Bari
More informationRheological Properties
Rheological Properties Shear Stress the force experienced by the liquid is proportional to the area of the upper plate and is defined as the shear stress, σ = F/A (Units = Nm -2 or Pascals [Pa]) Shear
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationSIMULATIONS OF THE 2004 LAVA FLOW AT ETNA VOLCANO BY THE MAGFLOW CELLULAR AUTOMATA MODEL
SIMULATIONS OF THE 2004 LAVA FLOW AT ETNA VOLCANO BY THE MAGFLOW CELLULAR AUTOMATA MODEL Ciro Del Negro 1, Luigi Fortuna 2, Alexis Herault 1,3, Annamaria Vicari 1 1 Istituto Nazionale di Geofisica e Vulcanologia
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 informationDuality methods for variational inequalities and Non-Newtonian fluid mechanics
Duality methods for variational inequalities and Non-Newtonian fluid mechanics Enrique Fernández-Nieto, Paul Vigneaux Dpto. Matemática Aplicada I, Universidad de Sevilla UMPA, Ecole Normale Supérieure
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 informationCHAPTER 1 Fluids and their Properties
FLUID MECHANICS Gaza CHAPTER 1 Fluids and their Properties Dr. Khalil Mahmoud ALASTAL Objectives of this Chapter: Define the nature of a fluid. Show where fluid mechanics concepts are common with those
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
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 informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationOn a variational inequality of Bingham and Navier-Stokes type in three dimension
PDEs for multiphase ADvanced MATerials Palazzone, Cortona (Arezzo), Italy, September 17-21, 2012 On a variational inequality of Bingham and Navier-Stokes type in three dimension Takeshi FUKAO Kyoto University
More informationModelling the Rheology of Semi-Concentrated Polymeric Composites
THALES Project No 1188 Modelling the Rheology of Semi-Concentrated Polymeric Composites Research Team Evan Mitsoulis (PI), Professor, NTUA, Greece Costas Papoulias (Research Student), NTUA, Greece Souzanna
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More information1. The Properties of Fluids
1. The Properties of Fluids [This material relates predominantly to modules ELP034, ELP035] 1.1 Fluids 1.1 Fluids 1.2 Newton s Law of Viscosity 1.3 Fluids Vs Solids 1.4 Liquids Vs Gases 1.5 Causes of viscosity
More informationChapter 6 Molten State
Chapter 6 Molten State Rheology ( 流變學 ) study of flow and deformation of (liquid) fluids constitutive (stress-strain) relation of fluids shear flow shear rate ~ dγ/dt ~ velocity gradient dv 1 = dx 1 /dt
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 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 informationPetroleum Engineering Dept. Fluid Mechanics Second Stage Dr. Ahmed K. Alshara
Continents Chapter 1. Fluid Mechanics -Properties of fluids -Density, specific gravity, specific volume and Viscosity -Newtonian and non Newtonian fluids -Surface tension Compressibility -Pressure -Cavitations
More informationProcess Development for a High-Throughput Fine Line Metallization Approach Based on Dispensing Technology
Process Development for a High-Throughput Fine Line Metallization Approach Based on Dispensing Technology Rheology, CFD-Simulation and Print-Head Design M. Pospischil, M. Klawitter, M. Kuchler, J. Specht,
More informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
More informationCSCI1950V Project 4 : Smoothed Particle Hydrodynamics
CSCI1950V Project 4 : Smoothed Particle Hydrodynamics Due Date : Midnight, Friday March 23 1 Background For this project you will implement a uid simulation using Smoothed Particle Hydrodynamics (SPH).
More informationContents. Preface XIII. 1 General Introduction 1 References 6
VII Contents Preface XIII 1 General Introduction 1 References 6 2 Interparticle Interactions and Their Combination 7 2.1 Hard-Sphere Interaction 7 2.2 Soft or Electrostatic Interaction 7 2.3 Steric Interaction
More informationKELVIN-HELMHOLTZ INSTABILITY BY SPH
II International Conference on Particle-based Methods - Fundamentals and Applications PARTICLES 2011 E. Oñate and D.R.J. Owen (Eds) KELVIN-HELMHOLTZ INSTABILITY BY SPH M. S. Shadloo, M. Yildiz Faculty
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationPharmaceutics I. Unit 6 Rheology of suspensions
Pharmaceutics I اينالديصيدلينيات 1 Unit 6 Rheology of suspensions 1 Rheology, the science of the flow or deformation of matter (liquid or soft solid) under the effect of an applied force. It addresses
More information(Refer Slide Time: 2:14)
Fluid Dynamics And Turbo Machines. Professor Dr Shamit Bakshi. Department Of Mechanical Engineering. Indian Institute Of Technology Madras. Part A. Module-1. Lecture-3. Introduction To Fluid Flow. (Refer
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 informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationAPPENDIX A USEFUL EQUATIONS (METRIC AND IMPERIAL SYSTEMS) THE DEFINITION OF VISCOSITY RHEOLOGICAL (VISCOUS BEHAVIOR) PROPERTIES OF FLUIDS
APPENDIX A USEFUL EQUATIONS (METRIC AND IMPERIAL SYSTEMS) THE DEFINITION OF VISCOSITY RHEOLOGICAL (VISCOUS BEHAVIOR) PROPERTIES OF FLUIDS A APPENDIX A APPENDIX A Flow vs. velocity Specific gravity vs.
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 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 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 informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationSmooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520
Smooth Particle Hydrodynamic (SPH) Presented by: Omid Ghasemi Fare Nina Zabihi XU Zhao Miao Zhang Sheng Zhi EGEE 520 OUTLINE Ø Introduction and Historical Perspective: Ø General Principles: Ø Governing
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
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 informationLiquid fuels viscosity (non-newtonian fluids)
Liquid fuels viscosity (non-newtonian fluids) Introduction Viscosity is an important parameter when preparing liquid fuels for combustion as it determines how much drag the fuel experiences when passing
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 informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationThe Physics of Non-Newtonian Fluids
The Physics of Non-Newtonian Fluids M. Grogan 10. April 2006 [1][4] R. Fox, A. McDonald, P. Pritchard, Introduction to Fluid Mechanics, 6 th Ed. 2004 Wiley, United States of America. p3, p27. [2],[3] B
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.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 informationPharmaceutical compounding I Colloidal and Surface-Chemical Aspects of Dosage Forms Dr. rer. nat. Rebaz H. Ali
University of Sulaimani School of Pharmacy Dept. of Pharmaceutics Pharmaceutical Compounding Pharmaceutical compounding I Colloidal and Surface-Chemical Aspects of Dosage Forms Dr. rer. nat. Rebaz H. Ali
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationNUMERICAL SIMULATION OF THE FLOW OF A POWER LAW FLUID IN AN ELBOW BEND. A Thesis KARTHIK KANAKAMEDALA
NUMERICAL SIMULATION OF THE FLOW OF A POWER LAW FLUID IN AN ELBOW BEND A Thesis by KARTHIK KANAKAMEDALA Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationME 262 BASIC FLUID MECHANICS Assistant Professor Neslihan Semerci Lecture 4. (Buoyancy and Viscosity of water)
ME 262 BASIC FLUID MECHANICS Assistant Professor Neslihan Semerci Lecture 4 (Buoyancy and Viscosity of water) 16. BUOYANCY Whenever an object is floating in a fluid or when it is completely submerged in
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationIntroduction to Viscometry and Rheology, Basics, Rotational Testing. Basic Seminar Applied Rheology
Introduction to Viscometry and Rheology, Basics, Rotational Testing Basic Seminar Applied Rheology Contents Definition of basic rheological parameters Viscosity and elasticity Deformation, shear stress
More informationLiquids and solids are essentially incompressible substances and the variation of their density with pressure is usually negligible.
Properties of Fluids Intensive properties are those that are independent of the mass of a system i.e. temperature, pressure and density. Extensive properties are those whose values depend on the size of
More informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationPetroleum Engineering Department Fluid Mechanics Second Stage Assist Prof. Dr. Ahmed K. Alshara
Continents Petroleum Engineering Department Fluid Mechanics Second Stage Assist Prof. Dr. Ahmed K. Alshara Chapter 1. Fluid Mechanics -Properties of fluids -Density, specific gravity, specific volume and
More informationPolymerization Technology Laboratory Course
Polymerization Technology Laboratory Course Viscometry/Rheometry Tasks 1. Comparison of the flow behavior of polystyrene- solution and dispersion systems 2. Determination of the flow behaviour of polyvinylalcohol
More informationPart II Fundamentals of Fluid Mechanics By Munson, Young, and Okiishi
Part II Fundamentals of Fluid Mechanics By Munson, Young, and Okiishi WHAT we will learn I. Characterization of Fluids - What is the fluid? (Physical properties of Fluid) II. Behavior of fluids - Fluid
More informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
More informationFluid Mechanics Abdusselam Altunkaynak
Fluid Mechanics Abdusselam Altunkaynak 1. Unit systems 1.1 Introduction Natural events are independent on units. The unit to be used in a certain variable is related to the advantage that we get from it.
More informationRheometer: Procedure: Part A: Viscosity v Time
Rheometer A fluid is defined as a substance that deforms continuously under the action of a shear stress, no matter how small the shear stress may be. Without shear stress, there will be no deformation.
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More informationCH5716 Processing of Materials
CH5716 Processing of Materials Ceramic Thick Film Processing Lecture MC5 Slurry Characterisation Specific Surface Area Powder size & specific surface area (area per unit wt) closely related As particle
More informationIntroduction to Geology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.001 Introduction to Geology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. RHEOLOGICAL MODELS Rheology
More informationCHAPTER (2) FLUID PROPERTIES SUMMARY DR. MUNZER EBAID MECH.ENG.DEPT.
CHAPTER () SUMMARY DR. MUNZER EBAID MECH.ENG.DEPT. 08/1/010 DR.MUNZER EBAID 1 System Is defined as a given quantity of matter. Extensive Property Can be identified when it is Dependent on the total mass
More informationPhysics-Based Animation
CSCI 5980/8980: Special Topics in Computer Science Physics-Based Animation 13 Fluid simulation with grids October 20, 2015 Today Presentation schedule Fluid simulation with grids Course feedback survey
More informationNotes 4: Differential Form of the Conservation Equations
Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationSmoothed Particle Hydrodynamics
Smoothed Particle Hydrodynamics for astronomical model demonstration (2010) Thomas A. Grønneløv Figure 1: Collision between protoplanets e-mail:tag@greenleaf.dk This text can currently (april 2010) be
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationNumerical Simulation of Newtonian and Non-Newtonian Flows in Bypass
Numerical Simulation of Newtonian and Non-Newtonian Flows in Bypass Vladimír Prokop, Karel Kozel Czech Technical University Faculty of Mechanical Engineering Department of Technical Mathematics Vladimír
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationFUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL
Annual Journal of Hydraulic Engineering, JSCE, Vol.54, 2010, February FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL How Tion PUAY1 and Takashi HOSODA2 1 Member of JSCE, Phd Student,
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationSmoothed Particle Hydrodynamics (SPH) 4. May 2012
Smoothed Particle Hydrodynamics (SPH) 4. May 2012 Calculating density SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationSmoothed Particle Hydrodynamics (SPH) Huamin Wang
Smoothed Particle Hydrodynamics (SPH) Huamin Wang Fluid Representation Fluid is represented using a set of particles. Particle i has position x i, velocity v i, and mass m i. Animation Example Using 10
More informationWe may have a general idea that a solid is hard and a fluid is soft. This is not satisfactory from
Chapter 1. Introduction 1.1 Some Characteristics of Fluids We may have a general idea that a solid is hard and a fluid is soft. This is not satisfactory from scientific or engineering point of view. In
More information