Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

Meso-Cenozoic Evolution of Earth Surface System under the East Asian Tectonic Superconvergence

The East Asian geological setting has a long duration related to the superconvergence of the Paleo-Asian, Tethyan and Paleo-Pacific tectonic domains. The Triassic Indosinian Movement contributed to an unified passive continental margin in East Asia. The later ophiolites and I-type granites associated with subduction of the Paleo-Pacific Plate in the Late Triassic, suggest a transition from passive to active continental margins. With the presence of the ongoing westward migration of the Paleo-Pacific Subduction Zone, the sinistral transpressional stress field could play an important role in the intraplate deformation in East Asia during the Late Triassic to Middle Jurassic, being characterized by the transition from the E-W-trending structural system controlled by the Tethys and Paleo-Asian oceans to the NE-trending structural system caused by the Paleo-Pacific Ocean subduction. The continuously westward migration of the subduction zones resulted in the transpressional stress field in East Asia marked by the emergence of the Eastern North China Plateau and the formation of the Andean-type active continental margin from late Late Jurassic to Early Cretaceous （160-135 Ma）, accompanied by the development of a small amount of adakites. In the Late Cretaceous （135-90 Ma）, due to the eastward retreat of the Paleo-Pacific Subduction Zone, the regional stress field was replaced from sinistral transpression to transtension. Since a large amount of late-stage adakites and metamorphic core complexes developed, the Andean-type active continental margin was destroyed and the Eastern North China Plateau started to collapse. In the Late Cretaceous, the extension in East Asia gradually decreased the eastward retreat of the Paleo-Pacific subduction zones. Futhermore, a significant topographic inversion had taken place during the Cenozoic that resulted from a rapid uplift of the Tibet Plateau resulting from the India-Eurasian collision and the formation of the Bohai Bay Basin and other basins in the East Asian continental margin. The inversion caused a remarkable eastward migration of deformation, basin formation and magmatism. Meanwhile, the basins that mainly developed in the Paleogene resulted in a three-step topography which typically appears to drop eastward in altitude. In the Neogene, the basins underwent a rapid subsidence in some depressions after basin-controlled faulting, as well as the intracontinental extensional events in East Asia, and are likely to be a contribution to the uplift of the Tibetan Plateau.

APPLICATION OF SUPERCONVERGENCE TO A MODEL FOR COMPRESSIBLE MISCIBLE DISPLACEMENT

A modification of a finite element method of Douglas and Roberts for approximating the solution of the equations describing compressible miscible displacement in a porous medium is proposed and analyzed. The pressure is treated by a parabolic mixed finite element method using a Raviart-Thomas space of index rover a quasiregular partition, An extension of the Darcy velocity along Gauss lines is used in the evaluation of the coefficients in the Galerkin procedure for the concentration. A simple computational procedure allows the superconvergence property of the fluid velocity to be retained in our total algorithm.

Superconvergence analysis of Wilson element on anisotropic meshes

The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain the superclose and global superconvergence on anisotropic meshes. Numerical example is also given to confirm our theoretical analysis.

Superconvergence of nonconforming finite element penalty scheme for Stokes problem using L^2 projection method

A modified penalty scheme is discussed for solving the Stokes problem with the Crouzeix-Raviart type nonconforming linear triangular finite element. By the L^2 projection method, the superconvergence results for the velocity and pressure are obtained with a penalty parameter larger than that of the classical penalty scheme. The numerical experiments are carried out to confirm the theoretical results.

SUPERCONVERGENCE OF THE FULL-DISCRETE F.E.M. FOR COMPRESSIBLE MISCIBLE DISPLACEMENT: THE FULL TENSOR CASE

An efficient time stepping procedure is proposed to treat the system describing compressible miscible displacement in a porous medium by employing a mixed finite element method to approximate the pressure and the fluid velocity and a standard Galerkin method to approximate the concentration. An extension of the Darcy velocity along Gauss lines is used in the evaluation of the coefficients in the Galerkin procedure for the concentration. These results show that the total algorithm has the superconvergence property of the fluid velocity.

