Nonstandard Finite Difference Research Articles

Nonstandard finite difference methods for a convective predator-prey pursuit and evasion model

One standard and two nonstandard finite difference methods are used to solve a convective predator-prey model consisting of cross-diffusion terms for which no exact solution is known. The model involves a system of strongly coupled nonlinear parabolic partial differential equations. Through some numerical experiments, we observe that reasonable numerical solutions are obtained when the temporal step size satisfies k ≤ 0.0005 with spatial step size h = 0.1 . We construct a nonstandard method using ideas from Anguelov et al. and the scheme is denoted by NSFD1. We obtain four conditions for NSFD1 to replicate the positivity property of the continuous model, but these conditions consist of many parameters. Therefore, an unconditionally positivity-preserving scheme is considered which is denoted as NSFD2. However, we found that NSFD2 is in general not consistent. We obtain conditions under which NSFD2 is both consistent and bounded and some numerical results are presented when these conditions are satisfied.

Development of a Compartment Model to Study the Pharmacokinetics of Medical THC after Oral Administration

The therapeutic potential of delta9-tetrahydrocannabinol (THC), a primary cannabinoid in the cannabis plant, has led to its development into oral medical products for treating various conditions. However, THC, being a psychoactive substance, can lead to addiction if taken in inappropriate amounts. Thus, studying the pharmacokinetics of THC is crucial for understanding how the drug behaves in the body after administration. This study aims to develop a multi-compartmental model to investigate the pharmacokinetics of medical THC and its metabolites after oral administration. Using the law of mass action, the model was converted into ordinary differential equations (ODEs) to describe the rate of concentration changes of THC and its metabolites in each compartment. The nonstandard finite difference (NSFD) method was then applied to construct numerical solution schemes, which were implemented in MATLAB along with estimated pharmacokinetic rate constants. The results demonstrate that the simulation curves depicting the plasma concentration–time profiles of THC and 11-hydroxy-THC (THC-OH) closely resemble actual data samples, indicating the model’s accuracy. Moreover, the model predicts the pharmacokinetics of THC and its metabolites in various tissues. Consequently, this model serves as a valuable tool for enhancing our understanding of the pharmacokinetics of THC and its metabolites, guiding dosage adjustments, and determining administration durations for oral medical THC.

Analysis of a diffusive chemical reaction model in three space dimensions

This article proposes an implicit operator splitting nonstandard finite difference (OS-NSFD) scheme for numerical treatment of two species in three space dimensions reaction-diffusion glycolysis model. Since, the unknown state variables exhibiting the concentrations of species in glycolysis models and they cannot be negative and obtaining their positive solutions is a challenging task. The established theoretical result ensures that our proposed OS-NSFD scheme is unconditionally convergent at equilibrium point and fulfills the condition of positivity of solutions on contrary to other methods. Further, we analyze the existence and uniqueness of the solution obtained for the underlying system. To highlight the effectiveness of OS-NSFD scheme we compare the simulation results of OS-NSFD scheme with three well-known existing operators splitting finite difference (FD) schemes, namely, forward Euler explicit, backward Euler implicit and Crank Nicolson splitting schemes. Many existing techniques provide with the restricted positive solutions which do not work always. These techniques are only applicable if certain conditions on the discretized parameters are considered otherwise; they produce negative solutions, which is not the physical feature of the real system. The current work bridges this gap by catering the unconditional positive solutions to the reaction diffusion models.

An efficient numerical approach for singularly perturbed time delayed parabolic problems with two-parameters

ObjectivesThe main objective of this work is to design an efficient numerical scheme is proposed for solving singularly perturbed time delayed parabolic problems with two parameters.ResultsThe scheme is constructed via the implicit Euler and non-standard finite difference method to approximate the time and space derivatives, respectively. Besides, to enhance the accuracy and order of convergence of the method Richardson extrapolation technique is employed. Intensive numerical experimentation has been done on some model examples. Further, the layer behavior of the solutions is presented using graphs and observed to agree with the existing theories. Finally, error analysis of the scheme is done and observed that the proposed method is parameter uniform convergent with the order of convergence (Δt)2+h2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\left( (\\Delta t )^2+h^2\\right) $$\\end{document}

