Loading [MathJax]/extensions/TeX/boldsymbol.js
Research article

Exploring analytical results for (2+1) dimensional breaking soliton equation and stochastic fractional Broer-Kaup system

  • Received: 04 February 2024 Revised: 13 March 2024 Accepted: 14 March 2024 Published: 25 March 2024
  • MSC : 34A05, 35C08, 60G22

  • This paper introduces a pioneering exploration of the stochastic (2+1) dimensional breaking soliton equation (SBSE) and the stochastic fractional Broer-Kaup system (SFBK), employing the first integral method to uncover explicit solutions, including trigonometric, exponential, hyperbolic, and solitary wave solutions. Despite the extensive application of the Broer-Kaup model in tsunami wave analysis and plasma physics, existing literature has largely overlooked the complexity introduced by stochastic elements and fractional dimensions. Our study fills this critical gap by extending the traditional Broer-Kaup equations through the lens of stochastic forces, thereby offering a more comprehensive framework for analyzing hydrodynamic wave models. The novelty of our approach lies in the detailed investigation of the SBSE and SFBK equations, providing new insights into the behavior of shallow water waves under the influence of randomness. This work not only advances theoretical understanding but also enhances practical analysis capabilities by illustrating the effects of noise on wave propagation. Utilizing MATLAB for visual representation, we demonstrate the efficiency and flexibility of our method in addressing these sophisticated physical processes. The analytical solutions derived here mark a significant departure from previous findings, contributing novel perspectives to the field and paving the way for future research into complex wave dynamics.

    Citation: Faeza Lafta Hasan, Mohamed A. Abdoon, Rania Saadeh, Ahmad Qazza, Dalal Khalid Almutairi. Exploring analytical results for (2+1) dimensional breaking soliton equation and stochastic fractional Broer-Kaup system[J]. AIMS Mathematics, 2024, 9(5): 11622-11643. doi: 10.3934/math.2024570

    Related Papers:

    [1] Hasib Khan, Jehad Alzabut, Dumitru Baleanu, Ghada Alobaidi, Mutti-Ur Rehman . Existence of solutions and a numerical scheme for a generalized hybrid class of n-coupled modified ABC-fractional differential equations with an application. AIMS Mathematics, 2023, 8(3): 6609-6625. doi: 10.3934/math.2023334
    [2] Shakeel Muhammad, Obaid J. Algahtani, Sayed Saifullah, Amir Ali . Theoretical and numerical aspects of the Malaria transmission model with piecewise technique. AIMS Mathematics, 2023, 8(12): 28353-28375. doi: 10.3934/math.20231451
    [3] Rahat Zarin, Abdur Raouf, Amir Khan, Aeshah A. Raezah, Usa Wannasingha Humphries . Computational modeling of financial crime population dynamics under different fractional operators. AIMS Mathematics, 2023, 8(9): 20755-20789. doi: 10.3934/math.20231058
    [4] Elkhateeb S. Aly, Mohammed A. Almalahi, Khaled A. Aldwoah, Kamal Shah . Criteria of existence and stability of an n-coupled system of generalized Sturm-Liouville equations with a modified ABC fractional derivative and an application to the SEIR influenza epidemic model. AIMS Mathematics, 2024, 9(6): 14228-14252. doi: 10.3934/math.2024691
    [5] Nadeem Khan, Amjad Ali, Aman Ullah, Zareen A. Khan . Mathematical analysis of neurological disorder under fractional order derivative. AIMS Mathematics, 2023, 8(8): 18846-18865. doi: 10.3934/math.2023959
    [6] Mdi Begum Jeelani, Kamal Shah, Hussam Alrabaiah, Abeer S. Alnahdi . On a SEIR-type model of COVID-19 using piecewise and stochastic differential operators undertaking management strategies. AIMS Mathematics, 2023, 8(11): 27268-27290. doi: 10.3934/math.20231395
    [7] Muhammed Jamil, Rahmat Ali Khan, Kamal Shah, Bahaaeldin Abdalla, Thabet Abdeljawad . Application of a tripled fixed point theorem to investigate a nonlinear system of fractional order hybrid sequential integro-differential equations. AIMS Mathematics, 2022, 7(10): 18708-18728. doi: 10.3934/math.20221029
    [8] Hongyan Wang, Shaoping Jiang, Yudie Hu, Supaporn Lonapalawong . Analysis of drug-resistant tuberculosis in a two-patch environment using Caputo fractional-order modeling. AIMS Mathematics, 2024, 9(11): 32696-32733. doi: 10.3934/math.20241565
    [9] Saleh S. Redhwan, Maoan Han, Mohammed A. Almalahi, Maryam Ahmed Alyami, Mona Alsulami, Najla Alghamdi . Piecewise implicit coupled system under ABC fractional differential equations with variable order. AIMS Mathematics, 2024, 9(6): 15303-15324. doi: 10.3934/math.2024743
    [10] Gulalai, Shabir Ahmad, Fathalla Ali Rihan, Aman Ullah, Qasem M. Al-Mdallal, Ali Akgül . Nonlinear analysis of a nonlinear modified KdV equation under Atangana Baleanu Caputo derivative. AIMS Mathematics, 2022, 7(5): 7847-7865. doi: 10.3934/math.2022439
  • This paper introduces a pioneering exploration of the stochastic (2+1) dimensional breaking soliton equation (SBSE) and the stochastic fractional Broer-Kaup system (SFBK), employing the first integral method to uncover explicit solutions, including trigonometric, exponential, hyperbolic, and solitary wave solutions. Despite the extensive application of the Broer-Kaup model in tsunami wave analysis and plasma physics, existing literature has largely overlooked the complexity introduced by stochastic elements and fractional dimensions. Our study fills this critical gap by extending the traditional Broer-Kaup equations through the lens of stochastic forces, thereby offering a more comprehensive framework for analyzing hydrodynamic wave models. The novelty of our approach lies in the detailed investigation of the SBSE and SFBK equations, providing new insights into the behavior of shallow water waves under the influence of randomness. This work not only advances theoretical understanding but also enhances practical analysis capabilities by illustrating the effects of noise on wave propagation. Utilizing MATLAB for visual representation, we demonstrate the efficiency and flexibility of our method in addressing these sophisticated physical processes. The analytical solutions derived here mark a significant departure from previous findings, contributing novel perspectives to the field and paving the way for future research into complex wave dynamics.



    Dynamical modeling systems and their computational analysis have been the interest of researchers in science and engineering. Recently, fractional operators have been receiving a lot of attention for their ability to generalize the scope of the classical works. Singular and non-singular kernels and local and non-local kernels, are studied in fractional calculus (FC). Readers can explore the impressive works in the field of FC. We highlight some useful works as follows: a fractional order COVID-19 model and simulations [1], fractional order sine-Gordon equations [2], fractional order nonlinear models [3], and imperfect testing disease models [4,5], which can be studied in the literature.

    The infectious disease known as tuberculosis (TB) is typically brought on by the bacteria known as Mycobacterium TB. Although TB most frequently infects the lungs, the disease can manifest in other places of the body as well. Latent TB is the term used when an infection does not manifest any symptoms. This occurs in the majority of cases. About 10% are reinfected if treatment is not received, and approximately half of those who are afflicted will pass away due to the active disease that was caused by the latent infection that eventually became active. Active TB often presents itself with a cough that brings up blood-tinged mucus, night sweats, fever, and a reduction in body weight. Other symptoms may accompany these as well. Because of the significant weight loss that was associated with the disease in the past, it was historically known as consumption. There is a wide variety of symptoms that can be brought on by an infection in another organ. People who have active TB in their lungs can pass the disease on to others through the air when they cough, spit, speak, or sneeze. This is how the disease is transmitted. Latent TB does not contribute to the spread of the infection. Active infection happens more frequently in population who already have HIV or AIDS, or are smokers. Chest X-rays, as well as microscopic analysis and the culture of fluids in a body, are utilized in the diagnostic process for active tuberculosis. The tuberculin skin test, also known as the TST, or a blood is used to diagnose latent TB. One of the most interconnected important public health issues is TB. TB is a leading cause of deaths in HIV-positive individuals, according to [6]. Those who are infected with Mycobacterium TB and also have HIV or AIDS have a tenfold increased risk of developing TB compared to those who are HIV-uninfected [7]. Because of how the disease expresses itself in a unique way, TB is usually challenging to diagnose [8]. Additionally, the presence of TB makes HIV infection worse through a number of processes, some of which include an enhanced viral structure and an CD4/CD8 cell growth [9]. Because of the drug complications, close toxic effects characteristics and immune recovery inflammatory condition, and a high medication burden that may affect compliance, co-infection of TB and HIV, makes it difficult to handle. This is because of the similarities in toxicity traits between the two conditions.

    Moreover, there has been indication that prompt diagnosis of TB reduces the risk of disease and death in individuals who have been infected with the condition, and successful treatments are regularly used. One of the most significant challenges to global health that the world's nations have been confronted with over the passed 30 years has been the spread of infection of various groups that have the HIV-TB diseases. These illnesses are two of the leading causes of death around the world, and one of the reasons for this is that HIV infection weakens the immune system of humans.

    Even though the mathematical modeling of biological processes is not a novel concept in the realm of scientific fields, the COVID-19 pandemic which has significantly affected the entire world has enhanced its usage. This is the case despite the fact that the mathematical modeling of biological processes has been an interest of researchers for quite some time [10]. Co-infection of HIV and TB disorders is typical and had been mathematically estimated and examined by applied researchers in numerous works. As [11] shows, some of the research included the dynamics of HIV-TB co-infection as well as integrated treatment strategies. Using such a model, different combination therapy strategies were examined, simulating HIV-TB co-infection. They also tried combining different medicines that were given at different periods. Mycobacterium TB has been shown to worsen the clinical record in HIV infected individuals; ultimately, dealing with TB in HIV infected individuals may have a substantial impact on the disease spread. A model of a host with co-infections of TB and HIV was constructed by Kirschner et al. [12]. It was the first time anyone had attempted to understand how TB impacts the characteristics of HIV-infected people. It is believed that TB, which is the most prevalent and severe opportunistic infection that affects HIV-positive people, is the primary cause of more than half of the cases that develop into AIDS syndrome in underdeveloped countries.

    Zhao et al. [13], considered the following TB model with consideration of the age factor in the susceptible class:

    dS1dt=Ad1S1ϑ1S1Im1S1,dS2dt=m1S1ϑ2S2Id2S2m2S2,dS3dt=m2S2d3S3ϑ3S3I,dEdt=(1p)[ϑ1S1+ϑ2S2I+ϑ3S3I]vEdE,dIdt=p(ϑ1S1I+ϑ2S2I+ϑ3S3I)(d+γ+μ)I+ηR+vE,dRdt=γIdRηR, (1.1)

    with

    {S1(0)=S010,S2(0)=S020,S2(0)=S020,E(0)=E00,I(0)=I00,R(0)=R00.

    In this model there are four compartments including the susceptible group "S(t)", "E(t)" which represents the exposed but not yet actively infected individuals, "I(t)" which represents the infected people among the population and the recovered people which are grouped in the class "R(t)". For better interpretation of the age factor in the infection spread, the susceptible portion has been divided into three subclasses: child "S1(t)", middle-aged people "S2(t)", and elders "S3(t)". The death rate of the recovered people is denoted by "d". Death caused by TB is expressed in terms of the rate of "μ" which is added to "d". The annual birth in the whole population is "A", "m1" and "m2" are the rates at which the child susceptibles become middle age susceptibles and the susceptibles from the middle age group transition to the elderly, respectively, "ϑ1", "ϑ2" and "ϑ3", are the death rates for the child, middle-aged people and elderly susceptible groups, respectively "p" is the rate of fast-spreading infection cases; "v" is the rate at which the dormant organisms are re-activated in TB patients, d1,d2,d3 are the death rates for the child, the middle-aged and the subclass of the elderly susceptible, "d" is the natural death rate, "μ" is the rate of death due to TB, "γ" represents the rate of regaining health, and "η" is the rate of TB treatment.

    FC is used in so many fields of science and engineering because of its theoretical outcomes, numerical computations and experimental results. The use of FC-theory chaotic problems for the investigation of previously unknown and novel forms of chaotic behavior in dynamical systems can be seen in [14,15]. In mathematical physics, FC has numerous applications [16,17]. Bio-mathematics [18,19,20] and engineering [21,22] problems have recently become target applications for FC. Regarding dynamical systems, smoking models can be seen in previous works [23,24]. Non-oscillatory solutions for nabla forced equations were studied in [25]. Asymptotic and oscilatory criteria for nonlinear differential equations were considered in [26]. Impulsive fractional differential equations (FDEs) with piecewise Caputo's FDEs were investigated for the existence and uniqueness of solutions (EUS) with applications in [27,28,29].

    Further, Atangana and Araz studied a new subject of operators known as piecewise operators [30]. The exponential and Mittag-Lefler kernels fail to address the crossover behavior of a case problem. Scientists are utilizing this approach to analyze the crossover behaviors in dynamical systems. In the study of disease dynamics, this type of piecewise operator is extensively applied. For example, Dengue virus transmission within populations was considered for piecewise operators by Ahmad et al. [31]. One can see more related and useful works on piecewise operators in [32,33,34]. For more details on numerical methods and their applications the readers may study in [35,36].

    In light of these merits, we give the notion of the piecewise modified Atangana-Baleanu-Caputo (mABC) fractional derivative and apply it to a TB model. This novel operator is a combination of the classical derivative and the recently developed modified Atangana-Baleanu operator in the Caputo's sense. For this combination, we have considered the splitting of an interval [0,t2] for t2R+, such that, the classical derivative is applied in the first portion [0,t1] while the second differential operator is applied in the interval [t1,t2]. As a result, we have the piecewise mABC operator. Its corresponding integral is also given accordingly. This new operator is then applied to a TB model for the study of crossover behavior. The existence of solutions and stability of solutions are studied for the nonlinear piecewise mABC TB model. A numerical scheme for the simulations is presented and applied to the available data. We shall explore the model (1.1) in the following fashion, employing the piecewise differential operator in the case of the classical and mABC derivatives.

    PCC0DϖtS1=Am1S1d1S1ϑ1S1I,PCC0DϖtS2=m1S1m2S2d2S2ϑ2S2I,PCC0DϖtS3=m2S2d3S3ϑ3S3I,PCC0DϖtE=(1p)[ϑ1S1+ϑ2S2I+ϑ3S3I]vEdE,PCC0DϖtI=p(ϑ1S1I+ϑ2S2I+ϑ3S3I)(γ+μ+d)I+ηR+vE,PCC0DϖtR=γIdRηR. (1.2)

    In the piecewise TB model (1.2), the PCC represents piecewise the classical and mABC derivatives.

    In order to make the calculations easier, we rewrite TB model (1.2) in piecewise form, as follows

    PCC0Dϖt(S1)={Dϖt(S1(t))=G1(S1,S2,S3,E,I,R,t), 0<tt1,mABC0Dϖt(S1(t))=G1(S1,S2,S3,E,I,R,t),  t1<tT,,PCC0Dϖt(S2)={Dϖt(S2(t))=G2(S1,S2,S3,E,I,R,t), 0<tt1,mABC0Dϖt(S2(t))=G2(S1,S2,S3,E,I,R,t), t1<tT,,PCC0Dϖt(S3(t))={Dϖ0(S3(t))=G3(S1,S2,S3,E,I,R,t), 0<tt1,mABC0Dϖt(S3(t))=G3(S1,S2,S3,E,I,R,t),  t1<tT,,PCC0Dϖt(E(t))={Dϖt(E(t))=G4(S1,S2,S3,E,I,R,t), 0<tt1,mABC0Dϖt(E(t))=G4(S1,S2,S3,E,I,R,t),  t1<tT,,PCC0Dϖt(I(t))={Dϖt(I(t))=G5(S1,S2,S3,E,I,R,t),  0<tt1,mABC0Dϖ0(I(t))=G5(S1,S2,S3,E,I,R,t),  t1<tT,,PCC0Dϖt(R(t))={Dϖt(R(t))=G6(S1,S2,S3,E,I,R,t),  0<tt1,mABC0Dϖt(R(t))=G6(S1,S2,S3,E,I,R,t),  t1<tT, (1.3)

    where, G1=Ad1S1m1S1ϑ1S1I,G2=m1S1d2S2ϑ2S2Im2S2,G3=m2S2d3S3ϑ3S3I,G4=(1p)[ϑ1S1+ϑ2S2I+ϑ3S3I]vEdE,G5=p(ϑ1S1I+ϑ2S2I+ϑ3S3I)(d+γ+μ)I+ηR+vE,G6=γIdRηR.

    In this section, the fundamental definitions associated with the mABC calculus are explained. Additionally, the piecewise operators are defined with respect to both the classical and mABC operators. These definitions and findings will be put to use in the future to obtain existence results and implement a numerical scheme for the simulations of the dynamics of the piecewise model (1.2)

    Definition 2.1. [37,38] For ϱ(0,1), and fL1(0,T), the mABC derivative is given as follows

    mABCDϱ0f(t)=B(ϱ)1ϱ[f(t)Eϱ(μϱtϱ)f(0)μϱt0(ts)ϱ1Eϱ,ϱ(μϱ(ts)ϱ)f(s)ds].

    One can easily obtain that mABCDϱ0C=0 [37].

    Definition 2.2. [37,38] The modified AB-integral for ϱ(0,1), and fL1(0,T), is defined by

    mABIϱ0f(t)=B(1ϱ)B(ϱ)[f(t)f(0)]+[RLIϱ0(f(t)μϱf(0))]. (2.1)

    Lemma 2.3. [37] For fL1(0,), and ϱ(0,1), we have

    mABIϱ0mABCDϱ0f(t)=f(t)f(0). (2.2)

    Definition 2.4. [30] The piecewise integral for a differentiable h(t), is given by:

    PF0Ith(t)={t0h(s)ds,  0<tt1,tt1h(s)ds  t1<tt2.

    Definition 2.5. Consider a differentiable function h(t); then, the piecewise derivative (PD) of h(t) is

    PF0Dϖth(t)={h(t),  0<tt1,mABC0Dϖth(t)  t1<tt2.

    Definition 2.6. For a differentiable function h(t), the piecewise integration of the classical and mAB-integral is given by

    PF0Ith(t)={t0h(s)ds,  0<tt1,B(1ϖ)B(ϖ)[h(t)h(0)]+[RLIϖ0(h(t)μϖh(0))],  t1<tt2,

    where PF0Ith(t) is integer order integration in 0<tt1 and in t1<tt2; it is an mAB fractional integral.

    Lemma 2.7. The piecewise DE for t(0,T] is given by

    PFC0Dϖth(t)=G(t,h(t)),

    has the following solution

    h(t)={h0+t0h(s)ds,0<tt1,h0+1ϖB(ϖ)G(t,h(t))+ϖΓ(ϖ)B(ϖ)t0(ts)ϖ1G(s,h(s))ds1ϱB(ϖ)G(0,h(0))(1+γϖΓ(ϖ+1)tϖ),t1<tt2.

    In this section, we focus on EUS for the piecewise model (1.2). To accomplish this, we can make use of Lemma 2.7. Furthermore, in order to provide further clarification, we can also write the aforementioned model in the following form; for which, one can be benefit from [39,40,41].

    The solution of

    PFC0Dϖtg(t)=G(t,g(t)),  0<tT, (3.1)

    is

    g(t)={g0+t0g(s)ds,0<tt1,g0+1ϖB(ϖ)G(t,g(t))+ϖΓ(ϖ)B(ϖ)t0(ts)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ),t1<tt2,

    where

    g(t)={S1(t)S2(t)S3(t)E(t)I(t)R(t),   g0={S1(0)S2(0)S3(0)E(0)I(0)R(0),   gt1={S1(t1)S2t1S3(t1)Et1It1R(t1),   G(t,g(t))={G1(S1,S2,S3,E,I,R),G2(S1,S2,S3,E,I,R),G3(S1,S2,S3,E,I,R),G4(S1,S2,S3,E,I,R),G5(S1,S2,S3,E,I,R),G6(S1,S2,S3,E,I,R). (3.2)

    Let us take >t2t>t1>0 and the set E=C6[0,T] as a Banach space with

    g=maxt[0,T]maxi=1,...6|gi|,

    for g(t)=(g1,g2,,g6).

    The presumptions are as follows:

    (C1) Lg>0; G, ˉgE we arrive at

    |G(t,g)G(t,ˉg)|LG|gˉg|.

    (C2) CG>0 & MG>0,;

    |G(t,g(t))|CG|g|+MG.

    Theorem 3.1. With the assumption of (C2) and the piecewise function G, the TB model (3.1) has a solution.

    Proof. Considering a closed subset B of E,

    B={gE: gR1,2, R>0}.

    Assume that T:BB and from (5.1), we have the following

    T(g)={g0+t0g(s)ds,0<tt1,g0+1ϖB(ϖ)G(t,g(t))+ϖΓ(ϖ)B(ϖ)t0(ts)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ),t1<tt2.

    For the gB, we have

    |T(g)(t)|{|g0|+t10|G(ˉδ,g(ˉδ))|dˉδ,|gt1|+1ϖB(ϖ)|G(t,g(t))|+ϖΓ(ϖ)B(ϖ)t0(ts)ϖ1|G(s,g(s))|ds+1ϱB(ϖ)|G(0,g(0))(1+γϖΓ(ϖ+1)tϖ),{|g0|+t1(CG|g|+MG),|gt1|+1ϖB(ϖ)(CG|g|+MG)+1Γ(ϖ)B(ϖ)tϖ2(CG|g|+MG)|+1ϱB(ϖ)|(CG|g|+MG)(1+γϖΓ(ϖ+1)tϖ2),{R1, 0<tt1,R2, t1<tt2.

    Hence, T(B)B. Thus, we have that T is a closed and complete. Now, to show that the operator T is completely continuous, we assume that tn>tm[0,t1], which gives us the following:

    |T(g)(tn)T(g)(tm)|=|tn0G(s,g(s))dstm0G(s,g(s))ds|tn0|G(s,g(s))|dstm0|G(s,g(s))|ds[tn0(CG|g|+MG)tm0(CG|g|+MG)(CGg+MG)[tntm]. (3.3)

    Thanks to (3.3), when tmtn, we have that

    |T(g)(tn)T(g)(tm)|0, as tmtn.

    This implies equicontinuity of T in [0,t1]. For ti,tj[t1,T], we have

    |T(g)(tn)T(g)(tm)|=|g0+1ϖB(ϖ)G(tn,g(tn))+ϖΓ(ϖ)B(ϖ)tn0(tns)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖn)[g0+1ϖB(ϖ)G(tm,g(tm))+ϖΓ(ϖ)B(ϖ)tm0(tms)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖm)]1ϖB(ϖ)|G(tn,g(tn))G(tm,g(tm))|+ϖΓ(ϖ)B(ϖ)tn0[(tns)ϖ1(tms)ϖ1]|G(s,g(s))|ds+1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)|tϖntϖm|)1ϖB(ϖ)|G(tn,g(tn))G(tm,g(tm))|+|tϖntϖm|Γ(ϖ)B(ϖ)(CGg+MG)+1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)|tϖntϖm|). (3.4)

    If tntm, then

    |T(g)(tn)T(g)(tm)|0, as tmtn.

    This demonstrates that T is equally continuous over the interval [t1,t2]. As a result, we proved the equicontinuity of the operator T. Thus, by the Arzel'a-Ascoli and Schauder theorems TB model (3.1) has a solution.

    Theorem 3.2. Let us assume that (C1) holds true and Lg(1ϖB(ϖ)+1B(ϖ)(t2t1)ϖΓ(ϖ))<1; then, the TB system with the piecewise operator given by (3.1) has a unique solution.

    Proof. We already discussed the continuity of T:BB. Furthermore, for g,ˉgB on [0,t1], we have

    T(g)T(ˉg)=maxt[0,t1]|t10G(s,g(s))dst10G(s,ˉg(s))ds|t1LGgˉg. (3.5)

    From (3.5), we have

    T(g)T(ˉg)t1LGgˉg. (3.6)

    This implies that T is a contraction and by Banach's contraction principle, (3.1) has a unique solution. Moreover, for t[t1,t2], we have

    T(g)T(ˉg)=maxt[t1,t2]|g0+1ϖB(ϖ)G(tn,g(tn))+ϖΓ(ϖ)B(ϖ)t2t1(t2s)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ)[ˉg0+1ϖB(ϖ)G(t,ˉg(t))+ϖΓ(ϖ)B(ϖ)t2t1(t2s)ϖ1G(s,ˉg(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ)]|1ϖB(ϖ)Lg|gˉg|+1B(ϖ)Lg|gˉg|(t2t1)ϖΓ(ϖ)=Lg(1ϖB(ϖ)+1B(ϖ)(t2t1)ϖΓ(ϖ))|gˉg|. (3.7)

    From (3.7), we have

    T(g)T(ˉg)Lg(1ϖB(ϖ)+1B(ϖ)(t2t1)ϖΓ(ϖ))|gˉg|. (3.8)

    So, T is a contraction. By the help of Banach's theorem, the model has a unique solution.

    Here, we give definition for the Hyers-Ulam (HU) stability and derive it for the proposed model.

    Definition 4.1. The piecewise model (1.2) is HU-stable if for α>0, and

    |PCCDϖtg(t)F(t,g(t),PCCDϖtg(t))|<α,forall,tB, (4.1)

    there is ¯gZ with H>0, a constant, satisfying that

    ||g¯g||ZHα,forall,tB. (4.2)

    Also, for a nondecreasing function Φ:[0,)R+, if we have

    ||g¯g||ZHΦ(α),atevery,tB, (4.3)

    where Φ(0)=0, then the model (1.2) is generalized HU-stable.

    Lemma 4.2. Consider the function

    PCC0Dϖtg(t)=F(t,g(t)),  0<ϖ1. (4.4)

    The solution of (4.4) is

    g(t)={g0+Φ(g)+t10F(s,g(s))ds+t10ϕ(s)ds, 0<tt1,g0+1ϖB(ϖ)G(t,g(t))+ϖΓ(ϖ)B(ϖ)t0(ts)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ),t1<tt2

    It is easy to obtain the following

    ||F(g)F(¯g)||{(t11Lf1Mf)α,tB1,[Lgα(1ϱB(ϖ)+1Γ(ϖ)B(ϖ))]α=Λα,tB2. (4.5)

    Theorem 4.3. With the help of lemma 4.2 and Lf1Mf<1, the solution of (1.2) is HU-stable as well as gHU- stable.

    Proof. Let us assume that g is a solution of the model (1.2) and consider an approaximate solution ¯g of (1.2), for tB, implies that

    ||g¯g||=suptB|g(g+Φ(¯g)+t10F(s,¯g(s))dˉs+t10g(s)ds)|suptB|g(g+Φ(¯g)+t10F(s,¯g(s))ds)|+suptB|Φ(g)Φ(¯g)+t10F(s,g(s))dˉs+t10F(s,¯g(s))ds|t1α+Lf1Mf||g¯g||Z.

    Upon further simplification, we get

    ||g¯g||(t11Lf1Mf)α. (4.6)

    Case 2: for tB, we have

    ||g¯g||=suptB|g0+1ϖB(ϖ)G(t,g(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G(s,g(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ)[g0+1ϖB(ϖ)G(t,ˉg(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G(s,ˉg(s))ds1ϱB(ϖ)G(0,g(0))(1+γϖΓ(ϖ+1)tϖ)]|suptB|1ϱB(ϖ)Lggˉg|+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G(s,g(s)G(s,ˉg(s)dsLg(1ϱB(ϖ)+1Γ(ϖ)B(ϖ))gˉg=ˉδggˉgα, (4.7)

    for ˉδg=Lgα(1ϱB(ϖ)+1Γ(ϖ)B(ϖ)). Since, g,ˉg are solutions of (1.2), this implies that these are the fixed points (FPs) of T, or mathematically, T(g)=g, and T(ˉg)=ˉg. Thus, from (4.7), we have

    T(g)T(ˉg)ˉδggˉg. (4.8)

    Also, we have

    |g¯g|=gT(g)+T(ˉg)¯ggT(g)+T(g)T(¯g). (4.9)

    With the use of (4.7), (4.8), and (4.10), we have

    g¯gα1ˉδg=Hα. (4.10)

    For H=11ˉδg. Thus, (1.2) is HU-stable. Also, with a replacement of α by Φ(α), (4.4), implies that

    ||g¯g||ZHΦ(α),ateachtB.

    From Φ(0)=0, we have that (1.2) is gHU-stable.

    Here, we aim to produce a numerical Scheme for (1.2). The scheme has been applied to obtain computational results. To demonstrate this, consider the following:

    S1(t)={S1(0)+t10G1(s,S1)ds,  0<tt1,S1(t1)+1ϖB(ϖ)G1(t,S1(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G1(s,S1(s))ds1ϱB(ϖ)G1(0,S1(0))(1+γϖΓ(ϖ+1)tϖ),  t1<tt2,,
    S2(t)={S2(0)+t10G2(s,S2)ds,  0<tt1,S2(t1)+1ϖB(ϖ)G2(t,S2(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G2(s,S2(s))ds1ϱB(ϖ)G2(0,S2(0))(1+γϖΓ(ϖ+1)tϖ),  t1<tt2,,S3(t)={S3(0)+t10G3(s,S3)ds,  0<tt1,S3(t1)+1ϖB(ϖ)G1(t,S3(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G3(s,S3(s))ds1ϱB(ϖ)G3(0,S3(0))(1+γϖΓ(ϖ+1)tϖ) t1<tt2,,
    E(t)={E0+t10G4(s,E)ds,  0<tt1,E(t1)+1ϖB(ϖ)G4(t,E(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G4(s,E(s))ds1ϱB(ϖ)G4(0,E(0))(1+γϖΓ(ϖ+1)tϖ) t1<tt2,,I(t))={I0+t10G5(s,I)ds,  0<tt1,I(t1)+1ϖB(ϖ)G5(t,I(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G5(s,I(s))ds1ϱB(ϖ)G5(0,I(0))(1+γϖΓ(ϖ+1)tϖ),  t1<tt2,,R(t)={R0+t10G6(s,R)ds,  0<tt1,R(t1)+1ϖB(ϖ)G6(t,R(t))+ϖΓ(ϖ)B(ϖ)t20(t2s)ϖ1G6(s,R(s))ds1ϱB(ϖ)G6(0,R(0))(1+γϖΓ(ϖ+1)tϖ),  t1<tt2.,

    At t=tn+1, we deduce for (5.1) that

    s(tn+1))={S0+t10G1(U,ˉδ)dˉδ,0<tt1,s(t1)+1Γ(ϖ)tn+1t1(tˉδ)ϖ1G1(U,ˉδ)dˉδ,  t1<tt2. (5.1)

    Now, (5.1) can be expressed by using Lagrange's interpolation polynomial (LIP), we have

    S1(tn+1)={S1(0)+{iK=2[512G1(U2,tK2)ˉδt43G1(U1,tK1)ˉδt+G1(U,tK)],S1(t1)+{+1ϖB(ϖ)G1(tk,S1(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G1(tk,S1(tk))((n+1k)ϖ(nk+ϖ+2)(nk)ϖ(n+2k+2ϖ))G1(tk1,S1k1)((n+1k)ϖ+1(nk)ϖ(n+1k+ϖ))]1ϖB(ϖ)G1(0,S1(0))(1+γϖΓ(ϖ+1)(kh)ϖ),}, (5.2)
    S2(tn+1)={S2(0)+{iK=2[512G2(U2,tK2)ˉδt43G2(U1,tK1)ˉδt+G2(U,tK)],S2(t1)+{+1ϖB(ϖ)G2(tk,S2(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G1(tk,S1(tk))((n+1k)ϖ(n+2k+ϖ)(nk)ϖ(n+2k+2ϖ))G2(tk1,S2k1)((n+1k)ϖ+1(nk)ϖ(n+1k+ϖ))]1ϖB(ϖ)G2(0,S2(0))(1+γϖΓ(ϖ+1)(kh)ϖ),}, (5.3)
    S3(tn+1)={S3(0)+{iK=2[512G3(U2,tK2)ˉδt43G3(U1,tK1)ˉδt+G3(U,tK)],S3(t1)+{+1ϖB(ϖ)G3(tk,S3(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G3(tk,S3(tk))((n+1k)ϖ(n+2k+ϖ)(nk)ϖ(n+2k+2ϖ))G3(tk1,S3k1)((n+1k)ϖ+1(nk)ϖ(n+1+ϖk))]1ϖB(ϖ)G3(0,S3(0))(1+γϖΓ(ϖ+1)(kh)ϖ)}, (5.4)
    E(tn+1)={E0+{iK=2[512G4(U2,tK2)ˉδt43G4(U1,tK1)ˉδt+G4(U,tK)],E(t1)+{+1ϖB(ϖ)G4(tk,E(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G4(tk,E(tk))((n+1k)ϖ(n+2k+ϖ)(nk)ϖ(n+2k+2ϖ))G4(tk1,Ek1)((n+1k)ϖ+1(nk)ϖ(n+1+ϖk))]1ϖB(ϖ)G4(0,E(0))(1+γϖΓ(ϖ+1)(kh)ϖ)}, (5.5)
    I(tn+1)={I0+{iK=2[512G5(U2,tK2)ˉδt43G5(U1,tK1)ˉδt+G5(U,tK)],I(t1)+{+1ϖB(ϖ)G5(tk,I(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G5(tk,I(tk))((n+1k)ϖ(n+2k+ϖ)(nk)ϖ(n+2k+2ϖ))G5(tk1,Ik1)((n+1k)ϖ+1(n+1k+ϖ)(nk)ϖ)]1ϖB(ϖ)G5(0,I(0))(1+γϖΓ(ϖ+1)(kh)ϖ),}, (5.6)
    R(tn+1)={R(0)+{iK=2[512G6(U2,tK2)ˉδt43G6(U1,tK1)ˉδt+G6(U,tK)],R(t1)+{+1ϖB(ϖ)G6(tk,R(tk))+ϱ1hϖΓ(ϖ+2)Σnk=1[G6(tk,R(tk))((n+1k)ϖ(n+2k+ϖ)(nk)ϖ(n+2k+2ϖ))G6(tk1,Rk1)((n+1k)ϖ+1(n+1k+ϖ)(nk)ϖ)]1ϖB(ϖ)G6(0,R(0))(1+γϖΓ(ϖ+1)(kh)ϖ).}. (5.7)

    Mycobacterium TB affects the brain, kidneys, and spine, but mostly the lungs. TB is one of the biggest infectious disease killers globally, and it spreads through coughs and sneezes. TB symptoms might be moderate and may not present until the disease progresses. Coughing, chest discomfort, fever, and weariness are common. Night sweats, weight loss, and appetite loss may occur. Untreated TB can cause lung damage, organ failure, and death [42]. A chest X-ray, TST, and sputum cultures can diagnose TB. TB treatment requires months of medication. Even if symptoms improve, patients must finish the medications to prevent the development of antibiotic-resistant germs. Early diagnosis, treatment, immunization, and infection control measures like hospital isolation and masks can prevent TB from spreading. TB remains a public health issue, especially in low- and middle-income nations, despite advances. The COVID-19 pandemic has disrupted TB detection and treatment, increasing infections in some nations. TB impacts millions of individuals globally. Early diagnosis, treatment, and prevention are crucial to minimizing the disease's burden, and more research is needed to find better medicines and vaccines to fight this global health issue.

    Here, we present the application of the computational scheme for the (1.2) with the parametric values given by A=1.623107,μ=0.0025,p=0.05,γ=0.496,η=0.00341,ϕ=0.9,d=0.007,d1=0.0017,d2=0.0023,,ν=6,ϑ2=5.111010,ϑ1=1.2681010,ϑ3=2.553109,d3=0.0367,m1=0.079,m2=0.0067, and ,ξ=0.51. For the fractional orders φ=1.0,φ=0.99,φ=0.98,φ=0.97 and ,φ=0.96, with the initial conditions given by S1(0)=26504104,S2(0)=94197104,S3(0)=10055104,E(0)=9702790.121,I(0)=1259308, and R(0)=776223 as given in [13]. Figure 1 shows the dynamics of the susceptible group S1(t).

    Figure 1.  The dynamics of the youngest group of susceptible individuals in the piecewise TB model (1.2) with t1=35.

    Here we describe the simulations that are shown in the figures. Thus, the differential operator and integral operator that we propose have been derived from a combination of classical and mABC operators. As such, we have assumed a study period of 100 days, that is, I=[0,100]. This interval was divided into two sub-intervals. The first one is I1=[0,35] and the second one is [35,100]. Thus, t1=35 and t2=100. Figure 1 presents the results of simulations for the S1(t) dynamics for the crossover behaviour is presented in which is for the younger most susceptible population, which exhibited a decrease. The role of the modified derivative in the graph is apparent. Because it is closer to the classical case that is 1.

    There is a slight higher rate of decrease. And beyond t1, the cross over behaviour is well established. We have presumed the orders of the derivative to be 1.0,0.99,0.98,0.97,0.96. This class is a bit sensitive to the cross over behaviour as compared to the results in Figure 2, and Figure 3. Both of these graphs show a sudden decrease of the population of the susceptible classes. This decrease in population indicates a transition of individuals to the exposed group, the results of which are given in the Figure 4. In this class, the role of the derivative is very much apparent. In the early days there is a rise in the dynamics of the class and then we have a slight and gradual decrease in the population. This decrease is clearly reflected in the change in population of the infected group, as shown in Figure 5. There is an increase in the population of the recovered class in the [t1,t2], as shown in Figure 6. Here, as the order decreases, there is an increase in the population of the group. The cross over behavior of the recovered class can be seen in Figure 6.

    Figure 2.  The dynamics of the middle-aged susceptible individuals for the piecewise TB model (1.2) with t1=35.
    Figure 3.  The dynamics of the elderly susceptible individuals for the piecewise TB model (1.2) with t1=35.
    Figure 4.  The dynamics of exposed individuals for the piecewise TB model (1.2) with t1=35.
    Figure 5.  The dynamics of infected individuals for the piecewise TB model (1.2) with t1=35.
    Figure 6.  The recovered individuals for the piecewise TB model (1.2) with t1=35.

    In this paper, we introduced a piecewise differential operator which is a combination of the classical derivative and the mABC fractional derivative. Its corresponding integral is also given. The new operator was applied to a TB model for the study of its dynamics. We have observed very nice interpretations of the dynamics in the two sub intervals. For the numerical simulations, we considered the the interval of study as [0,100] in days. This interval was divided into sub-intervals as I1=[0,t1) and [t1,t2], where t1=25 and t2=100. In the I1 interval, we have considered the classical case while in the second interval I2, we have the mABC operators. This work is an interpretation of the dynamics of crossover behavior of TB under the two types of operators.

    The EUS with a piecewise derivative was investigated for the TB-infection model. HU-stability based on the nonlinear analysis was used to show the stability of the solutions. To approximate the solution to the stated issue, the piecewise LIP was used. The computations for the proposed infection model have been given for various fractional orders. In the simulation aspect, we have observed that beyond t1, the cross over behaviour is well established. We have presumed the orders of the derivative to be 1.0,0.99,0.98,0.97,0.96. Figures 2 and 3 show sudden decreases in the population of the susceptible classes. These decreases in population are included as a component of the population that is exposed, and its details are presented in Figure 4. Within this context, the function of the derivative is made abundantly clear. After an initial period in which there is an increase in the population of the exposed group, there is a minor and steady drop in the overall population of the class. This drop is likely the one that was reflected in the infected group shown in Figure 5. There is also increase in the number of people belonging to the recovered class during the time period [t1,t2]. In this situation, the population of the group is growing at a rate that is disproportionately faster than the order's decline.

    The authors declare that they have not used Artificial Intelligence tools in the creation of this article.

    Hasib Khan and Jehad Alzabut express their sincere thanks to Prince Sultan University. Jehad Alzabut, also expresses his gratitude to OSTİM Technical University for their endless support. The fifth author would like to thank Azarbaijan Shahid Madani University.

    The authors declare no conflict of interest.

    Hasib Khan: Conceptualization, Methodology, Writing- Original draft preparation, Software; Jehad Alzabut: Conceptualization, Methodology, Software, Writing-Original draft preparation; J.F. Gómez Aguilar: Conceptualization, Methodology, Writing- Original draft preparation, Software, Supervision; Praveen Agarwal: Methodology, Software, Writing-Original draft preparation. All authors read and approved the final manuscript.



    [1] R. Rosen, Dynamical system theory in biology, NY: Wiley and Sons Inc., 1970.
    [2] P. Waltman, Deterministic threshold models in the theory of epidemics, New York: Springer-Verlag, 1974.
    [3] D. Šiljak, Large scale dynamic systems: stability and structure, New York: Noth-Holland, 1978.
    [4] A. Coddington, N. Levinson, Theory of ordinary differential equations, New York: McGraw Hill, 1955.
    [5] P. Hartman, Ordinary differential equations, New York: John Wiley, 1964.
    [6] V. Lakshmikantham, S. Leela, Differential equations, Academic Press, I– II (1969).
    [7] G. Ladde, S. Ladde, Dynamic processes under random environment, Bull. Marathwada Math. Soc., 8 (2007), 96–123.
    [8] N. Katugampola, A new approach to generalized fractional derivatives, Bull. Math. Anal. Appl., 6 (2014), 1–15. https://doi.org/10.48550/arXiv.1106.0965 doi: 10.48550/arXiv.1106.0965
    [9] A. Mohamed, H. Faeza, T. Nidal, Computational technique to study analytical solutions to the fractional modified KDV-Zakharov-Kuznetsov equation, Abstr. Appl. Anal., 2022 (2022), 2162356. https://doi.org/10.1155/2022/2162356 doi: 10.1155/2022/2162356
    [10] H. Faeza, A. Mohamed, The generalized (2+1) and (3+1)- dimensional with advanced analytical wave solutions via computational applications, Int. J. Nonlinear Anal., 12 (2021), 1213–1241.
    [11] A. Mohamed, H. Faeza, Advantages of the differential equations for solving problems in mathematical physics with symbolic computation, Math. Model. Eng. Probl., 9 (2022), 268–276. https://doi.org/10.18280/mmep.090133 doi: 10.18280/mmep.090133
    [12] S. Alshammari, W. W. Mohammed, S. K. Samura, S. Faleh, The analytical solutions for the stochastic-fractional Broer-Kaup equations, Math. Probl. Eng., 2022 (2022). https://doi.org/10.1155/2022/6895875 doi: 10.1155/2022/6895875
    [13] M. Alquran, M. Khaled, H. Ananbeh, New soliton solutions for systems of nonlinear evolution equations by the rational sine-cosine method, Stud. Math. Sci., 3 (2011), 1–9. https://doi.org/10.3968/j.sms.1923845220110301.105 doi: 10.3968/j.sms.1923845220110301.105
    [14] M. A. Abdoon, R. Saadeh, M. Berir, F. E. Guma, Analysis, modeling and simulation of a fractional-order influenza model, Alex. Eng. J., 74 (2023), 231–240. https://doi.org/10.1016/j.aej.2023.05.011
    [15] M. S. Osman, On multi-soliton solutions for the (2+1)-dimensional breaking soliton equation with variable coefficients in a graded-index waveguide, Comput. Math. Appl., 75 (2018), 1–6. https://doi.org/10.1016/j.camwa.2017.08.033 doi: 10.1016/j.camwa.2017.08.033
    [16] C. M. Szpilka, R. L. Kolar, Numerical analogs to Fourier and dispersion analysis: Development, verification and application to the shallow water equations, Adv. Water Resour., 26 (2003), 649–662. https://doi.org/10.1016/S0309-1708(03)00028-9
    [17] H. Yasmin, A. S. Alshehry, A. H. Ganie, A. M. Mahnashi, R. Shah, Perturbed Gerdjikov-Ivanov equation: Soliton solutions via Backlund transformation, Optik, 298 (2024), 171576. https://doi.org/10.1016/j.ijleo.2023.171576 doi: 10.1016/j.ijleo.2023.171576
    [18] A. H. Ganie, F. Aljuaydi, Z. Ahmad, E. Bonyah, N. Khan, N. S. Alharthi, et al., A fractal-fractional perspective on chaotic behavior in 4D memristor-nonlinear system, AIP Adv., 14 (2024). https://doi.org/10.1063/5.0187218
    [19] L. Yang, M. Ur Rahman, M. A. Khan, Complex dynamics, sensitivity analysis and soliton solutions in the (2+1)-dimensional nonlinear Zoomeron model, Results Phys., 56 (2024), 107261. https://doi.org/10.1016/j.rinp.2023.107261
    [20] U. Ali, A. H. Ganie, I. Khan, F. Alotaibi, K. Kamran, S. Muhammad, et al., Traveling wave solutions to a mathematical model of fractional order (2+1)-dimensional breaking soliton equation, Fractals, 30 (2022), 2240124. https://doi.org/10.1142/S0218348X22401247
    [21] A. Qazza, R. Hatamleh, The existence of a solution for semi-linear abstract differential equations with infinite B-chains of the characteristic sheaf, Int. J. Appl. Math., 31 (2018), 611. https://doi.org/10.12732/ijam.v31i5.7 doi: 10.12732/ijam.v31i5.7
    [22] W. W. Mohammed, M. El-Morshedy, A. Moumen, E. E. Ali, M. Benaissa, A. E. Abouelregal, Effects of M-truncated derivative and multiplicative noise on the exact solutions of the breaking soliton equation, Symmetry, 15 (2023), 288. https://doi.org/10.3390/sym15020288 doi: 10.3390/sym15020288
    [23] F. M. Al-Askar, W. W. Mohammed, A. M. Albalahi, M. El-Morshedy, The impact of the Wiener process on the analytical solutions of the stochastic (2+1)-dimensional breaking soliton equation by using tanh-coth method, Mathematics, 10 (2022), 817. https://doi.org/10.3390/math10050817
    [24] K. Hosseini, F. Alizadeh, K. Sadri, E. Hinçal, A. R. Z. U. Akbulut, H. M. Alshehri, et al., Lie vector fields, conservation laws, bifurcation analysis, and Jacobi elliptic solutions to the Zakharov-Kuznetsov modified equal-width equation, Opt. Quant. Electron., 56 (2024), 506. https://doi.org/10.1007/s11082-023-06086-9
    [25] Z. Feng, The first-integral method to study the Burgers-Korteweg-de Vries equation, J. Phys. A-Math. Gen., 35 (2002), 343. https://doi.org/10.1088/0305-4470/35/2/312 doi: 10.1088/0305-4470/35/2/312
    [26] A. M. Wazwaz, Integrable (2+1)-dimensional and (3+1)-dimensional breaking soliton equations, Phys. Scr., 81 (2010), 035005. https://doi.org/10.1088/0031-8949/81/03/035005 doi: 10.1088/0031-8949/81/03/035005
    [27] R. Saadeh, M. Abu-Ghuwaleh, A. Qazza, E. Kuffi, A fundamental criteria to establish general formulas of integrals, J. Appl. Math., 2022 (2022). https://doi.org/10.1155/2022/6049367 doi: 10.1155/2022/6049367
    [28] M. A. Akbar, N. H. M. Ali, S. T. Mohyud-Din, Assessment of the further improved (G'/G)-expansion method and the extended tanh-method in probing exact solutions of nonlinear PDEs, SpringerPlus, 2 (2013), 1-9. https://doi.org/10.1186/2193-1801-2-1
    [29] E. Atilgan, M. Senol, A. Kurt, O. Tasbozan, New wave solutions of time-fractional coupled Boussinesq-Whitham-Broer-Kaup equation as a model of water waves, China Ocean Eng., 33 (2019), 477-483. https://doi.org/10.1007/s13344-019-0045-1 doi: 10.1007/s13344-019-0045-1
    [30] R. Saadeh, O. Ala'yed, A. Qazza, Analytical solution of coupled Hirota-Satsuma and KdV equations, Fractal Fract., 6 (2022), 694. https://doi.org/10.3390/fractalfract6120694 doi: 10.3390/fractalfract6120694
    [31] I. Podlubny, Fractional differential equations, Academic Press, 1999.
    [32] S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional integrals and derivatives: Theory and applications, Gordon and Breach Science Publishers, 1993.
    [33] M. M. Meerschaert, H. P. Scheffler, Limit distributions for sums of independent random vectors: Heavy tails in theory and practice, John Wiley & Sons, 2004.
    [34] M. Abu-Ghuwaleh, R. Saadeh, A. Qazza, General master theorems of integrals with applications, Mathematics, 10 (2022), 3547. https://doi.org/10.3390/math10193547 doi: 10.3390/math10193547
    [35] S. B. G. Karakoc, K. K. Ali, D. Y. Sucu, A new perspective for analytical and numerical soliton solutions of the Kaup-Kupershmidt and Ito equations, J. Comput. Appl. Math., 421 (2023), 114850. https://doi.org/10.1016/j.cam.2022.114850 doi: 10.1016/j.cam.2022.114850
    [36] M. A. Akbar, A. M. Wazwaz, F. Mahmud, D. Baleanu, R. Roy, H. K. Barman, et al., Dynamical behavior of solitons of the perturbed nonlinear Schrödinger equation and microtubules through the generalized Kudryashov scheme, Results Phys., 43 (2022), 106079. https://doi.org/10.1016/j.rinp.2022.106079
    [37] F. S. Alshammari, M. Asif, M. F. Hoque, A. Aldurayhim, Bifurcation analysis on ion sound and Langmuir solitary waves solutions to the stochastic models with multiplicative noises, Heliyon, 9 (2023). https://doi.org/10.1016/j.heliyon.2023.e16570
    [38] Y. Q. Chen, Y. H. Tang, J. Manafian, H. Rezazadeh, M. S. Osman, Dark wave, rogue wave and perturbation solutions of Ivancevic option pricing model, Nonlinear Dynam., 105 (2021), 2539-2548. https://doi.org/10.1007/s11071-021-06642-6
    [39] M. M. Hossain, M. A. N. Sheikh, M. M. Roshid, M. A. Taher, New soliton solutions and modulation instability analysis of the regularized long-wave equation in the conformable sense, Part. Differ. Equ. Appl. Math., 9 (2024), 100615. https://doi.org/10.1016/j.padiff.2024.100615 doi: 10.1016/j.padiff.2024.100615
    [40] Z. Islam, M. A. N. Sheikh, H. O. Roshid, M. A. Hossain, M. A. Taher, A. Abdeljabbar, Stability and spin solitonic dynamics of the HFSC model: Effects of neighboring interactions and crystal field anisotropy parameters, Opt. Quant. Electron., 56 (2024), 190. https://doi.org/10.1007/s11082-023-05739-z doi: 10.1007/s11082-023-05739-z
    [41] S. Malik, M. S. Hashemi, S. Kumar, H. Rezazadeh, W. Mahmoud, M. S. Osman, Application of new Kudryashov method to various nonlinear partial differential equations, Opt. Quant. Electron., 55 (2023), 8. https://doi.org/10.1007/s11082-022-04261-y doi: 10.1007/s11082-022-04261-y
    [42] H. U. Rehman, G. S. Said, A. Amer, H. Ashraf, M. M. Tharwat, M. Abdel-Aty, et al., Unraveling the (4+1)-dimensional Davey-Stewartson-Kadomtsev-Petviashvili equation: Exploring soliton solutions via multiple techniques, Alex. Eng. J., 90 (2024), 17-23. https://doi.org/10.1016/j.aej.2024.01.058
    [43] M. S. Ullah, H. O. Roshid, M. Z. Ali, New wave behaviors and stability analysis for the (2+1)-dimensional Zoomeron model, Opt. Quant. Electron., 56 (2024), 240. https://doi.org/10.1007/s11082-023-05804-7 doi: 10.1007/s11082-023-05804-7
    [44] J. Liu, W. Wei, W. Xu, An averaging principle for stochastic fractional differential equations driven by fBm involving impulses, Fractal Fract., 6 (2022), 256. https://doi.org/10.3390/fractalfract6050256 doi: 10.3390/fractalfract6050256
    [45] E. Salah, R. Saadeh, A. Qazza, R. Hatamleh, Direct power series approach for solving nonlinear initial value problems, Axioms, 12 (2023), 111. https://doi.org/10.3390/axioms12020111 doi: 10.3390/axioms12020111
    [46] A. B. M. Alzahrani, R. Saadeh, M. A. Abdoon, M. Elbadri, M. Berir, A. Qazza, Effective methods for numerical analysis of the simplest chaotic circuit model with Atangana-Baleanu Caputo fractional derivative, J. Eng. Math., 144 (2024), 9. https://doi.org/10.1007/s10665-023-10319-x doi: 10.1007/s10665-023-10319-x
    [47] M. Elbadri, M. A. Abdoon, M. Berir, D. K. Almutairi, A numerical solution and comparative study of the symmetric rossler attractor with the generalized Caputo fractional derivative via two different methods, Mathematics, 11 (2023), 2997. https://doi.org/10.3390/math11132997 doi: 10.3390/math11132997
    [48] M. A. Abdoon, R. Saadeh, M. Berir, F. E. Guma, Analysis, modeling and simulation of a fractional-order influenza model, Alex. Eng. J., 74 (2023), 231-240. https://doi.org/10.1016/j.aej.2023.05.011
    [49] R. Saadeh, M. A. Abdoon, A. Qazza, M. Berir, A numerical solution of generalized Caputo fractional initial value problems, Fractal Fract., 7 (2023), 332. https://doi.org/10.3390/fractalfract7040332 doi: 10.3390/fractalfract7040332
    [50] M. Elbadri, M. A. Abdoon, M. Berir, D. K. Almutairi, A symmetry chaotic model with fractional derivative order via two different methods, Symmetry, 15 (2023), 1151. https://doi.org/10.3390/sym15061151 doi: 10.3390/sym15061151
    [51] A. B. M. Alzahrani, M. A. Abdoon, M. Elbadri, M. Berir, D. E. Elgezouli, A comparative numerical study of the symmetry chaotic jerk system with a hyperbolic sine function via two different methods, Symmetry, 15 (2023), 1991. https://doi.org/10.3390/sym15111991 doi: 10.3390/sym15111991
    [52] A. Qazza, M. Abdoon, R. Saadeh, M. Berir, A new scheme for solving a fractional differential equation and a chaotic system, Eur. J. Pure Appl. Math., 16 (2023), 1128-1139. https://doi.org/10.29020/nybg.ejpam.v16i2.4769 doi: 10.29020/nybg.ejpam.v16i2.4769
    [53] F. Hasan, M. A. Abdoon, R. Saadeh, M. Berir, A. Qazza, A new perspective on the stochastic fractional order materialized by the exact solutions of Allen-Cahn equation, Int. J. Math. Eng. Manag. Sci., 8 (2023), 912. http://dx.doi.org/10.33889/IJMEMS.2023.8.5.052 doi: 10.33889/IJMEMS.2023.8.5.052
    [54] M. M. Al-Sawalha, S. Mukhtar, R. Shah, A. H. Ganie, K. Moaddy, Solitary waves propagation analysis in nonlinear dynamical system of fractional coupled Boussinesq-Whitham-Broer-Kaup equation, Fractal Fract., 7 (2023), 889. https://doi.org/10.3390/fractalfract7120889 doi: 10.3390/fractalfract7120889
    [55] M. A. Iqbal, M. M. Miah, H. S. Ali, N. H. M. Shahen, A. Deifalla, New applications of the fractional derivative to extract abundant soliton solutions of the fractional order PDEs in mathematics physics, Part. Differ. Equ. Appl. Math., 9 (2024), 100597. https://doi.org/10.1016/j.padiff.2023.100597 doi: 10.1016/j.padiff.2023.100597
  • This article has been cited by:

    1. Kamal Shah, Muhammad Sarwar, Thabet Abdeljawad, On rotavirus infectious disease model using piecewise modified ABC fractional order derivative, 2024, 19, 1556-1801, 214, 10.3934/nhm.2024010
    2. Nana Pan, Akbar Zada, Ioan-Lucian Popa, Fairouz Tchier, On the (k,φ)-Hilfer Langevin fractional coupled system having multi point boundary conditions and fractional integrals, 2024, 20904479, 103111, 10.1016/j.asej.2024.103111
    3. Aziz Khan, Kamal Shah, Thabet Abdeljawad, Inas Amacha, Fractal fractional model for tuberculosis: existence and numerical solutions, 2024, 14, 2045-2322, 10.1038/s41598-024-62386-4
    4. Hasib Khan, Jehad Alzabut, D. K. Almutairi, Wafa Khalaf Alqurashi, The Use of Artificial Intelligence in Data Analysis with Error Recognitions in Liver Transplantation in HIV-AIDS Patients Using Modified ABC Fractional Order Operators, 2024, 9, 2504-3110, 16, 10.3390/fractalfract9010016
    5. Hasib Khan, Jehad Alzabut, Mohamed Tounsi, Dalal Khalid Almutairi, AI-Based Data Analysis of Contaminant Transportation with Regression of Oxygen and Nutrients Measurement, 2025, 9, 2504-3110, 125, 10.3390/fractalfract9020125
    6. Ayesha Saleem, Mati ur Rahman, Salah Boulaaras, Rafik Guefaifia, Dumitru Baleanu, Exploring the dynamics of HIV and HCV co-infection through piecewise modified Mittag-Leffler fractional derivatives, 2025, 33, 2769-0911, 10.1080/27690911.2025.2478038
    7. Mukhtiar Khan, Nadeem Khan, Ibad Ullah, Kamal Shah, Thabet Abdeljawad, Bahaaeldin Abdalla, A novel fractal fractional mathematical model for HIV/AIDS transmission stability and sensitivity with numerical analysis, 2025, 15, 2045-2322, 10.1038/s41598-025-93436-0
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(905) PDF downloads(56) Cited by(2)

Figures and Tables

Figures(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog