A Helmholtz--Hodge projection method using an iterative gauge computation to solve the 3D generalized Stokes problem

Abdou Garba Pierre Haldenwang 1
1 O&I - Ondes et Imagerie
LMA - Laboratoire de Mécanique et d'Acoustique [Marseille]
Abstract : The generalized Stokes problem (GSP) is broadly recognized as a keystone in implicit or semi-implicit discretizations of the Navier--Stokes equations (NSE), either incompressible or of low Mach number type, i.e., for which there exists a constraint on the vector field. The GSP is also known as the steady Darcy--Brinkman model for flows in porous media or in Hele--Shaw configurations. Up to now, only pressure- (preconditioned) Uzawa methods claimed to solve the three-dimensional (3D) GSP exactly (i.e., without introducing an error due to some time stepping). In the present article, we present another exact 3D solver that is developed in the particular context of the Helmholtz--Hodge projection in $(H^1)^d$. Instead of working iteratively with the pressure operator (i.e., the Schur complement), we are interested here in using an alternative scalar field, which performs the projection and is called the gauge. It turns out that iterative computation of the gauge presents efficient properties in terms of operator conditioning. To prove the relevance of the gauge method, we demonstrate several mathematical properties of the operator acting on the gauge, especially when preconditioned in an appropriate way. In particular, we prove that the gauge operator can be solved efficiently with any gradient method. The convergence properties of the continuous preconditioned operator are quantitatively estimated in the context of a semiperiodic two-dimensional (2D) domain. We then study the discrete gauge operator implemented in the framework of Chebyshev 2D-3D pseudospectral methods. It exhibits the same favorable properties as those found for the continuous operator. Finally, we perform numerical experiments that corroborate our analyses.
Type de document :
Article dans une revue
SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2013, 35 (3), pp.A1560-A1583. <10.1137/110860902>
Liste complète des métadonnées

https://hal.archives-ouvertes.fr/hal-00905812
Contributeur : Elena Rosu <>
Soumis le : lundi 18 novembre 2013 - 16:50:56
Dernière modification le : mardi 10 mai 2016 - 11:45:42

Identifiants

Citation

Abdou Garba, Pierre Haldenwang. A Helmholtz--Hodge projection method using an iterative gauge computation to solve the 3D generalized Stokes problem. SIAM Journal on Scientific Computing, Society for Industrial and Applied Mathematics, 2013, 35 (3), pp.A1560-A1583. <10.1137/110860902>. <hal-00905812>

Partager

Métriques

Consultations de la notice

181