伪双曲型方程的一个H1-Galerkin非协调混合元格式

讨论了一类伪双曲型方程的一个H1-Galerkin非协调混合有限元方法.利用插值算子的特殊性质,在半离散和全离散格式下,得到了与传统混合有限元相同的误差估计且不需要满足LBB条件.     

一类非线性双曲型方程的非协调有限元方法

研究了一类非线性双曲型方程的非协调有限元方法,在不需要传统的Ritz投影的情况下,得到了半离散格式下的误差估计及超收敛结果.     

特征值问题的Lagrange型各向异性有限元方法

在各向异性网格下首先研究了二阶椭圆特征值问题算子谱逼近的若干抽象结果.然后将这些结果具体应用于线性和双线性Lagrange型协调有限元,得到了与传统有限元网格剖分下相同的最优误差估计,从而拓宽了已有的成果.     

双曲型方程的非协调变网格有限元方法

采用变网格的思想讨论了双曲型方程在各向异性网格下的Crouzeix-Raviart型非协调有限元逼近.在不需要引入传统分析中Riesz投影的情况下,得到了相应最优误差估计.     

带弱奇异核的抛物型积分微分方程的非协调有限元方法

研究了带弱奇异核的抛物型积分微分方程的非协调有限元方法,在不需要Ritz-Volterra投影的情况下,在半离散和全离散的格式下分别得到了与协调有限元方法相同的误差估计.     

抛物型方程基于POD方法的有限元格式

将特征正交分解(Proper Orthogonal Decomposition, 简记为POD)方法应用于抛物型方程通常的有限元格式, 简化其为一个计算量很少但 具有足够高精度的POD有限元格式, 并给出简化的POD有限元解的误差分析. 数值例子表明在简化的POD有限元解和通常的有限元解之间的误差足够小的情形下, POD有限元格式比通常的有限元格式大大地节省计算量, 从而验证POD方法的有效性.     

重叠型非匹配网格抛物型方程的有限元法及收敛性分析

基于单位分解技术,本文介绍了重叠型非匹配网格抛物型方程初边值问题的有限元方法,分别给出了半离散有限元、全离散有限元格式和收敛性分析的结果,并给出了数值例子.     

抛物型方程一个新的非协调混合元超收敛性分析及外推

本文研究了抛物型方程在新混合元格式下的非协调混合有限元方法. 在抛弃传统有限元分析的必要工具-Ritz 投影算子的前提下,直接利用单元的插值性质,运用高精度分析和对时间t的导数转移技巧,借助于插值后处理技术,分别导出了关于原始变量u的H¹-模和通量p=▽u在L²-模下的O（h²）阶超逼近性质和整体超收敛. 进一步,通过构造合适的辅助问题,运用Richardson 外推格式,得到了具有更高精度O（h³）阶的外推结果. 最后,给出了一些数值结果验证了理论分析的正确性.     

几乎不可压线弹性问题的新的Uzawa型自适应有限元方法

本文针对几乎不可压线弹性问题设计新的Uzawa型自适应有限元方法,该方法可克服"闭锁"现象·通过引入"压力"变量将弹性问题转化为一个鞍点系统,对该系统将Uzawa型迭代法和自适应有限元方法相结合,建立了Uzawa型自适应有限元方法,并给出了该算法的收敛性.该算法采用低阶协调有限元通近空间变量,选取的有限元空间对无需满足离散的BB条件.最后,数值算例验证了理论结果的正确性.     

非线性抛物型方程的全离散有限元方法

§1.引言关于非线性抛物型方程有限元方法的研究已有许多工作.但所讨论的方程关于梯度是线性的,仅得到全离散有限元逼近的L_2误差估计及较弱意义下(两层平均)的H~1误差估计.本文讨论关于梯度亦是非线性的抛物型方程.得到了最佳L_2、H~1(较强意义下)、L.及时间导数的误差估计. 考虑下述抛物型方程的混合问题:     

一种高精度的裂纹奇异单元

基于广义伽辽金法的多变量有限元算法,增加了连续体力学有限元模型建立的灵活性.本文利用它,通过数值试验的对比建立了一种高精度的含奇异性的裂纹单元,并对多变量奇异元的构成进行了探讨.     

