A benchmark study on mantle convection in a 3-D spherical shell using CitcomS

Shijie Zhong, Allen McNamara, Eh Tan, Louis Moresi, Michael Gurnis

Research output: Contribution to journalArticle

136 Citations (Scopus)

Abstract

As high-performance computing facilities and sophisticated modeling software become available, modeling mantle convection in a three-dimensional (3-D) spherical shell geometry with realistic physical parameters and processes becomes increasingly feasible. However, there is still a lack of comprehensive benchmark studies for 3-D spherical mantle convection. Here we present benchmark and test calculations using a finite element code CitcomS for 3-D spherical convection. Two classes of model calculations are presented: the Stokes' flow and thermal and thermochemical convection. For Stokes' flow, response functions of characteristic flow velocity, topography, and geoid at the surface and core-mantle boundary (CMB) at different spherical harmonic degrees are computed using CitcomS and are compared with those from analytic solutions using a propagator matrix method. For thermal and thermochemical convection, 24 cases are computed with different model parameters including Rayleigh number (7 × 103 or 105) and viscosity contrast due to temperature dependence (1 to 107). For each case, time-averaged quantities at the steady state are computed, including surface and CMB Nussult numbers, RMS velocity, averaged temperature, and maximum and minimum flow velocity, and temperature at the midmantle depth and their standard deviations. For thermochemical convection cases, in addition to outputs for thermal convection, we also quantified entrainment of an initially dense component of the convection and the relative errors in conserving its volume. For nine thermal convection cases that have small viscosity variations and where previously published results were available, we find that the CitcomS results are mostly consistent with these previously published with less than 1% relative differences in globally averaged quantities including Nussult numbers and RMS velocities. For other 15 cases with either strongly temperature-dependent viscosity or thermochemical convection, no previous calculations are available for comparison, but these 15 test calculations from CitcomS are useful for future code developments and comparisons. We also presented results for parallel efficiency for CitcomS, showing that the code achieves 57% efficiency with 3072 cores on Texas Advanced Computing Center's parallel supercomputer Ranger.

Original languageEnglish (US)
Article numberQ10017
JournalGeochemistry, Geophysics, Geosystems
Volume9
Issue number10
DOIs
StatePublished - Oct 2008

Fingerprint

mantle convection
spherical shells
Earth mantle
convection
shell
free convection
thermal convection
viscosity
core-mantle boundary
flow velocity
Stokes flow
temperature
Viscosity
Rayleigh number
Flow velocity
spherical harmonics
geoid
entrainment
supercomputers
modeling

Keywords

  • Benchmark
  • Mantle convection

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geophysics

Cite this

A benchmark study on mantle convection in a 3-D spherical shell using CitcomS. / Zhong, Shijie; McNamara, Allen; Tan, Eh; Moresi, Louis; Gurnis, Michael.

In: Geochemistry, Geophysics, Geosystems, Vol. 9, No. 10, Q10017, 10.2008.

Research output: Contribution to journalArticle

Zhong, Shijie ; McNamara, Allen ; Tan, Eh ; Moresi, Louis ; Gurnis, Michael. / A benchmark study on mantle convection in a 3-D spherical shell using CitcomS. In: Geochemistry, Geophysics, Geosystems. 2008 ; Vol. 9, No. 10.
@article{44649bffd98d4536a444f2d74bb3c126,
title = "A benchmark study on mantle convection in a 3-D spherical shell using CitcomS",
abstract = "As high-performance computing facilities and sophisticated modeling software become available, modeling mantle convection in a three-dimensional (3-D) spherical shell geometry with realistic physical parameters and processes becomes increasingly feasible. However, there is still a lack of comprehensive benchmark studies for 3-D spherical mantle convection. Here we present benchmark and test calculations using a finite element code CitcomS for 3-D spherical convection. Two classes of model calculations are presented: the Stokes' flow and thermal and thermochemical convection. For Stokes' flow, response functions of characteristic flow velocity, topography, and geoid at the surface and core-mantle boundary (CMB) at different spherical harmonic degrees are computed using CitcomS and are compared with those from analytic solutions using a propagator matrix method. For thermal and thermochemical convection, 24 cases are computed with different model parameters including Rayleigh number (7 × 103 or 105) and viscosity contrast due to temperature dependence (1 to 107). For each case, time-averaged quantities at the steady state are computed, including surface and CMB Nussult numbers, RMS velocity, averaged temperature, and maximum and minimum flow velocity, and temperature at the midmantle depth and their standard deviations. For thermochemical convection cases, in addition to outputs for thermal convection, we also quantified entrainment of an initially dense component of the convection and the relative errors in conserving its volume. For nine thermal convection cases that have small viscosity variations and where previously published results were available, we find that the CitcomS results are mostly consistent with these previously published with less than 1{\%} relative differences in globally averaged quantities including Nussult numbers and RMS velocities. For other 15 cases with either strongly temperature-dependent viscosity or thermochemical convection, no previous calculations are available for comparison, but these 15 test calculations from CitcomS are useful for future code developments and comparisons. We also presented results for parallel efficiency for CitcomS, showing that the code achieves 57{\%} efficiency with 3072 cores on Texas Advanced Computing Center's parallel supercomputer Ranger.",
keywords = "Benchmark, Mantle convection",
author = "Shijie Zhong and Allen McNamara and Eh Tan and Louis Moresi and Michael Gurnis",
year = "2008",
month = "10",
doi = "10.1029/2008GC002048",
language = "English (US)",
volume = "9",
journal = "Geochemistry, Geophysics, Geosystems",
issn = "1525-2027",
publisher = "American Geophysical Union",
number = "10",

}