A fitted parameter convergent finite difference scheme for two-parameter singularly perturbed parabolic differential equation

The objective of this paper is to develop a numerical scheme that is uniform in its parameters for a specific type of time-dependent parabolic problem with two perturbation parameters. The existence of these two parameters in the terms with the highest-order derivatives results in the formation of boundary layer(s) in the solution of such problems. Solving these model problems using classical methods does not yield satisfactory results due to the layer behavior. Therefore, nonstandard finite difference schemes have been developed as a means to obtain numerical solutions for these problems. To develop the scheme, we employ the Crank-Nicolsondiscretization on a uniform time mesh and apply a fitted operator method with a uniform spatial mesh. We have established the stability and convergence of the proposed scheme. The proposed scheme exhibits uniform convergence of second order in the temporal direction and first order in the spatial direction. However, temporal mesh refinements is employed to enhance the order to two in both directions. Model examples are provided to validate the practicality of the proposed numerical scheme.

Some novel exact and non-standard finite difference schemes for a class of diffusion–advection–reaction equation

In this paper, a generalized diffusion–advection–reaction equation (GDARE) is considered and by using the solitary wave solution two new class of exact finite difference (EFD) schemes are proposed under certain functional relationship between the temporal and space step sizes. Suitable denominator functions are identified from these EFD schemes and a new non-standard finite difference (NSFD) scheme is developed. This NSFD scheme preserves positivity and boundedness properties of the solution of GDARE. The L 2 , L ∞ errors and CPU times are obtained for the proposed NSFD and EFD schemes to show the proficiency and accuracy of the schemes.

Two-dimensional distributed order Cable equation with non-singular kernel: A nonstandard implicit compact finite difference approach

Two-dimensional distributed order Cable equation with non-singular kernel: A nonstandard implicit compact finite difference approach

A comprehensive analysis of COVID-19 nonlinear mathematical model by incorporating the environment and social distancing

This research conducts a detailed analysis of a nonlinear mathematical model representing COVID-19, incorporating both environmental factors and social distancing measures. It thoroughly analyzes the model’s equilibrium points, computes the basic reproductive rate, and evaluates the stability of the model at disease-free and endemic equilibrium states, both locally and globally. Additionally, sensitivity analysis is carried out. The study develops a sophisticated stability theory, primarily focusing on the characteristics of the Volterra–Lyapunov (V–L) matrices method. To understand the dynamic behavior of COVID-19, numerical simulations are essential. For this purpose, the study employs a robust numerical technique known as the non-standard finite difference (NSFD) method, introduced by Mickens. Various results are visually presented through graphical representations across different parameter values to illustrate the impact of environmental factors and social distancing measures.

A multi-population approach to epidemiological modeling of Listeriosis transmission dynamics incorporating food and environmental contamination

Listeriosis is a food-borne disease that mainly affects pregnant women and newborns. We propose and analyze a deterministic model of Listeriosis by considering three groups of individuals: newborns, pregnant women, and others. Mathematical analysis of the model is performed, and equilibrium points are determined. The model has three equilibria, namely, the disease-free equilibrium, the bacteria-free equilibrium, and the endemic equilibrium. We use Castillo-Chavez theorem to establish the global stability of the disease-free equilibrium when the basic reproduction number is less than 1. The local asymptotic stability of the bacteria-free, and endemic equilibria are also established using the sign of the eigenvalues of the Jacobian matrix. We use the non-standard finite difference scheme and carried numerical simulations to confirm the theoretical results. We further show the impact of specific parameters on the dynamics of infectious individuals and observe that intervention is required in all the sub-populations by reducing the contact rate and vertical transmission to mininmize the number of infectious.

Variable order fractional diabetes models: numerical treatment

