Analytical and numerical solutions of the Fitzhugh–Nagumo equation and their multistability behavior

In this paper, we propose an analytical method and a modification of explicit exponential finite difference method (EEFDM) for analytical and numerical solutions of the Fitzhugh–Nagumo (FN) and Newell–Whitehead (NW) equations. The method is improved computationally by using the Padé approximation technique. Furthermore, multistability behavior of traveling wave solutions of the FN and NW equations are examined in presence of external forcing. It is observed that there exist coexisting periodic and quasiperiodic orbits for the FN equation, where as only quasiperiodic orbits is observed in case of NW equation.     

Convergence and Mean-Square Stability of Exponential Euler Method for Semi-Linear Stochastic Delay Integro-Differential Equations

In this paper, the numerical methods for semi-linear stochastic delay integro-differential equations are studied. The uniqueness, existence and stability of analytic solutions of semi-linear stochastic delay integro-differential equations are studied and some suitable conditions for the mean-square stability of the analytic solutions are also obtained. Then the numerical approximation of exponential Euler method for semi-linear stochastic delay integro-differential equations is constructed and the convergence and the stability of the numerical method are studied. It is proved that the exponential Euler method is convergent with strong order $\frac{1}{2}$ and can keep the mean-square exponential stability of the analytical solutions under some restrictions on the step size. In addition, numerical experiments are presented to confirm the theoretical results.     

A Novel Numerical Simulations for Fornberg-Whitham and Modified Fornberg-Whitham Equations with Nonhomogeneous Boundary Conditions

In this study, the numerical solutions of the Fornberg-Whitham (FW) equation modeling the qualitative behavior of wave refraction and the modified Fornberg-Whitham (mFW) equation describing the solitary wave and peakon waves with a discontinuous first derivative at the peak have been obtained. To obtain numerical results, the collocation finite element method has been combined with quintic B-spline bases. Although there are solutions to these equations by semi-analytical and analytical methods in the literature, there are very few studies using numerical methods. The stability analysis of the applied method is examined by the von-Neumann Fourier series method. We have considered four test problems with nonhomogeneous boundary conditions that have analytical solutions to show the performance of the method. The numerical results of the two problems are compared with some studies in the literature. Additionally, peakon wave solutions and some new numerical results of the mFW equation, which are not available in the literature, are given in the last two problems. No comparison has been made since there are no numerical results in the literature for the last two problems. The error norms $L_{2}$ and $L_{\infty }$ are calculated to demonstrate the presented numerical scheme''s accuracy and efficiency. The advantage of the scheme is that it produces accurate and reliable solutions even for modest values of space and time step lengths, rather than small values that cause excessive data storage in the computation process. In general, large step lengths in the space and time directions result in smaller matrices. This means less storage on the computer and results in faster outcomes. In addition, the present method gives more accurate results than some methods given in the literature.     

Polynomial and rational wave solutions of Kudryashov-Sinelshchikov equation and numerical simulations for its dynamic motions

Polynomial and rational wave solutions of Kudryashov-Sinelshchikov equation and numerical simulations for its dynamic motions are investigated. Conservation flows of the dynamic motion are obtained utilizing multiplier approach. Using the unified method, a collection of exact solitary and soliton solutions of Kudryashov-Sinelshchikov equation is presented. Collocation finite element method based on quintic B-spline functions is implemented to the equation to evidence the accuracy of the proposed method by test problems. Stability analysis of the numerical scheme is studied by employing von Neumann theory. The obtained analytical and numerical results are in good agreement.     

Numerical Solution of Partial Differential Equations in Arbitrary Shaped Domains Using Cartesian Cut-Stencil Finite Difference Method. Part I: Concepts and Fundamentals

A new finite difference (FD) method, referred to as "Cartesian cut-stencil FD", is introduced to obtain the numerical solution of partial differential equations on any arbitrary irregular shaped domain. The 2^nd-order accurate two-dimensional Cartesian cut-stencil FD method utilizes a 5-point stencil and relies on the construction of a unique mapping of each physical stencil, rather than a cell, in any arbitrary domain to a generic uniform computational stencil. The treatment of boundary conditions and quantification of the solution accuracy using the local truncation error are discussed. Numerical solutions of the steady convection-diffusion equation on sample complex domains have been obtained and the results have been compared to exact solutions for manufactured partial differential equations (PDEs) and other numerical solutions.     

