非定常的热传导-对流问题的混合有限元-Galerkin交替方向有限元法

1 引 言设Ω　R2为具有光滑边界的有界区域,考虑非定常的,无量纲化的,而且带有热传导的粘性不可压缩流体力学问题:问题（Ⅰ）:求u=（u1,u2）,p,T满足:     

一种有限元-边界元耦合分域算法

提出了一种有限元-边界元耦合分域算法．该算法将所分析问题的区域分解成有限元和边界元子域,在满足两子域界面上位移和面力协调连续的条件下,通过迭代求解得到问题的解．在迭代求解过程中,引入动态松弛系数,使收敛得以加速．该方法在两子域界面上有限单元结点和边界单元结点的位置相互独立,无需协调一致,对诸如裂纹扩展过程的模拟具有独特的优势．用所提出的耦合算法分析算例,得到的结果与有限元法、边界元法和另一种耦合算法的数值计算结果一致,验证了这种算法的正确性和可行性．     

三维热传导型半导体问题的交替方向特征有限元方法及理论分析

１．引言 三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述 ［１,２］．工程研究中一般考虑绝流边条件,由于绝流条件可以看作一反射条件来处理、为了数值分析方便,我们在此考虑三维周期问题： 其中,　＝［０,１］３,未知函数是电子位势　;电子,空穴浓度ｅ,ｐ;温度函数Ｔ．方程（１,１）－（１.４）中出现的系数均有正的上下界,且是　周期的． ａ＝Ｑ／ε,Ｑ,ε分别表示电子负荷和介电系数,均为正常数．Ｎ（ｘ）是给定的函数．Ｄｓ（ｘ）为扩散系数,μｓ（ｘ）为迁移率,ｓ＝ｅ,Ｐ．Ｒ（ｅ,ｐ,Ｔ）…     

裂缝-孔隙介质中地下水污染问题的Galerkin交替方向有限元方法

考虑裂缝 孔隙介质中地下水污染问题均匀化模型的数值模拟．对压力方程采用混合元方法,对浓度方程采用Ｇａｌｅｒｋｉｎ交替方向有限元方法,对吸附浓度方程采用标准Ｇａｌｅｒｋｉｎ方法,证明了交替方向有限元格式具有最优犔２ 和犎１ 模误差估计．     

Multiscale methods for elliptic homogenization problems

In this article we study two families of multiscale methods for numerically solving elliptic homogenization problems. The recently developed multiscale finite element method [Hou and Wu, J Comp Phys 134 (1997), 169–189] captures the effect of microscales on macroscales through modification of finite element basis functions. Here we reformulate this method that captures the same effect through modification of bilinear forms in the finite element formulation. This new formulation is a general approach that can handle a large variety of differential problems and numerical methods. It can be easily extended to nonlinear problems and mixed finite element methods, for example. The latter extension is carried out in this article. The recently introduced heterogeneous multiscale method [Engquist and Engquist, Comm Math Sci 1 (2003), 87–132] is designed for efficient numerical solution of problems with multiscales and multiphysics. In the second part of this article, we study this method in mixed form (we call it the mixed heterogeneous multiscale method). We present a detailed analysis for stability and convergence of this new method. Estimates are obtained for the error between the homogenized and numerical multiscale solutions. Strategies for retrieving the microstructural information from the numerical solution are provided and analyzed. Relationship between the multiscale finite element and heterogeneous multiscale methods is discussed. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

拟协调元的位移函数及节点误差

直接从拟协调元的应变关系式出发,构造具有明确物理意义的幂级数形式的位移函数,从而得出拟协调元的常应变和线性应变系数是唯一确定的,它只能收敛到常应变的结论;刚性位移项可采用多种构造方法,不同的方法得出的节点参数与单元的本身的节点参数存在不同阶次的误差,这与常规位移法有限元不同。     

一种新的声椭球无限单元

