Processing math: 68%
Research article

Relation-theoretic almost ϕ-contractions with an application to elastic beam equations

  • In this article, we prove some results on existence and uniqueness of fixed points for an almost ϕ-contraction mapping defined on a metric space endowed with an amorphous relation. Our results generalize and improve several well known fixed point theorems of the existing literature. To substantiate the credibility of our results, we construct some examples. We also apply our results to determine a unique solution of a boundary value problem associated with nonlinear elastic beam equations.

    Citation: Ebrahem A. Algehyne, Nifeen Hussain Altaweel, Mounirah Areshi, Faizan Ahmad Khan. Relation-theoretic almost ϕ-contractions with an application to elastic beam equations[J]. AIMS Mathematics, 2023, 8(8): 18919-18929. doi: 10.3934/math.2023963

    Related Papers:

    [1] Yongyi Gu, Najva Aminakbari . Two different systematic methods for constructing meromorphic exact solutions to the KdV-Sawada-Kotera equation. AIMS Mathematics, 2020, 5(4): 3990-4010. doi: 10.3934/math.2020257
    [2] Jin Hyuk Choi, Hyunsoo Kim . Exact traveling wave solutions of the stochastic Wick-type fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation. AIMS Mathematics, 2021, 6(4): 4053-4072. doi: 10.3934/math.2021240
    [3] Jianhong Zhuang, Yaqing Liu, Ping Zhuang . Variety interaction solutions comprising lump solitons for the (2+1)-dimensional Caudrey-Dodd-Gibbon-Kotera-Sawada equation. AIMS Mathematics, 2021, 6(5): 5370-5386. doi: 10.3934/math.2021316
    [4] Saima Noor, Azzh Saad Alshehry, Asfandyar Khan, Imran Khan . Innovative approach for developing solitary wave solutions for the fractional modified partial differential equations. AIMS Mathematics, 2023, 8(11): 27775-27819. doi: 10.3934/math.20231422
    [5] Humaira Yasmin, Aljawhara H. Almuqrin . Efficient solutions for time fractional Sawada-Kotera, Ito, and Kaup-Kupershmidt equations using an analytical technique. AIMS Mathematics, 2024, 9(8): 20441-20466. doi: 10.3934/math.2024994
    [6] Huiqing Wang, Md Nur Alam, Onur Alp İlhan, Gurpreet Singh, Jalil Manafian . New complex wave structures to the complex Ginzburg-Landau model. AIMS Mathematics, 2021, 6(8): 8883-8894. doi: 10.3934/math.2021515
    [7] Noorah Mshary . Exploration of nonlinear traveling wave phenomena in quintic conformable Benney-Lin equation within a liquid film. AIMS Mathematics, 2024, 9(5): 11051-11075. doi: 10.3934/math.2024542
    [8] Naveed Iqbal, Muhammad Bilal Riaz, Meshari Alesemi, Taher S. Hassan, Ali M. Mahnashi, Ahmad Shafee . Reliable analysis for obtaining exact soliton solutions of (2+1)-dimensional Chaffee-Infante equation. AIMS Mathematics, 2024, 9(6): 16666-16686. doi: 10.3934/math.2024808
    [9] A. K. M. Kazi Sazzad Hossain, M. Ali Akbar . Solitary wave solutions of few nonlinear evolution equations. AIMS Mathematics, 2020, 5(2): 1199-1215. doi: 10.3934/math.2020083
    [10] Abdulghani R. Alharbi . Traveling-wave and numerical solutions to a Novikov-Veselov system via the modified mathematical methods. AIMS Mathematics, 2023, 8(1): 1230-1250. doi: 10.3934/math.2023062
  • In this article, we prove some results on existence and uniqueness of fixed points for an almost ϕ-contraction mapping defined on a metric space endowed with an amorphous relation. Our results generalize and improve several well known fixed point theorems of the existing literature. To substantiate the credibility of our results, we construct some examples. We also apply our results to determine a unique solution of a boundary value problem associated with nonlinear elastic beam equations.



    Since the dawn of mankind, humans have shown a keen interest in comprehending natural phenomena, starting with fire, lightning, thunder, earthquakes and volcanoes, extending all the way down to the nano particle. Many scientists have expressed interest in researching these events from the start of the 18th century. Until the middle of the 19th century, these studies are still being developed. The study of various scientific fields has benefited greatly from the use of PDEs, which have been used to describe a wide variety of events. Since the discovery of partial differential equations and their types, methods for solving them, and methods for determining their precise solutions, the world has been shrouded in the darkness of ignorance for a very long time, this darkness has started to dissipate. This has made it easier to comprehend these models and how to use them for global advancement, for keeping up with development, and for keeping up with the scientific renaissance.

    Since then, mathematicians have investigated the forms of nonlinear partial differential equations. Scientists and academics have devised several methods for solving PDEs in an effort to comprehend these phenomena. The precise and approximate solutions to these models can be found using a variety of techniques. The key findings in identifying explicit solutions to nonlinear partial differential equations (NLPDEs) were deduced in [1]. In the scientific and social sciences, NLPDEs are frequently used to describe complicated events. Modern science and engineering have produced several well-known models to demonstrate the dynamics of nonlinearity, including the Kortewegde Vries (KdV) equation [2], the KdV Burgers equation [3], the modified KdV (mKdV) equation [4,5], the modified KdV Kadomtsev Petviashvili (mKdVKP) equation [6] and others.

    The SK equation was created by Wazwaz [7] from the mKdV equation, and a few multiple solitons results were found using the condensed Hirota method. Later, utilizing the sine-cosine function and Kudryashov approaches, many researchers explored the mKdV model and discovered several kink, brilliant and periodic solutions in singular form. Dual wave solutions were created by the authors using modified Kudryashov and auxiliary equation approaches. The SK equation is a special instance of the mKdV equation [8]. In optical fibre, the variable-coefficients in NLPDEs can indicate group velocity dispersion, Kerr nonlinearity, third-order dispersion, time-delay related to the cubic term, tidal and wave energy and external potential [9,10,11,12]. To the best of our knowledge, the (1+1)-dimensional time-dependent variable-coefficient SK equation has not been researched before. The literature shows that there is potential for more research into the SK equation using different approaches, as well as demonstrating how their physical explanation works. The SK equation is one of the well-known partial differential equations. It is a significant mathematical model that is applied to numerous physical situations, including quantum physics, nonlinear optics and the travel of long waves in shallow water under gravity. Sawada and Kotera produced the (1+1)-dimensional SK equation, which established the presence of multiple solutions in [13]. The classical (1+1) dimensional SK equation is given as [14]

    ut+45u2ur15ururr15uurrr+urrrrr=0, (1.1)

    which represents a nonlinear physical model equation. Different approaches were examined with this important model equation[13,15,16,17,18,19,20]. The SK equation has a vast range from classical to bidirectional and dual-wave forms, which are used by a lot of researchers in their articles. Zeng used the Hirota bilinear method to construct the N-soliton solution of the (1+1)-dimensional variable-coefficient Sawada–Kotera equation[21], lump-type solutions to the bidirectional Sawada–Kotera equation[22], modified Kudryashov and new auxiliary equation methods to investigate the abundant solutions to the two-mode Sawada-Kotera equation[23], the trial function method to find the solitary wave solutions of bidirectional SK and variable-coefficient SK equations[24], exp(ψ(z))-expansion method to find meromorphic solutions to the KdV-SK equation[25], exact traveling wave solutions of the CDGSK equation using subequations method[26], Lie symmetry classification of the generalized PDE[27] and the lump solution of the CDGSK using bilinear operator[28].

    Soliton theory is essential for finding these nonlinear models' exact solutions because it can reveal important underlying details about phenomena like optical fibres, plasma physics, fluid mechanics, biophysics, solid-state physics, control theory, systems identification and many more. In numerous branches of applied research and engineering, the age of analytical solutions is becoming increasingly important to gain a deeper understanding of nonlinear processes and their potential practical uses. To realize the wave structures more effectively, the traveling wave structures plays a significant role.

    Because of this, eminent researchers of this decade are more eager to create any cutting-edge analytical answers for such models. Both fundamental theory and applications of nonlinear science heavily utilize NLPDE solitons. The soliton and solitary wave solutions give physical details and enhance the problem's mathematical aspects[29,30,31,32,33,34]. The extended direct algebraic method has been used previously for finding solutions of many equations, like to find optical solitons for quadratic-cubic nonlinearity[35], to find the solutions of cubic nonlinear Schrodinger equation[36] and for finding the solitary wave solutions of nonlinear evolution equations[37]. The extended direct algebraic method has not yet been applied to the SK equation. By displaying the discovered solutions in two- and three-dimensional graphs with various parameter values, this approach offers insightful information regarding the internal formation of the traveling wave phenomenon.

    According to the literature, there is need for additional research on the SK equation in order to develop new dual-wave solutions using the modified Kudryashov method (mKM) and the auxiliary equation method (AEM), as well as to derive their physical clarifications. The modified auxiliary equation method is an interesting special case of auxiliary equation method. In this study, a modest effort has been made to create some soliton solutions to the SK equation using the modified auxiliary equation method and the extended direct algebraic method. This effort has been prompted by the current literature. The suggested methods are among the best direct methods for creating precise solitary wave solutions to nonlinear partial differential equations. The strategy in this paper is summarized as follows: Section 2 offers a description of MAEM. In Section 3, the MAEM technique is explained in detail. The extended direct algebraic method is given in Section 4. Section 5 provides an explanation of the graphical representations. Section 6 contains the closing comments.

    Let u=u(r,t) satisfy the general equation given as

    G(u,ur,ut,urr,urt,utt,....)=0, (2.1)

    which can be converted to an ODE of the form

    H(u,u,u,...)=0. (2.2)

    The provided expansion is present in the general solution

    u(ζ)=α0+kj=1[αj(ϕg)j+βj(ϕg)j], (2.3)

    where α and β are arbitrary constants.

    The nature of the general solution g(ζ) to Eq (2.1) is determined by using an auxiliary equation given as

    g(ζ)=a+bϕj+cϕjlnϕ, (2.4)

    where ϕ>0,ϕ1. Also, a, b, c are arbitrary values. The roots of Eq (2.4) are given below.

    For Eq (2.4), the following situations occurred as

    ● When  b24ac<0,c0,

    ϕg(ζ)=b+4acb2 tan(4acb22)ζ2c,

    or

    ϕg(ζ)=b+4acb2 cot(4acb22)ζ2c.

    ● When  b24ac>0,c0,

    ϕg(ζ)=b+b24ac tanh(b24ac2)ζ2c,

    or

    ϕg(ζ)=b+b24ac coth(b24ac2)ζ2c.

    ● When  b24ac=0,c0,

    ϕg(ζ)=2+bζ2cζ.

    Consider the wave transformation for the SK equation

    u(r,t)=u(ζ),ζ=rd t, (3.1)

    where d is a constant. Substituting Eq (3.1) into Eq (1.1) and integrating, we get the ODE

    d u+15u315uu+u=0, (3.2)

    where u=dudζ and the constants of integration are equal to zero. Balancing the highest order derivative u and the highest order nonlinear term u3 in Eq (3.2), we get k=2 in accordance with the homogeneous balance principle. As a result, (2.3) produces the differential equation:

    u(ζ)=α0+α1ϕg+β1ϕg+α2ϕ2g+β2ϕ2g. (3.3)

    The following set of algebraic equations is constructed using Eqs (3.3) and (2.4) in Eq (3.2), and balancing the coefficients of each power of ϕg to zero.

    {(ϕg)6:15β2(8a46a2β2+β22)=0,(ϕg)5:3(8a4β1+8a2(14ab5β1)β2+5(10ab+3β1)β22)=0,(ϕg)4:15(abβ1(4a213β2)+β21(2a2+3β2)+β2(2a2(11b2+8ac3α0)+(4(b2+2ac)+3α0)β2))=0,(ϕg)3:5(9abβ21+3β31+β1(2a2(5b2+4ac3α0)3(5(b2+2ac)6α0)β2)+β2(2a(13b3+44abc15bα09aα1)+9(2bc+α1)β2))=0,(ϕg)2:15(b2+2ac3α0)β21+15β1(a(b3+4abc3bα02aα1)+(7bc+6α1)β2)β2(d8(2b4+29ab2c+17a2c2)+15(4(b2+2ac)α03α20+11abα1+α2(8a23β2)+2c2β2))=0,(ϕg)1:β1(d+b4+22ab2c+16a2c215((b2+2ac)α03α20+4a2α2+α1(4ab3β1)+bcβ1))+15(5(b2+2ac)α1+6α0(bc+α1)+2α2(8ab+3β1)+2c(b3+4abccβ1))β2=0,(ϕg)0:dα0+15α30+ab3α1+8a2bcα115abα0α1+14a2b2α2+16a3cα230a2α0α2+b3cβ1+8abc2β115bcα0β130b2α1β160acα1β190α0α1β1135abα2β1+45α2β21+14b2c2β2+16ac3β230c2α0β2135bcα1β2+45α21β2120b2α2β2240acα2β2+90α0α2β2=0,(ϕg)1:15α21(ab3β1)+15α2(2ab(b2+4ac)5(b2+2ac)β1+6α0(ab+β1)16bcβ2)+α1(d+b4+22ab2c+16a2c215((b2+2ac)α03α20+2α2(b23β2)+4c(bβ1+cβ2)))=0,(ϕg)2:15(b2+2ac3α0)α21+15α1(3acα0+α2(7ab+6β1)+r(b3+4abc2cβ1))α2(d8(2b4+29ab2c+17a2c2)+15(4(b2+2ac)α03α20+α2(2a23β2)+c(11bβ1+8cβ2)))=0,(ϕg)3:5(9bcα21+3α31+α1(2c2(5b2+4ac3α0)3(5(b2+2ac)6α0)α2)+α2(30bcα0+9α2(2ab+β1)+2c(13b3+44abc9cβ1)))=0,(ϕg)4:15(bcα1(4c213α2)+α21(2c2+3α2)+α2(2c2(11b2+8ac3α0)+(4(b2+2ac)+3α0)α2))=0,(ϕg)5:3(8c4α1+8c2(14bc5α1)α2+5(10bc+3α1)α22)=0,(ϕg)6:15α2(8c46c2α2+α22)=0.

    The system yields the following families.

    Family 1

    When α0=130(5b2+40ac54d+b48ab2c+16a2c2),α1=2bc,β1=0,α2=2c2,β2=0, the following cases are obtained.

    For b24ac>0 and n0,

    u1,1(r,t)=130(5b2+40ac54d+b48ab2c+16a2c2)+2bc(b+b24actanh(b24acζ2)2c)+2c2(b+b24actanh(b24acζ2)2c)2, (3.4)

    or

    u1,2(r,t)=130(5b2+40ac54d+b48ab2c+16a2c2)+2bc(b+b24accoth(b24acζ2)2c)+2c2(b+b24accoth(b24acζ2)2c)2, (3.5)

    which are the hyperbolic solutions.

    For b24ac<0 and n0,

    u1,3(r,t)=130(5b2+40ac54d+b48ab2c+16a2c2)+2bc(b+4acb2tan(4acb2ζ2)2c)+2c2(b+4acb2tan(4acb2ζ2)2c)2, (3.6)

    or

    u1,4(r,t)=130(5b2+40ac54d+b48ab2c+16a2c2)+2bc(b+4acb2cot(4acb2ζ2)2c)+2c2(b+4acb2cot(4acb2ζ2)2c)2, (3.7)

    which are trigonometric solutions.

    Family 2

    When α0=130(5b2+40ac+54d+b48ab2c+16a2c2),α1=2bc,β1=0,α2=2c2,β2=0, the following cases occurred.

    For b24ac>0 and n0,

    u2,1(r,t)=130(5b2+40ac+54d+b48ab2c+16a2c2)+2bc(b+b24actanh(b24acζ2)2c)+2c2(b+b24actanh(b24acζ2)2c)2, (3.8)

    or

    u2,2(r,t)=130(5b2+40ac+54d+b48ab2c+16a2c2)+2bc(b+b24accoth(b24acζ2)2c)+2c2(b+b24accoth(b24acζ2)2c)2, (3.9)

    which are the hyperbolic solutions.

    For b24ac<0 and n0,

    u2,3(r,t)=130(5b2+40ac+54d+b48ab2c+16a2c2)+2bc(b+4acb2tan(4acb2ζ2)2c)+2c2(b+4acb2tan(4acb2ζ2)2c)2, (3.10)

    or

    u2,4(r,t)=130(5b2+40ac+54d+b48ab2c+16a2c2)+2bc(b+4acb2cot(4acb2ζ2)2c)+2c2(b+4acb2cot(4acb2ζ2)2c)2, (3.11)

    which are trigonometric solutions.

    Family 3

    When α0=130(5b2+40ac54d+b48ab2r+16a2c2),α1=0,β1=2ab,α2=0,β2=2a2, the following cases are made.

    For b24ac>0 and n0

    u3,1(r,t)=130(5b2+40ac54d+b48ab2r+16a2c2)+2ab(b+b24actanh(b24acζ2)2c)1+2a2(b+b24actanh(b24acζ2)2c)2, (3.12)

    or

    u3,2(r,t)=130(5b2+40ac54d+b48ab2r+16a2c2)+2ab(b+b24accoth(b24acζ2)2c)1+2a2(b+b24accoth(b24acζ2)2c)2, (3.13)

    which are the hyperbolic solutions.

    For b24ac<0 and n0,

    u3,3(r,t)=130(5b2+40ac54d+b48ab2r+16a2c2)+2ab(b+4acb2tan(4acb2ζ2)2c)1+2a2(b+4acb2tan(4acb2ζ2)2c)2, (3.14)

    or

    u3,4(r,t)=130(5b2+40ac54d+b48ab2r+16a2c2)+2ab(b+4acb2cot(4acb2ζ2)2c)1+2a2(b+4acb2cot(4acb2ζ2)2c)2, (3.15)

    which are the trigonometric solutions.

    Family 4

    When α0=130(5b2+40ac+54d+b48ab2r+16a2c2),α1=0,β1=2ab,α2=0,β2=2a2, the following cases are obtained.

    For b24ac>0 and n0,

    u4,1(r,t)=130(5b2+40ac+54d+b48ab2r+16a2c2)+2ab(b+b24actanh(b24acζ2)2c)1+2a2(b+b24actanh(b24acζ2)2c)2, (3.16)

    or

    u4,2(r,t)=130(5b2+40ac+54d+b48ab2r+16a2c2)+2ab(b+b24accoth(b24acζ2)2c)1+2a2(b+b24accoth(b24acζ2)2c)2, (3.17)

    which are the hyperbolic solutions.

    For b24ac<0 and n0,

    u4,3(r,t)=130(5b2+40ac+54d+b48ab2r+16a2c2)+2ab(b+4acb2tan(4acb2ζ2)2c)1+2a2(b+4acb2tan(4acb2ζ2)2c)2, (3.18)

    or

    u4,4(r,t)=130(5b2+40ac+54d+b48ab2r+16a2c2)+2ab(b+4acb2cot(4acb2ζ2)2c)1+2a2(b+4acb2cot(4acb2ζ2)2c)2, (3.19)

    which are the trigonometric solutions.

    Assume a general NLPDE of the type

    X(u,ut,ur,utt,urr,urt,....)=0, (4.1)

    which can be transformed as

    F(R,R,R,...)=0, (4.2)

    by using following transformation as

    U(r,t)=U(ζ)eιΦ, (4.3)

    where ζ=rdt and Φ=ζ1r+ζ2t, and the order of differentiation for various properties in Eq (4.2) is represented by the prime sign. Assume that Eq (4.2) has a solution of the form

    U(ζ)=lj=1al(P(ζ))l, (4.4)

    where

    P(ζ)=lnχ(ν+μP(ζ)+λ(P(ζ))2),χ0,1. (4.5)

    In addition, λ, μ and ν are real constants. The general form of the solutions of Eq (4.5) w.r.t real parameters λ, μ and ν are as follows:

    1) For μ24νλ<0 and λ0,

    P1(ζ)=μ2λ+δ2λtanχ(δ2ζ), (4.6)
    P2(ζ)=μ2λδ2λcotχ(δ2ζ), (4.7)
    P3(ζ)=μ2λ+δ2λ(tanχ(δζ)±mnsecχ(δζ)),
    P4(ζ)=μ2λ+δ2λ(cotχ(δζ)±mncscχ(δζ)),
    P5(ζ)=μ2λ+δ4λ(tanχ(δ4ζ)cotχ(δ4ζ)).

    where δ=μ24νλ.

    2) For μ24νλ>0 and λ0,

    P6(ζ)=μ2λδ2λtanhχ(δ2ζ),
    P7(ζ)=μ2λδ2λcothχ(δ2ζ),
    P8(ζ)=μ2λ+δ2λ(tanhχ(δζ)±ιmnχ(δζ)),
    P9(ζ)=μ2λ+δ2λ(cothχ(δζ)±mnχ(δζ)),
    P10(ζ)=μ2λδ4λ(tanhχ(δ4ζ)+cothχ(δ4ζ)).

    where δ=μ24νλ.

    3) For νλ>0 and μ=0,

    P11(ζ)=νλtanχ(νλζ),
    P12(ζ)=νλcotχ(νλζ),
    P13(ζ)=νλ(tanχ(2νλζ)±mnsecχ(2νλζ)),
    P14(ζ)=νλ(cotχ(2νλζ)±mncscχ(2νλζ)),
    P15(ζ)=12νλ(tanχ(νλ2ζ)+cotχ(νλ2ζ)).

    4) For νλ<0 and μ=0,

    P16(ζ)=νλtanhχ(νλζ),
    P17(ζ)=νλcothχ(νλζ),
    P18(ζ)=νλ(tanhχ(2νλζ)±ιmnχ(2νλζ)),
    P19(ζ)=νλ(cothχ(2νλζ)±mnχ(2νλζ)),
    P20(ζ)=12νλ(tanhχ(νλ2ζ)+cothχ(νλ2ζ)).

    5) For μ=0 and ν=λ,

    P21(ζ)=tanχ(νζ),
    P22(ζ)=cotχ(νζ),
    P23(ζ)=tanχ(2νζ)±mnsecχ(2νζ),
    P24(ζ)=cotχ(2νζ)±mncscχ(2νζ),
    P25(ζ)=12(tanχ(ν2ζ)cotχ(ν2ζ)).

    6) For μ=0 and λ=ν,

    P26(ζ)=tanhχ(νζ),
    P27(ζ)=cothχ(νζ),
    P28(ζ)=tanhχ(2νζ)±ιmnχ(2νζ),
    P29(ζ)=cothχ(2νζ)±mnχ(2νζ),
    P30(ζ)=12(tanhχ(ν2ζ)+cothχ(ν2ζ)).

    7) For μ2=4νλ,

    P31(ζ)=2ν(μζln(χ)+2)μ2ζln(χ).

    8) For μ=p, ν=pq, (q0), and λ=0,

    P32(ζ)=χpζq.

    9) For μ=λ=0,

    P33(ζ)=νζln(χ).

    10) For μ=ν=0,

    P34(ζ)=1λζln(χ).

    11) For μ0 and ν=0,

    P35(ζ)=mμλ(coshχ(μζ)sinhχ(μζ)+m),
    P36(ζ)=μ(sinhχ(μζ)+coshχ(μζ))ξ(sinhχ(μζ)+coshχ(μζ)+n).

    12) For μ=p, λ=pq, (q0 and ν=0),

    P37(ζ)=mχpζmqnχpζ.

    sinhχ(ζ)=mχζnχζ2, coshχ(ζ)=mχζ+nχζ2,

    tanhχ(ζ)=mχζnχζmχζ+nχζ,

    cscχ(ζ)=2mχζnχζ, secχ(ζ)=2mχζ+nχζ,

    cothχ(ζ)=mχζ+nχζmχζnχζ,

    sinχ(ζ)=mχιζnχιζ2ι, cosχ(ζ)=mχιζ+nχιζ2,

    tanχ(ζ)=ιmχιζnχιζmχιζ+nχιζ,

    cscχ(ζ)=2ιmχζnχζ, secχ(ζ)=2mχζ+nχζ,

    cotχ(ζ)=ιmχιζ+nχιζmχιζnχιζ,

    where the terms "deformation parameters" m and n refer to arbitrary constants with values larger than 0.

    We employ an ansatz in Eq (1.1) to determine the solution

    U(r,t)=U(ζ), (4.8)

    and ζ=rdt, where the coefficient of r is the wave number and d is the speed of soliton. By imposing the considered transformation in Eq (4.8) on Eq (1.1), we produce the following real parts:

    dU+15U315UU+U=0. (4.9)

    The homogeneous balancing constant in Eq (4.9) is 2, and thus the solution is given as follows:

    U(ζ)=α0+α1P(ζ)+α2(P(ζ))2, (4.10)

    where

    P(ζ)=ln(χ)(ν+μP+λ(P(ζ))2). (4.11)

    By entering in the solution of Eq (4.10) along with Eqs (4.9)–(4.11) and by calculating the coefficients of powers of P(ζ), we get the following system of equations. This system of algebraic equations was solved using Wolfram Mathematica 12 software. The results are as follows:

    {α0=13λμ2Δ+83νΔ,α1=4μΔ,α2=4λΔ, (4.12)

    where Δ=13λlogχ2.

    We get the general solution of Eq (1.1) by plugging Eq (4.12) into Eq (4.10):

    U(r,t)=Δ3λ(μ2+8λν+12λ(μPi(ζ)+λPi(ζ)2)). (4.13)

    It should be seen that by choosing different values of Pj for the above equations, we get many solutions.

    (1) For μ24νλ<0 and λ0,

    U1(r,t)=83λνΔ2μ2Δ3λδΔλtan2χ(δ2ζ),
    U2(r,t)=23λμ2Δ+83νΔΔδλcot2χ(δ2ζ),
    U3(r,t)=2μ2Δ3λ+8νΔ3δΔλ(tanχ(δζ)±mnsecχ(δζ))2, (4.14)
    U4(r,t)=2μ2Δ3λ+8νΔ3δΔλ(cotχ(δζ)±mncscχ(δζ))2,
    U5(r,t)=2μ2Δ3λ+8νΔ3δΔ4λ(tanχ(δ4ζ)cotχ(δ4ζ))2.

    (2) For μ24νλ>0 and λ0,

    U6(r,t)=2μ2Δ3λ+8νΔ3+δΔλtanh2χ(δ2ζ),
    U7(r,t)=2μ2Δ3λ+8νΔ3+δΔλcoth2χ(δ2ζ),
    U8(r,t)=2μ2Δ3λ+8νΔ3+δΔλ(tanhχ(δζ)±ιmnsechχ(δζ))2,
    U9(r,t)=2μ2Δ3λ+8νΔ3+δΔλ(cothχ(δζ)±ιmncschχ(δζ))2,
    U10(r,t)=2μ2Δ3λ+8νΔ3+δΔλ(tanhχ(δ4ζ)+cothχ(δ4ζ))2.

    (3) For νλ>0 and μ=0,

    U11(r,t)=13λμ2Δ+83νΔ+4μΔνλtanχ(νλζ)+4λνtan2χ(νλζ),
    U12(r,t)=13λμ2Δ+83νΔ4μΔνλcotχ(νλζ)+4λνcot2χ(νλζ),
    U13(r,t)=Δ(μ23λ+8ν3+2μ2νλ(tanχ(2νλζ)±mnsecχ(2νλζ))+4ν(tanχ(2νλζ)±mnsecχ(2νλζ))2),
    U14(r,t)=μ2Δ3λ+8νΔ3+4μΔνλ(cotχ(2νλζ)±mncscχ(2νλζ))+4νΔ(cotχ(2νλζ)±mncscχ(2νλζ))2,
    U15(r,t)=μ2Δ3λ+8νΔ3+2μΔνλ(tanχ(νλ2ζ)+cotχ(νλ2ζ))+2νΔ(tanχ(νλ2ζ)+cotχ(νλ2ζ))2.

    (4) For νλ<0 and μ=0,

    U16(r,t)=13λμ2Δ+83νΔ4μΔ(νλtanhχ(νλζ))4νΔ(tanhχ(νλζ))2,
    U17(r,t)=13λμ2Δ+83νΔ4μΔ(νλcothχ(νλζ))4νΔ(cothχ(νλζ))2,
    U18(r,t)=13λμ2Δ+83νΔ+4μΔνλ(tanhχ(2νλζ)±ιmnsechχ(2νλζ))4νΔ(tanhχ(2νλζ)±ιmnsechχ(2νλζ))2,
    U19(r,t)=13λμ2Δ+83νΔ+4μΔνλ(cothχ(2νλζ)±ιmncschχ(2νλζ))4νΔ(cothχ(2νλζ)±ιmncschχ(2νλζ))2,
    U20(r,t)=13λμ2Δ+83νΔ+4μΔ(12νλ(tanhχ(νλ2ζ)+cothχ(νλ2ζ)))+4λΔ(12νλ(tanhχ(νλ2ζ)+cothχ(νλ2ζ)))2.

    (5) For μ=0 and ν=λ,

    U21(r,t)=13λμ2Δ+83νΔ+4μΔ(tanχ(νζ))+4λΔ(tanχ(νζ))2,
    U22(r,t)=13λμ2Δ+83νΔ+4μΔ(cotχ(νζ))+4λΔ(cotχ(νζ))2,
    U23(r,t)=13λμ2Δ+83νΔ+4μΔ(tanχ(2νζ)±mnsecχ(2νζ))+4λΔ(tanχ(2νζ)±mnsecχ(2νζ))2,
    U24(r,t)=13λμ2Δ+83νΔ+4μΔ(cotχ(2νζ)±mncscχ(2νζ))+4λΔ(cotχ(2νζ)±mncscχ(2νζ))2,
    U25(r,t)=13λμ2Δ+83νΔ+4μΔ(12(tanχ(ν2ζ)cotχ(ν2ζ)))+4λΔ(12(tanχ(ν2ζ)cotχ(ν2ζ)))2.

    (6) For μ=0 and ν=λ,

    U26(r,t)=13λμ2Δ+83νΔ+4μΔ(tanhχ(νζ))+4λΔ(tanhχ(νζ))2,
    U27(r,t)=13λμ2Δ+83νΔ+4μΔ(cothχ(νζ))+4λΔ(cothχ(νζ))2,
    \begin{eqnarray} \label{E4.80} U_{28}(r,t) = \frac{1}{{3\lambda }}{\mu ^2}\Delta + \frac{8}{3}\nu \Delta + 4\mu \Delta ( - {\tanh _\chi }\left( {2\nu \zeta } \right) \pm \iota \sqrt {mn} {\coth _\chi }\left( {2\nu \zeta } \right)) \\ + 4\lambda \Delta {( - {\tanh _\chi }\left( {2\nu \zeta } \right) \pm \iota \sqrt {mn} {\coth _\chi }\left( {2\nu \zeta } \right))^2}, \end{eqnarray}
    \begin{eqnarray} \label{E4.81} U_{29}(r,t) = \frac{1}{{3\lambda }}{\mu ^2}\Delta + \frac{8}{3}\nu \Delta + 4\mu \Delta ( - {\coth _\chi }\left( {2\nu \zeta } \right) \pm \iota \sqrt {mn} {\csc h _\chi }\left( {2\nu \zeta } \right)) \\ + 4\lambda \Delta {( - {\coth _\chi }\left( {2\nu \zeta } \right) \pm \iota \sqrt {mn} {\csc h _\chi }\left( {2\nu \zeta } \right))^2}, \end{eqnarray}
    \begin{eqnarray} \label{E4.82} U_{30}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta(-\frac{1}{{2}}\left(\tanh_{\chi}\left( \frac{\nu}{2}\zeta \right) + \coth_{\chi}\left(\frac{\nu}{2}\zeta \right)\right)) \\+4\lambda \Delta(-\frac{1}{{2}}\left(\tanh_{\chi}\left( \frac{\nu}{2}\zeta \right) + \coth_{\chi}\left( \frac{\nu}{2}\zeta \right)\right))^2. \end{eqnarray}

    (7) For \mu^2 = 4\nu\lambda,

    \begin{equation} \label{E4.83} U_{31}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta(\frac{-2\nu(\mu \zeta \ln(\chi)+2)}{\mu^2\zeta \ln(\chi)})+4\lambda \Delta(\frac{-2\nu(\mu\zeta \ln(\chi)+2)}{\mu^2\zeta \ln(\chi)})^2.\nonumber \end{equation}

    (8) For \mu = p, \nu = pq, (q\neq0), and\; \lambda = 0,

    \begin{equation} \label{E4.84} U_{32}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta( \chi^{p\zeta}-q)+4\lambda \Delta( \chi^{p\zeta}-q)^2.\nonumber \end{equation}

    (9) For \mu = \lambda = 0,

    \begin{equation} \label{E4.85} U_{33}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta(\nu \zeta\ln(\chi))+4\lambda \Delta(\nu \zeta\ln(\chi))^2.\nonumber \end{equation}

    (10) For \mu = \nu = 0,

    \begin{equation} \label{E4.86} U_{34}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta(\frac{-1}{\lambda \zeta\ln(\chi)})+4\lambda \Delta(\frac{-1}{\lambda \zeta\ln(\chi)})^2.\nonumber \end{equation}

    (11) For \mu\neq0 and \nu = 0,

    \begin{eqnarray} \label{E4.87} U_{35}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta(-\frac{m\mu}{\lambda\left(\cosh_{\chi}(\mu \zeta)- \sinh_{\chi}(\mu \zeta)+ m\right)}) \\+4\lambda\Delta(-\frac{m\mu}{\lambda\left(\cosh_{\chi}(\mu \zeta)- \sinh_{\chi}(\mu \zeta)+ m\right)})^2, \end{eqnarray}
    \begin{eqnarray} \label{E4.88} U_{36}(r,t) = \frac{1}{3 \lambda}\mu^2 \Delta+ \frac{8}{3} \nu \Delta+4\mu \Delta (-\frac{\mu\left(\sinh_{\chi}(\mu \zeta) + \cosh_{\chi}(\mu \zeta)\right)}{\lambda\left(\sinh_{\chi}(\mu \zeta) + \cosh_{\chi}(\mu \zeta)+ n\right)}) \\+4\lambda \Delta(-\frac{\mu\left(\sinh_{\chi}(\mu \zeta) + \cosh_{\chi}(\mu\zeta)\right)}{\lambda\left(\sinh_{\chi}(\mu \zeta) + \cosh_{\chi}(\mu \zeta)+ n\right)})^2. \end{eqnarray}

    (12) For \mu = p , \lambda = pq , ( q\neq0 and \nu = 0 ),

    \begin{equation} \label{E4.89} U_{37}(r,t) = \frac{{{\mu ^2}\Delta }}{{3\lambda }} + \frac{{8\nu \Delta }}{3} + (\frac{{4\mu \Delta m{\chi ^{p\zeta }}}}{{qn{\chi ^{p\zeta }} - m}}) + {(\frac{{2\sqrt {\lambda \Delta } m{\chi ^{p\zeta }}}}{{m - qn{\chi ^{p\zeta }}}})^2}.\nonumber \end{equation}

    The SK equation is a nonlinear equation. The wave behaviour connected to the SK equation is shown using the MAEM, which provides the equation's solitary wave solutions. Several of the answers that have been found are shown graphically in figures to more fully understand the dynamics of the solutions created. Figure 1 depicts the solitary waves associated with {|u_{1, 1}|} . The 3D plot shows that the variation in the fractional order has a significant influence on the exchange of amplitude and soliton nature. It demonstrates the spatial and temporal behavior of U(r, t) in the situation of a fractional order system. The behavior indicates that the solitary wave solution's wave nature varies as the fractional parameter is increased.

    Figure 1.  3D and 2D plots for {|u_{1, 1}|} when d = 1 , a = 0.1 , b = 1 , c = 0.1 .

    Figure 2 shows the solitary waves associated with {|u_{1, 3}|} with parametric values d = 1 , a = 0.1 , b = 0.01 , c = 0.1 , while Figure 3 illustrates the solution for {|u_{2, 1}|} having parametric values d = 1 , a = 0.1 , b = 1 , c = 0.1 .

    Figure 2.  3D and 2D plots for {|u_{1, 3}|} when d = 1 , a = 0.1 , b = 0.01 , c = 0.1 .
    Figure 3.  3D and 2D plots for {|u_{2, 1}|} when d = 1 , a = 0.1 , b = 1 , c = 0.1 .

    Figure 4 depicts the solitary waves associated with {|u_{2, 3}|} having parametric values d = 1 , a = 0.1 , b = 0.01 , c = 0.1 . Figure 5 shows the solitary waves associated with {|u_{3, 1}|} with parametric values d = 1 , a = 0.1 , b = 1 , c = 0.1 , while Figure 6 illustrates the solution for {|u_{4, 1}|} having parametric values d = 1 , a = 0.1 , b = 1 , c = 0.1 .

    Figure 4.  3D and 2D plots for {|u_{2, 3}|} when d = 1 , a = 0.1 , b = 0.01 , c = 0.1 .
    Figure 5.  3D and 2D plots for {|u_{3, 1}|} when d = 1 , a = 0.1 , b = 1 , c = 0.1 .
    Figure 6.  3D and 2D plots for {|u_{4, 1}|} when d = 1 , a = 0.1 , b = 1 , c = 0.1 .

    The propagation of waves and the contour graphs are also shown in the results. Solutions obtained by the extended direct algebraic method are also shown by 3D, 2D and contour graphs. For Eq (4.9), the real part of the solution is presented in Figure 7 by using the parameters \chi = 3, m = 9, n = 1, a = 0.1, b = 1, \lambda = 0.3, \nu = 0.2, \zeta_1 = 2, \mu = 0.3, \xi = 5 , for \alpha = 0.002 , \alpha = 0.08 and \alpha = 0.2 respectively. For Eq (4.9), the imaginary part of the solution is presented in Figure 8 by using the parameters \chi = 3, m = 9, n = 1, a = 0.1, b = 1, \lambda = 0.3, \nu = 0.2, \zeta_1 = 2, \mu = 0.3, \xi = 5 , for \alpha = 0.002 , \alpha = 0.08 and \alpha = 0.2 respectively. For Eq (4.9), the imaginary part of the solution is presented in Figure 9 by using the parameters \chi = 3, m = 9, n = 1, a = 0.1, b = 1, \lambda = 0.3, \nu = 0.2, \zeta_1 = 2, \mu = 0.3, \xi = 5 , for \alpha = 1 , \alpha = 0.5 and \alpha = 0.7 respectively.

    Figure 7.  For the real part of solution { U_3 }, 3D, contour and 2D plots with different values of \alpha .
    Figure 8.  For the imaginary part of solution { U_3 }, 3D, contour and 2D plots with different values of \alpha .
    Figure 9.  For the imaginary part of solution { U_3 }, 3D, contour and 2D plots with different values of \alpha .

    In this study, we used the modified auxiliary equation method and the extended direct algebraic method to locate distinct generalized solitary solutions to the SK equations. These results are validated for accuracy using symbolic computation tools on a computer. In contrast to the majority of the already used methodologies, the modified auxiliary equation method can produce creative and more general solutions. This exemplifies the value and immense potential of the modified auxiliary equation method in resolving difficult solitary wave issues that arise in mathematical physics. Using a modified auxiliary equation method, the wave solutions to the Sawada-Kotera problem are examined. Brilliant periodic waves, steady bright soliton waves and symmetric waves all have physical qualities that have been identified. Graphics are used to show the effects of the parameters.

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

    The author Magda Abd El-Rahman extends their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through large group Research Project under grant number RGP2/69/44. The author Muhammad Bilal Riaz thankful to Ministry of Education, Youth and Sports of the Czech Republic for their support through the e–INFRA CZ (ID: 90254).

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



    [1] F. E. Browder, On the convergence of successive approximations for nonlinear functional equations, Indag. Math., 30 (1968), 27–35.
    [2] D. W. Boyd, J. S. W. Wong, On nonlinear contractions, Proc. Amer. Math. Soc., 20 (1969), 158–164.
    [3] J. Matkowski, Integrable solutions of functional equations, Warszawa: Instytut Matematyczny Polskiej Akademi Nauk, 1975.
    [4] V. Berinde, Iterative approximation of fixed points, Springer Berlin, Heidelberg, 2007. https://doi.org/10.1007/978-3-540-72234-2
    [5] V. Berinde, Approximating fixed points of weak contractions using the Picard iteration, Nonlinear Anal. Forum, 9 (2004), 43–53.
    [6] V. Berinde, M. Pǎcurar, Fixed points and continuity of almost contractions, Fixed Point Theor., 9 (2008), 23–34.
    [7] M. Berinde, V. Berinde, On a general class of multi-valued weakly Picard mappings, J. Math. Anal. Appl., 326 (2007), 772–782. https://doi.org/10.1016/j.jmaa.2006.03.016 doi: 10.1016/j.jmaa.2006.03.016
    [8] M. A. Alghamdi, V. Berinde, N. Shahzad, Fixed points of non-self almost contractions, Carpathian J. Math., 30 (2014), 7–14.
    [9] G. V. R. Babu, M. L. Sandhya, M. V. R. Kameshwari, A note on a fixed point theorem of Berinde on weak contractions, Carpathian J. Math., 24 (2008), 8–12.
    [10] A. Alam, M. Imdad, Relation-theoretic contraction principle, J. Fixed Point Theory Appl., 17 (2015), 693–702. https://doi.org/10.1007/s11784-015-0247-y doi: 10.1007/s11784-015-0247-y
    [11] A. Alam, M. Imdad, Relation-theoretic metrical coincidence theorems, Filomat, 31 (2017), 4421–4439. https://doi.org/10.2298/FIL1714421A doi: 10.2298/FIL1714421A
    [12] A. Alam, M. Imdad, Nonlinear contractions in metric spaces under locally T-transitive binary relations, Fixed Point Theor., 19 (2018), 13–24. https://doi.org/10.24193/fpt-ro.2018.1.02 doi: 10.24193/fpt-ro.2018.1.02
    [13] M. Arif, M. Imdad, A. Alam, Fixed point theorems under locally T-transitive binary relations employing Matkowski contractions, Miskolc Math. Notes, 23 (2022), 71–83. https://doi.org/10.18514/MMN.2022.3220 doi: 10.18514/MMN.2022.3220
    [14] E. A. Algehyne, M. S. Aldhabani, F. A. Khan, Relational contractions involving (c)-comparison functions with applications to boundary value problems, Mathematics, 11 (2023), 1127. https://doi.org/10.3390/math11061277 doi: 10.3390/math11061277
    [15] F. A. Khan, Almost contractions under binary relations, Axioms, 11 (2022), 441. https://doi.org/10.3390/axioms11090441 doi: 10.3390/axioms11090441
    [16] A. Alam, R. George, M. Imdad, Refinements to relation-theoretic contraction principle, Axioms, 11 (2022), 316. https://doi.org/10.3390/axioms11070316 doi: 10.3390/axioms11070316
    [17] K. Sawangsup, W. Sintunavarat, A. F. Roldán-López-de-Hierro, Fixed point theorems for F_\mathcal{R}-contractions with applications to solution of nonlinear matrix equations, J. Fixed Point Theory Appl., 19 (2017), 1711–1725. https://doi.org/10.1007/s11784-016-0306-z doi: 10.1007/s11784-016-0306-z
    [18] H. H. Al-Sulami, J. Ahmad, N. Hussain, A. Latif, Relation-theoretic (\theta, \mathcal{R})-contraction results with applications to nonlinear matrix equations, Symmetry, 10 (2018), 767. https://doi.org/10.3390/sym10120767 doi: 10.3390/sym10120767
    [19] S. Shukla, N. Dubey, Some fixed point results for relation theoretic weak \varphi-contractions in cone metric spaces equipped with a binary relation and application to the system of Volterra type equation, Positivity, 24 (2020), 1041–1059. https://doi.org/10.1007/s11117-019-00719-8 doi: 10.1007/s11117-019-00719-8
    [20] B. S. Choudhury, N. Metiya, S. Kundu, Existence, well-posedness of coupled fixed points and application to nonlinear integral equations, CUBO, 23 (2021), 171–190.
    [21] C. Zhai, M. Hao, Fixed point theorems for mixed monotone operators with perturbation and applications to fractional differential equation boundary value problems, Nonlinear Anal., 75 (2012), 2542–2551. https://doi.org/10.1016/j.na.2011.10.048 doi: 10.1016/j.na.2011.10.048
    [22] M. Subaşi, S. I. Araz, Numerical regularization of optimal control for the coefficient function in a wave equation, Iran J. Sci. Technol. Trans. Sci., 43 (2019), 2325–2333. https://doi.org/10.1007/s40995-019-00690-9 doi: 10.1007/s40995-019-00690-9
    [23] A. Atangana, S. I. Araz, An accurate iterative method for ordinary differential equations with classical and Caputo-Fabrizio derivatives, 2023, hal-03956673.
    [24] J. Wu, Y. Liu, Fixed point theorems for monotone operators and applications to nonlinear elliptic problems, Fixed Point Theory Appl., 134 (2013), 134. https://doi.org/10.1186/1687-1812-2013-134 doi: 10.1186/1687-1812-2013-134
    [25] M. Jleli, V. C. Rajić, B. Samet, C. Vetro, Fixed point theorems on ordered metric spaces and applications to nonlinear elastic beam equations, J. Fixed Point Theory Appl., 12 (2012), 175–192. https://doi.org/10.1007/s11784-012-0081-4 doi: 10.1007/s11784-012-0081-4
    [26] S. Lipschutz, Schaum's outlines of theory and problems of set theory and related topics, McGraw-Hill, 1998.
    [27] B. Kolman, R. C. Busby, S. C. Ross, Discrete mathematical structures, 6 Eds, Pearson/Prentice Hall, 2009.
  • This article has been cited by:

    1. Sajawal Abbas Baloch, Muhammad Abbas, Farah Aini Abdullah, Mawahib Elamin, Mona Elmahi, M-shaped rational, homoclinic breather, kink-cross rational, multi-wave and interactional soliton solutions to the fifth-order Sawada-Kotera equation, 2024, 11, 26668181, 100919, 10.1016/j.padiff.2024.100919
  • Reader Comments
  • © 2023 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(1351) PDF downloads(62) Cited by(9)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog