Numerical investigation of fractional model of Phytoplankton–Toxic Phytoplankton–Zooplankton system with convergence analysis

In this paper, a fractional order model of the phytoplankton–toxic phytoplankton–zooplankton system with Caputo fractional derivative is investigated via three computational methods, namely, residual power series method (RPSM), homotopy perturbation Sumudu transform method (HPSTM) and the homotopy analysis Sumudu transform method (HASTM). This model is constituted by three components: phytoplankton, toxic phytoplankton and zooplankton. Phytoplankton species are self-feeding members of the plankton community and play a very significant role in ecosystems. A wide range of sea creatures get food through phytoplankton. This paper focuses on the implementation of the three above-mentioned computational methods for a nonlinear time-fractional phytoplankton–toxic phytoplankton–zooplankton (PTPZ) model with a perception to study the dynamics of a model. This study shows that the solutions obtained by employing the suggested computational methods are in good agreement with each other. The computational procedures reveal that the HASTM solution generates a more general solution as compared to RPSM and HPSTM and incorporates their results as a special case. The numerical results presented in the form of graphs authenticate the accuracy of computational schemes. Hence, the implemented methods are very appropriate and relevant to handle nonlinear fractional models. In addition, the effect of variation of fractional order of a time derivative and time [Formula: see text] on populations of phytoplankton, toxic–phytoplankton and zooplankton has also been studied through graphical presentations. Moreover, the uniqueness and convergence analyses of HASTM solution have also been discussed in view of the Banach fixed-point theory.

Free terminal time optimal control of a fractional-order model for the HIV/AIDS epidemic

In this paper, we present a fractional model for the HIV/AIDS epidemic and incorporate into the model control parameters of pre-exposure prophylaxis (PrEP), behavioral change and antiretroviral therapy (ART) aimed at controlling the spread of diseases. We prove the local and global asymptotic stability of disease-free and endemic equilibria of the model. We present a general fractional optimal control problem (FOCP) with free terminal time and develop the Adapted Forward-Backward Sweep method for numerical solving of the FOCP. Necessary conditions for a state/control/terminal time triplet to be optimal are obtained. The results show that the use of all controls increases the life expectancy of HIV-treated patients with ART and remarkably increases the number of people undergoing PrEP and changing their sexual habits. Also, when the derivative order [Formula: see text] ([Formula: see text]) limits to 1, the value of optimal terminal time increases while the value of objective functional decreases.

Impact of fear on delayed three species food-web model with Holling type-II functional response

This paper deals with the investigation of the three species food-web model. This model includes two logistically growing interaction species, namely [Formula: see text] and [Formula: see text], and the third species [Formula: see text] behaves as the predator and also host for [Formula: see text]. The species [Formula: see text] predating on the species [Formula: see text] with the Holling type-II functional response, while the first species [Formula: see text] is benefited from the third species [Formula: see text]. Further, the effect of fear is incorporated in the growth rate of species [Formula: see text] due to the predator [Formula: see text] and time lag in [Formula: see text] due to the gestation process. We explore all the biologically possible equilibrium points, and their local stability is analyzed based on the sample parameters. Next, we investigate the occurrence of Hopf-bifurcation around the interior equilibrium point by taking the value of the fear parameter as a bifurcation parameter for the non-delayed system. Moreover, we verify the local stability and existence of Hopf-bifurcation for the corresponding delayed system. Also, the direction and stability of the bifurcating periodic solutions are determined using the normal form theory and the center manifold theorem. Finally, we perform extensive numerical simulations to support the evidence of our analytical findings.

Traveling wave solutions in predator–prey models with competition

This paper studies the minimal wave speed of traveling wave solutions in predator–prey models, in which there are several groups of predators that compete among different groups. We investigate the existence and nonexistence of traveling wave solutions modeling the invasion of predators and coexistence of these species. When the positive solution of the corresponding kinetic system converges to the unique positive steady state, a threshold that is the minimal wave speed of traveling wave solutions is obtained. To finish the proof, we construct contracting rectangles and upper–lower solutions and apply the asymptotic spreading theory of scalar equations. Moreover, multiple propagation thresholds in the corresponding initial value problem are presented by numerical examples, and one threshold may be the minimal wave speed of traveling wave solutions.

