Nonstandard Finite Difference Research Articles

A reaction–diffusion model of major emerging infectious diseases in a spatially heterogeneous environment and case study

During the outbreak of major emerging infectious diseases, virus droplets can survive in the environment as aerosols for hours to days, and their impact on human infectious diseases is often overlooked. In addition, the speed of transmission of infectious diseases is often closely related to transportation. Therefore, studying the impact of environmental viruses and transportation on disease development is significant for effective infectious disease prevention and control. We proposed a degenerate reaction–diffusion infectious disease model ([Formula: see text]) considering environmental virus interference and established a well-posedness and threshold system for this model. We have obtained the system solution approaches the disease-free equilibrium ([Formula: see text]) when the basic reproduction number [Formula: see text]. The system has at least one positive steady state solution (PSS) when [Formula: see text]. This paper used a non-standard finite difference method discretization model while data fitting and parameter estimation were performed based on data provided by the Health Commission. Further sensitivity analysis was conducted on [Formula: see text]. At the same time, we also discussed the impact of various parameters in the early stages of the outbreak of major emerging infectious diseases on the development of the disease. Research found that even if the contact rate between people is controlled at a shallow level, the disease may persist. In the early stages of major emerging outbreak of infectious diseases, immediately reducing the use of transportation can effectively reduce the speed of disease spread.

MODELING THE EFFECTS OF PESTS AND PESTICIDE ON CROP YIELDS IN A MULTIPLE CROPPING SYSTEM

Pest infestation poses a significant threat to agricultural crop yields, and to control it, farmers spray chemical pesticides. The persistent use of these chemical agents not only leads to pesticide residues within crops but also exerts collateral damage on the beneficial pest population. In this research work, we formulate a nonlinear mathematical model to assess the impacts of pesticide on crop yields within a multiple cropping system. Model analysis illustrates that crop consumption rates destabilize, and the spraying rate of pesticide stabilizes the system. Furthermore, we determine conditions for the global stability of the coexisting equilibrium and conduct a global sensitivity analysis to identify model parameters that significantly influence pest population density. Our findings emphasize that, for effective pest population control and enhanced crop yields, farmers should choose either pesticides with a high pest abatement rate or those with a higher pesticide uptake rate. Considering the spraying rate of pesticide as time-dependent, we also suggest an optimal control strategy to minimize the pest population and associated costs. We provide analytical results backed by numerical simulations implemented through the non-standard finite difference scheme to support our findings.

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.

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.

Exploring local and global stability of COVID-19 through numerical schemes

Respiratory sensitivity and pneumonia are possible outcomes of the coronavirus (COVID-19). Surface characteristics like temperature and sunshine affect how long the virus survives. This research article analyzes COVID-19 mathematical model behavior based on symptomatic and non-symptomatic individuals. In the reproductive model, the best result indicates the intensity of the epidemic. Our model remained stable at a certain point under controlled conditions after we evaluated a specific element. This approach is in place of traditional approaches such as Euler’s and Runge–Kutta’s. An unusual numerical approach known as the non-standard finite difference (NSFD) scheme is used in this article. This numerical approach gives us positivity. A dependable numerical analysis allowed us to evaluate different approaches and verify our theoretical results. Unlike the widely used Euler and RK4 approaches, we investigated the benefits of implementing NSFD schemes. By numerically simulating COVID-19 in a variety of scenarios, we demonstrated how our theoretical concepts work. The simulation findings support the usefulness of both approaches.

Nonstandard Finite Difference Time Domain Methodology to Simulate Light Propagation in Nonlinear Materials

We extend the nonstandard (NS) finite difference time domain (FDTD) methodology, originally developed to solve Maxwell’s equations in linear materials, to nonlinear ones. We validate it by computing harmonics generation in a nonlinear dielectric and comparing with theory. The methodology also applies to the quantum electrodynamics that describes the interaction of charged particles with electromagnetic fields, and also to the Ginzburg-Landau model of superconductivity.