三维光滑区域上Navier-Stokes方程的有限元方法

§1.引言设Ω是R~3中的有界区域,且属于C~2.v是正常数.已知:当f∈(L~2(Ω))~3.且||f||0充分小时,Navier—Stokes方程的解存在且唯一.另外,u∈(H~2(Ω)∩H_1~0(Ω))~3,P∈H~1(Ω)\R. 最近,Bernardi讨论三维多面体区域Ω上Stokes方程的有限元解法.有限元空间由分片多项式(关于u为分片特殊三次多项式,关于p为分片常数)构成,进行误差估计时.要求u∈(H~2(Ω)∩H_0~1(Ω))~3,P∈H~1(Ω)\R.当Ω为二维区域上的凸多角形时,Stokes     

一类带弱奇异核非线性偏积分微分方程的全离散有限元

1引言我们将研究下面一类带弱奇异核非线性偏积分微分方程的数值解:u_t-▽·(a(u)▽u)-integral from n=0 to tβ(t-s)△u(s)ds=f(u),x∈Ω,t∈(?),(1.1) u(·,t)=0,x∈(?)Ω,t∈J,(1.2) u(·,0)=v(x),x∈Ω,(1.3)其中Ω为平面上的凸角域,J=(0,T],α和f为R上的光滑函数,满足0相似文献   

关于不协调有限元空间最大模的不等式

§1.引言 对于二维协调元有限元空间的元素v_h,已有下述估计 ||v_h||_(0,∞,Ω)≤C|lnh|~(1/2)|v_h||(1,Z,Ω) (1)它为有限元解的L~∞收敛性及超收敛性研究提供了工具.本文试图把上述估计推广到一类包括非协调元、杂交元和拟协调元空间的有限元空间.     

关于特征值问题有限元近似解的最大模估计及其校正加速估计

本文记通常CooeB空间W_p~m(Ω)(P=2时记为H~m(Ω))的模为||·||_(m,p),m=0,即空间L_p(Ω)的模简记为||·||_p。 设Ω为平面有界区域,边界?Ω分段光滑,则(1)的解存在,且特征值可列:     

位移障碍下一个四阶变分不等式的某些强间断非协调元逼近

考虑[1]中四阶变分不等式问题:其中为非空闭凸集,而障碍函数φ∈C~2(Ω),φ<0,在?Ω上.关于解的性质,有下述结果:当Ω?R~2是具有光滑边界?Ω的有界凸区域且f∈L~2(Ω)时,问题(1)存在唯     

平面负位势相关的Schr(o)dinger方程的非负广义解

设β是复平面上圆盘Ωa={z ||z|＜a}内的一个零容紧致集.考虑Ωβα=Ωα\β上的定常Schrodinger方程(-A+μ)u=0,其中位势μ≤0是Kato类Radon测度.将方程在广义函数意义下的在{z||z|=a}上取极限值0的非负连续解族记为μH+.对Ωβα的Kerekjato-Stoilow意义下的理想边界β的任一点ζ,本文通过定义μH+→μH+的线性算子πζ,引入Martin函数Kζ,证明了μH+=Hβ Pβ,其中Hβ={u∈μH+|πζ(u)=0,vζβ},Pβ={u∈μH+|u=∞∑i=i ciKζi,ζi∈β,ci≥0}.     

平面负位势相关的Schrödinger方程的非负广义解

设β是复平面上圆盘Ωa={z ||z|＜a}内的一个零容紧致集.考虑Ωβα=Ωα\β上的定常Schrodinger方程(-A+μ)u=0,其中位势μ≤0是Kato类Radon测度.将方程在广义函数意义下的在{z||z|=a}上取极限值0的非负连续解族记为μH+.对Ωβα的Kerekjato-Stoilow意义下的理想边界β的任一点ζ,本文通过定义μH+→μH+的线性算子πζ,引入Martin函数Kζ,证明了μH+=Hβ Pβ,其中Hβ={u∈μH+|πζ(u)=0,vζβ},Pβ={u∈μH+|u=∞∑i=i ciKζi,ζi∈β,ci≥0}.     

关于不协调有限元空间最大模的不等式

§1.引言 对于二维协调元有限元空间的元素v_h,已有下述估计 ||v_h||_(0,∞,Ω)≤C|lnh|~(1/2)|v_h||(1,Z,Ω) (1)它为有限元解的L~∞收敛性及超收敛性研究提供了工具.本文试图把上述估计推广到一类包括非协调元、杂交元和拟协调元空间的有限元空间.     

关于“分片线性有限元最优L∞—误差估计”的注

