Simulating general flows in curvilinear domains with Dedalus

Dedalus is an open-source Python framework for solving general partial differential equations at scale with modern spectral methods. Here we will describe recent additions to the code for simulating general fluid flows in curvilinear domains. These additions are based on new methods for discretizing arbitrary tensorial equations in disks, cylinders, spherical shells, and balls. These methods enable Dedalus to directly solve both incompressible and compressible hydrodynamical models with full spectral accuracy in these domains, as no reductions via scalar decompositions are necessary. We will discuss the new capabilities of the code, detail our new interface for vector-invariant model specification, and look at a range of geophysical and astrophysical applications. We will also describe ongoing development to support more complex geometries, nonlocal boundary conditions, and data-driven subgrid closures.