提出了一种新的声椭球无限单元．这种声无限单元基于一种新的声压表达式,这种声压表达式能够更准确地代表着椭球声场的声传播模式．这种新方法的形函数类似于Burnett方法,而权函数定义为形函数和一个附加因子的乘积．因为仅需要一维的数值积分,这种新方法的代码生成十分容易,就像处理一维单元一样．耦合标准的有限元程序,这种声无限单元理论上能够高效地求解任何形状的声源的声辐射和声散射现象．简要地推导了这种新方法,并给出了这种方法详尽的推导结果．为更有效地检验该无限元方法的可行性,文中例子仅考虑无限元求解的精度,而不包括相应的有限元．使用这种新方法,精确地推导出了摆动球的理论计算公式．而长旋转椭球的例子则表明了这种方法优于边界元方法和其他声椭球无限元方法．这些例子表明了这种新方法是切实可行的．     

特征值问题的自适应反迭代有限元算法

1引言设Ω∈R~2为Lipschitz单连通的有界闭区域,X为定义在Ω的Sobolev空间,a(·,·)和b(·,·)为X×X→C的有界双线性或半双线性泛函,考虑变分特征值问题:求(λ,u≠0)∈C×X使得a(u,v)=λb(u,u),(?)u∈X,其中a(·,·)满足X上的"V-强制性"条件或者连续的inf-sup条件,设M_h为Q区域上的正则三角形剖分,X_h∈X为定义在M_h有限元子空间,上述变分问题对应的有限元离散问题为:求(λ_h,u_h)∈R×X,u_h≠0使得     

A Characteristic Mixed Finite Element Two-Grid Method for Compressible Miscible Displacement Problem

A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium. The concentration equation is treated by a mixed finite element method with characteristics (CMFEM) and the pressure equation is treated by a parabolic mixed finite element method (PMFEM). Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations. Moreover, the $L^q$ error estimates are conducted for the pressure, Darcy velocity and concentration variables in the two-grid solutions. Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.     

Development and Analysis of Higher Order Finite Volume Methods over Rectangles for Elliptic Equations

Currently used finite volume methods are essentially low order methods. In this paper, we present a systematic way to derive higher order finite volume schemes from higher order mixed finite element methods. Mostly for convenience but sometimes from necessity, our procedure starts from the hybridization of the mixed method. It then approximates the inner product of vector functions by an appropriate, critical quadrature rule; this allows the elimination of the flux and Lagrange multiplier parameters so as to obtain equations in the scalar variable, which will define the finite volume method. Following this derivation with different mixed finite element spaces leads to a variety of finite volume schemes. In particular, we restrict ourselves to finite volume methods posed over rectangular partitions and begin by studying an efficient second-order finite volume method based on the Brezzi–Douglas–Fortin–Marini space of index two. Then, we present a general global analysis of the difference between the solution of the underlying mixed finite element method and its related finite volume method. Then, we derive finite volume methods of all orders from the Raviart–Thomas two-dimensional rectangular elements; we also find finite volume methods to associate with BDFM ₂ three-dimensional rectangles. In each case, we obtain optimal error estimates for both the scalar variable and the recovered flux.     

弱有限元方法在线弹性问题中的应用

本文考虑弱有限元（简称WG）方法在线弹性问题中的应用.WG方法是传统有限元方法的推广,用于偏微分方程的数值求解.和传统有限元一样,它的基本思想源于变分原理.WG方法的特点是使用在剖分单元内部和剖分单元边界上分别有定义的分片多项式函数（即弱函数）作为近似函数来逼近真解,并针对弱函数定义相应的弱微分算子代入数值格式进行计算.除此之外,WG方法允许在数值格式中引进稳定子以实现近似函数的弱连续性.WG方法具有允许使用任意多边形或多面体剖分,数值格式与逼近函数构造简单,易于满足相应的稳定性条件等优点.本文考虑WG方法在求解线弹性问题中的应用.围绕线弹性问题数值求解中常见的三个问题,即：数值格式的强制性,闭锁性,应力张量的对称性介绍WG方法在线弹性问题求解中的应用.     

