一类带参数的非线性奇异摄动问题的自适应移动网格算法

本文讨论一类带参数的非线性奇异摄动问题的自适应移动网格方法.首先,在任意非均匀网格下,利用向后欧拉公式对方程进行离散,并给出相应的局部截断误差.然后,基于局部截断误差和网格等分布原理,利用精确解的弧长函数,证明半离散格式下自适应移动网格算法是一阶收敛的.同时,基于近似的弧长控制函数,给出易于实现的网格生成算法,并给出全离散格式下的后验误差估计.最后,数值实验结果验证了本文所给出的理论结果.     

一类强耦合的奇异摄动对流扩散方程组的移动网格方法（英文）

本文研究一类强耦合的奇异摄动对流扩散方程组的移动网格方法.首先,利用迎风有限差分格式对方程组进行离散.然后,推出数值解的后验误差估计,并以此设计出相应的自适应网格生成算法.同时,证明数值解具有一阶一致收敛性.最后,数值实验验证了本文移动网格方法的一致收敛性.     

二维空间中分数阶扩散方程的变网格有限元方法

针对二维空间分数阶偏微分方程,给出了一个变网格全离散有限元格式,并得到了相应最优误差估计.其主要思想是对空间变量采用有限元离散,对时间交量采用差分,但不同时刻的有限元网格可以不同.这对于没计相应的自适应算法是十分有益的.     

对流扩散方程的时空间断有限元方法

研究对流扩散方程的时空间断Galerkin有限元方法,该方法采用时,空两个变量都允许间断的基函数,更适用于移动网格,自适应算法以及并行计算.本文利用拉格朗日欧拉方法,采用F.Brezzi数值流通量,给出对流扩散方程的间断时空有限元离散格式,并证明格式的相容性,强制性,稳定性,解的存在唯一性,以及总体误差估计.     

固定利率抵押贷款模型的基于自适应网格的有限差分法

研究求解固定利率抵押贷款模型的基于自适应网格的有限差分策略.采用中心差分格式来离散微分算子的空间变量导数项,构造分片一致的自适应网格,使得与离散算子相应的系数矩阵为M-阵,以保证所构造的差分策略是在无穷模意义下稳定的.通过区分不同网格点集,在相应的网格点集上应用极大模原理来直接导出误差估计.此有限差分策略对于任意波动率和任意利率都是稳定的,并且是关于标的资产价格二阶收敛的.数值实验证实了理论结果的准确性.     

中子输运方程误差估计及自适应计算

 本文研究了全离散方法求解二维中子输运方程的有限元自适应算法, 角度变量用离散纵坐标方法展开, 空间变量用间断元方法求解. 基于间断元方法给出了空间离散的残量型后验误差估计. 在后验误差估计的基础上, 我们设计了自适应有限元算法.由残量型后验估计可以给出局部加密网格的自适应算法. 最后, 我们给出了数值算例来验证我们的理论结果.     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

Sobolev方程的一类各向异性非协调有限元逼近

在各向异性网格下,分别讨论了Sobolev方程在半离散和全离散格式下的一类非协调有限元逼近,得到了与传统有限元方法相同的误差估计和一些超逼近性质.同时在半离散格式下,通过构造具有各向异性特征的插值后处理算子得到了整体超收敛结果.     

粒子输运方程的子网格平衡格式的稳定性和收敛性

本文研究四边形网格上求解粒子输运方程的有限体积格式,其中角方向变量采用离散纵标(Sn)方法,空间离散采用子网格平衡(SCB)格式.利用能量估计方法,证明了在正交网格上该格式的稳定性和离散解的收敛性.数值实验结果验证了格式的稳定性和离散解的收敛性.     

自适应隐式有限差分方法求解美式期权定价问题

本文针对美式期权的定价问题设计了基于有限差分方法的预估-校正数值算法.该算法采用显式离散格式先对自由边界条件进行预估,再对经过变量替换后的关于期权价格的偏微分方程采用隐式格式离散,并用Fourier方法分析了此离散格式的稳定性.接下来,引入基于Richardson外推法的后验误差指示子.这个后验误差指示子能够在给定的误差阈值范围内,针对期权价格和自由边界找到合适的网格划分.最后,通过设计多组数值实验并与Fazio^[1]采用显式离散格式算得的数值结果相比较,验证了所提算法的有效性,稳定性和收敛性.     

An accurate a posteriori error estimator for the Steklov eigenvalue problem and its applications

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine mesh and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.

     

Steklov eigenvalue problem and its applications An accurate a posteriori error estimator for the

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine meshes and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.     

A Regularized Iterative Scheme for Solving Nonlinear Ill-Posed Problems

