170 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
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
Optimized Schwarz Waveform Relaxation for Advection Reaction Diffusion Equations in Two Dimensions
Optimized Schwarz Waveform Relaxation methods have been developed over the
last decade for the parallel solution of evolution problems. They are based on
a decomposition in space and an iteration, where only subproblems in space-time
need to be solved. Each subproblem can be simulated using an adapted numerical
method, for example with local time stepping, or one can even use a different
model in different subdomains, which makes these methods very suitable also
from a modeling point of view. For rapid convergence however, it is important
to use effective transmission conditions between the space-time subdomains, and
for best performance, these transmission conditions need to take the physics of
the underlying evolution problem into account. The optimization of these
transmission conditions leads to a mathematically hard best approximation
problem of homographic type. We study in this paper in detail this problem for
the case of linear advection reaction diffusion equations in two spatial
dimensions. We prove comprehensively best approximation results for
transmission conditions of Robin and Ventcel type. We give for each case closed
form asymptotic values for the parameters, which guarantee asymptotically best
performance of the iterative methods. We finally show extensive numerical
experiments, and we measure performance corresponding to our analysisComment: 42 page
A new Algorithm Based on Factorization for Heterogeneous Domain Decomposition
Often computational models are too expensive to be solved in the entire
domain of simulation, and a cheaper model would suffice away from the main zone
of interest. We present for the concrete example of an evolution problem of
advection reaction diffusion type a heterogeneous domain decomposition
algorithm which allows us to recover a solution that is very close to the
solution of the fully viscous problem, but solves only an inviscid problem in
parts of the domain. Our new algorithm is based on the factorization of the
underlying differential operator, and we therefore call it factorization
algorithm. We give a detailed error analysis, and show that we can obtain
approximations in the viscous region which are much closer to the viscous
solution in the entire domain of simulation than approximations obtained by
other heterogeneous domain decomposition algorithms from the literature.Comment: 23 page
Schwarz waveform relaxation with adaptive pipelining
Schwarz waveform relaxation (SWR) methods have been developed to solve a wide range of diffusion-dominated and reaction-dominated equations. The appeal of these methods stems primarily from their ability to use nonconforming space-time discretizations; SWR methods are consequently well-adapted for coupling models with highly varying spatial and time scales. The efficacy of SWR methods is questionable, however, since in each iteration, one propagates an error across the entire time interval. In this manuscript, we introduce an adaptive pipeline approach wherein one subdivides the computational domain into space-time blocks, and adaptively selects the waveform iterates which should be updated given a fixed number of computational workers. Our method is complementary to existing space and time parallel methods, and can be used to obtain additional speedup when the saturation point is reached for other types of parallelism. We analyze these waveform relaxation with adaptive pipelining (WRAP) methods to show convergence and the theoretical speedup that can be expected. Numerical experiments on solutions to the linear heat equation, the advection-diffusion equation, and a reaction-diffusion equation illustrate features and efficacy of WRAP methods for various transmission conditions
- …