Viscoelasticity in mantle convection
|
|
- Deborah Adams
- 5 years ago
- Views:
Transcription
1 Mgr. et Mgr. Vojtěch Patočka Supervised by: RNDr. Ondřej Čadek, CSc. 5th June 2015
2 Content Method of implementation Model Method Testing the method Thermal convection Mechanical deformation Results Numerically instaneous behaviour Time step decoupling Elastic episodes Construction of the constitutve equation 3D analogs Entropy production maximization
3 I2ELVIS - Taras Gerya staggered grid, markers regional studies, benchmarks 3D version not elastic yet - I3VIS Ellipsis - Moresi, Muhlhaus finite elements, markers used as integration points regional studies StagYY - Paul Tackley staggered grid, global studies, 3D, spherical shell not elastic yet
4 Model Extended Boussinesq - nondimensionalised equation of motion + continuity + thermal equation rheology p + σ = Ra T e z v = 0 T t = 2 T v T σ + D σ = 2ηD η(p, T ) = exp( A(T T ref ) + B(p p ref ))
5 Method Constitutive equation σ + D σ = 2ηD D(p, T ) := τ rel κ h 2 = ηκ Gh 2 σ := σ + ( v ) σ + (σw Wσ) a(dσ + σd) t advection co-rotation co-deformation Discretization (mixed euler, 1st order accurate) σ n = dt dt + D 2ηDn + D dt + D σn 1 + D dt dt + D (Wσ σw+a(dσ+σd))
6 Method σ n = dt dt + D 2ηDn + D ( ) σ n 1 + dt (Wσ σw + a (Dσ + σd)) dt + D σ n 1 := σ n 1 + dt ( Wσ σw + a ( Dσ + σd ) ) n 1/n σ n = 2η n num D n + Dn dt + D n σn 1 ; η num := dt dt + D η Viscoelasticity parameter Z Z := dt dt + D ; σn = 2 Z η n D n + (1 Z) σ n 1
7 Method dt n 1, n = stability( v n 1, D, d); dt a ; dt el T n = advection-diffusion (T n 1, v n 1, dt a ) σ n 1 = σ n 1 + deadvection (σ n 1, v n 1, dt a ) η n = η(t n ); η n num = η n num(η n, dt el, D n ) rhs n = ( D n dt el +D n σ n 1 ) Ra T n e z no decoupling: dt a = dt n 1, n dt el = dt n 1, n Stokes(η n num, v n, p n ) = rhs n σ n = σ n ( σ n 1, η n num, v n )
8 Thermal convection Harder - Journal of Non-Newtonian Fluid Mechanics, 39 (1991) 67-88
9 Thermal convection Co-rotational derivative (Jaumann, a = 0)
10 Thermal convection Upper convected derivative (a = 1)
11 Thermal convection Viscous case
12 Mechanical deformation Bending of an elastic slab (under gravity for 20 Kyr) ρ 1 = 4000; η 1 = 1 e27; G 1 = 1 e10; ρ 2 = 1; η 2 = 1 e21; G 2 = 1 e20;
13 Numerically instaneous behaviour Onset of convection is related to numerically instaneous elastic behaviour for Deborah number above certain critical value. Ra = 10 4 D critical =
14 Numerically instaneous behaviour Example of velocity singularity
15 Numerically instaneous behaviour Jaumann - parametric study, failed runs
16 Time step decoupling dt n 1, n = stability( v n 1, D, d); dt a ; dt el T n = advection-diffusion (T n 1, v n 1, dt a ) σ n 1 = σ n 1 + deadvection (σ n 1, v n 1, dt a ) η n = η(t n ); η n num = η n num(η n, dt el, D n ) Gerya / Moresi decoupling: multiple advection rhs n = ( D n dt el +D n σ n 1 ) Ra T n e z Stokes(η n num, v n, p n ) = rhs n σ n = σ n ( σ n 1, η n num, v n )
17 Time step decoupling dt n 1, n = stability( v n 1, D, d); dt a ; dt el T n = advection-diffusion (T n 1, v n 1, dt a ) σ n 1 = σ n 1 + deadvection (σ n 1, v n 1, dt a ) η n = η(t n ); η n num = η n num(η n, dt el, D n ) rhs n = ( D n dt el +D n σ n 1 ) Ra T n e z convection decoupling: interpolation needed Stokes(η n num, v n, p n ) = rhs n σ n = σ n ( σ n 1, η n num, v n )
18 Time step decoupling Realistic case: separate viscous and viscoelastic parts of the domain Z := dtel dt el + D ; σn = 2 Z η n D n + (1 Z) σ n 1 bulid-up relaxation Desired behaviour: viscoelastic material does not move across many cells within elastic time step (viscous material can - elastic time step does not play any role in viscous part of the domain). D = ηd 0 ; dt el = 0.01 τ rel κ h 2 = 0.01 D 0; Z = η ; 0.01
19 Time step decoupling Jaumann - failed runs, decoupled
20 Time step decoupling Upper convected - failed runs, decoupled
21 Time step decoupling Jaumann - not trusted runs, decoupled
22 Time step decoupling Critical Deborah, time step decoupled
23 Time step decoupling Critical Deborah, time step decoupled: elastic strain rate
24 Time step decoupling Critical Deborah, time step decoupled: upper convected
25 Elastic episodes
26 Elastic episodes Undercritical Deborah number - time step/scheme convergence
27 Elastic episodes Overcritical Deborah number - time step convergence
28 Elastic episodes Overcritical Deborah number - time scheme convergence
29 Elastic episodes Time step analysis - overcritical Deborah number
30 Four ways to derive our rheological equation 1D analogs objective rate not specified Entropy production maximization dissipation needs to be prescribed objective rate specified 3D analogs - is it new? objective rate specified Microscopical model e.g. dumbbells diluted in newtonian fluid We need to know how the equation was derived for interpretation: elastic strain rate, elastic energy, dissipation...
31 3D analogs 3D analogs Geodynamical texts often use σ ij = 2η ε ij - incorrect Four possible strain tensors: eulerian/lagrangian, full/linearized ε lin := 1 2 ( u + ( u)t ); ξ lin := 1 2 (GradU + (GradU)T ) ε := 1 2 (I F T F 1 ); ξ := 1 2 (FT F I) which have following material time derivatives ε lin = D 1 2 (( u)l + LT ( u) T ) ε = D (εl + L T ε) ξ lin = D 1 2 (( u)t L T + L( u)) ξ = F T D F
32 3D analogs Maxwell rheology is based on this idea: Elastic and viscous deformations add together, but stress is the same for both elastic and viscous parts of the deformation τ = p + σ; D = D vis + D el σ = 2ηD vis = 2Gε el Choice of ε or ξ or ε lin or ξ lin gives us the relation ε el D el and thus specifies the rate used in our rheology (only for ε we get an objective rate - lower convected one)
33 Rajagopal Entropy production maximization Motivation - interpretation of physical quantities: what is dissipation, elastic strain rate, elastic energy... σ : v = σ : D = σ : σ 2η + σ : σ 2G Dissipation? Rate of change of elastic energy?
34 Rajagopal Oldroyd B model - ICC Thermodynamic considerations e(δ, ρ, B Kp(t)) = e 0 (δ, ρ) + G 2ρ (trb K p(t) 3 log det B Kp(t)) ζ = (T G(B Kp(t) I)) d : D d + G(C Kp(t) I) : D Kp(t) Prescribed dissipation Resulting model ζ = 2η 2 D 2 + 2ηD Kp(t) C Kp(t) : D Kp(t) T = pi + 2η 2 D + G(B Kp(t) I) B Kp(t) = G η (B K p(t) I)
35 Rajagopal After setting η 2 = 0, we get Maxwell model T = pi + G(B Kp(t) I) B Kp(t) = G η (B K p(t) I) T = pi + σ σ + τ rel σ = 2ηD But now we now the dissipation of our model, we prescribed it So the question is ζ = 2ηD Kp(t) C Kp(t) : D Kp(t) 2ηD Kp(t) C Kp(t) : D Kp(t)? = σ : σ 2η
36 Rajagopal By using σ = 2ηD τ rel σ and σ = G(BKp(t) I), this question reduces to (kinematic consideration regarding the natural configuration needed): (F Kp(t) D Kp(t) F T K p(t) ) : (F K p(t) D Kp(t) F T K p(t) )? = which is clearly not satisfied? = (D Kp(t) F T K p(t) F K p(t)) : D Kp(t)
37 Rajagopal Conclusions We used a grid based method to implement viscoelastic rheology into a simple thermal convection model Time step decoupling which could be trusted in realistic cases was proposed Simulations show that locally high stresses are produced by viscoelasticity, which means that viscoelasticity could trigger plasticity in visco-elasto-plastic code
38 Rajagopal Thank you
Vít Průša (joint work with K.R. Rajagopal) 30 October Mathematical Institute, Charles University
On models for viscoelastic fluid-like materials that are mechanically incompressible and thermally compressible or expansible and their Oberbeck Boussinesq type approximations Vít Průša (joint work with
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 informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
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 informationRobust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems
Physics of the Earth and Planetary Interiors 163 (2007) 83 105 Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems Taras V. Gerya a David A. Yuen b a
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 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 informationSubduction II Fundamentals of Mantle Dynamics
Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011 Rheology Elasticity vs. viscous deformation η
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
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 informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
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 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 informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More informationThree-dimensional numerical simulations of thermo-chemical multiphase convection in Earth s mantle Takashi Nakagawa a, Paul J.
Three-dimensional numerical simulations of thermo-chemical multiphase convection in Earth s mantle Takashi Nakagawa a, Paul J. Tackley b a Department of Earth and Planetary Sciences, University of Tokyo,
More informationOn the analysis for a class of thermodynamically compatible viscoelastic fluids with stress diffusion
On the analysis for a class of thermodynamically compatible viscoelastic fluids with stress diffusion Josef Málek Nečas Center for Mathematical Modeling and Charles University, Faculty of Mathematics and
More informationThe finite difference code (fully staggered grid) includes a marker-in-cell Lagrangian marker
GSA Data Repository 2018289 Ruh et al., 2018, Shale-related minibasins atop a massive olistostrome in an active accretionary wedge setting: Two-dimensional numerical modeling applied to the Iranian Makran:
More informationContinuum Mechanics and Theory of Materials
Peter Haupt Continuum Mechanics and Theory of Materials Translated from German by Joan A. Kurth Second Edition With 91 Figures, Springer Contents Introduction 1 1 Kinematics 7 1. 1 Material Bodies / 7
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 informationEffects of elasticity on the Rayleigh Taylor instability: implications for large-scale geodynamics
Geophys. J. Int. (2006) doi: 10.1111/j.1365-246X.2006.03201.x Effects of elasticity on the Rayleigh Taylor instability: implications for large-scale geodynamics Boris J. P. Kaus and Thorsten W. Becker
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 informationMechanical properties of polymers: an overview. Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore
Mechanical properties of polymers: an overview Suryasarathi Bose Dept. of Materials Engineering, IISc, Bangalore UGC-NRCM Summer School on Mechanical Property Characterization- June 2012 Overview of polymer
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 informationSimulations of viscoelastic droplet deformation through a microfluidic contraction
Advances in Fluid Mechanics VI 69 Simulations of viscoelastic droplet deformation through a microfluidic contraction D. J. E. Harvie 1,M.R.Davidson 1 & J. J. Cooper-White 2 1 Department of Chemical and
More informationA Three-Dimensional Anisotropic Viscoelastic Generalized Maxwell Model for Ageing Wood Composites
A Three-Dimensional Anisotropic Viscoelastic Generalized for Ageing Wood Composites J. Deteix, G. Djoumna, A. Fortin, A. Cloutier and P. Blanchet Society of Wood Science and Technology, 51st Annual Convention
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 informationarxiv: v1 [cs.na] 21 Jul 2010
DEVELOPMENT OF THREE DIMENSIONAL CONSTITUTIVE THEORIES BASED ON LOWER DIMENSIONAL EXPERIMENTAL DATA Satish Karra, K. R. Rajagopal, College Station arxiv:1007.3760v1 [cs.na] 1 Jul 010 Abstract. Most three
More informationComputational Modelling of Mechanics and Transport in Growing Tissue
Computational Modelling of Mechanics and Transport in Growing Tissue H. Narayanan, K. Garikipati, E. M. Arruda & K. Grosh University of Michigan Eighth U.S. National Congress on Computational Mechanics
More informationMicrostructural studies for rheology. Chapter 7: Microstructural studies for rheology. Separation of length scales. Micro & macro views
Chapter 7: Microstructural studies for rheology Microstructural studies for rheology To calculate the flow of complex fluids, need governing equations, in particular, the constitutive equation relating
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 informationBenchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003
Introduction Benchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003 In early October 2002 a group of researchers met at the University of Michigan at Ann Arbor,
More information1 Scaling analysis and non-dimensional numbers
1 Scaling analysis and non-dimensional numbers While this is a textbook on numerical analysis, it is crucial to keep the nature of the physical processes which we would like to model in mind. This will
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 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 informationInitiation of localized shear zones in viscoelastoplastic rocks
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb003652, 2006 Initiation of localized shear zones in viscoelastoplastic rocks Boris J. P. Kaus 1,2 and Yuri Y. Podladchikov 3 Received 28 January
More information1 2D Stokes equations on a staggered grid using primitive variables
Figure 1: Staggered grid definition. Properties such as viscosity and density inside a control volume (gray) are assumed to be constant. Moreover, a constant grid spacing in x and -direction is assumed.
More informations ij = 2ηD ij σ 11 +σ 22 +σ 33 σ 22 +σ 33 +σ 11
Plasticity http://imechanicaorg/node/17162 Z Suo NONLINEAR VISCOSITY Linear isotropic incompressible viscous fluid A fluid deforms in a homogeneous state with stress σ ij and rate of deformation The mean
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary information: Our 1080 km x 360 km model setup includes a 7 km thick oceanic crust adjacent to a 60 km thick, 250 km wide orogenic crust. Both the oceanic and the
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 informationarxiv: v3 [physics.flu-dyn] 13 Dec 2016
ON THERMODYNAMICS OF VISCOELASTIC RATE TYPE FLUIDS WITH TEMPERATURE DEPENDENT MATERIAL COEFFICIENTS JAROSLAV HRON, VOJTĚCH MILOŠ, VÍT PRŮŠA, ONDŘEJ SOUČEK, AND KAREL TŮMA arxiv:1612.01724v3 [physics.flu-dyn]
More informationMantle convection models with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling
Mantle convection models with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling Louis Moresi, Hans Mühlhaus, Frédéric Dufour CSIR Exploration and Mining, P Box 437, Nedlands,
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 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 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 informationA Nonlinear Generalized Standard Solid Model for Viscoelastic Materials
A Nonlinear Generalized Standard Solid Model for Viscoelastic Materials A Nonlinear Generalized Standard Solid Model for Viscoelastic Materials Marc Delphin MONSIA From: Département de Physique Université
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
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 informationInterfacial dynamics
Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic
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 informationDETAILS ABOUT THE TECHNIQUE. We use a global mantle convection model (Bunge et al., 1997) in conjunction with a
DETAILS ABOUT THE TECHNIQUE We use a global mantle convection model (Bunge et al., 1997) in conjunction with a global model of the lithosphere (Kong and Bird, 1995) to compute plate motions consistent
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 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 informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationINTRODUCTION TO NUMERICAL GEODYNAMIC MODELLING
INTRODUCTION TO NUMERICAL GEODYNAMIC MODELLING Until now, numerical modelling of geodynamic processes has been the domain of highly trained mathematicians with long experience of numerical and computational
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
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 informationCourse Business. Today: isostasy and Earth rheology, paper discussion
Course Business Today: isostasy and Earth rheology, paper discussion Next week: sea level and glacial isostatic adjustment Email did you get my email today? Class notes, website Your presentations: November
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More informationRheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep
Rheology III Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Ideal materials fall into one of the following categories:
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 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 informationFor an imposed stress history consisting of a rapidly applied step-function jump in
Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For
More informationWorkshop on Geophysical Data Analysis and Assimilation. 29 October - 3 November, Modelling Mantle Convection and Plate Tectonics
2373-16 Workshop on Geophysical Data Analysis and Assimilation 29 October - 3 November, 2012 Modelling Mantle Convection and Plate Tectonics Paul J. Tackley ETH Zuerich Switzerland Modelling Mantle Convection
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 informationSimulation of Thermomechanical Couplings of Viscoelastic Materials
Simulation of Thermomechanical Couplings of Viscoelastic Materials Frank Neff 1, Thomas Miquel 2, Michael Johlitz 1, Alexander Lion 1 1 Institute of Mechanics Faculty for Aerospace Engineering Universität
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 informationWhy does the Nazca plate slow down since the Neogene? Supplemental Information
GSA DATA REPOSITORY 2013004 Why does the Nazca plate slow down since the Neogene? Supplemental Information Javier Quinteros a,b,, Stephan V. Sobolev a,c a Deutsches GeoForschungsZentrum GFZ, Telegrafenberg,
More informationElements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004
Elements of Polymer Structure and Viscoelasticity David M. Parks Mechanics and Materials II 2.002 February 18, 2004 Outline Elements of polymer structure Linear vs. branched; Vinyl polymers and substitutions
More informationTheoretical Seismology. Astrophysics and Cosmology and Earth and Environmental Physics. Anelasticity. Fabio ROMANELLI
Theoretical Seismology Astrophysics and Cosmology and Earth and Environmental Physics Anelasticity Fabio ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it Intrinsic
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationSimulation of the Three-Dimensional Flow of Blood Using a Shear-Thinning Viscoelastic Fluid Model
Math. Model. Nat. Phenom. Vol. 6, No. 5, 2011, pp. 1-24 DOI: 1051/mmnp/20116501 Simulation of the Three-Dimensional Flow of Blood Using a Shear-Thinning Viscoelastic Fluid Model T. Bodnár 1, K.R. Rajagopal
More informationLinear viscoelastic behavior
Harvard-MIT Division of Health Sciences and Technology HST.523J: Cell-Matrix Mechanics Prof. Ioannis Yannas Linear viscoelastic behavior 1. The constitutive equation depends on load history. 2. Diagnostic
More informationCreep. Creep behavior of viscoelastic polymeric materials
B1 Version: 2.2_EN Date: 15. March 2018. BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS FACULTY OF MECHANICAL ENGINEERING DEPARTMENT OF POLYMER ENGINEERING Creep Creep behavior of viscoelastic polymeric
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 informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationRHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA
Communications to SIMAI Congress, ISSN 187-915, Vol. (7) DOI: 1.1685/CSC6157 RHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA A. CIANCIO, V. CIANCIO Department of Mathematics, University
More informationLarge time-step numerical modelling of the flow of Maxwell materials
Large time-step numerical modelling of the flow of Maxwell materials R.C. Bailey Geology and Physics Depts., University of Toronto, Toronto, Ontario M5S 1A7, Canada August 1, 2003 Maxwell viscoelastic
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 informationSplash singularity for a free-boundary incompressible viscoelastic fluid model
Splash singularity for a free-boundary incompressible viscoelastic fluid model Pierangelo Marcati (joint work with E.Di Iorio, S.Spirito at GSSI) Workshop 2016 Modeling Computation of Shocks Interfaces
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 informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationThe numerical method used for experiments is based on an explicit finite element
Bialas 1 Model Supplementary Data The numerical method used for experiments is based on an explicit finite element technique similar to the Fast Lagrangian Analysis of Continua (FLAC) method (Cundall,
More informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationUniversity Graz / Austria Institut für Chemie Volker Ribitsch
University Graz / Austria Institut für Chemie Volker Ribitsch 1 Rheology Oscillatory experiments Dynamic experiments Deformation of materials under non-steady conditions in the linear viscoelastic range
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
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 informationBeall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/ /g
GSA Data Repository 2018160 Beall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/10.1130/g39943.1. Numerical Modeling Methodology Stokes
More informationThe Hopf equation. The Hopf equation A toy model of fluid mechanics
The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationF7. Characteristic behavior of solids
F7. Characteristic behavior of solids F7a: Deformation and failure phenomena: Elasticity, inelasticity, creep, fatigue. à Choice of constitutive model: Issues to be considered è Relevance? Physical effect
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More informationEngineering Solid Mechanics
}} Engineering Solid Mechanics 1 (2013) 1-8 Contents lists available at GrowingScience Engineering Solid Mechanics homepage: www.growingscience.com/esm Impact damage simulation in elastic and viscoelastic
More informationAn introduction to implicit constitutive theory to describe the response of bodies
An introduction to implicit constitutive theory to describe the response of bodies Vít Průša prusv@karlin.mff.cuni.cz Mathematical Institute, Charles University in Prague 3 July 2012 Balance laws, Navier
More informationSEISMOLOGY. Master Degree Programme in Physics - UNITS Physics of the Earth and of the Environment ANELASTICITY FABIO ROMANELLI
SEISMOLOGY Master Degree Programme in Physics - UNITS Physics of the Earth and of the Environment ANELASTICITY FABIO ROMANELLI Department of Mathematics & Geosciences University of Trieste romanel@units.it
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 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 informationAnalysis and numerics of the mechanics of gels
Analysis and numerics of the mechanics of gels A DISSERTATION SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Brandon Michael Chabaud IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
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 informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More information