Analysis of a non-standard mixed finite element method with applications to superconvergence

We show that a non-standard mixed finite element method proposed by Barrios and Gatica in 2007, is a higher order perturbation of the least-squares mixed finite element method. Therefore, it is also superconvergent whenever the least-squares mixed finite element method is superconvergent. Superconvergence of the latter was earlier investigated by Brandts, Chen and Yang between 2004 and 2006. Since the new method leads to a non-symmetric system matrix, its application seems however more expensive than applying the least-squares mixed finite element method. Dedicated to Ivan Hlaváček on the occasion of his 75th birthday     

二维非均匀水沙模型的多步特征-混合元数值模拟

通过对方程的对流部分采用沿着特征线方向向后两步差分格式进行离散,而对扩散部分采用混合有限元格式进行离散,从而利用多步特征-混合有限元方法对平面非均匀水沙模型进行了数值模拟,给出了相应的误差分析及数值算例.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

A space‐time mixed‐hybrid finite element method for the damped wave equation

A space‐time finite element method is introduced to solve the linear damped wave equation. The scheme is constructed in the framework of the mixed‐hybrid finite element methods, and where an original conforming approximation of H(div;Ω) is used, the latter permits us to obtain an upwind scheme in time. We establish the link between the nonstandard finite difference scheme recently introduced by Mickens and Jordan and the scheme proposed. In this regard, two approaches are considered and in particular we employ a formulation allowing the solution to be marched in time, i.e., one only needs to consider one time increment at a time. Numerical results are presented and compared with the analytical solution illustrating good performance of the present method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

A New Reduced-Order fve Algorithm Based on POD Method for Viscoelastic Equations

A proper orthogonal decomposition (POD) technique is used to reduce the finite volume element (FVE) method for two-dimensional (2D) viscoelastic equations. A reduced-order fully discrete FVE algorithm with fewer degrees of freedom and sufficiently high accuracy based on POD method is established. The error estimates of the reduced-order fully discrete FVE solutions and the implementation for solving the reduced-order fully discrete FVE algorithm are provided. Some numerical examples are used to illustrate that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order fully discrete FVE algorithm is one of the most effective numerical methods by comparing with corresponding numerical results of finite element formulation and finite difference scheme and that the reduced-order fully discrete FVE algorithm based on POD method is feasible and efficient for solving 2D viscoelastic equations.     

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

分析了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²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

Each averaging technique yields reliable a posteriori error control in FEM on unstructured grids. Part I: Low order conforming, nonconforming, and mixed FEM

Averaging techniques are popular tools in adaptive finite element methods for the numerical treatment of second order partial differential equations since they provide efficient a posteriori error estimates by a simple postprocessing. In this paper, their reliablility is shown for conforming, nonconforming, and mixed low order finite element methods in a model situation: the Laplace equation with mixed boundary conditions. Emphasis is on possibly unstructured grids, nonsmoothness of exact solutions, and a wide class of averaging techniques. Theoretical and numerical evidence supports that the reliability is up to the smoothness of given right-hand sides.

     

An algebraic multigrid method for finite element systems on criss-cross grids

In this paper, we design and analyze an algebraic multigrid method for a condensed finite element system on criss-cross grids and then provide a convergence analysis. Criss-cross grid finite element systems represent a large class of finite element systems that can be reduced to a smaller system by first eliminating certain degrees of freedoms. The algebraic multigrid method that we construct is analogous to many other algebraic multigrid methods for more complicated problems such as unstructured grids, but, because of the specialty of our problem, we are able to provide a rigorous convergence analysis to our algebraic multigrid method. Dedicated to Professor Charles A. Micchelli on the occasion of his 60th birthday The work was supported in part by NSAF(10376031) and National Major Key Project for basic researches and by National High-Tech ICF Committee in China.