Simple finite element method in vorticity formulation for incompressible flows

A very simple and efficient finite element method is introduced for two and three dimensional viscous incompressible flows using the vorticity formulation. This method relies on recasting the traditional finite element method in the spirit of the high order accurate finite difference methods introduced by the authors in another work. Optimal accuracy of arbitrary order can be achieved using standard finite element or spectral elements. The method is convectively stable and is particularly suited for moderate to high Reynolds number flows.

     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

Stokes问题Q_2-P_1混合元外推方法

考虑Stokes问题的有限元解与精确解插值的Q2-P1混合元的渐近误差展开和分裂外推.首先利用积分恒等式技巧确定了微分方程精确解与有限元插值之间积分式的主项,其次再借助插值后处理和分裂外推技术,得到了比通常的误差估计提高两阶的收敛速度.     

A wavelet-based stochastic finite element method of thin plate bending

A wavelet-based stochastic finite element method is presented for the bending analysis of thin plates. The wavelet scaling functions of spline wavelets are selected to construct the displacement interpolation functions of a rectangular thin plate element and the displacement shape functions are expressed by the spline wavelets. A new wavelet-based finite element formulation of thin plate bending is developed by using the virtual work principle. A wavelet-based stochastic finite element method that combines the proposed wavelet-based finite element method with Monte Carlo method is further formulated. With the aid of the wavelet-based stochastic finite element method, the present paper can deal with the problem of thin plate response variability resulting from the spatial variability of the material properties when it is subjected to static loads of uncertain nature. Numerical examples of thin plate bending have demonstrated that the proposed wavelet-based stochastic finite element method can achieve a high numerical accuracy and converges fast.     

A posteriori error analysis for linearization of nonlinear elliptic problems and their discretizations

The paper is devoted to a posteriori quantitative analysis for errors caused by linearization of non-linear elliptic boundary value problems and their finite element realizations. We employ duality theory in convex analysis to derive computable bounds on the difference between the solution of a non-linear problem and the solution of the linearized problem, by using the solution of the linearized problem only. We also derive computable bounds on differences between finite element solutions of the nonlinear problem and finite element solutions of the linearized problem, by using finite element solutions of the linearized problem only. Numerical experiments show that our a posteriori error bounds are efficient.     

Conforming and nonconforming harmonic finite elements

ABSTRACT

Instead of using the full polynomial space, a conforming and a nonconforming finite element methods are designed where only harmonic polynomials (a much smaller space) are employed in the computation. The conforming quadratic harmonic polynomial finite element is defined only on a special triangular grid. The nonconforming quadratic harmonic finite element is defined on general triangular grids. The optimal order of convergence is proved for both finite element methods, and confirmed by numerical computations. In addition, numerical comparisons with the standard conforming and nonconforming finite elements are presented.     

An operator splitting method with FDM and FEM for the convection- diffusion equation using boundary-fitted coordinate system

An operator splitting method combining finite difference method and finite element method is proposed in this paper by using boundary-fitted coordinate system. The governing equation is split into advection and diffusion equations and solved by finite difference method using boundary-fitted coordinate system and finite element method respectively. An example for which analytic solution is available is used to verified the proposed methods and the agreement is very good. Numerical results show that it is very efficient.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Stability analysis of frame structures by bubble-function method

In the stability analysis of frame structures, the results by conventional finite element method (FEM) in which one member is taken as one element are sometimes unavailable. This paper took a new basic function system with bubble functions as the shape function of a bar element to develop a bubble function finite element method (BFEM), in which the bending and the geometric stiffness matrices were derived from the principle of virtual work. Bubble functions are finite element modes that are located entirely within a single element and are zero on boundaries of the element, but are nonzero at the other points. BFEM is as concise as conventional bar FEM but has better accuracy, and is adaptable to the buckling analysis of all kinds of frame structures. The use of bubble functions significantly improves the convergence of finite element analysis, and efficiently reduces the computation cost for the buckling analysis of frame structures. Numerical results show that using bubble functions in finite element for the stability analysis of structures is very efficient, especially for high-rise and large-scale frame structures.