SUPERCONVERGENCE OF GENERALIZED DIFFERENCE METHOD FOR ELLIPTIC BOUNDARY VALUE PROBLEM

Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.

Superconvergence analysis of bi-k-degree rectangular elements for two-dimensional time-dependent Schrodinger equation

Superconvergence has been studied for long, and many different numerical methods have been analyzed. This paper is concerned with the problem of superconvergence for a two-dimensional time-dependent linear Schr?dinger equation with the finite element method. The error estimate and superconvergence property with order O(hk+1)in the H1norm are given by using the elliptic projection operator in the semi-discrete scheme. The global superconvergence is derived by the interpolation post-processing technique. The superconvergence result with order O(hk+1+ τ2) in the H1norm can be obtained in the Crank-Nicolson fully discrete scheme.

Superconvergence Study of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations

In this paper,we apply the Fourier analysis technique to investigate superconvergence properties of the direct disontinuous Galerkin(DDG)method(Liu and Yan in SIAM J Numer Anal 47(1):475-698,2009),the DDG method with the interface correction(DDGIC)(Liu and Yan in Commun Comput Phys 8(3):541-564,2010),the symmetric DDG method(Vidden and Yan in Comput Math 31(6):638-662,2013),and the nonsymmetric DDG method(Yan in J Sci Comput 54(2):663-683,2013).We also include the study of the interior penalty DG(IPDG)method,due to its close relation to DDG methods.Error estimates are carried out for both P2 and P3 polynomial approximations.By investigating the quantitative errors at the Lobatto points,we show that the DDGIC and symmetric DDG methods are superior,in the sense of obtaining(k+2)th superconvergence orders for both P2 and P3 approximations.Superconvergence order of(k+2)is also observed for the IPDG method with P3 polynomial approximations.The errors are sensitive to the choice of the numerical flux coefficient for even degree P2 approximations,but are not for odd degree P3 approxi-mations.Numerical experiments are carried out at the same time and the numerical errors match well with the analytically estimated errors.

Superconvergence analysis of lower order anisotropic finite element

The convergence analysis of the lower order nonconforming element pro- posed by Park and Sheen is applied to the second-order elliptic problem under anisotropic meshes. The corresponding error estimation is obtained. Moreover, by using the interpo- lation postprocessing technique, a global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is derived. Numerical results are also given to verify the theoretical analysis.

Superconvergence and the Numerical Flux: a Study Using the Upwind-Biased Flux in Discontinuous Galerkin Methods

One of the beneficial properties of the discontinuous Galerkin method is the accurate wave propagation properties.That is,the semi-discrete error has dissipation errors of order 2k+1(≤Ch2k+1)and order 2k+2 for dispersion(≤Ch2k+2).Previous studies have concentrated on the order of accuracy,and neglected the important role that the error constant,C,plays in these estimates.In this article,we show the important role of the error constant in the dispersion and dissipation error for discontinuous Galerkin approximation of polynomial degree k,where k=0,1,2,3.This gives insight into why one may want a more centred flux for a piecewise constant or quadratic approximation than for a piecewise linear or cubic approximation.We provide an explicit formula for these error constants.This is illustrated through one particular flux,the upwind-biased flux introduced by Meng et al.,as it is a convex combination of the upwind and downwind fluxes.The studies of wave propagation are typically done through a Fourier ansatz.This higher order Fourier information can be extracted using the smoothness-increasing accuracy-conserving(SIAC)filter.The SIAC filter ties the higher order Fourier information to the negative-order norm in physical space.We show that both the proofs of the ability of the SIAC filter to extract extra accuracy and numerical results are unaffected by the choice of flux.

Superconvergence analysis of the finite element method for nonlinear hyperbolic equations with nonlinear boundary condition

This paper discusses the semidiscrete finite element method for nonlinear hyperbolic equations with nonlinear boundary condition. The superclose property is derived through interpolation instead of the nonlinear H^1 projection and integral identity technique. Meanwhile, the global superconvergence is obtained based on the interpolated postprocessing techniques.

