求解Oseen流的交替线松弛多重网格方法

利用Riemann解的通量差分分裂法——Godunov方法对Oseen流控制方程进行离散,得到了基于一阶上迎风格式的离散方程,并给出了使用多重网格方法求解该离散方程的V-循环算法和W-循环算法的收敛性分析.通过局部Fourier分析方法,对获得的离散方程的聚对称交替线GaussSeidel松弛的光滑性质进行了研究.结果表明:使用多重网格的两层网格及三层网格算法求解具有不同Reynolds数的Oseen流,即便是在高Reynolds数情况下,聚对称交替线Gauss-Seidel松弛具有很好的光滑性质,多重网格W-循环算法收敛性比V-循环算法好.     

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

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

无限元多重网格算法

在求偏微分方程数值解时,往往需要解一个规模很大的代数方程组,而多重网格是一种十分有效的迭代方法.大量数值试验证明,它具有很高的收敛速度.理论分析表明,这种迭代法的收敛速度并不随网格的加密而降低,这一突出优点是其它迭代方法望尘莫及的. 在使用有限元多重网格算法时,如果区域边界的角点使解具有奇性,理论分析会遇到     

用自适应MG法,FMG法求解Helmholtz方程

本文针对Ｖ循环、Ｗ循环和多重网格法中最优光滑次数及循环体个数难以确定的缺点，以Ｈｅｌｍｈｏｌｔｚ方程为例给出自适应的多重网格算法和自适应的完全多重网格算法。     

Wilson非协调元的V循环多重网格法 *

非协调有限元V循环多重网格法的收敛性至今仍是一个没有很好解决的问题 .给出了Wilson非协调有限元的两类V循环多重网格法的收敛性证明 .     

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

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

等几何分析的多重网格共轭梯度法

提高NURBS基函数阶数可以提高等几何分析的精度,同时也会降低多重网格迭代收敛速度.将共轭梯度法与多重网格方法相结合,提出了一种提高收敛速度的方法,该方法用共轭梯度法作为基础迭代算法,用多重网格进行预处理.对Poisson(泊松)方程分别用多重网格方法和多重网格共轭梯度法进行了求解,计算结果表明:等几何分析中采用高阶NURBS基函数处理三维问题时,多重网格共轭梯度法比多重网格法的收敛速度更快.     

曲边区域上的多边形网格间断有限元离散及其多重网格算法

本文利用多边形网格上的间断有限元方法离散二阶椭圆方程,在曲边区域上,采用多条直短边逼近曲边的以直代曲的策略,实现了高阶元在能量范数下的最优收敛.本文还将这一方法用于带曲边界面问题的求解,同样得到高阶元的最优收敛.此外我们还设计并分析了这一方法的\linebreakW-cycle和Variable V-cycle多重网格预条件方法,证明当光滑次数足够多时,多重网格预条件算法一致收敛.最后给出了数值算例,证实该算法的可行性并验证了理论分析的结果.     

几乎不可压缩线性弹性问题的多重网格Uzawa型混合有限元方法

该文针对几乎不可压缩弹性问题,设计了多重网格Uzawa型混合有限元方法,成功克服了"闭锁"现象.通过引入"压力"变量p将弹性问题转化为一个鞍点型系统,对该系统将Uzawa型迭代法和多重网格方法相结合,建立了多重网格和套迭代多重网格Uzawa型混合有限元方法,并给出了该算法的收敛性.数值算例验证了方法的有效性和稳定性.     

旋转Q1非协调元的V循环多重网格法

