A hybrid domain decomposition and multigrid method for the acceleration of compressible viscous flow calculations on unstructured triangular meshes

Victorita Dolean, Stéphane Lanteri

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.

LanguageEnglish
Pages287-304
Number of pages18
JournalInternational Journal of Computational Fluid Dynamics
Volume14
Issue number4
DOIs
Publication statusPublished - 1 Jan 2001

Fingerprint

multigrid methods
viscous flow
Viscous flow
mesh
Decomposition
decomposition
Linear systems
linear systems
Domain decomposition methods
Steady flow
Airfoils
formulations
Navier Stokes equations
airfoils
steady flow
Agglomeration
matrices
agglomeration
Navier-Stokes equation
Fluxes

Keywords

  • domain decomposition method
  • finite elements
  • finite volumes
  • multigrid algorithm
  • Navier-Stokes equations
  • parallel computing
  • triangular meshes

Cite this

@article{3434964ce13d4769bf6d90e98c02fa18,
title = "A hybrid domain decomposition and multigrid method for the acceleration of compressible viscous flow calculations on unstructured triangular meshes",
abstract = "This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.",
keywords = "domain decomposition method, finite elements, finite volumes, multigrid algorithm, Navier-Stokes equations, parallel computing, triangular meshes",
author = "Victorita Dolean and St{\'e}phane Lanteri",
year = "2001",
month = "1",
day = "1",
doi = "10.1080/10618560108940730",
language = "English",
volume = "14",
pages = "287--304",
journal = "International Journal of Computational Fluid Dynamics",
issn = "1061-8562",
number = "4",

}

TY - JOUR

T1 - A hybrid domain decomposition and multigrid method for the acceleration of compressible viscous flow calculations on unstructured triangular meshes

AU - Dolean, Victorita

AU - Lanteri, Stéphane

PY - 2001/1/1

Y1 - 2001/1/1

N2 - This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.

AB - This paper is concerned with the formulation and the evaluation of a hybrid solution method that makes use of domain decomposition and multigrid principles for the calculation of two-dimensional compressible viscous flows on unstructured triangular meshes. More precisely, a non-overlapping additive domain decomposition method is used to coordinate concurrent subdomain solutions with a multigrid method. This hybrid method is developed in the context of a flow solver for the Navier-Stokes equations which is based on a combined finite element/finite volume formulation on unstructured triangular meshes. Time integration of the resulting semi-discrete equations is performed using a linearized backward Euler implicit scheme. As a result, each pseudo time step requires the solution of a sparse linear system. In this study, a non-overlapping domain decomposition algorithm is used for advancing the solution at each implicit time step. Algebraically, the Schwarz algorithm is equivalent to a Jacobi iteration on a linear system whose matrix has a block structure. A substructuring technique can be applied to this matrix in order to obtain a fully implicit scheme in terms of interface unknowns. In the present approach, the interface unknowns are numerical fluxes. The interface system is solved by means of a full GMRES method. Here, the local system solves that are induced by matrix-vector products with the interface operator, are performed using a multigrid by volume agglomeration method. The resulting hybrid domain decomposition and multigrid solver is applied to the computation of several steady flows around a geometry of NACA0012 airfoil.

KW - domain decomposition method

KW - finite elements

KW - finite volumes

KW - multigrid algorithm

KW - Navier-Stokes equations

KW - parallel computing

KW - triangular meshes

UR - http://www.scopus.com/inward/record.url?scp=0035634180&partnerID=8YFLogxK

U2 - 10.1080/10618560108940730

DO - 10.1080/10618560108940730

M3 - Article

VL - 14

SP - 287

EP - 304

JO - International Journal of Computational Fluid Dynamics

T2 - International Journal of Computational Fluid Dynamics

JF - International Journal of Computational Fluid Dynamics

SN - 1061-8562

IS - 4

ER -