Overview of the PyLith Finite Element Code
|
|
|
- Alexander Banks
- 5 years ago
- Views:
Transcription
1 Overview of the PyLith Finite Element Code Charles Williams (GNS Science) Brad Aagaard (USGS) Matt Knepley (University of Chicago)
2 CIG Software
3 PyLith is a Tool for Crustal Deformation Modeling Typical CDM workflow
4 PyLith Problem Types Quasi-static modeling associated with earthquakes: Strain accumulation associated with interseismic deformation: What is the stressing rate on faults X and Y? Where is strain accumulating in the crust? Coseismic stress changes and fault slip: What was the slip distribution in earthquake A? How did earthquake A change the stresses on faults X and Y? Postseismic relaxation of the crust: What rheology is consistent with observed postseismic deformation? Can aseismic creep or afterslip explain the deformation?
5 PyLith Problem Types (cont.) Dynamic modeling associated with earthquakes: Modeling of strong ground motions: Forecasting the amplitude and spatial variation in ground motion for scenario earthquakes. Coseismic stress changes and fault slip: How did earthquake A change the stresses on faults X and Y? Earthquake rupture behavior: What fault constitutive models/parameters are consistent with the observed rupture propagation in earthquake A?
6 PyLith Problem Types (cont.) Volcanic deformation associated with magma chambers and/or dikes/sills: Magma chamber inflation/deflation: What is the geometry and pressure associated with the magma chamber? How close to failure is the material surrounding the chamber? Dike/sill intrusions: What is the orientation and amount of opening associated with the intrusion?
7 Usage Examples Frictional fault behavior near a subducting seamount SUTRA fluid flow model PyLith frictional model
8 Usage Examples (cont.) Modeling of Slow Slip Events
9 Usage Examples (cont.) Generation of Green s Functions for Darfield Earthquake Elastic properties from NZ-wide velocity model Surface response due to unit Slip on Greendale fault
10 Usage Examples (cont.) Dynamic rupture simulation of Darfield earthquake
11 Usage Examples (cont.) Stress changes due to dike intrusion at Mt. Asama volcano Maximum principal stress magnitude Maximum principal stress orientation
12 Ingredients for Running PyLith Finite-element mesh Mesh exported from LaGriT Mesh exported from CUBIT Mesh constructed by hand (PyLith mesh ASCII format) Simulation parameters Command-line Parameters in.cfg file Spatial databases for physical properties, boundary conditions, and rupture parameters SCEC CVM-H, USGS Bay Area Velocity Model Simple ASCII files
13 Examples Spatial Databases Uniform value for Dirichlet (0D) Piecewise linear traction variation for Neumann BC (1D) SCEC CVM-H seismic velocity model (3D) Generally independent of problem discretization Available spatial databases UniformDB Optimized for uniform value SimpleDB Simple ASCII files (0D, 1D, 2D, or 3D) SCECCVMH SCEC CVM-H seismic velocity model v5.3 SimpleGridDB Gridded data (1D, 2D, or 3D) ZeroDispDB Special case of UniformDB
14 PyLith: Focus on Geodynamics Leverage packages developed by computational scientists
15 Fault Implementation
16 Implementing Fault Slip with Advantages Lagrange Multipliers Fault implementation is local to cohesive cell Solution includes tractions generating slip (Lagrange multipliers) Retains block structure of matrix, including symmetry Offsets in mesh mimic slip on natural faults Disadvantages Cohesive cells require adjusting topology of finiteelement mesh Scalable preconditioner/solver is more complex
17 PyLith 1.9 Features Time integration schemes and formulations Implicit for quasistatic problems (neglect inertial terms) Infinitesimal strains Small strains (~30-40%) Explicit for dynamic problems Infinitesimal strains Small strains (~30-40%) Numerical damping via viscosity
18 PyLith 1.9 Features (cont.) Bulk constitutive models Elastic model (1D, 2D, and 3D) Linear Maxwell viscoelastic models (2D and 3D) Generalized Maxwell viscoelastic models (2D and 3D) Power-law viscoelastic model (2D and 3D) Drucker-Prager elastoplastic model (2-D and 3-D)
19 PyLith 1.9 Features (cont.) Boundary and interface conditions Time-dependent Dirichlet BC Time-dependent Neumann (traction) BC Absorbing BC Kinematic (prescribed slip) faults w/multiple ruptures Dynamic (friction) fault interfaces Time-dependent point forces Gravitational body forces
20 PyLith 1.9 Features (cont.) Fault constitutive models Kinematic (specified slip) Static friction Linear slip-weakening friction Linear time-weakening friction Dieterich-Ruina rate and state friction w/aging law
21 PyLith 1.9 Features (cont.) Automatic and user-controlled time stepping Ability to specify initial stress/strain state Importing meshes LaGriT: GMV/Pset CUBIT: Exodus II ASCII: PyLith mesh ASCII format (for toy problems)
22 PyLith 1.9 Features (cont.) Output: VTK and HDF5 files (parallel I/O for HDF5) Solution over volume Solution over surface boundary Solution interpolated to user-specified points State variables (e.g., stress and strain) for each material Fault information (e.g., slip and tractions)
23 PyLith 1.9 Features (cont.) Automatic conversion of units for all parameters Ability to use geographic projections (using Proj.4) Automatic generation of static Green s functions Parallel uniform global refinement PETSc linear and nonlinear solvers Custom preconditioner with algebraic multigrid solver
24 Short-term priorities PyLith Development Under-the-hood improvements New finite-element data structures [done] Support higher order basis functions [in progress] Provides much higher resolution for a given mesh Prepare for multi-physics [done] Multi-cycle earthquake modeling Resolve interseismic, coseismic, and postseismic deformation Elastic/viscoelastic/plastic rheologies Coseismic slip, afterslip, and creep Long-term priorities Multiphysics: Elasticity + Fluid flow + Heat flow Scaling to 1000 processors
25 PyLith Development: Planned v2.0 (Summer 2013) Releases New finite-element data structures Support for higher order basis functions v2.1 (Spring 2014) Coupling of quasi-static and dynamic simulations Moment tensor point sources v2.2 (Fall 2014) Support for incompressible elasticity Heat and fluid flow coupled to elastic deformation v2.x Support for finite-element integrations on GPUs
26 Where to Get More Info Mini-workshop this afternoon Simple example including topography, a spherical magma chamber, and a dike Instructors Charles Williams Adrian Shelley CIG website Binaries available for Mac OS X, Linux, and Windows Source code + installer available Extensive manual + example problems Subscribe to cig-short mailing list Low/moderate traffic list to report problems, get help, etc.
Earthquake and Volcano Deformation
Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c
Introduction: Advancing Simulations of Sequences of Earthquakes and Aseismic Slip (SEAS)
Introduction: Advancing Simulations of Sequences of Earthquakes and Aseismic Slip (SEAS) Brittany Erickson (Portland State University) Junle Jiang (University of California, San Diego) SCEC DR-SEAS Workshop,
Seismic 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
Surface changes caused by erosion and sedimentation were treated by solving: (2)
GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)
CAEFEM v9.5 Information
CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
Surface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes
Surface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes Brad Aagaard May 21, 2009 Surface Rupture from Kinematic Rupture Models Objective Holistic approach to calculating surface
Defmod, 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
Data Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
A 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
Can geodetic strain rate be useful in seismic hazard studies?
Can geodetic strain rate be useful in seismic hazard studies? F. Riguzzi 1, R. Devoti 1, G. Pietrantonio 1, M. Crespi 2, C. Doglioni 2, A.R. Pisani 1 Istituto Nazionale di Geofisica e Vulcanologia 2 Università
Fault 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,
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological
DETAILS 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
GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise
GPS Strain & Earthquakes Unit 5: 2014 South Napa earthquake GPS strain analysis student exercise Strain Analysis Introduction Name: The earthquake cycle can be viewed as a process of slow strain accumulation
Modeling Approaches That Reproduce a Range of Fault Slip Behaviors: What We Have and What We Need Nadia Lapusta. California Institute of Technology
Modeling Approaches That Reproduce a Range of Fault Slip Behaviors: What We Have and What We Need Nadia Lapusta California Institute of Technology Modeling Approaches That Reproduce a Range of Fault Slip
Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction
Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose
USER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
The Mechanics of Earthquakes and Faulting
The Mechanics of Earthquakes and Faulting Christopher H. Scholz Lamont-Doherty Geological Observatory and Department of Earth and Environmental Sciences, Columbia University 2nd edition CAMBRIDGE UNIVERSITY
Rheology 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:
Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation
Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences
Development of a Predictive Simulation System for Crustal Activities in and around Japan - II
Development of a Predictive Simulation System for Crustal Activities in and around Japan - II Project Representative Mitsuhiro Matsu'ura Graduate School of Science, The University of Tokyo Authors Mitsuhiro
USER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
Variability of earthquake nucleation in continuum models of rate-and-state faults and implications for aftershock rates
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005154, 2008 Variability of earthquake nucleation in continuum models of rate-and-state faults and implications
Multiphysics and multiscale earthqauke modeling
Multiphysics and multiscale earthqauke modeling C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 29, 2017 Multiphysics Poro(visco-)elasticity, (quasi-)static and dynamic processes. ground overburden reservoir
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University
Slow Slip and Tremor Along San Andreas fault system
Slow Slip and Tremor Along San Andreas fault system Slow slip in Upper Crust Aseismic creep, afterslip and spontaneous slow slip events on some faults in upper 15 km Mostly aseismic but accompanied by
Elizabeth H. Hearn modified from W. Behr
Reconciling postseismic and interseismic surface deformation around strike-slip faults: Earthquake-cycle models with finite ruptures and viscous shear zones Elizabeth H. Hearn hearn.liz@gmail.com modified
Introduction to Displacement Modeling
Introduction to Displacement Modeling Introduction Deformation on the Earth surface informs us about processes and material properties below surface Observation tools: GPS (static-dynamic) InSAR (static)
Toward a SCEC Community Rheology Model: TAG Kickoff and Workshop SCEC Workshop Proposal Final Report
Toward a SCEC Community Rheology Model: TAG Kickoff and Workshop SCEC Workshop Proposal 17206 Final Report A one-day Community Rheology Model workshop was held at the Palm Springs Hilton on the Saturday
Two 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
Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L19305, doi:10.1029/2008gl035270, 2008 Predicted reversal and recovery of surface creep on the Hayward fault following the 1906 San Francisco earthquake D. A. Schmidt
Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller
Lecture 20: Slow Slip Events and Stress Transfer GEOS 655 Tectonic Geodesy Jeff Freymueller Slow Slip Events From Kristine Larson What is a Slow Slip Event? Slip on a fault, like in an earthquake, BUT
Fault-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
USER S MANUAL. 1D Seismic Site Response Analysis Example. University of California: San Diego.
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 2, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
Geophysical Journal International
Geophysical Journal International Geophys. J. Int. (2011) 186, 1389 1403 doi: 10.1111/j.1365-246X.2011.05117.x Shallow slip deficit due to large strike-slip earthquakes in dynamic rupture simulations with
Constitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
Friction can increase with hold time. This happens through growth and increasing shear strength of contacts ( asperities ).
Friction can increase with hold time. This happens through growth and increasing shear strength of contacts ( asperities ). If sliding speeds up, the average lifespan of asperities decreases This means
SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation L.L.A. Vermeersen (1), H. Schotman (1), M.-W. Jansen (1), R. Riva (1) and R. Sabadini (2) (1) DEOS,
Deformation cycles of great subduction earthquakes in a viscoelastic Earth
Deformation cycles of great subduction earthquakes in a viscoelastic Earth Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada School of Earth and Ocean Science, University of Victoria????
2D Embankment and Slope Analysis (Numerical)
2D Embankment and Slope Analysis (Numerical) Page 1 2D Embankment and Slope Analysis (Numerical) Sunday, August 14, 2011 Reading Assignment Lecture Notes Other Materials FLAC Manual 1. 2. Homework Assignment
The 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,
Spectral 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
Fundamentals 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
March 2017 SCEC Rupture Dynamics Code Validation Workshop
Presentation for March 1, 2017 SCEC Boardroom Los Angeles, CA March 2017 SCEC Rupture Dynamics Code Validation Workshop Ruth A. Harris (U.S. Geological Survey) Ralph J. Archuleta (UC Santa Barbara) INTRODUCTION
Verification of the asperity model using seismogenic fault materials Abstract
Verification of the asperity model using seismogenic fault materials Takehiro Hirose*, Wataru Tanikawa and Weiren Lin Kochi Institute for Core Sample Research/JAMSTEC, JAPAN * Corresponding author: hiroset@jamstec.go.jp
Chapter 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
Generic Earthquake Simulator
Generic Earthquake Simulator by Terry E. Tullis, Keith Richards-Dinger, Michael Barall, James H. Dieterich, Edward H. Field, Eric Heien, Louise H. Kellogg, Fred Pollitz, John Rundle, Michael Sachs, Donald
Christine M. Puskas. Objective: A geophysical position where a background in geophysical modeling, data processing, and GPS studies can be applied.
Christine M. Puskas 777 E South Temple 9C Salt Lake City, UT 84102 www.uusatrg.utah.edu/people/christine (801) 581-7856 Office (801) 521-4107 Home c.puskas@utah.edu Education Ph.D. in Geophysics, University
Earthquake and Volcano Clustering at Mono Basin (California)
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Earthquake and Volcano Clustering at Mono Basin (California) D. La Marra *,1, A. Manconi 2,3 and M. Battaglia 1 1 Dept of Earth Sciences,
Earthquake nucleation. Pablo Ampuero Caltech Seismolab
Earthquake nucleation Pablo Ampuero Caltech Seismolab How do earthquakes start? Do small and large earthquakes start differently? Predictive value of earthquake onset and foreshock sequences? Seismological
Large-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)
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
Analytic and Numeric Tests of Fourier Deformation Model (Copyright 2003, Bridget R. Smith and David T. Sandwell)
Analytic and Numeric Tests of Fourier Deformation Model (Copyright 2003, Bridget R. Smith and David T. Sandwell) Although the solutions of our Fourier deformation model have been checked using computer
Secondary Project Proposal
Secondary Project Proposal Post-seismic deformation of Chi-chi earthquake Yunyue (Elita) Li 11:, Wednesday, June 2, 21 Li 2 Secondary project proposal Personal prospective MOTIVATION My interests for earthquake
ETH Swiss Federal Institute of Technology Zürich
Swiss Federal Institute of Technology Zürich Earthquake Statistics using ZMAP Recent Results Danijel Schorlemmer, Stefan Wiemer Zürich, Swiss Seismological Service, Switzerland Contributions by: Matt Gerstenberger
An Introduction to Virtual Quake
An Introduction to Virtual Quake Kasey Schultz and John Wilson University of California, Davis Department of Physics - Rundle Group CIG Webinar December 3, 2015 Outline for Today! What is Virtual Quake
GEOL 10: Environmental Geology Mid Term II Study Guide
(1) 10 pts. Use the seismographs below, along with the travel time plot for California, and the location map for seismographs in this region to do the following: a. Determine the origin time for the earthquake
Afterslip, slow earthquakes and aftershocks: Modeling using the rate & state friction law
Afterslip, slow earthquakes and aftershocks: Modeling using the rate & state friction law Agnès Helmstetter (LGIT Grenoble) and Bruce Shaw (LDE0 Columbia Univ) Days after Nias earthquake Cumulative number
2D Liquefaction Analysis for Bridge Abutment
D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering
Megathrust earthquakes: How large? How destructive? How often? Jean-Philippe Avouac California Institute of Technology
Megathrust earthquakes: How large? How destructive? How often? Jean-Philippe Avouac California Institute of Technology World seismicity (data source: USGS) and velocities relative to ITRF1997 (Sella et
DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
Annual Report for Research Work in the fiscal year 2005
JST Basic Research Programs C R E S T (Core Research for Evolutional Science and Technology) Annual Report for Research Work in the fiscal year 2005 Research Area : High Performance Computing for Multi-scale
When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress?
When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress? How could you exert a non-zero shear stress on the ground? Hydrostatic Pressure (fluids)
LETTER Earth Planets Space, 65, , 2013
LETTER Earth Planets Space, 65, 1041 1046, 2013 Numerical simulation of crustal deformation using a three-dimensional viscoelastic crustal structure model for the Japanese islands under east-west compression
seismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model
seismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model properties inversion of seismic observables for 3D Earth structure
ENGN 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
Scenario Earthquake and Tsunami Simulations for a Pacific Rim GNSS Tsunami Early Warning System
Scenario Earthquake and Tsunami Simulations for a Pacific Rim GNSS Tsunami Early Warning System Kasey Schultz February 2, 2015 Abstract The first operational prototype for a Pacific Rim Tsunami Early Warning
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
SCEC INCITE OLCF PROJECT PROGRESS REPORT,
SCEC INCITE OLCF PROJECT PROGRESS REPORT, 2015 Title: PI and Co-PI(s): Thomas H. Jordan, Jacobo Bielak, Kim Olsen, Yifeng Cui, Po Chen, Ricardo Taborda, Philip Maechling Applying Institution/Organization:
Spectral 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,
MHA042 - 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
Variations in Tremor Activity and Implications for Lower Crustal Deformation Along the Central San Andreas Fault
Variations in Tremor Activity and Implications for Lower Crustal Deformation Along the Central San Andreas Fault David R. Shelly USGS, Menlo Park Shelly and Hardebeck, GRL, 2010 Collaborators: Jeanne Hardebeck
Rate and State-Dependent Friction in Earthquake Simulation
Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity
Fonctionnalités de Code_Aster version 12
Fonctionnalités de Code_Aster version 12 Outline Version 12 in a nutshell Fracture mechanics Non-linear constitutive laws Linear and non-linear dynamics Numerical methods Architecture, ergonomics, performances
A 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,
MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
Differential 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
Qualitative modeling of earthquakes and aseismic slip in the Tohoku-Oki area. Nadia Lapusta, Caltech Hiroyuki Noda, JAMSTEC
Qualitative modeling of earthquakes and aseismic slip in the Tohoku-Oki area Nadia Lapusta, Caltech Hiroyuki Noda, JAMSTEC Constitutive law on the fault: Rate-and-state friction at low slip rates + Potential
Elasto-plastic FEM models explain the emplacement of shallow magma intrusions in volcanic complexes
Elasto-plastic FEM models explain the emplacement of shallow magma intrusions in volcanic complexes Andrea Bistacchi University of Milano Bicocca - Structural Geology Lab Excerpt from the Proceedings of
1 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
The 3 rd SCEC CSM workshop
The 3 rd SCEC CSM workshop Welcome on behalf of the organizers Jeanne Hardebeck Brad Aagaard David Sandwell Bruce Shaw John Shaw Thorsten Becker Thanks for playing! SCEC Community Stress Model (CSM) Community
Jihoon 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,
Numerical Methods in Geophysics. Introduction
: Why numerical methods? simple geometries analytical solutions complex geometries numerical solutions Applications in geophysics seismology geodynamics electromagnetism... in all domains History of computers
Numerical approximation of rate-and-state friction problems
Preprint No. 1 Numerical approximation of rate-and-state friction problems Elias Pipping Ralf Kornhuber Matthias Rosenau Onno Oncken April, 2015 Project B01 This research has been funded by Deutsche Forschungsgemeinschaft
Dynamic 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.
Chapter 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
Syllabus and Course Description Geophysical Geodesy Fall 2013 GPH 411/611
Syllabus and Course Description Geophysical Geodesy Fall 2013 GPH 411/611 Course Location: LME 422 Course Time: Tuesday & Thursday 12:00-1:15 Units: 3 Instructor Name: Bill Hammond Office Location: SEM
Multibody System Dynamics: MBDyn Hydraulics Modeling
Multibody System Dynamics: MBDyn Hydraulics Modeling Pierangelo Masarati Politecnico di Milano Dipartimento di Scienze e Tecnologie Aerospaziali Outline 2 Introduction Modeling
Source 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
ALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis
ALASKA ENERGY AUTHORITY AEA11-022 ENGINEERING FEASIBILITY REPORT Appendix B8 Finite Element Analysis Susitna-Watana Hydroelectric Project Alaska Energy Authority FERC Project No. 14241 December 2014 Seismic
High-Resolution Finite Volume Methods and Adaptive Mesh Refinement
High-Resolution Finite Volume Methods and Adaptive Mesh Refinement Randall J. LeVeque Department of Applied Mathematics University of Washington CLAWPACK and TsunamiClaw Software http://www.amath.washington.edu/~claw
Buckling, Postbuckling, and Collapse Analysis with Abaqus. Abaqus 2017
Buckling, Postbuckling, and Collapse Analysis with Abaqus Abaqus 2017 About this Course Course objectives Upon completion of this course you will be able to: Perform linear eigenvalue buckling analysis
Dynamics Manual. Version 7
Dynamics Manual Version 7 DYNAMICS MANUAL TABLE OF CONTENTS 1 Introduction...1-1 1.1 About this manual...1-1 2 Tutorial...2-1 2.1 Dynamic analysis of a generator on an elastic foundation...2-1 2.1.1 Input...2-1
Physics and Chemistry of the Earth and Terrestrial Planets
MIT OpenCourseWare http://ocw.mit.edu 12.002 Physics and Chemistry of the Earth and Terrestrial Planets Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much
Mechanics 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
Frictional rheologies have a wide range of applications in engineering
A liquid-crystal model for friction C. H. A. Cheng, L. H. Kellogg, S. Shkoller, and D. L. Turcotte Departments of Mathematics and Geology, University of California, Davis, CA 95616 ; Contributed by D.
The Non-volcanic tremor observation in Northern Cascadia. Hsieh Hsin Sung 3/22
The Non-volcanic tremor observation in Northern Cascadia Hsieh Hsin Sung 3/22 Reference Kao, H., S. J. Shan, H. Dragert, and G. Rogers (2009), Northern Cascadia episodic tremor and slip: A decade of observations