Méthode de décomposition de domaine et conditions aux limites artificielles en mécanique des fluides: méthode Optimisée d'Orde 2.

Ce travail a pour objet le developpement et l'etude d'une methode de decomposition de domaine, la methode Optimisee d'Ordre 2 (OO2), pour la resolution de l'equation de convection-diffusion. Son atout principal est de permettre d'utiliser un decoupage quelconque du domaine, sans savoir a l'avance ou sont situes les phenomenes physiques tels que les couches limites ou les zones de recirculation. La methode OO2 est une methode de decomposition de domaine sans recouvrement, iterative, parallelisable. Le domaine de calcul est divise en sous-domaines, et on resout le probleme de depart dans chaque sous-domaine, avec des conditions de raccord specifiques sur les interfaces des sous-domaines. Ce sont des conditions differentielles d'ordre 1 dans la direction normale et d'ordre 2 dans la direction tangente a l'interface qui approchent, par une procedure d'optimisation, les Conditions aux Limites Artificielles (CLA). L'utilisation des CLA en decomposition de domaine permet de definir des algorithmes stables. Une reformulation de la methode de Schwarz conduit a un probleme d'interface. Celui-ci est resolu par une methode iterative de type Krylov (BICG-STAB, GMRES, GCR). La methode est appliquee a un schema aux differences finies decentre, puis a un schema volumes finis. Un preconditionneur ``basses frequences'' est ensuite introduit et etudie, dans le but d'avoir une convergence independante du nombre de sous-domaines. Ce preconditionneur est une extension aux problemes non-symetriques d'un preconditionneur utilise pour des problemes symetriques. Enfin, l'utilisation de conditions differentielles d'ordre 2 le long de l'interface necessite d'ajouter des conditions de raccord aux points de croisement des sous-domaines. Une etude est menee a ce sujet, qui permet de montrer que les problemes dans chaque sous-domaine sont bien poses.

[1]  A Schwarz auditive method with high order interface conditions and nonoverlapping subdomains , 1998 .

[2]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[3]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[4]  Christopher C. Borel,et al.  High performance parallelized implicit Euler solver for the analysis of unsteady aerodynamic flows , 1992 .

[5]  Laurence Halpern,et al.  Absorbing boundary conditions for the discretization schemes of the one-dimensional wave equation , 1982 .

[6]  C. Farhat,et al.  The two-level FETI method for static and dynamic plate problems Part I: An optimal iterative solver for biharmonic systems , 1998 .

[7]  G. Golub,et al.  Iterative solution of linear systems , 1991, Acta Numerica.

[8]  P. Tallec Domain decomposition methods in computational mechanics , 1994 .

[9]  C. Farhat,et al.  The two-level FETI method. Part II: Extension to shell problems, parallel implementation and performance results , 1998 .

[10]  Alfio Quarteroni,et al.  Adaptive Domain Decomposition Methods for Advection-Diffusion Problems , 1995 .

[11]  P. Grisvard Singularities in Boundary Value Problems , 1992 .

[12]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[13]  Charbel Farhat,et al.  Parallel Implementation of the Two-level Feti Method , .

[14]  Thomas Hagstrom,et al.  Numerical Experiments on a Domain Decomposition Algorithm for Nonlinear Elliptic Boundary Value Problems , 1988 .

[15]  Patrick Joly,et al.  Mathematical and Numerical Aspects of Wave Propagation Phenomena , 1991 .

[16]  V. Venkatakrishnan,et al.  Viscous computations using a direct solver , 1990 .

[17]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[18]  F. Nataf On the Use of Open Boundary Conditions in Block Gauss-Seidel Methods for the Convection-Di usion Equation , 1993 .

[19]  C. Farhat,et al.  Optimal convergence properties of the FETI domain decomposition method , 1994 .