In this article, we consider a regularized iterative scheme for solving nonlinear ill-posed problems. The convergence analysis and error estimates are derived by choosing the regularization parameter according to both a priori and a posteriori methods. The iterative scheme is stopped using an a posteriori stopping rule, and we prove that the scheme converges to the solution of the well-known Lavrentiev scheme. The salient features of the proposed scheme are: (i) convergence and error estimate analysis require only weaker assumptions compared to standard assumptions followed in literature, and (ii) consideration of an adaptive a posteriori stopping rule and a parameter choice strategy that gives the same convergence rate as that of an a priori method without using the smallness assumption, the source condition. The above features are very useful from theory and application points of view. We also supply the numerical results to illustrate that the method is adaptable. Further, we compare the numerical result of the proposed method with the standard approach to demonstrate that our scheme is stable and achieves good computational output.     

A posteriori error analysis of the discontinuous finite element methods for first order hyperbolic problems

In this paper, we consider the a posteriori error analysis of discontinuous Galerkin finite element methods for the steady and nonsteady first order hyperbolic problems with inflow boundary conditions. We establish several residual-based a posteriori error estimators which provide global upper bounds and a local lower bound on the error. Further, for nonsteady problem, we construct a fully discrete discontinuous finite element scheme and derive the a posteriori error estimators which yield global upper bound on the error in time and space. Our a posteriori error analysis is based on the mesh-dependent a priori estimates for the first order hyperbolic problems. These a posteriori error analysis results can be applied to develop the adaptive discontinuous finite element methods.     

An augmented mixed finite element method with Lagrange multipliers: A priori and a posteriori error analyses

In this paper, we provide a priori and a posteriori error analyses of an augmented mixed finite element method with Lagrange multipliers applied to elliptic equations in divergence form with mixed boundary conditions. The augmented scheme is obtained by including the Galerkin least-squares terms arising from the constitutive and equilibrium equations. We use the classical Babuška–Brezzi theory to show that the resulting dual-mixed variational formulation and its Galerkin scheme defined with Raviart–Thomas spaces are well posed, and also to derive the corresponding a priori error estimates and rates of convergence. Then, we develop a reliable and efficient residual-based a posteriori error estimate and a reliable and quasi-efficient Ritz projection-based one, as well. Finally, several numerical results illustrating the performance of the augmented scheme and the associated adaptive algorithms are reported.     

Numerical Analysis of a Nonlinear Singularly Perturbed Delay Volterra Integro-Differential Equation on an Adaptive Grid

In this paper, we study a nonlinear first-order singularly perturbed Volterra integro-differential equation with delay. This equation is discretized by the backward Euler for differential part and the composite numerical quadrature formula for integral part for which both an a priori and an a posteriori error analysis in the maximum norm are derived. Based on the a priori error bound and mesh equidistribution principle, we prove that there exists a mesh gives optimal first order convergence which is robust with respect to the perturbation parameter. The a posteriori error bound is used to choose a suitable monitor function and design a corresponding adaptive grid generation algorithm. Furthermore, we extend our presented adaptive grid algorithm to a class of second-order nonlinear singularly perturbed delay differential equations. Numerical results are provided to demonstrate the effectiveness of our presented monitor function. Meanwhile, it is shown that the standard arc-length monitor function is unsuitable for this type of singularly perturbed delay differential equations with a turning point.     

Locking-free adaptive discontinuous Galerkin FEM for linear elasticity problems

An adaptive discontinuous Galerkin finite element method for linear elasticity problems is presented. We develop an a posteriori error estimate and prove its robustness with respect to nearly incompressible materials (absence of volume locking). Furthermore, we present some numerical experiments which illustrate the performance of the scheme on adaptively refined meshes.

     

An adaptive time stepping method with efficient error control for second-order evolution problems

This work is concerned with time stepping fnite element methods for abstract second order evolution problems.We derive optimal order a posteriori error estimates and a posteriori nodal superconvergence error estimates using the energy approach and the duality argument.With the help of the a posteriori error estimator developed in this work,we will further propose an adaptive time stepping strategy.A number of numerical experiments are performed to illustrate the reliability and efciency of the a posteriori error estimates and to assess the efectiveness of the proposed adaptive time stepping method.     

An a posteriori error estimator for anisotropic refinement

Summary. Besides an algorithm for local refinement, an a posteriori error estimator is the basic tool of every adaptive finite element method. Using information generated by such an error estimator the refinement of the grid is controlled. For 2nd order elliptic problems we present an error estimator for anisotropically refined grids, like -d cuboidal and 3-d prismatic grids, that gives correct information about the size of the error; additionally it generates information about the direction into which some element has to be refined to reduce the error in a proper way. Numerical examples are presented for 2-d rectangular and 3-d prismatic grids. Received March 15, 1994 / Revised version received June 3, 1994     

A posteriori error estimate for discontinuous Galerkin finite element method on polytopal mesh

In this work, we derive a posteriori error estimates for discontinuous Galerkin finite element method on polytopal mesh. We construct a reliable and efficient a posteriori error estimator on general polygonal or polyhedral meshes. An adaptive algorithm based on the error estimator and DG method is proposed to solve a variety of test problems. Numerical experiments are performed to illustrate the effectiveness of the algorithm.