Preprints
https://doi.org/10.5194/gmd-2022-252
https://doi.org/10.5194/gmd-2022-252
Submitted as: development and technical paper
10 Nov 2022
Submitted as: development and technical paper | 10 Nov 2022
Status: this preprint is currently under review for the journal GMD.

A Comparison of 3-D Spherical Shell Thermal Convection results at Low to Moderate Rayleigh Number using ASPECT (version 2.2.0) and CitcomS (version 3.3.1)

Grant Thomas Euen1, Shangxin Liu1,2, Rene Gassmöller2, Timo Heister3, and Scott David King1 Grant Thomas Euen et al.
  • 1Department of Geosciences, 926 West Campus Drive, Virginia Tech, Blacksburg, VA, USA
  • 2Department of Geological Sciences, 241 Williamson Hall, University of Florida, Gainesville, FL, USA
  • 3Mathematical and Statistical Sciences, O-110 Martin Hall, Clemson University, Clemson, SC, USA

Abstract. Due to the increasing availability of high-performance computing over the past decades numerical models have become an important tool for research in geodynamics. Several generations of mantle convection software have been developed, but due to their differing methods and increasing complexity it is important to evaluate the accuracy of each new model generation to ensure published geodynamic research is reliable and reproducible. We here explore the accuracy of the open-source, finite-element codes ASPECT and CitcomS as a function of mesh spacing using low to moderate Rayleigh number models in steady-state, thermal convection. ASPECT (Advanced Solver for Problems in Earth’s ConvecTion) is a new generation mantle convection code that enables modeling global mantle convection with realistic parameters and complicated physical processes using adaptive mesh refinement (Kronbichler et al., 2012; Heister et al., 2017). We compare the ASPECT results with calculations from the finite element code CitcomS (Zhong et al., 2000; Tan et al., 2006; Zhong et al., 2008), which has a long history of use in the geodynamics community. We find that the globally-averaged quantities: RMS velocity, mean temperature, and Nusselt number at the top and bottom of the shell, agree to within 1 %, and often much better, for calculations with sufficient mesh resolution. We also show that there is excellent agreement of the time-evolution of both the RMS velocity and the Nusselt numbers between the two codes for otherwise identical parameters. Based on our results we are optimistic that similar agreement would be achieved for calculations performed at convective vigor expected for Earth, Venus, and Mars.

Grant Thomas Euen et al.

Status: open (until 05 Jan 2023)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Grant Thomas Euen et al.

Grant Thomas Euen et al.

Viewed

Total article views: 228 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
184 38 6 228 4 2
  • HTML: 184
  • PDF: 38
  • XML: 6
  • Total: 228
  • BibTeX: 4
  • EndNote: 2
Views and downloads (calculated since 10 Nov 2022)
Cumulative views and downloads (calculated since 10 Nov 2022)

Viewed (geographical distribution)

Total article views: 214 (including HTML, PDF, and XML) Thereof 214 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 06 Dec 2022
Download
Short summary
Due to the increasing availability of high-performance computing over the past decades numerical models have become an important tool for research. Here we test two geodynamic codes that produce such models: ASPECT, a newer code, and CitcomS, an older one. We show that they produce solutions that are extremely close. As methods and codes become more complex over time, showing reproducibility allows us to seamlessly link previously-known information to modern methodologies.