Ravi Shanker Dubey,Dumitru Baleanu,Manvendra Narayan Mishra,★and Pranay Goswami
1Department of Mathematics,AMITY School of Applied Science,AMITY University Rajasthan,Jaipur,302022,India
2Department of Mathematics and Computer Sciences,Cankaya University,Balgat,06530,Turkey
3Institute of Space Sciences,Magurele-Bucharest,R76900,Romania
4School of Liberal Studies,Dr.B.R.Ambedkar University Delhi,Delhi,110006,India
ABSTRACT Diabetes is a burning issue in the whole world.It is the imbalance between body glucose and insulin.The study of this imbalance is very much needed from a research point of view.For this reason,Bergman gave an important model named-Bergman minimal model.In the present work,using Caputo-Fabrizio(CF)fractional derivative,we generalize Bergman’s minimal blood glucose-insulin model.Further,we modify the old model by including one more component known as diet D(t),which is also essential for the blood glucose model.We solve the modified model withthe help of Sumudu transform and fixed-point iteration procedures.Also,using the fixed point theorem,we examine the existence and uniqueness of the results along with their numerical and graphical representation.Furthermore,the comparison between the values of parameters obtained by calculating different values of t with experimental data is also studied.Finally,we draw the graphs of G(t),X(t),I(t),and D(t)for different values of τ.It is also clear from the obtained results and their graphical representation that the obtained results of modified Bergman’s minimal model are better than Bergman’s model.
KEYWORDS Bergman minimal model; blood glucose; Caputo-Fabrizio fractional derivative; uniqueness and existence;fractional calculus
These days, Mathematical modeling is turning out to be a vital tool in mathematical science.Since it translates real-world problems into mathematical language and after applying the necessary methods, we again translate the results into real-world languages to forecast the objective.In every field, modeling is being used to achieve or predict the future prospective.In recent years, modeling is being used with another important but less interactive field of mathematics—Fractional Calculus (see [1-8]).We have studied several problems and their solutions by ordinary calculus methods, but sometimes fractional calculus gives us better results to describe the model than the old one.Fractional calculus has several real-world applications.A few of them are Fractional conservation of mass, Groundwater flow problem, Time-space fractional diffusion equation models, Acoustical wave equations, Fractional Schrdinger equation in quantum theory, etc.(see [9-22]).
Nowadays, diabetes is turning out to be a fatal disease.The imbalance between body glucose and insulin causes diabetes.It is categorized into two types.Type-I diabetes is more severe than Type-II.According to a survey, Type-I and Type-II diabetes sufferers are one ratio nine.So, the study of this imbalance was very much needed from a research point of view.For this reason,Bergman gave an important model named-Bergman minimal model (see [23-27]).
Riemann-Liouville and Caputo introduced the initial concept of fractional calculus.But they use a singular kernel in their definition.However, Caputo and Fabrizio recently found specific systems related to material heterogeneities that cannot be well connected with Riemann-Liouville or Caputo derivative.Because of that, Caputo and Fabrizio introduced a new fractional derivative operator involving a non-singular kernel.Due to its memory and non-singularity property, this operator is used to study different models of engineering, science, and bio-mathematical fields.In general, we get more reliable results while using this operator than others for more details of applications of Caputo-Fabrizio derivative (see [28-45]).
Apart from changing the model to fractional-order, we also include the Diet factor in the model that describes the effect of meals on the glucose level.The generalize of Bergman’s minimal blood glucose-insulin model gives us a more precise and detailed prediction about the problem.
Definition 1:Supposeh∈H1(a1,b1),b1>a1,β∈[0,1] hence the Caputo-Fabrizio differential coefficient of fractional order:
M(β)is function of normalization such asM(0)=M(1)=1.
If,h/∈H1(a1,b1), then the derivative is
also,
One thing is to be noted here that the Caputo-Fabrizio operator has an exponent differential coefficient.We know that the exponential function has no singularity which, means the derivative exists at every point of the territory.So CF operator gives better results in its domain, see [49-56].
Definition 2:The fractional integral of the function h(t) of orderβ (0<β <1), is
Remark 1:From the above equation, we conclude that the fractional integral of orderβ(0<β <1), is the mean of h and its anti-derivative.So Nieto et al.gave this condition,
Based on the above relation, Nieto and Losada defined the following operator;
Supposef(t)be a function whose Caputo-Fabrizio derivative occurs, hence the Sumudu Transform of Caputo-Fabrizio fractional differential coefficient of f(t) is defined below:
The whole paper is divided into six sections.Section first introduces fractional calculus and mathematical modeling with its brief history.The second section describes the Bergman model’s fractional exponent and an overview of the fractional modified minimal model.The third section discusses the existence and uniqueness of the modified mathematical model with the help of the Caputo-Fabrizio operator.Section 4 is devoted to the solution of the model by using the Sumudu Transform operator.Section 5 deals with numerical solutions along with graphical representations of the modified Bergman model.In section 6, we have concluded our findings.
Here, we are going to discuss the Bergman model of fractional order.In this model, we consider a glucose chamber and plasma insulin which is inherent to carry out across the isolated chamber to control the final glucose intake.This model, asserts that this is good to satisfy specific validation criteria with minimum parameters.This system is efficient in describing the gesture of interactivity of blood sugar and insulin.In the presented model, we suppose thatG(t)is the imbalance of plasma glucose cluster andI(t)is the free plasma insulin cluster, from their initial values.The model was
with starting conditionsG(0)=G0,X(0)=X0, andI(0)=I0.
Here, we redefined the structure of the previous model.We made some changes in the model and added few parameters.As we have, ‘Diet’having significant effect on the glucose level.Now the changed structured Bergman modified system has been described in the following equations-
Here one thing is to be noted that 0<τ≤1 and initial restrictions are,G(0)=G0,X(0)=X0,I(0)=I0andD(0)=D0.
This model can be a tool in search of an artificial pancreas.But, unfortunately, it also adopts the problems with the glucose-minimal model.Here,G(t)—blood glucose cluster,X(t)—aftermath of effective insulin,I(t)—blood insulin cluster,D(t)—infusion of exogenous glucose,u(t)—insulin distribution function,Gb—initial blood glucose cluster,Ib—initial blood insulin cluster,V1—insulin issuance volume,n—fractional fading rate of insulin,p1—insulin free glucose dispensation rate,p2—active insulin dispensation rate andp3—the increase in absorption capacity caused by insulin.
DefineK1,K2,K3,K4and their relations with variables.
Proof:Since the system is
Here one thing is to be noted that 0<τ≤1.Then converting the above system into a system of integral equations
Then by definition stated by Nieto, we get
Now let us suppose the kernels are given as
Show thatK1,K2,K3andK4satisfies the Lipchitz condition.
Proof:At first, we shall show this forK1.SupposeGandG1are any functions, so
Similarly, we can have
and
Now consider the recursive formula
Now suppose that the deviation amid two successive terms is
Now
ButK1satisfies Lipchitz condition so,
Similarly, we have
Establish that Bergman Minimal Model with Fractional-order is a minimum system of sugar insulin dynamics.
Proof:Using the recursive technique, we have
Hence the existence of results is verified, which is continuous too.Now we get
wherePn,Qn,RnandSnare residues of series solution.
So,
Similarly, we have
From the above, it is clear that,
Now,
Now takingn→∞we have
On behalf of the above equations, we can state that the solution of the system exists.
To show the uniqueness of results, we suppose that other sets of results exist for the system specified from Eqs.(46) to (49) such that
taking norm both sides, we get
using Lipchitz condition, we obtain
which is valid for alln, so
similarly
Hence, it claims uniqueness of the system.
Considering the system has various equations so it can be challenging to find the exact results.For this, we will adopt the iterative technique together with the Sumudu Transform.Now take Sumudu Transform either sides, side of Eq.(15)
or
At this moment, take inverse Sumudu Transform on both ends, we have-
Similarly, we have from Eqs.(16)-(18)
Then we get the following recurrent form from the above
The result is obtained by
and
For numerical solution, we use the values given below by the experiment defined in [59,60].
Table 1:Table with initial value and parameters
Table 1(continued)
Table 2:Comparison between the values of G for different values of t with experimental data
Table 3:Comparison between the values of I for different values of t with experimental data
By using the above-defined values, we can easily determine the numerical results for the defined model.In that analysis, we found the numerical result by using the sumudu transform.We found the result for some differentτvalues and formed a table that shows the comparison of the obtained result with the experimental data (see Tab.1).In Tab.2, we discussed G, and in Tab.3,we discussed I and found from these analyses that the fractional model gives less error than the integer model.We also drew four figures to show the numerical outcomes.In Fig.1 we showed the values ofGat different values ofτ, i.e.,τ=1,τ=0.95,τ=0.9,τ=0.8, andτ=0.7.In Fig.2, we showed the values ofIat different values ofτ.In the same way, we drew the Figs.3 and 4 to show the numerical results of X and D at different values ofτ.From the numerical outcome, it is clear that the diet is an essential component of glucose level.
Figure 1:Graph of blood glucose cluster (G) with respect to time t for different values of τ
Figure 2:Graph of blood insulin cluster (I) with respect to time t for different values of τ
Furthermore, it has been observed from Tabs.2 and 3, that the values obtained for G and I with linear order have more error than fractional order.So we can see that the fractional model with Caputo-Fabrizio operator gives better results in contrast to linear model, and our model defined the real-world problem in a better manner.
Figure 3:Graph of aftermath of effective insulin (X) with respect to time t for different values of τ
Figure 4:Graph of infusion of exogenous glucose (D) with respect to time t for different values of τ
The presented work strives to explain the presence and oneness of the Modified Bergman Minimal Model that is stretched out by Caputo-Fabrizio fractional differential coefficient in the frame of reference of sugar and insulin quantity in blood.From the obtained results, it is clear that the fractional model error reduces compared to integer order.Thus, we get nearby results of a set-up that displays the aftermath of time on the concentrationsG(t),X(t),I(t), andD(t).As in future work, we can generalize the model or compare the results using various differential operators.
Acknowledgement:The authors wish to express their appreciation to the reviewers for their helpful suggestions which greatly improved the presentation of this paper.
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.
Computer Modeling In Engineering&Sciences2021年9期