Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation

R. F. Katz; M. G. Knepley; B. Smith; M. Spiegelman; E. T. Coon

doi:10.1016/j.pepi.2007.04.016

Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation

R. F. Katz, M. G. Knepley, B. Smith, M. Spiegelman, E. T. Coon

Research output: Contribution to journal › Article › peer-review

58 Scopus citations

Abstract

Geodynamics simulations are characterized by rheological nonlinearity, localization, three-dimensional effects, and separate but interacting length scales. These features represent a challenge for computational science. We discuss how a leading software framework for advanced scientific computing (the Portable Extensible Toolkit for Scientific Computation, PETSc) can facilitate the development of geodynamics simulations. To illustrate our use of PETSc, we describe simulations of (i) steady-state, non-Newtonian passive flow and thermal structure beneath a mid-ocean ridge, (ii) magmatic solitary waves in the mantle, and (iii) the formation of localized bands of high porosity in a two-phase medium being deformed under simple shear. We highlight two supplementary features of PETSc, structured storage of application parameters and self-documenting output, that are especially useful for geodynamics simulations.

Original language	English
Pages (from-to)	52-68
Number of pages	17
Journal	Physics of the Earth and Planetary Interiors
Volume	163
Issue number	1-4
DOIs	https://doi.org/10.1016/j.pepi.2007.04.016
State	Published - Aug 15 2007
Externally published	Yes

Funding

We thank Argonne National Laboratory and the Systems Group in the Mathematics and Computer Science Division for access to and management of the Blue Gene/L and Jazz systems. We are also grateful to Satish Balay for PETSc support and to Suzanne Carbotte for MOR consultations. This work was supported by a U.S. Department of Energy Computational Science Graduate Fellowship and NSF International Research Postdoctoral Fellowship grant 0602101 for Katz, as well as a U.S. Department of Energy Computational Science Graduate Fellowship for Coon. Spiegelman is supported by NSF grants 0525973 and 0610138. Contributions of Knepley and Smith are supported by the Mathematical, Information, and Computational Sciences Division subprogram of the Office of Advanced Scientific Computing Research, Office of Science, U.S. Department of Energy, under Contract DE-AC02-06CH11357.

Keywords

Geodynamics
Mid-ocean ridge
PETSc
Parallel computation
Scientific computation

Access to Document

10.1016/j.pepi.2007.04.016

Cite this

@article{cafd63d7843f4c2887dacf7d39cb6afa,

title = "Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation",

abstract = "Geodynamics simulations are characterized by rheological nonlinearity, localization, three-dimensional effects, and separate but interacting length scales. These features represent a challenge for computational science. We discuss how a leading software framework for advanced scientific computing (the Portable Extensible Toolkit for Scientific Computation, PETSc) can facilitate the development of geodynamics simulations. To illustrate our use of PETSc, we describe simulations of (i) steady-state, non-Newtonian passive flow and thermal structure beneath a mid-ocean ridge, (ii) magmatic solitary waves in the mantle, and (iii) the formation of localized bands of high porosity in a two-phase medium being deformed under simple shear. We highlight two supplementary features of PETSc, structured storage of application parameters and self-documenting output, that are especially useful for geodynamics simulations.",

keywords = "Geodynamics, Mid-ocean ridge, PETSc, Parallel computation, Scientific computation",

author = "Katz, {R. F.} and Knepley, {M. G.} and B. Smith and M. Spiegelman and Coon, {E. T.}",

year = "2007",

month = aug,

day = "15",

doi = "10.1016/j.pepi.2007.04.016",

language = "English",

volume = "163",

pages = "52--68",

journal = "Physics of the Earth and Planetary Interiors",

issn = "0031-9201",

publisher = "Elsevier",

number = "1-4",

}

TY - JOUR

T1 - Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation

AU - Katz, R. F.

AU - Knepley, M. G.

AU - Smith, B.

AU - Spiegelman, M.

AU - Coon, E. T.

PY - 2007/8/15

Y1 - 2007/8/15

N2 - Geodynamics simulations are characterized by rheological nonlinearity, localization, three-dimensional effects, and separate but interacting length scales. These features represent a challenge for computational science. We discuss how a leading software framework for advanced scientific computing (the Portable Extensible Toolkit for Scientific Computation, PETSc) can facilitate the development of geodynamics simulations. To illustrate our use of PETSc, we describe simulations of (i) steady-state, non-Newtonian passive flow and thermal structure beneath a mid-ocean ridge, (ii) magmatic solitary waves in the mantle, and (iii) the formation of localized bands of high porosity in a two-phase medium being deformed under simple shear. We highlight two supplementary features of PETSc, structured storage of application parameters and self-documenting output, that are especially useful for geodynamics simulations.

AB - Geodynamics simulations are characterized by rheological nonlinearity, localization, three-dimensional effects, and separate but interacting length scales. These features represent a challenge for computational science. We discuss how a leading software framework for advanced scientific computing (the Portable Extensible Toolkit for Scientific Computation, PETSc) can facilitate the development of geodynamics simulations. To illustrate our use of PETSc, we describe simulations of (i) steady-state, non-Newtonian passive flow and thermal structure beneath a mid-ocean ridge, (ii) magmatic solitary waves in the mantle, and (iii) the formation of localized bands of high porosity in a two-phase medium being deformed under simple shear. We highlight two supplementary features of PETSc, structured storage of application parameters and self-documenting output, that are especially useful for geodynamics simulations.

KW - Geodynamics

KW - Mid-ocean ridge

KW - PETSc

KW - Parallel computation

KW - Scientific computation

UR - http://www.scopus.com/inward/record.url?scp=34547864593&partnerID=8YFLogxK

U2 - 10.1016/j.pepi.2007.04.016

DO - 10.1016/j.pepi.2007.04.016

M3 - Article

AN - SCOPUS:34547864593

SN - 0031-9201

VL - 163

SP - 52

EP - 68

JO - Physics of the Earth and Planetary Interiors

JF - Physics of the Earth and Planetary Interiors

IS - 1-4

ER -

Numerical simulation of geodynamic processes with the Portable Extensible Toolkit for Scientific Computation

Abstract

Funding

Keywords

Access to Document

Other files and links

Fingerprint

Cite this