Dedalus: A flexible framework for numerical simulations with spectral methods
Top Cited Papers
Open Access
- 23 April 2020
- journal article
- research article
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 2 (2) , 023068
Abstract
Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical method that enables this translation, describes the design and implementation of the codebase, and illustrates its capabilities with a variety of example problems. The numerical method is a first-order generalized tau formulation that discretizes equations into banded matrices. This method is implemented with an object-oriented design. Classes for spectral bases and domains manage the discretization and automatic parallel distribution of variables. Discretized fields and mathematical operators are symbolically manipulated with a basic computer algebra system. Initial value, boundary value, and eigenvalue problems are efficiently solved using high-performance linear algebra, transform, and parallel communication libraries. Custom analysis outputs can also be specified in plain text and stored in self-describing portable formats. The performance of the code is evaluated with a parallel scaling benchmark and a comparison to a finite-volume code. The features and flexibility of the codebase are illustrated by solving several examples: the nonlinear Schrödinger equation on a graph, a supersonic magnetohydrodynamic vortex, quasigeostrophic flow, Stokes flow in a cylindrical annulus, normal modes of a radiative atmosphere, and diamagnetic levitation.Keywords
All Related Versions
Funding Information
- National Science Foundation (NSF PHY-1748958, 1122374, AST10-09802, DGE 1106400)
- University of Sydney
- U.S. Department of Energy
- Australian Research Council (DE140101960)
- University of Colorado
- National Aeronautics and Space Administration (NNX16AC92G, 80NSSC19K0026)
This publication has 122 references indexed in Scilit:
- HIGH-LUNDQUIST NUMBER SCALING IN THREE-DIMENSIONAL SIMULATIONS OF PARKER'S MODEL OF CORONAL HEATINGThe Astrophysical Journal, 2012
- Charney isotropy and equipartition in quasi-geostrophic turbulenceJournal of Fluid Mechanics, 2010
- A numerical comparison of Chebyshev methods for solving fourth order semilinear initial boundary value problemsJournal of Computational and Applied Mathematics, 2010
- Magnetic Buoyancy Instabilities of a Shear‐generated Magnetic LayerThe Astrophysical Journal, 2008
- Generalized quasi-geostrophy for spatially anisotropic rotationally constrained flowsJournal of Fluid Mechanics, 2006
- The effects of a hot outer atmosphere on acoustic-gravity wavesThe Astrophysical Journal, 1994
- Anisotropy and Coherent Vortex Structures in Planetary TurbulenceScience, 1994
- The equations of reduced magnetohydrodynamicsJournal of Plasma Physics, 1992
- The numerical solution of linear differential equations in Chebyshev seriesMathematical Proceedings of the Cambridge Philosophical Society, 1957
- The general circulation of the atmosphere: A numerical experimentQuarterly Journal of the Royal Meteorological Society, 1956