Stabilized Crouzeix‐Raviart element for the Darcy‐Stokes problem

We stabilize the nonconforming Crouzeix‐Raviart element for the Darcy‐Stokes problem with terms motivated by a discontinuous Galerkin approach. Convergence of the method is shown, also in the limit of vanishing viscosity. Finally, some numerical examples verifying the theoretical predictions are presented. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 21, 2005.     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Notice of retraction and apology

It has come to the attention of the editors and publisher that an article published in Numerical Methods and Partial Differential Equations, Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations, by Mohamed Bensaada, Driss Esselaoui, and Pierre Saramito, Numer Methods Partial Differential Eq 21(6) (2005), 1099–1121, Numerical Methods for Partial Differential Equations Numer Methods Partial Differential Eq(2005)21(6)1099 included large portions that were copied from the following paper without proper citation: Convergence and nonlinear stability of the Lagrange‐Galerkin method for the Navier‐Stokes equations, Endre Suli, Numerische Mathematik, Vol. 53, No. 4, pp. 459–486 (July, 1988). We have retracted the paper and apologize to Dr. Suli.     

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.     

Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations (retracted article)

It has come to the attention of the editors and publisher that an article published in Numerical Methods and Partial Differential Equations, “Second‐order Galerkin‐Lagrange method for the Navier‐Stokes equations,” by Mohamed Bensaada, Driss Esselaoui, and Pierre Saramito, Numer Methods Partial Differential Eq 21(6) (2005), 1099–1121 included large portions that were copied from the following paper without proper citation: “Convergence and nonlinear stability of the Lagrange‐Galerkin method for the Navier‐Stokes equations,” Endre Suli, Numerische Mathematik, Vol. 53, No. 4, pp. 459–486 (July, 1988). We have retracted the paper and apologize to Dr. Suli Numer Methods Partial Differential Eq (2007)23(1)211 .     

Generalization of Runge‐Kutta discontinuous Galerkin method to LWR traffic flow model with inhomogeneous road conditions

The discussed model is characterized by changeable lane numbers and free flow velocities that give rise to the spatially varying flux function in conservation equation. Accordingly a new numerical flux and a new limiter for the Runge‐Kutta Discontinuous Galerkin method are considered, which are compared with a natural but simple extension. It is verified that the new generalization is of high‐resolution and has wider stable and convergent ranges. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Multilevel numerical solutions of convection‐dominated diffusion problems by spline wavelets

In this article, we utilize spline wavelets to establish an adaptive multilevel numerical scheme for time‐dependent convection‐dominated diffusion problems within the frameworks of Galerkin formulation and Eulerian‐Lagrangian localized adjoint methods (ELLAM). In particular, we shall use linear Chui‐Quak semi‐orthogonal wavelets, which have explicit expressions and compact supports. Therefore, both the diffusion term and boundary conditions in the convection‐diffusion problems can be readily handled. Strategies for efficiently implementing the scheme are discussed and numerical results are interpreted from the viewpoint of nonlinear approximation. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A cell boundary element method for elliptic problems

An elementary analysis on the cell boundary element (CBEM) was given by Jeon and Sheen. In this article we improve the previous results in various aspects. First of all, stability and convergence analysis on the rectangular grids are established. Moreover, error estimates are improved. Our improved analysis was possible by recasting of the CBEM in a Petrov‐Galerkin setting. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Time accurate fast wavelet‐Taylor Galerkin method for partial differential equations

We introduce the concept of fast wavelet‐Taylor Galerkin methods for the numerical solution of partial differential equations. In wavelet‐Taylor Galerkin method discretization in time is performed before the wavelet based spatial approximation by introducing accurate generalizations of the standard Euler, θ and leap‐frog time‐stepping scheme with the help of Taylor series expansions in the time step. We will present two different time‐accurate wavelet schemes to solve the PDEs. First, numerical schemes taking advantage of the wavelet bases capabilities to compress the operators and sparse representation of functions which are smooth, except for in localized regions, up to any given accuracy are presented. Here numerical experiments deal with advection equation with the spiky solution in one dimension, two dimensions, and nonlinear equation with a shock in solution in two dimensions. Second, our schemes deal with more regular class of problems where wavelets are not efficient procedure for data compression but we can use the good approximation properties of wavelet. Here time‐accurate schemes lead to consistent mass matrix in an explicit time stepping, which can be solved by approximate factorization techniques. Numerical experiment deals with more regular class of problems like heat equation as well as coupled linear system in two dimensions. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Convergence analysis of the streamline diffusion and discontinuous Galerkin methods for the Vlasov‐Fokker‐Planck system

