This paper reviews the basic mathematical ideas and convergence analysis of domain decomposition methods. These are parallel and scalable iterative methods for the efficient numerical solution of partial differential equations. Two examples are then presented showing the application of domain decomposition methods to large-scale numerical simulations in computational mechanics and electrocardiology.
We present a parallel solution algorithm for the transient heat equation in one and two spatial dimensions. The problem is discretized in space by the lowest-order conforming finite element method. Further, a one-step time integration scheme is used for the numerical solution of the arising system of ordinary differential equations. For the latter, the parareal method decomposing the time interval into subintervals is employed. It leads to parallel solution of smaller time-dependent problems. At...
We consider mixed finite element discretizations of second order elliptic boundary value problems. Emphasis is on the efficient iterative solution by multilevel techniques with respect to an adaptively generated hierarchy of nonuniform triangulations. In particular, we present two multilevel solvers, the first one relying on ideas from domain decomposition and the second one resulting from mixed hybridization. Local refinement of the underlying triangulations is done by efficient and reliable a...
In this paper a black-box solver based on combining the unknowns aggregation with smoothing is suggested. Convergence is improved by overcorrection. Numerical experiments demonstrate the efficiency.
We consider linear elliptic systems which arise in coupled elastic continuum mechanical models. In these systems, the strain tensor εP := sym (P-1∇u) is redefined to include a matrix valued inhomogeneity P(x) which cannot be described by a space dependent fourth order elasticity tensor. Such systems arise naturally in geometrically exact plasticity or in problems with eigenstresses. The tensor field P induces a structural change of the elasticity equations. For such a model the FETI-DP method is...
We consider linear elliptic systems which arise in coupled elastic continuum mechanical models. In these systems, the strain tensor εP := sym (P-1∇u) is redefined to include a matrix valued inhomogeneity P(x) which cannot be described by a space dependent fourth order elasticity tensor. Such systems arise naturally in geometrically exact plasticity or in problems with eigenstresses. The tensor field P induces a structural change of the elasticity equations. For such a model the FETI-DP method is...
We propose and analyze a domain decomposition method on non-matching grids for partial differential equations with non-negative characteristic form. No weak or strong continuity of the finite element functions, their normal derivatives, or linear combinations of the two is imposed across the boundaries of the subdomains. Instead, we employ suitable bilinear forms defined on the common interfaces, typical of discontinuous Galerkin approximations. We prove an error bound which is optimal with respect...
We propose and analyze a domain decomposition method on non-matching grids for partial differential equations with non-negative characteristic form. No weak or strong continuity of the finite element functions, their normal derivatives, or linear combinations of the two is imposed across the boundaries of the subdomains. Instead, we employ suitable bilinear forms defined on the common interfaces, typical of discontinuous Galerkin approximations. We prove an error bound which is optimal with respect...
The smoothed aggregation method has became a widely used tool for solving the linear systems arising by the discretization of elliptic partial differential equations and their singular perturbations. The smoothed aggregation method is an algebraic multigrid technique where the prolongators are constructed in two steps. First, the tentative prolongator is constructed by the aggregation (or, the generalized aggregation) method. Then, the range of the tentative prolongator is smoothed by a sparse linear...
We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...
We consider high order finite difference approximations to the Helmholtz equation in an exterior domain. We include a simplified absorbing boundary condition to approximate the Sommerfeld radiation condition. This yields a large, but sparse, complex system, which is not self-adjoint and not positive definite. We discretize the equation with a compact fourth or sixth order accurate scheme. We solve this large system of linear equations with a Krylov subspace iterative method. Since the method converges...