Dynamics in a diffusive predator–prey system with double Allee effect and modified Leslie–Gower scheme

In this paper, we investigate a diffusive modified Leslie–Gower predator–prey system with double Allee effect on prey. The global existence, uniqueness and a priori bound of positive solutions are determined. The existence and local stability of constant steady–state solutions are analyzed. Next, we induce the nonexistence of nonconstant positive steady–state solutions, which indicates the effect of large diffusivity. Furthermore, we discuss the steady–state bifurcation and the existence of nonconstant positive steady–state solutions by the bifurcation theory. In addition, Hopf bifurcations of the spatially homogeneous and inhomogeneous periodic orbits are studied. Finally, we make some numerical simulations to validate and complement the theoretical analysis. Our results demonstrate that the dynamics of the system with double Allee effect and modified Leslie–Gower scheme are richer and more complex.

An efficient method for oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics

In this paper, a novel numerical technique, the first-order Smooth Composite Chebyshev Finite Difference method, is presented. Imposing a first-order smoothness of the approximation polynomial at the ends of each subinterval is originality of the method. Both round-off and truncation error analyses of the method are performed beside the convergence analysis. Diffusion of oxygen in a spherical cell including nonlinear uptake kinetics is solved by using the method. The obtained results are compared with the existing methods in the literature and it is observed that the proposed method gives more reliable results.

Dynamical Aspects of Corona Virus Infection

In this paper, we consider five mathematical models of corona virus infection. The first model is a mathematical model of corona virus entry. The second model is a mathematical model for interactions of virus N-protein and viral RNA. Here, we prove that phosphorylated N protein increases the affinity of viral RNA. The third model is a mathematical model of virion assembly. It is a six-dimensional model. But there is an invariant three-dimensional submodel, which we can prove has a positive stable equilibrium. The fourth model is a model of an enzyme inhibitor that blocks viral replication. The fifth model is a model of a virus and a vaccine.

Numerical study of Zika model as a mosquito-borne virus with non-singular fractional derivative

Zika virus infection, which is closely related to dengue, is becoming a global threat to human society. The transmission of the Zika virus from one human to another is spread by bites of Aedes mosquitoes. Recent studies also reveal the fact that the Zika virus can be transmitted by sexual interaction. In this paper, we use the fractional derivative with Mittag–Leffler non-singular kernel to study Zika virus transmission dynamics. This fractional is also known as the Atangana–Baleanu Caputo (ABC) derivative which is employed for the resulting system of ordinary differential equations. We investigate the proposed Zika virus model by using the Legendre spectral method. The model parameters are estimated and validated numerically by investigating the effect of fractional order exponent on various cases such as Susceptible human, infected human, asymptomatic carrier, exposed human, and recovered human. Numerical results obtained with the proposed method have been compared with exact solutions, showing in all parameters a very satisfactory agreement.

An efficient numerical simulation and mathematical modeling for the prevention of tuberculosis

The main purpose of this research is to use a fractional-mathematical model including Atangana–Baleanu derivatives to explore the clinical associations and dynamical behavior of the tuberculosis. Herein, we used a lately introduced fractional operator having Mittag-Leffler kernel. The existence and inimitability problems to the relevant model were examined through the fixed-point theory. To verify the significance of the arbitrary fractional-order derivative, numerical outcomes were explored from the biological and mathematical viewpoints using the values of model parameters. The graphical simulations show the comparison of the predictor–corrector method (PCM) and Caputo method (CM) for different fractional orders and the results indicated the significant preference of PCM over CM.

Dynamical behavior of recurrent neural networks with different external inputs

The main aim of this paper is to study the dynamics of a recurrent neural networks with different input currents in terms of asymptotic point. Under certain circumstances, we studied the existence, the uniqueness of bounded solutions and their homoclinic and heteroclinic motions of the considered system with rectangular currents input. Moreover, we studied the unpredictable behavior of the continuous high-order recurrent neural networks and the discrete high-order recurrent neural networks. Our method was primarily based on Banach’s fixed-point theorem, topology of uniform convergence on compact sets and Gronwall inequality. For the demonstration of theoretical results, we give examples and their numerical simulations.