We prove stability estimates and derive optimal convergence rates for the streamline diffusion and discontinuous Galerkin finite element methods for discretization of the multi‐dimensional Vlasov‐Fokker‐Planck system. The focus is on the theoretical aspects, where we deal with construction and convergence analysis of the discretization schemes. Some related special cases are implemented in M. Asadzadeh [Appl Comput Meth 1(2) (2002), 158–175] and M. Asadzadeh and A. Sopasakis [Comput Meth Appl Mech Eng 191(41–42) (2002), 4641–4661]. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Mixed finite volume method for nonlinear elliptic problems

In this article we construct and analyze a mixed finite volume method for second‐order nonlinear elliptic problems employing H(div; Ω)‐conforming approximations for the vector variable and completely discontinuous approximations for the scalar variable. The main attractive feature of our method is that, although the vector variable is H(div; Ω)‐conforming, one can eliminate it in a local manner to obtain a discontinuous Galerkin method for the scalar variable. Optimal error estimates will be established for both vector and scalar variables. We also present a fully discrete version of this method that is more convenient for computational purposes. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

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     

Inf‐sup stability of Petrov‐Galerkin immersed finite element methods for one‐dimensional elliptic interface problems

In this article, we analyze the Petrov‐Galerkin immersed finite element method (PG‐IFEM) when applied to one‐dimensional elliptic interface problems. In the PG‐IFEM (T. Hou, X. Wu and Y. Zhang, Commun. Math. Sci., 2 (2004), 185‐205, and S. Hou and X. Liu, J. Comput. Phys., 202 (2005), 411‐445), the classic immersed finite element (IFE) space was taken as the trial space while the conforming linear finite element space was taken as the test space. We first prove the inf‐sup condition of the PG‐IFEM and then show the optimal error estimate in the energy norm. We also show the optimal estimate of the condition number of the stiffness matrix. The results are extended to two dimensional problems in a special case.     

Spectral approximation for drainage of an Elastico‐viscous liquid and error analysis

A Fourier‐Galerkin spectral method is proposed and used to analyze a system of quasilinear partial differential equations governing the drainage of liquids of the Oldroyd four‐constant type. It is shown that, Fourier‐Galerkin approximations are convergent with spectral accuracy. An efficient and accurate algorithm based on the Fourier‐Galerkin approximations to the system of quasilinear partial differential equations are developed and implemented. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 492–505, 2012     

Analysis and computational method based on quadratic B‐spline FEM for the Rosenau–Burgers equation

L_∞‐error estimates for B‐spline Galerkin finite element solution of the Rosenau–Burgers equation are considered. The semidiscrete B‐spline Galerkin scheme is studied using appropriate projections. For fully discrete B‐spline Galerkin scheme, we consider the Crank–Nicolson method and analyze the corresponding error estimates in time. Numerical experiments are given to demonstrate validity and order of accuracy of the proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 877–895, 2016     

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 posteriori error estimates for the mixed finite element method with Lagrange multipliers

We consider the mixed finite element method with Lagrange multipliers as applied to second‐order elliptic equations in divergence form with mixed boundary conditions. The corresponding Galerkin scheme is defined by using Raviart‐Thomas spaces. We develop a posteriori error analyses yielding a reliable and efficient estimate based on residuals, and a reliable and quasi‐efficient estimate based on local problems, respectively. Several numerical results illustrate the suitability of these a posteriori estimates for the adaptive computation of the discrete solutions. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An Eulerian‐Lagrangian discontinuous Galerkin method for transient advection‐diffusion equations

We develop an Eulerian‐Lagrangian discontinuous Galerkin method for time‐dependent advection‐diffusion equations. The derived scheme has combined advantages of Eulerian‐Lagrangian methods and discontinuous Galerkin methods. The scheme does not contain any undetermined problem‐dependent parameter. An optimal‐order error estimate and superconvergence estimate is derived. Numerical experiments are presented, which verify the theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Stokes' first problem for a Newtonian fluid in a non‐Darcian porous half‐space using a Laguerre–Galerkin method

A Laguerre–Galerkin method is proposed and analysed for the Stokes' first problem of a Newtonian fluid in a non‐Darcian porous half‐space on a semi‐infinite interval. It is well known that Stokes' first problem has a jump discontinuity on boundary which is the main obstacle in numerical methods. By reformulating this equation with suitable functional transforms, it is shown that the Laguerre–Galerkin approximations are convergent on a semi‐infinite interval with spectral accuracy. An efficient and accurate algorithm based on the Laguerre–Galerkin approximations of the transformed equations is developed and implemented. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented. Copyright © 2007 John Wiley & Sons, Ltd.