基于当地笛卡尔架构的无网格方法

提出了一种新的无网格方法,该方法是自动地在每一样点建立一个局部笛卡尔架构并选取相应的邻近点,然后运用全导数公式构造该样点的所有导数,它不需要任何网格单元,所以是彻底的无网格方法．数值算例表明,该方法具有很高的精度．     

解两点边值问题的一类修改的三次有限体积元法

构造了求解两点边值问题的一类修改的Lagrange型三次有限体积元法.试探函数空间取以四次Lobatto多项式的零点作为插值节点的Lagrange型三次有限元空间.将插值多项式的导数超收敛点(应力佳点)作为对偶单元的节点,检验函数空间取相应于对偶剖分的分片常数函数空间.证明了新方法具有最优的H1模和L2模收敛阶,讨论了在应力佳点导数的超收敛性,并通过数值实验验证了理论分析结果.     

一个值得关注的高考新考点——导数

导数的综合应用是多方面的.如求曲线在某点处切线的斜率,判断函数的单调性,求单调区间以及求函数的极值与最值等.而且导数知识可直接跟函数、数列、不等式、向量、解几、立几等重要知识块产生密切联系,表现得非常活跃.现在高考命题十分强调“能力立意”,注     

解Hamilton-Jacobi方程的MUSCL格式

本文利用单调数值通量和分片线性重构导数的方法构造了一种求HJ方程数值解的有限差分格式:MUSCL格式,并证明该格式具有TVB稳定性.数值实验表明该格式具有二阶精度,能避免产生伪振荡,尤其在类似"角点"的间断处有较好的分辩率.     

投影型插值算子的超收敛性质及其应用

本文首先将证明矩形剖分单元上的Ｌｏｂａｔｔｏ点,Ｇａｕｓｓ点和拟Ｌｏｂａｔｔｏ点分别是二维投影型插值算子函数,梯度和二阶导数的逼近佳点;然后考虑了二阶椭圆边值问题的有限元近似．通过建立投影型插值算子各种形式的超收敛基本估计,证明了投影型插值算子的各类...     

抛物型积分-微分方程有限元近似的超收敛性质

1　引　　言有限元超收敛性质在有限元方法的研究中占有重要的地位 .利用超收敛性不仅可提高有限元实际计算的精度 ,而且还可得到后验误差估计 .对于椭圆问题有限元超收敛性质的研究目前已有了较丰富的结果 [1 - 3] ,而对于近年来引起广泛关注的发展型积分 -微分方程[4- 6] ,这方面的研究尚不成熟 .本文将研究一维抛物型积分 -微分方程半离散有限元近似的超收敛性质 ,证明了剖分单元上的 Lobatto点、Gauss点和拟 Lobatto点分别是函数、一阶和二阶导数逼近的超收敛点 ;并且在一定条件下证明了强超收敛二择一定理 ;在每个单元上 ,单元中点或…     

基于层次T网格的分片孔斯曲面重构

本文提出一种基于任意层次T网格的多项式(PHT)样条空间$S(3,3,1,1,T)$的一个新的曲面重构算法.该算法由分片插值于层次T网格上每个小矩形单元对应4个顶点的16个参数的孔斯曲面形式给出.对于一个给定的T网格和相应基点处的几何信息(函数值,两个一阶偏导数和混合导数值),可得到与$S(3,3,1,1,T)$的PHT样条曲面相同的结果,且曲面表达形式更简单,同时,在离散数据点的曲面拟合中,我们给出了自适应的曲面加细算法.数值算例显示,该自适应算法能够有效的拟合离散数据点.     

关于“分段点”可导性判定的一个定理

如何判断分段函数在分段点处可导性,并求出导数?通常的作法(1)先判断连续性,若不连续,必不可导.(2)如果连续,再按导数的定义求导,由于在分段点两侧,函数表达式可能不同,则一般要通过计算分段点处左右导数来判断.实际上,在函数连续的基础上,可借助导函数在分段点处的极限,来判定并求出分段点的导数.这是因为有如下的定理:     

求解瞬时温度场的有限元显式算法

在一公共节点为中心的各单元中,对于线性形函数,实际计算和数字实验表明,温度在单元各节点上的时间导数用它在中心节点上的时间导数表示是可取和合理的。由此可在对微分方程用有限元法进行空间离散的基础上得到单个节点温度的时间导数与其周围节点温度的关系,建立温度场的显式计算格式。它具有计算简捷的特点。用最大值原理对稳定性的分析导出了与稳式算法类似的稳定性条件。     

一维有限元后处理的EEP法的数学分析

利用一维投影型插值与有限元超收敛基本估计,对一类两点边值问题,严格证明了袁驷等人由单元能量投影(EEP)法获得的节点恢复导数,当有限元空间的次数不超过4时,具有最佳阶超收敛．理论分析圆满地解释了已有的数值结果．     

Finite element derivative interpolation recovery technique and superconvergence

A new finite element derivative recovery technique is proposed by using the polynomial interpolation method. We show that the recovered derivatives possess superconvergence on the recovery domain and ultraconvergence at the interior mesh points for finite element approximations to elliptic boundary problems. Compared with the well-known Z-Z patch recovery technique, the advantage of our method is that it gives an explicit recovery formula and possesses the ultraconvergence for the odd-order finite elements. Finally, some numerical examples are presented to illustrate the theoretical analysis.     

导数小片插值恢复技术与超收敛性

１．引言 有限元超收敛的研究自七十年代起至今方兴未艾．现有的研究工作基本遵循两种途径：一是找出有限元插值逼近的超收敛点,然后再利用插值弱估计等手段导出有限元解本身所具有的超收敛性质［１,２］;二是利用各种后处理技术,如平均技术,投影技术,外插技术和插值有限元技术等[3,6],来导出经过后处理的有限元解的超收敛性．近年来一种新的超收敛后处理技术,即所谓的“Ｚ-Ｚ导数小片恢复技术”,得到众多的研究［７-１１］,并被 Ｂａｂｕｓｋａ等人认为是用于渐进准确的后验误差估计效果最好的技术之一［１２］．这种技术是利用…     

The Derivative Patch Interpolating Recovery Technique for Finite Element Approximations

A derivative patch interpolating recovery technique is analyzed for the finite element approximation to the second order elliptic boundary value problems in two dimensional case.It is shown that the convergence rate of the recovered gradient admits superc onvergence on the recovered subdomain, and is two order higher than the optimal global convergence rate (ultracovergence) at an internal node point when even order finite element spaces and local uniform meshes are used.     

Mathematical analysis of Zienkiewicz—Zhu's derivative patch recovery technique

Zienkiewicz-Zhu's derivative patch recovery technique is analyzed for general quadrilateral finite elements. Under certain regular conditions on the meshes, the arithmetic mean of the absolute error of the recovered gradient at the nodal points is superconvergent for the second-order elliptic operators. For rectangular meshes and the Laplacian, the recovered gradient is superconvergent in the maximum norm at the nodal points. Furthermore, it is proved for a model two-point boundary-value problem that the recovery technique results in an “ultra-convergent” derivative recovery at the nodal points for quadratic finite elements when uniform meshes are used. © 1996 John Wiley & Sons, Inc.     

Recovery type superconvergence and a posteriori error estimates for control problems governed by Stokes equations

In this paper, we derive recovery type superconvergence analysis and a posteriori error estimates for the finite element approximation of the distributed optimal control governed by Stokes equations. We obtain superconvergence results and asymptotically exact a posteriori error estimates by applying two recovery methods, which are the patch recovery technique and the least-squares surface fitting method. Our results are based on some regularity assumption for the Stokes control problems and are applicable to the first order conforming finite element method with regular but nonuniform partitions.     

Pantograph-Catenary Dynamics: An Analysis of Models and Simulation Techniques

In this paper we analyze models and simulation techniques for the interaction of pantograph and catenary. Detailed models for catenary and pantograph and the propagation of waves are first investigated. Next, the semi discretization by the finite element method and the time integration are described. In this context numerical techniques like GGL-stabilization and superconvergent patch recovery are applied. The latter yields an error estimation for the finite element grid and shows the critical points of the system.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

A BOUNDARY ELEMENT PROCEDURE FOR THE SIGNORINI PROBLEM IN THREE-DIMENSIONAL ELASTICITY

In this paper, a new numerical method for the Signorini problem in three-dimensional elasticity is presented. The problem is reduced to a boundary variational inequality based on a new representation of the derivative of the doublelayer potential. Furthermore, a boundary element procedure is described for the numerical approximation of its solution and an abstract error estimate is given.     

New construction and ultraconvergence of derivative recovery operator for odd-degree rectangular elements

In this paper the ultraconvergence of the derivative for odd-degree rectangular elements is addressed. A new, discrete least-squares patch recovery technique is proposed to post-process the solution derivatives. Such recovered derivatives are shown to possess ultraconvergence by using projection type interpolation.