Numerical modeling of dynamic rupture propagation
|
|
- Allyson Mills
- 5 years ago
- Views:
Transcription
1 IX. Slovenská geofyzikálna konferencia, jún 2011 Fakulta matematiky, fyziky a informatiky UK, Bratislava Numerical modeling of dynamic rupture propagation Martin Galis Peter Moczo Jozef Kristek Christian Pelties Comenius University Bratislava, Slovakia Slovak Academy of Sciences, Bratislava, Slovakia Ludwig-Maximiliams-Universität München, Germany
2 introduction FEM (Finite-Element Method) TSN (Traction at Split Nodes) approach and its implementation in FEM
3 introduction FEM (Finite-Element Method) TSN (Traction at Split Nodes) approach and its implementation in FEM ADER-DGM (Arbitrarily high-order DERivative Discontinuous Galerkin Method) implementation of the frictional boundary conditions in ADER-DGM
4 introduction FEM (Finite-Element Method) TSN (Traction at Split Nodes) approach and its implementation in FEM ADER-DGM (Arbitrarily high-order DERivative Discontinuous Galerkin Method) implementation of the frictional boundary conditions in ADER-DGM 3D numerical comparison FEM-ASA vs ADER-DGM effect of the numerical nucleation zone on the rupture propagation
5 FEM (Finite-Element Method)
6 displacement formulation of equation of motion - D-EqM ( ) ρu = σ, σ = λδ u, + μ( u, + u, ) i ij j ij ij k k i j j i
7 displacement formulation of equation of motion - D-EqM ( ) ρu = σ, σ = λδ u, + μ( u, + u, ) i ij j ij ij k k i j j i choose an element with n nodes e.g. tetrahedron with 4 nodes node
8 displacement formulation of equation of motion - D-EqM ( ) ρu = σ, σ = λδ u, + μ( u, + u, ) i ij j ij ij k k i j j i choose an element with n nodes e.g. tetrahedron with 4 nodes node choose shape functions and approximation to displacement in the element k k ( ) ( ) ( ) u x, t = s x U t ; k = 1,..., n i j j i k unique displacements U i at nodes continuity of displacement u i at a contact of elements
9 multiply D-EqM by the shape functions k ρus = σ, s ; k= 1,..., n i ij j k
10 multiply D-EqM by the shape functions k ρus = σ, s ; k= 1,..., n i ij j k integrate over an element Ω e i k ρusdv = Ω e σ, ij j k sdv
11 multiply D-EqM by the shape functions ρus = σ, s ; k= 1,..., n i k ij j k integrate over an element Ω e i k ρusdv = Ω e σ, ij j k sdv integrate the r.h.s. by parts e i k k k e ij j e Ω Ω Ω i ρusdv = σ s, dv + Ts ds
12 multiply D-EqM by the shape functions ρus = σ, s ; k= 1,..., n i k ij j k integrate over an element Ω e i k ρusdv = Ω e σ, ij j k sdv integrate the r.h.s. by parts e i k k k e ij j e Ω Ω Ω i ρusdv = σ s, dv + Ts ds e e e M u = r + surface-traction boundary term local mass matrix local vector of nodal displacements local vector of restoring forces
13 e e e M u = r + surface-traction boundary term assemble all elements covering volume Ω closed by surface Ω= ΩFS Ω Ω Dirichlet M u = r + 0 global mass matrix global vector of nodal displacements global vector of restoring forces
14 e e e M u = r + surface-traction boundary term Ω e k Ts ds i assemble all elements covering volume Ω closed by surface Ω= ΩFS Ω Ω Dirichlet M u = r + 0 global mass matrix global vector of nodal displacements global vector of restoring force the theoretical boundary term vanishes - at a contact of two elements due to traction continuity - at the free surface due to zero traction
15 e e e M u = r + surface-traction boundary term Ω e k Ts ds i assemble all elements covering volume Ω closed by surface Ω= ΩFS Ω Ω Dirichlet M u = r + 0 global mass matrix global vector of nodal displacements global vector of restoring force the theoretical boundary term vanishes - at a contact of two elements due to traction continuity - at the free surface due to zero traction in fact, however, the final discretization does not give - the traction continuity at a contact of two elements - zero traction at the free surface and thus the zero boundary term in the global equation they are just low-order approximated
16 TSN (Traction at Split Nodes) approach and its implementation in FEM
17 elastic halfspaces H - and H + do not interact H - H traction-free surface traction-free surface
18 elastic halfspaces H - and H + do not interact H - H traction-free surface traction-free surface halfspaces coupled by a constraint surface traction n c T { } { } a F () t AT c = + ( n,) t m - + fault + + a F () t AT c + = ( n,) t m
19 enforcement of the frictional boundary condition on the fault if T () t S() t ct sh then c ct T () t = T () t no slip
20 enforcement of the frictional boundary condition on the fault if T () t S() t ct sh then c ct T () t = T () t no slip if T () t > S() t ct sh then T c f () t T = () t sh sh c ct T () t = T () t n n + m + m ct f Dv ( t + dt 2) = A T T dt + mm
21 we apply our adaptive smoothing algorithm (ASA) to the trial traction before it is used for updating the slip rate in order to reduce spurious high-frequency oscillations of the slip rate
22 ADER-DGM (Arbitrarily high-order DERivative Discontinuous Galerkin Method)
23 velocity-stress formulation of equation of motion - VS-EqM σ λv, δ μ( v, v, ) = 0 ij k k ij i j j i ρv i σ, = 0 ij j
24 velocity-stress formulation of equation of motion - VS-EqM σ λv, δ μ( v, v, ) = 0 ρv ij k k ij i j j i i σ, = 0 ij j define a vector of unknown variables ( σxx, σyy, σzz, σxy, σyz, σzx, x, y, z) Q = v v v T
25 velocity-stress formulation of equation of motion - VS-EqM σ λv, δ μ( v, v, ) = 0 ρv ij k k ij i j j i i σ, = 0 ij j define a vector of unknown variables ( σxx, σyy, σzz, σxy, σyz, σzx, x, y, z) Q = v v v T VS-EqMin thematrixform Q + A Q, + B Q, + C Q, = 0 p pq q x pq q y pq q z A, B, C space-dependent matrices include material properties
26 velocity-stress formulation of equation of motion - VS-EqM σ λv, δ μ( v, v, ) = 0 ρv ij k k ij i j j i i σ, = 0 ij j define a vector of unknown variables ( σxx, σyy, σzz, σxy, σyz, σzx, x, y, z) Q = v v v T VS-EqMin thematrixform Q + A Q, + B Q, + C Q, = 0 p pq q x pq q y pq q z A, B, C space-dependent matrices include material properties tetrahedral element (e.g.) ( Q ) = Qˆ ( t) Φ ( x ) h p pk k j polynomial basis functions of an optional degree
27 multiply VS-EqM by a test function and integrate over an element volume Ω e ( ) Q Φ dv A Q, + B Q, + C Q, Φ dv = 0 p k e pq q x pq q y pq q z k Ω
28 multiply VS-EqM by a test function and integrate over an element volume Ω e ( ) Q Φ dv A Q, + B Q, + C Q, Φ dv = 0 p k e pq q x pq q y pq q z k Ω integrate the 2nd integral by parts Ω e (,,, ) Q Φ dv Φ A + Φ B + Φ C Q dv p k e k x pq k y pq k z pq q Ω + Φ F ds = Ω e k p 0 numerical flux introduced Q h because may be discontinuous at an element boundary
29 multiply VS-EqM by a test function and integrate over an element volume Ω e ( ) Q Φ dv A Q, + B Q, + C Q, Φ dv = 0 p k e pq q x pq q y pq q z k Ω integrate the 2nd integral by parts Ω e (,,, ) Q Φ dv Φ A + Φ B + Φ C Q dv p k e k x pq k y pq k z pq q Ω + Φ F ds = Ω e k p 0 numerical flux introduced because Q h may be discontinuous at an element boundary Riemann problem an evolution physically continuous problem with initial discontinuous approximation of unknowns across an interface to find a flux such that continuity of particle velocity and traction at an element boundary is assured
30 implementation of the frictional boundary conditions in ADER-DGM
31 imposing frictional boundary condition at S (the 2D example) frictional traction G = + imposed at S σ σ δσ xy xy xy
32 imposing frictional boundary condition at S (the 2D example) frictional traction G = + imposed at S σ σ δσ xy xy xy consequently G cs Q ( S ± ) = Q + 0, 0, δσ xy, 0, ± δσ xy μ and thus G cs G v y ± = vy ± ( σ xy σ xy) μ T
33 imposing frictional boundary condition at S (the 2D example) frictional traction G = + imposed at S σ σ δσ xy xy xy consequently G cs Q ( S ± ) = Q + 0, 0, δσ xy, 0, ± δσ xy μ and thus G cs G v y ± = vy ± ( σ xy σ xy) μ the slip rate is ( ) 2cs G y xy xy Δ v = σ σ μ T
34 imposing frictional boundary condition at S (the 2D example) frictional traction G = + imposed at S σ σ δσ xy xy xy consequently G cs Q ( S ± ) = Q + 0, 0, δσ xy, 0, ± δσ xy μ and thus G cs G v y ± = vy ± ( σ xy σ xy) μ the slip rate is ( ) 2cs G y xy xy Δ v = σ σ μ T all this means that the imposed frictional traction, G σ xy different from, instantly and locally generates an imposed slip rate parallel to the fault σ xy
35 the problem of the spurious high-frequency oscillations [m/s] [s] 7.5
36 the problem of the spurious high-frequency oscillations FEM + f ct mm Tsh = T Da A m + ( + m ) sh - collinearity n T + Dv = v v + fault ADER-DG σ G μ = σ + Δv 2c xy xy y s analogous to the analytical formulation in BIEM μ Δ v 2c s y - instantaneous (high-frequency) stress response to the slip rate + antiparallelism n T + Dv = v v + fault
37 3D numerical comparison FEM-ASA vs ADER-DGM planar fault receivers nucleation zone R1 R3 R2 7.5 km 6 km direction of initial shear traction friction law: linear slip weakening
38 3D numerical comparison FEM-ASA vs ADER-DGM planar fault receivers nucleation zone R1 R3 R2 7.5 km 6 km direction of initial shear traction friction law: linear slip weakening configuration sub-rayleigh supershear initial traction τ MPa 70.0 MPa static friction τ s MPa 73.5 MPa dynamic friction τ d MPa 63.0 MPa critical distance 0.4 m m strength par. S = ( τ τ ) ( τ τ ) s 0 0 d square nucleation zone, side length 3 km 3 km
39 3D numerical comparison FEM-ASA vs ADER-DGM order # of integ. points element size element type ADER-DG O m TET ADER-DG O m TET ADER-DG O m TET FEM-ASA O m HEX FEM-ASA O m HEX FEM-ASA O m HEX
40 sub-rayleigh rupture - slip-rate receiver R1 (antiplane) x-component receiver R2 (inplane) x-component receiver R3 x-component receiver R3 y-component
41 supershear rupture - slip-rate receiver R1 (antiplane) x-component receiver R2 (inplane) x-component receiver R3 x-component receiver R3 y-component
42 effect of the numerical nucleation zone on the rupture propagation simulations by FEM-ASA nucleation zones compared square with an abrupt change in traction for example SCEC ellipse with a smooth change in traction τ s + overshoot τ s τ 0 L
43 no rupture super-shear square nucleation zone sub-rayleigh 3.00 dimensionless size of nucleation zone L / L A strength parameter S = ( τ τ ) ( τ τ ) s 0 0 d
44 no rupture no rupture super-shear square nucleation zone super-shear elliptic nucleation zone sub-rayleigh sub-rayleigh 3.00 dimensionless size of nucleation zone L / L A strength parameter S = ( τ s τ0) ( τ0 τd )
45 no rupture no rupture super-shear square nucleation zone super-shear elliptic nucleation zone sub-rayleigh sub-rayleigh 3.00 dimensionless size of nucleation zone L / L A L L L D A / L / L A A L / U& R A strength parameter S = ( τ s τ0) ( τ0 τd )
46 thank you for your attention
Aspects of a Discontinuous Galerkin Approach for 3D Dynamic Rupture Modeling in the Case of a Complex Fault System
AGU Fall Meeting 2011 S54C-08 Aspects of a Discontinuous Galerkin Approach for 3D Dynamic Rupture Modeling in the Case of a Complex Fault System Christian Pelties Josep de la Puente (BSC), Jean Paul Ampuero
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2015) 200, 888 907 GJI Seismology doi: 10.1093/gji/ggu436 On the initiation of sustained slip-weakening ruptures by localized stresses M. Galis, 1 C.
More informationLokálne chyby numerických schém
Slovenská geofyzikálna konferencia 211 Lokálne chyby numerických schém Peter Moczo Jozef Kristek Martin Galis Emmanuel Chaljub, Vincent Etienne FMFI UK Bratislava GFÚ SAV Bratislava ISTerre Université
More informationA new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities
A new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities Setare Hajarolasvadi Ahmed E. Elbanna https://www.youtube.com/watch?v=bltx92tuwha MOTIVATION
More informationFault Representation Methods for Spontaneous Dynamic Rupture Simulation. Luis A. Dalguer
Fault Representation Methods for Spontaneous Dynamic Rupture Simulation Luis A. Dalguer Computational Seismology Group Swiss Seismological Service (SED) ETH-Zurich July 12-18, 2011 2 st QUEST Workshop,
More informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationDefmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More information4. Numerical algorithms for earthquake rupture dynamic modeling
Research Signpost 37/661 (2), Fort P.O. Trivandrum-695 023 Kerala, India The Mechanics of Faulting: From Laboratory to Real Earthquakes, 2012: 93-124 ISBN: 978-81-308-0502-3 Editors: Andrea Bizzarri and
More informationThree-Dimensional Dynamic Rupture Simulation with a High-order Discontinuous Galerkin Method on Unstructured Tetrahedral Meshes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, Three-Dimensional Dynamic Rupture Simulation with a High-order Discontinuous Galerkin Method on Unstructured Tetrahedral Meshes Christian Pelties
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
More informationComparison of finite difference and boundary integral solutions to three-dimensional spontaneous rupture
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2005jb003813, 2005 Comparison of finite difference and boundary integral solutions to three-dimensional spontaneous rupture Steven M. Day and Luis
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 informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 9, 21 Sep. 2017 www.geosc.psu.edu/courses/geosc508 Rate and State Friction Velocity step test to measure RSF parameters SHS test to measure RSF parameters
More informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
More informationEML4507 Finite Element Analysis and Design EXAM 1
2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationRupture dynamics with energy loss outside the slip zone
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003191, 2005 Rupture dynamics with energy loss outside the slip zone D. J. Andrews U.S. Geological Survey, Menlo Park, California, USA Received
More informationHow to promote earthquake ruptures: Different nucleation strategies. in a dynamic model with slip weakening friction
How to promote earthquake ruptures: Different nucleation strategies in a dynamic model with slip weakening friction Andrea Bizzarri Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna, Italy
More informationBake, shake or break - and other applications for the FEM. 5: Do real-life experimentation using your FEM code
Bake, shake or break - and other applications for the FEM Programming project in TMA4220 - part 2 by Kjetil André Johannessen TMA4220 - Numerical solution of partial differential equations using the finite
More informationDeformable Materials 2 Adrien Treuille. source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes.
Deformable Materials 2 Adrien Treuille source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes. Goal Overview Strain (Recap) Stress From Strain to Stress Discretization Simulation Overview Strain
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationSurface Waves and Free Oscillations. Surface Waves and Free Oscillations
Surface waves in in an an elastic half spaces: Rayleigh waves -Potentials - Free surface boundary conditions - Solutions propagating along the surface, decaying with depth - Lamb s problem Surface waves
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationDynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement
Emily Morton Geop 523: Theoretical Seismology 29 April 2011 Introduction Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Earthquakes can be triggered by other earthquakes.
More informationPrediction of geometric dimensions for cold forgings using the finite element method
Journal of Materials Processing Technology 189 (2007) 459 465 Prediction of geometric dimensions for cold forgings using the finite element method B.Y. Jun a, S.M. Kang b, M.C. Lee c, R.H. Park b, M.S.
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 informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationSource parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop
Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for
More informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationDEFORMATION PATTERN IN ELASTIC CRUST
DEFORMATION PATTERN IN ELASTIC CRUST Stress and force in 2D Strain : normal and shear Elastic medium equations Vertical fault in elastic medium => arctangent General elastic dislocation (Okada s formulas)
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
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 informationHeterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21306, doi:10.1029/2008gl035614, 2008 Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationFault-Zone Properties arid Earthquake Rupture Dynamics
Fault-Zone Properties arid Earthquake Rupture Dynamics Eiichi Fukuyama National Research Institute for Earth Science,and Disaster Prevention Tsukuba, Japan i ':-> i' ' -':'.." \,' " '' L VH ELSEVIER AMSTERDAM
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 information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationSeismic and aseismic processes in elastodynamic simulations of spontaneous fault slip
Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Most earthquake simulations study either one large seismic event with full inertial effects or long-term slip history
More informationfriction friction a-b slow fast increases during sliding
µ increases during sliding faster sliding --> stronger fault --> slows sliding leads to stable slip: no earthquakes can start velocity-strengthening friction slow fast µ velocity-strengthening friction
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationEarthquake modelling with complex geometries by curved grid finite-difference method
FSEF2017, Shenzhen, China Nov 27 Dec. 1, 2017 Earthquake modelling with complex geometries by curved grid finite-difference method Zhenguo Zhang Department of Earth and Space Sciences, Southern University
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 informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 16
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 16 In the last lectures, we have seen one-dimensional boundary value
More informationShear stresses around circular cylindrical openings
Shear stresses around circular cylindrical openings P.C.J. Hoogenboom 1, C. van Weelden 1, C.B.M. Blom 1, 1 Delft University of Technology, the Netherlands Gemeentewerken Rotterdam, the Netherlands In
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationFREE VIBRATION OF A THERMO-PIEZOELECTRIC PLATE
Inter national Journal of Pure and Applied Mathematics Volume 113 No. 11 217, 217 225 ISSN: 1311-88 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu FREE VIBRATION
More informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationSpectral element modeling of spontaneous earthquake rupture on rate and state faults: Effect of velocity-strengthening friction at shallow depths
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005553, 2008 Spectral element modeling of spontaneous earthquake rupture on rate and state faults: Effect of velocity-strengthening friction
More informationECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017
ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM Prof. Peter Bermel February 6, 2017 Outline Static Equilibrium Dynamic Equilibrium Thermal transport mechanisms Thermal transport modeling
More informationNumerical study on multi-scaling earthquake rupture
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L266, doi:1.129/23gl1878, 24 Numerical study on multi-scaling earthquake rupture Hideo Aochi Institut de Radioprotection et de Sûreté Nucléaire, France Satoshi Ide
More informationCalculation of the local rupture speed of dynamically propagating earthquakes
ANNALS OF GEOPHYSICS,,,, S; doi:./ag-79 Calculation of the local rupture speed of dynamically propagating earthquakes Andrea Bizzarri Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna,
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationSPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES
2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationTwo ways to think about the dynamics of earthquake ruptures
Two ways to think about the dynamics of earthquake ruptures (1) In terms of friction (2) In terms of fracture mechanics Scholz describes conditions for rupture propagation (i.e. instability) via energy
More informationDEFORMATION PATTERN IN ELASTIC CRUST
DEFORMATION PATTERN IN ELASTIC CRUST Stress and force in 2D Strain : normal and shear Elastic medium equations Vertical fault in elastic medium => arctangent General elastic dislocation (Okada s formulas)
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationJournalofGeophysicalResearch: SolidEarth
JournalofGeophysicalResearch: SolidEarth RESEARCH ARTICLE Key Points: Relations between source dimensions and corner frequencies are established Stress drop estimates are affected by differences in source
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 57 AA242B: MECHANICAL VIBRATIONS Dynamics of Continuous Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations: Theory
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
More informationFinite-amplitude shear wave in pre-stressed thin elastomers
Wave Motion 43 2005) 20 28 Finite-amplitude shear wave in pre-stressed thin elastomers Sia Nemat-Nasser, Alireza V. Amirkhizi Center of Excellence for Advanced Materials, Mechanical and Aerospace Engineering,
More informationOn the diminishing of spurious oscillations in explicit finite element analysis of linear and non-linear wave propagation and contact problems
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16315 On the diminishing of spurious oscillations
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More informationEarthquake stress drop estimates: What are they telling us?
Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationStress Integration for the Drucker-Prager Material Model Without Hardening Using the Incremental Plasticity Theory
Journal of the Serbian Society for Computational Mechanics / Vol. / No., 008 / pp. 80-89 UDC: 59.74:004.0 Stress Integration for the Drucker-rager Material Model Without Hardening Using the Incremental
More informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More information16.20 Techniques of Structural Analysis and Design Spring Instructor: Raúl Radovitzky Aeronautics & Astronautics M.I.T
16.20 Techniques of Structural Analysis and Design Spring 2013 Instructor: Raúl Radovitzky Aeronautics & Astronautics M.I.T February 15, 2013 2 Contents 1 Stress and equilibrium 5 1.1 Internal forces and
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
More informationLarge-scale multi-physics earthquake scenarios with the ADER-DG method on modern supercomputers
Largescale multiphysics earthquake scenarios with the ADERDG method on modern supercomputers Stephanie Wollherr, Dr. AliceAgnes Gabriel, Dr. Betsy Madden, Thomas Ulrich LudwigsMaximilians Universität (LMU)
More informationDynamic stress field of a kinematic earthquake source model with k-squared slip distribution
Geophys. J. Int. (2007) doi: 10.1111/j.1365-246X.2007.03548.x Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution Jan Burjánek and Jiří Zahradník Charles University
More informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationScale Dependence in the Dynamics of Earthquake Rupture Propagation: Evidence from Geological and Seismological Observations
Euroconference of Rock Physics and Geomechanics: Natural hazards: thermo-hydro-mechanical processes in rocks Erice, Sicily, 25-30 September, 2007 Scale Dependence in the Dynamics of Earthquake Rupture
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationSustained Petascale Performance of Seismic Simulations with SeisSol
SIAM EX Workshop on Exascale Applied Mathematics Challenges and Opportunities Sustained Petascale Performance of Seismic Simulations with SeisSol M. Bader, A. Breuer, A. Heinecke, S. Rettenberger C. Pelties,
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationA FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA
A FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA Eren Ucar a, Eirik Keilegavlen a, Inga Berre a,b, Jan Martin Nordbotten a,c a Department of Mathematics, University of Bergen, Bergen,
More informationIntroduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0
Introduction to Finite Element Analysis Using Pro/MECHANICA Wildfire 5.0 Randy H. Shih Oregon Institute of Technology SDC PUBLICATIONS Schroff Development Corporation www.schroff.com Better Textbooks.
More information