求解奇异摄动问题混合有限元的最优一致收敛的统一方法

给出了在些Shiskin型网格[21,23,19,18]上,利用一个任意次的混合有限元方法在L2一模下得到奇异摄动问题解的最优一致收敛阶的一个统一方法,通过研究一个四阶问题,定常和不定常问题,我们显示了这个方法的一般性,结果显示非传统Shiskin型网格上的误差估计比传统Shiskin型网格上的误差估计更容易得到,但两种网格给出的误差估计是相容的,它们证明了Roos的猜想[21]是合理的。     

关于半离散Galerkin近似解的L_∞估计的一个注记

在《计算数学》和《高等学校计算数学学报》上最近发表的文章[1]和[2]中,分别讨论了抛物和二阶双曲方程半离散Galerkin近似解(分片线性函数情形)的L_∞估计。文章作者采用正则Green函数方法证明了阶为h~2ln(1/h)的误差估计式。值得指出,[1]和[2]中所给出的估计式的一个不足之处就是它们所需要的精确解的正则性过于强。在这个注记里,我们将说明如下事实,利用熟知的半离散Galerkin近似解的超收敛估计和有限元函数空间的一个弱嵌入性质,可以证明得到阶也是h~2ln(1/h)的误差估计式,然而对解的正则性的要求则较[1]和[2]中估计式所需要的弱得多。 先讨论抛物问题,文[1]讨论的是热传导问题     

定常Navier-Stokes方程流函数形式两重网格算法的残量型后验误差估计

运用七种两重网格协调元方法得出了不可压Navier-Stokes方程流函数形式的残量型后验误差估计．对比标准有限元方法的后验误差估计,两重网格算法的后验误差估计多了一些额外项(三线性项)．说明了这些额外项在误差估计中对研究离散解渐近性的重要性,推出了对于最优网格尺寸,这些额外项的收敛阶不高于标准离散解的收敛阶．     

混合有限元法的误差分析

关于混合变分问题有限元方法的研究工作,见[1]—[4].其中已得出混合法的最优误差估计.本文讨论抽象混合有限元法的误差并证明一些超收敛估计,然后将其应用到具体问题上,即应用到一个四阶边值问题和一个二阶边值问题.     

非协调Wilson有限元的多重网格方法

§1.引言 非协调Wilson有限元[1—3]对解弹性力学方程有实用价值,在工程上有用。本文分析Wilson元的多重网格法,给出用多重网格方法求得的近似解按L~2模和能量模的最佳收敛阶误差估计。对于W-循环,可以证明其计算量与离散空间的维数为同一量级O(N_k)。 考虑二阶椭圆Dirchlet边值问题:     

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

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

四阶障碍问题的稳定化混合有限元方法

本文讨论了四阶障碍问题的稳定化混合有限元方法.首先,引入网格依赖范数,通过加罚方法得到了与四阶障碍问题的等价的混合变分形式.随后给出了基于C~0协调有限元空间(W_h,V_h)的混合有限元逼近,例如P_k-P_k三角形有限元.在网格依赖范数下,(W_h,V_h)满足离散的inf-sup条件.最后,我们在不同的假设下,得到了一些误差估计.     

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

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

轴对称Stokes流的无限元逼近(Ⅱ)

[1]中讨论了无界区域上轴对称Stokes绕流的无限元方法,我们利用转移矩阵X以及组合刚度矩阵K_z将问题归结为一个有限阶代数方程组。[1]又给出了两种计算K_z的迭代方法,并证明了迭代方法的收敛性。最后证明了无限元解收敛于精确解,估计了误差的阶。这个方法的优点是:无穷远边界条件自然,计算规模小,边界形状不受限制,程序通用,并且理论基础比较完整。 本文是[1]的继续。我们将迭代格式作了一些简化,使之更便于计算;并且利用这种     

p—型有限元方法的L^2模误差估计的一个补充

本文讨论Kenneth Erikssion提出的模型问题的p-型有限元方法,解决了文[1]定理2后提出的问题,并给出提高误差收敛阶的一个方法。     

Graded Galerkin methods for the high‐order convection‐diffusion problem

