Analysis of a mathematical model for cancer treatment by the nonstandard finite difference methods

We construct two nonstandard finite difference schemes and use them to study a mathematical model of cancer therapy. Several recent studies show various aspects of the immune response against the cancer. Our discrete models emphasize the role of antibodies in any form of therapy by taking into account the development of anticancer therapies (chemotherapy, immunotherapy, radiation therapy). The nonstandard finite difference models are implemented by using Matlab. Numerical simulations show the existence of a separation line between the basins of attraction of cancerous cell-free and the highest equilibrium cancerous cell.     

A nonstandard finite difference method for solving a mathematical model of HIV-TB co-infection

We design and analyze an unconditionally convergent nonstandard finite-difference method to study transmission dynamics of a mathematical model of HIV-TB co-infection. The dynamics of this model are studied using the qualitative theory of dynamical systems. These qualitative features of the continuous model are preserved by the numerical method that we propose in this paper. This method also preserves positivity of the solution which is one of the essential requirements when modelling epidemic diseases. Furthermore, we show that the numerical method is unconditionally stable. Competitive numerical results confirming theoretical investigations are provided. Comparisons are also made with other conventional approaches that are routinely used to solve these types of problems.     

A nonstandard finite difference scheme for contaminant transport with kinetic langmuir sorption

This study focuses on a contaminant transport model with Langmuir sorption under nonequilibrium conditions. The numerical instabilities of the standard finite difference schemes including the upwind method are investigated. By using the nonstandard finite difference method, a better finite difference model is constructed. The numerical simulation on a specific system configuration proves the advantages of the new finite difference model. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 767–785, 2011     

Exact finite difference scheme for linear differential equation with constant coefficients

An exact finite difference equation for the n-th order linear differential equation with real, constant coefficients is constructed. The exact finite difference scheme is expressed differently but equivalent to that given by Potts [3].     

Positive numerical solution for a nonarbitrage liquidity model using nonstandard finite difference schemes

In this article, we construct a numerical method based on a nonstandard finite difference scheme to solve numerically a nonarbitrage liquidity model with observable parameters for derivatives. This nonlinear model considers that the parameters involved are observable from order book data. The proposed numerical method use a exact difference scheme in the linear convection‐reaction term, and the spatial derivative is approximated using a nonstandard finite difference scheme. It is shown that the proposed numerical scheme preserves the positivity as well as stability and consistence. To illustrate the accuracy of the method, the numerical results are compared with those produced by other methods. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 210‐221, 2014     

Bifurcation analysis of a non-standard finite difference scheme for a time-delayed model of asset prices

In this paper, we apply a non-standard finite difference scheme to a time-delayed model of speculative asset markets and discuss the effect of time delay on the dynamics of asset prices. Firstly, the stability of the positive equilibrium of the system is investigated by analysing the characteristic equation. By choosing the time delay as a bifurcation parameter, we prove that Hopf bifurcations occur when the delay passes a sequence of critical values. Then, the explicit algorithm for determining the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions are derived. Finally, some numerical simulations are given to verify the theoretical analysis.     

A Newton linearized compact finite difference scheme for one class of Sobolev equations

In this article, a Newton linearized compact finite difference scheme is proposed to numerically solve a class of Sobolev equations. The unique solvability, convergence, and stability of the proposed scheme are proved. It is shown that the proposed method is of order 2 in temporal direction and order 4 in spatial direction. Moreover, compare to the classical extrapolated Crank‐Nicolson method or the second‐order multistep implicit–explicit methods, the proposed scheme is easier to be implemented as it only requires one starting value. Finally, numerical experiments on one and two‐dimensional problems are presented to illustrate our theoretical results.     

Asymptotical stability of partial difference equations with variable coefficients

