GeoFramework: Example of using frameworks in geophysics

1
GeoFramework Example of using frameworks in
geophysics
Mike Gurnis Caltech

• Tectonic Modeling Workshop
• June, 2005

2
A hierarchy of software components
?
Superstructure
?
Geodynamic Specific
?
Infrastructure
?
Library
3
Superstructure and Geodynamic Specific Layers

• Simulation controller
• Monitor Simulation
• Couple Fluid to Solid

Superstructure

• Component A
• Component B
• .

Geodynamic Specific
4
Pyre Framework (Science Neutral)

• The integration framework is a set of
  co-operating abstract services

python
package
abstract class
specialization
FORTRAN/C/C
Developed by Michael Aivazis
5
GeoFramework TECTON Solver
6
Virtual Test Facility (VTF) Example From Pyre
7
CitcomS
main
Parser
Mesher
T solver
Stokes solver
Output
8
Pyre version of CitcomS
regionalcitcoms.py
Components (Python)
Python bindings
Mesher
T solver
Stokes solver
Output
RegionalModule.so (a shared object)
9
CitcomS.py as a single component
10
Coupling With Pyre
11
Regional and Global Mantle Flow Coupled with Pyre
CitcomS.py, Eh Tan
12
Regional CitcomS coupled to full CitcomS
CitcomS.py, Eh Tan
13
Predictions Hot-spot track synthetic tomography
14
Example Mantle convection coupled to lithosphere

• Monitor Simulation
• Couple Fluid to Solid
• Visualization

Superstructure

• Self-contained geophysics (single phyiscs)
• Rheology modules

Geodynamic Specific

• Mesher Solid Fluid
• Solver Solid Fluid

Infrastructure
Library PETSc, BLAS, MPI
15
Example of Geodynamic Specific Infrastructure
Layers

• SNARK -- particle based FEM with implicit solver
• SNAC -- Lagragian explicit FEM

Geodynamic Specific

• StGermain A framework with entry points
  plugins for building meshes, advecting
  particles, calling solvers, I/O..

Infrastructure
Library PETSc, MPI
16
Example of Geodynamic Specific Infrastructure
Layers
Geodynamic Specific SNARK Monash group
(Moresi) SNAC Caltech/Texas groups (Gurnis,
Lavier)

• SNARK -- particle based FEM with implicit solver
• SNAC -- Lagragian explicit FEM

Infrastructure Software Engineers at
the Victorian Partnership for Advanced
Computing (Steve Quenette Team Leader)

• StGermain A framework with entry points
  plugins for building meshes, advecting
  particles, calling solvers, I/O..

17
SNARK SNAC
Moresi et al.
Choi et al.
StGermain A framework with entry points
plugins for building meshes, advecting
particles, calling solvers, I/O..
Quenette et al.
18
SNAC
19
Examples of coupling codes with Pyre
(superstructure framework) GeoFramework
Pyre
CitcomS
SNAC
pHMelts
a geophysics solver
Exchanger
20
SNAC CitcomS coupling (Crust-Mantle Interaction)
Eun-seo Choi et al.