二维不可压缩流的谱方法

本文以二维涡度方程为模型,介绍了谱方法和拟谱方法以及它们与差分方法和有限元法相结合的混合解法．这些方法可推广应用于其它一些类似的非线性问题．本文还给出了这些方法的某些数值例子和误差估计结果     

Alternating group explicit method for the diffusion equation

A new alternating group explicit method is presented for the finite difference solution of the diffusion equation. The new method uses stable asymmetric approximations to the partial differential equation which, when coupled in groups of two adjacent points on the grid, result in implicit equations which can be easily converted to explicit form and which offer many advantages. By judicious alternation of this strategy on the grid points of the domain an algorithm which possesses unconditional stability is obtained. This approach also results in more accurate solutions because of truncation error cancellations. The stability, consistency, convergence and truncation error of the new method are briefly discussed and the results of numerical experiments presented.     

Galerkin finite element methods for two-dimensional RLW and SRLW equations

In this paper, Galerkin finite methods for two-dimensional regularized long wave and symmetric regularized long wave equation are studied. The discretization in space is by Galerkin finite element method and in time is based on linearized backward Euler formula and extrapolated Crank–Nicolson scheme. Existence and uniqueness of the numerical solutions have been shown by Brouwer fixed point theorem. The error estimates of linearlized Crank–Nicolson method for RLW and SRLW equations are also presented. Numerical experiments, including the error norms and conservation variables, verify the efficiency and accuracy of the proposed numerical schemes.     

The transport dynamics in complex systems governing by anomalous diffusion modelled with Riesz fractional partial differential equations

In this paper, numerical solutions of fractional Fokker–Planck equations with Riesz space fractional derivatives have been developed. Here, the fractional Fokker–Planck equations have been considered in a finite domain. In order to deal with the Riesz fractional derivative operator, shifted Grünwald approximation and fractional centred difference approaches have been used. The explicit finite difference method and Crank–Nicolson implicit method have been applied to obtain the numerical solutions of fractional diffusion equation and fractional Fokker–Planck equations, respectively. Numerical results are presented to demonstrate the accuracy and effectiveness of the proposed numerical solution techniques. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite difference method combined with differential quadrature method for numerical computation of the modified equal width wave equation

The aim of this study is to improve the numerical solution of the modified equal width wave equation. For this purpose, finite difference method combined with differential quadrature method with Rubin and Graves linearizing technique has been used. Modified cubic B‐spline base functions are used as base function. By the combination of two numerical methods and effective linearizing technique high accurate numerical algorithm is obtained. Three main test problems are solved for various values of the coefficients. To observe the performance of the present method, the error norms of the single soliton problem which has analytical solution are calculated. Besides these error norms, three invariants are reported. Comparison of the results displays that our algorithm produces superior results than those given in the literature.     

Smoothed Particle ElectroMagnetics: A mesh-free solver for transients

In this paper an advanced mesh-free particle method for electromagnetic transient analysis, is presented. The aim is to obtain efficient simulations by avoiding the use of a mesh such as in the most popular grid-based numerical methods. The basic idea is to obtain numerical solutions for partial differential equations describing the electromagnetic problem by using a set of particles arbitrarily placed in the problem domain. The mesh-free smoothed particle hydrodynamics method has been adopted to obtain numerical solution of time domain Maxwell's curl equations. An explicit finite difference scheme has been employed for time integration. Details about the numerical treatment of electromagnetic vector fields components are discussed. Two case studies in one and in two dimensions are reported. In order to validate the new proposed methodology, named as Smoothed Particle ElectroMagnetics, a comparison with the standard finite difference time domain method results is performed. The intrinsic adaptive capability of the proposed method, has been exploited by introducing irregular particles distribution.     

