半线性抛物方程的时空有限元方法

讨论半线性抛物方程的连续Galerkin时空有限元方法,利用有限元方法和有限差分方法相结合的技巧,证明了弱解的存在唯一性,给出了时间最大模,空间L~2模,即L~∞(L~2)模误差估计.并给出数值算例证明了连续时空有限元方法对于半线性抛物方程的有效性.     

一类三维变动区域上的非线性抛物型方程之全离散有限元形式

1 引言 1986年,L.Cermak和M.Zlamal研究了半导体器件中杂质的重新分布,对具有活动边界的二维非线性扩散问题。给出在时间方向上是一阶精度的全离散有限元格式。证明了格式最优的H~1模和次最优的L~2模估计。1989年.P.Lesaint和R.Touzani对一维变动区域上的热传导方程。经过坐标变换,给出了在固定区域上的全离散有限元格式和最优的L~2模估计。1990年,梁国平和陈志明利用时空有限元,给出了变动区域上线性抛物型的方程的全离散变网格有限元格式。证明了最优的L~2收敛性。本文考虑了一类具有活动边界的三维     

对流扩散问题流线扩散法的最优误差估计

在本文中,我们考虑对流占优扩散问题流线扩散双线性有限元方法。原先的文献在ε≤h~2的条件下,得到了L~2-模最优误差估计,而本文则在ε≤h的条件下得到了相同估计。     

热传导方程一个新的H~1-Galerkin非协调混合有限元格式

对热传导方程提出了一个新的H~1-Galerkin非协调混合有限元格式,其逼近空间不需满足LBB相容性条件,且在不引进传统的Rutz投影的情况下,得到了与以往协调有限元方法相同的L~2-模和H~1-模的误差估计.     

带弱奇异核非线性积分微分方程的有限元分析

讨论了带弱奇异核的非线性抛物积分微分方程的Hermite型各向异性矩形元逼近.在各向异性网格下导出了关于Riesz投影的L~2和H~1模的误差估计.在半离散和向后欧拉全离散格式下,基于Riesz投影的性质并利用平均值技巧,分别得到了L~2模意义下的最优误差估计.     

反应扩散方程的连续时空有限元方法

本文研究了反应扩散方程的连续时空有限元方法.首先建立了其连续时空有限元格式并证明了有限元解的存在唯一性及稳定性.然后通过引入时空投影算子在没有时空网格限制的条件下给出其近似解在节点处的L~2,H~1最优范数估计以及全局L~2(L~2),L~2(H~1)最优范数估计.最后给出两个数值算例来验证方法的有效性与灵活性并说明结论的正确性.     

页岩气藏渗流模型的特征混合有限元数值模拟

将特征有限元方法和混合有限元方法进行耦合,对页岩气藏渗流模型进行了数值模拟,给出了详细的误差分析,得到了最优的L~2模误差估计,并用数值实验验证了方法的有效性.     

电报方程H~1-Galerkin非协调混合有限元分析

主要研究一类电报方程的H~1-Galerkin非协调混合有限元方法,在任意四边形网格剖分下,其逼近空间分别取为类Wilson元与双线性Q_1元,在不需要满足LBB相容性条件及不采用传统的Ritz投影的情况下,得到了与常规有限元方法相同的L~2-模和H~1-模的误差估计,进一步拓展了H~1-Galerkin混合有限元和类Wilson元的应用范围.     

期权定价问题的数值方法

本文研究以股票为标的资产的美式看跌期权定价问题的数值方法,即有限元方法。通过将所考虑的问题转化为等价的变分不等式,并利用积分恒等式与超逼近分析技术,得到了半离散有限元方法的最优L~2-模与L~∞-模的误差估计。     

多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.     

The Crouzeix-Raviart type nonconforming finite element method for the nonstationary Navier-Stokes equations on anisotropic meshes

This paper is devoted to study the Crouzeix-Raviart （C-R） type nonconforming linear triangular finite element method （FEM） for the nonstationary Navier-Stokes equations on anisotropic meshes. By intro- ducing auxiliary finite element spaces, the error estimates for the velocity in the L2-norm and energy norm, as well as for the pressure in the L2-norm are derived.     

A symmetric finite volume scheme for selfadjoint elliptic problems

Based on a linear finite element space, a symmetric finite volume scheme for a self-adjoint elliptic boundary-value problem is proposed. Error estimates in L²-norm, H¹-norm, and L^∞-norm are derived. Some post-processing techniques are also provided.     

Convergence and stability of a stabilized finite volume method for the stationary Navier-Stokes equations

