Search

Abstract: Finite Element Tearing and Interconnecting (FETI) methods are a powerful approach to designing solvers for large-scale problems in computational mechanics. The numerical simulation problem is subdivided into a number of independent sub-problems, which are then coupled in appropriate ways. NURBS- (Non-Uniform Rational B-spline) based isogeometric analysis (IGA) applied to complex geometries requires to represent the computational domain as a collection of several NURBS geometries. Since there is a natural decomposition of the computational domain into several subdomains, NURBS-based IGA is particularly well suited for using FETI methods. This paper proposes the new IsogEometric Tearing and Interconnecting (IETI) method, which combines the advanced solver design of FETI with the exact geometry representation of IGA. We describe the IETI framework for two classes of simple model problems (Poisson and linearized elasticity) and discuss the coupling of the subdomains along interfaces (both for matching interfaces and for interfaces with T-joints, i.e. hanging nodes). Special attention is paid to the construction of a suitable preconditioner for the iterative linear solver used for the interface problem. We report several computational experiments to demonstrate the performance of the proposed IETI method. [Copyright &y& Elsevier]

Abstract: Finite element tearing and interconnecting (FETI) methods and boundary element tearing and interconnecting (BETI) methods are special iterative substructuring methods with Lagrange multipliers. For elliptic boundary value problems on bounded domains, the condition number of these methods can be rigorously bounded by , where H is the subdomain diameter and h the mesh size. The constant C is independent of H, h and possible jumps in the coefficients of the partial differential equation. In certain situations, e.g., in electromagnetic field computations, instead of imposing artificial boundary conditions one may be interested in modelling the real physical behaviour in an exterior domain with a radiation condition. In this work we analyze one-level BETI methods for such unbounded domains and show explicit condition number estimates similar to the one above. Our theoretical results are confirmed in numerical experiments. [Copyright &y& Elsevier]

In this article, we discuss domain decomposition parallel iterative solvers for highly heterogeneous problems of flow and transport in porous media. We are particularly interested in highly unstructured coefficient variation where standard periodic or stochastic homogenization theory is not applicable. When the smallest scale at which the coefficient varies is very small, it is often necessary to scale up the equation to a coarser grid to make the problem computationally feasible. Standard upscaling or multiscale techniques require the solution of local problems in each coarse element, leading to a computational complexity that is at least linear in the global number N of unknowns on the subgrid. Moreover, except for the periodic and the isotropic random case, a theoretical analysis of the accuracy of the upscaled solution is not yet available. Multilevel iterative methods for the original problem on the subgrid, such as multigrid or domain decomposition, lead to similar computational complexity (i.e. O(N)) and are therefore a viable alternative. However, previously no theory was available guaranteeing the robustness of these methods to large coefficient variation. We review a sequence of recent papers where simple variants of domain decomposition methods, such as overlapping Schwarz and one-level FETI, are proposed that are robust to strong coefficient variation. Moreover, we also extend the theoretical results, for the first time, to the dual-primal variant of FETI. [ABSTRACT FROM AUTHOR]

The present paper deals with the numerical solution of boundary value problems for nonlinear potential equations of the form –∇· [ν(|∇ u|) ∇ u] = f arising e. g. in 2D magnetic field computations. We apply Newton's method to the nonlinear variational formulation and solve the linearized problems with coupled Finite and Boundary Element Tearing and Interconnecting (FETI/BETI) methods. The spectrum of the subdomain Jacobi matrices may show high variation due to the nonlinear behavior of the reluctivity on ferromagnetic subdomains. A special preconditioner is proposed to overcome these problems. In order to get a good initial guess for the Newton iteration, we use the nested iteration strategy based on a hierarchy of nested grids. Finally, we discuss our first numerical results obtained from the numerical solution of some 2D magnetostatic model problems. [ABSTRACT FROM AUTHOR]