Convergence rate of some domain decomposition methods for overlapping and nonoverlapping subdomains

Summary. Three iterative domain decomposition methods are considered: simultaneous updates on all subdomains (Additive Schwarz Method), flow directed sweeps and double sweeps. By using some techniques of formal language theory we obtain a unique criterion of convergence for the three methods. The convergence rate is a function of the criterion and depends on the algorithm. Received October 24, 1994 / Revised version received November 27, 1995     

Unsteady convection and convection-diffusion problems via direct overlapping domain decomposition methods

In solving unsteady problems,domain decomposition methods may be used either for iterative preconditioning each global implicit time-step or directly (noniteratively) within a blockwise implicit time-stepping procedure, in the latter case, the inner boundary values for the subproblems are generated by explicit time-extrapolation. The overlapping variants of this method have been proved to be efficient tools for solving parabolic and first-order hyperbolic problems on modern parallel computers, because they require global communication only once per time-step. The mechanism making this possible is the exponential decay in space of the time-discrete Green's function. We investigate several model problems of convection and convection-diffusion. Favorable optimal and far-reaching estimates of the overlap required have been established in the case of exemplary standard upwind finite-difference schemes. In particular, it has been shown that the overlap for the convection-diffusion problem is the additive function of overlaps for the corresponding convection and diffusion problem to be considered independently. These results have been confirmed with several numerical test examples. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 387–406, 1998     

Automatic domain decomposition for semiclassical Bohmian mechanics using k-means clustering

We propose a method to automatically decompose domains in the context of semiclassical Bohmian mechanics. The algorithm is based on the approximate quantum potential method and the technique of k-means clustering. Two numerical examples, static analysis of quantum forces for a Pearson Type IV distribution and temporal analysis of the scattering on the Eckart barrier, are presented to show the viability of the method. The first example demonstrates that approximate quantum forces using our domain decomposition technique achieves convergence as the number of domains increases. In the second example, it is demonstrated that the domains constructed from k-means clustering has well adapted themselves to the evolving wave packet, providing coverage to both transmission and reflection waves. We also confirm that the use of multiple domains improves the evolution of the wave packet by comparing the result with the quantum mechanical solution, previously obtained. The computational cost remains manageable even with a naive implementation of time-consuming summation routines, but development of more sophisticated methodology is recommended for large scale, multidimensional calculations.     

Convergence rate estimate for a domain decomposition method

Summary We provide a convergence rate analysis for a variant of the domain decomposition method introduced by Gropp and Keyes for solving the algebraic equations that arise from finite element discretization of nonsymmetric and indefinite elliptic problems with Dirichlet boundary conditions in ². We show that the convergence rate of the preconditioned GMRES method is nearly optimal in the sense that the rate of convergence depends only logarithmically on the mesh size and the number of substructures, if the global coarse mesh is fine enough.This author was supported by the National Science Foundation under contract numbers DCR-8521451 and ECS-8957475, by the IBM Corporation, and by the 3M Company, while in residence at Yale UniversityThis author was supported by the Applied Mathematical Sciences subprogram of the Office of Energy Research, U.S. Department of Energy under Contract W-31-109-Eng-38This author was supported by the National Science Foundation under contract number ECS-8957475, by the IBM Corporation, and by the 3M Company     

Substructured two-grid and multi-grid domain decomposition methods

Two-level Schwarz domain decomposition methods are very powerful techniques for the efficient numerical solution of partial differential equations (PDEs). A two-level domain decomposition method requires two main components: a one-level preconditioner (or its corresponding smoothing iterative method), which is based on domain decomposition techniques, and a coarse correction step, which relies on a coarse space. The coarse space must properly represent the error components that the chosen one-level method is not capable to deal with. In the literature, most of the works introduced efficient coarse spaces obtained as the span of functions defined on the entire space domain of the considered PDE. Therefore, the corresponding two-level preconditioners and iterative methods are defined in volume. In this paper, we use the excellent smoothing properties of Schwarz domain decomposition methods to define, for general elliptic problems, a new class of substructured two-level methods, for which both Schwarz smoothers and coarse correction steps are defined on the interfaces (except for the application of the smoother that requires volumetric subdomain solves). This approach has several advantages. On the one hand, the required computational effort is cheaper than the one required by classical volumetric two-level methods. On the other hand, our approach does not require, like classical multi-grid methods, the explicit construction of coarse spaces, and it permits a multilevel extension, which is desirable when the high dimension of the problem or the scarce quality of the coarse space prevents the efficient numerical solution. Numerical experiments demonstrate the effectiveness of the proposed new numerical framework.

     