In this paper, a new stabilized finite volume method is studied and developed for the stationary Navier-Stokes equations. This method is based on a local Gauss integration technique and uses the lowest equal order finite element pair P ₁–P ₁ (linear functions). Stability and convergence of the optimal order in the H ¹-norm for velocity and the L ²-norm for pressure are obtained. A new duality for the Navier-Stokes equations is introduced to establish the convergence of the optimal order in the L ²-norm for velocity. Moreover, superconvergence between the conforming mixed finite element solution and the finite volume solution using the same finite element pair is derived. Numerical results are shown to support the developed convergence theory.     

New a posteriori L∞(L2) and L2(L2)-error estimates of mixed finite element methods for general nonlinear parabolic optimal control problems

We study new a posteriori error estimates of the mixed finite element methods for general optimal control problems governed by nonlinear parabolic equations. The state and the co-state are discretized by the high order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions. We derive a posteriori error estimates in L^∞(J; L²Ω)-norm and L²(J; L²Ω)-norm for both the state, the co-state and the control approximation. Such estimates, which seem to be new, are an important step towards developing a reliable adaptive mixed finite element approximation for optimal control problems. Finally, the performance of the posteriori error estimators is assessed by two numerical examples.     

Nonconforming mixed finite element approximation to the stationary Navier-Stokes equations on anisotropic meshes

The main aim of this paper is to study the error estimates of a rectangular nonconforming finite element for the stationary Navier-Stokes equations under anisotropic meshes. That is, the nonconforming rectangular element is taken as approximation space for the velocity and the piecewise constant element for the pressure. The convergence analysis is presented and the optimal error estimates both in a broken H¹-norm for the velocity and in an L²-norm for the pressure are derived on anisotropic meshes.     

A coupling method based on new MFE and FE for fourth-order parabolic equation

In this article, a coupling method of new mixed finite element (MFE) and finite element (FE) is proposed and analyzed for fourth-order parabolic partial differential equation. First, the fourth-order parabolic equation is split into the coupled system of second-order equations. Then, an equation is solved by finite element method, the other equation is approximated by the new mixed finite element method, whose flux belongs to the square integrable space replacing the classical H(div;Ω) space. The stability for fully discrete scheme is derived, and both semi-discrete and fully discrete error estimates are obtained. Moreover, the optimal a priori error estimates in L ² and H ¹-norm for both the scalar unknown u and the diffusion term γ and a priori error estimate in (L ²)²-norm for its flux σ are derived. Finally, some numerical results are provided to validate our theoretical analysis.     

A class of nonconforming quadrilateral finite elements for incompressible flow

This paper focuses on the low-order nonconforming rectangular and quadrilateral finite elements approximation of incompressible flow.Beyond the previous research works,we propose a general strategy to construct the basis functions.Under several specific constraints,the optimal error estimates are obtained,i.e.,the first order accuracy of the velocities in H1-norm and the pressure in L2-norm,as well as the second order accuracy of the velocities in L2-norm.Besides,we clarify the differences between rectangular and quadrilateral finite element approximation.In addition,we give several examples to verify the validity of our error estimates.     

Error estimates of the finite element method for the exterior Helmholtz problem with a modified DtN boundary condition

A priori error estimates in the H¹- and L²-norms are established for the finite element method applied to the exterior Helmholtz problem, with modified Dirichlet-to-Neumann (MDtN) boundary condition. The error estimates include the effect of truncation of the MDtN boundary condition as well as that of discretization of the finite element method. The error estimate in the L²-norm is sharper than that obtained by the author [D. Koyama, Error estimates of the DtN finite element method for the exterior Helmholtz problem, J. Comput. Appl. Math. 200 (1) (2007) 21-31] for the truncated DtN boundary condition.     

Optimal convergence analysis of an immersed interface finite element method

We analyze an immersed interface finite element method based on linear polynomials on noninterface triangular elements and piecewise linear polynomials on interface triangular elements. The flux jump condition is weakly enforced on the smooth interface. Optimal error estimates are derived in the broken H ¹-norm and L ²-norm.     

Convergence and superconvergence analysis for nonlinear delay reaction–diffusion system with nonconforming finite element

In this article, we propose and analyze several numerical methods for the nonlinear delay reaction–diffusion system with smooth and nonsmooth solutions, by using Quasi-Wilson nonconforming finite element methods in space and finite difference methods (including uniform and nonuniform L1 and L2-1_σ schemes) in time. The optimal convergence results in the senses of L²-norm and broken H¹-norm, and H¹-norm superclose results are derived by virtue of two types of fractional Grönwall inequalities. Then, the interpolation postprocessing technique is used to establish the superconvergence results. Moreover, to improve computational efficiency, fast algorithms by using sum-of-exponential technique are built for above proposed numerical schemes. Finally, we present some numerical experiments to confirm the theoretical correctness and show the effectiveness of the fast algorithms.