Parallel split least‐squares mixed finite element method for parabolic problem

Based on overlapping domain decomposition, a new class of parallel split least‐squares (PSLS) mixed finite element methods is presented for solving parabolic problem. The algorithm is fully parallel. In the overlapping domains, the partition of unity is applied to distribute the corrections reasonably, which makes that the new method only needs one or two iteration steps to reach given accuracy at each time step while the classical Schwarz alternating methods need many iteration steps. The dependence of the convergence rate on the spacial mesh size, time increment, iteration times, and subdomains overlapping degree is analyzed. Some numerical results are reported to confirm the theoretical analysis.     

Overlapping Schwarz methods on unstructured meshes using non-matching coarse grids

Summary. We consider two level overlapping Schwarz domain decomposition methods for solving the finite element problems that arise from discretizations of elliptic problems on general unstructured meshes in two and three dimensions. Standard finite element interpolation from the coarse to the fine grid may be used. Our theory requires no assumption on the substructures that constitute the whole domain, so the substructures can be of arbitrary shape and of different size. The global coarse mesh is allowed to be non-nested to the fine grid on which the discrete problem is to be solved, and neither the coarse mesh nor the fine mesh need be quasi-uniform. In addition, the domains defined by the fine and coarse grid need not be identical. The one important constraint is that the closure of the coarse grid must cover any portion of the fine grid boundary for which Neumann boundary conditions are given. In this general setting, our algorithms have the same optimal convergence rate as the usual two level overlapping domain decomposition methods on structured meshes. The condition number of the preconditioned system depends only on the (possibly small) overlap of the substructures and the size of the coarse grid, but is independent of the sizes of the subdomains. Received March 23, 1994 / Revised version received June 2, 1995     

A general approach to analyse preconditioners for two‐by‐two block matrices

Two‐by‐two block matrices arise in various applications, such as in domain decomposition methods or when solving boundary value problems discretised by finite elements from the separation of the node set of the mesh into ‘fine’ and ‘coarse’ nodes. Matrices with such a structure, in saddle point form arise also in mixed variable finite element methods and in constrained optimisation problems. A general algebraic approach to construct, analyse and control the accuracy of preconditioners for matrices in two‐by‐two block form is presented. This includes both symmetric and nonsymmetric matrices, as well as indefinite matrices. The action of the preconditioners can involve element‐by‐element approximations and/or geometric or algebraic multigrid/multilevel methods. Copyright © 2011 John Wiley & Sons, Ltd.     

Fast parallel solvers for symmetric boundary element domain decomposition equations

Summary. The boundary element method (BEM) is of advantage in many applications including far-field computations in magnetostatics and solid mechanics as well as accurate computations of singularities. Since the numerical approximation is essentially reduced to the boundary of the domain under consideration, the mesh generation and handling is simpler than, for example, in a finite element discretization of the domain. In this paper, we discuss fast solution techniques for the linear systems of equations obtained by the BEM (BE-equations) utilizing the non-overlapping domain decomposition (DD). We study parallel algorithms for solving large scale Galerkin BE–equations approximating linear potential problems in plane, bounded domains with piecewise homogeneous material properties. We give an elementary spectral equivalence analysis of the BEM Schur complement that provides the tool for constructing and analysing appropriate preconditioners. Finally, we present numerical results obtained on a massively parallel machine using up to 128 processors, and we sketch further applications to elasticity problems and to the coupling of the finite element method (FEM) with the boundary element method. As shown theoretically and confirmed by the numerical experiments, the methods are of algebraic complexity and of high parallel efficiency, where denotes the usual discretization parameter. Received August 28, 1996 / Revised version received March 10, 1997     

Local and parallel finite element algorithm based on the partition of unity method for the incompressible MHD flow

Based on the partition of unity method (PUM), a local and parallel finite element method is designed and analyzed for solving the stationary incompressible magnetohydrodynamics (MHD). The key idea of the proposed algorithm is to first solve the nonlinear system on a coarse mesh, divide the globally fine grid correction into a series of locally linearized residual problems on some subdomains derived by a class of partition of unity, then compute the local subproblems in parallel, and obtain the globally continuous finite element solution by assembling all local solutions together by the partition of unity functions. The main feature of the new method is that the partition of unity provide a flexible and controllable framework for the domain decomposition. Finally, the efficiency of our theoretical analysis is tested by numerical experiments.     

利用自然边界归化求解平面弹性方程外边值问题的SCHWARZ算法