ABSTRACT In this paper, we present two numerical approximation solutions for the updated nonlinear variable order fractional differential equations of diabetes disease. The variable order fractional differential operator is based on the Caputo fractional derivative concept. The first used method is developed using the nonstandard finite difference process, such that a new quick approximation is applied to tackle the nonlinear terms. By applying some reasonable assumptions to the given data, we can prove that this scheme maintains the positivity and boundedness of the discretized solutions. The second method is based on approximation with shifted Jacobi polynomials. The properties of Jacobi polynomials, together with the shifted Jacobi-Gauss-Lobatto nodes were utilized to reduce the proposed system into a set of algebraic nonlinear equations. Furthermore, some numerical experiments are conducted to compare the performance of the introduced schemes. Our numerical simulations show that the nonstandard finite difference approach requires less computational time compared to the Jacobi-Gauss-Lobatto spectral collocation. While the last one is more accurate when used for solving the variable order fractional diabetes model.

Deterministic and Stochastic Nonlinear Model for Transmission Dynamics of COVID-19 with Vaccinations Following Bayesian-Type Procedure

We develop a mathematical model for the SARAS-CoV-2 double variant transmission characteristics with variant 1 vaccination to address this novel aspect of the disease. The model is theoretically examined, and adequate requirements are derived for the stability of its equilibrium points. The model includes the single variant 1 and variant 2 endemic equilibria in addition to the endemic and disease-free equilibria. Various approaches are used for the global and local stability of the model. For both strains, we determine the basic reproductive numbers R1 and R2. To investigate the occurrence of the layers (waves), we expand the model to include some analysis based on the second-order derivative. The model is then expanded to its stochastic form, and numerical outcomes are computed. For numerical purposes, we use the nonstandard finite difference method. Some error analysis is also recorded.

Analysis of a fractional-order model for dengue transmission dynamics with quarantine and vaccination measures

A comprehensive mathematical model is proposed to study two strains of dengue virus with saturated incidence rates and quarantine measures. Imperfect dengue vaccination is also assumed in the model. Existence, uniqueness and stability of the proposed model are proved using the results from fixed point and degree theory. Additionally, well constructed Lyapunov function candidates are also applied to prove the global stability of infection-free equilibria. It is also demonstrated that the model is generalized Ulam-Hyers stable under some appropriate conditions. The model is fitted to the real data of dengue epidemic taken from the city of Espirito Santo in Brazil. For the approximate solution of the model, a non-standard finite difference(NSFD) approach is applied. Sensitivity analysis is also carried out to show the influence of different parameters involved in the model. The behaviour of the NSFD is also assessed under different denominator functions and it is observed that the choice of the denominator function could influence the solution trajectories. Different scenario analysis are also assessed when the reproduction number is below or above one. Furthermore, simulations are also presented to assess the epidemiological impact of dengue vaccination and quarantine measures for infected individuals.

Analytical Properties and Numerical Preservation of an Age-Group Susceptible-Infected-Recovered Model: Application to the Diffusion of Information

Abstract This paper analyzes an age-group susceptible-infected-recovered (SIR) model. Theoretical results concerning the conservation of the total population, the positivity of the analytical solution, and the final size of the epidemic are derived. Since the model is a nonlinear system of ordinary differential equations (ODEs), a numerical approximation is considered, based on Standard and non-Standard Finite Difference methods, and on a Modified Patankar-Runge–Kutta (MPRK) method. The numerical preservation of the qualitative properties of the analytical solution is studied. The obtained results are applied to the diffusion of information in social networks, and the effectiveness of the different numerical approaches is shown through several numerical tests on real data.

Dynamical behavior of a hepatitis B epidemic model and its NSFD scheme