仿样有限条U变换逼近法的确解及其收敛性

仿样有限条法(spline finite strip method)是分析等截面结构最流行的数值方法之一.在以往的研究中,与一些基准问题的解析结果相比较,论证了该方法数值结果的有效性和收敛性,但至今未对该方法的精确解和显式误差项进行过数学推导,解析地论证过其收敛性.该文在对平板的分析中,使用酉变换(简称U变换)逼近法,导出了仿样有限条法精确的数学解,这是首次在公开文献中给出的精确解.和常规的仿样有限条法相比较,总矩阵方程的集成及其数值解都不同,U变换法的总矩阵方程,减少为仅含有2个未知量的方程,然后导出仿样有限条法显式的精确解.精确解按Taylor级数展开,导出误差项和收敛率,并和其他数值方法直接比较.在这一点上可以发现,仿样有限条法收敛速度和非协调有限元相同时,包含的未知量少得多,收敛率比常规的有限差分法快得多.     

Radial basis collocation method and quasi‐Newton iteration for nonlinear elliptic problems

This work presents a radial basis collocation method combined with the quasi‐Newton iteration method for solving semilinear elliptic partial differential equations. The main result in this study is that there exists an exponential convergence rate in the radial basis collocation discretization and a superlinear convergence rate in the quasi‐Newton iteration of the nonlinear partial differential equations. In this work, the numerical error associated with the employed quadrature rule is considered. It is shown that the errors in Sobolev norms for linear elliptic partial differential equations using radial basis collocation method are bounded by the truncation error of the RBF. The combined errors due to radial basis approximation, quadrature rules, and quasi‐Newton and Newton iterations are also presented. This result can be extended to finite element or finite difference method combined with any iteration methods discussed in this work. The numerical example demonstrates a good agreement between numerical results and analytical predictions. The numerical results also show that although the convergence rate of order 1.62 of the quasi‐Newton iteration scheme is slightly slower than rate of order 2 in the Newton iteration scheme, the former is more stable and less sensitive to the initial guess. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Nonstandard perturbation approximation and travelling wave solutions of nonlinear reaction diffusion equations

This paper deals with the construction of a nonstandard numerical method to compute the travelling wave solutions of nonlinear reaction diffusion equations at high wave speeds. Related general properties are studied using the perturbation approximation. At high wave speed the perturbation parameter approaches to zero and the problem exhibits a multiscale character. That is, there are thin layers where the solution varies rapidly, while away from these layers the solution behaves regularly and varies slowly. Most of the conventional methods fail to capture this layer behavior. Thus, the quest for some new numerical techniques that may handle the travelling wave solutions at high wave speeds earns relevance. In this paper, one such parameter robust nonstandard numerical scheme is constructed, in the sense that its numerical solution converges in the maximum norm to the exact solution uniformly well for all finite wave speeds. To overcome the difficulty due to the nonlinearity, the problem is linearized using the quasilinearization process followed by nonstandard finite difference discretization. An extensive amount of analysis is carried out which uses a suitable decomposition of the error into smooth and singular component and a comparison principle combined with appropriate barrier functions. The error estimates are obtained, which ensures uniform convergence of the method. A set of numerical experiment is carried out in support of the predicted theory that validates computationally the theoretical results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Mobile and immobile gaseous transport: Embedded analytical solutions to finite volume methods

The motivation is driven by deposition processes based on chemical vapor problems. The underlying model problem is based on coupled transport–reaction equations with mobile and immobile areas. We deal with systems of ordinary and partial differential equations. Such equation systems are delicate to solve and we introduce a novel solver method, that takes into account ways to solve analytically parts of the transport and reaction equations. The main idea is to embed the analytical and semianalytical solutions, which can then be explicitly given to standard numerical schemes of higher order. The numerical scheme is based on flux‐based characteristic methods, which is a finite volume method. Such a method is an attractive alternative to the standard numerical schemes, which fully discretize the full equations. We instead reduce the computational time while embedding fast computable analytical parts. Here, we can accelerate the solver process, with a priori explicitly given solutions. We will focus on the derivation of the analytical solutions for general and special solutions of the characteristic methods that are embedded into a finite volume method. In the numerical examples, we illustrate the higher‐order method for different benchmark problems. Finally, the method is verified with realistic results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

