MLegS: A modern mapped Legendre spectral method solver for unbounded domains with parallelization
Abstract
MLegS is an open-source Modern Fortran code package for solving linear and nonlinear time-dependent partial differential equations on radially unbounded domains. It builds on an earlier in-house prototype, recasting it as a modular solver with a unified interface for fields, transforms, and spectral operators. Within this framework, the numerical formulation maps 0≤ r<∞ onto a finite interval, where associated Legendre functions provide the radial basis, while Fourier series represent the periodic azimuthal and axial directions. The mapped radial basis represents radial infinity without imposing an artificial boundary condition at finite radius while incorporating the regularity required at r=0 according to azimuthal order. The implementation distributes the corresponding field layouts using the Message-Passing Interface (MPI), with strong- and weak-scaling measurements extending to 512 cores across four compute nodes. Beyond parallel performance, three examples assess the solver across diffusion, reaction-diffusion, and vortex dynamics: closed-form diffusion solutions establish spectral spatial convergence and the design orders of both semi-implicit integrators; a Fisher-Kolmogorov-Petrovsky-Piskunov front matches the linearized solution and propagates at the classical speed within its logarithmic correction; and vortex-pair simulations based on toroidal-poloidal decomposition yield a measured Crow growth rate 2.7% above the theoretical prediction and demonstrate the reconnection of the vortex pair.
Related
- Essential spectral geometry of the Maxwell system in unbounded domains
- On Simply Connected Domains Supporting an Unbounded Analytic Function with Bounded Derivative
- Global Well-Posedness for Planar Compressible Magnetohydrodynamics with Degenerate Viscosity on Unbounded Domains
- On the geometry of unbounded wandering domains
- Global Existence and Stability of 3D Stochastic NSEs in Bounded and Unbounded Domains Driven by a Special Multiplicative Wiener Process
BibTeX
@misc{lee2026mlegs,
title = {{MLegS: A modern mapped Legendre spectral method solver for unbounded domains with parallelization}},
author = {Sangjoon Lee and Jinge Wang},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.18026},
eprint = {2609.18026},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.18026},
}
RIS
TY - UNPB
TI - MLegS: A modern mapped Legendre spectral method solver for unbounded domains with parallelization
AU - Sangjoon Lee
AU - Jinge Wang
PY - 2026
DA - 2026/09/16
PB - arXiv
DO - 10.48550/arxiv.2609.18026
UR - https://arxiv.org/abs/2609.18026
AB - MLegS is an open-source Modern Fortran code package for solving linear and nonlinear time-dependent partial differential equations on radially unbounded domains. It builds on an earlier in-house prototype, recasting it as a modular solver with a unified interface for fields, transforms, and spectral operators. Within this framework, the numerical formulation maps 0≤ r<∞ onto a finite interval, where associated Legendre functions provide the radial basis, while Fourier series represent the periodic azimuthal and axial directions. The mapped radial basis represents radial infinity without imposing an artificial boundary condition at finite radius while incorporating the regularity required at r=0 according to azimuthal order. The implementation distributes the corresponding field layouts using the Message-Passing Interface (MPI), with strong- and weak-scaling measurements extending to 512 cores across four compute nodes. Beyond parallel performance, three examples assess the solver across diffusion, reaction-diffusion, and vortex dynamics: closed-form diffusion solutions establish spectral spatial convergence and the design orders of both semi-implicit integrators; a Fisher-Kolmogorov-Petrovsky-Piskunov front matches the linearized solution and propagates at the classical speed within its logarithmic correction; and vortex-pair simulations based on toroidal-poloidal decomposition yield a measured Crow growth rate 2.7% above the theoretical prediction and demonstrate the reconnection of the vortex pair.
ER -