103 research outputs found
A robust adaptive algebraic multigrid linear solver for structural mechanics
The numerical simulation of structural mechanics applications via finite
elements usually requires the solution of large-size and ill-conditioned linear
systems, especially when accurate results are sought for derived variables
interpolated with lower order functions, like stress or deformation fields.
Such task represents the most time-consuming kernel in commercial simulators;
thus, it is of significant interest the development of robust and efficient
linear solvers for such applications. In this context, direct solvers, which
are based on LU factorization techniques, are often used due to their
robustness and easy setup; however, they can reach only superlinear complexity,
in the best case, thus, have limited applicability depending on the problem
size. On the other hand, iterative solvers based on algebraic multigrid (AMG)
preconditioners can reach up to linear complexity for sufficiently regular
problems but do not always converge and require more knowledge from the user
for an efficient setup. In this work, we present an adaptive AMG method
specifically designed to improve its usability and efficiency in the solution
of structural problems. We show numerical results for several practical
applications with millions of unknowns and compare our method with two
state-of-the-art linear solvers proving its efficiency and robustness.Comment: 50 pages, 16 figures, submitted to CMAM
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
Hybrid multigrid methods for high-order discontinuous Galerkin discretizations
The present work develops hybrid multigrid methods for high-order
discontinuous Galerkin discretizations of elliptic problems. Fast matrix-free
operator evaluation on tensor product elements is used to devise a
computationally efficient PDE solver. The multigrid hierarchy exploits all
possibilities of geometric, polynomial, and algebraic coarsening, targeting
engineering applications on complex geometries. Additionally, a transfer from
discontinuous to continuous function spaces is performed within the multigrid
hierarchy. This does not only further reduce the problem size of the
coarse-grid problem, but also leads to a discretization most suitable for
state-of-the-art algebraic multigrid methods applied as coarse-grid solver. The
relevant design choices regarding the selection of optimal multigrid coarsening
strategies among the various possibilities are discussed with the metric of
computational costs as the driving force for algorithmic selections. We find
that a transfer to a continuous function space at highest polynomial degree (or
on the finest mesh), followed by polynomial and geometric coarsening, shows the
best overall performance. The success of this particular multigrid strategy is
due to a significant reduction in iteration counts as compared to a transfer
from discontinuous to continuous function spaces at lowest polynomial degree
(or on the coarsest mesh). The coarsening strategy with transfer to a
continuous function space on the finest level leads to a multigrid algorithm
that is robust with respect to the penalty parameter of the SIPG method.
Detailed numerical investigations are conducted for a series of examples
ranging from academic test cases to more complex, practically relevant
geometries. Performance comparisons to state-of-the-art methods from the
literature demonstrate the versatility and computational efficiency of the
proposed multigrid algorithms
- …