Roland Masson

V. Dolean, M. J. Gander, F. Kwok et al.

For linear problems, domain decomposition methods can be used directly as iterative solvers, but also as preconditioners for Krylov methods. In practice, Krylov acceleration is almost always used, since the Krylov method finds a much better residual polynomial than the stationary iteration, and thus converges much faster. We show in this paper that also for non-linear problems, domain decomposition methods can either be used directly as iterative solvers, or one can use them as preconditioners for Newton's method. For the concrete case of the parallel Schwarz method, we show that we obtain a preconditioner we call RASPEN (Restricted Additive Schwarz Preconditioned Exact Newton) which is similar to ASPIN (Additive Schwarz Preconditioned Inexact Newton), but with all components directly defined by the iterative method. This has the advantage that RASPEN already converges when used as an iterative solver, in contrast to ASPIN, and we thus get a substantially better preconditioner for Newton's method. The iterative construction also allows us to naturally define a coarse correction using the multigrid full approximation scheme, which leads to a convergent two level non-linear iterative domain decomposition method and a two level RASPEN non-linear preconditioner. We illustrate our findings with numerical results on the Forchheimer equation and a non-linear diffusion problem.

Jerome Droniou, Julian Hennicker, Roland Masson

We present a new model for two phase Darcy flows in fractured media, in which fractures are modelled as submanifolds of codimension one with respect to the surrounding domain (matrix). Fractures can act as drains or as barriers, since pressure discontinuities at the matrix-fracture interfaces are permitted. Additionally, a layer of damaged rock at the matrix-fracture interfaces is accounted for. The numerical analysis is carried out in the general framework of the Gradient Discretisation Method. Compactness techniques are used to establish convergence results for a wide range of possible numerical schemes; the existence of a solution for the two phase flow model is obtained as a byproduct of the convergence analysis. A series of numerical experiments conclude the paper, with a study of the influence of the damaged layer on the numerical solution.

K. Brenner, J. Hennicker, R. Masson et al.

We investigate the discretization of Darcy flow through fractured porous media on general meshes. We consider a hybrid dimensional model, invoking a complex network of planar fractures. The model accounts for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. The numerical analysis is performed in the general framework of gradient discretizations which is extended to the model under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to satisfy the gradient scheme framework, which yields, in particular, convergence. Numerical tests confirm the theoretical results. Gradient Discretization; Darcy Flow, Discrete Fracture Networks, Finite Volume