Merge pull request #58 from SeisSol/fabian/ttpv-benchmark-quakeworx

Merge PR after incorporating feedback from Thomas and Alice.

Author: fabian-kutschera

.DS_Store

@@ -1 +1,3 @@
 *.ipynb_checkpoints/
+.DS_Store
+../.DS_Store

.gitignore

@@ -0,0 +1,45 @@
+!Switch
+[mu_s, mu_d, d_c]: !ConstantMap
+      map:
+        mu_s:        0.575
+        mu_d:        0.450
+        d_c:         0.18
+[forced_rupture_time, cohesion, T_n, T_s, T_d]: !LuaMap
+  returns: [forced_rupture_time, cohesion, T_n, T_s, T_d]
+  function: |
+     function f (x)
+      xh = 17387.0
+      yh = 0.0
+      zh = -4658.0
+      r_crit = 4000.0
+      r = math.sqrt((x["x"]-xh)^2 + (x["y"]-yh)^2 + (x["z"]-zh)^2 )
+      Vs = 3464.0
+      if (r <= r_crit) then
+        forced_rupture_time = r/(0.7*Vs)+(0.081*r_crit/(0.7*Vs))*(1.0/(1.0-(r/r_crit)^2)-1.0)
+      else 
+        forced_rupture_time = 1000000000.0
+      end
+
+      pi = 4.0 * math.atan(1.0)
+      depth = math.abs(x["z"])
+      w = depth/math.sin(15.0*pi/180.)
+
+      if (w < 8000) then
+        cohesion = -500.0 * math.max(8000.0-w,0)
+      else 
+        cohesion = 0.0
+      end
+      
+      return {
+        forced_rupture_time = forced_rupture_time,
+        cohesion = cohesion,
+        T_n = -4240.0 * w,
+        T_s = 0.0,
+        T_d = 2120.0 * w,
+      }
+     end
+[Tnuc_n, Tnuc_s, Tnuc_d]: !ConstantMap
+    map:
+        Tnuc_n: 0
+        Tnuc_s: 0
+        Tnuc_d: 0

quakeworx/earthquake-tsunami/fault.yaml

@@ -0,0 +1,17 @@
+!Switch
+[rho, mu, lambda]: !Any
+  components:
+    - !GroupFilter
+      groups: [0]
+      components: !ConstantMap
+        map: 
+          rho:    1000
+          mu:     0.0
+          lambda: 2.250e9
+    - !GroupFilter
+      groups: [1]
+      components: !ConstantMap
+        map:
+          rho: 2670
+          mu: 32038120320
+          lambda: 32043759360

quakeworx/earthquake-tsunami/figures/Figure01.png

@@ -0,0 +1,99 @@
+&equations
+MaterialFileName = 'material_WL.yaml'
+Plasticity = 0
+!numflux = 'godunov'
+!numfluxnearfault = 'godunov'
+/
+
+&IniCondition
+/
+
+&DynamicRupture
+FL = 16                                        ! Friction law
+ModelFileName = 'fault.yaml'
+
+!characteristic time of the forced nucleation
+t_0 = 0.5
+
+!reference vector for defining strike and dip direction
+XRef = -0.1
+YRef = 0.0
+ZRef = -1.0
+refPointMethod = 1
+
+OutputPointType = 4                            ! 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
+                            ! 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)
+/
+
+!see: https://seissol.readthedocs.io/en/latest/fault-output.html
+!parameterize paraview file output
+&Elementwise
+printtimeinterval_sec = 1.0                    ! Time interval at which output will be written
+OutputMask = 1 1 1 1 1 1 1 1 1 1 1             ! turn on and off fault outputs
+refinement_strategy = 2
+refinement = 1
+vtkorder = -1 ! positive values enable the high-order VTKHDF format for the elementwise fault output
+/
+
+! parameterize ascii fault file outputs
+&Pickpoint
+printtimeinterval = 100       ! Index of printed info at timesteps
+OutputMask = 1 1 1 1 1 1 1 1 1 1 1 1        ! turn on and off fault outputs
+!pickdt = 0.005 ! not yet implemented
+PPFileName = 'onfault_receivers.dat'
+/
+
+&SourceType
+/
+            
+&MeshNml
+MeshFile = 'mesh_ttpv1_small_18k.puml.h5'
+!MeshFile = 'mesh_ttpv2_225k.puml.h5'
+meshgenerator = 'PUML'                         ! Name of meshgenerator (Gambit3D-fast, Netcdf or PUML)
+pumlboundaryformat='i32'                       ! the boundary data type for PUML files ('auto', 'i32', 'i64', 'i32x4')
+/
+
+&Discretization
+CFL = 0.4                                      ! CFL number (<=1.0)
+ClusteredLTS = 2                               ! Enable local time-stepping
+/
+
+&Output
+OutputFile = 'outputs/ttpv1'
+!OutputFile = 'outputs/ttpv2' 
+WavefieldOutput = 0        ! disable/enable volumne wavefield output (right now, format=6 needs to be set as well)
+Format = 10                                     ! Format (10=no output, 6=hdf5 output)
+!             |stress     |vel
+iOutputMask = 1 1 1 1 1 1 1 1 1
+TimeInterval = 1.0                           ! Index of printed info at time
+refinement = 1
+
+! Free surface output
+SurfaceOutput = 1
+SurfaceOutputRefinement = 1
+SurfaceOutputInterval = 1.0
+
+ReceiverOutput = 0 !1
+pickdt = 0.01                       ! Pickpoint Sampling
+RFileName = ''      ! Record Points in extra file
+ReceiverOutputInterval = 10.0
+ReceiverComputeStrain = 0
+
+xdmfWriterBackend = 'hdf5' ! (optional) The backend used in fault, wavefield,
+! and free-surface output. The HDF5 backend is only supported when SeisSol is compiled with
+! HDF5 support.
+
+EnergyOutput = 1 ! Computation of energy, written in csv file
+EnergyTerminalOutput = 1 ! Write energy to standard output
+EnergyOutputInterval = 1.0
+
+wavefieldvtkorder = -1
+surfacevtkorder = -1
+
+Checkpoint = 0
+/
+
+&AbortCriteria
+EndTime = 240.0
+/