非协调Wilson有限元的多重网格方法

§1.引言 非协调Wilson有限元[1—3]对解弹性力学方程有实用价值,在工程上有用。本文分析Wilson元的多重网格法,给出用多重网格方法求得的近似解按L~2模和能量模的最佳收敛阶误差估计。对于W-循环,可以证明其计算量与离散空间的维数为同一量级O(N_k)。 考虑二阶椭圆Dirchlet边值问题:     

混合有限元非嵌套网格型预处理法

1 引言 众所周知,二阶椭圆型问题混合有限元离散以后的矩阵是不定的,所以对混合法很难形成一种有效的区域分解法,在文[9]、[10]、[11]中提出了一些混合有限元方法的区域分解法,但在实际计算中有很多局限性。最近Chen对混合有限元法提出一种全新的解释并把它应用到多重网格法中,他的基本思想是混合有限元离散的代数系统实际上等价于某个非协调有限元离散的代数系统,这样可把一个不定问题转化为一个正定问题,本文将基于这种思想考虑混合有限元的区域分解法。 若按传统的Dryia-widlund两水平加性Schwarz方法,要求两层网格间具有嵌套关系,这样在应用中将带来很大的不便。本文将不要求粗网格嵌入细网格中,减少两层网格间的     

线性抛物型方程的多重网格方法

线性抛物型方程的多重网格方法蔚喜军（北京应用物理与计算数学研究所，计算物理实验室）ＭＵＬＴＩＧＲＩＤＭＥＴＨＯＤＦＯＲＴＨＥＬＩＮＥＡＲＰＡＲＡＢＯＬＩＣＰＲＯＢＬＥＭ￥ＹｕＸｉ－ｊｕｎ（ＬａｂｏｒｕｔｏｒｙｏｆＣｏｍｐｕｔａｔｉｏｎａｌＰｈｙｓｉｃ...     

非嵌套网格上的Morley元两水平加性Schwarz方法

1．引言设0是RZ中的有界多边形区域，其边界为Rfl．考虑下面的重调和Dirichlet问题：（1．1）的变分形式为：求。EHI（fi）使得对？／EL‘（m，问题（1．幻的唯一可解性可由冯（m上的M线性型的强制性和连续性以及La。Mlgram定理得出（of［4］）．令人一｛丸）是n的一个三角剖分，并且满足最小角条件，其中h是它的网格参数．设Vh为Money元空间［41．问题（1．2）的有限元离散问题为：求。eVh使得当有限元参数人很小时，这个方程组很大，而且矩阵A的条件数变得非常大，直接求解，存贮量及计算量都很大．如果B可逆，则方程组（1．4）等…     

双参数Morley有限元多重网格法

1、引言 多重网格方法是求解偏微分方程的高效快速算法,在实际中得到广泛应用．[2][6]中考察了Morley元的多重网格方法,并用于双调和方程问题。     

二次Lagrangian有限元方程的几何多重网格法

为了构造快速求解二次Lagrangian有限元方程的几何多重网格法,在选择二次Lagrangian有限元空间和一系列线性Lagrangian有限元空间分别作为最细网格层和其余粗网格层以及构造一种新限制算子的基础上,提出了一种新的几何多重网格法,并对它的计算量进行了估计.数值实验结果,与通常的几何多重网格法和AMG01法相比,表明了新算法计算量少且稳健性强.     

抛物问题的mortar有限元逼近的瀑布型多重网格法

用瀑布型多重网格法解决椭圆、抛物问题，已有不少研究工作[1-2]，本文对抛物问题的mortar有限元的全离散格式提出瀑布型多重网格法，证明了该方法是最优的，即具有最优精确度和复杂度．     

半线性椭圆型问题Mortar有限元逼近的瀑布型多重网格法

Mortar有限元法作为一个非协调的区域分解技术已得到许多研究者的关注(如文献[2]、[5]等)。本文对半线性椭圆型问题的Mortar有限元逼近提出了瀑布型多重网格法,并给出了此法的误差估计和计算复杂度估计定理。     

The Mortar Element Method with Locally Nonconforming Elements

