Skip to main content
Log in

Efficient implementation of characteristic-based schemes on unstructured triangular grids

  • Published:
Computational and Applied Mathematics Aims and scope Submit manuscript

Abstract

Using characteristics to treat advection terms in time-dependent PDEs leads to a class of schemes, e.g., semi-Lagrangian and Lagrange–Galerkin schemes, which preserve stability under large Courant numbers, and may therefore be appealing in many practical situations. Unfortunately, the need of locating the feet of characteristics may cause a serious drop of efficiency in the case of unstructured space grids, and thus prevent the use of large time-step schemes on complex geometries. In this paper, we perform an in-depth analysis of the main recipes available for characteristic location, and propose a technique to improve the efficiency of this phase, using additional information related to the advecting vector field. This results in a clear improvement of execution times in the unstructured case, thus extending the range of applicability of large time-step schemes.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12
Fig. 13
Fig. 14
Fig. 15
Fig. 16
Fig. 17

Similar content being viewed by others

References

  • Bonaventura L, Calzola E, Carlini E, Ferretti R (2020) A fully semi-Lagrangian method for the Navier–Stokes equations in primitive variables. In: van Brummelen H, Corsini A, Perotto S, Rozza G (eds) Numerical Methods for Flows. Lecture Notes in Computational Science and Engineering, vol 132. Springer, Cham

  • Bonaventura L, Calzola E, Carlini E, Ferretti R (2021) Second order fully semi-Lagrangian discretizations of advection–diffusion–reaction systems. J Sci Comp 88:1–29

    Article  MathSciNet  Google Scholar 

  • Boscheri W (2020) A space-time semi-Lagrangian advection scheme on staggered Voronoi meshes applied to free surface flows. Comput Fluids 202

  • Boscheri W, Dumbser M, Righetti M (2013) A semi-implicit scheme for 3D free surface flows with high-order velocity reconstruction on unstructured Voronoi meshes. Int J Num Meth Fluids 72:607–631

    Article  MathSciNet  Google Scholar 

  • Devillers O, Pion S, Teillaud M (2001) Walking in a triangulation. In: Proceedings of the seventeenth annual symposium on Computational geometry, 106–114

  • Douglas J, Russell TF (1982) Numerical methods for convection-dominated diffusion problems based on combining the method of characteristics with finite element or finite difference procedures. SIAM J Num Anal 19:871–885

    Article  MathSciNet  Google Scholar 

  • Falcone M, Ferretti R (2013) Semi-Lagrangian approximation schemes for linear and Hamilton-Jacobi equations. SIAM, Philadelphia

    Book  Google Scholar 

  • Ferretti R, Mehrenberger M (2020) Stability of Semi-Lagrangian schemes of arbitrary odd degree under constant and variable advection speed. Math Comp 89:1783–1805

    Article  MathSciNet  Google Scholar 

  • Finkel RA, Bentley JL (1974) Quad trees a data structure for retrieval on composite keys. Acta Informatica 4:1–9

    Article  Google Scholar 

  • Giraldo FX (1998) The Lagrange–Galerkin spectral element method on unstructured quadrilateral grids. J Comp Phys 147:114–146

    Article  MathSciNet  Google Scholar 

  • Giraldo FX (2000) The Lagrange–Galerkin method for the two-dimensional shallow water equations on adaptive grids. Int J Num Meth Fluids 33:789–832

    Article  MathSciNet  Google Scholar 

  • Pironneau O (1982) On the transport-diffusion algorithm and its application to the Navier–Stokes equations. Num Math 38:309–332

    Article  MathSciNet  Google Scholar 

  • Restelli M, Bonaventura L, Sacco R (2006) A semi-Lagrangian discontinuous Galerkin method for scalar advection by incompressible flows. J Comp Phys 216:195–215

    Article  MathSciNet  Google Scholar 

  • Russell TF, Celia MA (2002) An overview of research on Eulerian–Lagrangian localized adjoint methods (ELLAM). Adv Water Res 25:1215–1231

    Article  Google Scholar 

  • Triangle, A two-dimensional quality mesh generator and delaunay triangulator. https://www.cs.cmu.edu/~quake/triangle

Download references

Acknowledgements

This work has been partially supported by the PRIN 2017 project “Innovative Numerical Methods for Evolutionary Partial Differential Equations and Applications”, by the INdAM–GNCS project “Approssimazione numerica di problemi di natura iperbolica ed applicazioni” and by Roma Tre University. We thank Dr. Beatrice Beco and Dr. Lorenzo Della Cioppa for taking part in the first steps of this work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to R. Ferretti.

Additional information

Communicated by Raphaèle Herbin.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cacace, S., Ferretti, R. Efficient implementation of characteristic-based schemes on unstructured triangular grids. Comp. Appl. Math. 41, 19 (2022). https://doi.org/10.1007/s40314-021-01716-y

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s40314-021-01716-y

Keywords

Mathematics Subject Classification

Navigation

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy