Front. Phys. Frontiers in Physics Front. Phys. 2296-424X Frontiers Media S.A. 1114319 10.3389/fphy.2023.1114319 Physics Original Research An investigation of a closed-form solution for non-linear variable-order fractional evolution equations via the fractional Caputo derivative Ali et al. 10.3389/fphy.2023.1114319 Ali Umair 1 * Naeem Muhammad 2 Alahmadi Reham 3 Abdullah Farah Aini 4 Khan Muhammad Asim 5 Ganie Abdul Hamid 3 * 1 Department of Applied Mathematics and Statistics, Institute of Space Technology, Islamabad, Pakistan 2 Department of Mathematics of Applied Sciences, Umm-Al-Qura University, Makkah, Saudi Arabia 3 Department of Basic Sciences, College of Science and Theoretical Studies, Saudi Electronic University, Riyadh, Saudi Arabia 4 School of Mathematical Sciences, Universiti Sains Malaysia, Penang, Malaysia 5 Department of Mathematical Sciences, Faculty of Sciences, Universiti Teknologi Malaysia, Johor Bahru, Malaysia

Edited by: Gangwei Wang, Hebei University of Economics and Business, China

Reviewed by: Weipeng Hu, Northwestern Polytechnical University, China

Zhonglong Zhao, North University of China, China

 *Correspondence: Umair Ali, umairkhanmath@gmail.com; Abdul Hamid Ganie, a.ganie@seu.edu.sa

This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics

 14 02 2023 2023 11 1114319 02 12 2022 23 01 2023 Copyright © 2023 Ali, Naeem, Alahmadi, Abdullah, Khan and Ganie. 2023 Ali, Naeem, Alahmadi, Abdullah, Khan and Ganie

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Determining the non-linear traveling or soliton wave solutions for variable-order fractional evolution equations (VO-FEEs) is very challenging and important tasks in recent research fields. This study aims to discuss the non-linear space–time variable-order fractional shallow water wave equation that represents non-linear dispersive waves in the shallow water channel by using the Khater method in the Caputo fractional derivative (CFD) sense. The transformation equation can be used to get the non-linear integer-order ordinary differential equation (ODE) from the proposed equation. Also, new exact solutions as kink- and periodic-type solutions for non-linear space–time variable-order fractional shallow water wave equations were constructed. This confirms that the non-linear fractional variable-order evolution equations are natural and very attractive in mathematical physics.

 space-time variable-order fractional shallow water wave equation variable-order Caputo fractional derivative Khater method closed-form solution graphical representation 1 Introduction

Fractional calculus is a generalization of traditional integer-order integration and derivation actions onto non-integer order. The idea of fractional calculus is as old as classical calculus; it was discussed for the first time by Leibniz and L’Hospital in 1965. The fractional- and variable-order VO fractional models gained more attention because these models describe the physical phenomenon properly as compared to integer-order differential models. The non-linear FEEs define different phenomena in various areas, such as signal preparation, medication, biology, and organic framework [1, 2]. Many strategies have been produced to solve integer/fractional-order problems. Various fractional-order literature works directed that the memory and/or non-locality of the system may change with time, space, or other conditions. So, here our focus is on VO fractional differential models, which describe the physical models that vary with time or space or space–time. For example, Akgül et al. [3] solved the VO FPDE numerically and presented numerical experiments to confirm the efficiency and feasibility. Katsikadelis [4] developed a numerical method for linear and non-linear VO FPDEs in the Caputo sense. The resultant numerical values demonstrated the accuracy of the proposed method. Sahoo et al. [5] reviewed the VO operator definitions and properties. They discussed the new transfer function and investigated the model of a dynamic viscoelastic oscillator. Sing et al. [6] suggested an SEIR model that modeled the 2014–2015 outbreak of the Ebola virus in Africa. They discussed the system of VO FDEs and estimated its parameters for one or more variables. Semary et al. [7] approximated the solution of Liouville–Caputo VO FPDEs with 0 < α t ≤ 1 based on the Chebyshev function and discussed many linear and non-linear non-integer-order PDEs. Taghipour and Aminikhah [8] proposed the ADI numerical scheme for the fractional-order model and discussed the theoretical analysis. Other related studies can be seen in [9–16]. The effective analytical and closed-form solutions are studied in the recent literature. For example, Uddin et al. [17] considered the two important fractional-order models, namely, equal width and generalized equal width that describe the dispersive waves. They used the fractional derivative in the Riemann–Liouville sense and the G ′ G , 1 G expansion approach has been used, and they confirmed that the proposed approach is powerful, very convenient, and computationally efficient. Barman et al. [18] worked on a generalized Kudryashov method to provide a generic and inclusive closed-form solution. The proposed approach confirmed various shapes of waveform solutions such as kink-shaped, bell-shaped, singular, and flat in a 3D form. In another study, Barman et al. [19] proposed the same technique for Konopelchenko–Dubrovsky and Landau–Ginzburg–Higgs models. They obtained various varieties of analytical solutions for different parameters. The solutions are obtained in 2D and 3D forms, which demonstrated the efficiency and reliability of the proposed method. Roy et al. [20] solved the two significant types of models and implemented the new generalized G ′ / G expansion method. They constructed the solution in trigonometric, hyperbolic, and rational forms with different parameters. Kumar et al. [21] found out the exact solution for the higher-dimensional Fokas and breaking soliton models by the generalized exponential function method. The authors observed that the suggested method is effective and powerful. Ali et al. [22] investigated the exact solution for the VO fractional modified equal width equation based on the exp − ϕ ξ method. The fractional derivative is obtained in the Caputo sense, and the obtained exact solution is new and somewhat natural in mathematical physics. Akhtar et al. [23] constructed exact and traveling wave solutions for the Konopelchenko–Dubrovsky model and used two types of integration schemes. The resultant solutions are dark, single, anti-kink forms having a wide range of applications in applied sciences. Islam et al. [24] worked on analytical techniques and found the solution for the fractional-order foam drainage equation and SRLW equation. They used the G ′ / G expansion method and investigated the traveling wave solution for the proposed models. Mamun et al. [25] discussed the double G ′ G , 1 G expansion approach for the breaking soliton and the (1 + 1)-dimensional classical Boussinesq equations and obtained different soliton solutions, such as kink, multi-periodic, single soliton, and periodic wave solutions for different values of parameters. The comprehensive study can be found in [26–35].

The aforementioned cited literature reported that so far only numerical studies have been discussed for VO models and no attempt has been made to find the closed form for such types of VO-FEEs. The objective of this paper is to discuss the closed-form solution of the non-linear VO-FEEs. Here, we solve the non-linear VO fractional shallow water wave equation with CFD using the Khater method. The VO fractional problems are more complex computationally than a constant fractional order, and the evolution of a system can be furthermore clearly and accurately described. This contribution seems natural and simple and models many systems with VO [36]. The traveling wave solutions for the VO physical models are not known to the authors.

 2 The outline of the Khater method