Application of domain decomposition methods to isogeometric analysis

During the past decade various new spatial discretization techniques have been developed. In particular, the usage of NURBS based shape functions, well known to the CAD community, has been adapted to finite element technology. In the present work we use the mortar finite element method for the coupling of nonconforming discretized sub-domains in the framework of nonlinear elasticity. We show that the method can be applied to isogeometric analysis with little effort, once the framework of NURBS based shape functions has been implemented. Furthermore, a specific coordinate augmentation technique allows the design of an energy-momentum scheme for the constrained mechanical system under consideration. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

General theory of domain decomposition: Indirect methods

According to a general theory of domain decomposition methods (DDM), recently proposed by Herrera, DDM may be classified into two broad categories: direct and indirect (or Trefftz‐Herrera methods). This article is devoted to formulate systematically indirect methods and apply them to differential equations in several dimensions. They have interest since they subsume some of the best‐known formulations of domain decomposition methods, such as those based on the application of Steklov‐Poincaré operators. Trefftz‐Herrera approach is based on a special kind of Green's formulas applicable to discontinuous functions, and one of their essential features is the use of weighting functions which yield information, about the sought solution, at the internal boundary of the domain decomposition exclusively. A special class of Sobolev spaces is introduced in which boundary value problems with prescribed jumps at the internal boundary are formulated. Green's formulas applicable in such Sobolev spaces, which contain discontinuous functions, are established and from them the general framework for indirect methods is derived. Guidelines for the construction of the special kind of test functions are then supplied and, as an illustration, the method is applied to elliptic problems in several dimensions. A nonstandard method of collocation is derived in this manner, which possesses significant advantages over more standard procedures. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 296–322, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10008     

Convergence of alternating domain decomposition schemes for kinetic and aerodynamic equations

A domain decomposition scheme linking linearized kinetic and aerodynamic equations is considered. Convergence of the alternating scheme is shown. Moreover the physical correctness of the obtained coupled solutions is discussed.     

The multipliers‐free domain decomposition methods

This article concludes the development and summarizes a new approach to dual‐primal domain decomposition methods (DDM), generally referred to as “the multipliers‐free dual‐primal method.” Contrary to standard approaches, these new dual‐primal methods are formulated without recourse to Lagrange‐multipliers. In this manner, simple and unified matrix‐expressions, which include the most important dual‐primal methods that exist at present are obtained, which can be effectively applied to floating subdomains, as well. The derivation of such general matrix‐formulas is independent of the partial differential equations that originate them and of the number of dimensions of the problem. This yields robust and easy‐to‐construct computer codes. In particular, 2D codes can be easily transformed into 3D codes. The systematic use of the average and jump matrices, which are introduced in this approach as generalizations of the “average” and “jump” of a function, can be effectively applied not only at internal‐boundary‐nodes but also at edges and corners. Their use yields significant advantages because of their superior algebraic and computational properties. Furthermore, it is shown that some well‐known difficulties that occur when primal nodes are introduced are efficiently handled by the multipliers‐free dual‐primal method. The concept of the Steklov–Poincaré operator for matrices is revised by our theory and a new version of it, which has clear advantages over standard definitions, is given. Extensive numerical experiments that confirm the efficiency of the multipliers‐free dual‐primal methods are also reported here. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

General theory of domain decomposition: Beyond Schwarz methods

Recently, Herrera presented a general theory of domain decomposition methods (DDM). This article is part of a line of research devoted to its further development and applications. According to it, DDM are classified into direct and indirect, which in turn can be subdivided into overlapping and nonoverlapping. Some articles dealing with general aspects of the theory and with indirect (Trefftz–Herrera) methods have been published. In the present article, a very general direct‐overlapping method, which subsumes Schwarz methods, is introduced. Also, this direct‐overlapping method is quite suitable for parallel implementation. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 495–517, 2001     

Preconditioning for radial basis functions with domain decomposition methods

In our previous work, an effective preconditioning scheme that is based upon constructing least-squares approximation cardinal basis functions (ACBFs) from linear combinations of the RBF-PDE matrix elements has shown very attractive numerical results. The preconditioner costs O(N²) flops to set up and O(N) storage. The preconditioning technique is sufficiently general that it can be applied to different types of different operators. This was applied to the 2D multiquadric method, with c~1/√N on the Poisson test problem, the preconditioned GMRES converges in tens of iterations. In this paper, we combine the RBF methods and the ACBF preconditioning technique with the domain decomposition method (DDM). We studied different implementations of the ACBF-DDM scheme and provide numerical results for N > 10,000 nodes. We shall demonstrate that the efficiency of the ACBF-DDM scheme improves dramatically as successively finer partitions of the domain are considered.     

