Hybrid Modeling For Fault Rupture and Wave Propagation ' Metadata and other religious issues

1
Hybrid Modeling For Fault Rupture and Wave
Propagation.Metadata and other religious
issues
Hybrid Numerical Modeling for Fault Rupture and
Wave Propagation

• Geoffrey Ely
• Bernard Minster
• Steve Day
• SCEC ITR AHM
• Jan 26 2004

2
Computational Pathways
3
(No Transcript)
4
(No Transcript)
5
Teragrid AWM Simulation

• Output Statistics
• Velocity vector field
• 3000 time steps
• 3GB per time step
• Stress tensor field
• 300 time steps
• 6GB per time step
• 10 Terabytes total storage
• 1.8 million files
• Transferred from GPFS to the SRB in Aug 03
• Transferred at 6Mb/s 3 weeks
• Current performance 100Mb/s (?) a few hours

• Computation
• Parallelization
• Domain decomposition using MPI
• Run on SDSC Blue Horizon
• 512 Processors
• 12 hours
• 6000 CPU Hours
• Output written to GPFS parallel file system
• Throughput - 230 Mb/s

6
Hybrid Earthquake Model
Rupture Dynamics Model (RDM) performs fault plane
and free surface calculations. Handles irregular
geometry needed for including realistic
topography and dipping faults. Anelastic Wave
Model (AWM) used away from the boundaries were
the grid can be made regular. Is much more
efficient than the RDM and uses higher order
(4th) approximations. Also includes anelastic
corrections important for accurate simulation of
wave traveling long distances.
Fault trace
7
Coupling Scheme Interpolate and exchange
velocity field variables at the boundary
RDM/AWM boundary
AWM
RDM
Vx
Vy
Vz
RDM V interpolated from AWM Vx, Vy Vz
V
AWM Vx, Vy Vz interpolated from RDM V
8
(No Transcript)
9
Support Operator Method

• Structured mesh
• Interpolate mesh node locations and function
  values with Lagrange polynomials
• Differentiate and evaluate at cell centers to
  derive prime operator (e.g. grad)
• Use a discrete analog of this integral identity
  to derive the support operator (e.g. div)div
  is the negative adjoint of grad

10
Animations
11
Metadata Standards Structure
12
SCEC Simulation Metadata
ISO 19107
13
Metadata schemadesign criteria

• Searchable through SQL techniques
• Machine/program interface fast searches
• Keep users happyintuitive organization
• No duplication of information. Maintain
  integrity
• Functionally complete description of metadata
• Use relational schema, try for DBMS independence
  (e.g. Oracle, MySQL, SyBase, Postgress)
• Maintain SRB compatibility
• Basic hub object The Run

14
SCEC Simulation Metadata
ISO 19107
15
(No Transcript)
16
Tables
17
Formats. This is the religious part.

• Recommend that seismograms be translated on
  demand into relevant popular formats (e.g.
  COSMOS, SAC, SEED, SEG-Y, CSS, etc) This could
  be done through a web service at user request.
• Same for grid snapshots. Preliminary web service
  to translate grid snapshots from KOF format to
  HDF-5
• Can further down-translate to HDF-4

18
Conversion service