MARS: A method for the adaptive removal of stiffness in PDEs

Laurent Duchemin*, Jens Eggers

*Corresponding author for this work

Research output: Contribution to journalArticle (Academic Journal)peer-review


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.
Original languageEnglish
Article number11624
Pages (from-to)1-14
Number of pages14
JournalJournal of Computational Physics
Early online date24 Sept 2022
Publication statusPublished - 15 Dec 2022


Dive into the research topics of 'MARS: A method for the adaptive removal of stiffness in PDEs'. Together they form a unique fingerprint.

Cite this