1引言平面弹性方程在水利土建等工程技术领域有着广泛应用．其中,孔边应力集中等问题,都是无界区域问题．我们可以通过各种实验手段研究上述问题．而随着计算机和有限元技术的迅猛发展,数值解法提供了一种研究上述问题的有效途径．对于有界区域上的平面弹性方程,我们可以直接利用有限元方法求解,对于其中的大规模问题可以利用区域分解和并行技术求解．但这些方法难以处理无界区域问题．虽然对于某些典型区域上的外问题（例如,圆孔外区域和＿些规则形状裂纹）可以针对具体情况利用复变函数论方法予以解决,但对于一般的无界区域问题广…     

关于解椭圆型问题的两个子区域不重叠区域分解算法

关于解椭圆型问题的两个子区域不重叠区域分解算法顾金生，胡显承（清华大学）ＯＮＴＨＥＤＯＭＡＩＮＤＥＣＯＭＰＯＳＩＴＩＯＮＭＥＴＨＯＤＳＦＯＲＥＬＬＩＰＴＩＣＰＲＯＢＬＥＭＳＷＩＴＨＴＷＯＳＵＢＳＴＲＵＣＴＵＲＥＳ￥ＧｕＪｉｎ－ｓｈｅｎｇ；ＨｕＸｉａｎ...     

无界区域上基于自然边界归化的一种重叠型区域分解法及其离散化

1．引言由于科学技术的迅猛发展，人们遇到许多大规模科学和工程计算问题．随着并行计算机的出现和应用，并行技术越来越得到人们的重视和研究．区域分解法成为并行计算和处理这类问题的主要方法之一．但是，对于无界区域上的椭圆边值问题，因进行区域分解后至少有一个区域仍为无界区域，故仅应用通常的区域分解算法求解是不够的．由于边界归化是处理无界区域问题的有效手段，通常采用边界元和有限元耦合的方法求解此类问题IZ，6。8。121．或片什适当的人工边界并在此边界上加近似边界条件，再在有限区域应用有限元方法求解【人习．近年来…     

不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

无界区域上基于自然边界归化的双调和方程的一种重叠型区域分解法

1．引言双调和方程边值问题的一个力学背景是薄板弯曲问题．对于有界区域上的双调和方程，可以直接利用协调元和非协调元求解18，10］．我们考虑双调和方程Dirichlet外边值问题其中0是充分光滑闭曲线ro之外的无界区域，naro关于0的单位外法向量．引理1．且卜］．若。0EH’／‘（几），gEH‘／‘（fo），则问题（1．1）在W0z（fi）中有唯一解·这里由迹定理，可找到一个具有紧支集的函数识。。，g）E护（炉）满足可。0；g）【F。＝。0和则问题（1．1）等价于如下齐次边值问题其中f—一面‘B有紧支集．边值问题（工．2）又可转化为变分问…     

SOME ESTIMATES WITH NONCONFORMING ELEMENTS IN DOMAIN DECOMPOSITION ANALYSIS

1.IntroductionFOrsimplicityoftheexposition,weconsidertheellipticboundaryvalueproblemonaboundedopenpolygonaldomainfiCEZwhereItiswell--knownthat(1.1)hasauniquesolutionueH'(fl)(of.[7,15,16]).Supposethatfib~{e}isaquajsi--uniformmeshoffi,i.e.,fibsatisfieswhere…     

A conforming finite element method for overlapping and nonmatching grids

In this paper we propose a finite element method for nonmatching overlapping grids based on the partition of unity. Both overlapping and nonoverlapping cases are considered. We prove that the new method admits an optimal convergence rate. The error bounds are in terms of local mesh sizes and they depend on neither the overlapping size of the subdomains nor the ratio of the mesh sizes from different subdomains. Our results are valid for multiple subdomains and any spatial dimensions.

     

Upper bound of the constant in the strengthened C.B.S. inequality for FEM 2D elasticity equations

The basic theory of the strengthened Cauchy–Buniakowskii–Schwarz (C.B.S.) inequality is the main tool in the convergence analysis of the recently proposed algebraic multilevel iterative methods. An upper bound of the constant γ in the strengthened C.B.S. inequality for the case of the finite element solution of 2D elasticity problems is obtained. It is assumed that linear triangle finite elements are used, the initial mesh consisting of right isosceles triangles and the mesh refinement procedure being uniform. For the resulting linear algebraic systems we have proved that γ²<0.75 uniformly on the mesh parameter and on Poisson's ratio ν ? (0, 1/2). Furthermore, the presented numerical tests show that the same relation holds for arbitrary initial right triangulations, even in the case of degeneracy of triangles. The theoretical results obtained are practically important for successful implementation of the finite element method to large-scale modeling of complicated structures. They allow us to construct optimal order algebraic multilevel iterative solvers for a wide class of real–life elasticity problems.     

