We implement a second-order exponential integrator for semidiscretized advection-diffusion-reaction equations, obtained by coupling exponential-like Euler and Midpoint integrators, and computing the relevant matrix exponentials by polynomial interpolation at Leja points. Numerical tests on 2D models discretized in space by finite differences or finite elements, show that the Leja-Euler-Midpoint (LEM) exponential integrator can be up to 5 times faster than a classical second-order implicit solver.
The LEM exponential integrator for advection-diffusion-reaction equations
BERGAMASCHI, LUCA;VIANELLO, MARCO
2007
Abstract
We implement a second-order exponential integrator for semidiscretized advection-diffusion-reaction equations, obtained by coupling exponential-like Euler and Midpoint integrators, and computing the relevant matrix exponentials by polynomial interpolation at Leja points. Numerical tests on 2D models discretized in space by finite differences or finite elements, show that the Leja-Euler-Midpoint (LEM) exponential integrator can be up to 5 times faster than a classical second-order implicit solver.File in questo prodotto:
File | Dimensione | Formato | |
---|---|---|---|
cvb07.pdf
accesso aperto
Tipologia:
Preprint (submitted version)
Licenza:
Accesso libero
Dimensione
641.35 kB
Formato
Adobe PDF
|
641.35 kB | Adobe PDF | Visualizza/Apri |
1-s2.0-S0377042706006479-main.pdf
non disponibili
Tipologia:
Published (publisher's version)
Licenza:
Accesso privato - non pubblico
Dimensione
582.32 kB
Formato
Adobe PDF
|
582.32 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.