Goals Learn to simulate dynamic rupture model Explore the effects - - PowerPoint PPT Presentation

goals
SMART_READER_LITE
LIVE PREVIEW

Goals Learn to simulate dynamic rupture model Explore the effects - - PowerPoint PPT Presentation

Practice : dynamic rupture in 2D and 2.5D Huihui Weng Jean-Paul Ampuero Goals Learn to simulate dynamic rupture model Explore the effects of fault heterogeneities on dynamic ruptures Discussion: seismology problems 2D and 2.5D


slide-1
SLIDE 1

Practice:

dynamic rupture in 2D and 2.5D

Huihui Weng Jean-Paul Ampuero

slide-2
SLIDE 2

Goals

  • Learn to simulate dynamic rupture model
  • Explore the effects of fault heterogeneities
  • n dynamic ruptures
  • Discussion: seismology problems
slide-3
SLIDE 3

2D and 2.5D numerical codes

Download from GitHub

slide-4
SLIDE 4

2D and 2.5D numerical codes

Download from GitHub

slide-5
SLIDE 5

MATLAB code semlab

  • ---- run it directly

Fortran code sem2dpack

  • ---- installation is simple
  • ---- simulation is fast
slide-6
SLIDE 6

Model parameters

Frictional fault LY LX or LX/2 Discretized size S wave speed Density

Static friction coeffjcient Dynamic friction coeffjcient Normal stress Shear stress Slip-weakening distance

τ Slip τ0 τs τd Dc

Stations

slide-7
SLIDE 7

Seismogenic width W

x2 x1 x3

Slip profile

Seismogenic width W

Aseismic region Aseismic region Aseismic region Aseismic region

Sin(k3x3)

(3 equations) Reduce to 1 equation

Slip approximation

slide-8
SLIDE 8

Tutorial of sem2dpack Tutorial of semlab

slide-9
SLIDE 9

cd ${work_dir}/semlab-master/SEMLAB matlab &

  • pen this file from matlab:

sem2d_eq_dyn_tutorial.m

How to run semlab?

slide-10
SLIDE 10

Parameter setup in semlab

slide-11
SLIDE 11

Result presentation (figure1)

slide-12
SLIDE 12

Result presentation (figure2)

slide-13
SLIDE 13

Result presentation (figure3)

Tip Cohesive tail Stopping phase Stopping phase

slide-14
SLIDE 14

Result presentation (figure3)

Tip Cohesive tail Stopping phase Stopping phase

Slip Slip rate Stress Tail of cohesive zone Tip

Definition of tip: Slip rate > threshold Definition of cohesive tail: Slip = Dc

slide-15
SLIDE 15

Result presentation (figure4)

slide-16
SLIDE 16

STF and its spectrum (final)

slide-17
SLIDE 17

STF and spectra

slide-18
SLIDE 18

Add a heterogeneity

slide-19
SLIDE 19

STF comparison

slide-20
SLIDE 20

Models to test

Uniform fault Barrier

Model 1

Uniform fault Barrier

Model 2

abrupt stop Rough fault Barrier

Model 3

abrupt stop gradually stop

Ø Heterogeneity can be initial stress or Dc Ø Question: what are the differences of STF and spectrum? Ø Discussion: what are the differences between 2D and 2.5D models?

slide-21
SLIDE 21

Discussion

Uniform fault Barrier

Model 1

Uniform fault Barrier

Model 2

abrupt stop Rough fault Barrier

Model 3

abrupt stop gradually stop

Ø model 1 vs. model 2 Ø model 1 vs. model 3 Ø 2D vs. 2.5D

Apple to apple

slide-22
SLIDE 22

Discussion

Ø model 1 vs. model 2

slide-23
SLIDE 23

Discussion

Ø model 1 vs. model 3

slide-24
SLIDE 24

Discussion

Ø 2D vs. 2.5D

slide-25
SLIDE 25

Models to test

Ø Heterogeneity is initial stress Ø Question: do these heterogeneities increase high-frequency radiation? Ø Discussion: what parameters may control the high- frequency radiation?

Uniform fault

Model 4

Uniform fault

Model 5

Previous ruptured region Steps Initial stress Dc

slide-26
SLIDE 26

Discussion

Ø model 1 vs. model 4 Ø model 1 vs. model 5 Ø 2D vs. 2.5D

Apple to apple

Uniform fault

Model 4

Uniform fault

Model 5

Previous ruptured region Steps Initial stress Dc

slide-27
SLIDE 27

Discussion

Ø model 1 vs. model 4

slide-28
SLIDE 28

Discussion

Ø model 1 vs. model 5

slide-29
SLIDE 29

Tutorial of sem2dpack Tutorial of SEMLAB

slide-30
SLIDE 30

Installation of sem2dpack

slide-31
SLIDE 31

cd ${work_dir}/sem2dpack- 25D/EXAMPLES/2.5D_inplane

  • pen this file by vim or other method:

vi Par.inp

How to run sem2dpack?

slide-32
SLIDE 32

Parameter setup in sem2dpack

slide-33
SLIDE 33

Parameter setup in sem2dpack

slide-34
SLIDE 34

cd ${work_dir}/sem2dpack-25D/POST

  • pen the script by MATLAB:

create_movie.m run it by MATLAB

How to present results?

slide-35
SLIDE 35
slide-36
SLIDE 36

Discussion

ØFind seismic phases, such as P wave front, S wave, Rayleigh wave, etc. ØWhat parameters may control the formation of supershear rupture?

slide-37
SLIDE 37
slide-38
SLIDE 38
slide-39
SLIDE 39
slide-40
SLIDE 40
slide-41
SLIDE 41