Nonstandard finite difference schemes for a class of generalized convection–diffusion–reaction equations

In this work, a class of nonstandard finite difference (NSFD) schemes are proposed to approximate the solutions of a class of generalized convection–diffusion–reaction equations. First, in the case of no diffusion, two exact finite difference schemes are presented using the method of characteristics. Based on these two exact schemes, a class of exact schemes are presented by introducing a parameter α. Second, since the forms of these exact schemes are so complicated that they are not convenient to use, a class of NSFD schemes are derived from the exact schemes using numerical approximations. It follows that, under certain conditions about denominator function of time‐step sizes, these NSFD schemes are elementary stable and the solutions are positive and bounded. Third, by means of the Mickens' technique of subequations, a new class of implicit NSFD schemes are constructed for the full convection–diffusion–reaction equations. It is shown that, under certain parameters set, these NSFD schemes are capable of preserving the non‐negativity and boundedness of the analytical solutions. Finally, some numerical simulations are provided to verify the validity of our analytical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1288–1309, 2015     

A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Nonstandard finite difference schemes for Michaelis–Menten type reaction‐diffusion equations

We compare and investigate the performance of the exact scheme of the Michaelis–Menten (M–M) ordinary differential equation with several new nonstandard finite difference (NSFD) schemes that we construct using Mickens' rules. Furthermore, the exact scheme of the M–M equation is used to design several dynamically consistent NSFD schemes for related reaction‐diffusion equations, advection‐reaction equations, and advection‐reaction‐diffusion equations. Numerical simulations that support the theory and demonstrate computationally the power of NSFD schemes are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Numerical solutions of Burgers–Huxley equation by exact finite difference and NSFD schemes

In this paper, numerical solution of the Burgers–Huxley (BH) equation is presented based on the nonstandard finite difference (NSFD) scheme. At first, two exact finite difference schemes for BH equation obtained. Moreover an NSFD scheme is presented for this equation. The positivity, boundedness and local truncation error of the scheme are discussed. Finally, the numerical results of the proposed method with those of some available methods compared.     

Non-standard finite difference method applied to an initial boundary value problem describing hepatitis B virus infection

In this paper, two non-standard finite difference (NSFD) schemes are proposed for a mathematical model of hepatitis B virus (HBV) infection with spatial dependence. The dynamic properties of the obtained discretized systems are completely analyzed. Relying on the theory of M-matrix, we prove that the proposed NSFD schemes is unconditionally positive. Furthermore, we establish that the NSFD method used preserves all constant steady states of the corresponding continuous initial boundary value problem (IBVP) model. We prove that the conditions for those equilibria to be asymptotically stable are consistent with the continuous IBVP model independently of the numerical grid size. The global asymptotical properties of the HBV-free equilibrium of the proposed NSFD schemes are derived via the construction of a suitable discrete Lyapunov function, and coincides with the continuous system. This confirms that the discretized models are dynamically consistent since they maintain essential properties of the corresponding continuous IBVP model. Finally, numerical simulations are performed from which it is demonstrated that the proposed NSFD method is advantageous over the standard finite difference (SFD) method.     

Analysis of and numerical schemes for a mouse population model in Hantavirus epidemics

This paper considers a non-linear system of ordinary differential equations, which arises in the study of hantavirus epidemics. The system has the property that the total population obeys the logistic equation. For this system, we use linearization and the dynamical properties of the logistic equation to analyze the dynamics of the subpopulation system. In view of the usual numerical instabilities associated with standard finite difference methods used for simulating such systems, we construct non-standard finite difference (NSFD) schemes, which preserve the dynamic properties of the system, and may therefore be used for its simulation.     

Analysis of numerical schemes of a mathematical model for sickle cell depolymerization

The purpose of this paper is to present and discuss numerical schemes for a mathematical model that describes carbon monoxide mediated sickle cell polymer melting. Two Runge-Kutta methods are analyzed and shown to be unstable by calculating the first failure value of step size and displaying the bifurcation diagram of RK4. Two nonstandard finite difference (NSFD) schemes are proposed and analyzed; one is shown to be stable subject to a predictable bound on step size, while the second one is unconditionally stable.     

Lyapunov direct method for investigating stability of nonstandard finite difference schemes for metapopulation models

