A hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shell
KAUST Grant NumberKUK-C1-013-04
Permanent link to this recordhttp://hdl.handle.net/10754/597291
MetadataShow full item record
AbstractA novel hybrid spectral method that combines radial basis function (RBF) and Chebyshev pseudospectral methods in a "2 + 1" approach is presented for numerically simulating thermal convection in a 3-D spherical shell. This is the first study to apply RBFs to a full 3-D physical model in spherical geometry. In addition to being spectrally accurate, RBFs are not defined in terms of any surface-based coordinate system such as spherical coordinates. As a result, when used in the lateral directions, as in this study, they completely circumvent the pole issue with the further advantage that nodes can be "scattered" over the surface of a sphere. In the radial direction, Chebyshev polynomials are used, which are also spectrally accurate and provide the necessary clustering near the boundaries to resolve boundary layers. Applications of this new hybrid methodology are given to the problem of convection in the Earth's mantle, which is modeled by a Boussinesq fluid at infinite Prandtl number. To see whether this numerical technique warrants further investigation, the study limits itself to an isoviscous mantle. Benchmark comparisons are presented with other currently used mantle convection codes for Rayleigh number (Ra) 7 × 10<sup>3</sup> and 10<sup>5</sup>. Results from a Ra = 10<sup>6</sup> simulation are also given. The algorithmic simplicity of the code (mostly due to RBFs) allows it to be written in less than 400 lines of MATLAB and run on a single workstation. We find that our method is very competitive with those currently used in the literature. Copyright 2010 by the American Geophysical Union.
CitationWright GB, Flyer N, Yuen DA (2010) A hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shell. Geochem Geophys Geosyst 11: n/a–n/a. Available: http://dx.doi.org/10.1029/2009GC002985.
SponsorsG. B. Wright was supported by NSF grants ATM-0801309 and DMS-0934581. N. Flyer was supported by NSF grants ATM-0620100 and DMS-0934317. N. Flyer's work was in part carried out while she was a Visiting Fellow at Oxford Centre for Collaborative Applied Mathematics (OCCAM) under support provided by award KUK-C1-013-04 to the University of Oxford, UK, by King Abdullah University of Science and Technology (KAUST). D. A. Yuen was supported by NSF grant DMS-0934564.
PublisherAmerican Geophysical Union (AGU)