The dynamical systems considered have scalar state, are multivariate, linear, time-discrete, and time-variable and are described by an initial value problem for a class of evolutionary partial difference equations. The time set is the nonnegative part of the integer lattice in several dimensions. Parts of the asymptotical stability set in the parameter space spanned by the time-variable coefficients are explicitly found. To assess the quality of the sufficient stability criteria, a comparison with the exact stability set is made in an example.     

A new positivity‐preserving nonstandard finite difference scheme for the DWE

An improved positivity‐preserving nonstandard finite difference scheme for the linear damped wave equation is presented. Unlike an earlier such scheme developed by the authors, the new scheme involves three time levels and is therefore able to include the effects of the equation's relaxation coefficient. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential, 2005     

Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition

An essential feature of nonstandard finite difference schemes for differential equations is the precise manner in which the discretization of derivatives is made. We demonstrate, for differential equations modeling systems where the solutions satisfy a positivity condition, that procedures can be formulated to calculate the so‐called denominator functions that appear in the discrete derivatives. These procedures are applied to a number of both ordinary and partial model differential equations to illustrate their use. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Oscillation and global asymptotic stability in a discrete epidemic model

In this paper, we study the oscillation, global asymptotic stability, and other properties of the positive solutions of the difference equation

     

A monotone scheme for Hamilton–Jacobi equations via the nonstandard finite difference method

A usual way of approximating Hamilton–Jacobi equations is to couple space finite element discretization with time finite difference discretization. This classical approach leads to a severe restriction on the time step size for the scheme to be monotone. In this paper, we couple the finite element method with the nonstandard finite difference method, which is based on Mickens' rule of nonlocal approximation. The scheme obtained in this way is unconditionally monotone. The convergence of the new method is discussed and numerical results that support the theory are provided. Copyright © 2009 John Wiley & Sons, Ltd.     

A positivity‐preserving nonstandard finite difference scheme for the damped wave equation

A positivity‐preserving nonstandard finite difference scheme is constructed to solve an initial‐boundary value problem involving heat transfer described by the Maxwell‐Cattaneo thermal conduction law, i.e., a modified form of the classical Fourier flux relation. The resulting heat transport equation is the damped wave equation, a PDE of hyperbolic type. In addition, exact analytical solutions are given, special cases are mentioned, and it is noted that the positivity condition is equivalent to the usual linear stability criteria. Finally, solution profiles are plotted and possible extensions to a delayed diffusion equation and nonlinear reaction‐diffusion systems are discussed. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 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.     

Nonstanard finite difference scheme for a scalar reaction-convection PDE

This paper is dedicated to my friend and colleague, professor Greey Ladas, in celebration of his sixtieth birthday. May his outstanding contributions to mathematics continue for many years.

We construct a nonstrandard finite difference shceme for a scalar reaction-convection partial differntial equation (PDE) having a symmetric, nonlinear cubic source term in the dependen variable. Is mathematical properies are studied and compared with those of the original PDE. We also construct the exact scheme for this equations.     

Global stability and a non-standard finite difference scheme for a diffusion driven HBV model with capsids

A diffusion driven model for hepatitis B virus (HBV) infection, taking into account the spatial mobility of both the HBV and the HBV DNA-containing capsids is presented. The global stability for the continuous model is discussed in terms of the basic reproduction number. The analysis is further carried out on a discretized version of the model. Since the standard finite difference (SFD) approximation could potentially lead to numerical instability, it has to be restricted or eliminated through dynamic consistency. The latter is accomplished by using a non-standard finite difference (NSFD) scheme and the global stability properties of the discretized model are studied. The results are numerically illustrated for the dynamics and stability of the various populations in addition to demonstrating the advantages of the usage of NSFD method over the SFD scheme.     

A finite difference scheme for solving the heat transport equation at the microscale

Heat transport at the microscale is of vital importance in microtechnology applications. In this study, we develop a finite difference scheme of the Crank‐Nicholson type by introducing an intermediate function for the heat transport equation at the microscale. It is shown by the discrete energy method that the scheme is unconditionally stable. Numerical results show that the solution is accurate. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 697–708, 1999