22,145 research outputs found
Time-parallel iterative solvers for parabolic evolution equations
We present original time-parallel algorithms for the solution of the implicit
Euler discretization of general linear parabolic evolution equations with
time-dependent self-adjoint spatial operators. Motivated by the inf-sup theory
of parabolic problems, we show that the standard nonsymmetric time-global
system can be equivalently reformulated as an original symmetric saddle-point
system that remains inf-sup stable with respect to the same natural parabolic
norms. We then propose and analyse an efficient and readily implementable
parallel-in-time preconditioner to be used with an inexact Uzawa method. The
proposed preconditioner is non-intrusive and easy to implement in practice, and
also features the key theoretical advantages of robust spectral bounds, leading
to convergence rates that are independent of the number of time-steps, final
time, or spatial mesh sizes, and also a theoretical parallel complexity that
grows only logarithmically with respect to the number of time-steps. Numerical
experiments with large-scale parallel computations show the effectiveness of
the method, along with its good weak and strong scaling properties
A Class of Stable, Globally Noniterative, Nonoverlapping Domain Decomposition Algorithms for the Simulation of Parabolic Evolutionary Systems.
Parabolic systems are governed by time dependent partial differential equations. To obtain a high simulation quality that captures important features of a parabolic system requires solving the governing equation to an adequate accuracy, which necessitates a large sampling size in the spatial and temporal dimensions, and hence a large amount of simulation data and high computing cost. Domain decomposition is an effective method of divide-and-conquer paradigm that divides the problem domain into several subdomains, reducing the original problem into several smaller interdependent problems which can be solved in parallel. In this dissertation, we propose a class of stabilized explicit-implicit time marching (SEITM) domain decomposition algorithms for parabolic equations. Explicit-implicit time marching (EITM) algorithms are globally non-iterative nonoverlapping domain decomposition methods, which, when compared with Schwartz algorithm based parabolic solvers, are both computationally and communicationally efficient for each time step simulation but suffer from small time step size restrictions due to conditional stability. The proposed stabilization techniques in the SEITM algorithms retain the time-stepwise efficiency in computation and communication of the EITM algorithms but free the algorithms from small time step size restrictions, rendering SEITM algorithms excellent candidates for large scale parallel simulation problems. Three algorithms of the SEITM class are presented in this dissertation, which are mathematically analyzed and experimentally tested to show excellent numerical stability, computation and communication efficiencies, and high parallel speedup and scalability
Hybrid PDE solver for data-driven problems and modern branching
The numerical solution of large-scale PDEs, such as those occurring in
data-driven applications, unavoidably require powerful parallel computers and
tailored parallel algorithms to make the best possible use of them. In fact,
considerations about the parallelization and scalability of realistic problems
are often critical enough to warrant acknowledgement in the modelling phase.
The purpose of this paper is to spread awareness of the Probabilistic Domain
Decomposition (PDD) method, a fresh approach to the parallelization of PDEs
with excellent scalability properties. The idea exploits the stochastic
representation of the PDE and its approximation via Monte Carlo in combination
with deterministic high-performance PDE solvers. We describe the ingredients of
PDD and its applicability in the scope of data science. In particular, we
highlight recent advances in stochastic representations for nonlinear PDEs
using branching diffusions, which have significantly broadened the scope of
PDD.
We envision this work as a dictionary giving large-scale PDE practitioners
references on the very latest algorithms and techniques of a non-standard, yet
highly parallelizable, methodology at the interface of deterministic and
probabilistic numerical methods. We close this work with an invitation to the
fully nonlinear case and open research questions.Comment: 23 pages, 7 figures; Final SMUR version; To appear in the European
Journal of Applied Mathematics (EJAM
Algorithms and Data Structures for Multi-Adaptive Time-Stepping
Multi-adaptive Galerkin methods are extensions of the standard continuous and
discontinuous Galerkin methods for the numerical solution of initial value
problems for ordinary or partial differential equations. In particular, the
multi-adaptive methods allow individual and adaptive time steps to be used for
different components or in different regions of space. We present algorithms
for efficient multi-adaptive time-stepping, including the recursive
construction of time slabs and adaptive time step selection. We also present
data structures for efficient storage and interpolation of the multi-adaptive
solution. The efficiency of the proposed algorithms and data structures is
demonstrated for a series of benchmark problems.Comment: ACM Transactions on Mathematical Software 35(3), 24 pages (2008
Multigrid waveform relaxation for the time-fractional heat equation
In this work, we propose an efficient and robust multigrid method for solving
the time-fractional heat equation. Due to the nonlocal property of fractional
differential operators, numerical methods usually generate systems of equations
for which the coefficient matrix is dense. Therefore, the design of efficient
solvers for the numerical simulation of these problems is a difficult task. We
develop a parallel-in-time multigrid algorithm based on the waveform relaxation
approach, whose application to time-fractional problems seems very natural due
to the fact that the fractional derivative at each spatial point depends on the
values of the function at this point at all earlier times. Exploiting the
Toeplitz-like structure of the coefficient matrix, the proposed multigrid
waveform relaxation method has a computational cost of
operations, where is the number of time steps and is the number of
spatial grid points. A semi-algebraic mode analysis is also developed to
theoretically confirm the good results obtained. Several numerical experiments,
including examples with non-smooth solutions and a nonlinear problem with
applications in porous media, are presented
- …