Two-level Schwarz method for unilateral variational inequalities

The numerical solution of variational inequalities of obstacletype associated with second-order elliptic operators is considered.Iterative methods based on the domain decomposition approachare proposed for discrete obstacle problems arising from thecontinuous, piecewise linear finite element approximation ofthe differential problem. A new variant of the Schwarz methodology,called the two-level Schwarz method, is developed offering thepossibility of making use of fast linear solvers (e.g., linearmultigrid and fictitious domain methods) for the genuinely nonlinearobstacle problems. Namely, by using particular monotonicityresults, the computational domain can be partitioned into (mesh)subdomains with linear and nonlinear (obstacle-type) subproblems.By taking advantage of this domain decomposition and fast linearsolvers, efficient implementation algorithms for large-scalediscrete obstacle problems can be developed. The last part ofthe paper is devoted to illustrate numerical experiments.     

椭圆型方程的重叠型区域分裂混合元方法

本文研究椭圆型方程的重叠型区域分解混合元方法，对第一边值和第二边值问题，分别给出了离散形式的区域分解混合元格式；证明了区域分裂格式解的存在唯一性和算法的收敛性，并给出数值算例．     

Parallel characteristic finite element method for time‐dependent convection–diffusion problem

Based on the overlapping‐domain decomposition and parallel subspace correction method, a new parallel algorithm is established for solving time‐dependent convection–diffusion problem with characteristic finite element scheme. The algorithm is fully parallel. We analyze the convergence of this algorithm, and study the dependence of the convergent rate on the spacial mesh size, time increment, iteration times and sub‐domains overlapping degree. Both theoretical analysis and numerical results suggest that only one or two iterations are needed to reach to optimal accuracy at each time step. Copyright © 2010 John Wiley & Sons, Ltd.     

A parallel finite element variational multiscale method based on fully overlapping domain decomposition for incompressible flows

Based on fully overlapping domain decomposition and a recent variational multiscale method, a parallel finite element variational multiscale method for convection dominated incompressible flows is proposed and analyzed. In this method, each processor computes a local finite element solution in its own subdomain using a global mesh that is locally refined around its own subdomain, where a stabilization term based on two local Gauss integrations is adopted to stabilize the numerical form of the Navier–Stokes equations. Using the technical tool of local a priori estimate for the finite element solution, error bounds of the discrete solution are estimated. Algorithmic parameter scalings are derived. Numerical tests are also given to verify the theoretical predictions and demonstrate the effectiveness of the method. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 856–875, 2015     

SUBSTRUCTURE PRECONDITIONERS FOR NONCONFORMING PLATE ELEMENTS

1.IntroductionInthispaper,wegeneralizetheBPSalgorithm[1]tononconformingelementfproximationsofthebiharmonicequation.WeconstructapreconditionerforMor:elementbysubstructuringonthebasisofafunctiondecompositionfordiscretebibmonicfunctions.Thefunctiondecomposit…     

Mesh and solver co‐adaptation in finite element methods for anisotropic problems

Mesh generation and algebraic solver are two important aspects of the finite element methodology. In this article, we are concerned with the joint adaptation of the anisotropic triangular mesh and the iterative algebraic solver. Using generic numerical examples pertaining to the accurate and efficient finite element solution of some anisotropic problems, we hereby demonstrate that the processes of geometric mesh adaptation and the algebraic solver construction should be adapted simultaneously. We also propose some techniques applicable to the co‐adaptation of both anisotropic meshes and linear solvers. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

基于子结构法构造用非协调元解椭圆型问题的预处理器（Ⅰ）

基于子结构法构造用非协调元解椭圆型问题的预处理器（Ⅰ）顾金生，胡显承（清华大学应用数学系）ＴＨＥＣＯＮＳＴＲＵＣＴＩＯＮＯＦＰＲＥＣＯＮＤＩＴＩＯＮＥＲＳＦＯＲＥＬＬＩＰＴＩＣＰＲＯＢＬＥＭＳＤＩＳＣＲＥＴＩＺＥＤＢＹＮＯＮＣＯＮＦＯＲＭＩＮＧＦＩＮ...