Hepatitis is inflammation of the liver, and one of its types, hepatitis B, is a contagious infection. Using mathematical models, the nature of the spread of the Hepatitis B virus can be predicted. In the present paper, a hepatitis B epidemic model with a Beddington–DeAngelis type incidence rate and a constant vaccination rate is considered. Some dynamical properties of this model, such as non-negativity, boundedness character, the basic reproduction number R0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {R}_0$$\\end{document}, stability nature, and the bifurcation phenomenon, are investigated. By the Bendixson theorem, it is demonstrated that the disease-free equilibrium is globally asymptotically stable. It is shown that a transcritical bifurcation phenomenon occurs when R0=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {R}_0=1$$\\end{document}. It is concluded that the endemic equilibrium is globally asymptotically stable when R0>1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {R}_0>1$$\\end{document}, by utilizing Dulac’s criteria. Also, a discrete system of difference equations is obtained by constructing a non-standard finite difference (NSFD) scheme for the continuous model. It is shown that the solutions of this discrete system are dynamically consistent for all finite step sizes. The theoretical results obtained are also supported and visualized by numerical simulations. These simulations also demonstrate that the NSFD scheme produces much more efficient results than the Euler or RK4 schemes, as shown in the theoretical results obtained.

Numerical method for second order singularly perturbed delay differential equations with fractional order in time via fitted computational method

The numerical solution of time fractional parabolic differential equations with singular perturbations and delay is the subject of this article. An arbitrarily small perturbation parameter ɛ(0<ɛ<<1) is multiplied the highest order derivative term. The problem’s solution shows an exponential boundary layer on the right side of the spatial domain for small values of ɛ. The solution and its derivatives are discussed together with their properties and bounds. The Crank–Nicolson approach for time discretization and the fitted operator non-standard finite difference method for space discretization are used to solve the time fractional singularly perturbed time delay problem under consideration. The comparison principle and solution bound are used to examine and analyze the scheme’s stability. The scheme’s uniform convergence is examined and demonstrated. The proposed scheme has a uniform convergence order of Mx−1+(Δt)2−γ. Two numerical illustrations for various values of ɛ are taken into consideration to validate the theoretical analysis of the scheme. The constructed scheme provides a better order of convergence and more accurate results than methods found in the literature.

Comparative analysis of a fractional co-infection model using nonstandard finite difference and two-step Lagrange polynomial methods

In this work, Caputo fractional model for co-dynamics of two-strains (Alpha and Delta variants) of SARS-CoV-2 and tuberculosis (TB) is formulated. We investigates the basic mathematical analysis, existence of unique solution, stability analysis, sensitivity analysis and simulations of the proposed model. In the basic mathematical analysis, it is proved that the model’s solutions are non-negative, bounded and locally asymptotically stable at disease free equilibrium (DFE) point. The existence of unique solution is proved by using Banach fixed point theorem (BFPT). Stability analysis is carried out by employing Hyers–Ulam stability criteria. Sensitivity analysis is performed to assess the impact of the model parameters on the general dynamics of the model by implementing partial rank correlation coefficient (PRCC) technique. Simulations of the proposed model obtained from non-standard finite difference scheme (NSFDs) are compared with the simulations obtained from two-step Lagrange polynomial method (TLPM). It is concluded that NSFDs gives more accurate and realistic results than TLPM. To assess some control measures necessary for reducing the co-spread of both diseases, comparative simulations are carried out for the infected classes. Also, numerical assessments showed that preventive efforts against SARS-CoV-2 variants could results in the reduction of TB prevalence and the co-infection of both diseases.

Global Properties of a Discrete SARS-CoV-2/HIV Co-Dynamics Model

Coronavirus disease 2019 (COVID-19), which is caused by the virus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a respiratory disease. In this paper, we analyze the global stability of a discrete SARS-CoV-2/HIV co-dynamics model. We create the discrete model by applying a nonstandard finite difference (NSFD) method. We demonstrate that NSFD retains essential solution properties, including positivity and boundedness. We determine the fixed points and identify their existence conditions. We investigate the global stability of these fixed points through the application of the Lyapunov method. To complement our analytical findings, we present numerical simulations.

Analyzing the dynamic patterns of COVID-19 through nonstandard finite difference scheme

