Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

A numerical method for the one‐dimensional sine‐Gordon equation

A numerical method based on a predictor–corrector (P‐C) scheme arising from the use of rational approximants of order 3 to the matrix‐exponential term in a three‐time level recurrence relation is applied successfully to the one‐dimensional sine‐Gordon equation, already known from the bibliography. In this P‐C scheme a modification in the corrector (MPC) has been proposed according to which the already evaluated corrected values are considered. The method, which uses as predictor an explicit finite‐difference scheme arising from the second order rational approximant and as corrector an implicit one, has been tested numerically on the single and the soliton doublets. Both the predictor and the corrector schemes are analyzed for local truncation error and stability. From the investigation of the numerical results and the comparison of them with other ones known from the bibliography it has been derived that the proposed P‐C/MPC schemes at least coincide in terms of accuracy with them. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A space‐time mixed‐hybrid finite element method for the damped wave equation

A space‐time finite element method is introduced to solve the linear damped wave equation. The scheme is constructed in the framework of the mixed‐hybrid finite element methods, and where an original conforming approximation of H(div;Ω) is used, the latter permits us to obtain an upwind scheme in time. We establish the link between the nonstandard finite difference scheme recently introduced by Mickens and Jordan and the scheme proposed. In this regard, two approaches are considered and in particular we employ a formulation allowing the solution to be marched in time, i.e., one only needs to consider one time increment at a time. Numerical results are presented and compared with the analytical solution illustrating good performance of the present method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Finite analytic numerical method for solving two‐dimensional quasi‐Laplace equation

A typical power series analytic solution of quasi‐Laplace equation in the infinitesimal angle domain around the singular point of the square cells is provided in this article. Toward the singular point, the gradient of the potential variable will tend to infinity, which is described by the first term of the power series solution. Based on this analytic solution, three finite analytic numerical methods are proposed. These methods are analogous and are constructed, respectively, when considering different numbers of the terms or using different schemes to determine the relevant parameters in the power series. Numerical examples show that all of the three finite analytic numerical methods proposed can provide rather accurate solutions than the traditional numerical methods. In contrast, when using the traditional numerical schemes to solve the quasi‐Laplace equation in a strong heterogeneous medium, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result. In practical applications, subdividing each origin cell into 2 × 2 or 3 × 3 subcells is enough for the finite analytical numerical methods to get relatively accurate results. The finite analytical numerical methods are also convenient to construct the flux field with high accuracy.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1755–1769, 2014     

Influence of Gauss and Gauss‐Lobatto quadrature rules on the accuracy of a quadrilateral finite element method in the time domain

In this article, we examine the influence of numerical integration on finite element methods using quadrilateral or hexahedral meshes in the time domain. We pay special attention to the use of Gauss‐Lobatto points to perform mass lumping for any element order. We provide some theoretical results through several error estimates that are completed by various numerical experiments. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A numerical method based on extended Raviart–Thomas (ER‐T) mixed finite element method for solving damped Boussinesq equation

In this paper, we present the approximate solution of damped Boussinesq equation using extended Raviart–Thomas mixed finite element method. In this method, the numerical solution of this equation is obtained using triangular meshes. Also, for discretization in time direction, we use an implicit finite difference scheme. In addition, error estimation and stability analysis of both methods are shown. Finally, some numerical examples are considered to confirm the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

An equation decomposition method for the numerical solution of a fourth‐order elliptic singular perturbation problem

In this article, we propose a tailored finite point method (TFPM) for the numerical solution of a type of fourth‐order singular perturbation problem in two dimensions based on the equation decomposition technique. Our finite point method has been tailored based on the local exponential basis functions. Furthermore, our TFPM satisfies the discrete maximum principle automatically. Our numerical examples show that our method has second order convergence rate in energy norm as $\varepsilon\to0$ . © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

On the numerical solution of the Klein‐Gordon equation

A predictor–corrector (P–C) scheme based on the use of rational approximants of second‐order to the matrix‐exponential term in a three‐time level reccurence relation is applied to the nonlinear Klein‐Gordon equation. This scheme is accelerated by using a modification (MPC) in which the already evaluated values are used for the corrector. Both the predictor and the corrector scheme are analyzed for local truncation error and stability. The proposed method is applied to problems possessing periodic, kinks and single, double‐soliton waves. The accuracy as well as the long time behavior of the proposed scheme is discussed. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A numerical method for spatial diffusion in age‐structured populations

We propose a method to approximate numerically the diffusion of an age‐structured population in a spatial environment. We integrate separately the age and time variables by finite differences and we discretize the space variable by finite elements. The method is implicit in time and, inside each time step, implicit in age. We provide stability and convergence results and we illustrate our approach with some numerical result. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Solution of some plate bending problems using the boundary element method

Some fundamental aspects of the boundary element method of the Kirchhoff theory of thin plate flexure are given. The direct boundary integral equation method with higher conforming properties (using first-order Hermitian interpolation for plate displacement ω, and zero-order Hermitian interpolation for angle of rotation θ, the moment M andthe equivalent shear V) are used for several computational examples. They are: square plate with simply-supported or clamped edges, the same square plate with square central opening and the cantilevered triangular plates. The results of computation as compared with some known experimental and theoritical results showed that the numerical schemes seemed to be satisfactory for the practical applications.     

Space‐time spectral collocation method for the one‐dimensional sine‐Gordon equation

In this article, we introduce a new space‐time spectral collocation method for solving the one‐dimensional sine‐Gordon equation. We apply a spectral collocation method for discretizing spatial derivatives, and then use the spectral collocation method for the time integration of the resulting nonlinear second‐order system of ordinary differential equations (ODE). Our formulation has high‐order accurate in both space and time. Optimal a priori error bounds are derived in the L²‐norm for the semidiscrete formulation. Numerical experiments show that our formulation have exponential rates of convergence in both space and time. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 670–690, 2015     

An optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows

This article is concerned about an optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows. Using the method, an classical domain decomposition problem is transformed into a constrained minimization problem for which the objective functional is chosen to measure the jump in the dependent variables across the common interfaces between subdomains. The Lagrange multiplier rule is used to transform the constrained optimization problem into an unconstrained one and that rule is applied to derive an optimality system from which optimal solutions may be obtained. The optimality system is also derived using “sensitivity” derivatives instead of the Lagrange multiplier rule. We consider a gradient‐type approach to the solution of domain decomposition problem. The results of some numerical experiments are presented to demonstrate the feasibility and applicability of the algorithm developed in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A pseudo‐spectral method that uses an overlapping multidomain technique for the numerical solution of sine‐Gordon equation in one and two spatial dimensions

In this article, we study an explicit scheme for the solution of sine‐Gordon equation when the space discretization is carried out by an overlapping multidomain pseudo‐spectral technique. By using differentiation matrices, the equation is reduced to a nonlinear system of ordinary differential equations in time that can be discretized with the explicit fourth‐order Runge–Kutta method. To achieve approximation with high accuracy in large domains, the number of space grid points must be large enough. This yields very large and full matrices in the pseudo‐spectral method that causes large memory requirements. The domain decomposition approach provides sparsity in the matrices obtained after the discretization, and this property reduces storage for large matrices and provides economical ways of performing matrix–vector multiplications. Therefore, we propose a multidomain pseudo‐spectral method for the numerical simulation of the sine‐Gordon equation in large domains. Test examples are given to demonstrate the accuracy and capability of the proposed method. Numerical experiments show that the multidomain scheme has an excellent long‐time numerical behavior for the sine‐Gordon equation in one and two dimensions. Copyright © 2013 John Wiley & Sons, Ltd.     

Some existence results for an elliptic equation of Kirchhoff‐type with changing sign data and a logarithmic nonlinearity

