“…PDE solvers based on such an approach maintain the desirable qualities of spectrally-accurate numerical solvers for time-dependent systems: they can produce accurate solutions by means of relatively coarse discretizations and, importantly, they carry minimal numerical diffusion or dispersion errors (manifesting as artificial amplitude decay or period elongation). FC-based PDE solvers have been successfully constructed for a variety of physical equations including classical wave equations [17,18], non-linear Burgers systems [19], Euler equations [20,21], Navier-Stokes equations [22,23,24], radiative transfer equations [25], Navier elastodynamics equations [26,27] and 1D fluid-structure hemodynamics equations [28]. The numerical methodology introduced in this work represents a fully-realized 3D convection-diffusion solver based on this Fast Fourier Transform (FFT)-speed Fourier continuation approximation procedure for both Dirichlet and Neumann boundary conditions.…”