Rayleigh is a 3-D convection code designed for the study of dynamo behavior in spherical shell geometry. It evolves the incompressible and anelastic MHD equations in spherical geometry using a pseudo-spectral approach.
Rayleigh employs spherical harmonics in the horizontal direction and Chebyshev polynomials in the radial direction. The code has undergone extensive accuracy testing using the Christensen et al. (2001) Boussinesq benchmarks and the Jones et al. (2011) anelastic benchmarks. Rayleigh has been developed with NSF support through the Computational Infrastructure for Geodynamics (CIG).
Featherstone, N. (2018), Rayleigh 0.9.1, Available from: geodynamics.org, doi: 10.5281/zenodo.1236565, url: https://doi.org/10.5281/zenodo.1236565
Featherstone, N.A.; Hindman, B.W. (2016), The Spectral Amplitude Of Stellar Convection And Its Scaling In The High-Rayleigh-Number Regime, The Astrophysical Journal, 818 (1) , 32, doi: 10.3847/0004-637X/818/1/32, url: http://stacks.iop.org/0004-637X/818/i=1/a=32?key=crossref.a90f82507dd0eeb7a6e7562d1e4b0210
Please also cite: