高阶偏微分方程局部间断Galerkin方法的最优误差估计（英文）

本文针对含三阶和四阶空间导数的高阶偏微分方程,得到了基于广义交替数值通量局部间断Galerkin方法的最优L~2-模误差估计.主要技术是基于有关辅助变量的能量方程和最新提出的整体Gauss-Radau投影.数值实验验证了理论结果.     

具转向点常微分方程混合边值问题的一致收敛差分格式

本文得到了具转向点的二阶常微分方程混合边值问题解的导数估计,提出了Il′in型差分格式,证明了此差分格式按L~1模关于小参数ε的一阶一致收敛性。最后,给出了一个数值例子,计算结果与理论分析一致。     

各向异性网格下抛物方程一个新的非协调混合元收敛性分析

本文将 Crouzeix-Raviart 型非协调线性三角形元应用到抛物方程,建立了一个新的混合元格式.在抛弃传统有限元分析的必要工具 Ritz 投影算子的前提下,直接利用单元的插值性质和导数转移技巧, 分别得到了各向异性剖分下关于原始变量u 的H^-1-模和积分意义下L²-模以及通量p=-▽u 在L²-模下的最优阶误差估计.数值结果与我们的理论分析是相吻合的.     

时间分数阶扩散方程双线性元的高精度分析

针对具有Caputo导数的二维时间分数阶扩散方程进行高精度有限元分析.首先,基于双线性元和L1逼近建立了一个全离散格式,并证明其在H~1模意义下的无条件稳定性;其次,借助Riesz投影和分数阶导数的技巧得到了L~2模意义下的最优误差估计,结合该元插值算子与Riesz投影算子之间的高精度结果和插值后处理技术,导出了H~1意义下的超逼近性质和超收敛结果.该结果是单独利用双线性插值算子和Riesz投影算子均无法得到的.最后,利用数值算例验证了理论分析的正确性.     

一阶双曲方程的Legendre—Tau方法

文章采用Legendre—tau方法对一阶双曲方程进行数值求解,此方法可以被有效实施,且可以得到L^2模意义下的最优误差估计,将以往对此类问题的收敛阶估计由O（N^1-τ）提高到O（N^-τ）,改进了原有的理论分析结果,数值算例证实了此方法的有效性．     

非线性分数阶常微分方程的分段线性插值多项式方法

通过分段线性插值多项式方法构造了一类含有Hadamard有限部分积分的非线性常微分方程的数值离散格式.在时间方向上, 利用分段线性插值多项式方法对分数阶导数项进行近似, 并通过二阶向后差分格式来离散整数阶导数项.经过详细的证明, 得到了收敛精度为O(τmin{1+α,1+β})的误差估计结果.最后,通过数值算例和理论结果的对比直观地说明了理论分析的正确性.     

一维非线性抛物问题两层网格有限体积元逼近

该文主要研究一维非线性抛物问题两层网格有限体积元逼近.对一维非线性抛物问题有限体积元解的存在性进行了讨论,给出了最优阶L~2-模和H~1-模误差估计结果,并研究了其两层网格算法.证明了当粗细网格步长满足h=O(H~2)时两层网格算法具有最优阶H~1-模误差估计.数值算例验证了理论结果.     

电报方程的类Wilson非协调有限元分析

本文在矩形网格上讨论了半离散和全离散格式下电报方程的类Wilson非协调有限元逼近.利用该元在H1模意义下O(h2)阶的相容误差结果,平均值理论和关于时间t的导数转移技巧得到了超逼近性.进而,借助于插值后处理方法导出了超收敛结果.又由于该元在H1模意义下的相容误差可以达到O(h3)阶,构造了新的外推格式,给出了比传统误差估计高两阶的外推估计.最后,对于给出的全离散逼近格式得到了最优误差估计.     

一类具耗散项双曲系统解的存在性问题

讨论具有非线性耗散项双曲系统的初值问题,对初值的模不加小性假设,而要求其一阶导数适当小情形下,证明其光滑解的整体存在性,并用经典解的特征线法获得解的模估计,同时应用极值原理得到解的偏导数的一致估计.     

Cahn-Hilliard方程的有限元分析

建立了求解非线性发展型Cahn-Hilliard方程的有限元方法,借助于一个双调和问题的有限元投影逼近,给出了最优阶L_2模误差估计。特别对于3次Hermite型有限元,导出了L_∞模和W_∞~1模的最优阶误差估计和导数逼近的超收敛结果。     

Galerkin alternating-direction method for a kind of second-order quasi-linear hyperbolic problems

A kind of second-order quasi-linear hyperbolic equation is firstly transformed into a first-order system of equations, then the Galerkin alternating-direction procedure for the system is derived. The optimal order estimates in H¹ norm and L² norm of the procedure are obtained respectively by using the theory and techniques of priori estimate of differential equations. The numerical experiment is also given to support the theoretical analysis. Comparing the results of numerical example with the theoretical analysis, they are uniform.     

Analysis of a two‐stage least‐squares finite element method for the planar elasticity problem

A new first‐order formulation for the two‐dimensional elasticity equations is proposed by introducing additional variables which, called stresses here, are the derivatives of displacements. The resulted stress–displacement system can be further decomposed into two dependent subsystems, the stress system and the displacement system recovered from the stresses. For constructing finite element approximations to these subsystems with appropriate boundary conditions, a two‐stage least‐squares procedure is introduced. The analysis shows that, under suitable regularity assumptions, the rates of convergence of the least‐squares approximations for all the unknowns are optimal both in the H¹‐norm and in L²‐norm. Also, numerical experiments with various Poisson's ratios are examined to demonstrate the theoretical estimates. Copyright © 1999 John Wiley & Sons, Ltd.     

A splitting positive definite mixed finite element method for two classes of integro-differential equations

In this article, we establish a new mixed finite element procedure to solve the second-order hyperbolic and pseudo-hyperbolic integro-differential equations, in which the mixed element system is symmetric positive definite without requiring the LBB consistency condition. Convergence analysis shows that the method yields the approximate solutions with optimal accuracy in L ²(??) norm for u and in H(?div;??) norm for the flux???. Numerical experiments are given to verify the theoretical results.     

Stabilized mixed finite element methods for the Navier‐Stokes equations with damping

In this paper, the stabilized mixed finite element methods are presented for the Navier‐Stokes equations with damping. The existence and uniqueness of the weak solutions are proven by use of the Brouwer fixed‐point theorem. Then, optimal error estimates for the H¹‐norm and L²‐norm of the velocity and the L²‐norm of the pressure are derived. Moreover, on the basis of the optimal L²‐norm error estimate of the velocity, a stabilized two‐step method is proposed, which is more efficient than the usual stabilized methods. Finally, two numerical examples are implemented to confirm the theoretical analysis.     

Superconvergence analysis of finite element method for Poisson–Nernst–Planck equations

This article concerns with the superconvergence analysis of bilinear finite element method (FEM) for nonlinear Poisson–Nernst–Planck (PNP) equations. By employing high accuracy integral identities together with mean value technique, the superclose estimates in H¹‐norm are derived for the semi‐discrete and the backward Euler fully‐discrete schemes, which improve the suboptimal error estimate in L²‐norm in the previous literature. Furthermore, the global superconvergence results in H¹‐norm are obtained through interpolation postprocessing approach. Finally, a numerical example is provided to confirm the theoretical analysis.     

Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L²-norm. The numerical examples are given to illustrate the theoretical results.     

Numerical approximation of solution for the coupled nonlinear Schrödinger equations

In this article, a compact finite difference scheme for the coupled nonlinear Schrödinger equations is studied. The scheme is proved to conserve the original conservative properties. Unconditional stability and convergence in maximum norm with order O(τ² + h⁴) are also proved by the discrete energy method. Finally, numerical results are provided to verify the theoretical analysis.     

Numerical analysis of a Stokes interface problem based on formulation using the characteristic function

Numerical analysis of a model Stokes interface problem with the homogeneous Dirichlet boundary condition is considered. The interface condition is interpreted as an additional singular force field to the Stokes equations using the characteristic function. The finite element method is applied after introducing a regularization of the singular source term. Consequently, the error is divided into the regularization and discretization parts which are studied separately. As a result, error estimates of order h^1/2 in H¹ × L² norm for the velocity and pressure, and of order h in L² norm for the velocity are derived. Those theoretical results are also verified by numerical examples.     

Superconvergent estimates of conforming finite element method for nonlinear time‐dependent Joule heating equations

In this article, we study the superconvergence analysis of conforming bilinear finite element method (FEM) for nonlinear Joule heating equations. Based on the rigorous estimates together with high accuracy analysis of this element, mean value technique and interpolation postprocessing approach, the superclose and superconvergent estimates about the related variables in H¹‐norm are derived for semidiscrete and a linearized backward Euler fully discrete schemes, which extends the results of optimal estimates obtained for conforming FEMs in the previous literature. At last, a numerical experiment is performed to verify the theoretical analysis.     

Conjugate gradient-boundary element solution using multiple reciprocity method for distributed elliptic optimal control problems

This paper deals with the numerical solution of optimal control problems, where the state equations are given by the fourth order elliptic partial differential equations. An iterative algorithm for this class of problems is developed. This new proposal is obtained by combining the Conjugate Gradient Method (CGM) with the Boundary Element Method (BEM) and Multiple Reciprocity Method (MRM). The local error estimates based on the stability of this scheme in the H² norm, L² norm and L^∞ norm are obtained. Finally, the numerical results on a test case show that this method is correct and feasible.