Water and Gender in Recreating Family Life with Maa Ganga: The Confluence of Nature and Culture in a North Indian River Pilgrimage

1.   Introduction

The investigation of exact solutions to nonlinear fractional differential equations (NLFDEs) plays an active role in nonlinear sciences in the perspective of the physical phenomena of real world can effectively be modeled by the making use of the theory of derivatives and integrals of fractional order. NLFDEs have recently proved to be valuable tools to the modeling of many physical phenomena and have expanded the attention of many studies due to their frequent entrance in various applications, such as optical fibers, plasma physics, solid state physics, control theory, chemical kinematics, signal processing, fractional dynamics, fluid flow and other areas. The NLFDEs are also used in modeling of many chemical progressions, mathematical biology and many other problems in physics and engineering. In the last two decades, numerous techniques have been suggested to examine the exact solution to NLFDEs. As a consequence, a variety of methods have been established to get the exact solution to NLFDEs such as the variational iteration method ^[1,2,3,4,5], the Adomian's decomposition method ^[6,7,8], the differential transformation method ^[9,10], the Laplace perturbation method ^[11], the homotopy perturbation method ^[12,13,14] the finite element method ^[15], the $(G'/G)$-expansion method ^{[16,17,18,19]}, the improve Bernoulli sub-equation method ^[20,21,22], the Hirota bilinear method ^[23], the extended sinh-Gordon expansion method ^[24], the exp-function method ^[25,26,27], the fractional sub-equation method ^[28,29,30], the first integral method ^[31,32], the two variable $(G'/G, 1/G)$-expansion method ^[33,34,35] and so forth ^[36,37,38].

Alzaidy ^[39] derived the exact and analytical solution to space-time fractional mBBM and space-time fractional ZKBBM equation by using fractional sub-equation method. Ege et al. ^[40] established the analytical and approximate solution to the space time fractional mBBm and space time fractional potential Kadomtsev-Petviashvili equation by the modified Kudryashov method. Recently, Bekir et al. ^[41] studied the exact solution to the space time fractional mKdV and space time fractional mBBM equation by the first integral method. On the other hand, Song et al. ^[42] investigated the exact traveling wave solution to the space time fractional ZKBBM equation by Bifurcation method. Aksoy et al. ^[43] derived the analytical solution to the time fractional Fokas equation, fractional ZKBBM equation, and fractional couple Burgers equation by means of the exponential rational function method. In addition, Ekici et al. ^[44] studied the solitary wave solution and periodic wave solution to the space-time fractional ZKBBM by a fractional sub-equation method. To our comprehension, the space-time fractional mBBM and fractional ZKBBM equation have not been studied by means of recently established two variables $(G'/G, 1/G)$-expansion method. So, the goal of this study is to establish new solutions to the above-mentioned equation by means of the proposed method. It is noteworthy to observe that, the two variable $(G'/G, 1/G)$-expansion method is the generalization of $(G'/G)$-expansion method. The main idea of this method is that the solutions to NLFDEs are represented as a polynomial in two variables $(G'/G)$ and $\left({1/G} \right)$, wherein $G = G\left(\xi \right)$ satisfies the second order ODE $G''\left(\xi \right) + \lambda G\left(\xi \right) = \mu$, where $G''\left(\xi \right) + \lambda G\left(\xi \right) = \mu$ and $\mu$ are constants. The objective of this article is to establish further general and some fresh close form solitary wave solution to the space-time fractional mBBM equation and the space-time fractional ZKBBM equation by means of the two variable $(G'/G, 1/G)$-expansion method. The results obtained in this article have been associated with those existing in the literature and shown that the achieved solutions are new and further general. We have to end with determined that the solution might bring up their prominence throughout the contribution and be recorded in the literature.

The rest of this article is organized as follows: In section 2, we have introduced preliminaries and fundamental tools, in section 3, we have described the two variable $(G'/G, 1/G)$-expansion method, in section 4, we have established the exact solution to the space-time fractional mBBM equation and the space time fractional ZKBBM equation by the proposed method. In section 5, we investigate results and discussion and in the last section, the conclusions are drawn.

2.   Preliminaries and fundamental tools

Let $f:\left[ {0, \infty } \right) \to \mathbb{R}$, be a function. The $\alpha$-order "conformable derivative" of $f$ is defined by ^[45]:

for all t > 0, α∈(0, 1). If f is α-differentiable in some (0, a), a > 0, and ${\lim _{t \to 0}} + {f^{\left(\alpha \right)}}\left(t \right)$ exists, then define ${f^{\left(\alpha \right)}}\left(0 \right) = {\lim _{t \to 0}} + {f^{\left(\alpha \right)}}\left(t \right)$. The following theorems point out some properties which are satisfied with the conformable derivative.

Theorem 1. Let $\alpha \in \left({0, 1} \right]$ and suppose $f, g$ be $\alpha$-differentiable at a point $t > 0$. Then

• ${T_\alpha }\left({cf + dg} \right) = c{T_\alpha }\left(f \right) + d{T_\alpha }\left(g \right)$, for all $c, d \in \mathbb{R}$.

• ${T_\alpha }\left({{t^p}} \right) = p{t^{p - \alpha }}$, for all $p \in \mathbb{R}$.

• ${T_\alpha }\left(c \right) = 0$, for all constant function $f\left(t \right) = c$.

• ${T_\alpha }\left({fg} \right) = f{T_\alpha }\left(g \right) + g{T_\alpha }\left(f \right)$.

• ${T_\alpha }\left({\frac{f}{g}} \right) = \frac{{g{T_\alpha }\left(f \right) - f{T_\alpha }\left(g \right)}}{{{g^2}}}$.

• In addition, if $f$ is differentiable, then ${T_\alpha }\left(f \right)\left(t \right) = {t^{1 - \alpha }}\frac{{df}}{{dt}}$.

Some more properties involving the chain rule, Gronwall's inequality, some integration techniques, Laplace transform, Tailor series expansion and exponential function with respect to the conformable fractional derivative are explained in the study ^[46].

Theorem 2. Let $f$ be an $\alpha$-differentiable function in conformable differentiable and suppose that $g$ is also differentiable and defined in the range of $f$. Then

3.   The two variable $(G'/G, 1/G)$-expansion method

Suppose the second order differential equation

and let us consider the following relations

Thus, it is obtained

The solution to the mentioned equation in (3.1) depends on $\lambda$ for which its values as $\lambda < 0, \lambda > 0$ and $\lambda = 0$.

Remark 1: If $\lambda < 0$, the general solution to equation (3.1) is:

where ${A_1}$ and ${A_2}$ are arbitrary constants. Accordingly, it yields

where $\sigma = A_1^2 - A_2^2$.

Remark 2: If $\lambda > 0$, the general solution to equation (3.1) is:

where ${A_1}$ and ${A_2}$ are arbitrary constants. Consequently, it provides

where $\sigma = A_1^2 + A_2^2$.

Remark 3: If $\lambda = 0$, the general solution to equation (3.1) is:

where ${A_1}$ and ${A_2}$ are arbitrary constants. As a result, it delivers

Let us suppose the general nonlinear fractional differential equation (NLFDE) of the form

where $u$ is an unidentified function of spatial derivative x and temporal derivative t and p is a polynomial of u and its partial fractional derivatives wherein the maximum number of derivatives and the nonlinear terms are involved.

Step 1: We suppose the traveling wave transformation in the subsequent

where c and k are nonzero arbitrary constants.

Using this wave transformation, we can rewrite the equation (3.10) as:

where prime stands for the derivative of u with respect to $\xi$.

Step 2: Assume that the solution of (3.12) can be revealed as a polynomial in two variables $\phi$ and $\psi$ as the following form:

where in ${a_i}\left({i = 0, 1, 2, \ldots } \right)$ and ${b_i}\left({i = 0, 1, 2, \ldots } \right)$ are constants to be evaluated afterword.

Step 3: The homogeneous balancing between the highest order derivatives and the nonlinear terms appearing in equation (3.12) determined the positive integer N.

Step 4: Setting (3.13) into (3.12) together with (3.3) and (3.5), the equation (3.12) can be reduced into a polynomial in $\phi$ and $\psi$, where the degree of $\psi$ is not greater than one. Equating the coefficients of this polynomial of like power to zero give a system of algebraic equations which can be solved by using the computer algebra, like Maple or Mathematica yields the values of ${a_i}, {b_i}, \mu, {A_1}, {A_2}$ and $\lambda$ where $\lambda < 0$, which provide hyperbolic function solutions.

Step 5: Similarly, we examine the values of ${a_i}, {b_i}, \mu, {A_1}, {A_2}$ and $\lambda$ when $\lambda > 0$ and $\lambda = 0$, yield the trigonometric and rational function solutions respectively.

4.   Formulation of exact solutions

4.1.   The exact wave solutions to the space-time fractional mBBM equation

In this subsection, we investigate more general and new closed form wave solutions to the space time fractional mBBM equation by means of the two variables $(G'/G, 1/G)$-expansion method. We introduce the space time fractional mBBM equation ^[39] as follows:

where v is a nonzero constant. If spatial derivative term $D_x^\alpha u$ is removed then the equation (4.1) transformed into mKdV equation. This equation was first derived to designate an approximation for surface long waves in nonlinear dispersive media. It is also describe the acoustic waves in inharmonic crystals, the hydro-magnetic waves in cold plasma and acoustic gravity waves in compressible fluids.

For the mBBM equation (4.1), we present the following transformation:

where c is the speed of the traveling wave. Using the transformation (4.2), the equation (4.1) reduces to the following integer order ordinary differential equation (ODE):

Integrating equation (4.3) with zero constant (see ^[47] for details), we attain

Balancing the highest order derivative ${u''}$ and the nonlinear term ${{u^3}}$, yields $N = 1$. So, the solution to (4.3) converted to the following form:

where ${a_0}, {a_1}$ and ${b_1}$ are constants to be determined.

Case 1: For $\lambda < 0$, setting equation (4.5) into (4.4) along with (3.3) and (3.5) yields a set of algebraic equations and solving these equations by using symbolic computation software Maple, we obtained the following results:

Substituting these values into (4.5), we find the solution to the mBBM equation as the form:

where $\xi = \pm {b_1}\sqrt { - \frac{{2v\lambda }}{{3{\mu ^2} + 3{\lambda ^2}\sigma }}} \frac{{{x^\alpha }}}{\alpha } \mp \frac{1}{3}\frac{{\sqrt { - \frac{{2v\lambda }}{{3{\mu ^2} + 3{\lambda ^2}\sigma }}} \times {b_1}\left({3{\mu ^2} - vb_1^2{\lambda ^2} + 3{\lambda ^2}\sigma } \right)}}{{{\mu ^2} + {\lambda ^2}\sigma }}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 - A_2^2$.

Since ${{A_1}}$ and ${{A_2}}$ are integral constants, one may choose arbitrarily their values. Therefore, if we choose ${A_1} = 0$, ${A_2} \ne 0$ and $\mu = 0$ in (4.6), we obtain the following solitary wave solution

Again if ${A_1} \ne 0$, ${A_2} = 0$ and $\mu = 0$ in (4.6), we obtain the solitary wave solution.

Case 2: In a similar manner, substituting (4.5) into (4.4) along side with (3.3) and (3.7) yield a set of algebraic equations for ${a_0}, {a_1}$, c and k and solving these equations, we obtain the following solutions:

When we substituted these values into equation (4.5), we found the following solution to mBBM equation:

where $\xi = {b_1}\sqrt {\frac{{2v\lambda }}{{3{\lambda ^2}\sigma - 3{\mu ^2}}}} \frac{{{x^\alpha }}}{\alpha } \pm \frac{1}{3}\frac{{\sqrt {\frac{{2v\lambda }}{{3{\lambda ^2}\sigma - 3{\mu ^2}}}} \times {b_1}\left({3{\mu ^2} - vb_1^2{\lambda ^2} - 3{\lambda ^2}\sigma } \right)}}{{{\mu ^2} - {\lambda ^2}\sigma }}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 + A_2^2$.

If we set ${A_1} = 0$, ${A_2} \ne 0$ (or ${A_1} \ne 0$, ${A_2} = 0$) and $\mu = 0$ in (4.9), we obtain the solitary wave solution

Case 3: In the similar fashion when $\lambda = 0$, using(4.5) and (4.4) along with(3.3) and (3.9), we will find a group of algebraic equations, whose solutions are as follows:

Substituting these values into equation (4.5), we attained the solution to equation (4.3)

It is remarkable to see that the traveling wave solution ${u_{{1_1} - {u_{{1_7}}}}}$ of the space-time fractional mBBM equation are new and further general and have not been investigated in the previous study. These solutions forces are convenient to characterize the hydromagnetic waves in cold plasma, acoustic waves in inharmonic crystals and acoustic-gravity waves incompressible fluids.

4.2.   The exact wave solutions to the space-time fractional ZKBBM equation

In this subsection, we establish the general and some fresh solutions to the space-time fractional ZKBBM equation over the two variable $(G'/G, 1/G)$-expansion method. This equation is significant in many physical phenomena and it arises as an explanation of gravity water waves in the long-wave regime. We consider the space-time fractional ZKBBM equation ^[39] as follows:

where a and b are arbitrary constants. If we removed the spatial derivative term $D_x^\alpha u$, the equation (4.13) converted to BBM equation.

For the space-time fractional ZKBBM equation, we present the following transformation:

where c is the velocity of the traveling wave. Using traveling wave variable (4.14) equation (4.13) reduces to the following ODE for $u = u\left(\xi \right)$:

Integrating (4.15) once with respect to $\xi$ with zero constant, we obtain

Balancing the highest order derivatives and nonlinear terms, we obtained $N = 2$. Therefore, the solution to equation (4.16) is of the following form:

Case 1: For $\lambda < 0$, inserting equation (4.17) into (4.16) along with (3.3) and (3.5) yields a set of algebraic equations whose solutions from the symbolic computation software Maple are as follows:

where k is an arbitrary constant.

For Set 1, we obtain the following general solution to the space time fractional ZKBBM equation:

where $\xi = k\frac{{{x^\alpha }}}{\alpha } - \frac{k}{{b{k^2}\lambda + 1}}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 - A_2^2$.

Here ${{A_1}}$ and ${{A_2}}$ are constants of integration. Therefore, one can freely select their values. If we choose ${A_1} = 0$, ${A_2} \ne 0$ and $\mu = 0$ in (4.18), we attain the solitary wave solution

Again if we set ${A_1} \ne 0, {A_2} = 0$ and $\mu = 0$ in (4.18), we obtain the solitary wave solution

Again for Set 2, we obtain the following general solution to the space time fractional ZKBBM equation:

where $\xi = k\frac{{{x^\alpha }}}{\alpha } + \frac{k}{{b{k^2}\lambda + 1}}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 - A_2^2$.

If we choose ${A_1} = 0$, ${A_2} \ne 0$ (or ${A_1} \ne 0$, ${A_1} \ne 0$) and $\mu = 0$ in (4.21), we obtain the solitary wave solution

Case 2: In the similar manner, when $\lambda > 0$ inserting (4.17) into (4.16) along with (3.3) and (3.7) yields a set of algebraic equation for ${a_0}, {a_1}, {b_1}, {a_2}, {b_2}, c, k$ and solving these equations we have achieved the following results:

Now, we find the following exact solution to the space time fractional ZKBBM equation (4.13) for Set 1:

where $\xi = k\frac{{{x^\alpha }}}{\alpha } + \frac{k}{{b{k^2}\lambda + 1}}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 + A_2^2$.

As ${A_1}$ and ${A_2}$ are integral constants, if we select ${A_1} = 0$, ${A_2} \ne 0$(or ${A_1} \ne 0$, ${A_2} = 0$) and $\mu = 0$ into (4.24), we attain the solitary wave solution

Again, we find the following exact solution to the space time fractional ZKBBM equation for Set 2:

where $\xi = k\frac{{{x^\alpha }}}{\alpha } - \frac{k}{{b{k^2}\lambda + 1}}\frac{{{t^\alpha }}}{\alpha }$ and $\sigma = A_1^2 + A_2^2$.

If we set ${A_1} = 0$, ${A_2} \ne 0$ (or ${A_1} \ne 0$, ${A_2} = 0$) and $\mu = 0$ into (4.27), we find the solitary wave solution:

Case 3: Finally when $\lambda = 0$, substituting equation (4.17) into (4.16) along with (3.3) and (3.9) yields a set of algebraic equations for ${a_0}, {a_1}, {b_1}, {a_2}, {b_2}, c, k$ and whose solution are as follows:

Inserting these values into equation (4.16), we attain the rational function solution to the space time fractional ZKBBM equation (4.13) as follows:

5.   Results and discussion

It is remarkable to observe that some of the obtained solutions demonstrate good similarity with earlier established solutions. A comparison between Mohyud-Din et al. ^[48] solutions and our obtained solutions is presented in the following Table 1:

Table 1. Comparison between Mohyud-Din et al. ^[48] solutions and the obtained solutions.

The trigonometric function solutions referred to the above table is similar and if we set definite values of the arbitrary constants they are identical. It is substantial to understand that the traveling wave solutions ${u_{{2_2}}}\left({x, t} \right)$, ${u_{{2_3}}}\left({x, t} \right)$, ${u_{{2_5}}}\left({x, t} \right)$, ${u_{{2_6}}}\left({x, t} \right)$, and ${u_{{2_2}}}\left({x, t} \right)$ of the fractional ZKBBM equation are all new and very much important which were not originate in the previous work. This diffusion equation is significant in various physical phenomena. Itarises as an explanation of gravity water waves in the long-wave regime that creates outstanding model in physics and engineering.

6.   Conclusion

In this study, we have obtained some new and further general solitary wave solutions to two nonlinear space time fractional differential equation, namely, the time fractional mBBM equation and the fractional ZKBBM equation in terms of hyperbolic, trigonometric and rational function solution involving parameters. It is remarkable to see that our achieve solutions through the suggested method are more new and further general compared to the existing literature. The obtained solutions to these equations are capable to investigate the mathematical model of gravity water waves in the long-wave regime, the acoustic waves in inharmonic crystals, hydromagnetic waves in cold plasma and acoustic-gravity waves incompressible fluids. The competence of the two variables $(G'/G, 1/G)$-expansion method is consistent, reliable and very much attractive. Since each nonlinear equation has its own inconsistent characteristic, the future research might be how to recommend method is compatible for revealing the solutions to other NLFDEs.

Acknowledgement

The authors would like to express their deepest appreciation to the reviewers for their valuable suggestions and comments to improve the article.

Conflict of interest

The authors declare no conflict of interest.