Fourth-order time-stepping for stiff PDEs on the sphere release_lvzwma2d5fd4xdnnsnkrpz5mpe

by Hadrien Montanelli, Yuji Nakatsukasa

Released as a article .

2017  

Abstract

We present in this paper algorithms for solving stiff PDEs on the unit sphere with spectral accuracy in space and fourth-order accuracy in time. These are based on a variant of the double Fourier sphere method in coefficient space with multiplication matrices that differ from the usual ones, and implicit-explicit time-stepping schemes. Operating in coefficient space with these new matrices allows one to use a sparse direct solver, avoids the coordinate singularity and maintains smoothness at the poles, while implicit-explicit schemes circumvent severe restrictions on the time-steps due to stiffness. A comparison is made against exponential integrators and it is found that implicit-explicit schemes perform best. Implementations in MATLAB and Chebfun make it possible to compute the solution of many PDEs to high accuracy in a very convenient fashion.
In text/plain format

Archived Files and Locations

application/pdf  704.3 kB
file_ljmmwrz26fdppbbshil6p6d5tu
arxiv.org (repository)
web.archive.org (webarchive)
Read Archived PDF
Preserved and Accessible
Type  article
Stage   submitted
Date   2017-01-21
Version   v1
Language   en ?
arXiv  1701.06030v1
Work Entity
access all versions, variants, and formats of this works (eg, pre-prints)
Catalog Record
Revision: 1d4cd9ff-3a4b-4d2d-8076-8b894988b23f
API URL: JSON