粘性不可压缩对流占优Oseen方程的非协调子格稳定化有限元法分析

<正>1引言本文考虑粘性不可压缩对流占优Oseen方程,(?)(1)其中:Ω(?)R~d(d=2)为具有Lipschitz连续边界的有界开集,β∈W~(1,∞)(Ω)且▽·β=0,μ、σ为常数,f∈L~2(Ω).当采用通常的混合有限元方法(MFEM)求解时,一般会遇到以下两个困难:·为保证速度和压力数值解稳定,要求有限元空间满足inf-sup(or Babuska-Brezzi)条件.·当对流占优,即0μ《||β||_(L~∞(Ω))时,数值解会产生伪振荡.     

对流占优的Sobolev方程的投影稳定化有限元方法

对于对流占优的Sobolev方程,提出了一种新的投影稳定化有限元方法,建立了半离散和全离散的投影稳定化格式,给出了解的稳定性和收敛性分析.该方法能够有效克服对流占优,与内罚方法相比,投影格式更简单,计算量更小,且得到的C—N格式是无条件稳定的,时间精度达到了二阶.最后,通过实验证明,数值结果与理论结果完全一致.     

非定常对流扩散问题的非协调局部投影有限元方法

本文将近年来基于协调有限元逼近提出的涡旋粘性法推广应用到非协调有限元逼近,对非定常的对流占优扩散问题,空间采用非协调Crouzeix-Raviart元逼近,时间用Crank-Nicolson差分离散格式,提出了Crank-Nicolson差分-局部投影法稳定化有限元格式,我们对稳定性和误差估计给出了详细的分析,得出了最...     

对流扩散方程的稳定化流线扩散法最小二乘非协调有限元分析

将最小二乘法和稳定化的流线扩散法相结合,研究了对流扩散方程的非协调有限元格式,用矩形EQ_1~(rot)元和零阶R-T元分别来逼近位移和应力,利用单元本身的特殊性质,证明了离散格式解的存在惟一性,得到了位移H~1-模和应力H(div)-模的最优误差估计.     

对流占优扩散方程的特征有限体积元方法

考虑对流占优扩散方程初边值问题的特征有限体积元方法,并给出特征有限体积元解的误差分析.理论分析表明特征有限体积元解具有最优阶L~2和H~1模误差估计.数值算例说明此方法是有效的.     

一类对流占优扩散方程的MMOCAAFVEM方法及分析

本文用具有调整对流的特征线修正方法(MMOCAA)与有限体积元方法相结合,构造出一种新的守恒型计算格式-MMOCAAFVEM,这种方法综合了特征线方法和有限体积元方法的主要优点.通过对对流项进行调整,以很小的额外计算量获取了问题的质量守恒性质,并且证明该方法具有最优阶H1误差估计.     

线性对流占优扩散方程的后验误差估计

对于线性对流占优扩散方程,采用特征线有限元方法离散时间导数项和对流项,用分片线性有限元离散空间扩散项,并给出了一致的后验误差估计,其中估计常数不依赖与扩散项系数。     

非定常线性化Navier-Stokes方程的子格粘性非协调有限元方法

本文结合子格粘性法的思想,空间采用非协调Crouzeix-Raviart元逼近,时间采用Crank-Nicolson差分离散,对非定常线性化Navier-Stokes方程建立了全离散的子格粘性非协调有限元格式.对稳定性和误差估计作出了详细的分析, 得出了最优的误差估计.最后, 通过数值算例进一步验证了该方法的稳定性和收敛性.     

线性对流占优Sobolev方程的流线扩散有限元方法

Abastract. In this paper,a streamline-diffusion F. E. M. for linear Sobolev equations with con-vection-dominated term is given. According to the range of space-time F. E mesh parameter h,two choices for artifical diffusion parameter are presented,and for the corresponding computa-tion schemes the stability and error estimates in suitable norms are estabilished.     

对流占优扩散方程的分裂特征混合有限元方法

利用修正的特征线方法,构建一类求解对流占优扩散方程的分裂特征混合有限元算法.在新的算法中,混合系统的系数矩阵对称正定,且原未知函数u与流函数σ=-ε▽u可分离求解.推导了加权能量模意义下的最优阶误差估计,并给出数值算例验证理论上的分析结果.     

Residual-based a posteriori error estimates of nonconforming finite element method for elliptic problems with Dirac delta source terms

Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.     

基于Riesz-表示算子非协调稳定化有限元法

Following the framework of the finite element methods based on Riesz-representing operators developed by Duan Huoyuan in 1997,through discrete Riesz representing-operators on some virtual(non-) conforming finite-dimensional subspaces,a stabilization formulation is presented for the Stokes problem by employing nonconforming elements.This formulation is uniformly coercive and not subject to the Babu ka-Brezzi condition,and the resulted linear algebraic system is positive definite with the spectral condition number O(h-2). Quasi-optimal error bounds are obtained,which is consistent with the interpolation properties of the finite elements used.     

二维发展型对流占优扩散方程的FD-SD法的后验误差估计

引言 对流占优扩散问题是流体力学中一个典型的模型问题,对其数值求解始终是众多学者相当关心的课题．［１１］中指出,即使对于线性问题,通常其解在外流边界附近也会产生剧烈变化．倘若在内流边界上所给出的边值函数存在不连续点时,则在沿过此不连续点的特征线（流线）附近会出现断层．因此在数值求解对流占优扩散问题时,尽管标准有限元法具有高阶精度,但常产生数值剧烈振荡Ｓ而古典人工粘性Ｇａｌｅｒｋｉｎ法虽具有较好的稳定性,但仅具有一阶精度．流线扩散法（Ｓｔｒｅａｍｌｉｎｅ　 Ｄｉｆｆｕｓｉｏｎ Ｍｅｔｈｏｄ,简称 ＳＤ…     

各向异性网格下具有数值积分的非协调有限元逼近

In this paper we mainly discuss the nonconforming finite element method for second order elliptic boundary value problems on anisotropic meshes.By changing the discretization form(i.e.,by use of numerical quadrature in the procedure of computing the left load),we obtain the optimal estimate O(h),which is as same as in the traditional finite element analysis when the load f∈H~1(Ω)∩C~0(Ω)which is weaker than the previous studies.The results obtained in this paper are also valid to the conforming triangular element and nonconforming Carey's element.     

A multilevel preconditioner for the C-R FEM for elliptic problems with discontinuous coefficients

In this paper, we propose a multilevel preconditioner for the Crouzeix-Raviart finite element approximation of second-order elliptic partial differential equations with discontinuous coefficients. Since the finite element spaces are nonnested, weighted intergrid transfer operators, which are stable under the weighted L2 norm, are introduced to exchange information between different meshes. By analyzing the eigenvalue distribution of the preconditioned system, we prove that except a few small eigenvalues, all the other eigenvalues are bounded below and above nearly uniformly with respect to the jump and the mesh size. As a result, we get that the convergence rate of the preconditioned conjugate gradient method is quasi-uniform with respect to the jump and the mesh size. Numerical experiments are presented to confirm our theoretical analysis.