1．引言近年来,多重网格法已成为行之有效的偏微分方程数值解法,而对非协调元的多重网格法也有众多的研究,例在[1,3］中,作者研究了非协P1元的w循环多重网格法,[10］中,作者研究了＊11s。n非协调元的V循环多重网格法．此外在K豆1,12］中,作者研究了板问题非协调有限元的多重网格法．最近,Rannacher和Turek同构造了所谓的QI非协调元,并用该元离散StokeS问题．而在问中,利用该元来计算晶体,数值效果非常好．同时在同中,作者给出了该元的误差估计和超收敛分析．最近,Chen和oswald同又讨论了该元的多重网格法,并证明了W循环…     

MULTIGRID ALGORITHM FOR THE COUPLING SYSTEM OF NATURAL BOUNDARY ELEMENT METHOD AND FINITE ELEMENT METHOD FOR UNBOUNDED DOMAIN PROBLEMS

In this paper, some V-cycle multigrid algorithms are presented for the coupling system arising from the discretization of the Dirichlet exterior problem by coupling the natural boundary element method and finite element method. The convergence of these multigrid algorithms is obtained even with only one smoothing on all levels. The rate of convergence is found uniformly bounded independent of the number of levels and the mesh sizes of all levels, which indicates that these multigrid algorithms are optimal. Some numerical results are also reported.     

A new type of full multigrid method for the elasticity eigenvalue problem

This paper introduces a new type of full multigrid method for the elasticity eigenvalue problem. The main idea is to avoid solving large scale elasticity eigenvalue problem directly by transforming the solution of the elasticity eigenvalue problem into a series of solutions of linear boundary value problems defined on a multilevel finite element space sequence and some small scale elasticity eigenvalue problems defined on the coarsest correction space. The involved linear boundary value problems will be solved by performing some multigrid iterations. Besides, some efficient techniques such as parallel computing and adaptive mesh refinement can also be absorbed in our algorithm. The efficiency and validity of the multigrid methods are verified by several numerical experiments.     

MONOLITHIC MULTIGRID FOR REDUCED MAGNETOHYDRODYNAMIC EQUATIONS

In this paper,the monolithic multigrid method is investigated for reduced magnetohydrodynamic equations.We propose a diagonal Braess-Sarazin smoother for the finite element discrete system and prove the uniform convergence of the MMG method with respect to mesh sizes.A multigrid-preconditioned FGMRES method is proposed to solve the magnetohydrodynamic equations.It turns out to be robust for relatively large physical parameters.By extensive numerical experiments,we demonstrate the optimality of the monolithic multigrid method with respect to the number of degrees of freedom.     

Multigrid for the mortar element method for P1 nonconforming element

In this paper, a multigrid algorithm is presented for the mortar element method for P1 nonconforming element. Based on the theory developed by Bramble, Pasciak, Xu in [5], we prove that the W-cycle multigrid is optimal, i.e. the convergence rate is independent of the mesh size and mesh level. Meanwhile, a variable V-cycle multigrid preconditioner is constructed, which results in a preconditioned system with uniformly bounded condition number. Received May 11, 1999 / Revised version received April 1, 2000 / Published online October 16, 2000     

Multigrid fictitious boundary and grid deformation methods for flow-induced rotation of wing

Numerical simulations of flow-induced rotation of wing by multigrid fictitious boundary and grid deformation methods are presented. The flow is computed by a special ALE formulation with a multigrid finite element solver. The solid wing is allowed to move freely through the computational mesh which is adaptively aligned by a special mesh deformation method. The advantage of this approach is that no expensive remeshing has to be performed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An efficient algebraic multigrid method for quadratic discretizations of linear elasticity problems on some typical anisotropic meshes in three dimensions

The quality of the mesh used in the finite element discretizations will affect the efficiency of solving the discreted linear systems. The usual algebraic solvers except multigrid method do not consider the effect of the grid geometry and the mesh quality on their convergence rates. In this paper, we consider the hierarchical quadratic discretizations of three‐dimensional linear elasticity problems on some anisotropic hexahedral meshes and present a new two‐level method, which is weakly independent of the size of the resulting problems by using a special local block Gauss–Seidel smoother, that is LBGS_v iteration when used for vertex nodes or LBGS_m iteration for midside nodes. Moreover, we obtain the efficient algebraic multigrid (AMG) methods by applying DAMG (AMG based on distance matrix) or DAMG‐PCG (PCG with DAMG as a preconditioner) to the solution of the coarse level equation. The resulting AMG methods are then applied to a practical example as a long beam. The numerical results verify the efficiency and robustness of the proposed AMG algorithms. Copyright © 2015 John Wiley & Sons, Ltd.     

Convergence analysis of HSS-multigrid methods for second-order nonselfadjoint elliptic problems

In this paper, the multigrid methods using Hermitian/skew-Hermitian splitting (HSS) iteration as smoothers are investigated. These smoothers also include the modified additive and multiplicative smoothers which result from subspace decomposition. Without full elliptic regularity assumption, it is shown that the multigrid methods with these smoothers converge uniformly for second-order nonselfadjoint elliptic boundary value problems if the mesh size of the coarsest grid is sufficiently small (but independent of the number of the multigrid levels). Numerical results are reported to confirm the theoretical analysis.     

Comparison of multigrid algorithms for high‐order continuous finite element discretizations

We present a comparison of different multigrid approaches for the solution of systems arising from high‐order continuous finite element discretizations of elliptic partial differential equations on complex geometries. We consider the pointwise Jacobi, the Chebyshev‐accelerated Jacobi, and the symmetric successive over‐relaxation smoothers, as well as elementwise block Jacobi smoothing. Three approaches for the multigrid hierarchy are compared: (1) high‐order h‐multigrid, which uses high‐order interpolation and restriction between geometrically coarsened meshes; (2) p‐multigrid, in which the polynomial order is reduced while the mesh remains unchanged, and the interpolation and restriction incorporate the different‐order basis functions; and (3) a first‐order approximation multigrid preconditioner constructed using the nodes of the high‐order discretization. This latter approach is often combined with algebraic multigrid for the low‐order operator and is attractive for high‐order discretizations on unstructured meshes, where geometric coarsening is difficult. Based on a simple performance model, we compare the computational cost of the different approaches. Using scalar test problems in two and three dimensions with constant and varying coefficients, we compare the performance of the different multigrid approaches for polynomial orders up to 16. Overall, both h‐multigrid and p‐multigrid work well; the first‐order approximation is less efficient. For constant coefficients, all smoothers work well. For variable coefficients, Chebyshev and symmetric successive over‐relaxation smoothing outperform Jacobi smoothing. While all of the tested methods converge in a mesh‐independent number of iterations, none of them behaves completely independent of the polynomial order. When multigrid is used as a preconditioner in a Krylov method, the iteration number decreases significantly compared with using multigrid as a solver. Copyright © 2015 John Wiley & Sons, Ltd.     

PRECONDITIONING HIGHER ORDER FINITE ELEMENT SYSTEMS BY ALGEBRAIC MULTIGRID METHOD OF LINEAR ELEMENTS

We present and analyze a robust preconditioned conjugate gradient method for the higher order Lagrangian finite element systems of a class of elliptic problems. An auxiliary linear element stiffness matrix is chosen to be the preconditioner for higher order finite elements. Then an algebraic multigrid method of linear finite element is applied for solving the preconditioner. The optimal condition number which is independent of the mesh size is obtained. Numerical experiments confirm the efficiency of the algorithm.     

Uniform convergence of multigrid V-cycle on adaptively refined finite element meshes for second order elliptic problems

In this paper we prove the uniform convergence of the standard multigrid V-cycle algorithm with the Gauss-Seidel relaxation performed only on the new nodes and their "immediate" neighbors for discrete elliptic problems on the adaptively refined finite element meshes using the newest vertex bisection algorithm. The proof depends on sharp estimates on the relationship of local mesh sizes and a new stability estimate for the space decomposition based on the Scott-Zhang interpolation operator. Extensive numerical results are reported, which confirm the theoretical analysis.