Lagrange四边形单位分解有限元法的最优误差分析

用构造最优局部逼近空间的方法对Lagrange型四边形单位分解有限元法进行了最优误差分析.单位分解取Lagrange型四边形上的标准双线性基函数,构造了一个特殊的局部多项式逼近空间,给出了具有2阶再生性的Lagrange型四边形单位分解有限元插值格式,从而得到了高于局部逼近阶的最优插值误差.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

节点应力连续的四边形单元

节点应力连续的四边形单元Q4-CNS是一种基于单位分解理论的混合的有限元无网格法.Q4-CNS可以视作FE-LSPIM QUAD4的发展.Q4-CNS形函数的导数在节点处是连续的,因此可以自然的得到节点应力,而不需要使用节点应力磨平算法.数值实验表明,与传统四边形单元(QUAD4)相比,Q4-CNS具有更好的计算精度和更高的收敛速度.在扭曲网格下,Q4-CNS也能取得满意的数值精度.然而,QUAD4的数值精度则会随着网格的扭曲明显的变差.基于Kirchhoff-Love假设的非协调板单元计算中,不仅要求形函数在单元的交界面上要保持C0连续性,而且要求形函数在节点处具有C1连续性,所以在任意的四边形单元上构造满足插值条件的非协调板单元形函数较为困难.Q4-CNS形函数的导数在节点处是连续的,所以Q4-CNS在求解基于Kirchhoff-Love假设的板单元问题中具有潜在的应用价值.     

基于插值方法构造对称的矩形板单元

1　引　　言在有限元方法中,构造多项式类有限元的问题可以归结为多元多项式插值问题.在多元多项式插值的情形中,插值条件与插值多项式空间之间存在所谓的“匹配”(correct)问题,参见deBoor的论文[9].在构造非协调有限元时,关键的是设计适当的有限元形参数和适当的有限元形函数空间,再设计与之匹配的形参数.在经典的有限元构造方法中,例如,基于位移假设的板元的构造,就是首先定义好形函数空间.在这种情况下,所设计形参数必须满足两个条件:(1)形参数作为插值条件,必须与事先给定的形函数空间是插值匹配的,也即,由形参数定义的插值条件在形函…     

Navier-Stokes方程的局部投影稳定化方法

将Matthies,Skrzypacz和Tubiska的思想从线性的Oseen方程拓展到了非线性的Navier-Stokes方程,针对不可压缩的定常Navier-Stokes方程,提出了一种局部投影稳定化有限元方法．该方法既克服了对流占优,又绕开了inf-sup条件的限制．给出的局部投影空间既可以定义在两种不同网格上,又可以定义在相同网格上．与其他两级方法相比,定义在同一网格空间上的局部投影稳定化格式更紧凑．在同一网格上,除了给出需要bubble函数来增强的逼近空间外,还特别考虑了两种不需要用bubble函数来增强的新的空间．基于一种特殊的插值技巧,给出了稳定性分析和误差估计．最后,还列举了两个数值算例,进一步验证了理论结果的正确性．     

有限元分层快速高精度算法

1.引言 多重网格法和区域分解法实质是有限元空间的分解,在子空间上实行逐次校正迭代或并行校正迭代。[1]对一维有限元空间,利用正交化过程,消去单元内结点,修改单元角点的基函数,提出了所谓快速高精度算法。实例表明,这一算法十分有效。本文对一般区域Ω R~d(d=1,2,3)上有限元空间进行分层正交分解,提出所谓分层快速高精度算法。     

半直线上三阶方程的Legendre-Laguerre耦合谱元法

针对建立在半直线上的三阶微分方程,提出Legendre-Laguerre耦合谱元法.通过构造满足试探函数空间和检验函数空间的基函数,分解得到的线性系统的系数矩阵是稀疏的,可以有效地进行求解.数值例子验证了方法的有效性和高精度.     

一阶双曲型方程组的时空间断全离散有限元的收敛性

利用能量方法和单元正交分析方法，构造了特殊的Radau型单元正交展开和张量积分解，简明论证了一阶双曲方程组时空间断有限元的收敛性，得到了丰满阶的整体误差估计．数值实验证实了这些理论结果．     

一种h型自适应有限单元

ｈ型自适应有限单元在网格局部细划时．会产生非常规节点，从而破坏了一般意义上的单元连续性假定．本文利用参照节点对非常规单元进行坐标和位移插值．为保证单元之间坐标和位移的连续性，本文提出了一组修正的形函数，常用的形函数是它的一个特例．本方法应用于有限元程序时，除形函数外无须做任何改动．算例表明水文的方法具有方法简单、精度高、自由度少、计算量小等优点．     

W空间中有界线性算子的最佳逼近及其应用

１．引言 对线性算子的有限秩算子逼近是最经典的问题．并且它的应用极广．如数值积分公式、函数的逼近、数值原函数、方程的数值解法等．１９８６年,在文［１］中,首次给出了在再生核空间中函数的最佳逼近算子（恒等算子的有限秩算子逼近）．之后;在文［２］中给出了数值原函数．又在文［３］、［５］、［６］等中利用有限秩算子逼近（并非是最佳逼近）给出了一些方程的数值解法．但这些讨论都是在一元函数空间上只对特殊算子进行的．１９９７年,虽然在文［４］中给出了完备的二元再生核空间及二元函数的最佳逼近插值算子．但是对多元…     