In this paper nonstandard finite difference (NSFD) schemes of two metapopulation models are constructed. The stability properties of the discrete models are investigated by the use of the Lyapunov stability theorem. As a result of this we have proved that the NSFD schemes preserve essential properties of the metapopulation models (positivity, boundedness and monotone convergence of the solutions, equilibria and their stability properties). Especially, the basic reproduction number of the continuous models is also preserved. Numerical examples confirm the obtained theoretical results of the properties of the constructed difference schemes. The method of Lyapunov functions proves to be much simpler than the standard method for studying stability of the discrete metapopulation model in our very recent paper.     

Stability properties of a nonstandard finite difference scheme for a hantavirus epidemic model

This paper studies the stability properties of a nonstandard finite difference (NSFD) scheme used to simulate the dynamics of a mouse population model in hantavirus epidemics. It is shown that this difference scheme and the underlying system of differential equations have the same dynamics. The proof uses the fact that the total population obeys the logistic equation, as well as techniques from calculus, graphical analysis, and dynamical systems.     

SYMPLECTIC SCHEMES FOR QUASILINEAR WAVE EQUATIONS OF KLEIN-GORDON AND SINE-GORDON TYPE

1. IntroductionThere has been much discussion recently aboat deswi symplectic numrical schemesfor both finite and idste dimensional Hamilonian systems ([1] - [8]). A dass of syInPlecticsMes for linear wav equations has been sUggested [1], and numerical exPerimellt8 for tlieseschemes have been conducted [8, 9]. Results shOw that the syInPlectic method8 are inhereotlyfree from arttheial dissipatinn and all kinds of nonHtalltoulan pollutions, and are thu8 Of highquality and resolution.In thi…     

Parameterized families of finite difference schemes for the wave equation

This article is devoted to an analysis of simple families of finite difference schemes for the wave equation. These families are dependent on several free parameters, and methods for obtaining stability bounds as a function of these parameters are discussed in detail. Access to explicit stability bounds such as those derived here may, it is hoped, lead to optimization techniques for so‐called spectral‐like methods, which are difference schemes dependent on many free parameters (and for which maximizing the order of accuracy may not be the defining criterion). Though the focus is on schemes for the wave equation in one dimension, the analysis techniques are extended to two dimensions; implicit schemes such as ADI methods are examined in detail. Numerical results are presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 463–480, 2004.     

Nonstandard finite difference method revisited and application to the Ebola virus disease transmission dynamics

We provide effective and practical guidelines on the choice of the complex denominator function of the discrete derivative as well as on the choice of the nonlocal approximation of nonlinear terms in the construction of nonstandard finite difference (NSFD) schemes. Firstly, we construct nonstandard one-stage and two-stage theta methods for a general dynamical system defined by a system of autonomous ordinary differential equations. We provide a sharp condition, which captures the dynamics of the continuous model. We discuss at length how this condition is pivotal in the construction of the complex denominator function. We show that the nonstandard theta methods are elementary stable in the sense that they have exactly the same fixed-points as the continuous model and they preserve their stability, irrespective of the value of the step size. For more complex dynamical systems that are dissipative, we identify a class of nonstandard theta methods that replicate this property. We apply the first part by considering a dynamical system that models the Ebola Virus Disease (EVD). The formulation of the model involves both the fast/direct and slow/indirect transmission routes. Using the specific structure of the EVD model, we show that, apart from the guidelines in the first part, the nonlocal approximation of nonlinear terms is guided by the productive-destructive structure of the model, whereas the choice of the denominator function is based on the conservation laws and the sub-equations that are associated with the model. We construct a NSFD scheme that is dynamically consistent with respect to the properties of the continuous model such as: positivity and boundedness of solutions; local and/or global asymptotic stability of disease-free and endemic equilibrium points; dependence of the severity of the infection on self-protection measures. Throughout the paper, we provide numerical simulations that support the theory.     

A note on an NSFD scheme for a mathematical model of respiratory virus transmission

We construct a non-standard finite difference (NSFD) scheme for an SIRS mathematical model of respiratory virus transmission. This discretization is in full compliance with the NSFD methodology as formulated by Mickens. By use of an exact conservation law satisfied by the SIRS differential equations, we are able to determine the corresponding denominator function for the discrete first-order time derivatives. Our scheme is dynamically consistent with the SIRS differential equations, since the conservation laws are preserved. Furthermore, the scheme is shown to satisfy a positivity condition for its solutions for all values of the time step size.     

Exact finite-difference schemes for first order differential equations having three distinct fixed-points

We construct nonstandard finite-difference (NSFD) schemes that provide exact numerical methods for a first-order differential equation having three distinct fixed-points. An explicit, but also nonexact, NSFD scheme is also constructed. It has the feature of preserving the critical properties of the original differential equation such as the positivity of the solutions and the stability behavior of the three fixed-points.     

