1,124 research outputs found
Space-time domain decomposition for advection-diffusion problems in mixed formulations
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{\'e} 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
A Domain Decomposition method based on iterative Operator Splitting method.
In this article a new approach is proposed for constructing domain decomposition methods based on iterative operator splitting methods. We study the convergence properties of such a method. The main feature of the method is the decoupling the space and time dimension. We confirm with two numerical applications the effectiveness of the proposed iterative operator splitting method in comparison with classical Schwarz waveform relaxation method as a standard method for domain decomposition. We provide improved results and convergence rates. The efficiency of considering the whole domain in the case of the iterative operator splitting method allows more accurate results
Optimized Schwarz waveform relaxation and discontinuous Galerkin time stepping for heterogeneous problems
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
Optimized Schwarz waveform relaxation for Primitive Equations of the ocean
In this article we are interested in the derivation of efficient domain
decomposition methods for the viscous primitive equations of the ocean. We
consider the rotating 3d incompressible hydrostatic Navier-Stokes equations
with free surface. Performing an asymptotic analysis of the system with respect
to the Rossby number, we compute an approximated Dirichlet to Neumann operator
and build an optimized Schwarz waveform relaxation algorithm. We establish the
well-posedness of this algorithm and present some numerical results to
illustrate the method
Nonlinear Preconditioning: How to use a Nonlinear Schwarz Method to Precondition Newton's Method
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
- …