CMES CMES CMES Computer Modeling in Engineering & Sciences 1526-1506 1526-1492 Tech Science Press USA 23059 10.32604/cmes.2023.023059 Article The Influence of Saturated and Bilinear Incidence Functions on the Dynamical Behavior of HIV Model Using Galerkin Scheme Having a Polynomial of Order Two The Influence of Saturated and Bilinear Incidence Functions on the Dynamical Behavior of HIV Model Using Galerkin Scheme Having a Polynomial of Order Two The Influence of Saturated and Bilinear Incidence Functions on the Dynamical Behavior of HIV Model Using Galerkin Scheme Having a Polynomial of Order Two Attaullah1attaullah@bkuc.edu.pk ZebKamil1 MohamedAbdullah2 Department of Mathematics & Statistics, Bacha Khan University, Charsadda, 24461, Pakistan Research Centre, Future University in Egypt, New Cairo, 11835, Egypt Corresponding Author: Attaullah. Email: attaullah@bkuc.edu.pk 2023 4 2 2023 136 2 1661 1685 0742022 2692022 © 2023 Attaullah et al. 2023 Attaullah et al. This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Mathematical modelling has been extensively used to measure intervention strategies for the control of contagious conditions. Alignment between different models is pivotal for furnishing strong substantiation for policymakers because the differences in model features can impact their prognostications. Mathematical modelling has been widely used in order to better understand the transmission, treatment, and prevention of infectious diseases. Herein, we study the dynamics of a human immunodeficiency virus (HIV) infection model with four variables: S (t), I (t), C (t), and A (t) the susceptible individuals; HIV infected individuals (with no clinical symptoms of AIDS); HIV infected individuals (under ART with a viral load remaining low), and HIV infected individuals (with two different incidence functions (bilinear and saturated incidence functions). A novel numerical scheme called the continuous Galerkin-Petrov method is implemented for the solution of the model. The influence of different clinical parameters on the dynamical behavior of S (t), I (t), C (t) and A (t) is described and analyzed. All the results are depicted graphically. On the other hand, we explore the time-dependent movement of nanofluid in porous media on an extending sheet under the influence of thermal radiation, heat flux, hall impact, variable heat source, and nanomaterial. The flow is considered to be 2D, boundary layer, viscous, incompressible, laminar, and unsteady. Sufficient transformations turn governing connected PDEs into ODEs, which are solved using the proposed scheme. To justify the envisaged problem, a comparison of the current work with previous literature is presented.

HIV/AIDS Galerkin technique bilinear and saturated incidence functions
Introduction

The study of epidemic models is a powerful tool for the dynamics of different infectious diseases in real-world phenomena. For the transmission dynamics of infectious diseases in a population, mathematicians and biologists used various epidemic models . There are innovative scientific advances and significant health intervention measures in the globe, yet HIV/AIDS remains one of humanity’s graves devastating diseases. Many countries are still seriously afflicted by this disease. Currently, the global spread of HIV infection is influencing the occurrence of other infectious diseases such as tuberculosis (TB) . HIV is a virus that causes HIV infection and is transferred during sexual activity, breastfeeding, and sharing injectable drug gear such as needles with HIV positive people. AIDS, the most severe stage of HIV infection, is triggered by the HIV pathogens. In 2018, the number of individuals living with HIV/AIDS and the number of deaths worldwide is expected to hit 37.9 million and 1.2 million, respectively. Approximately 62% of those infected were confirmed and started on antiretroviral therapy (ART) . Many therapies have been proposed to improve the quality of life of HIV patients, including antiretroviral therapy , chemotherapy, and stem cell therapy. Antiretroviral therapy, which is the most commonly used combination of drugs to treat HIV infection, has many side effects . Stem cell therapy is very limited due to the high cost of the procedure as well as the difficulty of obtaining healthy and consistent donors. A mathematical model is a mathematically based description of a dynamical system. It is essential for evaluating and controlling the HIV/AIDS infectious disease. Several assumptions and factors have substantial effects on the construction of a model, which may be changed employing controlling functions. Thus, using the idea of optimal control theory, a mathematical model of the HIV/AIDS pandemic can be reconstructed, and the disease’s regulating mechanisms may be studied. This theory contains several useful concepts that explain how disease, whether epidemic or pandemic, may be managed via biological controls. This concept has been adopted by many authors in order to control infection. Several HIV models have been developed in recent years to better understand the dynamics of HIV infection, disease progression, and the interaction of the immune system with HIV in the area of HIV infection of CD4+T cells.

Naresh et al.  presented a nonlinear HIV/AIDS mathematical model. They claimed that HIV infection has been reduced significantly because of increased awareness of HIV infectives as identified by screening and contact tracing, but that the illness remains prevalent due to immigration and the lack of contact tracking. Finally, they believe that the most effective way to minimize the disease burden is to spread awareness about HIV/AIDS. Nyabadza et al.  investigated a deterministic HIV/AIDS model that describes condom usage, HIV counselling and testing (HCT), and therapy. They examined the concept because HCT practice is still in its early stages. According to the model, this campaign has very little impact on reducing HIV endemicity. A mathematical model for HIV/AIDS dynamics was proposed by Mushanyu . He looked at the effects of HIV late diagnosis on the disease’s spread. His numerical findings show that early HIV/AIDS treatment motivation and improved HIV self-testing schedules offer more undiagnosed people the knowledge they need to know their HIV status, reducing HIV transmission. Ullah et al.  established an optimal control model for the COVID-19 pandemic. They used real-world data to quantitatively evaluate the model. They proved that the proposed method could control the disease. Geffen et al.  proposed a mathematical model of the hepatitis B virus, including isolation, treatment, and vaccine technology. Alrabaiah et al.  used the Galerkin method to solve the HIV infection model. They used a method called residual correction. The purpose of this technique is to reduce the error rate of the solution. Yüzbaşı et al.  used the cGP (2) and “LWCM” to approximate the solution of the proposed model. Furthermore, they solved the model using the traditional RK4-method. Finally, they compared the results obtained from the RK4-method to those acquired from the proposed schemes in order to ensure their validity. Sohaib  came up with a completely new way to think about the HIV pandemic. This model allows for a lot of new people to get infected. They examined the impact of public health education initiatives on the prevalence of the condition and found that they had no effect. In order to define the control and determine the best system, they employed “Pontryagin’s maximal principle”. Seatlhodi  developed and tested HIV/AIDS models using Caputo-fractional derivatives as a medical therapy. They first established Caputo-fractional order HIV/AIDS models with switching parameters and studied their dynamics using the Lyapunov–Razumikhin approach, based on the fractional derivative order linked to memory and genetic effects, and considering that the model coefficients are time-varying parameters. Wang et al.  proposed an SIR model with long-range temporal memory. The proposed model consists of delayed differential equations. They considered that the susceptible individual is following the logistic form, in which the incidence term takes the saturated form. The existence of steady states and the stability of those states are also examined. The concept of the Lyapunov function is used to figure out a new set of conditions that keep the steady states stable. The dynamical behavior of various infectious diseases are described using the idea of mathematical modeling (see  for detail information). Attaullah et al.  established a mathematical model for the dynamics of Human Immunodeficiency Virus (HIV) infection. They implemented the continuous Galerkin Petrov time discretization scheme and a fourth-order Runge-Kutta (RK4)-method to illustrate the dynamical behavior of the model, as well as a detailed description of the effects of different physical parameters of interest, which are depicted graphically and discussed how the level of healthy, infected CD T-cells, and free HIV particles varies related to the emerging parameters in the model. Sabir et al.  considered a novel designed prevention class in the HIV nonlinear model and solved numerically. Amin et al.  used the Haar wavelet approach to estimate the solution of the mathematical model of HIV infection CD4+T-Cells.

In this manuscript, we implemented a new method, namely the continuous Galerkin-Petrov scheme for finding the approximate solution of the non-linear model for HIV infection presented by Mehdi et al. . The proposed model is divided into four different compartments. We presented the impact of saturated and bilinear incidence functions and different clinical parameters (the parameter φ shows the default treatment rate for the individual, ζ is concerned with HIV treatment rate for the individual) on the dynamical behavior of the model. The present findings are conducive in the field of mathematical modeling of HIV infection of CD4+T-cells. This will be used to analyzed the population dynamics of CD4+T-cells in the presence and absence of HIV, helpful to observe the symptoms of AIDS seen clinically and valuable to hold back the disease. Medically, it provides sufficient information to clinicians to reduce the viral load of the infection.

Mathematical Description of the Model

In the absence of HIV, it is important to understand the population of T-cells produced by the bone marrow. Therefore, the premature cells shifted to an organ called the thymus, which is present in the chest sternum for further maturation and conversion into immune component T-cells. In humans, at the time of puberty, the thymus secretion in adults has minimal consequences, despite the thymus being in full operation and the fever lymphocytes performing as precursors of T-cells and immune component T-cells. The progression chain can be calculated by the number of T-cells, which shows us the initial symptoms. Enormous models have been developed for HIV infection. Mehdi et al.  suggested that the HIV infection model consists of four variables as follows: dSdt=λγS(t)ψ(S(t),I(t))I(t), dIdt=ψ(S(t),I(t))I(t)(φ+ζ+γ)I(t)+εA(t)+σC(t), dCdt=ζI(t)(σ+γ)C(t), dAdt=φI(t)(ε+γ+η)A(t)initial conditions are given as follows: S(0)=S00,I(0)=I00,C(0)=C00,A(0)=A00,the unknows S(t),I(t),C(t),A(t) illustrate the susceptible individuals, HIV infected individuals (with no clinical symptoms of AIDS), HIV infected individuals (under ART treatment with a viral load remaining low) and HIV infected individuals, respectively. The detailed description of the parameters used in the model are presented in Table 1.

Explanation of parameters involved in HIV infection model [<xref ref-type="bibr" rid="ref-1">1</xref>]. Unit: <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:msup><mml:mrow><mml:mtext>mm</mml:mtext></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mtext>days</mml:mtext></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>
Parameters Meaning Values
N(0) Initial population 23023935
So Population at risk 0.999
Io HIV-positive population 8.68×108
Co HIV-positive population (under ART treatment) 0
Ao Clinical signs HIV infected people with AIDS 140×107
λ Recruitment rate 2.19γ
β The rate of HIV infection 0.75
ζ The rate of HIV infection for I individuals 1.0
φ Default treatment rate for I individuals 0.1
η AIDS treatment rate 0.33
σ Default treatment rate of C individual 0.09
ε AIDS induced death rate 1.0

The general incidence function ψ(S,I) which is assumed to depend on the effective contact rate ρ>0. Then,ψ can take many forms. Table 2 is composed of some special incidence functions.

Some special incidence functions
Incidence functions ψ(S,I) References
Bilinear ρS 
Saturated ρS1+ε1S or ρS1+ε2I [22,23]
The Continuous Galerkin Petrov Technique

The Galerkin technique is an effective tool for numerically investigating critical challenges. This approach is commonly employed for complicated problems and is capable of dealing with nonlinear system and complicated problems.

This section is focused on the application and implementation of the suggested technique to the aforementioned model. For simplicity some assumptions are given, i.e., u1(t)=S(t),u2(t)=I(t),u3(t)=C(t),u4(t)=A(t), initially at t=0, u1(0)=S(0)=ρ1,u2(0)=I(0)=ρ2,u3(0)=C(0)=ρ3,u4(0)=A(0)=ρ4.

Find u:J=[0,T]V. Here J = [0,T] describes the time interval, for function

u:J×V and tT dtu(t)=ψ(t,u(t))tJ=[0,T]u(0)=u0,

Here dtu(t) refers to the time derivative of u(t). The derivative of u(t) w.r.t tu(t)= (u1(0),u2(0),u3(0))V represents u(t) at t = 0, and ψ=(ψ1,ψ2,ψ3) and describes as ψ:J×VV. The formulation of Eq. (6) is: find uX such that u(0)=uo and Jdtu(t),ϑ(t)dt=Jψ(t,u(t),ϑ(t)dtforallϑY,where X, and Y represent the solution space and the test space respectively. To describe function tu(t), we consider the space E(,V)=EO(J,V) as the space of contionuous functions u:JV equipped with norm ||u||E(,V)=EO(,V)=supt∈=[0,T]||u(t)||V.

M2(J,V) represents the space of square integrable function M2(,V) by containing discontinuous functions, which is expressed in the form as M2(,V)={u:[0,T]V:||u||M2(,V)<} with ||u||M2(,V)=(J||u(t)||V2dt)12

We divide the time interval J into N subintervals for the Galerkin time discretization. Jn=[tn1,tn], where n=1,2,3,,N, and 0 = t0<t1<tN1<tn=T. The symbol τ denotes the time discretization parameter, which will be used to determine the maximum time step size τ=max1nNτn, where τn=tntn1, which is the length of each Jn. Now we will approximate u:JV using a function uτ:JV (see [14,26,27,30,31] for details). Then, we will find the space for Xτl={uE(JV):u|JnHl(Jn,V)for all JnGτ},where Hl(Jn,V)={u:JnV,u(t)=j=0lUjtj,tJn,UjV,forallj}, and the test space for uτ is Yτk, illustrated as: Yτl={VM2(J,V):V|JnHl1(Jn,V)for all JnGτ},where Yτl consists of piecewise polynomials of order l1, which are discontinuous at the end points of the time intervals. By taking a test function ϑτYτl and multiplying it by Eq. (6) , and integrate over J (see [14,26,27,30] for details). Find uτXτl such that uτ(0)=0 and Jdtuτ(t),ϑτ(t)dt=Jψ(t,uτ(t)),ϑτ(t)dtϑτYτl.

This discretization is called the exact cGP-technique of order l. (see [14,2231] for details). Now, to find u|JnHl(Jn,V) such that Jndtuτ(t),ϑφ(t)dt=Jnφ(t,uτ(t)),ϑφ(t)dtϑVandφH(Jn),with the initial condition uτ|Jn(tn1)=uτ|Jn1(tn1), for n2 and uτ|Jn(tn1)=u0 for n=1. As a general case of a nonlinear function ψ.,., we approximate the integral on the right-hand side of Eq. (9) by the (l+1)-points Gauss-Lobatto quadrature formula (see [14,2231] for details). Find

u|JnHl(Jn,V), such that uτ(tn1)=un1 j=0lwjdtuτ(tn,j)φ(tn,j)=j=0lwjψ(tn,j,uτ(tn,j))φ(tn,j)φHk1(Jn),where wj are the weights.

To determine uτ|Jn, we represent it by a polynomial ansatz uτ(t)=j=0lUnjn,j(t)tJnwhere the coefficient Unj is the elements of V and the real valued function n,jH(Jn) are the Lagrange basis functions with respect to (l+1) suitable nodal points tn,jJn satisfying the conditions. n,j(tn,j)=δi,j,i,j=0,1,2,,l,where δi,j is the Kronecker delta that is δi,j={1:ifi=j0:ifij.

For the choice of initial conditions, we set tn,0 = tn1, which implies that the initial conditions for Eq. (9) is given as Un0=uτ|Jn(tn1),ifn2, Un0=u0,ifn=1.

The other points tn,1,tn,2,,ttn,l are chosen as the quadrature points of the l-points Gauss-Lobatto formula on the time interval Jn. For representation (11), for dtuτ, we get dtuτ=j=0lUnjn,j(t),tJn,

Using Eq. (14) in Eq. (9), we get Jndtuτ(t),ϑφ(t)dt=Jnj=0lUnj,ϑj(t)φ(t)dt.

This implies that Jndtuτ(t),ϑφ(t)dt=j=0lUnj,ϑJnj(t)φ(t)dt.

We define the basis functions n,jHk(Jn) via the affine reference transformation

ϖn: J^Jn, where J^=[1,1] and t=ϖnt^=tn+tn12+τn2t^Jnt^J^,n=1,2,3,,N.

Let ^jHk(J^), j=0,1,2,,l, denote the basis functions satisfying the conditions ^(t^i)=δi,j,i,j=0,1,2,,l,where t^0=1andt^i,i=1,2,3,,l, are the quadrature points for the interval J^. Then, we define the basis functions on the original time interval Jn by the mapping n,j(t)=^j(t^)witht^=ϖn1(t)=2τn(t+tn1tn2)J^.

Similarly, we define the test basis functions φn,i by suitable reference basis functions

φ^iHl1(J^), i.e., φn,i(t)=φ^i(ϖn1(t)tJn,i=1,2,3,,l.

Now, we transform the integral into a reference interval J^=[1,1] and approximate it by the (l+1) point Gauss-Lobatto quadrature formula which leads, for each test basis function φHl1 and for all ϑV, to the following: J^nj=0lUnj,ϑ^j(t^)φ^(t^)dt^=τn2J^nψ(ωn(t^),j=0lUnj(t^)),ϑφ^(t^)dt^ϑV.

This implies that μ=0lϖμ^j=0lUnj,ϑ^j(t^μ)φ^(t^μ)=τn2μ=0lϖμ^ψ(ωnt^μ,j=0lUnj(t^μ)),ϑφ^(t^μ).

Here ϖμ^ are the weights and t^μ [1, −1] are the integration points with the t^0=1 and t^l=1

φ^i(t^μ)=(ϖμ^)1δi,μi,μ={1,2,3,,l}.

Then we get the special form of the numerically integrated J^n-problem of cGP (k), find the l unknown coefficients UnjV where j=1,2,3,,l, such that j=olϖ^i,jUnj=τn2{ψ(tn,i,Unj)+σiψ(tn,0,Uno)}i=1,2,3,,lwhere Uno=Un1lforn>1andU1o=u0forn=1, indicated initial values and zi,j and σi are defined as ϖ^i,j=^j(t^μ)+σi^j(t^μ),tn,μ=ϖ^n(t^μ)andσi=ϖ^0φ^i(t^μ).

The cGP(2)-Method

Here, we apply the Gauss-Lobatto formula (Simpson rule) along the points tn,0=tn1, tn,1=(tn+tn12),tn,2=tn and the weights ϖ0^=ϖ2^=13, ϖ1^=43. Then we get the coefficients ϖ^i,j=(54114242),σi=(121),i=1,2j=0,1,2.with respect to the time interval J^n=[tn1,tn], we have to compute the system for the two unknowns ϖ^

Unj=uτ(tn,j) with tn,j=ϖn(t^) for j=(1,2). The couple (2 × 2) block-system for Un1,Un2V, is as follows: ϖ^1,1Un1+ϖ^1,2Un2=ϖ^1,0Un0+τn2{ψ(tn,1,Un1)+σ1ψ(tn,0,Uno)}, ϖ^2,1Un1+ϖ^2,2Un2=ϖ^2,0Un0+τn2{ψ(tn,2,Un2)+σ2ψ(tn,0,Uno)},

Un0 indicates the initial value at the time interval J^n obtaining from the time interval or the initial value u0.

Results and Discussions

In this study, we consider the global dynamics of a HIV infection model with general incidence rate. The aforesaid model represents the dynamical behavior of four different compartments. In the given model, we insert two types of incidence function to observe distinct variations through the graphical representation. Figs. 14 show the dynamical behavior for the bilinear incidence function ψ1(S,I) and saturated incidence function ψ2(S,I) initially. In Fig. 1, ψ1(S,I) and ψ2(S,I) show graphical results at time t, but after some intervals of time, we observe a sharp decrease in ψ1(S,I) and a gradual decrease in ψ2(S,I). In Fig. 2, at the initial stage, it is noticed that ψ1(S,I) and ψ2(S,I) are in dynamical equilibrium. However, a sudden rise occurs in the graph of ψ1(S,I), while gradual rises are observed in ψ2(S,I). Fig. 3 shows the graphical results for HIV infected individuals C(t). At the start, the curves for ψ1(S,I) and ψ2(S,I) having the same behavior and after some time a rapid increase occurs in ψ1(S,I) and a slow increase in ψ2(S,I). In Fig. 4, by increasing the concentration of HIV infected individuals A(t), the curves for both functions remain the same for some time, but after increase in time, the graph of ψ1(S,I) reaches the top with higher rate, while the graph of ψ2(S,I) reaches with lower rate. Figs. 58 present the bilinear incidence function ρS to illustrate the graphical results for S(t),C(t),I(t), and A(t), respectively. It could be seen in Fig. 5 that by varying the parameter values of (default treatment rate), the population of S(t) could be observed with no change at the beginning, but with the increase in default rate, a significant variation is observed in the population of S(t). Fig. 6 depicts that the increase in default treatment rate, i.e., I(t) population exhibits slow to rapid behavior. In Fig. 7, the graphical simulation of bilinear incidence function ρS for HIV infected individuals C(t) by varing φ have been examined. The population of C(t) is not affected initially, but after changing the value of φ from 0.1 to 0.4, a gradual influence occurs in the corresponding curves. The same variation in φ is considered for HIV infected individuals A(t) in Fig. 8. Figs. 912 show the influence of the saturated incidence function ρS1+ε1S for HIV susceptible individuals S(t), HIV infected individuals I(t), HIV infected individuals C(t), HIV infected individuals A(t), respectively. The effect of saturated incidence functions on the population dynamics of HIV susceptible individual S(t) are depicted in Fig. 9.

Population dynamics of the HIV susceptible individuals <inline-formula id="ieqn-189"><mml:math id="mml-ieqn-189"><mml:mi>S</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>Population dynamics of the HIV infected individuals <inline-formula id="ieqn-190"><mml:math id="mml-ieqn-190"><mml:mi>I</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>Population dynamics of HIV infected individuals <inline-formula id="ieqn-191"><mml:math id="mml-ieqn-191"><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>Population dynamics of the HIV infected individuals <inline-formula id="ieqn-192"><mml:math id="mml-ieqn-192"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>Graphical simulations of bilinear incidence function for <inline-formula id="ieqn-193"><mml:math id="mml-ieqn-193"><mml:mi>S</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-194"><mml:math id="mml-ieqn-194"><mml:mi>φ</mml:mi></mml:math></inline-formula>Graphical simulations of bilinear incidence function for <inline-formula id="ieqn-195"><mml:math id="mml-ieqn-195"><mml:mi>I</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-196"><mml:math id="mml-ieqn-196"><mml:mi>φ</mml:mi></mml:math></inline-formula>Graphical simulations of bilinear incidence function for <inline-formula id="ieqn-197"><mml:math id="mml-ieqn-197"><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-198"><mml:math id="mml-ieqn-198"><mml:mi>φ</mml:mi></mml:math></inline-formula>Graphical simulations of bilinear incidence function for <inline-formula id="ieqn-199"><mml:math id="mml-ieqn-199"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-200"><mml:math id="mml-ieqn-200"><mml:mi>φ</mml:mi></mml:math></inline-formula>Influence of saturated incidence function for <inline-formula id="ieqn-201"><mml:math id="mml-ieqn-201"><mml:mi>S</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-202"><mml:math id="mml-ieqn-202"><mml:mi>φ</mml:mi></mml:math></inline-formula>Influence of saturated incidence function for <inline-formula id="ieqn-203"><mml:math id="mml-ieqn-203"><mml:mi>I</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-204"><mml:math id="mml-ieqn-204"><mml:mi>φ</mml:mi></mml:math></inline-formula>Influence of saturated incidence function for <inline-formula id="ieqn-205"><mml:math id="mml-ieqn-205"><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-206"><mml:math id="mml-ieqn-206"><mml:mi>φ</mml:mi></mml:math></inline-formula>Influence of saturated incidence function for <inline-formula id="ieqn-207"><mml:math id="mml-ieqn-207"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-208"><mml:math id="mml-ieqn-208"><mml:mi>φ</mml:mi></mml:math></inline-formula>

The graph shows that raising the value of φ raises the susceptible population consistently. However, as time passes, the concentration of this population falls. The same phenomenon is observed for HIV infected individuals I(t), HIV infected individuals C(t), HIV infected individuals A(t) in Figs. 1012. Figs. 1316 demonstrate the impact of bilinear incidence function for S(t),I(t),C(t), and A(t) populations. From the graph, it is concluded that when the value of parameter ζ (HIV treatment rate) increases, the strength of population also increases. Initially, the susceptible population exhibits almost similar behavior; but, when the value of ζ is increased, a diminution in the population graphing could be noticed. The Fig. 13 displays the impact of HIV treatment rate on healthy population S(t). In the beginning, the concentration level of S(t) increases in uniform succession. However, increasing the HIV treatment rate ζ, a decline is observed in the graph of healthy cells population. Fig. 14 represents a minor change in the concentration of I(t), but after some time, the growth rate of the infected population increases by varying ζ. The similar behavior is stated in the Figs. 15 and 16 for C(t), and A(t), respectively.

Impact of bilinear incidence function for <inline-formula id="ieqn-209"><mml:math id="mml-ieqn-209"><mml:mi>S</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-210"><mml:math id="mml-ieqn-210"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Impact of bilinear incidence function <inline-formula id="ieqn-211"><mml:math id="mml-ieqn-211"><mml:mrow><mml:mi mathvariant="normal">ρ</mml:mi></mml:mrow><mml:mrow><mml:mtext>S</mml:mtext></mml:mrow></mml:math></inline-formula> for <inline-formula id="ieqn-212"><mml:math id="mml-ieqn-212"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-213"><mml:math id="mml-ieqn-213"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Impact of bilinear incidence function for <inline-formula id="ieqn-214"><mml:math id="mml-ieqn-214"><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-215"><mml:math id="mml-ieqn-215"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Impact of bilinear incidence function <inline-formula id="ieqn-216"><mml:math id="mml-ieqn-216"><mml:mrow><mml:mi mathvariant="normal">ρ</mml:mi></mml:mrow><mml:mrow><mml:mtext>S</mml:mtext></mml:mrow></mml:math></inline-formula> for <inline-formula id="ieqn-217"><mml:math id="mml-ieqn-217"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-218"><mml:math id="mml-ieqn-218"><mml:mi>ζ</mml:mi></mml:math></inline-formula>

Figs. 1720 demonstrate the dependence of the saturated incidence function ρS1+ε1S on each population one by one. The curves for the susceptible population follow the uniform trajectory initially in Fig. 17. However, after changing the parameter values of ζ from 1.1to1.4, a regular decrease is noticed in the graph. Fig. 18 states that the infected population I(t) is initially in decline, but after varying the HIV treatment rate ζ a regular increase is noticed in the curves. The same behavior is observed in the concentration of the C(t)andA(t) individuals.

Dependence of saturated incidence function <inline-formula id="ieqn-225"><mml:math id="mml-ieqn-225"><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="normal">ρ</mml:mi></mml:mrow><mml:mrow><mml:mtext>S</mml:mtext></mml:mrow></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msub><mml:mrow><mml:mi mathvariant="normal">ε</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mtext>S</mml:mtext></mml:mrow></mml:mrow></mml:mfrac></mml:mstyle></mml:math></inline-formula> for <inline-formula id="ieqn-226"><mml:math id="mml-ieqn-226"><mml:mi>S</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-227"><mml:math id="mml-ieqn-227"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Dependence of saturated incidence function for <inline-formula id="ieqn-228"><mml:math id="mml-ieqn-228"><mml:mi>I</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-229"><mml:math id="mml-ieqn-229"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Dependences of saturated incidence function for <inline-formula id="ieqn-230"><mml:math id="mml-ieqn-230"><mml:mi>C</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-231"><mml:math id="mml-ieqn-231"><mml:mi>ζ</mml:mi></mml:math></inline-formula>Dependences of saturated incidence function for <inline-formula id="ieqn-232"><mml:math id="mml-ieqn-232"><mml:mi>A</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> by varying <inline-formula id="ieqn-233"><mml:math id="mml-ieqn-233"><mml:mi>ζ</mml:mi></mml:math></inline-formula>
Mathematical Formulation

Consider an infinite horizontal parallel disk in which incompressible laminar flow of the nanomaterial fluid is examined. The plates are porous and rotating with angular velocity Ω=(0,Ω1,0) along the yaxis. The inspiration of Hall, and radiation is also considered. The nan-fluid model characterizes the Brownian motion and thermophoresis. The framework of the Cartesian coordinate system has been used to determine the problem’s geometry, in which the xaxis is analogous to the platters, the yaxis is the normal, and the zaxis is oblique to the xy-plane. The flow is represented in Fig. 21 as a schematic diagram. When a consistent magnetic flux is pragmatic, the fluid becomes electrically conduct. The lower plate is elongated linearly uw=cx in the x-direction at y=0, while the top plate is at y = h. The lower plate suctions fluid at v=v0(v0>0 correlates to suction, while v0<0 refers to injection). The imprint of the temperature-thermal conductivity, varied heat generation/absorption mixed with chemical response, and stimulation energy enhance heat and mass transmission. The model equations of nano-fluid stream are :

Geometry of the problem

The basic flow equations of Nano liquid are : xu~+yv~=0, u~xu~+v~yu~+2Ω1w~=1ρxρ+ν(xxu~+yyu~)+σ1B02ρ(1+m2)(u~mw~), u~xu~+v~yv~=1ρyρ+ν(xxv~+yyv~), u~xw~+v~yw~2Ω1u~=ν(xxw~+yyw~)σ1B02ρ(1+m2)(mu~+w~), u~xT~+v~yT~=1ρCPy[k(T)yT~]+τ[DByC~yT~+DTT~l(YT~)2]1ρcpyqr+K(t)uwxvρcp[D(T~wT~l)f+H(T~T~l)], u~xC~+v~yC~=DB[(xxC~+yyC~]+DTT~l(xxT~+yyT~)kr2(T~T~l)exp(EakT~)

With bounded conditions [34,36,37,40] u~|y=0=uw=cx,v~|y=0=v0,w~|y=0=0,T~|y=0=T~w,C~|y=0=C~wu~|y=h=0,v~|y=h=0,w~|y=h=0,T~|y=h=T~l,C~|y=h=C~l.

The radioactive heat flux is expressed by the following relation [34,35]: qr=43ϱ¯k¯yT4,whereT4=4Tl3T3Tl4

The varying thermal conductivity given in Eq. (29) is explored as [37,39,41,42]: K=k0(1+d(TTlTwTl))

By using Eqs. (32) and (33), Eq. (29) becomes u~xT~+v~yT~=1ρCPy[k(t)yT~]+τ[DByC~yT~+DTT~l(YT~)2]+163ϱ¯k¯yyT+k0(1+dθ)uwxvρcρ[D(T~wT~l)f+H(T~T~l)],

By proper conversion [34,35,38] u~=cxf(η),v~=chf(η),w~=cxj(η),θ(ζ)=T~T~lT~wT~l,(η)=C~C~lC~wC~l,ζ=YH.

By using the upstairs conversion Eq. (25) is triflingly equated. However, Eqs. (26)(28), (30), (31), and (34) yield the system d4fdη4=Re(dfdηd2fdη2fd3fdη3)+2α1djdηHa2(1+m2)(d2fdη2mdjdη), d2jdη2=Re(fdjdηfdfdηj)2α1dfdη+Ha2(1+m2)(j+mdfdη), ((1+dθ)+43Rd)Re)d2θdη2=Pr(fdθdη+Re(Nbdθdηdϕdη+Nt(dθdη)2))(1+dθ)(Ddfdη+Hdθdη)Red(dθdη)2, d2dη2=NtNbd2θdη2+Sc(δ(1+αθ)nexp(E(1+αθ))f.Redϕdη).

The boundary constraints gross the form: Lower plate:dfdη(0)=1,f(0)=K,j(0)=0,θ(0)=1,ϕ(0)=1, Upper plate:dfdη(1)=0,f(1)=0,j(1)=0,θ(1)=0,ϕ(1)=0.

The physical quantities of interest are as follows: Cf,lower=μρuw2yu|y=0 Cf,upper=μρuw2yu|y=h Nulower=hQwk0(TwTl),Qw=k(t)yT+qr|y=0 Nuupper=hQwk0(TwTl),Qw=k(t)yT+qr|y=h Shlower=hQmDB(CwCl),Qm=DByC|y=0 Shlower=hQmDB(CwCl),Qm=DByC|y=h.

Using Eqs. (35), (41)(46) are transformed into (RehCf)flower=d2fdη2|ζ=0,(RehCf)upper=d2fdη2|η=1 (Nu)lower=(1+(43).Rd1+dθ)dθdη|η=0,(Nu)upper=(1+(43).Rd1+dθ)dθdη|η=1 (Sh)lower=dφdη|η=0,(Sh)upper=dφdη|η=1.

Numerical Solutions

This segment is devoted to enclosing a well-known Galerkin scheme to handle the aforementioned nonlinear problems. For the validity of the obtained solution, we compared the present solution with those exist in the published literature. Table 3 provides a comparative analysis of the current study with Tlili et al. . There is a significant correlation between the outcomes.

Validation of the present work with published work for temperature and concentration profiles
η Tili et al.  Current Tili et al.  Galerkin scheme
θ(η) ϕ(η)
0 1 1 1 1
0.1 0.7568230 0.7568230 0.6160560 0.6160560
0.3 0.6071153 0.6071151 0.3736382 0.3736381
0.5 0.4046126 0.4046123 0.1671527 0.1671527
0.7 0.1718404 0.1718401 0.0026777 0.0026778
1 0 0 0 0
Conclusions

In this paper, we considered the HIV infection model, which consists of four nonlinear ordinary differential equations. We applied an innovative numerical approach known as the continuous Galerkin-Petrov scheme to determine the solution of the model. In addition, we analyzed the dynamics of HIV-infected model with a different incidence rate. Assessed the effects of various clinical parameters on the dynamical behavior of distinct compartments. In the suggested model, we included two types of incidence functions (bilinear and saturated incidence functions) in order to observe visually distinctive fluctuations. By varying the values of various parameters, we observed the periodic rise and fall of the curves of various populations. The bilinear incidence function and the saturated incidence function initially exhibit identical dynamical behavior, as indicated in the illustrations. Nevertheless, with time, various graphical representations evolve. The aforementioned results highlight the importance for mathematical modelling of HIV infection. This will be performed to analyze the population dynamics of CD4+T-cells in the existence and exclusion of HIV, which will be beneficial in identifying clinical AIDS manifestations and in halting the epidemic. It enables physicians with enough information to minimize the viral burden of the disease. The aforementioned approach was employed to investigate a mathematical model for nano-material fluid flow between two analogous infinite disks. The findings are validated through comparison toward those reported in the literature.

Future Recommendations:

It is a well-known observation that fractional analysis has increasingly become a prominent research area. It has been demonstrated that fractional calculus is especially beneficial for imitating a number of legitimate situations. Employing fractional order derivatives and integrals, researchers have evaluated infectious maladies such as COVID-19, HIV, AIDS, and others. The future challenge will involve assessing the quantitative and qualitative aspects of our concept with various fractional order derivatives.

The authors would like to thank the editor and anonymous reviewers.

Funding Statement

The authors received no specific funding for this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest to report regarding the present study.