Legendre-Laguerre coupled spectral element methods for second- and fourth-order equations on the half line

Some Legendre spectral element/Laguerre spectral coupled methods are proposed to numerically solve second- and fourth-order equations on the half line. The proposed methods are based on splitting the infinite domain into two parts, then using the Legendre spectral element method in the finite subdomain and Laguerre method in the infinite subdomain. C⁰ or C¹-continuity, according to the problem under consideration, is imposed to couple the two methods. Rigorous error analysis is carried out to establish the convergence of the method. More importantly, an efficient computational process is introduced to solve the discrete system. Several numerical examples are provided to confirm the theoretical results and the efficiency of the method.     

Error estimates of CFVE method for fully nonlinear convection‐dominated diffusion problems

Finite volume method and characteristics finite element method are two important methods for solving the partial differential equations. These two methods are combined in this paper to establish a fully discrete characteristics finite volume method for fully nonlinear convection‐dominated diffusion problems. Through detailed theoretical analysis, optimal order H¹ norm error estimates are obtained for this fully discrete scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

Error estimates of H 1-Galerkin mixed finite element method for Schrödinger equation

An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.     

Convergence analysis of virtual element methods for semilinear parabolic problems on polygonal meshes

In this article, we discuss and analyze new conforming virtual element methods (VEMs) for the approximation of semilinear parabolic problems on convex polygonal meshes in two spatial dimension. The spatial discretization is based on polynomial and suitable nonpolynomial functions, and a Euler backward scheme is employed for time discretization. The discrete formulation of both the proposed schemes—semidiscrete and fully discrete (with time discretization) is discussed in detail, and the unique solvability of the resulted schemes is discussed. A priori error estimates for the proposed schemes (semidiscrete and fully discrete) in H¹‐ and L²‐norms are derived under the assumption that the source term f is Lipschitz continuous. Some numerical experiments are conducted to illustrate the performance of the proposed scheme and to confirm the theoretical convergence rates.     

求解两类Maxwell方程组棱元离散系统的快速算法和自适应方法

Maxwell方程组棱元离散系统的快速算法和自适应方法是当前计算电磁场中的研究热点和难点. 首先, 针对H(curl)椭圆方程组的棱元离散系统, 通过建立棱元空间的稳定性分解, 设计了相应的快速迭代法和高效预条件子, 并且证明了迭代算法的收敛率和预条件子的条件数均不依赖于模型参数和网格规模. 其次, 针对时谐Maxwell方程组的棱有限元方法, 利用离散的Helmholtz分解, 连续散度为零函数对离散散度为零函数的逼近性和对偶论证, 获得了在L²和H(curl)范数下的拟最优误差估计. 进而设计和分析了相应的两网格法. 最后, 分别针对变系数H(curl)椭圆方程组和不定时谐Maxwell方程组, 考虑了一种不需要标记振荡项和加密单元不需要满足“内节点” 性质的自适应棱有限元法(AEFEM), 并证明了AEFEM的收敛性. 进一步, 当初始网格和Dörfler标记策略参数满足一定的假设条件时, 利用AEFEM的收敛性、误差的整体下界和局部上界估计, 证明了AEFEM的拟最优复杂性.     

Error estimates for mixed finite element approximations of the viscoelasticity wave equation

This paper studies mixed finite element approximations to the solution of the viscoelasticity wave equation. Two new transformations are introduced and a corresponding system of first‐order differential‐integral equations is derived. The semi‐discrete and full‐discrete mixed finite element methods are then proposed for the problem based on the Raviart–Thomas–Nedelec spaces. The optimal error estimates in L²‐norm are obtained for the semi‐discrete and full‐discrete mixed approximations of the general viscoelasticity wave equation. Copyright © 2004 John Wiley & Sons, Ltd.     

Optimal l~∞ error estimates of finite difference methods for the coupled Gross-Pitaevskii equations in high dimensions

Due to the difficulty in obtaining the a priori estimate,it is very hard to establish the optimal point-wise error bound of a finite difference scheme for solving a nonlinear partial differential equation in high dimensions(2D or 3D).We here propose and analyze finite difference methods for solving the coupled GrossPitaevskii equations in two dimensions,which models the two-component Bose-Einstein condensates with an internal atomic Josephson junction.The methods which we considered include two conservative type schemes and two non-conservative type schemes.Discrete conservation laws and solvability of the schemes are analyzed.For the four proposed finite difference methods,we establish the optimal convergence rates for the error at the order of O(h~2+τ~2)in the l~∞-norm(i.e.,the point-wise error estimates)with the time stepτand the mesh size h.Besides the standard techniques of the energy method,the key techniques in the analysis is to use the cut-off function technique,transformation between the time and space direction and the method of order reduction.All the methods and results here are also valid and can be easily extended to the three-dimensional case.Finally,numerical results are reported to confirm our theoretical error estimates for the numerical methods.     

Optimal error estimates of mixed finite element methods for a fourth-order nonlinear elliptic problem

A mixed finite element method is developed for a nonlinear fourth-order elliptic problem. Optimal L² error estimates are proved by using a special interpolation operator on the standard tensor-product finite elements of order k?1. Then two iterative schemes are presented and proved to keep the same optimal error estimates. Three numerical examples are provided to support the theoretical analysis.     

Superconvergence analysis of splitting positive definite nonconforming mixed finite element method for pseudo-hyperbolic equations

In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.     

