Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

L p estimates of boundary integral equations for some nonlinear boundary value problems

Summary Recently, Galerkin and collocation methods have been analyzed for boundary integral equation formulations of some potential problems in the plane with nonlinear boundary conditions, and stability results and error estimates in theH ^1/2-norm have been proved (Ruotsalainen and Wendland, and Ruotsalainen and Saranen). We show that these results extend toL ^p setting without any extra conditions. These extensions are proved by studying the uniform boundedness of the inverses of the linearized integral operators, and then considering the nonlinear equations. The fact that inH ^1/2 setting the nonlinear operator is a homeomorphism with Lipschitz continuous inverse plays a crucial role. Optimal error estimates for the Galerkin and collocation method inL ^p space then follow.This research was performed while the second author was visiting professor at the University of Delaware, spring 1989     

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain are given. These procedures use implicit Galerkin method in the subdomains and simple explicit flux calculation on the interdomain boundaries by integral mean method or extrapolation method to predict the inner‐boundary conditions. Thus, the parallelism can be achieved by these procedures. These procedures are conservative both in the subdomains and across interboundaries. The explicit nature of the flux prediction induces a time‐step limitation that is necessary to preserve stability, but this constraint is less severe than that for a fully explicit method. L²‐norm error estimates are derived for these procedures. Compared with the work of Dawson and Dupont [Math Comp 58 (1992), 21–35], these L²‐norm error estimates avoid the loss of H^?1/2 factor. Experimental results are presented to confirm the theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Variational formulation for the stationary fractional advection dispersion equation

In this article a theoretical framework for the Galerkin finite element approximation to the steady state fractional advection dispersion equation is presented. Appropriate fractional derivative spaces are defined and shown to be equivalent to the usual fractional dimension Sobolev spaces H^s. Existence and uniqueness results are proven, and error estimates for the Galerkin approximation derived. Numerical results are included that confirm the theoretical estimates. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A Petrov–Galerkin spectral method for fourth‐order problems

In this paper, we consider the Petrov–Galerkin spectral method for fourth‐order elliptic problems on rectangular domains subject to non‐homogeneous Dirichlet boundary conditions. We derive some sharp results on the orthogonal approximations in one and two dimensions, which play important roles in numerical solutions of higher‐order problems. By applying these results to a fourth‐order problem, we establish the H²‐error and L²‐error bounds of the Petrov–Galerkin spectral method. Numerical experiments are provided to illustrate the high accuracy of the proposed method and coincide well with the theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Superconvergence analysis of a linearized three‐step backward differential formula finite element method for nonlinear Sobolev equation

A linearized three‐step backward differential formula (BDF) Galerkin finite element method (FEM) is developed for nonlinear Sobolev equation with bilinear element. Temporal error and spatial error are discussed through introducing a time‐discrete system. Solutions of the time‐discrete system are bounded in H²‐norm by the temporal error. Superconvergence results of order O(h² + τ³) in H¹‐norm for the original variable are deduced based on the spatial error. Some new tricks are utilized to get higher order of the temporal error and the spatial error. At last, two numerical examples are provided to support the theoretical analysis. Here, h is the subdivision parameter, and τ is the time step.     

A posteriori error estimates of constrained optimal control problem governed by convection diffusion equations

In this paper, we study a posteriori error estimates of the edge stabilization Galerkin method for the constrained optimal control problem governed by convection-dominated diffusion equations. The residual-type a posteriori error estimators yield both upper and lower bounds for control u measured in L ²-norm and for state y and costate p measured in energy norm. Two numerical examples are presented to illustrate the effectiveness of the error estimators provided in this paper.      

Characteristics—Galerkin and mixed finite element approximation of contamination by compressible nuclear waste-disposal in porous media

A compressible nuclear waste-disposal contamination in porous media is modeled by a coupled system of partial differential equations. The approximation of this system by a Galerkin method that makes use of a modified method of characteristics for the brine, radionuclides, and heat and by a mixed finite element method for the pressure and velocity are analyzed. Optimal order error estimates are obtained. This article improves upon previous works in two aspects. First, error analysis is given with no restriction on the diffusion-dispersion tensors. That is, we have included the effects of molecular diffusion and dispersion. Secondly, optimal error estimates in H¹ and L² are derived. © 1996 John Wiley & Sons, Inc.     

Interior penalty discontinuous Galerkin methods with implicit time‐integration techniques for nonlinear parabolic equations

We prove existence and numerical stability of numerical solutions of three fully discrete interior penalty discontinuous Galerkin methods for solving nonlinear parabolic equations. Under some appropriate regularity conditions, we give the l²(H¹) and l^∞(L²) error estimates of the fully discrete symmetric interior penalty discontinuous Galerkin–scheme with the implicit θ ‐schemes in time, which include backward Euler and Crank–Nicolson finite difference approximations. Our estimates are optimal with respect to the mesh size h. The theoretical results are confirmed by some numerical experiments. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Weak Galerkin finite element methods for Parabolic equations

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal‐order error estimates in both H¹ and L² norms are established. Numerical tests are performed and reported. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Higher order convergence results for the Rayleigh-Ritz method applied to eigenvalue problems: 2. Improved error bounds for eigenfunctions

The application of the Rayleigh-Ritz method for approximating the eigenvalues and eigenfunctions of linear eigenvalue problems in several dimensions is investigated. The object is to improve upon known error estimates for the approximate eigenfunctions. Results for the Galerkin approximation of the eigenfunctions are developed under varying assumptions on the boundary conditions and domain of definition of the eigenvalue problem. These results, coupled with a previous result relating Galerkin and Rayleigh-Ritz approximation of the eigenfunctions, are then used to obtain improved error estimates for the approximate eigenfunctions in theL ² and uniform norms.This research was supported in part by AEC Grant (11-1)-2075.     

