Uniform Convergence of Adaptive Multigrid Methods for Elliptic Problems and Maxwell's Equations

We consider the convergence theory of adaptive multigrid methods for second-order elliptic problems and Maxwell's equations. The multigrid algorithm only performs pointwise Gauss-Seidel relaxations on new degrees of freedom and their "immediate" neighbors. In the context of lowest order conforming finite element approximations, we present a unified proof for the convergence of adaptive multigrid V-cycle algorithms. The theory applies to any hierarchical tetrahedral meshes with uniformly bounded shape-regularity measures. The convergence rates for both problems are uniform with respect to the number of mesh levels and the number of degrees of freedom. We demonstrate our convergence theory by two numerical experiments.     

非自共轭和不定椭圆问题的有限体积元方法的一致收敛性

In this paper, we prove the existence, uniqueness and uniform convergence of the solution of finite volume element method based on the P1 conforming element for non-selfadjoint and indefinite elliptic problems under minimal elliptic regularity assumption.     

A Posteriori Error Estimators for Elliptic Equations with Discontinuous Coefficients

We consider linear elliptic equations with discontinuous coefficients in two and three space dimensions with varying boundary conditions. The problem is discretized with linear finite elements. An adaptive procedure based on a posteriori error estimators for the treatment of singularities is proposed. Within the class of quasi-monotonically distributed coefficients we derive a posteriori error estimators with bounds that are independent of the variation of the coefficients. In numerical test cases we confirm the robustness of the error estimators and observe that on adaptively refined meshes the reduction of the error is optimal with respect to the number of unknowns.     

Local Multilevel Method on Adaptively Refined Meshes for Elliptic Problems with Smooth Complex Coefficients

In this paper, a local multilevel algorithm is investigated for solving linear systems arising from adaptive finite element approximations of second order elliptic problems with smooth complex coefficients. It is shown that the abstract theory for local multilevel algorithm can also be applied to elliptic problems whose dominant coefficient is complex valued. Assuming that the coarsest mesh size is sufficiently small, we prove that this algorithm with Gauss-Seidel smoother is convergent and optimal on the adaptively refined meshes generated by the newest vertex bisection algorithm. Numerical experiments are reported to confirm the theoretical analysis.     

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

<正>Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.     

双障碍问题的多重网格法

本文研究双障碍问题的多重网格法,提出了两类算法,证明了其收敛性及对贴合分量的有限步收敛性,同时对其中一种算法的特款提出了一个 k无关收敛性定理。     

溃坝问题的间断有限元方法

本文研究90年代初提出的Runge-Kutta间断Galerkin有限元方法，给出该方法的精度分析，通过经典算例验证该方法处理间断问题、捕捉锐利波形的能力，并将其推广到求解浅水问题．针对坝底无摩擦，无坡度的理想情形进行讨论，给出方溃坝和圆溃坝问题的数值模拟结果．     

Preconditioning Mixed Finite Element Saddle-point Elliptic Problems

We consider saddle-point problems that typically arise from the mixed finite element discretization of second-order elliptic problems. By proper equivalent algebraic operations the considered saddle-point problem is transformed to another saddle-point problem. The resulting problem can then be efficiently preconditioned by a block-diagonal matrix or by a factored block-matrix (the blocks correspond to the velocity and pressure, respectively). Both preconditioners have a block on the main diagonal that corresponds to the bilinear form (δ is a positive parameter) and a second block that is equal to a constant times the identity operator. We derive uniform bounds for the negative and positive eigenvalues of the preconditioned operator. Then any known preconditioner for the above bilinear form can be applied. We also show some numerical experiments that illustrate the convergence properties of the proposed technique.     

Second Order Multigrid Methods for Elliptic Problems with Discontinuous Coefficients on an Arbitrary Interface,I: One Dimensional Problems

In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface. Second order accuracy for the first derivative is obtained as well. The method is based on the Ghost Fluid Method, making use of ghost points on which the value is defined by suitable interface conditions. The multi-domain formulation is adopted, where the problem is split in two sub-problems and interface conditions will be enforced to close the problem. Interface conditions are relaxed together with the internal equations (following the approach proposed in [10] in the case of smooth coefficients), leading to an iterative method on all the set of grid values (inside points and ghost points). A multigrid approach with a suitable definition of the restriction operator is provided. The restriction of the defect is performed separately for both sub-problems, providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient. Numerical tests will confirm the second order accuracy. Although the method is proposed in one dimension, the extension in higher dimension is currently underway [12] and it will be carried out by combining the discretization of [10] with the multigrid approach of [11] for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.     

半线性椭圆型问题Mortar有限元逼近的瀑布型多重网格法

