Skip to content

Duo/tpv6 7 #48

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 5 commits into
base: master
Choose a base branch
from
Open

Duo/tpv6 7 #48

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
179f0cd
add coarse mesh for training
Oct 15, 2020
f756342
fix tpv12
Sep 20, 2021
93763b0
Merge branch 'master' of github.com:SeisSol/Examples
Jul 1, 2022
b6ed390
fix fault mu_s setup
Jul 1, 2022
7458dad
I change the parameters which seems compare with scec better
Jul 1, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
19 changes: 17 additions & 2 deletions tpv6_7/fault.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,9 +22,24 @@
return 81600000;
}
return 70000000;
[mu_s, mu_d, d_c, cohesion]: !ConstantMap
[mu_s]: !FunctionMap
map:
x: return abs(x);
z: return abs(z);
components:
- !AxisAlignedCuboidalDomainFilter
limits:
x: [0, 15000]
z: [0, 15000]
components: !ConstantMap
map:
mu_s: 0.677
- !ConstantMap
map:
mu_s: 1000.0
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

should be 10000.

[mu_d, d_c, cohesion]: !ConstantMap
map:
mu_s: 0.677
mu_d: 0.525
d_c: 0.4
cohesion: 0.0

10 changes: 4 additions & 6 deletions tpv6_7/parameters.par
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -17,8 +17,8 @@ FL = 6 ! Friction law
!0: none, 16:LSW, 103: RS with strong velocity weakening
!yaml file defining spatial dependance of fault properties
ModelFileName = 'fault.yaml'
v_star = 1 ! Prakash-Clifton regularization parameter
L = 0.01
v_star = 0.1 ! Prakash-Clifton regularization parameter
L = 0.5
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These values are not prescribed by the Benchmark, right?
Why did you change them?



!reference vector for defining strike and dip direction
Expand All @@ -29,7 +29,7 @@ refPointMethod = 1

RF_output_on = 1 ! Rupture front ascii output
OutputPointType = 5 ! Type (0: no output, 3: ascii file, 4: paraview file, 5: 3+4)
SlipRateOutputType=1 ! 0: (smoother) slip rate output evaluated from the difference between the velocity on both side of the fault
SlipRateOutputType=0 ! 0: (smoother) slip rate output evaluated from the difference between the velocity on both side of the fault
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would not change that.

! 1: slip rate output evaluated from the fault tractions and the failure criterion (less smooth but usually more accurate where the rupture front is well developped)
/

Expand Down Expand Up @@ -58,7 +58,7 @@ PPFileName = 'tpv6_faultreceivers.dat'
/

&MeshNml
MeshFile = 'tpv5_f200m.puml.h5'
MeshFile = 'tpv5_f200_halfdom_sym'
meshgenerator = 'PUML' ! Name of meshgenerator (Gambit3D-fast, Netcdf or PUML)
/

Expand All @@ -75,8 +75,6 @@ OutputFile = 'output/tpv6'
Format = 10 ! Format (10= no output, 6=hdf5 output)
! |stress |vel
iOutputMask = 0 0 0 0 0 0 1 1 1
! |strain |eta
iPlasticityMask = 0 0 0 0 0 0 1
TimeInterval = 5. ! Index of printed info at time
refinement = 1

Expand Down