Approximate solutions of fractional Zakharov–Kuznetsov equations by VIM

This paper presents the approximate analytical solution of a fractional Zakharov–Kuznetsov equation with the help of the powerful variational iteration method. The fractional derivatives are described in the Caputo sense. Several examples are given and the results are compared to exact solutions. The results show that the variational iteration method is very effective, convenient and simple to use.     

A compact finite difference method for the solution of the generalized Burgers–Fisher equation

In this article, numerical solutions of the generalized Burgers–Fisher equation are obtained using a compact finite difference method with minimal computational effort. To verify this, a combination of a sixth‐order compact finite difference scheme in space and a low‐storage third‐order total variation diminishing Runge–Kutta scheme in time have been used. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. The approximate solutions to the equation have been computed without transforming the equation and without using linearization. Comparisons indicate that there is a very good agreement between the numerical solutions and the exact solutions in terms of accuracy. The present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Numerical solution of the Rosenau–KdV–RLW equation by operator splitting techniques based on B‐spline collocation method

In the present study, the operator splitting techniques based on the quintic B‐spline collocation finite element method are presented for calculating the numerical solutions of the Rosenau–KdV–RLW equation. Two test problems having exact solutions have been considered. To demonstrate the efficiency and accuracy of the present methods, the error norms L₂ and L_∞ with the discrete mass Q and energy E conservative properties have been calculated. The results obtained by the method have been compared with the exact solution of each problem and other numerical results in the literature, and also found to be in good agreement with each other. A Fourier stability analysis of each presented method is also investigated.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

A pseudospectral method for nonlinear Duffing equation involving both integral and non‐integral forcing terms

The Legendre pseudospectral method is developed for the numerical solution of nonlinear Duffing equation involving both integral and non‐integral forcing terms. By using differentiation matrix, the problem is reduced to the solution of a system of algebraic equations. The method is general, easy to implement, and yields very accurate results. Numerical experiments are presented to demonstrate the accuracy and the efficiency of the proposed computational procedure. Copyright © 2014 John Wiley & Sons, Ltd.     

A two-level finite element method for the Navier–Stokes equations based on a new projection

We consider a fully discrete two-level approximation for the time-dependent Navier–Stokes equations in two dimension based on a time-dependent projection. By defining this new projection, the iteration between large and small eddy components can be reflected by its associated space splitting. Hence, we can get a weakly coupled system of large and small eddy components. This two-level method applies the finite element method in space and Crank–Nicolson scheme in time. Moreover,the analysis and some numerical examples are shown that the proposed two-level scheme can reach the same accuracy as the classical one-level Crank–Nicolson method with a very fine mesh size h by choosing a proper coarse mesh size H. However, the two-level method will involve much less work.     

Explicit analytical solutions of the generalized Burger and Burger–Fisher equations by homotopy perturbation method

In this article, we apply the homotopy perturbation method (HPM) to obtain approximate analytical solutions of the generalized Burger and Burger‐Fisher (B–F) equations. Several numerical examples are given to illustrate the efficiency of the HPM. Comparison of the result obtained by the present method with exact solution reveals that the accuracy and fast convergence of the new method. It is predicted that the HPM can be found wide application in engineering problems. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A numerical approach for solving a class of the nonlinear Lane–Emden type equations arising in astrophysics

In this paper, a collocation method based on the Bessel polynomials is presented for the approximate solution of a class of the nonlinear Lane–Emden type equations, which have many applications in mathematical physics. The exact solution can be obtained if the exact solution is polynomial. In other cases, such as an increasing number of nodes, a good approximation can be obtained with applicable errors. In addition, the method is presented with error and stability analysis. The numerical results show the effectiveness of the method for this type of equations. Comparing the methodology with some known techniques shows that the present approach is relatively easy and highly accurate. Copyright © 2011 John Wiley & Sons, Ltd.     

Exponential finite-difference method applied to Korteweg–de Vries equation for small times

The Korteweg–de Vries equation is numerically solved by using the exponential finite-difference technique. The accuracy of computed solutions is examined by comparison with other numerical and analytical solutions using two examples. The close results agreement between the current results and the exact solutions confirms that the proposed finite-difference procedure is an effective technique for the solution of the Korteweg–de Vries equation at the small times.     

Exact travelling wave solutions for the generalized nonlinear Schrödinger (GNLS) equation with a source by Extended tanh–coth, sine–cosine and Exp-Function methods

The capability of Extended tanh–coth, sine–cosine and Exp-Function methods as alternative approaches to obtain the analytic solution of different types of applied differential equations in engineering mathematics has been revealed. In this study, the generalized nonlinear Schrödinger (GNLS) equation is solved by three different methods. To obtain the single-soliton solutions for the equation, the Extended tanh–coth and sine–cosine methods are used. Furthermore, for this nonlinear evolution equation the Exp-Function method is applied to derive various travelling wave solution. Results show that while the first two procedures easily provide a concise solution, the Exp-Function method provides a powerful mathematical means for solving nonlinear evolution equations in mathematical physics.     

A finite difference method for fractional partial differential equation

