The E(xplicit)I(implicit)N(null) method was developed recently to remove numerical instability from PDEs, adding and subtracting an operator D of arbitrary structure, treating the operator implicitly in one case, and explicitly in the other. Here we extend this idea by devising an adaptive procedure to find an optimal approximation for D. We propose a measure of the numerical error which detects numerical instabilities across all wavelengths, and adjust each Fourier component of D to the smallest value such that numerical instability is suppressed. We show that for a number of nonlinear and nonlocal PDEs, in one and two dimensions, the spectrum of D adapts automatically and dynamically to the theoretical result for marginal stability. The adaptive implicit part is diagonal in Fourier space, so that our method has the same stability properties as a fully implicit method, with minimal computational overhead coming only from performing the fast Fourier transform.
|Number of pages||14|
|Journal||Journal of Computational Physics|
|Early online date||24 Sep 2022|
|Publication status||Published - 15 Dec 2022|