TY - JOUR

T1 - A benchmark study on mantle convection in a 3-D spherical shell using CitcomS

AU - Zhong, Shijie

AU - McNamara, Allen

AU - Tan, Eh

AU - Moresi, Louis

AU - Gurnis, Michael

PY - 2008/10

Y1 - 2008/10

N2 - As high-performance computing facilities and sophisticated modeling software become available, modeling mantle convection in a three-dimensional (3-D) spherical shell geometry with realistic physical parameters and processes becomes increasingly feasible. However, there is still a lack of comprehensive benchmark studies for 3-D spherical mantle convection. Here we present benchmark and test calculations using a finite element code CitcomS for 3-D spherical convection. Two classes of model calculations are presented: the Stokes' flow and thermal and thermochemical convection. For Stokes' flow, response functions of characteristic flow velocity, topography, and geoid at the surface and core-mantle boundary (CMB) at different spherical harmonic degrees are computed using CitcomS and are compared with those from analytic solutions using a propagator matrix method. For thermal and thermochemical convection, 24 cases are computed with different model parameters including Rayleigh number (7 × 103 or 105) and viscosity contrast due to temperature dependence (1 to 107). For each case, time-averaged quantities at the steady state are computed, including surface and CMB Nussult numbers, RMS velocity, averaged temperature, and maximum and minimum flow velocity, and temperature at the midmantle depth and their standard deviations. For thermochemical convection cases, in addition to outputs for thermal convection, we also quantified entrainment of an initially dense component of the convection and the relative errors in conserving its volume. For nine thermal convection cases that have small viscosity variations and where previously published results were available, we find that the CitcomS results are mostly consistent with these previously published with less than 1% relative differences in globally averaged quantities including Nussult numbers and RMS velocities. For other 15 cases with either strongly temperature-dependent viscosity or thermochemical convection, no previous calculations are available for comparison, but these 15 test calculations from CitcomS are useful for future code developments and comparisons. We also presented results for parallel efficiency for CitcomS, showing that the code achieves 57% efficiency with 3072 cores on Texas Advanced Computing Center's parallel supercomputer Ranger.

AB - As high-performance computing facilities and sophisticated modeling software become available, modeling mantle convection in a three-dimensional (3-D) spherical shell geometry with realistic physical parameters and processes becomes increasingly feasible. However, there is still a lack of comprehensive benchmark studies for 3-D spherical mantle convection. Here we present benchmark and test calculations using a finite element code CitcomS for 3-D spherical convection. Two classes of model calculations are presented: the Stokes' flow and thermal and thermochemical convection. For Stokes' flow, response functions of characteristic flow velocity, topography, and geoid at the surface and core-mantle boundary (CMB) at different spherical harmonic degrees are computed using CitcomS and are compared with those from analytic solutions using a propagator matrix method. For thermal and thermochemical convection, 24 cases are computed with different model parameters including Rayleigh number (7 × 103 or 105) and viscosity contrast due to temperature dependence (1 to 107). For each case, time-averaged quantities at the steady state are computed, including surface and CMB Nussult numbers, RMS velocity, averaged temperature, and maximum and minimum flow velocity, and temperature at the midmantle depth and their standard deviations. For thermochemical convection cases, in addition to outputs for thermal convection, we also quantified entrainment of an initially dense component of the convection and the relative errors in conserving its volume. For nine thermal convection cases that have small viscosity variations and where previously published results were available, we find that the CitcomS results are mostly consistent with these previously published with less than 1% relative differences in globally averaged quantities including Nussult numbers and RMS velocities. For other 15 cases with either strongly temperature-dependent viscosity or thermochemical convection, no previous calculations are available for comparison, but these 15 test calculations from CitcomS are useful for future code developments and comparisons. We also presented results for parallel efficiency for CitcomS, showing that the code achieves 57% efficiency with 3072 cores on Texas Advanced Computing Center's parallel supercomputer Ranger.

KW - Benchmark

KW - Mantle convection

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

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

U2 - 10.1029/2008GC002048

DO - 10.1029/2008GC002048

M3 - Article

VL - 9

JO - Geochemistry, Geophysics, Geosystems

JF - Geochemistry, Geophysics, Geosystems

SN - 1525-2027

IS - 10

M1 - Q10017

ER -