An implicit unconditional stable difference scheme is presented for a kind of linear space–time fractional convection–diffusion equation. The equation is obtained from the classical integer order convection–diffusion equations with fractional order derivatives for both space and time. First-order consistency, unconditional stability, and first-order convergence of the method are proven using a novel shifted version of the classical Grünwald finite difference approximation for the fractional derivatives. A numerical example with known exact solution is also presented, and the behavior of the error is examined to verify the order of convergence.     

Local projection stabilization for advection–diffusion–reaction problems: One-level vs. two-level approach

Local projection stabilization (LPS) of finite element methods is a new technique for the numerical solution of transport-dominated problems. The main aim of this paper is a critical discussion and comparison of the one- and two-level approaches to LPS for the linear advection–diffusion–reaction problem. Moreover, the paper contains several other novel contributions to the theory of LPS. In particular, we derive an error estimate showing not only the usual error dependence on the mesh width but also on the polynomial degree of the finite element space. Based on this error estimate, we propose a definition of the stabilization parameter depending on the data of the solved problem. Unlike other papers on LPS methods, we observe that the consistency error may deteriorate the convergence order. Finally, we explain the relation between the LPS method and residual-based stabilization techniques for simplicial finite elements.     

Higher-order compact finite difference method for systems of reaction–diffusion equations

This paper is concerned with a compact finite difference method for solving systems of two-dimensional reaction–diffusion equations. This method has the accuracy of fourth-order in both space and time. The existence and uniqueness of the finite difference solution are investigated by the method of upper and lower solutions, without any monotone requirement on the nonlinear term. Three monotone iterative algorithms are provided for solving the resulting discrete system efficiently, and the sequences of iterations converge monotonically to a unique solution of the system. A theoretical comparison result for the various monotone sequences is given. The convergence of the finite difference solution to the continuous solution is proved, and Richardson extrapolation is used to achieve fourth-order accuracy in time. An application is given to an enzyme–substrate reaction–diffusion problem, and some numerical results are presented to demonstrate the high efficiency and advantages of this new approach.     

Homotopy perturbation method for the mixed Volterra–Fredholm integral equations

This article presents a numerical method for solving nonlinear mixed Volterra–Fredholm integral equations. The method combined with the noise terms phenomena may provide the exact solution by using two iterations only. Two numerical illustrations are given to show the pertinent features of the technique. The results reveal that the proposed method is very effective and simple.     

Exact and numerical solution of Black–Scholes matrix equation

In this paper a new method for solving Black–Scholes equation is proposed. The approach is based on the Mellin transform. A numerical procedure for the approximation of the solution is given.     

Composite generalized Laguerre–Legendre pseudospectral method for Fokker–Planck equation in an infinite channel

In this paper, we propose a composite generalized Laguerre–Legendre pseudospectral method for the Fokker–Planck equation in an infinite channel, which behaves like a parabolic equation in one direction, and behaves like a hyperbolic equation in other direction. We establish some approximation results on the composite generalized Laguerre–Legendre–Gauss–Radau interpolation, with which the convergence of proposed composite scheme follows. An efficient implementation is provided. Numerical results show the spectral accuracy in space of this approach and coincide well with theoretical analysis. The approximation results and techniques developed in this paper are also very appropriate for many other problems on multiple-dimensional unbounded domains, which are not of standard types.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.     

The convergence of Jacobi–Davidson iterations for Hermitian eigenproblems

Rayleigh quotient iteration is an iterative method with some attractive convergence properties for finding (interior) eigenvalues of large sparse Hermitian matrices. However, the method requires the accurate (and, hence, often expensive) solution of a linear system in every iteration step. Unfortunately, replacing the exact solution with a cheaper approximation may destroy the convergence. The (Jacobi‐) Davidson correction equation can be seen as a solution for this problem. In this paper we deduce quantitative results to support this viewpoint and we relate it to other methods. This should make some of the experimental observations in practice more quantitative in the Hermitian case. Asymptotic convergence bounds are given for fixed preconditioners and for the special case if the correction equation is solved with some fixed relative residual precision. A dynamic tolerance is proposed and some numerical illustration is presented. Copyright © 2002 John Wiley & Sons, Ltd.     

A level-2 reformulation–linearization technique bound for the quadratic assignment problem

This paper studies polyhedral methods for the quadratic assignment problem. Bounds on the objective value are obtained using mixed 0–1 linear representations that result from a reformulation–linearization technique (rlt). The rlt provides different “levels” of representations that give increasing strength. Prior studies have shown that even the weakest level-1 form yields very tight bounds, which in turn lead to improved solution methodologies. This paper focuses on implementing level-2. We compare level-2 with level-1 and other bounding mechanisms, in terms of both overall strength and ease of computation. In so doing, we extend earlier work on level-1 by implementing a Lagrangian relaxation that exploits block-diagonal structure present in the constraints. The bounds are embedded within an enumerative algorithm to devise an exact solution strategy. Our computer results are notable, exhibiting a dramatic reduction in nodes examined in the enumerative phase, and allowing for the exact solution of large instances.     

Application of Exp-function method to Dullin–Gottwald–Holm equation

In this paper, we adopt the Exp-function method and the traveling-wave transformation to study the so-called DGH equation, as a result a number of exact solutions of this equation have been found. The family of solution including some exact solutions such as solitary wave pattern, periodic traveling-wave solution, kink-wave solution and new bounded-wave solutions. And explained some of the solutions physical meaning.