21773 articles – 15587 Notices  [english version]
HAL : hal-00610690, version 1

Fiche détaillée  Récupérer au format
International Symposium FVCA6, Prague : Tchèque, République (2011)
A Spectacular Vector Penalty-Projection Method for Darcy and Navier-Stokes Problems
Philippe Angot ( ) 1, Jean-Paul Caltagirone 2, Pierre Fabrie 3
(06/2011)

We present a new fast vector penalty-projection method (VPP$_{\eps}$), issued from noticeable improvements of previous works [7, 3, 4], to efficiently compute the solution of unsteady Navier-Stokes/Brinkman problems governing incompressible multiphase viscous flows. The method is also efficient to solve anisotropic Darcy problems. The key idea of the method is to compute at each time step an accurate and curl-free approximation of the pressure gradient increment in time. This method performs a two-step approximate divergence-free vector projection yielding a velocity divergence vanishing as $\mathcal{O}(\eps\, \dt)$, $\dt$ being the time step, with a penalty parameter $\eps$ as small as desired until the machine precision, e.g. $\eps= 10^{−14}$, whereas the solution algorithm can be extremely fast and cheap. The method is numerically validated on a benchmark problem for two-phase bubble dynamics where we compare it to the Uzawa augmented Lagrangian (UAL) and scalar incremental projection (SIP) methods. Moreover, a new test case for fluid-structure interaction problems is also investigated. That results in a robust method running faster than usual methods and being able to efficiently compute accurate solutions to sharp test cases whatever the density, viscosity or anisotropic permeability jumps, whereas other methods crash.
1 :  Laboratoire d'Analyse, Topologie, Probabilités (LATP)
CNRS : UMR6632 – Université de Provence - Aix-Marseille I – Université Paul Cézanne - Aix-Marseille III
2 :  Transferts, écoulements, fluides, énergétique (TREFLE)
CNRS : UMR8508 – Université Sciences et Technologies - Bordeaux I – École Nationale Supérieure de Chimie et de Physique de Bordeaux (ENSCPB) – Arts et Métiers ParisTech
3 :  Institut de Mathématiques de Bordeaux (IMB)
CNRS : UMR5251 – Université Sciences et Technologies - Bordeaux I – Université Victor Segalen - Bordeaux II
Mathématiques/Analyse numérique

Sciences de l'ingénieur/Mécanique/Mécanique des fluides

Physique/Mécanique/Mécanique des fluides
Vector penalty-projection – Penalty method – Splitting method – Multiphase Navier-Stokes/Brinkman – Anisotropic Darcy problem – Incompressible flows
Liste des fichiers attachés à ce document : 
PDF
FVCA6_ACF11.pdf(1.9 MB)