Open AccessA community benchmark for 2‐D Cartesian compressible convection in the Earth's mantle
Author(s) -
King Scott D.,
Lee Changyeol,
Van Keken Peter E.,
Leng Wei,
Zhong Shijie,
Tan Eh,
Tosi Nicola,
Kameyama Masanori C.
Publication year - 2010
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2009.04413.x
Subject(s) - mantle convection , nusselt number , discretization , compressibility , convection , mechanics , rayleigh number , dissipation , physics , natural convection , mathematics , cartesian coordinate system , viscosity , geology , geophysics , thermodynamics , mathematical analysis , geometry , reynolds number , seismology , tectonics , lithosphere , turbulence
SUMMARY Benchmark comparisons are an essential tool to verify the accuracy and validity of computational approaches to mantle convection. Six 2‐D Cartesian compressible convection codes are compared for steady‐state constant and temperature‐dependent viscosity cases as well as time‐dependent constant viscosity cases. In general we find good agreement between all codes when comparing average flow characteristics such as Nusselt number and rms velocity. At Rayleigh numbers near 10 6 and dissipation numbers between 0 and 2, the results differ by approximately 1 per cent. Differences in discretization and use of finite volumes versus finite elements dominate the differences. There is a small systematic difference between the use of the anelastic liquid approximation (ALA) compared to that of the truncated ALA. In determining the onset of time‐dependence, there was less agreement between the codes with a spread in the Rayleigh number where the first bifurcation occurs ranging from 7.79 × 10 5 to 1.05 × 10 6 .