Benchmark results for testing adaptive finite element eigenvalue procedures

A discontinuous Galerkin method, with hp-adaptivity based on the approximate solution of appropriate dual problems, is employed for highly-accurate eigenvalue computations on a collection of benchmark examples. After demonstrating the effectivity of our computed error estimates on a few well-studied examples, we present results for several examples in which the coefficients of the partial-differential operators are discontinuous. The problems considered here are put forward as benchmarks upon which other adaptive methods for computing eigenvalues may be tested, with results compared to our own.     

Adaptive stabilized mixed finite volume methods for the incompressible flow

In this article, we study adaptive stabilized mixed finite volume methods for the incompressible flows approximated using the lower order elements. A residual type of a posteriori error estimator is designed and studied with the derivation of upper and lower bounds between the exact solution and the finite volume solution. A discrete local lower bound between two successive finite volume solutions is also obtained. Also, convergence of the adaptive stabilized mixed finite volume methods is established. The presented methods have three prominent features. First, it is of practical convenience in real applications with the same partitions for velocity and pressure. Second, less computational time is required by easily applying both the lower order elements and the local grid refinement necessary for the elements of interest. Third, compared with the standard finite element method, its analysis of H¹‐norm and L²‐norm for the velocity and pressure are usually derived without any high order regularity conditions on the exact solution. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1424–1443, 2015     

Superconvergence of triangular mixed finite elements for optimal control problems with an integral constraint

In this paper, we shall investigate the superconvergence property of quadratic elliptical optimal control problems by triangular mixed finite element methods. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite elements and the control is discretized by piecewise constant functions. We prove the superconvergence error estimate of h² in L²-norm between the approximated solution and the interpolation of the exact control variable. Moreover, by postprocessing technique, we find that the projection of the discrete adjoint state is superclose (in order h²) to the exact control variable.     

Attractors for Fully Discrete Finite Difference Scheme of Dissipative Zakharov Equations

A fully discrete finite difference scheme for dissipative Zakharov equations is analyzed. On the basis of a series of the time-uniform priori estimates of the difference solutions, the stability of the difference scheme and the error bounds of optimal order of the difference solutions are obtained in L²×H¹×H² over a finite time interval (0, T]. Finally, the existence of a global attractor is proved for a discrete dynamical system associated with the fully discrete finite difference scheme.     

Guaranteed velocity error control for the pseudostress approximation of the Stokes equations

The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L². Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

and Norms Error Estimates in Finite Element Method for Linear Parabolic Interface Problems

We derive new a priori error estimates for linear parabolic equations with discontinuous coefficients. Due to low global regularity of the solutions the error analysis of the standard finite element method for parabolic problems is difficult to adopt for parabolic interface problems. A finite element procedure is, therefore, proposed and analyzed in this paper. We are able to show that the standard energy technique of finite element method for non-interface parabolic problems can be extended to parabolic interface problems if we allow interface triangles to be curved triangles. Optimal pointwise-in-time error estimates in the L ²(Ω) and H ¹(Ω) norms are shown to hold for the semidiscrete scheme. A fully discrete scheme based on backward Euler method is analyzed and pointwise-in-time error estimates are derived. The interfaces are assumed to be arbitrary shape but smooth for our purpose.     

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.     

Optimal error estimates and modified energy conservation identities of the ADI-FDTD scheme on staggered grids for 3D Maxwell’s equations

This paper is concerned with the optimal error estimates and energy conservation properties of the alternating direction implicit finite-difference time-domain(ADI-FDTD) method which is a popular scheme for solving the 3D Maxwell’s equations.Precisely,for the case with a perfectly electric conducting(PEC) boundary condition we establish the optimal second-order error estimates in both space and time in the discrete H 1-norm for the ADI-FDTD scheme,and prove the approximate divergence preserving property that if the divergence of the initial electric and magnetic fields are zero,then the discrete L 2-norm of the discrete divergence of the ADI-FDTD solution is approximately zero with the second-order accuracy in both space and time.The key ingredient is two new discrete modified energy norms which are second-order in time perturbations of two new energy conservation laws for the Maxwell’s equations introduced in this paper.Furthermore,we prove that,in addition to two known discrete modified energy identities which are second-order in time perturbations of two known energy conservation laws,the ADI-FDTD scheme also satisfies two new discrete modified energy identities which are second-order in time perturbations of the two new energy conservation laws.This means that the ADI-FDTD scheme is unconditionally stable under the four discrete modified energy norms.Experimental results which confirm the theoretical results are presented.     

Two-grid methods for characteristic finite volume element solution of semilinear convection-diffusion equations

Two-grid methods for characteristic finite volume element solutions are presented for a kind of semilinear convection-dominated diffusion equations. The methods are based on the method of characteristics, two-grid method and the finite volume element method. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H). And the fine-grid solution (with grid size h) can be obtained by a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy H = O(h^1/3).     

Mortar element methods for parabolic problems

In this article a standard mortar finite element method and a mortar element method with Lagrange multiplier are used for spatial discretization of a class of parabolic initial‐boundary value problems. Optimal error estimates in L^∞(L²) and L^∞(H¹)‐norms for semidiscrete methods for both the cases are established. The key feature that we have adopted here is to introduce a modified elliptic projection. In the standard mortar element method, a completely discrete scheme using backward Euler scheme is discussed and optimal error estimates are derived. The results of numerical experiments support the theoretical results obtained in this article. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008