A new family of high regularity elements

In this article, we propose a new family of high regularity finite element spaces. The global approximation spaces are obtained in two steps. We first build an open cover of the computational domain and local approximation spaces on each patch of the cover. Then we construct partition of unity functions subordinate to the open cover depending on the regularity requirement. The basis functions of the global space is given by the products of the local basis functions and the corresponding partition of unity functions. The method can be used to construct finite element spaces of any desired regularity. Approximation properties and implementation details are discussed. Numerical examples for the biharmonic equation are presented to show the effectiveness of the proposed method. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 1–16, 2012     

A priori estimates for two multiscale finite element methods using multiple global fields to wave equations

We consider a scalar wave equation with nonseparable spatial scales. If the solution of the wave equation smoothly depends on some global fields, then we can utilize the global fields to construct multiscale finite element basis functions. We present two finite element approaches using the global fields: partition of unity method and mixed multiscale finite element method. We derive a priori error estimates for the two approaches and theoretically investigate the relation between the smoothness of the global fields and convergence rates of the approximations for the wave equation. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Explicit time integration with lumped mass matrix for enriched finite elements solution of time domain wave problems

We present a partition of unity finite element method for wave propagation problems in the time domain using an explicit time integration scheme. Plane wave enrichment functions are introduced at the finite elements nodes which allows for a coarse mesh at low order polynomial shape functions even at high wavenumbers. The initial condition is formulated as a Galerkin approximation in the enriched function space. We also show the possibility of lumping the mass matrix which is approximated as a block diagonal system. The proposed method, with and without lumping, is validated using three test cases and compared to an implicit time integration approach. The stability of the proposed approach against different factors such as the choice of wavenumber for the enrichment functions, the spatial discretization, the distortions in mesh elements or the timestep size, is tested in the numerical studies. The method performance is measured for the solution accuracy and the CPU processing times. The results show significant advantages for the proposed lumping approach which outperforms other considered approaches in terms of stability. Furthermore, the resulting block diagonal system only requires a fraction of the CPU time needed to solve the full system associated with the non-lumped approaches.     

A weak Galerkin finite element method for the stokes equations

This paper introduces a weak Galerkin (WG) finite element method for the Stokes equations in the primal velocity-pressure formulation. This WG method is equipped with stable finite elements consisting of usual polynomials of degree k≥1 for the velocity and polynomials of degree k?1 for the pressure, both are discontinuous. The velocity element is enhanced by polynomials of degree k?1 on the interface of the finite element partition. All the finite element functions are discontinuous for which the usual gradient and divergence operators are implemented as distributions in properly-defined spaces. Optimal-order error estimates are established for the corresponding numerical approximation in various norms. It must be emphasized that the WG finite element method is designed on finite element partitions consisting of arbitrary shape of polygons or polyhedra which are shape regular.     

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.     

OPTIMAL ERROR ESTIMATES OF THE PARTITION OF UNITY METHOD WITH LOCAL POLYNOMIAL APPROXIMATION SPACES

In this paper, we provide a theoretical analysis of the partition of unity finite elementmethod (PUFEM), which belongs to the family of meshfree methods. The usual erroranalysis only shows the order of error estimate to the same as the local approximations[12].Using standard linear finite element base functions as partition of unity and polynomials aslocal approximation space, in 1-d case, we derive optimal order error estimates for PUFEMinterpolants. Our analysis show that the error estimate is of one order higher than thelocal approximations. The interpolation error estimates yield optimal error estimates forPUFEM solutions of elliptic boundary value problems.     

Finite Element Modelling of Cracks Based on the Partition of Unity Method

In this paper a finite element model for the analysis of brittle materials in the post cracking regime is presented. The model allows the representation of failure zones several times smaller than the structure itself using relatively coarse finite element meshes. The formulation is based on the partition of unity method. Discontinuous shape functions are used to enrich the continuous approximation of the displacement field where a crack has opened [2]. The magnitude of the displacement jump is determined by extra degrees of freedom at existing nodes. The crack path is completely independent of the structure of the mesh and is continuous across element boundaries. To model inelastic deformations around the crack tip a cohesive crack model is used. A representative numerical example illustrates the performance of the proposed model.     

Generalizations of the Finite Element Method

This paper is concerned with the generalization of the finite element method via the use of non-polynomial enrichment functions. Several methods employ this general approach, e.g. the extended finite element method and the generalized finite element method. We review these approaches and interpret them in the more general framework of the partition of unity method. Here we focus on fundamental construction principles, approximation properties and stability of the respective numerical method. To this end, we consider meshbased and meshfree generalizations of the finite element method and the use of smooth, discontinuous, singular and numerical enrichment functions.     

非嵌套有限元的W循环多重网格方法

本文讨论了下述情形：１非嵌套网格；２曲边有限元；３非协调元；４拟协调元；５有限元的型函数有特殊性质，都能导致非嵌套的有限元空间．对一个包括上述情形的问题给出了非嵌套有限元的Ｗ循环多重网格方法，并证明了它的收敛性。