A hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shell

Handle URI:
http://hdl.handle.net/10754/597291
Title:
A hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shell
Authors:
Wright, G. B.; Flyer, N.; Yuen, D. A.
Abstract:
A 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.
Citation:
Wright 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.
Publisher:
Wiley-Blackwell
Journal:
Geochemistry, Geophysics, Geosystems
KAUST Grant Number:
KUK-C1-013-04
Issue Date:
Jul-2010
DOI:
10.1029/2009GC002985
Type:
Article
ISSN:
1525-2027
Sponsors:
G. 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.
Appears in Collections:
Publications Acknowledging KAUST Support

Full metadata record

DC FieldValue Language
dc.contributor.authorWright, G. B.en
dc.contributor.authorFlyer, N.en
dc.contributor.authorYuen, D. A.en
dc.date.accessioned2016-02-25T12:29:58Zen
dc.date.available2016-02-25T12:29:58Zen
dc.date.issued2010-07en
dc.identifier.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.en
dc.identifier.issn1525-2027en
dc.identifier.doi10.1029/2009GC002985en
dc.identifier.urihttp://hdl.handle.net/10754/597291en
dc.description.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.en
dc.description.sponsorshipG. 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.en
dc.publisherWiley-Blackwellen
dc.subjectIsoviscousen
dc.subjectMantle convectionen
dc.subjectPseudospectralen
dc.subjectRadial basis functionen
dc.subjectSpherical geometryen
dc.titleA hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shellen
dc.typeArticleen
dc.identifier.journalGeochemistry, Geophysics, Geosystemsen
dc.contributor.institutionDepartment of Mathematics; Boise State University; Boise; Idaho; 83725; USAen
dc.contributor.institutionInstitute for Mathematics Applied to Geosciences; National Center for Atmospheric Research; Boulder; Colorado; 80305; USAen
dc.contributor.institutionDepartment of Geology and Geophysics; University of Minnesota; Minneapolis; Minnesota; 55455; USAen
kaust.grant.numberKUK-C1-013-04en
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.