New constructions of domain decomposition methods for systems of PDEs

We propose new domain decomposition methods for systems of partial differential equations in two and three dimensions. The algorithms are derived with the help of the Smith factorization. This could also be validated by numerical experiments. To cite this article: V. Dolean et al., C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

Convergence analysis of domain decomposition algorithms with full overlapping for the advection-diffusion problems

The aim of this paper is to study the convergence properties of a time marching algorithm solving advection-diffusion problems on two domains using incompatible discretizations. The basic algorithm is first described, and theoretical and numerical results that illustrate its convergence properties are then presented.

     

Lagrangian formulation of domain decomposition methods: A unified theory

Domain decomposition methods based on one Lagrange multiplier have been shown to be very efficient for solving ill-conditioned problems in parallel. Several variants of these methods have been developed in the last ten years. These variants are based on an augmented Lagrangian formulation involving one or two Lagrange multipliers and on mixed type interface conditions between the sub-domains. In this paper, the Lagrangian formulations of some of these domain decomposition methods are presented both from a continuous and a discrete point of view.     

Contractivity of domain decomposition splitting methods for nonlinear parabolic problems

This work deals with the efficient numerical solution of nonlinear parabolic problems posed on a two-dimensional domain Ω. We consider a suitable decomposition of domain Ω and we construct a subordinate smooth partition of unity that we use to rewrite the original equation. Then, the combination of standard spatial discretizations with certain splitting time integrators gives rise to unconditionally contractive schemes. The efficiency of the resulting algorithms stems from the fact that the calculations required at each internal stage can be performed in parallel.     

Solution of the 2-D eddy current problem via the domain decomposition methods on serial and parallel computers

Domain decomposition algorithms are applied to the solution of a time harmonic two-dimensional eddy current problem. The system of differential equations describing this problem is considered as a singularly perturbed problem. An iterative domain decomposition algorithm suitable for parallelization is described, and convergence of this algorithm is established. The implementation on a shared memory multiprocessor is described, and numerical experiments are presented.     

Rate of Convergence for some constraint decomposition methods for nonlinear variational inequalities

Summary. Some general subspace correction algorithms are proposed for a convex optimization problem over a convex constraint subset. One of the nontrivial applications of the algorithms is the solving of some obstacle problems by multilevel domain decomposition and multigrid methods. For domain decomposition and multigrid methods, the rate of convergence for the algorithms for obstacle problems is of the same order as the rate of convergence for jump coefficient linear elliptic problems. In order to analyse the convergence rate, we need to decompose a finite element function into a sum of functions from the subspaces and also satisfying some constraints. A special nonlinear interpolation operator is introduced for decomposing the functions. Received December 13, 2001 / Revised version received February 19, 2002 / Published online June 17, 2002 This work was partially supported by the Norwegian Research Council under projects 128224/431 and SEP-115837/431.     

A relaxation procedure for domain decomposition methods using finite elements

Summary We present the convergence analysis of a new domain decomposition technique for finite element approximations. This technique was introduced in [11] and is based on an iterative procedure among subdomains in which transmission conditions at interfaces are taken into account partly in one subdomain and partly in its adjacent. No global preconditioner is needed in the practice, but simply single-domain finite element solvers are required. An optimal strategy for an automatic selection of a relaxation parameter to be used at interface subdomains is indicated. Applications are given to both elliptic equations and incompressible Stokes equations.     

Adaptive wavelet frame domain decomposition methods for nonlinear elliptic equations

In this article, we are concerned with the numerical treatment of nonlinear elliptic boundary value problems. Our method of choice is a domain decomposition strategy. Partially following the lines from (Cohen, Dahmen and deVore, SIAM J Numer Anal 41 (2003), 1785–1823; Kappei, Appl Anal J Sci 90 (2011), 1323–1353; Lui, SIAM J Sci Comput 21 (2000), 1506–1523; Stevenson and Werner, Math Comp 78 (2009), 619–644), we develop an adaptive additive Schwarz method using wavelet frames. We show that the method converges with an asymptotically optimal rate and support our theoretical results with numerical tests in one and two space dimensions. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013