Loading [MathJax]/jax/output/SVG/jax.js
 

Mixed vaccination strategy for the control of tuberculosis: A case study in China

  • This study first presents a mathematical model of TB transmission considering BCG vaccination compartment to investigate the transmission dynamics nowadays. Based on data reported by the National Bureau of Statistics of China, the basic reproduction number is estimated approximately as R0=1.1892. To reach the new End TB goal raised by WHO in 2015, considering the health system in China, we design a mixed vaccination strategy. Theoretical analysis indicates that the infectious population asymptotically tends to zero with the new vaccination strategy which is the combination of constant vaccination and pulse vaccination. We obtain that the control of TB is quicker to achieve with the mixed vaccination. The new strategy can make the best of current constant vaccination, and the periodic routine health examination provides an operable environment for implementing pulse vaccination in China. Numerical simulations are provided to illustrate the theoretical results and help to design the final mixed vaccination strategy once the new vaccine comes out.

    Citation: Siyu Liu, Yong Li, Yingjie Bi, Qingdao Huang. Mixed vaccination strategy for the control of tuberculosis: A case study in China[J]. Mathematical Biosciences and Engineering, 2017, 14(3): 695-708. doi: 10.3934/mbe.2017039

    Related Papers:

    [1] Saima Rashid, Rehana Ashraf, Fahd Jarad . Strong interaction of Jafari decomposition method with nonlinear fractional-order partial differential equations arising in plasma via the singular and nonsingular kernels. AIMS Mathematics, 2022, 7(5): 7936-7963. doi: 10.3934/math.2022444
    [2] Mohammed A. Almalahi, Satish K. Panchal, Fahd Jarad, Mohammed S. Abdo, Kamal Shah, Thabet Abdeljawad . Qualitative analysis of a fuzzy Volterra-Fredholm integrodifferential equation with an Atangana-Baleanu fractional derivative. AIMS Mathematics, 2022, 7(9): 15994-16016. doi: 10.3934/math.2022876
    [3] Kishor D. Kucche, Sagar T. Sutar, Kottakkaran Sooppy Nisar . Analysis of nonlinear implicit fractional differential equations with the Atangana-Baleanu derivative via measure of non-compactness. AIMS Mathematics, 2024, 9(10): 27058-27079. doi: 10.3934/math.20241316
    [4] Saima Rashid, Fahd Jarad, Fatimah S. Bayones . On new computations of the fractional epidemic childhood disease model pertaining to the generalized fractional derivative with nonsingular kernel. AIMS Mathematics, 2022, 7(3): 4552-4573. doi: 10.3934/math.2022254
    [5] Muhammad Farman, Ali Akgül, Kottakkaran Sooppy Nisar, Dilshad Ahmad, Aqeel Ahmad, Sarfaraz Kamangar, C Ahamed Saleel . Epidemiological analysis of fractional order COVID-19 model with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(1): 756-783. doi: 10.3934/math.2022046
    [6] Ahu Ercan . Comparative analysis for fractional nonlinear Sturm-Liouville equations with singular and non-singular kernels. AIMS Mathematics, 2022, 7(7): 13325-13343. doi: 10.3934/math.2022736
    [7] Hasib Khan, Jehad Alzabut, J.F. Gómez-Aguilar, Praveen Agarwal . Piecewise mABC fractional derivative with an application. AIMS Mathematics, 2023, 8(10): 24345-24366. doi: 10.3934/math.20231241
    [8] Kottakkaran Sooppy Nisar, Aqeel Ahmad, Mustafa Inc, Muhammad Farman, Hadi Rezazadeh, Lanre Akinyemi, Muhammad Mannan Akram . Analysis of dengue transmission using fractional order scheme. AIMS Mathematics, 2022, 7(5): 8408-8429. doi: 10.3934/math.2022469
    [9] Veliappan Vijayaraj, Chokkalingam Ravichandran, Thongchai Botmart, Kottakkaran Sooppy Nisar, Kasthurisamy Jothimani . Existence and data dependence results for neutral fractional order integro-differential equations. AIMS Mathematics, 2023, 8(1): 1055-1071. doi: 10.3934/math.2023052
    [10] Eman A. A. Ziada, Salwa El-Morsy, Osama Moaaz, Sameh S. Askar, Ahmad M. Alshamrani, Monica Botros . Solution of the SIR epidemic model of arbitrary orders containing Caputo-Fabrizio, Atangana-Baleanu and Caputo derivatives. AIMS Mathematics, 2024, 9(7): 18324-18355. doi: 10.3934/math.2024894
  • This study first presents a mathematical model of TB transmission considering BCG vaccination compartment to investigate the transmission dynamics nowadays. Based on data reported by the National Bureau of Statistics of China, the basic reproduction number is estimated approximately as R0=1.1892. To reach the new End TB goal raised by WHO in 2015, considering the health system in China, we design a mixed vaccination strategy. Theoretical analysis indicates that the infectious population asymptotically tends to zero with the new vaccination strategy which is the combination of constant vaccination and pulse vaccination. We obtain that the control of TB is quicker to achieve with the mixed vaccination. The new strategy can make the best of current constant vaccination, and the periodic routine health examination provides an operable environment for implementing pulse vaccination in China. Numerical simulations are provided to illustrate the theoretical results and help to design the final mixed vaccination strategy once the new vaccine comes out.


    Korteweg and de-Vries developed the classical KdV equation in 1895 as a nonlinear PDE to investigate the waves that occurs on the surfaces of shallow water. Many studies have been conducted on this exactly solvable model. Many scholars have proposed novel applications of the classical KdV equation, such as acoustic waves that produces in a plasma in ion form, and acoustic waves which are produces on a crystal lattice. In [1], the classical KdV equation is as follows

    Ut+β1Uxxx+β2UUx=0. (1.1)

    Different variations of Eq (1.1) have been published in the literature, including [2,3]. In literature [2], various efficient methods have been used to get solitons solutions of different kinds of KdV equations. Wang [3] used a quadratic function Ansatz to obtain lump solutions for the (2+1)-dimensional KdV equation. In [4,5,6], several systematic techniques were utiltized to investigate the different types of KdV equations. Wang and Kara introduced a new 2D-mKdV in 2019 [7]. The new (2+1)-dimensional mKdV equation is given by

    ft=6f2fx6f2fy+fxxxfyyy3fxxy+3fxyy. (1.2)

    After 1695, the non-integer order or fractional-order derivative (FOD) was described as a simple academic generalization of the classical derivative. A FOD is an operator that extends the order of differentiation from Natural numbers (N) to a set of real numbers (R) or even to a set of complex numbers (C). Fractional calculus has emerged as one of the most effective methods for describing long-memory processes over the last decade. Engineers and physicists, as well as pure mathematicians, are interested in such models. The models that are represented by differential equations with fractional-order derivatives are the most interesting [8,9,10]. Their evolutions are much more complicated than in the classical integer-order case, and deciphering the underlying principle is a difficult task. There are many fractional operators regarding to the kernels involved in the integration. The most popular fractional operator is Caputo-Liouville which is based on power-law kernel, but this kernel has issue about the singularity of the kernel. To tackle the limitation of derivative operators with power-law kernels, new types of nonlocal FOD have been introduced in recent literature. For example, in literature [11], Caputo and Fabrizio (CF) introduced a FOD that is focused on the exponential kernel. However, the CF derivative, on the other hand, has some issues with the kernel's locality. In 2016, Atangana and Baleanu constructed an updated version of FOD that is based on the Mittag-Leffler function [12]. This derivative solves the issues of locality and singularity. The FOD of Atangana-Baleanu in Caputo sense (ABC) accurately describes the memory. The ABC operator's most significant applications are available in [13,14,15,16,17].

    Certain important techniques have been used to solve fractional order differential equations (FODEs). Some of these includes the homotopy perturbation method (HPM), Laplace Adomian decomposition method (LADM) fractional operational matrix method (FOMM), homotopy analysis method (HAM) and many more [18,19,20,21,22]. In contrast to these techniques, the Laplace-Adomian decomposition method (LADM) is an important tool for solving non linear FODEs. The ADM and the Laplace transformation are two essential methods that are combined in LADM. Furthermore, unlike a Runge-Kutta process, LADM does not require a predefine size declaration. Every numerical or analytical approach has its own set of benefits and drawbacks. For instance, discretization of data is used in collocation techniques that require extra memory and a longer operation. Since the Laplace-Adomian approach has less parameters than all other methods, it is a useful tool that does not necessitate discretion or linearization [23]. With the aid of LADM, the smoke model was successfully solved in [24]. LADM was utilized by the authors to solve third order dispersive PDE defined by FOD in [22].

    Inspired by above literature, in this paper, we study Eq (1.2) under ABC-operator. We use LADM to solve the proposed equation. Consider Eq (1.2) under ABC-operator as

    ABCDΨtft=6f2fx6f2fy+fxxxfyyy3fxxy+3fxyy. (1.3)

    Definition 2.1. [12] Let 0<Ψ1 and f(t)H1. Then ABC FOD of order Ψ is expressed as

    ABCDΨtf(t)=c(Ψ)(1Ψ)t0f(η)EΨ((tη)ΨΨ(1Ψ))dη,

    where c(Ψ) is the normalization function such that c(0)=c(1)=1. The symbol EΨ denotes the Mittag-Leffler kernel which is defined as:

    EΨ(t)=k=0tkΓ(Ψk+1).

    Definition 2.2. [12] Let 0<Ψ1 and f(t)H1(0,T). Then AB fractional integral of order Ψ is defined as

    ABIΨtf(t)=(1Ψ)c(Ψ)f(t)+(Ψ)c(Ψ)Γ(Ψ)t0(tη)Ψ1f(η)dη.

    Definition 2.3. [12] The formula for the Laplace transform of ABC FOD of f(t) is defined by

    L[ABCDΨ0f(t)]=c(Ψ)sΨ(1Ψ)+Ψ[sΨL[f(t)]sΨ1u(0)].

    Theorem 2.4. [25] Let H be a Banach space and X:HHbe a mapping. Then X is said to be Picard's X-stable, if ξ,mH,

    XξXmaξXξ+bξm,

    where a0, and b[0,1]. Further, X has a fixedpoint.

    Theorem 3.1. [12] The ABC derivative holds the following Lipschitztype condition for 0<ϝ<.

    ABCDΨtf(t)ABCDΨtg(t)ϝf(t)g(t).

    Proof. By using definition of ABC, wehave

    ABCDΨtf(t)ABCDΨtg(t)=c(Ψ)(1Ψ)t0f(η)EΨ((tη)ΨΨ(1Ψ))dηc(Ψ)(1Ψ)t0g(η)EΨ((tη)ΨΨ(1Ψ))dη=c(Ψ)(1Ψ)[t0f(η)EΨ((tη)ΨΨ(1Ψ))dηt0g(η)EΨ((tη)ΨΨ(1Ψ))dη].

    Using the Lipschitz condition for the first order derivative, we can find a small positive constant ρ1 such that

    ABCDΨtf(t)ABCDΨtg(t)=c(Ψ)ρ1(1Ψ)EΨ(ΨtΨ(1Ψ))[t0f(η)dηt0g(η)dη]c(Ψ)ρ1(1Ψ)EΨ(ΨtΨ(1Ψ))t0f(η)dηt0g(η)dηϝt0f(η)dηt0g(η)dηϝf(t)g(t),

    where ϝ=c(Ψ)ρ1(1Ψ)EΨ(ΨtΨ(1Ψ)). Thus Lipschitz condition holds for ABC derivative.

    Let

    Φ(x,y,t;f)=6f2fx6f2fy+fxxxfyyy3fxxy+3fxyy, (4.1)

    Eq (1.3) can be written as

    ABCDΨtf(x,y,t)=Φ(x,y,t;f). (4.2)

    Applying the ABC FOI to Eq (4.2), we have

    f(x,y,t)f(x,y,0)=(1Ψ)c(Ψ)Φ(x,y,t;f)+(Ψ)c(Ψ)Γ(Ψ)t0(tη)Ψ1Φ(x,y,t;f)dη.

    First we have to verify that the Lipschitz condition holds for the kernel Φ(x,y,t;f). For this, let us take two bounded functions, f and g, i.e., fΔ1, and gΔ2 where Δ1,Δ2>0, and consider

    Φ(x,y,t;f)Φ(x,y,t;g)={(6f2fx6g2gx)(6f2fy6g2gy)+(fxxxgxxx)(fyyygyyy)(3fxxy3gxxy)+(3fxyy3gxyy)={2x(f3g3)2y(f3g3)+3x3(fg)3y3(fg)33yx2(fg)+33y2x(fg){2x(f3g3)+2y(f3g3)+3x3(fg)+3y3(fg)+33yx2(fg)+33y2x(fg).

    Since f and g are bounded functions. Therefore, their partial derivatives satisfy the Lipschitz conditions and there exists non-negative constants K,L,M,N,O,P such that

    Φ(x,y,t;f)Φ(x,y,t;g){2K(f3g3)+2L(f3g3)+M(fg)+N(fg)+3O(fg)+3P(fg){(2K+2L)(f2+fg+g2)(fg)+M(fg)+N(fg)+3O(fg)+3P(fg)=((2K+2L)(Δ21+Δ1Δ2+Δ22)+M+N+3O+3P)(fg)=κ(fg).

    Let

    κ=((2K+2L)(Δ21+Δ1Δ2+Δ22)+M+N+3O+3P),

    thus

    Φ(x,y,t;f)Φ(x,y,t;g)κfg.

    For further analysis, we make an iterative scheme as

    fξ+1(x,y,t)=(1Ψ)c(Ψ)Φ(x,y,t;fξ)+(Ψ)c(Ψ)Γ(Ψ)t0(tη)Ψ1Φ(x,y,η;fξ)dη,

    where f0(x,y,t)=f(x,y,0). Now the difference between two consecutive terms can be taken as

    eξ(x,y,t)=fξ(x,y,t)fξ1(x,y,t)=(1Ψ)c(Ψ)[Φ(x,y,t;fξ1)Φ(x,y,t;fξ2)]+(Ψ)c(Ψ)Γ(Ψ)t0(tη)Ψ1[Φ(x,y,η;fξ1)Φ(x,y,η;fξ2)]dη.

    Also, we have

    fξ(x,y,t)=ξk=0ek(x,y,t), (4.3)

    with f1=0.

    Theorem 4.1. Assume that f(x,y,t) is bounded a function. Then

    eξ(x,y,t)((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))ξf(x,y,0). (4.4)

    Proof. Consider

    eξ(x,y,t)=fξ(x,y,t)fξ1(x,y,t). (4.5)

    The Eq (4.5) gets the form under norm as

    eξ(x,y,t)=fξ(x,y,t)fξ1(x,y,t).

    To get the required result, we use the concept of mathematical induction. For ξ=1, one can get

    e1(x,y,t)=f1(x,y,t)f0(x,y,t)(1Ψ)c(Ψ)Φ(x,y,t;f0)Φ(x,y,t;f1)+Ψc(Ψ)Γ(Ψ)t0(tη)Ψ1Φ(x,y,η;f0)Φ(x,y,η;f1)dη(1Ψ)c(Ψ)κf0f1+Ψc(Ψ)Γ(Ψ)κt0(tη)Ψ1f0f1dη=(1Ψ)c(Ψ)κf(x,y,0)+Ψc(Ψ)Γ(Ψ)κf(x,y,0)t0(tη)Ψ1dη=(1Ψ)c(Ψ)κf(x,y,0)+Ψc(Ψ)Γ(Ψ)κf(x,y,0)t=((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))f(x,y,0).

    Now assume that the result is true for ξ1, i.e.,

    eξ1(x,y,t)((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))ξ1f(x,y,0). (4.6)

    Next, we have to show that

    eξ(x,y,t)((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))ξf(x,y,0). (4.7)

    To get the result (4.4), consider

    eξ(x,y,t)=fξ(x,y,t)fξ1(x,y,t)(1Ψ)c(Ψ)Φ(x,y,t;fξ1)Φ(x,y,t;fξ2)+Ψc(Ψ)Γ(Ψ)t0(tη)Ψ1Φ(x,y,η;fξ1)Φ(x,y,η;fξ2)dη(1Ψ)c(Ψ)κfξ1fξ2+Ψc(Ψ)Γ(Ψ)κt0(tη)Ψ1fξ1fξ2dη=(1Ψ)c(Ψ)κeξ1+κtΨ1c(Ψ)Γ(Ψ)eξ1=((1Ψ)c(Ψ)κ+κtΨ1c(Ψ)Γ(Ψ))eξ1=((1Ψ)c(Ψ)κ+κtΨ1c(Ψ)Γ(Ψ))((1Ψ)c(Ψ)κ+κtΨ1c(Ψ)Γ(Ψ))ξ1f(x,y,0)=((1Ψ)c(Ψ)κ+κtΨ1c(Ψ)Γ(Ψ))ξf(x,y,0).

    This ends the proof.

    Theorem 4.2. If the following relation holds at t=t00, where

    0((1Ψ)c(Ψ)κ+κtΨ10c(Ψ)Γ(Ψ))<1. (4.8)

    Then at least one solution of the new 2D KdV equation under the ABCfractional derivative exists.

    Proof. With the help of Eq (4.3), we have

    fξ(x,y,0)ξk=0ek(x,y,t)ξk=0(((1Ψ)c(Ψ)κ+κtΨ1c(Ψ)Γ(Ψ))kf(x,y,0)),

    for t=t0, we get

    fξ(x,y,0)f(x,y,0)ξk=0((1Ψ)c(Ψ)κ+κtΨ10c(Ψ)Γ(Ψ))k.

    From the above relation, we can say that

    limξfξ(x,y,0)f(x,y,0)limξξk=0((1Ψ)c(Ψ)κ+κtΨ10c(Ψ)Γ(Ψ))k.

    Since

    0((1Ψ)c(Ψ)κ+κtΨ10c(Ψ)Γ(Ψ))<1, (4.9)

    this implies that sequence fξ(x,y,t) is convergent and therefore the sequence is bounded for each ξ. Further, assume that

    Rξ(x,y,t)=f(x,y,t)fξ(x,y,t).

    Since fξ(x,y,t) is bounded. It follows that for λ>0, we have fξ(x,y,t)λ. After simple manipulation like we did in Theorems 4.1 and 4.2, we obtain

    Rξ(x,y,t)((1Ψ)c(Ψ)κ+κtΨ10c(Ψ)Γ(Ψ))ξ+1λ.

    Using Eq (4.9), one can get

    limξRξ(x,y,t)=0,

    it follows that limξfξ(x,y,t)=f(x,y,t). This finish the proof.

    Theorem 4.3. If the inequality (4.8) holds at t=t00.Then the unique solution of proposed equation exists.

    Proof. On contrary, suppose that there are two solutions f and g of the proposed equation such that fg. Now

    f(x,y,t)g(x,y,t)=(1Ψ)c(Ψ)[Φ(x,y,t;f)Φ(x,y,t;g)]+Ψc(Ψ)Γ(Ψ)×[t0(tη)Ψ1[Φ(x,y,η;f)Φ(x,y,η;g)]dη].

    Taking norm both side

    f(x,y,t)g(x,y,t)(1Ψ)c(Ψ)Φ(x,y,t;f)Φ(x,y,t;g)+Ψc(Ψ)Γ(Ψ)×[t0(tη)Ψ1Φ(x,y,η;f)Φ(x,y,η;g)dη](1Ψ)c(Ψ)κfg+Ψc(Ψ)Γ(Ψ)t0κ(tη)Ψ1fgdη((1Ψ)c(Ψ)κ+Ψc(Ψ)Γ(Ψ)κ)fgt0(tη)Ψ1dη=((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))fg,

    but

    0((1Ψ)c(Ψ)κ+κtΨc(Ψ)Γ(Ψ))<1.

    Using the above inequality, we achieve

    f(x,y,t)g(x,y,t)=0,

    thus, our supposition is wrong. Hence, the solution is unique.

    In this section, we briefly discuss the solution of the model by applying Ansatz method. For this purpose, we will consider a test function as

    f(x,y,t)=β0+β1sech(b1x+b2y+b3t). (5.1)

    By putting above equation into classical form of the model, we obtain

    6β20a16β20b2+b313b21b2+3b22b1b32b3=0,12β0β1b112β0b1b2=0,6β21b16β21b26b31+18b21b218b22b1+6b32=0,

    solution becomes as

    β0=0,β1=b1±b2,b3=b313b21b2+3b22b1b32,

    solution of classical model becomes as

    f1,2(x,y,t)=(b1±b2)sech(b1x+b2y+(b313b21b2+3b22b1b32)t).

    For b1=1 and b2=1, the above solution becomes

    f1,2(x,y,t)=4exp(xy+8t)1+exp(2x2y+16t). (5.2)

    In this section, to obtain analytic solution we applying Laplace transform (LT) on both sides of equations f(x,y,t) is the source term. Subject to the initial condition f(x,y,0)=f0(x,y,0). On utilizing LT, one can get

    L[ABCDΨtf]=L[6f2fx6f2fy+fxxxfyyy3fxxy+3fxyy],
    L[f(x,y,t)]=f(x,y,0)s+[s(1Ψ)+Ψc(Ψ)L[6f2fx6f2fy]+fxxxfyyy3fxxy+3fxyy]. (6.1)

    The approximate solution is represented by

    f(x,y,t)=ξ=0fξ(x,y,t), (6.2)

    and the nonlinear term is represented by Adomain polynomials, i.e., G(f)=f2=ξ=0Aξ, where Aξ is defined as follows for any ξ=0,1,2,

    Aξ=1Γ(ξ+1)dξdλξ[G(ξ=0(λξfξ))]λ=0.

    Using Eq (6.2), we obtain

    L[ξ=0fξ(x,y,t)]=f(x,y,0)s+s(1Ψ)+Ψc(Ψ)L[6ξ=0Aξξ=0xfξ6ξ=0Aξξ=0yfξ+ξ=0(3x3fξ3y3fξ33x2yfξ+33yx2fξ)].

    The following can be obtain by comparing terms

    L[f0(x,y,t)]=f(x,y,0)s,L[f1(x,y,t)]=sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6A0f0x6A0f0y+f0xxx+f0yyy3f0xxy+3f0xyy],L[f2(x,y,t)]=sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6A1f1x6A1f1y+f1xxx+f1yyy3f1xxy+3f1xyy],L[fξ+1(x,y,t)]=sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6Aξfξx6Aξfξy+fξxxx+fξyyy3fξxxy+3fξxyy].

    Applying L1, we get

    f0(x,y,t)=L1[f(x,y,0)s],f1(x,y,t)=L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6A0f0x6A0f0y+f0xxx+f0yyy3f0xxy+3f0xyy]]f2(x,y,t)=L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6A1f1x6A1f1y+f1xxx+f1yyy3f1xxy+3f1xyy]]fξ+1(x,y,t)=L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6Aξfξx6Aξfξy+fξxxx+fξyyy3fξxxy+3fξxyy]].

    The required series solution is given as

    f(x,y,t)=ξ=0fξ(x,y,t). (6.3)

    Here, we present two special cases of the proposed equation.

    For the first case, we take the initial condition as

    f(x,y,0)=4exp(x+y)1+exp(2(xy)).

    Using the detail procedure as discussed above, we achieve

    f0(x,y,t)=f(x,y,0)=4exp(x+y)1+exp(2(xy)),f1(x,y,t)=L1[sΨ(1Ψ)+ΨsΨc(Ψ)L[6f20xf06f20yf0+3x3f03y3f032x2(yf0)+3x(2y2f0)]],

    using Mathematica, we obtain

    f1(x,y,t)=(1Ψ+ΨtΨΓ(Ψ+1))(32texp(x+y)(exp(2x)+exp(2y))(exp(2x)+exp(2y))2c(Ψ)),

    similarly, other terms can be calculated with the help of Mathematica. The required series solution is given as

    {f(x,y,t)=4exp(x+y)1+exp(2(xy))+32t(1Ψ+ΨtΨΓ(Ψ+1))×(exp(x+y)(exp(2x)+exp(2y))(exp(2x)+exp(2y))2c(Ψ))+ (6.4)

    Remark 6.1. When we put Ψ=1 in Eq (6.4), the solution rapidly converges to the exact classical solution, i.e.,

    f(x,y,t)=4exp(xy+8t)1+exp(2(xy+8t)). (6.5)

    For the first case, we take the initial condition as

    f(x,y,0)=4exp(x+y)1+exp(2(xy)).

    Using the detail procedure as discussed above, we get

    f0(x,y,t)=f(x,y,0)=4exp(x+y)1+exp(2(xy)),f1(x,y,t)=L1[sΨ(1Ψ)+ΨsΨc(Ψ)L[6f20xf06f20yf0+3x3f03y3f032x2(yf0)+3x(2y2f0)]],

    using Mathematica, we obtain

    f1(x,y,t)=(1Ψ+ΨtΨΓ(Ψ+1))(32texp(x+y)(exp(2x)+exp(2y))(exp(2x)+exp(2y))2c(Ψ)),

    similarly, other terms can be calculated with the help of Mathematica. The required series solution is given as

    {f(x,y,t)=4exp(x+y)1+exp(2(xy))32t(1Ψ+ΨtΨΓ(Ψ+1))×(exp(x+y)(exp(2x)+exp(2y))(exp(2x)+exp(2y))2c(Ψ))+ (6.6)

    Remark 6.2. When we put Ψ=1 in Eq (6.6), the solution rapidly converges to the exact classical solution, i.e.,

    f(x,y,t)=4exp(xy+8t)1+exp(2(xy+8t)). (6.7)

    Here, we derive some results regarding to the convergence and stability of the proposed scheme with the help of functional analysis. The convergence of the proposed scheme of the is presented in the following theorem.

    Theorem 7.1. Let H be a Banach space and T:HHbe an operator. Suppose that f be the exact solution of the proposedequation. If ϖ such that 0ϖ<1 and fξ+1ϖfξ,ξN{0}, thenthe approximate solution ξ=0fξ convergesto the exact solution f.

    Proof. We construct a series as

    S0=f0,S1=f0+f1,S2=f0+f1+f2,Sξ=f0+f1+fξ.

    We want to prove that the sequence {Sξ}ξ=0 is a Cauchy sequence in H. Let us consider

    Sξ+1Sξ=fξϖfξϖ2fξ1ϖ3fξ2ϖξ+1f0.

    Now for every ξ,mN, we have

    SξSm=(SξSξ1)+(Sξ1Sξ2)++(Sm+1Sm)SξSξ1+Sξ1Sξ2++Sm+1Smϖξf0+ϖξ1f0++ϖm+1f0(ϖξ+1+ϖξ+2+)f0=ϖξ+11ϖf0.

    Now, limξ,mSξSm=0. This implies that {Sξ}ξ=0 is Cauchy sequence in H. So fH such that limξSξ=f. This ends the proof.

    Next, we present the Picard's X-stability of the proposed scheme in the following theorem.

    Theorem 7.2. Let X be a self-mapping which is definedas

    X(fξ(x,t))=fξ+1(x,t)=fξ(x,t)+L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[6f2ξfξx6f2ξfξy+fξxxx+fξyyy3fξxxy+3fξxyy]]. (7.1)

    The iteration is X-stable in L1(a,b), ifthe condition

    (2K+2L)(Φ21+Φ1Φ2+Φ2)1+M2N3+3O4+3P5<1, (7.2)

    is satisfied.

    Proof. With the help of Banach contraction theorem, first we show that the mapping X possesses a unique fixed point. For this, assume that the bounded iterations for (ξ,m)N×N. Let Φ1,Φ2>0 such that fξΦ1, and fmΦ2. Consider

    X(fξ(x,t))X(fm(x,t))=fξ(x,t)fm(x,t)+L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[2f2ξfξx2f2ξfξy+fξxxx+fξyyy3fξxxy+3fξxyy]]L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[2f2mfmx2f2mfmy+fmxxx+fmyyy3fmxxy+3fmxyy]]=fξ(x,t)fm(x,t)+L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[2K(f3ξf3m)6L(f3ξf3m)+M(fξfm)N(fξfm)3O(fξfm)+3P(fξfm)]].

    Now, using triangle inequality, we have

    X(fξ(x,t))X(fm(x,t))=fξ(x,t)fm(x,t)+L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[2K(f3ξf3m)2L(f3ξf3m)+M(fξfm)N(fξfm)3O(fξfm)+3P(fξfm)]].

    Using boundedness of fξ and fm, we have

    X(fξ(x,t))X(fm(x,t))[(2K+2L)(Φ21+Φ1Φ2+Φ2)1+M2N3+3O4+3P5]fξfm,

    where i,i=1,2,3,4,5, are functions obtained from L1[sΨ(1Ψ)+Ψ)sΨc(Ψ)L[]]. Using assumption (7.2), the mapping X fulfills the contraction condition. Hence by Banach fixed point result, X has a unique fixed point. Also, the mapping X satisfies the condition of Theorem 2.4 with

    a=0,b=(2K+2L)(Φ21+Φ1Φ2+Φ2)1+M2N3+3O4+3P5.

    Thus, the mapping X fulfills all conditions of Picard's X-stable. Hence our proposed scheme is Picard's X-stable.

    Thanks to the Mittag-Leffler kernel, which solves the singularity and locality problems with the Caputo and Caputo-Fabrizio FOD kernels. Since ABC-derivative is based on the Mittag-Leffler kernel, it has recently become popular for investigating the dynamics of a mathematical model that governs a physical process. We use ABC-derivative to investigate the soliton solution of the new modified KdV equation in (2+1) dimension in the current paper. Since the presence of a solution is essential for the study of a model, we have deduced some results using fixed point theory that guarantees at least one solution and the unique solution of the proposed equation. There are several techniques for solving FDEs, but among the analytical methods, LADM is the most effective and accurate. Under ABC-derivative, we used the LADM to obtain the solution of the proposed equation. Graphs of the solution are used to observe the method's convergence. The exact solution and the approximate solution obtained with the aid of LADM are in good agreement (See Figures 1 and 2). The dynamics of the solution under the ABC-derivative have been investigated. We can see from Figures 3 and 4 that the fractional-order solution curves are approaching the integer-order curve when fractional-order equals 1. Figure 5 shows that the ABC-derivative solution curve is much closer to the exact solution than the Caputo-Fabrizio derivative curve. As a result, the proposed model is better than Caputo-Fabrizio's. Figures 69 denote the graphical representation of solution obtained in the Case B. The considered equation will be studied under more generalized fractional operators in the next paper.

    Figure 1.  Two dimensional dynamics of exact and approximate solutions.
    Figure 2.  Dynamics of exact and approximate solutions in 3D.
    Figure 3.  Dynamics of approximate solution for different fractional orders in 2D.
    Figure 4.  Dynamics of approximate solution for different fractional orders in 3D.
    Figure 5.  Comparison between solution curves of exact solution, ABC-solution and Caputo-Fabrizio solution.
    Figure 6.  Dynamics of approximate solution in case B for different fractional orders in 2D.
    Figure 7.  Dynamics of approximate solution case B for different fractional orders in 3D.
    Figure 8.  Two dimensional dynamics of exact and approximate solutions in case B.
    Figure 9.  Dynamics of exact and approximate solutions in 3D in case B.

    The authors wish to express their sincere thanks to the honorable referees for their valuable comments and suggestions to improve the quality of the paper. In addition, authors would like to express their gratitude to the United Arab Emirates University, Al Ain, UAE for providing the financial support with Grant No. 12S005-UPAR 2020.

    In this research, there are no conflicts of interest.

    [1] [ J. P. Aparicio,C. Castillo-Chavez, Mathematical modelling of tuberculosis epidemics, Mathematical Biosciences and Engineering, 6 (2009): 209-237.
    [2] [ Z. Bai, Threshold dynamics of a time-delayed SEIRS model with pulse vaccination, Mathematical Biosciences, 269 (2015): 178-185.
    [3] [ S. M. Blower,T. Chou, Modeling the emergence of the 'hot zones': Tuberculosis and the amplification dynamics of drug resistance, Nature Medicine, 10 (2004): 1111-1116.
    [4] [ S. M. Blower,P. M. Small,P. C. Hopewell, Control strategies for tuberculosis epidemic: New models for old problems, Science, 273 (1996): 497-500.
    [5] [ S. Bowong,J. J. Tewa, Global analysis of a dynamical model for transmission of tuberculosis with a general contact rate, Communications in Nonlinear Science and Numerical Simulation, 15 (2010): 3621-3631.
    [6] [ H. Cao,Y. Zhou, The discrete age-structured SEIT model with application to tuberculosis transmission in China, Mathematical and Computer Modelling, 55 (2012): 385-395.
    [7] [ C. Castillo-Chavez,Z. Feng, To treat or not treat: The case of tuberculosis, Journal of Mathematical Biology, 35 (1997): 629-656.
    [8] [ C. Castillo-Chavez,Z. Feng, Global stability of an age-structure model for TB and its applications to optimal vaccination strategies, Mathematical Biosciences, 151 (1998): 135-154.
    [9] [ C. Castillo-Chavez,B. Song, Dynamical models of tuberculosis and their applications, Mathematical Biosciences and Engineering, 1 (2004): 361-404.
    [10] [ H. Chang, Quality monitoring and effect evaluation of BCG vaccination in neonatus, Occupation and Health, 9 (2013): 1109-1110.
    [11] [ O. Diekmann,J. A. P. Heesterbeek,J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, Journal of Mathematical Biology, 28 (1990): 365-382.
    [12] [ P. V. D. Driessche,J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Mathematical Biosciences, 180 (2002): 29-48.
    [13] [ M. J. Keeling,P. Rohani, null, Modeling Infectious Diseases in Humans and Animals, Princeton University Press, 2008.
    [14] [ J. Li, The spread and prevention of tuberculosis, Chinese Remedies and Clinics, 13 (2013): 482-483.
    [15] [ L. Liu,X. Zhao,Y. Zhou, A tuberculosis model with seasonality, Bulletin of Mathematical Biology, 72 (2010): 931-952.
    [16] [ J. S. Lopes, P. Rodrigues, S. T. Pinho, R. F. Andrade, R. Duarte and M. G. M. Gomes, Interpreting Measures of Tuberculosis Transmission: A Case Study on the Portuguese Population, BMC Infectious Diseases, 2014.
    [17] [ Z. Lu,X. Chi,L. Chen, The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission, Mathematical and Computer Modelling, 36 (2002): 1039-1057.
    [18] [ National Bureau of Statistics of China, Statistical Data of Tuberculosis 2004-2014. Available from: http://data.stats.gov.cn/easyquery.htm?cn=C01&zb=A0O0F01.
    [19] [ National Bureau of Statistics of China, China Statistical Yearbook 2014, Birth Rate, Death Rate and Natural Growth Rate of Population, 2014. Available from: http://www.stats.gov.cn/tjsj/ndsj/2014/indexch.htm.
    [20] [ National Technic Steering Group of the Epidemiological Sampling Survey for Tuberculosis, Report on fourth national epidemiological sampling survey of tuberculosis, Chinese Journal of Tuberculosis and Respiratory Diseases, 25 (2002), 3-7.
    [21] [ A. M. Samoilenko,N. A. Perestyuk, Periodic and almost-periodic solution of impulsive differential equations, Ukrainian Mathematical Journal, 34 (1982): 66-73,132.
    [22] [ A. M. Samoilenko and N. A. Perestyuk, Differential Equations with Impulse Effect, Visa Skola, Kiev, 1987.
    [23] [ O. Sharomi,C. N. Podder,A. B. Gumel,B. Song, Mathematical analysis of the transmission dynamics of HIV/TB coinfection in the presence of treatment, Mathematical Biosciences and Engineering, 5 (2008): 145-174.
    [24] [ M. Shen, Y. Xiao, W. Zhou and Z. Li, Global dynamics and applications of an epidemiological model for hepatitis C virus transmission in China, Discrete Dynamics in Nature and Society, 2015 (2015), Article ID 543029, 13pp.
    [25] [ B. Shulgin,L. Stone,Z. Agur, Pulse vaccination strategy in the SIR epidemic model, Bulletin of Mathematical Biology, 60 (1998): 1123-1148.
    [26] [ B. Song,C. Castillo-Chavez,J. P. Aparicio, Tuberculosis models with fast and slow dynamics: The role of close and casual contacts, Mathematical Biosciences, 180 (2002): 187-205.
    [27] [ Technical Guidance Group of the Fifth National TB Epidemiological Survey and The Office of the Fifth National TB Epidemiological Survey, The fifth national tuberculosis epidemiological survey in 2010, Chinese Journal of Antituberculosis, 34 (2012), 485–508.
    [28] [ E. Vynnycky,P. E. M. Fine, The long-term dynamics of tuberculosis and other diseases with long serial intervals: Implications of and for changing reproduction numbers, Epidemiology and Infection, 121 (1998): 309-324.
    [29] [ WHO, Global Tuberculosis Report 2015, World Health Organization, 2015. Available from: http://www.who.int/tb/publications/global_report/en/.
    [30] [ WHO, Tuberculosis Vaccine Development, 2015. Available from: http://www.who.int/immunization/research/development/tuberculosis/en/.
    [31] [ C. Xiong,X. Liang,H. Wang, A Systematic review on the protective efficacy of BCG against children tuberculosis meningitis and millet tuberculosis, Chinese Journal of Vaccines and Immunization, 15 (2009): 359-362.
    [32] [ B. Xu,Y. Hu,Q. Zhao,W. Wang,W. Jiang,G. Zhao, Molecular epidemiology of TB -Its impact on multidrug-resistant tuberculosis control in China, International Journal of Mycobacteriology, 4 (2015): 134.
    [33] [ Y. Yang,S. Tang,X. Ren,H. Zhao,C. Guo, Global stability and optimal control for a tuberculosis model with vaccination and treatment, Discrete and Continuous Dynamical Systems -Series B, 21 (2016): 1009-1022.
    [34] [ Y. Zhou,K. Khan,Z. Feng,J. Wu, Projection of tuberculosis incidence with increasing immigration trends, Journal of Theoretical Biology, 254 (2008): 215-228.
    [35] [ E. Ziv,C. L. Daley,S. M. Blower, Early therapy for latent tuberculosis infection, American Journal of Epidemiology, 153 (2001): 381-385.
  • This article has been cited by:

    1. Liu Xuan, Shabir Ahmad, Aman Ullah, Sayed Saifullah, Ali Akgül, Haidong Qu, Bifurcations, stability analysis and complex dynamics of Caputo fractal-fractional cancer model, 2022, 159, 09600779, 112113, 10.1016/j.chaos.2022.112113
    2. Mansour Mahmoud, Hanan Almuashi, On Some Asymptotic Expansions for the Gamma Function, 2022, 14, 2073-8994, 2459, 10.3390/sym14112459
    3. Umair Ali, Muhammad Naeem, Reham Alahmadi, Farah Aini Abdullah, Muhammad Asim Khan, Abdul Hamid Ganie, An investigation of a closed-form solution for non-linear variable-order fractional evolution equations via the fractional Caputo derivative, 2023, 11, 2296-424X, 10.3389/fphy.2023.1114319
    4. Hakeem Ullah, Mehreen Fiza, Ilyas Khan, Abd Allah A. Mosa, Saeed Islam, Abdullah Mohammed, Modifications of the Optimal Auxiliary Function Method to Fractional Order Fornberg-Whitham Equations, 2023, 136, 1526-1506, 277, 10.32604/cmes.2023.022289
    5. Asif Khan, Amir Ali, Shabir Ahmad, Sayed Saifullah, Kamsing Nonlaopon, Ali Akgül, Nonlinear Schrödinger equation under non-singular fractional operators: A computational study, 2022, 43, 22113797, 106062, 10.1016/j.rinp.2022.106062
    6. Fang Fang, Beibei Hu, Ling Zhang, Abdullahi Yusuf, Inverse Scattering Transform for the Generalized Derivative Nonlinear Schrödinger Equation via Matrix Riemann–Hilbert Problem, 2022, 2022, 1563-5147, 1, 10.1155/2022/3967328
    7. Ihtisham Ul Haq, Nigar Ali, Shabir Ahmad, Tayyaba Akram, Ibrahim Mahariq, A Hybrid Interpolation Method for Fractional PDEs and Its Applications to Fractional Diffusion and Buckmaster Equations, 2022, 2022, 1563-5147, 1, 10.1155/2022/2517602
    8. Tariq Mahmood, Fuad S. Al-Duais, Mei Sun, Dynamics of Middle East Respiratory Syndrome Coronavirus (MERS-CoV) involving fractional derivative with Mittag-Leffler kernel, 2022, 606, 03784371, 128144, 10.1016/j.physa.2022.128144
    9. Sajjad Ali, Aman Ullah, Shabir Ahmad, Kamsing Nonlaopon, Ali Akgül, Analysis of Kink Behaviour of KdV-mKdV Equation under Caputo Fractional Operator with Non-Singular Kernel, 2022, 14, 2073-8994, 2316, 10.3390/sym14112316
    10. Daba Meshesha Gusu, Megersa Danu, Abdullahi Yusuf, Existence of Solutions of Boundary Value Problem for Nonlinear One-Dimensional Wave Equations by Fixed Point Method, 2022, 2022, 1563-5147, 1, 10.1155/2022/5099060
    11. Naila Sajid, Zahida Perveen, Maasoomah Sadaf, Ghazala Akram, Muhammad Abbas, Thabet Abdeljawad, Manar A. Alqudah, Implementation of the Exp-function approach for the solution of KdV equation with dual power law nonlinearity, 2022, 41, 2238-3603, 10.1007/s40314-022-02047-2
    12. Suleyman Cetinkaya, Ali Demir, On the solution of time fractional initial value problem by a new method with ARA transform, 2023, 44, 10641246, 2693, 10.3233/JIFS-223237
    13. Shabir Ahmad, Aman Ullah, Ali Akgül, Fahd Jarad, A hybrid analytical technique for solving nonlinear fractional order PDEs of power law kernel: Application to KdV and Fornberg-Witham equations, 2022, 7, 2473-6988, 9389, 10.3934/math.2022521
    14. Shabir Ahmad, Sayed Saifullah, Analysis of the seventh-order Caputo fractional KdV equation: applications to the Sawada–Kotera–Ito and Lax equations, 2023, 75, 0253-6102, 085002, 10.1088/1572-9494/acded7
    15. Ihsan Ullah, Aman Ullah, Shabir Ahmad, Hijaz Ahmad, Taher A. Nofal, A survey of KdV-CDG equations via nonsingular fractional operators, 2023, 8, 2473-6988, 18964, 10.3934/math.2023966
    16. Panjaiyan Karthikeyann, Sadhasivam Poornima, Kulandhaivel Karthikeyan, Chanon Promsakon, Thanin Sitthiwirattham, Chang Phang, On Implicit Atangana–Baleanu–Caputo Fractional Integro-Differential Equations with Delay and Impulses, 2024, 2024, 2314-4785, 1, 10.1155/2024/5531984
    17. IBRAHIM SLIMANE, JUAN J. NIETO, SHABIR AHMAD, A FRACTIONAL-ORDER BOVINE BABESIOSIS EPIDEMIC TRANSMISSION MODEL WITH NONSINGULAR MITTAG-LEFFLER LAW, 2023, 31, 0218-348X, 10.1142/S0218348X23400339
    18. Yuqi Zhang, Peiluan Li, Changjin Xu, Xueqing Peng, Rui Qiao, Investigating the Effects of a Fractional Operator on the Evolution of the ENSO Model: Bifurcations, Stability and Numerical Analysis, 2023, 7, 2504-3110, 602, 10.3390/fractalfract7080602
    19. Ahmed E. Abouelregal, Faisal Alsharif, Hashem Althagafi, Yazeed Alhassan, Fractional heat transfer DPL model incorporating an exponential Rabotnov kernel to study an infinite solid with a spherical cavity, 2024, 9, 2473-6988, 18374, 10.3934/math.2024896
    20. Sivaranjani Ramasamy, Thangavelu Senthilprabu, Kulandhaivel Karthikeyan, Palanisamy Geetha, Saowaluck Chasreechai, Thanin Sitthiwirattham, Existence, uniqueness and controllability results of nonlinear neutral implicit ABC fractional integro-differential equations with delay and impulses, 2025, 10, 2473-6988, 4326, 10.3934/math.2025200
    21. Mariam Sultana, Uroosa Arshad, Ali Akgül, M. Khalid, Numerical Solutions of Fractional Systems Using Bessel Artificial Neural Network Based Integrated Intelligent Computing, 2025, 20, 1555-1415, 10.1115/1.4068236
  • Reader Comments
  • © 2017 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(5338) PDF downloads(707) Cited by(20)

Figures and Tables

Figures(5)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog