IDEAS home Printed from https://ideas.repec.org/a/gam/jmathe/v9y2021i18p2316-d638937.html
   My bibliography  Save this article

A Massively Parallel Hybrid Finite Volume/Finite Element Scheme for Computational Fluid Dynamics

Author

Listed:
  • Laura Río-Martín

    (Laboratory of Applied Mathematics, DICAM, University of Trento, Via Mesiano 77, 38123 Trento, Italy)

  • Saray Busto

    (Laboratory of Applied Mathematics, DICAM, University of Trento, Via Mesiano 77, 38123 Trento, Italy
    Departamento de Matemática Aplicada a la Ingeniería Industrial, Universidad Politécnica de Madrid, José Gutierrez Abascal 2, 28006 Madrid, Spain)

  • Michael Dumbser

    (Laboratory of Applied Mathematics, DICAM, University of Trento, Via Mesiano 77, 38123 Trento, Italy)

Abstract

In this paper, we propose a novel family of semi-implicit hybrid finite volume/finite element schemes for computational fluid dynamics (CFD), in particular for the approximate solution of the incompressible and compressible Navier-Stokes equations, as well as for the shallow water equations on staggered unstructured meshes in two and three space dimensions. The key features of the method are the use of an edge-based/face-based staggered dual mesh for the discretization of the nonlinear convective terms at the aid of explicit high resolution Godunov-type finite volume schemes, while pressure terms are discretized implicitly using classical continuous Lagrange finite elements on the primal simplex mesh. The resulting pressure system is symmetric positive definite and can thus be very efficiently solved at the aid of classical Krylov subspace methods, such as a matrix-free conjugate gradient method. For the compressible Navier-Stokes equations, the schemes are by construction asymptotic preserving in the low Mach number limit of the equations, hence a consistent hybrid FV/FE method for the incompressible equations is retrieved. All parts of the algorithm can be efficiently parallelized, i.e., the explicit finite volume step as well as the matrix-vector product in the implicit pressure solver. Concerning parallel implementation, we employ the Message-Passing Interface (MPI) standard in combination with spatial domain decomposition based on the free software package METIS. To show the versatility of the proposed schemes, we present a wide range of applications, starting from environmental and geophysical flows, such as dambreak problems and natural convection, over direct numerical simulations of turbulent incompressible flows to high Mach number compressible flows with shock waves. An excellent agreement with exact analytical, numerical or experimental reference solutions is achieved in all cases. Most of the simulations are run with millions of degrees of freedom on thousands of CPU cores. We show strong scaling results for the hybrid FV/FE scheme applied to the 3D incompressible Navier-Stokes equations, using millions of degrees of freedom and up to 4096 CPU cores. The largest simulation shown in this paper is the well-known 3D Taylor-Green vortex benchmark run on 671 million tetrahedral elements on 32,768 CPU cores, showing clearly the suitability of the presented algorithm for the solution of large CFD problems on modern massively parallel distributed memory supercomputers.

Suggested Citation

  • Laura Río-Martín & Saray Busto & Michael Dumbser, 2021. "A Massively Parallel Hybrid Finite Volume/Finite Element Scheme for Computational Fluid Dynamics," Mathematics, MDPI, vol. 9(18), pages 1-41, September.
  • Handle: RePEc:gam:jmathe:v:9:y:2021:i:18:p:2316-:d:638937
    as

    Download full text from publisher

    File URL: https://www.mdpi.com/2227-7390/9/18/2316/pdf
    Download Restriction: no

    File URL: https://www.mdpi.com/2227-7390/9/18/2316/
    Download Restriction: no
    ---><---

    References listed on IDEAS

    as
    1. Busto, S. & Río-Martín, L. & Vázquez-Cendón, M.E. & Dumbser, M., 2021. "A semi-implicit hybrid finite volume/finite element scheme for all Mach number flows on staggered unstructured meshes," Applied Mathematics and Computation, Elsevier, vol. 402(C).
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Busto, S. & Dumbser, M. & Río-Martín, L., 2023. "An Arbitrary-Lagrangian-Eulerian hybrid finite volume/finite element method on moving unstructured meshes for the Navier-Stokes equations," Applied Mathematics and Computation, Elsevier, vol. 437(C).
    2. Boscheri, Walter & Tavelli, Maurizio, 2022. "High order semi-implicit schemes for viscous compressible flows in 3D," Applied Mathematics and Computation, Elsevier, vol. 434(C).
    3. Han Veiga, Maria & Micalizzi, Lorenzo & Torlo, Davide, 2024. "On improving the efficiency of ADER methods," Applied Mathematics and Computation, Elsevier, vol. 466(C).

    Most related items

    These are the items that most often cite the same works as this one and are cited by the same works as this one.
    1. Boscheri, Walter & Tavelli, Maurizio, 2022. "High order semi-implicit schemes for viscous compressible flows in 3D," Applied Mathematics and Computation, Elsevier, vol. 434(C).
    2. Michel-Dansac, Victor & Thomann, Andrea, 2022. "TVD-MOOD schemes based on implicit-explicit time integration," Applied Mathematics and Computation, Elsevier, vol. 433(C).
    3. Saray Busto & Michael Dumbser & Laura Río-Martín, 2021. "Staggered Semi-Implicit Hybrid Finite Volume/Finite Element Schemes for Turbulent and Non-Newtonian Flows," Mathematics, MDPI, vol. 9(22), pages 1-38, November.
    4. Abgrall, Rémi & Busto, Saray & Dumbser, Michael, 2023. "A simple and general framework for the construction of thermodynamically compatible schemes for computational fluid and solid mechanics," Applied Mathematics and Computation, Elsevier, vol. 440(C).
    5. Busto, S. & Dumbser, M. & Río-Martín, L., 2023. "An Arbitrary-Lagrangian-Eulerian hybrid finite volume/finite element method on moving unstructured meshes for the Navier-Stokes equations," Applied Mathematics and Computation, Elsevier, vol. 437(C).
    6. Frolkovič, Peter & Žeravý, Michal, 2023. "High resolution compact implicit numerical scheme for conservation laws," Applied Mathematics and Computation, Elsevier, vol. 442(C).

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:gam:jmathe:v:9:y:2021:i:18:p:2316-:d:638937. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a bibliographic reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: MDPI Indexing Manager (email available below). General contact details of provider: https://www.mdpi.com .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.