Superconvergence of Energy-Conserving Discontinuous Galerkin Methods for Linear Hyperbolic Equations

In this paper,we study the superconvergence properties of the energy-conserving discontinuous Galerkin(DG)method in[18]for one-dimensional linear hyperbolic equations.We prove the approximate solution superconverges to a particular projection of the exact solution.The order of this superconvergence is proved to be k+2 when piecewise Pk polynomials with K≥1 are used.The proof is valid for arbitrary non-uniform regular meshes and for piecewise polynomials with arbitrary K≥1.Furthermore,we find that the derivative and function value approxi?mations of the DG solution are superconvergent at a class of special points,with an order of k+1 and R+2,respectively.We also prove,under suitable choice of initial discretization,a(2k+l)-th order superconvergence rate of the DG solution for the numerical fluxes and the cell averages.Numerical experiments are given to demonstrate these theoretical results.

Superconvergence Analysis of the Runge-Kutta Discontinuous Galerkin Method with Upwind-Biased Numerical Flux for Two-Dimensional Linear Hyperbolic Equation

In this paper,we shall establish the superconvergence properties of the Runge-Kutta dis-continuous Galerkin method for solving two-dimensional linear constant hyperbolic equa-tion,where the upwind-biased numerical flux is used.By suitably defining the correction function and deeply understanding the mechanisms when the spatial derivatives and the correction manipulations are carried out along the same or different directions,we obtain the superconvergence results on the node averages,the numerical fluxes,the cell averages,the solution and the spatial derivatives.The superconvergence properties in space are pre-served as the semi-discrete method,and time discretization solely produces an optimal order error in time.Some numerical experiments also are given.

Superconvergence of Finite Element Approximations to Parabolic and Hyperbolic Integro-Differential Equations

The object of this paper is to investigate the superconvergence properties of finite element approximations to parabolic and hyperbolic integro-differential equations. The quasi projection technique introduced earlier by Douglas et al. is developed to derive the O(h2r) order knot superconvergence in the case of a single space variable, and to show the optimal order negative norm estimates in the case of several space variables.

Application for Superconvergence of Finite Element Approximations for the Elliptic Problem by Global and Local L²-Projection Methods

Numerical experiments are given to verify the theoretical results for superconvergence of the elliptic problem by global and local L2-Projection methods.

Convergence and Superconvergence of Fully Discrete Finite Element for Time Fractional Optimal Control Problems

In this paper, we consider a fully discrete finite element approximation for time fractional optimal control problems. The state and adjoint state are approximated by triangular linear fi nite elements in space and L1 scheme in time. The control is obtained by the variational discretization technique. The main purpose of this work is to derive the convergence and superconvergence. A numerical example is presented to validate our theoretical results.

SUPERCONVERGENCE OF RITZ-VOLTERRA PROJECTION ON FINITE ELEMENT SPACES

The purpose of this paper is to study the superconvergence properties of Ritz-Volterra projection.Through construction a new type of Green function and making use of its properties and the principle of duality,the paper proves that the Ritz-Volterra projection defined on r-1 order finite element spaces of Lagrange type in one and two space variable cases possesses O(h2r^2)order and O(h4+1|Inh|)order nodal superconvergence,respectively,and the same type of superconver-gence results are demonstrated for the semidiscrete finite dement approximate solutions of Soboleve-quations.

Anisotropic Biquadratic Finite Element with Some Natural Superconvergence Results

The paper studies the convergence and the superconvergence of the biquadratic finite element for Poisson＇ problem on anisotropic meshes. By detailed analysis, it shows that the biquadratic finite element is anisotropically superconvergent at four Gauss points in the element. Key words：

Approximation and Superconvergence Analysis for a New Higher Order Wilson Element to Sobolev Type Equations

In this paper, a new higher order Wilson element is presented, and the convergence is proved. Then the interpolation postprocessing technique is used to obtain the global superconvergence and posterior error estimate of higher accuracy of this new element for the Sobolev type equations.