Space-time domain decomposition methods for mixed formulations of flow and transport problems in porous media
Séminaire Modèles et Algorithmes Déterministes: EDP-MOISE
13/02/2014 - 11:00 Mme Caroline Japhet (LAGA, Paris 13) Salle 1 - Tour IRMA
Flow and transport problems in porous media are well-known for their high computational cost. In the simulation of an underground nuclear waste disposal site, one has to work with extremely different length and time scales, and highly variable coefficients while satisfying strict accuracy requirements. One strategy for tackling these difficulties is to apply a non-overlapping domain decomposition method which allows local adaptation in both space and time and makes possible the use of parallel algorithms. We present two approaches. The first method uses a time-dependent Steklov-Poincaré operator and a generalized Neumann-Neumann preconditioner with weight matrices to handle the heterogeneities. The second method uses the optimized Schwarz waveform relaxation (OSWR). The OSWR algorithm uses more general (Robin or Ventcell) transmission operators in which coefficients can be optimized to improve convergence rates. For both approaches, an interface problem is derived on the space-time interface. Thus, different time steps can be used in different subdomains adapted to their physical properties. We use operator splitting to treat differently the advection and the diffunion. We show numerical experiments for various test cases, both academic and more realistic prototypes for nuclear waste disposal simulation. We finally study extensions of the two methods to the case in which the interface represents a discrete-fracture in a reduced fracture model for flow in a fractured porous medium.