We develop a Galerkin method using the Hermite spline on an admissible graded mesh for solving the high‐order singular perturbation problem of the convection‐diffusion type. We identify a special function class to which the solution of the convection‐diffusion problem belongs and characterize the approximation order of the Hermite spline for such a function class. The approximation order is then used to establish the optimal order of uniform convergence for the Galerkin method. Numerical results are presented to confirm the theoretical estimate.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

一类线性抛物型方程组的交替方向多步法及其理论分析

Efficient multistep procedure for time-stepping Galerkin method in which we use an alternating direction preconditioned iterative methods for approximately solving the linear equations arising at each timestep in a discrete Galerkin method for a class of linear parabolic systems is derived and analyzed. The optimal order error estimate is obtained. Numerical experiments show that the method has the characteristics of high efficiency and high accuracy.     

AD GALERKIN ANALYSIS FOR NONLINEAR PSEUDO—HYPERBOLIC EQUATIONS

AD(Alternating direction)Galerkin schemes for d-dimensional nonlinear pseudo-hyperbolic equations are studied.By using patch approximation technique,AD procedure is realized,and calculation,work is simplified.By using Galerkin approach,highly computational accuracy is kept.By using various priori estimate techniques for differential equations,difficulty coming form non-linearity is treated,and optimal H^1 and L^2 convergence prop-erties are demonstrated.Moreover,although all the existed AD Galerkin schemes using patch approximation are limited to have only one order accuracy in time increment,yet the schemes formulated in this paper have second order accuracy in it.This implies an essential advancement in AD Galerkin aualysis.     

Multilevel Monte Carlo method for parabolic stochastic partial differential equations

We analyze the convergence and complexity of multilevel Monte Carlo discretizations of a class of abstract stochastic, parabolic equations driven by square integrable martingales. We show under low regularity assumptions on the solution that the judicious combination of low order Galerkin discretizations in space and an Euler–Maruyama discretization in time yields mean square convergence of order one in space and of order 1/2 in time to the expected value of the mild solution. The complexity of the multilevel estimator is shown to scale log-linearly with respect to the corresponding work to generate a single path of the solution on the finest mesh, resp. of the corresponding deterministic parabolic problem on the finest mesh.     

A coupled continuous-discontinuous FEM approach for convection diffusion equations

In this article, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O((ε1 / 2 + h 1 / 2 )h k ) convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L 2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

A piecewise constant collocation method using cosine mesh grading for Symm's equation

Summary Here we study the piecewise constant collocation method using mesh grading to solve Symm's integral equation on [–1, 1]. We give a mesh grading for which this method achieves the optimal order of convergence even though the piecewise constant Galerkin method with the same mesh grading does not. Some numerical results are given.     

Convergence of a continuous Galerkin method with mesh modification for nonlinear wave equations

We consider space-time continuous Galerkin methods with mesh modification in time for semilinear second order hyperbolic equations. We show a priori estimates in the energy norm without mesh conditions. Under reasonable assumptions on the choice of the spatial mesh in each time step we show optimal order convergence rates. Estimates of the jump in the Riesz projection in two successive time steps are also derived.

     

Numerical analysis of a continuous Galerkin method for damped sine-Gordon equation

In this article, we discuss the numerical solution for the two-dimensional (2-D) damped sine-Gordon equation by using a space–time continuous Galerkin method. This method allows variable time steps and space mesh structures and its discrete scheme has good stability which are necessary for adaptive computations on unstructured grids. Meanwhile, it can easily get the higher-order accuracy in both space and time directions. The existence and uniqueness to the numerical solution are strictly proved and a priori error estimate in maximum-norm is given without any space–time grid conditions attached. Also, we prove that if the mesh in each time level is generated in a reasonable way, we can get the optimal order of convergence in both temporal and spatial variables. Finally, the convergence rates are presented and analyzed by some numerical experiments to illustrate the validity of the scheme.     

Max-Norm Estimates for Galerkin Approximations of One-Dimensional Elliptic,Parabolic and Hyperbolic Problems with Mixed Boundary Conditions

The Galerkin methods are studied for two-point boundary value problems and the related one-dimensional parabolic and hyperbolic problems. The boundary value problem considered here is of non-adjoint from and with mixed boundary conditions. The optimal order error estimate in the max-norm is first derived for the boundary problem for the finite element subspace. This result then gives optimal order max-norm error estimates for the continuous and discrete time approximations for the evolution problems described above.