Mortar有限元法作为一个非协调的区域分解技术已得到许多研究者的关注(如文献[2]、[5]等)。本文对半线性椭圆型问题的Mortar有限元逼近提出了瀑布型多重网格法,并给出了此法的误差估计和计算复杂度估计定理。     

Uniform Convergence of Trigonometric Series with Rarely Changing Coefficients

We consider the series and whose coefficients satisfy the condition for , where the sequence can be expressed as the union of a finite number of lacunary sequences. The following results are obtained. If as , then the series is uniformly convergent. If for all , then the sequence of partial sums of this series is uniformly bounded. If the series is convergent for and as , then this series is uniformly convergent. If the sequence of partial sums of the series for is bounded and for all , then the sequence of partial sums of this series is uniformly bounded. In these assertions, conditions on the rates of decrease of the coefficients of the series are also necessary if the sequence is lacunary. In the general case, they are not necessary.     

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

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

Cascadic Multigrid for Finite Volume Methods for Elliptic Problems

In this paper, some effective cascadic multigrid methods are proposed for solving the large scale symmetric or nonsymmetric algebraic systems arising from the finite volume methods for second order elliptic problems. It is shown that these algorithms are optimal in both accuracy and computational complexity. Numerical experiments are reported to support our theory.     

Newton-Multigrid for Biological Reaction-Diffusion Problems with Random Coefficients

An algebraic Newton-multigrid method is proposed in order to efficiently solve systems of nonlinear reaction-diffusion problems with stochastic coefficients. These problems model the conversion of starch into sugars in growing apples. The stochastic system is first converted into a large coupled system of deterministic equations by applying a stochastic Galerkin finite element discretization. This method leads to high-order accurate stochastic solutions. A stable and high-order time discretization is obtained by applying a fully implicit Runge-Kutta method. After Newton linearization, a point-based algebraic multigrid solution method is applied. In order to decrease the computational cost, alternative multigrid preconditioners are presented. Numerical results demonstrate the convergence properties, robustness and efficiency of the proposed multigrid methods.     

Low Regularity Error Analysis for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

This paper presents error estimates in both an energy norm and the $L^2$-norm for the weak Galerkin (WG) finite element methods for elliptic problems with low regularity solutions. The error analysis for the continuous Galerkin finite element remains same regardless of regularity. A totally different analysis is needed for discontinuous finite element methods if the elliptic regularity is lower than H-1.5. Numerical results confirm the theoretical analysis.     

一类四阶奇异非线性椭圆方程的加权模误差估计

利用有限元方法研究了一类四阶奇异非线性椭圆方程,利用有限元的逆性质,给出了考虑数值积分影响时的加权H2模误差估计.     

有限元网格的超收敛测度及其应用

给出线性有限元求解二阶椭圆问题的有限元网格超收敛测度及其应用.有限元超收敛经常是在具有一定结构的特殊网格条件下讨论的,而本文从一般网格出发,导出一种网格的范数用来描述超收敛所需要的网格条件以及超收敛的程度.并且通过对这种网格范数性质的考察,可以证明对于通常考虑的一些特殊网格的超收敛的存在性.更进一步,我们可以通过正则细分的方式在一般区域上也可以自动获得超收敛网格.最后给出相关的数值结果来验证本文的理论分析.     

Weak Convergence of a Sequence of Semimartingales to a Diffusion with Discontinuous Drift and Diffusion Coefficients

This paper presents a set of sufficient conditions for a sequence of semimartingales to converge weakly to a solution of a stochastic differential equation (SDE) with discontinuous drift and diffusion coefficients. This result is closely related to a well-known weak-convergence theorem due to Liptser and Shiryayev (see [27]) which proves the weak convergence to a solution of a SDE with continuous drift and diffusion coefficients in the Skorokhod–Lindvall J ₁-topology.The goal of this paper is to obtain a stronger result in order to solve outstanding problems in the area of large-scale queueing networks – in which the weak convergence of normalized queueing length is a solution of a SDE with discontinuous coefficients. To do this we need to make the stronger assumptions: (1) replacing the convergence in probability of the triplets of a sequence of semimartingales in the original Liptser and Shiryayev's theorem by stronger convergence in L ₂, (2) assuming the diffusion coefficient is coercive, and (3) assuming the discontinuity sets of the coefficients of the limit diffusion processs are of Lebesgue measure zero.     

Uniform Convergence Analysis for Singularly Perturbed Elliptic Problems with Parabolic Layers

In this paper,using Lin's integral identity technique,we prove the optimal uniform convergence θ(N_x~(-2)In~2N_x N_y~(-2)In~2N_y) in the L~2-norm for singularly per- turbed problems with parabolic layers.The error estimate is achieved by bilinear fi- nite elements on a Shishkin type mesh.Here N_x and N_y are the number of elements in the x- and y-directions,respectively.Numerical results are provided supporting our theoretical analysis.