A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems
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 th...
| Main Authors: | , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
2006
|
| _version_ | 1826292889200623616 |
|---|---|
| author | Spiegelman, M Katz, R |
| author_facet | Spiegelman, M Katz, R |
| author_sort | Spiegelman, M |
| collection | OXFORD |
| description | 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.© 2006 by the American Geophysical Union. |
| first_indexed | 2024-03-07T03:21:39Z |
| format | Journal article |
| id | oxford-uuid:b7a3ec35-e9f4-46f8-8246-33a008f97e4b |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-07T03:21:39Z |
| publishDate | 2006 |
| record_format | dspace |
| spelling | oxford-uuid:b7a3ec35-e9f4-46f8-8246-33a008f97e4b2022-03-27T04:50:04ZA semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problemsJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:b7a3ec35-e9f4-46f8-8246-33a008f97e4bEnglishSymplectic Elements at Oxford2006Spiegelman, MKatz, RWe 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.© 2006 by the American Geophysical Union. |
| spellingShingle | Spiegelman, M Katz, R A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title | A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title_full | A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title_fullStr | A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title_full_unstemmed | A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title_short | A semi-Lagrangian Crank-Nicolson algorithm for the numerical solution of advection-diffusion problems |
| title_sort | semi lagrangian crank nicolson algorithm for the numerical solution of advection diffusion problems |
| work_keys_str_mv | AT spiegelmanm asemilagrangiancranknicolsonalgorithmforthenumericalsolutionofadvectiondiffusionproblems AT katzr asemilagrangiancranknicolsonalgorithmforthenumericalsolutionofadvectiondiffusionproblems AT spiegelmanm semilagrangiancranknicolsonalgorithmforthenumericalsolutionofadvectiondiffusionproblems AT katzr semilagrangiancranknicolsonalgorithmforthenumericalsolutionofadvectiondiffusionproblems |