高阶非线性波动方程的有限差分方法

本文研究一类广泛的高阶非线性波动方程组初边值问题的有限差分格式，用离散泛函分析方法和先验估计的技巧得到了有限差分格式的收敛性。     

On Error Bounds for Approximation Schemes for Non-Convex Degenerate Elliptic Equations

In this paper we provide estimates of the rates of convergence of monotone approximation schemes for non-convex equations in one space-dimension. The equations under consideration are the degenerate elliptic Isaacs equations with x-depending coefficients, and the results applies in particular to certain finite difference methods and control schemes based on the dynamic programming principle. Recently, Krylov, Barles, and Jakobsen obtained similar estimates for convex Hamilton–Jacobi–Bellman equations in arbitrary space-dimensions. Our results are only valid in one space-dimension, but they are the first results of this type for non-convex second-order equations.     

On a new fractional-order Logistic model with feedback control

In this paper, we formulate and analyze a new fractional-order Logistic model with feedback control, which is different from a recognized mathematical model proposed in our very recent work. Asymptotic stability of the proposed model and its numerical solutions are studied rigorously. By using the Lyapunov direct method for fractional dynamical systems and a suitable Lyapunov function, we show that a unique positive equilibrium point of the new model is asymptotically stable. As an important consequence of this, we obtain a new mathematical model in which the feedback control variables only change the position of the unique positive equilibrium point of the original model but retain its asymptotic stability. Furthermore, we construct unconditionally positive nonstandard finite difference(NSFD) schemes for the proposed model using the Mickens' methodology. It is worth noting that the constructed NSFD schemes not only preserve the positivity but also provide reliable numerical solutions that correctly reflect the dynamics of the new fractional-order model. Finally, we report some numerical examples to support and illustrate the theoretical results. The results indicate that there is a good agreement between the theoretical results and numerical ones.     

Optimal l~∞ error estimates of finite difference methods for the coupled Gross-Pitaevskii equations in high dimensions

Due to the difficulty in obtaining the a priori estimate,it is very hard to establish the optimal point-wise error bound of a finite difference scheme for solving a nonlinear partial differential equation in high dimensions(2D or 3D).We here propose and analyze finite difference methods for solving the coupled GrossPitaevskii equations in two dimensions,which models the two-component Bose-Einstein condensates with an internal atomic Josephson junction.The methods which we considered include two conservative type schemes and two non-conservative type schemes.Discrete conservation laws and solvability of the schemes are analyzed.For the four proposed finite difference methods,we establish the optimal convergence rates for the error at the order of O(h~2+τ~2)in the l~∞-norm(i.e.,the point-wise error estimates)with the time stepτand the mesh size h.Besides the standard techniques of the energy method,the key techniques in the analysis is to use the cut-off function technique,transformation between the time and space direction and the method of order reduction.All the methods and results here are also valid and can be easily extended to the three-dimensional case.Finally,numerical results are reported to confirm our theoretical error estimates for the numerical methods.     

Properties of eigenvalues and spectral singularities for impulsive quadratic pencil of difference operators

In this paper, we investigate the spectral analysis of impulsive quadratic pencil of difference operators. We first present a boundary value problem consisting one interior impulsive point on the whole axis corresponding to the above mentioned operator. After introducing the solutions of impulsive quadratic pencil of difference equation, we obtain the asymptotic equation of the function related to the Wronskian of these solutions to be helpful for further works, then we determine resolvent operator and continuous spectrum. Finally, we provide sufficient conditions guarenteeing finiteness of eigenvalues and spectral singularities by means of uniqueness theorems of analytic functions. The main aim of this paper is demonstrating the impulsive quadratic pencil of difference operator is of finite number of eigenvalues and spectral singularities with finite multiplicities which is an uninvestigated problem proposed in the literature.     

A predator-prey model with Beddington-DeAngelis functional response: a non-standard finite-difference method

In this paper, we transform a continuous-time predator-prey model with Beddington–DeAngelis functional response into a discrete-time model by nonstandard finite difference scheme (NSFD). The NSFD model shows complete dynamic consistency with its continuous counterpart for any step size. However, the discrete model of same continuous system obtained by Euler forward method shows dynamic inconsistency for larger step size. Extensive numerical simulations have been done to compare the dynamics of NSFD system and Euler system. Our analysis reveals that dynamics of NSFD model is independent of the step-size, whereas the dynamics of the standard discrete model completely depends on the step-size and produces spurious dynamics like chaos.