The non-linear variable-order α x , y , … , t FPDE is given as H D x α x , y , . . , t Y , D x α x , y , . . , t D t α x , y , . . , t Y , D x α x , y , . . , t D y 2 α x , y , . . , t Y , … = 0 . where H is a polynomial for Y , Y t , Y x , D t α x , y , . . , t , D x α x , y , . . , t , D y α x , y , . . , t and D α x , y , . . , t represents Caputo fractional derivatives of the variable-order α x , y , . . . , t . The Caputo fractional derivative of the variable order for a function Y x , t of order γ ( x , t ) ϵ ( 0 , 1 is defined as follows [22]: D t γ x , t Y x , t = 1 Γ 1 + γ x , t ∫ 0 t Y ′ x , t Γ t − ξ γ x , t d ξ , 0 < γ x , t < 1 , Y ′ x , t , γ x , t = 1 .

Also, the important property is given as follows: D t γ x , t t β = Γ 1 − β Γ 1 − β + γ x , t t β − γ x , t , 0 < γ x , t < 1 .

Eq. 1 involved the linear and non-linear highest-order derivatives. A brief explanation of the proposed method is as follows [37]:

Convert the variable-order FPDE into an ordinary differential equation (ODE) by taking the transformation as Y x , y , t = y ξ , ξ = k x α x , y , . . , t Γ 1 + α x , y , . . , t + l y α x , y , . . , t Γ 1 + α x , y , . . , t + ⋯ − ω t α x , y , . . , t Γ 1 + α x , y , . . , t .

The obtained ODE is as follows: H y , ω y ′ , k y ″ , l y ‴ , ω l y ‴ , k l y ‴ , … = 0 , where k , l , m , and ω are constant parameters, if necessary, integrate Eq. 5. Next, we constructed a trial solution which can be expressed as y ξ = ∑ n = 0 M a n a n f ξ , where a n   n = 1 , 2 , … , M − 1 can be zero and a M ≠ 0 , and the function f ξ satisfies the following second-order linear equation: f ′ ξ = 1 ln a α a − f ξ + β + σ a f ξ .

The aforementioned equation has 27 possible solutions [33], which are derived by formulating various traveling wave solutions. Furthermore, the balancing principle is used to find M . Substituting Eq. 6 in Eq. 5 and Eq. 7, an equation involving the term a f ξ is obtained. In the obtained system of equations, the same power of a f ξ is equated to zero. The equations are solved simultaneously to find all unknown constants.

The solutions to Eq. 7:

When β 2 − 4 α σ < 0 and σ ≠ 0 , a f ξ = − β + 4 α σ − β 2 tan 4 α σ − β 2 2 ξ 2 σ

or a f ξ = − β − 4 α σ − β 2 cot 4 α σ − β 2 2 ξ 2 σ .

When β 2 − 4 α σ > 0 and σ ≠ 0 , a f ξ = − β − β 2 − 4 α σ tanh β 2 − 4 α σ 2 ξ 2 σ

or a f ξ = − β − β 2 − 4 α σ coth β 2 − 4 α σ 2 ξ 2 σ .

When β 2 + 4 α 2 < 0 , σ ≠ 0 , and σ = − p , a f ξ = β − − β 2 + 4 α 2 tan − β 2 + 4 α 2 2 ξ 2 α

or a f ξ = + β − β 2 + 4 α 2 cot − β 2 + 4 α 2 2 ξ 2 α .

When β 2 + 4 α 2 > 0 , σ ≠ 0 , and σ = − α , a f ξ = β + β 2 + 4 α 2 tanh q 2 + 4 p 2 2 ξ 2 α

or a f ξ = β + β 2 + 4 α 2 coth β 2 + 4 α 2 2 ξ 2 α .

When β 2 − 4   α 2 < 0 and σ = α , a f ξ = − β + − β 2 − 4 α 2 tanh − q 2 − 4 p 2 2 ξ 2 α

or a f ξ = − β − − β 2 − 4 α 2 coth − β 2 − 4 α 2 2 ξ 2 α .

When β 2 − 4 α 2 > 0 and σ = α , a f ξ = − β − β 2 − 4 α 2 tanh β 2 − 4 α 2 2 ξ 2 α

or a f ξ = − β − β 2 − 4 α 2 cot ⁡ h β 2 − 4 α 2 2 ξ 2 α .

When β 2 = 4 α σ , a f ξ = − 2 + β ξ ⁡ 2 σ ξ .

When σ α < 0 , β = 0 , and σ ≠ 0 , a f ξ = − − α σ tanh − σ α ξ

or a f ξ = − − α σ cot ⁡ h − σ α ξ .

When β = 0 and α = − σ , a f ξ = 1 + e − 2 σ ξ ⁡ − 1 + e − 2 σ ξ .

When α = σ = 0 , a f ξ = cosh β ξ + sinh β ξ .

When α = β = k and σ = 0 , a f ξ = e k ξ − 1 .

When β = σ = k and α = 0 , a f ξ = e k ξ 1 − e k ξ .

When β = α + σ , a f ξ = − 1 − α e α − σ ξ 1 − σ e α − σ ξ .

When β = − α + σ , a f ξ = α − e α − σ ξ σ − e α − σ ξ .

When α = 0 , a f ξ = β   e β ξ 1 − σ e β ξ .

When σ = β = α ≠ 0 , a f ξ = 1 2   3 tan 3 2 α ξ − 1 .

When σ = β = 0 , a f ξ = α   ξ .

When α = β = 0 , a f ξ = − 1 σ ξ .

When σ = α and β = 0 , a f ξ = tan α ξ .

When σ = 0 , a f ξ = e β ξ − α β .

The exact solutions for Eq. 1 are obtained by substituting unknown constants and Eq. 7 in Eq. 6.

 3 Formulation for the solutions of shallow water wave equations

Shallow water waves arise in the ocean when the waves move from the center of the ocean to the shore or beach known as shallow water waves. Most of the ocean waves are produced by wind, tsunamis, earthquakes, tides, etc. [38], which carry energy. Tsunamis and tides are both shallow water waves. The shallow water wave equation has been derived from the Navier–Stokes equations. Here, we apply the proposed method to study the non-linear space–time fractional VO shallow water wave equation and construct a traveling wave solution based on the Khater method.

3.1 The non-linear space–time variable-order fractional shallow water wave equation

We consider the space–time VO fractional shallow water wave equation as follows [39]: D t 2 γ x , t D x δ x , t Y + 3 D x δ x , t Y   D t δ x , t Y − D x δ x , t Y − D t γ x , t Y = 0 .

Using the wave variable ξ = k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t , Eq. 35 simplified to k ω 2 y ‴ − 3 k ω y ′ 2 − k y ′ + ω y ′ = 0 .

By balancing the highest-order non-linear term y ′ 2 and the highest-order linear term y ‴ , we obtain M = 1 . Therefore, the solution of Eq. 6 becomes y = a 0 + a 1 a f ξ .

Substituting Eq. 37 into Eq. 36 yields a polynomial equation for a f ξ . Equating the like powers of a f ξ n , we attain a system of algebraic equations given as a f ξ 0 : 2 α 2 k ω 2 σ a 1 2 − 3 α 2 k ω a 1 4 + α β 2 k ω 2 a 1 2 − α k a 1 2 + α ω a 1 2 = 0 , a f ξ 1 : 8 α β k ω 2 a 1 2 − 6 α β k ω a 1 4 + β 3 k ω 2 a 1 2 − β k a 1 2 + β ω a 1 2 = 0 , a f ξ 2 : 8 α k ω 2 σ 2 a 1 2 − 6 α k ω a 1 4 + 7 β 2 k ω a 1 4 − k σ a 1 2 + ω σ a 1 2 = 0 , a f ξ 3 : 12 β k ω 2 σ 2 a 1 2 − 6 α β k ω a 1 4 = 0 , a f ξ 4 : 6 k ω 2 σ 3 a 1 2 − 3 k ω σ 2 a 1 4 = 0 .

Solving the aforementioned system of algebraic equations by using computer algebra, we obtain S e t   1 : a 0 = a 0 , a 1 = 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ 4 α σ − β 2 k , ω = 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 2 4 α σ − β 2 k ,

where ω , k , α , β , and σ are arbitrary constants.

Substituting Eq. 38 into Eq. 37, we obtain y 1 = a 0 + 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ 4 α σ − β 2 k a f ξ .

Now, substituting the solutions of Eq. 7, we obtain the following 27 distinct traveling wave solutions for space–time fractional variable-order shallow water wave Eq. 35:

When β 2 − 4   α   σ < 0 and σ ≠ 0 , Y 1 = 4 α σ − β 2 tan 1 2 4 α σ − β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ 4 α σ − β 2 k + 2 a 0 σ 2 σ , Y 2 = − 4 α σ − β 2 cot 1 2 4 α σ − β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When β 2 − 4   α   σ > 0 and σ ≠ 0 , Y 3 = − − 4 α σ + β 2 tanh 1 2 − 4 α σ + β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ , Y 4 = − − 4 α σ + β 2 coth 1 2 − 4 α σ + β 2 k x δ x , y , t Γ 1 + δ x , y , t − ω t γ x , y , t Γ 1 + γ x , y , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When β 2 + 4   α 2 < 0 , σ ≠ 0 , and σ = − p , Y 5 = − − 4 α 2 − β 2 tan 1 2 − 4 α 2 − β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ , Y 6 = − − 4 α 2 − β 2 cot 1 2 − 4 α 2 − β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When β 2 + 4   α 2 > 0 , σ ≠ 0 , and σ = − p , Y 7 = 4 α 2 + β 2 tanh 1 2 4 α 2 + β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ , Y 8 = 4 α 2 + β 2 coth 1 2 4 α 2 + β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When β 2 − 4   α 2 < 0 and σ = α , Y 9 = 4 α 2 − β 2 tan 1 2 4 α 2 − β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ , Y 10 = − 4 α 2 − β 2 cot 1 2 4 α 2 − β 2 k x δ x , y , t Γ 1 + δ x , y , t − ω t γ x , y , t Γ 1 + γ x , y , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When q 2 − 4   p 2 > 0 and r = p , Y 11 = − − 4 α 2 + β 2 tanh 1 2 − 4 α 2 + β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ , Y 12 = − − 4 α 2 + β 2 coth 1 2 − 4 α 2 + β 2 k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − β − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + 2 a 0 σ 2 σ .

When q 2 = 4   p   r , Y 13 = − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k − β k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 2 + 2 a 0 σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t 2 σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t .

When σ   α < 0 , β = 0 , and σ ≠ 0 , Y 14 = − − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k − α σ tanh − α σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + a 0 , Y 15 = − − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k − α σ coth − α σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + a 0 .

When β = 0 and α = − σ , Y 16 = − e 2 σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k e 2 σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 + a 0 .

When α = σ = 0 , Y 17 = cos ⁡ h k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + sin ⁡ h k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + a 0 .

When α = β = k and σ = 0 , Y 18 = − e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + a 0 .

When β = σ = z and α = 0 , Y 19 = − − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 + a 0 .

When β = α + σ , Y 20 = a 0 − − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k .

When β = − α + σ , Y 21 = − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − e α k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t σ e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − e α k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + a 0 .

When α = 0 , Y 22 = − − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k β e β k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t σ e σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 + a 0 .

When σ = β = α ≠ 0 , Y 23 = a o + 1 2 − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k tan 3 2 α k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − 1 .

When σ = β = 0 , Y 24 = a o + − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k α k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t .

When p = q = 0 , Y 25 = a o σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t + − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k σ k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t .

When σ = α and β = 0 , Y 26 = a o + − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k tan α k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t .

When r = 0 , Y 27 = 1 β e β k x δ x , t Γ 1 + δ x , t − ω t γ x , t Γ 1 + γ x , t − α − 1 + 1 − 16 α k 2 σ + 4 β 2 k 2 σ − 4 α σ + β 2 k + a 0 .

4 Graphical representation

This section focuses on the graphical representation of some specific findings. Marwan and Aminah [40] solved the generalized shallow water equation by the (G′/G)-expansion and constructed a new exact solution for the proposed method. Bagchi et al. [41] extended the elliptic function method and found the traveling wave solution for the generalized shallow water wave equation. The obtained solutions are in the form of singular and periodic soliton solutions. Here, in this study, the graphical results obtained for different values of VO γ x , t and δ x , t are shown in Figures 1–5 for Eq. 35 in the form of 3D and 2D plots. Figure 1 and Figure 4 show the singleton soliton solution, and Figure 2, Figure 3, and Figure 5 represent the kink-shaped solution obtained using Maple 16 software.

Periodic solution for Eq. 35 for Y 12 at ω = 1 2 , k = 2 , a 0 = 1 , a 1 = − 1 2 , γ x , t = 2 ⁡ cos x t + x t 56 , δ x , t = 2 − sin x t 50 , α = 4 , β = 5 , σ = − 3 , x = 4.5 .

Kink-shaped solution for Eq. 35 for Y 5 at ω = 2 , k = 3 , a 0 = 5 , a 1 = 5 , γ x , t = 2 ⁡ cos x t + x t 56 , δ x , t = 2 − sin x t 50 , α = 4 , β = 5 , σ = − 30 , x = 4 .

Kink-shaped solution for Eq. 35 for Y 7 at ω = 1 , k = 1 , a 0 = 1 , a 1 = − 1 2 , γ x , t = 20 − x t 4 50 , δ x , t = 15 + sin 6 x t 50 , α = 4 , β = 5 , σ = − 3 , t = 1 .

Kink-shaped solution for Eq. 35 for Y 23 at ω = 1 , k = 1 , a 0 = 1 , a 1 = − 1 2 , γ x , t = 2 cos x t + x t 56 , δ x , t = 2 − sin x t 50 , α = 4 , β = 5 , σ = − 3 , e = 2.7128 , x = 3.5 .

Kink-shaped solution for Eq. 35 for Y 24 at ω = 1 , k = 1 , a 0 = 1 , a 1 = − 1 2 , γ x , t = 2 ⁡ cos x t + x t 56 , δ x , t = 2 − sin x t 50 , α = 4 , β = 5 , σ = − 3 , e = 2.7128 , z = 1 .

 5 Conclusion

In this paper, we solved the non-linear VO fractional evolution equation successfully in the Caputo fractional derivative sense and obtained new exact traveling wave solutions. The VO fractional evolution equation is discussed quite efficiently and accurately by using the Khater method. Here, 27 exact solutions having Kink and singular soliton-type solutions are obtained for different values of VO γ x , t and δ x , t for the proposed Caputo fractional VO equation. The different values of parameters examine different physical phenomena. This contribution is effective, instrumental, and evangelistic and seems more natural in the literature. This study can be extended to other types of VO FPDEs and can be solved by various analytical techniques.

 Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

 Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

 Funding

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work under grant code: 22UQU4310396DSR54.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Abbreviations

ADI, alternating direction implicit scheme; EW, equal width; FPDE, fractional partial differential equation; GEW, generalized equal width; ODE, ordinary differential equation; SRLW, symmetric regularized long wave; VO, variable order; 2D, two-dimensional; 3D, three-dimensional.

 References Bekir A Güner Ö . The G′G-expansion method using modified Riemann–Liouville derivative for some space-time fractional differential equations. Ain Shams Eng J (2014) 5(3):959–65. 10.1016/j.asej.2014.03.006 Bin Z . (G′/G)-expansion method for solving fractional partial differential equations in the theory of mathematical physics. Communic Theor Phys (2012) 58:623–30. 10.1088/0253-6102/58/5/02 Akgül A Baleanu D . On solutions of variable-order fractional differential equations. Int J Optimization Control Theories Appl (Ijocta) (2017) 7(1):112–6. 10.11121/ijocta.01.2017.00368 Katsikadelis JT . Numerical solution of variable order fractional differential equations (2018). Available at: https://arXiv.org/abs/preprint.arXiv:1802.00519 . Sahoo S Saha Ray S Das S Bera RK . The formation of dynamic variable order fractional differential equation. Int J Mod Phys C (2016) 27(07):1650074. 10.1142/s0129183116500741 Singh AK Mehra M Gulyani S . Learning parameters of a system of variable order fractional differential equations. Numer Methods Partial Differential Equations (2021) num.22796. 10.1002/num.22796 Semary MS Hassan HN Radwan AG . The minimax approach for a class of variable order fractional differential equation. Math Methods Appl Sci (2019) 42(8):2734–45. 10.1002/mma.5545 Taghipour M Aminikhah H . A new compact alternating direction implicit method for solving two-dimensional time fractional diffusion equation with Caputo-Fabrizio derivative. Filomat (2020) 34(11):3609–26. 10.2298/fil2011609t Ali U Sohail M Usman M Abdullah FA Khan I Nisar KS . Fourth-order difference approximation for time-fractional modified sub-diffusion equation. Symmetry (2020) 12(5):691. 10.3390/sym12050691 Ali U Kamal R Mohyud-Din ST . On nonlinear fractional differential equations. Int J Mod Math Sci (2012) 3:3., no. Zhao Z Yue J He L . New type of multiple lump and rogue wave solutions of the (2+ 1)-dimensional Bogoyavlenskii–Kadomtsev–Petviashvili equation. Appl Math Lett (2022) 133:108294. 10.1016/j.aml.2022.108294 Zubair T Usman M Ali U Mohyud-Din ST . Homotopy analysis method fa or system of partial differential equations. Int J Mod Eng Sci (2012) 1(2):67–79. Zhao Z He L . A new type of multiple-lump and interaction solution of the Kadomtsev–Petviashvili I equation. Nonlinear Dyn (2022) 109:1033–46. 10.1007/s11071-022-07484-6 Hu W Deng Z Han S Zhang W . Generalized multi-symplectic integrators for a class of Hamiltonian nonlinear wave PDEs. J Comput Phys (2013) 235:394–406. 10.1016/j.jcp.2012.10.032 Hu W Xu M Song J Gao Q Deng Z . Coupling dynamic behaviors of flexible stretching hub-beam system. Mech Syst Signal Process (2021) 151:107389. 10.1016/j.ymssp.2020.107389 Zhao Z He L . Lie symmetry, nonlocal symmetry analysis, and interaction of solutions of a (2+ 1)-dimensional KdV–mKdV equation. Theor Math Phys (2021) 206(2):142–62. 10.1134/s0040577921020033 Uddin MH Khatun MA Arefin MA Akbar MA . Abundant new exact solutions to the fractional nonlinear evolution equation via Riemann-Liouville derivative. Alexandria Eng J (2021) 60(6):5183–91. 10.1016/j.aej.2021.04.060 Barman HK Seadawy AR Roy R Akbar MA Raddadi MH . Rational closed form soliton solutions to certain nonlinear evolution equations ascend in mathematical physics. Results Phys (2021) 27:104450. 10.1016/j.rinp.2021.104450 Barman HK Akbar MA Osman MS Nisar KS Zakarya M Abdel-Aty AH Solutions to the Konopelchenko-Dubrovsky equation and the Landau-Ginzburg-Higgs equation via the generalized Kudryashov technique. Results Phys (2021) 24:104092. 10.1016/j.rinp.2021.104092 Roy R Akbar MA Seadawy AR Baleanu D . Search for adequate closed form wave solutions to space–time fractional nonlinear equations. Partial Differential Equations Appl Math (2021) 3:100025. 10.1016/j.padiff.2021.100025 Kumar S Niwas M Osman MS Abdou MA . Abundant different types of exact soliton solution to the (4+1)-dimensional Fokas and (2+1)-dimensional breaking soliton equations. Commun Theor Phys (2021) 73:105007. 10.1088/1572-9494/ac11ee Ali U Mastoi S Othman WAM Khater MM Sohail M . Computation of traveling wave solution for nonlinear variable-order fractional model of modified equal width equation. AIMS Math (2021) 6(9):10055–69. 10.3934/math.2021584 Akhtar J Seadawy AR Tariq KU Baleanu D . On some novel exact solutions to the time fractional (2+ 1) dimensional Konopelchenko–Dubrovsky system arising in physical science. Open Phys (2020) 18(1):806–19. 10.1515/phys-2020-0188 Islam T Akbar MA Azad AK . Traveling wave solutions to some nonlinear fractional partial differential equations through the rational (G′/G)-expansion method. J Ocean Eng Sci (2018) 3(1):76–81. 10.1016/j.joes.2017.12.003 Mamun Miah M Shahadat Ali HM Ali Akbar M Majid Wazwaz A . Some applications of the (G′/G, 1/G)-expansion method to find new exact solutions of NLEEs. The Eur Phys J Plus (2017) 132(6):252–15. 10.1140/epjp/i2017-11571-0 Islam MS Khan K Akbar MA Mastroberardino A . A note on improved F-expansion method combined with Riccati equation applied to nonlinear evolution equations. R Soc Open Sci (2014) 1(2):140038. 10.1098/rsos.140038 Hu W Huai Y Xu M Feng X Jiang R Zheng Y Mechanoelectrical flexible hub-beam model of ionic-type solvent-free nanofluids. Mech Syst Signal Process (2021) 159:107833. 10.1016/j.ymssp.2021.107833 Hu W Xu M Zhang F Xiao C Deng Z . Dynamic analysis on flexible hub-beam with step-variable cross-section. Mech Syst Signal Process (2022) 180:109423. 10.1016/j.ymssp.2022.109423 Ali U Ahmad H Abu-Zinadah H . Soliton solutions for nonlinear variable-order fractional Korteweg–de Vries (KdV) equation arising in shallow water waves. J Ocean Eng Sci (2022). 10.1016/j.joes.2022.06.011 Gulalai SA Rihan FA Ullah A Al-Mdallal QM Akgul A . Nonlinear analysis of a nonlinear modified KdV equation under Atangana Baleanu Caputo derivative. AIMS Math (2022) 7(5):7847–65. 10.3934/math.2022439 Hu W Wang Z Zhao Y Deng Z . Symmetry breaking of infinite-dimensional dynamic system. Appl Math Lett (2020) 103:106207. 10.1016/j.aml.2019.106207 Zayed Elsayed ME Amer Yasser A Shohib Reham MA . Exact traveling wave solution for nonlinear fractional partial differential equation using the improved (G’/G)-expansion method. Int J Engin (2014) 4(7):18–31. Ali U Ganie AH Khan I Alotaibi F Kamran K Muhammad S Traveling wave solutions to a mathematical model of fractional order (2+ 1)-dimensional breaking soliton equation. Fractals (2022). p. 2240124. Hu W Zhang C Deng Z . Vibration and elastic wave propagation in spatial flexible damping panel attached to four special springs. Commun Nonlinear Sci Numer Simulation (2020) 84:105199. 10.1016/j.cnsns.2020.105199 Ali U Ahmad H Baili J Botmart T Aldahlan MA . Exact analytical wave solutions for space-time variable-order fractional modified equal width equation. Results Phys (2022) 33:105216. 10.1016/j.rinp.2022.105216 Hu W Huai Y Xu M Deng Z . Coupling dynamic characteristics of simplified model for tethered satellite system. Acta Mechanica Sinica (2021) 37(8):1245–54. 10.1007/s10409-021-01108-9 Nayfeh AH Balachandran B . Applied nonlinear dynamics: Analytical, computational, and experimental methods. John Wiley & Sons (2008). Pushkarev AN Zakharov VEE . Nonlinear amplification of ocean waves in straits. Theor Math Phys (2020) 203(1):535–546. Rezazadeh H Korkmaz A Khater MM Eslami M Lu D Attia RA . New exact traveling wave solutions of biological population model via the extended rational sinh-cosh method and the modified Khater method. Mod Phys Lett B (2019) 33(28):1950338. 10.1142/s021798491950338x Clarkson PA Mansfield EL . On a shallow water wave equation. Nonlinearity (1994) 7(3):975–1000. 10.1088/0951-7715/7/3/012 Bagchi B Das S Ganguly A . New exact solutions of a generalized shallow water wave equation. Physica Scripta (2010) 82(2):025003. 10.1088/0031-8949/82/02/025003