A posteriori error control for finite element approximations of elliptic eigenvalue problems

We develop a new approach to a posteriori error estimation for Galerkin finite element approximations of symmetric and nonsymmetric elliptic eigenvalue problems. The idea is to embed the eigenvalue approximation into the general framework of Galerkin methods for nonlinear variational equations. In this context residual-based a posteriori error representations are available with explicitly given remainder terms. The careful evaluation of these error representations for the concrete situation of an eigenvalue problem results in a posteriori error estimates for the approximations of eigenvalues as well as eigenfunctions. These suggest local error indicators that are used in the mesh refinement process.     

Optimal L ∞(L 2)-Error Estimates for the DG Method Applied to Nonlinear Convection–Diffusion Problems with Nonlinear Diffusion

This article is concerned with the analysis of the discontinuous Galerkin finite element method (DGFEM) applied to the space semidiscretization of a nonstationary convection–diffusion problem with nonlinear convection and nonlinear diffusion. Optimal estimates in the L ^∞(L ²)-norm are derived for the symmetric interior penalty (SIPG) scheme in two dimensions. The error analysis is carried out for nonconforming triangular meshes under the assumption that the exact solution of the problem and the solution of a linearized elliptic dual problem are sufficiently regular.     

Multiscale discontinuous Petrov–Galerkin method for the multiscale elliptic problems

In this article, we present a new multiscale discontinuous Petrov–Galerkin method (MsDPGM) for multiscale elliptic problems. This method utilizes the classical oversampling multiscale basis in the framework of a Petrov–Galerkin version of the discontinuous Galerkin method, allowing us to better cope with multiscale features in the solution. MsDPGM takes advantage of the multiscale Petrov–Galerkin method (MsPGM) and the discontinuous Galerkin method (DGM). It can eliminate the resonance error completely and decrease the computational costs of assembling the stiffness matrix, thus, allowing for more efficient solution algorithms. On the basis of a new H² norm error estimate between the multiscale solution and the homogenized solution with the first‐order corrector, we give a detailed convergence analysis of the MsDPGM under the assumption of periodic oscillating coefficients. We also investigate a multiscale discontinuous Galerkin method (MsDGM) whose bilinear form is the same as that of the DGM but the approximation space is constructed from the classical oversampling multiscale basis functions. This method has not been analyzed theoretically or numerically in the literature yet. Numerical experiments are carried out on the multiscale elliptic problems with periodic and randomly generated log‐normal coefficients. Their results demonstrate the efficiency of the proposed method.     

One‐level and two‐level finite element methods for the Boussinesq problem

Previous works on the convergence of numerical methods for the Boussinesq problem were conducted, while the optimal L²‐norm error estimates for the velocity and temperature are still lacked. In this paper, the backward Euler scheme is used to discrete the time terms, standard Galerkin finite element method is adopted to approximate the variables. The MINI element is used to approximate the velocity and pressure, the temperature field is simulated by the linear polynomial. Under some restriction on the time step, we firstly present the optimal L² error estimates of approximate solutions. Secondly, two‐level method based on Stokes iteration for the Boussinesq problem is developed and the corresponding convergence results are presented. By this method, the original problem is decoupled into two small linear subproblems. Compared with the standard Galerkin method, the two‐level method not only keeps good accuracy but also saves a lot of computational cost. Finally, some numerical examples are provided to support the established theoretical analysis.     

Stability and error analysis for a spectral Galerkin method for the Navier‐Stokes equations with H2 or H1 initial data

In this article we consider a spectral Galerkin method with a semi‐implicit Euler scheme for the two‐dimensional Navier‐Stokes equations with H² or H¹ initial data. The H²‐stability analysis of this spectral Galerkin method shows that for the smooth initial data the semi‐implicit Euler scheme admits a large time step. The L²‐error analysis of the spectral Galerkin method shows that for the smoother initial data the numerical solution u exhibits faster convergence on the time interval [0, 1] and retains the same convergence rate on the time interval [1, ∞). © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

A CFL‐free explicit characteristic interior penalty scheme for linear advection‐reaction equations

We develop a CFL‐free, explicit characteristic interior penalty scheme (CHIPS) for one‐dimensional first‐order advection‐reaction equations by combining a Eulerian‐Lagrangian approach with a discontinuous Galerkin framework. The CHIPS method retains the numerical advantages of the discontinuous Galerkin methods as well as characteristic methods. An optimal‐order error estimate in the L² norm for the CHIPS method is derived and numerical experiments are presented to confirm the theoretical estimates. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Conform and mixed finite element schemes for the Dirichlet problem for the Bilaplacian in plane domains with corners

In plane domains with corners for the Bilaplacian a uniquely solvable conform variational principle is studied on weighted Sobolev spaces which is equivalent to the standard Dirichlet problem in the weak form. Clamped plates under point forces near corners are handled by this approach. With weighted Hsieh-Clough-Tocher elements on regular triangulations as conform C¹-finite elements a new error analysis is performed without higher regularity assumptions on the exact solution than given by the data and the boundary. The rate of convergence of the error depends on the eigenvalue with smallest imaginary part of a clamped infinite wedge since this eigenvalue describes the singularity of the exact solution in a sector with same angle. Using different spaces of trial and test functions in the standard Galerkin procedure it is shown that the error in the weighted energy norm does not pollute. For convex corners asymptotic error estimates, are proved yielding convergence for a mixed method in hydrodynamics where the solution of a system of 2^nd order and its Laplacian are approximated simultaneously by C⁰-finite elements being piecewise polynomials.