Caroline Japhet

Laurence Halpern, Jérémie Szeftel, Caroline Japhet

We design and analyze a Schwarz waveform relaxation algorithm for domain decomposition of advection-diffusion-reaction problems with strong heterogeneities. The interfaces are curved, and we use optimized Robin or Ventcell transmission conditions. We analyze the semi-discretization in time with Discontinuous Galerkin as well. We also show two-dimensional numerical results using generalized mortar finite elements in space.

Thi Thao Phuong Hoang, Jérôme Jaffré, Caroline Japhet et al.

This paper is concerned with global-in-time, nonoverlapping domain decomposition methods for the mixed formulation of the diffusion problem. Two approaches are considered: one uses the time-dependent Steklov-Poincaré operator and the other uses Optimized Schwarz Waveform Relaxation (OSWR) based on Robin transmission conditions. For each method, a mixed formulation of an interface problem on the space-time interfaces between subdomains is derived, and different time grids are employed to adapt to different time scales in the subdomains. Demonstrations of the well-posedness of the subdomain problems involved in each method and a convergence proof of the OSWR algorithm are given for the mixed formulation. Numerical results for 2D problems with strong heterogeneities are presented to illustrate the performance of the two methods.

Thi-Thao-Phuong Hoang, Caroline Japhet, Michel Kern et al.

In this paper we are interested in the "fast path" fracture and we aim to use global-in-time, nonoverlapping domain decomposition methods to model flow and transport problems in a porous medium containing such a fracture. We consider a reduced model in which the fracture is treated as an interface between the two subdomains. Two domain decomposition methods are considered: one uses the time-dependent SteklovPoincar{é} operator and the other uses optimized Schwarz waveform relaxation (OSWR) based on Ventcell transmission conditions. For each method, a mixed formulation of an interface problem on the space-time interface is derived, and different time grids are employed to adapt to different time scales in the subdomains and in the fracture. Demonstrations of the well-posedness of the Ventcell subdomain problems is given for the mixed formulation. An analysis for the convergence factor of the OSWR algorithm is given in the case with fractures to compute the optimized parameters. Numerical results for two-dimensional problems with strong heterogeneities are presented to illustrate the performance of the two methods.

Thi-Thao-Phuong Hoang, Caroline Japhet, Michel Kern et al.

This paper is concerned with the numerical solution of porous-media flow and transport problems , i. e. heterogeneous, advection-diffusion problems. Its aim is to investigate numerical schemes for these problems in which different time steps can be used in different parts of the domain. Global-in-time, non-overlapping domain-decomposition methods are coupled with operator splitting making possible the different treatment of the advection and diffusion terms. Two domain-decomposition methods are considered: one uses the time-dependent Steklov--Poincar{é} operator and the other uses optimized Schwarz waveform relaxation (OSWR) based on Robin transmission conditions. For each method, a mixed formulation of an interface problem on the space-time interface is derived, and different time grids are employed to adapt to different time scales in the subdomains. A generalized Neumann-Neumann preconditioner is proposed for the first method. To illustrate the two methods numerical results for two-dimensional problems with strong heterogeneities are presented. These include both academic problems and more realistic prototypes for simulations for the underground storage of nuclear waste.

Caroline Japhet, Yvon Maday, Frédéric Nataf

We design and analyze a new non-conforming domain decomposition method based on Schwarz type approaches that allows for the use of Robin interface conditions on non-conforming grids. The method is proven to be well posed, and the iterative solver to converge. The error analysis is performed in 2D piecewise polynomials of low and high order and extended in 3D for $P_1$ elements. Numerical results in 2D illustrate the new method.