Error analysis of approximate solutions of non-linear Burgers’ equation

This paper deals with application of the maximum principle for differential equations to the finite difference method for determining upper and lower approximate solutions of the non-linear Burgers’ equation and their error range. In term of mathematical architecture, the paper is based on the maximum principle for parabolic differential equations to establish monotonic residual relations of the Burgers’ equation; and in terms of numerical method, it applies the finite difference method to discretize the equation, followed by use of the proposed Residual Correction Method for obtaining its optimal solutions under constraint conditions for inequalities. Derived by using this approach, the upper and lower transient approximate solutions are not just useful in analyzing the range of the maximum possible error between them and the analytic solutions correctly, and the numerical validations also indicate good accuracy in mean values of the upper and lower approximate solutions.     

Multiscale asymptotic expansions and numerical algorithms for the wave equations of second order with rapidly oscillating coefficients

In this paper, we discuss the multiscale analysis and numerical algorithms for the wave equations of second order with rapidly oscillating coefficients. The formal multiscale asymptotic expansions of the solutions for these problems in four specific cases are presented. Higher order corrector methods are constructed and associated explicit convergence rates are obtained in some cases. A multiscale numerical method and a symplectic geometric scheme are introduced. Finally, some numerical results and unsolved problems are presented, and these numerical results support strongly the convergence theorem of this paper.     

Practical finite‐analytic method for solving differential equations by compact numerical schemes

Methodology for development of compact numerical schemes by the practical finite‐analytic method (PFAM) is presented for spatial and/or temporal solution of differential equations. The advantage and accuracy of this approach over the conventional numerical methods are demonstrated. In contrast to the tedious discretization schemes resulting from the original finite‐analytic solution methods, such as based on the separation of variables and Laplace transformation, the practical finite‐analytical method is proven to yield simple and convenient discretization schemes. This is accomplished by a special universal determinant construction procedure using the general multi‐variate power series solutions obtained directly from differential equations. This method allows for direct incorporation of the boundary conditions into the numerical discretization scheme in a consistent manner without requiring the use of artificial fixing methods and fictitious points, and yields effective numerical schemes which are operationally similar to the finite‐difference schemes. Consequently, the methods developed for numerical solution of the algebraic equations resulting from the finite‐difference schemes can be readily facilitated. Several applications are presented demonstrating the effect of the computational molecule, grid spacing, and boundary condition treatment on the numerical accuracy. The quality of the numerical solutions generated by the PFAM is shown to approach to the exact analytical solution at optimum grid spacing. It is concluded that the PFAM offers great potential for development of robust numerical schemes. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Compensated split-step balanced methods for nonlinear stiff SDEs with jump-diffusion and piecewise continuous arguments

This paper deals with numerical solutions of nonlinear stiff stochastic differential equations with jump-diffusion and piecewise continuous arguments. By combining compensated split-step methods and balanced methods, a class of compensated split-step balanced (CSSB) methods are suggested for solving the equations. Based on the one-sided Lipschitz condition and local Lipschitz condition, a strong convergence criterion of CSSB methods is derived. It is proved under some suitable conditions that the numerical solutions produced by CSSB methods can preserve the mean-square exponential stability of the corresponding analytical solutions. Several numerical examples are presented to illustrate the obtained theoretical results and the effectiveness of CSSB methods. Moreover, in order to show the computational advantage of CSSB methods, we also give a numerical comparison with the adapted split-step backward Euler methods with or without compensation and tamed explicit methods.

     

Numerical simulations of the improved Boussinesq equation

In this study, numerical simulations of the improved Boussinesq equation are obtained using two finite difference schemes and two finite element methods, based on the second‐and third‐order time discretization. The methods are tested on the problems of propagation of a soliton and interaction of two solitons. After the L_∞ error norm is used to measure differences between the exact and numerical solutions, the results obtained by the proposed methods are compared with recently published results. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010