Abstract
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 non-local 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.
Sign-in/Register to access full text options 
Free) 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