This paper presents a novel approach to analyzing the dynamics of COVID-19 using nonstandard finite difference (NSFD) schemes. Our model incorporates both asymptomatic and symptomatic infected individuals, allowing for a more comprehensive understanding of the epidemic's spread. We introduce an unconditionally stable NSFD system that eliminates the need for traditional Runge–Kutta methods, ensuring dynamical consistency and numerical accuracy. Through rigorous numerical analysis, we evaluate the performance of different NSFD strategies and validate our analytical findings. Our work demonstrates the benefits of using NSFD schemes for modeling infectious diseases, offering advantages in terms of stability and efficiency. We further illustrate the dynamic behavior of COVID-19 under various conditions using numerical simulations. The results from these simulations demonstrate the effectiveness of the proposed approach in capturing the epidemic's complex dynamics.

Ebola virus disease model with a nonlinear incidence rate and density-dependent treatment

This paper studies an Ebola epidemic model with an exponential nonlinear incidence function that considers the efficacy and the behaviour change. The current model also incorporates a new density-dependent treatment that catches the impact of the disease transmission on the treatment. Firstly, we provide a theoretical study of the nonlinear differential equations model obtained. More precisely, we derive the effective reproduction number and, under suitable conditions, prove the stability of equilibria. Afterwards, we show that the model exhibits the phenomenon of backward-bifurcation whenever the bifurcation parameter and the reproduction number are less than one. We find that the bi-stability and backward-bifurcation are not automatically connected in epidemic models. In fact, when a backward-bifurcation occurs, the disease-free equilibrium may be globally stable. Numerically, we use well-known standard tools to fit the model to the data reported for the 2018–2020 Kivu Ebola outbreak, and perform the sensitivity analysis. To control Ebola epidemics, our findings recommend a combination of a rapid behaviour change and the implementation of a proper treatment strategy with a high level of efficacy. Secondly, we propose and analyze a fractional-order Ebola epidemic model, which is an extension of the first model studied. We use the Caputo operator and construct the Grünwald-Letnikov nonstandard finite difference scheme, and show its advantages.

Fractional epidemic model of coronavirus disease with vaccination and crowding effects

Most of the countries in the world are affected by the coronavirus epidemic that put people in danger, with many infected cases and deaths. The crowding factor plays a significant role in the transmission of coronavirus disease. On the other hand, the vaccines of the covid-19 played a decisive role in the control of coronavirus infection. In this paper, a fractional order epidemic model (SIVR) of coronavirus disease is proposed by considering the effects of crowding and vaccination because the transmission of this infection is highly influenced by these two factors. The nonlinear incidence rate with the inclusion of these effects is a better approach to understand and analyse the dynamics of the model. The positivity and boundedness of the fractional order model is ensured by applying some standard results of Mittag Leffler function and Laplace transformation. The equilibrium points are described analytically. The existence and uniqueness of the non-integer order model is also confirmed by using results of the fixed-point theory. Stability analysis is carried out for the system at both the steady states by using Jacobian matrix theory, Routh–Hurwitz criterion and Volterra-type Lyapunov functions. Basic reproductive number is calculated by using next generation matrix. It is verified that disease-free equilibrium is locally asymptotically stable if R0<1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${R}_{0}<1$$\\end{document} and endemic equilibrium is locally asymptotically stable if R0>1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${R}_{0}>1$$\\end{document}. Moreover, the disease-free equilibrium is globally asymptotically stable if R0<1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${R}_{0}<1$$\\end{document} and endemic equilibrium is globally asymptotically stable if R0>1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${R}_{0}>1$$\\end{document}. The non-standard finite difference (NSFD) scheme is developed to approximate the solutions of the system. The simulated graphs are presented to show the key features of the NSFD approach. It is proved that non-standard finite difference approach preserves the positivity and boundedness properties of model. The simulated graphs show that the implementation of control strategies reduced the infected population and increase the recovered population. The impact of fractional order parameter α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} is described by the graphical templates. The future trends of the virus transmission are predicted under some control measures. The current work will be a value addition in the literature. The article is closed by some useful concluding remarks.