关于二阶问题的分片线性有限元近似解的L~∞—误差估计,已有不少出色的工作,Nitsche,Scott,Frehse和Rannacher,Schatz和Wahlbin等分别利用不同的方法,证明了在u∈W~(2,∞)(Ω)的条件下,有‖u-u_h‖_(0,∞,Ω=0(h~2|lnh|)。这一误差估计与分片线性插值的逼近误差阶相比,并不是丰满(或最优)的。能否得到与插值逼近相同阶的L~∞—误差     

三维带有衰减项的不可压缩磁流体力学方程组弱解与强解的研究

论文研究了带有衰减项的磁流体力学方程组的柯西问题.当β≥1及初值u_0,b_0∈L~2(R~3)时,采用Galerkin方法证明了方程组存在全局弱解.并且当初值u_0∈H_0~1∩L~(β+1)(R~3),b_0∈H_0~1(R~3)时,可以得到方程组存在唯一局部强解.     

A coupling of nonconforming and mixed finite element methods for Biot's consolidation model

In this article, we develop a nonconforming mixed finite element method to solve Biot's consolidation model. In particular, this work has been motivated to overcome nonphysical oscillations in the pressure variable, which is known as locking in poroelasticity. The method is based on a coupling of a nonconforming finite element method for the displacement of the solid phase with a standard mixed finite element method for the pressure and velocity of the fluid phase. The discrete Korn's inequality has been achieved by adding a jump term to the discrete variational formulation. We prove a rigorous proof of a‐priori error estimates for both semidiscrete and fully‐discrete schemes. Optimal error estimates have been derived. In particular, optimality in the pressure, measured in different norms, has been proved for both cases when the constrained specific storage coefficient c₀ is strictly positive and when c₀ is nonnegative. Numerical results illustrate the accuracy of the method and also show the effectiveness of the method to overcome the nonphysical pressure oscillations. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

受限制多项式插值及在构造形函数空间中的应用

1引言 G.Strang指出：有限元法的新思想在于试控函数的选择，目标是选择这样的分片多项式，它们被少数几个方面的节点值确定，并仍具有我们需要的连续性和逼近度。受限制多项式插值空间就是这样一类空间，在P.G.Ciarlet~[4]的书中有较多的介绍，采用的方法是通过约束条件来决定试验空间，但正如[1]中指出的，这样约束条件欠直观，且容易产生一些不确     

A SUPERCONVERGENT ESTIMATE OF WILSON-LIKE ELEMENTS

In this paper, the supercouvergence of a class of Wilson-like elements is considered, and a superconvergent estimate of Wilson-like elements is obtained for strong regular meshes.     

轴对称问题的Wilson四边形单元

给出了一族轴对称问题的改进Ｗｉｌｓｏｎ元，利用强分片检验证明了其收敛性，讨论了单元函数结构·从而给出了一种构造收敛的轴对称非协调元方法·     

Anisotropic interpolation and quasi-Wilson element for narrow quadrilateral meshes

In this paper an anisotropic interpolation theorem is presentedthat can be easily used to check the anisotropy of an element.A kind of quasi-Wilson element is considered for second-orderproblems on narrow quadrilateral meshes for which the usualregularity condition _K/h_K c₀ > 0 is not satisfied, whereh_K is the diameter of the element K and _K is the radius of thelargest inscribed circle in K. Anisotropic error estimates ofthe interpolation error and the consistency error in the energynorm and the L²-norm are given. Furthermore, we give a Poincaréinequality on a trapezoid which improves a result of eniek.     

A new rotated nonconforming pyramid element

In this paper, a new nonparametric nonconforming pyramid finite element is introduced. This element takes the five face mean values as the degrees of the freedom and the finite element space is a subspace of P₂. Different from the other nonparametric elements, the basis functions of this new element can be expressed explicitly without solving linear systems locally, which can be achieved by introducing a new reference pyramid. To evaluate the integration, a class of new quadrature formulae with only two/three equally weighted points on pyramid are constructed. We present the error estimation in the presence of quadrature formulae. Numerical results are shown to confirm the optimality of the convergence order for the second order elliptic problems.     

Nonconforming mixed finite elements for linear elasticity on simplicial grids

This paper introduces a new family of nonconforming mixed finite elements for solving the linear elasticity equations on simplicial grids. Besides, this paper describes the construction of the lowest order basis functions. The construction only involves simple computations due to the new explicit stress shape function spaces and the procedure applies for high order cases. Numerical experiments for four benchmark problems in mechanics indicate the robust locking‐free behavior and show that the lowest order nonconforming mixed method leads to smaller stress errors than the first and second order standard Galerkin methods for the nearly incompressible case.     

一类非协调元的收敛性分析

１．引言最近,［１］在求解Ｓｔｏｋｅｓ方程时提出了一类非协调元．这类非协调元定义在矩形网格上,形式简单,自由度少,比熟知的多线性元还少．例如,对于三维问题,它们只有六个自由度,是型函数在各个表面中点处的函数值或者表面上的平均值．与一些熟知的非协调元不同之处是,由于这类非协调元简单的形式及其非常少的自由度,它们不包含任何协调的部分,但是,已有的结果表明它们具有良好的数值表现［１］,［２］中已有一些理论分析,用它们来求解晶体的微结构方程,也取得较好的结果．本文将详细分析这类非协调元的收敛性,所用的工…     

A locally divergence-free nonconforming finite element method for the time-harmonic Maxwell equations

A new numerical method for computing the divergence-free part of the solution of the time-harmonic Maxwell equations is studied in this paper. It is based on a discretization that uses the locally divergence-free Crouzeix-Raviart nonconforming vector fields and includes a consistency term involving the jumps of the vector fields across element boundaries. Optimal convergence rates (up to an arbitrary positive ) in both the energy norm and the norm are established on graded meshes. The theoretical results are confirmed by numerical experiments.

     

A new quadratic nonconforming finite element on rectangles

A new quadratic nonconforming finite element on rectangles (or parallelograms) is introduced. The nonconforming element consists of P₂ ⊕ Span{x²y,xy²} on a rectangle and eight degrees of freedom. Our element is essentially of seven degrees of freedom since the degree of freedom associated with the integration on rectangle is essentially of bubble‐function nature. Global basis functions are constructed for both Dirichlet and Neumann type of problems; accordingly the corresponding dimensions are counted. The local and global interpolation operators are defined. Error estimates of optimal order are derived in both broken energy and L²(Ω) norms for second‐order of elliptic problems. Brief numerical results are also shown to confirm the optimality of the presented quadratic nonconforming element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006