We consider a discretization of linear elliptic boundary value problems in 2-D by the new version of the mortar finite element method which uses locally nonconforming Crouzeix-Raviart elements. We show that if a solution of the original differential problem belongs to the space H ²(), then an error is of the same order as in the standard nonconforming finite element method. We also propose an additive Schwarz method of solving the discrete problem and show that its rate of convergence is almost optimal.     

Some remarks on solutions to an overdetermined elliptic problem in divergence form in a ball

We prove the existence of a continuum of non-radial pairs (k,u) solutions to the following overposed problem div in B _r , u = 0 and on ∂B _r , where B _r is the Euclidean ball centered at zero of radius r in . Dedicato a Marco e Andrea Provera.     

Multigrid for second‐order ADN elliptic systems

This paper theoretically examines a multigrid strategy for solving systems of elliptic partial differential equations (PDEs) introduced in the work of Lee. Unlike most multigrid solvers that are constructed directly from the whole system operator, this strategy builds the solver using a factorization of the system operator. This factorization is composed of an algebraic coupling term and a diagonal (decoupled) differential operator. Exploiting the factorization, this approach can produce decoupled systems on the coarse levels. The corresponding coarse‐grid operators are in fact the Galerkin variational coarsening of the diagonal differential operator. Thus, rather than performing delicate coarse‐grid selection and interpolation weight procedures on the original strongly coupled system as often done, these procedures are isolated to the diagonal differential operator. To establish the theoretical results, however, we assume that these systems of PDEs are elliptic in the Agmon–Douglis–Nirenberg (ADN) sense and apply the factorization and multigrid only to the principal part of the system of PDEs. Two‐grid error bounds are established for the iteration applied to the complete system of PDEs. Numerical results are presented to illustrate the effectiveness of this strategy and to expose factors that affect the convergence of the methods derived from this strategy.     

Remarks on Multilevel Bases for Divergence-Free Finite Elements

The connection between the multilevel factorization method recently proposed by Sarin and Sameh for solving mixed discretizations of the Stokes equation using a divergence-free finite element formulation, and hierarchical basis preconditioners for the Poisson problem is established. For the 2D triangular Taylor–Hood element, a preconditioner is proposed that could be useful in fractional step methods.     

A Fully Nonlinear Boundary value Problem for the Laplace Equation in Dimension Two

We study the regularity up to the boundary of solutions to the boundary value problem:[math001] in D, ∣?u∣= g on &;pardD, where D is the unit disc. This problem finds its application in the study of geophysical and geomagnetic surveys. If g?C^,[math001](D) and is strictly positive, we prove that uis in the Holder class C1,α(D). An example shows that this is no longer true if g has some zeroes on ?D. In this case u isproved to be of class C¹(D)     

Semiregular Finite Elements in Solving Some Nonlinear Problems

In this paper, under the maximum angle condition, the finite element method is analyzed for nonlinear elliptic variational problem formulated in [4]. In [4] the analysis was done under the minimum angle condition.     

Finite Elements and characteristics for some parabolic-hyperbolic problems

This paper deals with the numerical integration of partial differential equations of the advection-diffusion type when the advection dominates the diffusion. It is shown that finite differencing of the total derivatives yields schemes which do not require upwinding. The method is numerically tested on three problems: the advection diffusion linear problem, the Navier-Stokes equation and the Euler equations.     

Uniform Convergence of Multigrid V-Cycle on Adaptively Refined Finite Element Meshes for Elliptic Problems with Discontinuous Coefficients

The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered. Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm, some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours. The multigrid V-cycle algorithm uses $\mathcal{O}(N)$ operations per iteration and is optimal.     

Finite element approximation of diffusion equations with convolution terms

Approximation of solutions to diffusion equations with memory represented by convolution integral terms is considered. Such problems arise from modeling of flows in fissured media. Convergence of the method is proved and results of numerical experiments confirming the theoretical results are presented. The advantages of implementation of the algorithm in a multiprocessing environment are discussed.

     

半线性问题的瀑布型多重网格法

本文提出了求解半线性椭圆问题的一类新的瀑布型多重网格法，在网格层数固定的条件下证明了此法的最优阶收敛性。