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.     

A non-standard finite difference scheme for a diffusive HBV infection model with capsids and time delay

In this paper, a non-standard finite difference (NSFD) scheme for a delayed diffusive hepatitis B virus (HBV) infection model with intracellular HBV DNA-containing capsids is proposed. Dynamic consistency of this NSFD scheme is achieved by showing that the scheme preserves the non-negativity and boundedness of the solutions and the global stability of the homogeneous steady states of the corresponding continuous model without any restriction on spatial and temporal grid sizes. We prove the global stability of the steady states by constructing suitable discrete Lyapunov functions.     

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.     

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.     

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.     

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.     

Square-root dynamics of a giving up smoking model

The paper presents a new giving up smoking model for which interaction term is square root of potential and occasional smokers of model presented in Zaman (2011) [15]. First, we will show formulation of the model. Then we will discuss local and global stability of the model and its general solutions. The non-standard finite difference method (NSFD) is used to solve the new giving up smoking model. Both non-negativity and conservative law for differential equations system are discussed. Numerical results are presented graphically and compared well with those obtained by Runge–Kutta fourth-order method (RK4) and ODE45.     

Global Aspects of Age-Structured Cigarette Smoking Model

Smoking impacts health and as a result creates several problems related to age which means smoking has a strong correlation with age. Keeping this problem in view, we consider the global asymptotic properties of age-structured smoking model. First, we formulate the model and present the existence and uniqueness of solution. Then we discuss the equilibrium points and construct the Lyapunov function to examine global stability of the free smoking and positive smoking equilibrium points. Finally, we fixed the age factor and use the non-standard finite difference (NSFD) scheme for numerical solutions and compare our results obtained with RK4 and ODE45 graphically with the help of MATLAB.     

Discretization of an eco-epidemiological model and its dynamic consistency

In this paper, we discretize a continuous-time eco-epidemiological model by non-standard finite difference (NSFD) scheme as well as standard Euler forward scheme. Dynamical properties of both the systems are explored and compared with their continuous-time model. We show that the solution of NSFD system remains positive for all positive initial values. Fixed points and their local stability properties are shown to be identical with that of the continuous model, indicating its dynamic consistency. Dynamics of the Euler model, however, depend on the step–size and therefore dynamically inconsistent. Solutions in this latter method may be negative and allows numerical instabilities, leading to chaos. Extensive numerical simulations have been performed to validate the theoretical results.     

Numerical treatment for nonlinear Brusselator chemical model

Nowadays, numerical models have great importance in every field of science, especially for solving the nonlinear differential equations, partial differential equations, biochemical reactions, etc. The total time evolution of the reactant species which interacts with other species is simulated by the Runge-Kutta of order four (RK4) and by Non-Standard finite difference (NSFD) method. A NSFD model has been constructed for the biochemical reaction problem and numerical experiments are performed for different values of discretization parameter h. The results are compared with the well-known numerical scheme, i.e. RK4. The developed scheme NSFD gives better results than RK4.     

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.     

Optimized composite finite difference schemes for atmospheric flow modeling

In this paper, we use some finite difference methods in order to solve an atmospheric flow problem described by an advection–diffusion equation. This flow problem was solved by Clancy using forward‐time central space (FTCS) scheme and is challenging to simulate due to large errors in phase and amplitude which are generated especially over long propagation times. Clancy also derived stability limits for FTCS scheme. We use Von Neumann stability analysis and the approach of Hindmarsch et al. which is an improved technique over that of Clancy in order to obtain the region of stability of some methods such as FTCS, Lax–Wendroff (LW), Crank–Nicolson. We also construct a nonstandard finite difference (NSFD) scheme. Properties like stability and consistency are studied. To improve the results due to significant numerical dispersion or numerical dissipation, we derive a new composite scheme consisting of three applications of LW followed by one application of NSFD. The latter acts like a filter to remove the dispersive oscillations from LW. We further improve the composite scheme by computing the optimal temporal step size at a given spatial step size using two techniques namely; by minimizing the square of dispersion error and by minimizing the sum of squares of dispersion and dissipation errors.     

Global dynamics of a discretized SIRS epidemic model with time delay

We derive a discretized SIRS epidemic model with time delay by applying a nonstandard finite difference scheme. Sufficient conditions for the global dynamics of the solution are obtained by improvements in discretization and applying proofs for continuous epidemic models. These conditions for our discretized model are the same as for the original continuous model.     

An NSFD scheme for a class of SIR epidemic models with vaccination and treatment

An non-standard finite difference scheme is employed to discuss a class of SIR epidemic model with vaccination and treatment. The dynamical properties of the discretized model are then analysed. The results demonstrate that the discretized epidemic model is dynamically consistent with the continuous model since it maintains essential properties of the corresponding continuous model, such as positivity property and boundness of solutions, equilibrium points and their local stability properties.     

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.     

An iterative approach to solve a nonlinear moving boundary problem describing the solvent diffusion within glassy polymers

This paper consists in studying a mathematical model of solvent diffusion through the glassy polymers as a one-dimensional moving boundary problem with kinetic undercooling. We establish an iterative variable time-step method based on a nonstandard finite difference (NSFD) scheme to solve the considered moving boundary problem. The monotonicity and positivity of the numerical solution are proved. The numerical approach is investigated for three test problems composed of constant and inconstant diffusion coefficients for different values of parameters to demonstrate the validity and ability of the method.     

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     

离散FitzHugh-Nagumo方程的整体吸引子和维数

该文对ＦｉｔｚＨｕｇｈ Ｎａｇｕｍｏ方程初边值问题用有限差分格式离散空间变量,证明了离散模型整体吸引子的存在性,并给出了与犿无关的Ｈａｕｓｄｏｒｆｆ维数和Ｆｒａｃｔａｌ维数上界估计。     

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.