The paper deals with some existence results for an elliptic equation of Kirchhoff‐type with changing sign data and a logarithmic nonlinearity by direct variational method, Galerkin approach, and sub‐super solutions method. Our results are natural extension of Boulaaras' work in (Math Methods Appl Sci; 41(13):5203‐5210).     

Error analysis for solving the Korteweg‐de Vries equation by a Legendre pseudo‐spectral method

A Legendre pseudo‐spectral method is proposed for the Korteweg‐de Vries equation with nonperiodic boundary conditions. Appropriate base functions are chosen to get an efficient algorithm. Error analysis is given for both semi‐discrete and fully discrete schemes. The numerical results confirm to the theoretical analysis. © (2000) John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 513–534, (2000)     

A mixed finite element method for a quasi‐Newtonian fluid flow

We propose a mixed formulation for quasi‐Newtonian fluid flow obeying the power law where the stress tensor is introduced as a new variable. Based on such a formulation, a mixed finite element is constructed and analyzed. This finite element method possesses local (i.e., at element level) conservation properties (conservation of the momentum and the mass) as in the finite volume methods. We give existence and uniqueness results for the continuous problem and its approximation and we prove error bounds. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

An efficient pseudo‐spectral Legendre–Galerkin method for solving a nonlinear partial integro‐differential equation arising in population dynamics

The pseudo‐spectral Legendre–Galerkin method (PS‐LGM) is applied to solve a nonlinear partial integro‐differential equation arising in population dynamics. This equation is a competition model in which similar individuals are competing for the same resources. It is a kind of reaction–diffusion equation with integral term corresponding to nonlocal consumption of resources. The proposed method is based on the Legendre–Galerkin formulation for the linear terms and interpolation operator at the Chebyshev–Gauss–Lobatto (CGL) points for the nonlinear terms. Also, the integral term, which is a kind of convolution, is directly computed by a fast and accurate method based on CGL interpolation operator, and thus, the use of any quadrature formula in its computation is avoided. The main difference of the PS‐LGM presented in the current paper with the classic LGM is in treating the nonlinear terms and imposing boundary conditions. Indeed, in the PS‐LGM, the nonlinear terms are efficiently handled using the CGL points, and also the boundary conditions are imposed strongly as collocation methods. Combination of the PS‐LGM with a semi‐implicit time integration method such as second‐order backward differentiation formula and Adams‐Bashforth method leads to reducing the complexity of computations and obtaining a linear algebraic system of equations with banded coefficient matrix. The desired equation is considered on one and two‐dimensional spatial domains. Efficiency, accuracy, and convergence of the proposed method are demonstrated numerically in both cases. Copyright © 2012 John Wiley & Sons, Ltd.     

An accurate numerical method for solving the linear fractional Klein–Gordon equation

In this article, an implementation of an efficient numerical method for solving the linear fractional Klein–Gordon equation (LFKGE) is introduced. The fractional derivative is described in the Caputo sense. The method is based upon a combination between the properties of the Chebyshev approximations and finite difference method (FDM). The proposed method reduces LFKGE to a system of ODEs, which is solved using FDM. Special attention is given to study the convergence analysis and deduce an error upper bound of the proposed method. Numerical example is given to show the validity and the accuracy of the proposed algorithm. Copyright © 2013 John Wiley & Sons, Ltd.     

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 .     

An iterative method based on equation decomposition for the fourth‐order singular perturbation problem

In this article, we propose an iterative method based on the equation decomposition technique ( 1 ) for the numerical solution of a singular perturbation problem of fourth‐order elliptic equation. At each step of the given method, we only need to solve a boundary value problem of second‐order elliptic equation and a second‐order singular perturbation problem. We prove that our approximate solution converges to the exact solution when the domain is a disc. Our numerical examples show the efficiency and accuracy of our method. Our iterative method works very well for singular perturbation problems, that is, the case of 0 < ε ? 1, and the convergence rate is very fast. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013