z-logo
open-access-imgOpen Access
A semi‐Lagrangian Crank‐Nicolson algorithm for the numerical solution of advection‐diffusion problems
Author(s) -
Spiegelman Marc,
Katz Richard F.
Publication year - 2006
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2005gc001073
Subject(s) - advection , interpolation (computer graphics) , numerical diffusion , mathematics , péclet number , diffusion , operator (biology) , limit (mathematics) , range (aeronautics) , numerical analysis , algorithm , computer science , mathematical analysis , physics , mechanics , thermodynamics , animation , biochemistry , chemistry , computer graphics (images) , materials science , repressor , transcription factor , gene , composite material
We present a hybrid method for the numerical solution of advection‐diffusion problems that combines two standard algorithms: semi‐Lagrangian schemes for hyperbolic advection‐reaction problems and Crank‐Nicolson schemes for purely diffusive problems. We show that the hybrid scheme is identical to the two end‐member schemes in the limit of infinite and zero Peclet number and remains accurate over a wide range of Peclet numbers. This scheme does not have a CFL stability criterion allowing the choice of time step to be decoupled from the spatial resolution. We compare numerical results with an analytic solution and test both an operator split version of our method and a combined version that solves advection and diffusion simultaneously. We also compare results of simple explicit and implicit numerical schemes and show that the semi‐Lagrangian Crank‐Nicolson (SLCN) scheme is both faster and more accurate on the same problem. We then apply the combined SLCN scheme to a more geologically relevant benchmark for calculating the thermal structure of a subduction zone. This problem demonstrates that the SLCN scheme can remain stable and accurate at large Courant numbers even in flows with highly curved streamlines. Finally, we introduce a variable order interpolation scheme for the semi‐Lagrangian schemes that reduces interpolation artifacts for sharp fronts without introducing numerical diffusion.

The content you want is available to Zendy users.

Already have an account? Click here to sign in.
Having issues? You can contact us here