Gauss‐Lobatto‐Legendre‐Birkhoff pseudospectral approximations for the multi‐term time fractional diffusion‐wave equation with Neumann boundaryconditions

A second‐order finite difference/pseudospectral scheme is proposed for numerical approximation of multi‐term time fractional diffusion‐wave equation with Neumann boundary conditions. The scheme is based upon the weighted and shifted Grünwald difference operators approximation of the time fractional calculus and Gauss‐Lobatto‐Legendre‐Birkhoff (GLLB) pseudospectral method for spatial discretization. The unconditionally stability and convergence of the scheme are rigorously proved. Numerical examples are carried out to verify theoretical results.     

Numerical experiments on a hybrid WENO5 filter for shock‐capturing

In the present paper, a hybrid filter is introduced for high accurate numerical simulation of shock‐containing flows. The fourth‐order compact finite difference scheme is used for the spatial discretization and the third‐order Runge–Kutta scheme is used for the time integration. After each time‐step, the hybrid filter is applied on the results. The filter is composed of a linear sixth‐order filter and the dissipative part of a fifth‐order weighted essentially nonoscillatory scheme (WENO5). The classic WENO5 scheme and the WENO5 scheme with adaptive order (WENO5‐AO) are used to form the hybrid filter. Using a shock‐detecting sensor, the hybrid filter reduces to the linear sixth‐order filter in smooth regions for damping high frequency waves and reduces to the WENO5 filter at shocks in order to eliminate unwanted oscillations produced by the nondissipative spatial discretization method. The filter performance and accuracy of the results are examined through several test cases including the advection, Euler and Navier–Stokes equations. The results are compared with that of a hybrid second‐order filter and also that of the WENO5 and WENO5‐AO schemes.     

Stabilized semi‐implicit spectral deferred correction methods for Allen–Cahn and Cahn–Hilliard equations

Stabilized semi‐implicit spectral deferred correction methods are constructed for the time discretization of Allen–Cahn and Cahn–Hilliard equations. These methods are unconditionally stable, lead to simple linear system to solve at each iteration, and can achieve high‐order time accuracy with a few iterations in each time step. Ample numerical results are presented to demonstrate the effectiveness of the stabilized semi‐implicit spectral deferred correction methods for solving the Allen–Cahn and Cahn–Hilliard equations. Copyright © 2013 John Wiley & Sons, Ltd.     

High order compact solution of the one‐space‐dimensional linear hyperbolic equation

In this article, we introduce a high‐order accurate method for solving one‐space dimensional linear hyperbolic equation. We apply a compact finite difference approximation of fourth order for discretizing spatial derivative of linear hyperbolic equation and collocation method for the time component. The main property of this method additional to its high‐order accuracy due to the fourth order discretization of spatial derivative, is its unconditionally stability. In this technique the solution is approximated by a polynomial at each grid point that its coefficients are determined by solving a linear system of equations. Numerical results show that the compact finite difference approximation of fourth order and collocation method produce a very efficient method for solving the one‐space‐dimensional linear hyperbolic equation. We compare the numerical results of this paper with numerical results of (Mohanty, 3 .© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Maximum norm error analysis of an unconditionally stable semi‐implicit scheme for multi‐dimensional Allen–Cahn equations

In this paper, a linearized finite difference scheme is proposed for solving the multi‐dimensional Allen–Cahn equation. In the scheme, a modified leap‐frog scheme is used for the time discretization, the nonlinear term is treated in a semi‐implicit way, and the central difference scheme is used for the discretization in space. The proposed method satisfies the discrete energy decay property and is unconditionally stable. Moreover, a maximum norm error analysis is carried out in a rigorous way to show that the method is second‐order accurate both in time and space variables. Finally, numerical tests for both two‐ and three‐dimensional problems are provided to confirm our theoretical findings.     

A fast singly diagonally implicit runge–kutta method for solving 1D unsteady convection‐diffusion equations

In this article, a fast singly diagonally implicit Runge–Kutta method is designed to solve unsteady one‐dimensional convection diffusion equations. We use a three point compact finite difference approximation for the spatial discretization and also a three‐stage singly diagonally implicit Runge–Kutta (RK) method for the temporal discretization. In particular, a formulation evaluating the boundary values assigned to the internal stages for the RK method is derived so that a phenomenon of the order of the reduction for the convergence does not occur. The proposed scheme not only has fourth‐order accuracy in both space and time variables but also is computationally efficient, requiring only a linear matrix solver for a tridiagonal matrix system. It is also shown that the proposed scheme is unconditionally stable and suitable for stiff problems. Several numerical examples are solved by the new scheme and the numerical efficiency and superiority of it are compared with the numerical results obtained by other methods in the literature. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 788–812, 2014     

Singly diagonally implicit runge‐kutta method for time‐dependent reaction‐diffusion equation

In this article, an efficient fourth‐order accurate numerical method based on Padé approximation in space and singly diagonally implicit Runge‐Kutta method in time is proposed to solve the time‐dependent one‐dimensional reaction‐diffusion equation. In this scheme, we first approximate the spatial derivative using the second‐order central finite difference then improve it to fourth‐order by applying Padé approximation. A three stage fourth‐order singly diagonally implicit Runge‐Kutta method is then used to solve the resulting system of ordinary differential equations. It is also shown that the scheme is unconditionally stable, and is suitable for stiff problems. Several numerical examples are solved by the scheme and the efficiency and accuracy of the new scheme are compared with two widely used high‐order compact finite difference methods. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1423–1441, 2011     

Instability of the finite‐difference split‐step method applied to the nonlinear Schrödinger equation. I. standing soliton

We consider numerical instability that can be observed in simulations of solitons of the nonlinear Schrödinger equation (NLS) by a split‐step method (SSM) where the linear part of the evolution is solved by a finite‐difference discretization. The von Neumann analysis predicts that this method is unconditionally stable on the background of a constant‐amplitude plane wave. However, simulations show that the method can become unstable on the background of a soliton. We present an analysis explaining this instability. Both this analysis and the features and threshold of the instability are substantially different from those of the Fourier SSM, which computes the linear part of the NLS by a spectral discretization. For example, the modes responsible for the numerical instability are not similar to plane waves, as for the Fourier SSM or, more generally, in the von Neumann analysis. Instead, they are localized at the sides of the soliton. This also makes them different from “physical” (as opposed to numerical) unstable modes of nonlinear waves, which (the modes) are localized around the “core” of a solitary wave. Moreover, the instability threshold for thefinite‐difference split‐step method is considerably relaxed compared with that for the Fourier split‐step. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1002–1023, 2016     

Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

An efficient numerical method for the resolution of the Kirchhoff‐Love dynamic plate equation

We solve numerically the Kirchhoff‐Love dynamic plate equation for an anisotropic heterogeneous material using a spectral method. A mixed velocity‐moment formulation is proposed for the space approximation allowing the use of classical Lagrange finite elements. The benefit of using high order elements is shown through a numerical dispersion analysis. The system resulting from this spatial discretization is solved analytically. Hence this method is particularly efficient for long duration experiments. This time evolution method is compared with explicit and implicit finite differences schemes in terms of accuracy and computation time. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Effect of numerical integration in the DGFEM for nonlinear convection‐diffusion problems

This paper is concerned with the effect of numerical integration applied to the discontinuous Galerkin finite element discretization of nonlinear convection‐diffusion problems in 2D. In the space semidiscretization the volume and line integrals are evaluated by numerical quadratures. Our goal is to estimate the error caused by the numerical integration and to show what numerical quadratures guarantee that the accuracy of the method with exact integration is preserved. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Minimum time‐step criteria for the Galerkin finite element methods applied to one‐dimensional parabolic partial differential equations

The finite element method has been well established for numerically solving parabolic partial differential equations (PDEs). Also it is well known that a too large time step should not be chosen in order to obtain a stable and accurate numerical solution. In this article, accuracy analysis shows that a too small time step should not be chosen either for some time‐stepping schemes. Otherwise, the accuracy of the numerical solution cannot be improved or can even be worsened in some cases. Furthermore, the so‐called minimum time step criteria are established for the Crank‐Nicolson scheme, the Galerkin‐time scheme, and the backward‐difference scheme used in the temporal discretization. For the forward‐difference scheme, no minimum time step exists as far as the accuracy is concerned. In the accuracy analysis, no specific initial and boundary conditions are invoked so that such established criteria can be applied to the parabolic PDEs subject to any initial and boundary conditions. These minimum time step criteria are verified in a series of numerical experiments for a one‐dimensional transient field problem with a known analytical solution. The minimum time step criteria developed in this study are useful for choosing appropriate time steps in numerical simulations of practical engineering problems. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A high‐order compact ADI method for solving three‐dimensional unsteady convection‐diffusion problems

We derive a high‐order compact alternating direction implicit (ADI) method for solving three‐dimentional unsteady convection‐diffusion problems. The method is fourth‐order in space and second‐order in time. It permits multiple uses of the one‐dimensional tridiagonal algorithm with a considerable saving in computing time and results in a very efficient solver. It is shown through a discrete Fourier analysis that the method is unconditionally stable in the diffusion case. Numerical experiments are conducted to test its high order and to compare it with the standard second‐order Douglas‐Gunn ADI method and the spatial fourth‐order compact scheme by Karaa. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Spectral analysis of a high-order Hermitian algorithm for structural dynamics

The spectral analysis of an efficient step-by-step direct integration algorithm for the structural dynamic equation is presented. The proposed algorithm is formulated in terms of two Hermitian finite difference operators of fifth-order local truncation error and it is unconditionally stable with no numerical damping presenting a fourth-order truncation error for period dispersion (global error). In addition, although it is in competition with higher-order algorithms presented in the literature, the computational effort is similar to that of the classical second-order Newmark’s method. The numerical application for nonlinear structural dynamic problems is also considered.     

High‐order numerical solution of second‐order one‐dimensional hyperbolic telegraph equation using a shifted Gegenbauer pseudospectral method

We present a high‐order shifted Gegenbauer pseudospectral method (SGPM) to solve numerically the second‐order one‐dimensional hyperbolic telegraph equation provided with some initial and Dirichlet boundary conditions. The framework of the numerical scheme involves the recast of the problem into its integral formulation followed by its discretization into a system of well‐conditioned linear algebraic equations. The integral operators are numerically approximated using some novel shifted Gegenbauer operational matrices of integration. We derive the error formula of the associated numerical quadratures. We also present a method to optimize the constructed operational matrix of integration by minimizing the associated quadrature error in some optimality sense. We study the error bounds and convergence of the optimal shifted Gegenbauer operational matrix of integration. Moreover, we construct the relation between the operational matrices of integration of the shifted Gegenbauer polynomials and standard Gegenbauer polynomials. We derive the global collocation matrix of the SGPM, and construct an efficient computational algorithm for the solution of the collocation equations. We present a study on the computational cost of the developed computational algorithm, and a rigorous convergence and error analysis of the introduced method. Four numerical test examples have been carried out to verify the effectiveness, the accuracy, and the exponential convergence of the method. The SGPM is a robust technique, which can be extended to solve a wide range of problems arising in numerous applications. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 307–349, 2016     

Convergence of Extrapolated BDF2 Finite Element Schemes for Unsteady Penetrative Convection Model

Extrapolated two-step backward difference (BDF2) in time and finite element in space discretization for the unsteady penetrative convection model is analyzed. Penetrative convection model employs a nonlinear equation of state making the problem more nonlinear. Optimal order error estimates are derived for the semi-discrete finite element spatial discretization. Two time discretization schemes based on linear extrapolation are proposed and analyzed, namely a coupled and a decoupled scheme. In particular, we show that although both schemes are unconditionally nonlinearly stable, the decoupled scheme converges unconditionally whereas coupled scheme requires that the time step be sufficiently small for convergence. These time discretization schemes can be implemented efficiently in practice, saving computational memory. Numerical computations and numerical convergence checks are presented to demonstrate the efficiency and the accuracy of the schemes.     

Space‐time spectral method for two‐dimensional semilinear parabolic equations

In this paper, a high‐order accurate numerical method for two‐dimensional semilinear parabolic equations is presented. We apply a Galerkin–Legendre spectral method for discretizing spatial derivatives and a spectral collocation method for the time integration of the resulting nonlinear system of ordinary differential equations. Our formulation can be made arbitrarily high‐order accurate in both space and time. Optimal a priori error bound is derived in the L²‐norm for the semidiscrete formulation. Extensive numerical results are presented to demonstrate the convergence property of the method, show our formulation have spectrally accurate in both space and time. John Wiley & Sons, Ltd.     

An IMEX‐BDF2 compact scheme for pricing options under regime‐switching jump‐diffusion models

In this paper, an implicit‐explicit two‐step backward differentiation formula (IMEX‐BDF2) together with finite difference compact scheme is developed for the numerical pricing of European and American options whose asset price dynamics follow the regime‐switching jump‐diffusion process. It is shown that IMEX‐BDF2 method for solving this system of coupled partial integro‐differential equations is stable with the second‐order accuracy in time. On the basis of IMEX‐BDF2 time semi‐discrete method, we derive a fourth‐order compact (FOC) finite difference scheme for spatial discretization. Since the payoff function of the option at the strike price is not differentiable, the results show only second‐order accuracy in space. To remedy this, a local mesh refinement strategy is used near the strike price so that the accuracy achieves fourth order. Numerical results illustrate the effectiveness of the proposed method for European and American options under regime‐switching jump‐diffusion models.     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Semi‐implicit schemes for transient Navier–Stokes equations and eddy viscosity models

This study presents two computational schemes for the numerical approximation of solutions to eddy viscosity models as well as transient Navier–Stokes equations. The eddy viscosity model is one example of a class of Large Eddy Simulation models, which are used to simulate turbulent flow. The first approximation scheme is a first order single step method that treats the nonlinear term using a semi‐implicit discretization. The second scheme employs a two step approach that applies a Crank–Nicolson method for the nonlinear term while also retaining the semi‐implicit treatment used in the first scheme. A finite element approximation is used in the spatial discretization of the partial differential equations. The convergence analysis for both schemes is discussed in detail, and numerical results are given for two test problems one of which is the two dimensional flow around a cylinder. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009