An efficient spectral Poisson solver for the nirvana-III code: the shearing-box case with vertical vacuum boundary conditions (2510.10070v1)

Abstract: The stability of a differentially rotating fluid subject to its own gravity is a problem with applications across wide areas of astrophysics--from protoplanetary discs (PPDs) to entire galaxies. The shearing box formalism offers a conceptually simple framework for studying differential rotation in the local approximation. Aimed at self-gravitating, and importantly, vertically stratified PPDs, we develop two novel methods for solving Poisson's equation in the framework of the shearing box with vertical vacuum boundary conditions (BCs). Both approaches naturally make use of multi-dimensional fast Fourier transforms for computational efficiency. While the first one exploits the linearity properties of the Poisson equation, the second, which is slightly more accurate, consists of finding the adequate discrete Green's function (in Fourier space) adapted to the problem at hand. To this end, we have revisited the method proposed by Vico et al. (2016) and have derived an analytical Green's function satisfying the shear-periodic BCs in the plane as well as vacuum BCs, vertically. Our spectral method demonstrates excellent accuracy, even with a modest number of grid points, and exhibits third-order convergence. It has been implemented in the NIRVANA-III code, where it exhibits good scalability up to 4096 CPU cores, consuming less than 6% of the total runtime. This was achieved through the use of P3DFFT, a fast Fourier Transform library that employs pencil decomposition, overcoming the scalability limitations inherent in libraries using slab decomposition. We have introduced two novel spectral Poisson solvers that guarantees high accuracy, performance, and intrinsically support vertical vacuum boundary conditions in the shearing-box framework. Our solvers enable high-resolution local studies involving self-gravity, such as MHD simulations of gravito-turbulence or gravitational fragmentation.