A new modified conformal expansion method to some nonlinear conformal partial differential equations

Alhakim, Lama; Moussa, Alaaeddin; Gasmi, Boubekeur; Mati, Yazid; Raihen, Md Nurul

doi:10.2478/ijmce-2025-0029

Full Article

1

Introduction

Various constraints encountered in many real-world phenomena can only be modeled by using nonlinear partial differential (NPD) equations. Some phenomena arising in ocean sciences include modeling short-distance seismic activities such as earthquakes and acoustic waves in nuclear physics. Other equations are used to model phenomena in physical science [1], finance [2], and optical science [3]. As a result, obtaining analytical solutions to nonlinear conformable partial differential (NCPD) equations allows for a more thorough examination of the phenomena governed by these equations.

Researchers have extensively studied the NPD equations to obtain exact and explicit soliton solutions using various approaches. Among the notable approaches are as follows: the reductive perturbation technique [4], the tanh(η) expansion method [5], the $(\frac{G^{'}}{G})$ \left({{{G'} \over G}} \right) expansion approach [6], the new extended auxiliary equation method [7], the improved Exp-function method [8], the double auxiliary equations approach [9, 10], the new generalized of e^(−Φ(ξ))-expansion approach [11], the coth_a (ξ) expansion approach [12], the generalized double auxiliary equation method [13,14,15], the F–expansion approach [16], the improved Jacobi elliptic function [17], the Cham method [18], the improved Cham method [19], the conformal fractional $(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})$ \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) -expansion approach [20], the Exp-function approach [21], the auxiliary equation approach [22], the extended auxiliary equation approach [23], the fractional sub-equation approach [24], the fractional Riccati equation rational expansion approach [25], the improved shooting method [26], and others [27,28,29,30].

This paper investigates the exact solutions to three different space-time nonlinear models with conformable such as the Whitham-Broer-Kaup equations [31] given as (1) $\begin{array}{l} D_{t}^{α} u + u D_{x}^{α} u + D_{x}^{α} v + β D_{x}^{2 α} u = 0, \\ D_{t}^{α} v + D_{x}^{α} (uv) - β D_{x}^{2 α} v + γ D_{x}^{3 α} u = 0, \end{array}$ \left\{{\matrix{{D_t^\alpha u + uD_x^\alpha u + D_x^\alpha v + \beta D_x^{2\alpha}u = 0,} \hfill \cr {D_t^\alpha v + D_x^\alpha \left({uv} \right) - \beta D_x^{2\alpha}v + \gamma D_x^{3\alpha}u = 0,} \hfill \cr}} \right. the Burgers equations [32] defined by (2) $\begin{array}{l} D_{t}^{α} h - D_{x}^{2 α} h + 2 h D_{x}^{α} h + p D_{x}^{α} (hg) = 0, \\ D_{t}^{α} g - D_{x}^{2 α} g + 2 g D_{x}^{α} g + q D_{x}^{α} (hg) = 0, \end{array}$ \left\{{\matrix{{D_t^\alpha h - D_x^{2\alpha}h + 2hD_x^\alpha h + pD_x^\alpha \left({hg} \right) = 0,} \hfill \cr {D_t^\alpha g - D_x^{2\alpha}g + 2gD_x^\alpha g + qD_x^\alpha \left({hg} \right) = 0,} \hfill \cr}} \right. and finally, the modified Korteweg-de Vries equations [33] read as (3) $\begin{array}{l} D_{t}^{α} ψ = \frac{1}{2} D_{x}^{3 α} ψ - 3 ψ^{2} D_{x}^{α} ψ + \frac{3}{2} D_{x}^{2 α} ψ + 3 D_{x}^{α} (ψ ϕ) - 3 δ D_{x}^{α} ψ, \\ D_{t}^{α} ϕ = - D_{x}^{3 α} ϕ - 3 ϕ D_{x}^{α} ϕ - 3 D_{x}^{α} ψ D_{x}^{α} ϕ + 3 ψ^{2} D_{x}^{α} ϕ + 3 δ D_{x}^{α} ϕ, \end{array}$ \left\{{\matrix{{D_t^\alpha \psi = {1 \over 2}D_x^{3\alpha}\psi - 3{\psi^2}D_x^\alpha \psi + {3 \over 2}D_x^{2\alpha}\psi + 3D_x^\alpha \left({\psi \phi} \right) - 3\delta D_x^\alpha \psi,} \hfill \cr {D_t^\alpha \phi = - D_x^{3\alpha}\phi - 3\phi D_x^\alpha \phi - 3D_x^\alpha \psi D_x^\alpha \phi + 3{\psi^2}D_x^\alpha \phi + 3\delta D_x^\alpha \phi,} \hfill \cr}} \right. where α ∈ (0, 1]; δ, p, and q are real constants.

The paper is divided into four sections. Section 2 provides an overview of Katugampola's conformal derivative and explains the proposed modified conformal expansion method (MCEM). Section 3 demonstrates the efficiency of our method on the three above mentioned equations and provides some physical interpretations of the results. Section 4 offers a conclusive overview of the work.

2

The proposed MCEM

The preliminaries and definitions of the proposed approach are given as follows.

2.1

The Katugampola's conformal derivative

Among the various derivative definitions, we use Katugampola's derivative [34] that fulfills the linearity property and the quotient, the power, the Leibniz, and the chain rules. For a given real function g(t) where t > 0, we can define the derivative of a given order α ∈ (0, 1) by: $D^{α} (g) (t) = lim_{ε \to 0} \frac{g (t e^{ε t^{- α}}) - g (t)}{ε}, \forall t > 0 .$ {D^\alpha}\left(g \right)\left(t \right) = \mathop {{\rm{lim}}}\limits_{\varepsilon \to 0} {{g\left({t{e^{\varepsilon {t^{- \alpha}}}}} \right) - g\left(t \right)} \over \varepsilon},\forall t > 0.

If the lim_t→0⁺ D^α (g)(t) exists and the function g is α–differentiable in (0, a), a > 0, the following value can be defined: (4) $D^{α} (g) (0) = lim_{t \to 0^{+}} D^{α} (g) (t) .$ {D^\alpha}\left(g \right)\left(0 \right) = \mathop {{\rm{lim}}}\limits_{t \to {0^ +}} {D^\alpha}\left(g \right)\left(t \right).

Let us consider two functions h, g that are α–differentiable, where α ∈ (0, 1]. Using the above definition, the following results can be stated [34]:

$D_{t}^{α} (ah + bg) (t) = a D_{t}^{α} (h (t)) + b D_{t}^{α} (g (t)), \forall a, b \in ℝ .$ D_t^\alpha \left({ah + bg} \right)\left(t \right) = aD_t^\alpha \left({h\left(t \right)} \right) + bD_t^\alpha \left({g\left(t \right)} \right),\forall a,b \in {\mathbb R}. .
$D_{t}^{α} (t^{n}) = n t^{n - α}, \forall n \in ℝ .$ D_t^\alpha \left({{t^n}} \right) = n{t^{n - \alpha }},\forall n \in {\mathbb R}.
$D_{t}^{α} (c) = 0, \forall h (t) = c .$ D_t^\alpha \left(c \right) = 0,\forall h\left(t \right) = c. .
$D_{t}^{α} (hg) (t) = g (t) D_{t}^{α} (h (t)) + h (t) D_{t}^{α} (g (t)) .$ D_t^\alpha \left({hg} \right)\left(t \right) = g\left(t \right)D_t^\alpha \left({h\left(t \right)} \right) + h\left(t \right)D_t^\alpha \left({g\left(t \right)} \right). .
$D_{t}^{α} (\frac{h}{g}) (t) = \frac{g (t) D_{t}^{α} (h (t)) - h (t) D_{t}^{α} (g (t))}{g {(t)}^{2}} .$ D_t^\alpha \left({{h \over g}} \right)\left(t \right) = {{g\left(t \right)D_t^\alpha \left({h\left(t \right)} \right) - h\left(t \right)D_t^\alpha \left({g\left(t \right)} \right)} \over {g{{\left(t \right)}^2}}}. .
$D_{t}^{α} (\frac{1}{α} t^{α}) = 1 .$ D_t^\alpha \left({{1 \over \alpha}{t^\alpha}} \right) = 1. .
$D_{t}^{α} (h \circ g) (t) = D_{g}^{α} h (g (t)) D_{t}^{α} g (t) .$ D_t^\alpha \left({h \circ g} \right)\left(t \right) = D_g^\alpha h\left({g\left(t \right)} \right)D_t^\alpha g\left(t \right). .

2.2

Describing of MCEM

Consider the NCPD equation in Eq.(5), where t and x are independent. (5) $F (u, u_{x}, u_{t}, D_{x}^{α} u, D_{t}^{α} u, \dots) = 0, α \in 0, 1],$ F\left({u,{u_x},{u_t},D_x^\alpha u,D_t^\alpha u, \cdots} \right) = 0,\;\;\;\;\alpha \in \left({0,1} \right], where u is an unknown function, $D_{x}^{α} u (resp . D_{t}^{α} u)$ D_x^\alpha u\,({\rm{resp}}.\,D_t^\alpha u) is the Katugampola's derivative of u of order α in respect to x (resp. t), and F is a polynomial in u.

The proposed modified conformal $(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})$ \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) expansion method to solve Eq.(5) is defined in five steps:

Step 1. Use the transformation given in Eq.(6) to transform the CPD equation into NCD equation. (6) $u (x, t) = u (ξ), ξ = \frac{k_{1}}{α} x^{α} + \frac{k_{2}}{α} t^{α},$ u\left({x,t} \right) = u\left(\xi \right),\;\;\;\;\;\;\;\;\xi = {{{k_1}} \over \alpha}{x^\alpha} + {{{k_2}} \over \alpha}{t^\alpha}, where k₁, k₂ ∈ ℝ will be determined later. Hence, Eq.(5) becomes a nonlinear ordinary differential equation with the new variable ξ. (7) $N (u, D_{ξ}^{α} u, D_{ξ}^{2 α} u, D_{ξ}^{3 α} u, \dots) = 0,$ N\left({u,D_\xi^\alpha u,D_\xi^{2\alpha}u,D_\xi^{3\alpha}u, \cdots} \right) = 0, where $D_{ξ}^{α}$ D_\xi^\alpha stands for differentiation with respect to ξ in the sense of Katugampola.
Step 2. Balance the term having the high-level derivative with nonlinear terms of equation (7) to find the value of m ∈ ℕ − {0} using the relation (8). (8) $Degree u^{p} {(\frac{d^{q} u}{d ξ^{q}})}^{s}] = mp + s (m + q) .$ Degree\left[ {{u^p}{{\left({{{{d^q}u} \over {d{\xi^q}}}} \right)}^s}} \right] = mp + s\left({m + q} \right).
Step 3. Suppose that solutions of equation (7) have the form given in Eq.(9): (9) $u (ξ) = \sum_{i = - m}^{m} l_{i} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{i},$ u(\xi) = \mathop \sum \limits_{i = - m}^m {l_i}{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^i}, where, l_i ∈ ℝ, ∀(i = −m, …, m) and $(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})$ \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) satisfies Eq.(10). (10) ${(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{'} = - {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{2} - λ (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) - μ .$ {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^{'}} = - {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^2} - \lambda \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) - \mu. Hence, Eq.(10) has several solutions defined as follows: (11) $\frac{D_{ξ}^{α} G (ξ)}{G (ξ)} = \begin{array}{l} - \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}); (λ^{2} - 4 μ) > 0, \\ - \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}); (λ^{2} - 4 μ) < 0, \\ - \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}}; (λ^{2} - 4 μ) = 0 . \end{array}$ {{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}} = \left\{{\matrix{{- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}{\rm{cosh}}\left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}{\rm{sinh}}\left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}{\rm{sinh}}\left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}{\rm{cosh}}\left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right);\;\;\;\;\left({{\lambda^2} - 4\mu} \right) > 0,} \hfill \cr {- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}{\rm{cos}}\left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}{\rm{sin}}\left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}{\rm{sin}}\left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}{\rm{cos}}\left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right);\;\;\;\;\left({{\lambda^2} - 4\mu} \right) < 0,} \hfill \cr {- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}};\;\;\;\;\left({{\lambda^2} - 4\mu} \right) = 0.\;\;\;\;} \hfill \cr}} \right.
Step 4. Find the system of equations for C₁,C₂, μ, λ, k₁, k₂, and l_i (i = −m, …, m) as follows: replace Eq.(9) in equation (7), use equation (10), and finally set coefficients of ${(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{i}$ {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^i} to zero.
Step 5. Resolve the obtained system, then replace the solutions of Eq.(10), all coefficients C₁,C₂, μ, λ, k₁, k₂, and l_i (i = −m, …, m) into Eq.(9), to get the solutions of the original Eq.(5).

3

Applications

The ability of the method to generate exact solutions is assessed using three different nonlinear fractional partial differential equations. Subsection 3.1 shows the performance of the method on the Whitham-Broer-Kaup equation. Subsection 3.2 presents the application of MCEM on the coupled Burgers equations. Finally, subsection 3.3 provides the efficiency of the method on the coupled modified Korteweg-de Vries (mKdv) equations.

3.1

The Whitham-Broer-Kaup equation

The governing equations in the conformal derivative are given by Eq.(12), in which β, γ are constants, u and v are functions of (x,t), and α ∈ (0, 1]. (12) $\begin{array}{l} D_{t}^{α} u + u D_{x}^{α} u + D_{x}^{α} v + β D_{x}^{2 α} u = 0, \\ D_{t}^{α} v + D_{x}^{α} (u v) - β D_{x}^{2 α} v + γ D_{x}^{3 α} u = 0 . \end{array}$ \left\{{\matrix{{D_t^\alpha u + uD_x^\alpha u + D_x^\alpha v + \beta D_x^{2\alpha}u = 0,} \hfill \cr {D_t^\alpha v + D_x^\alpha \left({uv} \right) - \beta D_x^{2\alpha}v + \gamma D_x^{3\alpha}u = 0.} \hfill \cr}} \right.

By employing the traveling wave transformation given in Eq.(13): (13) $v (x, t) = v (ξ), u (x, t) = u (ξ), ξ = \frac{k_{1} x^{α}}{α} + \frac{k_{2} t^{α}}{α},$ v\left({x,t} \right) = v\left(\xi \right),u\left({x,t} \right) = u\left(\xi \right),\xi = {{{k_1}{x^\alpha}} \over \alpha} + {{{k_2}{t^\alpha}} \over \alpha}, where k₁, k₂ ∈ ℝ, Eq.(12) becomes a nonlinear ordinary differential equation given in Eq.(14): (14) $\begin{array}{l} β k_{1}^{2} D_{ξ}^{2 α} u + k_{1} u D_{ξ}^{α} u + k_{1} D_{ξ}^{α} v + k_{2} D_{ξ}^{α} u = 0, \\ γ k_{1}^{3} D_{ξ}^{3 α} u - β k_{1}^{2} D_{ξ}^{2 α} v + k_{1} u D_{ξ}^{α} v + k_{1} v D_{ξ}^{α} u + k_{2} D_{ξ}^{α} v = 0 . \end{array}$ \left\{{\matrix{{\beta k_1^2D_\xi^{2\alpha}u + {k_1}uD_\xi^\alpha u + {k_1}D_\xi^\alpha v + {k_2}D_\xi^\alpha u = 0,} \hfill \cr {\gamma k_1^3D_\xi^{3\alpha}u - \beta k_1^2D_\xi^{2\alpha}v + {k_1}uD_\xi^\alpha v + {k_1}vD_\xi^\alpha u + {k_2}D_\xi^\alpha v = 0.} \hfill \cr}} \right.

By balancing the order of $D_{ξ}^{2 α} u$ D_\xi^{2\alpha}u and $u D_{ξ}^{α} u (resp . D_{ξ}^{3 α} u and u D_{ξ}^{α} v)$ uD_\xi^\alpha u\,({\rm{resp}}.\,D_\xi^{3\alpha}u\,{\rm{and}}\,uD_\xi^\alpha v) in the first (resp. second) equation of Eq.(14), we find n = 1 and m = 2. Hence, the solutions of Eq.(14) are given by Eq.(15), such that α_i ∈ ℝ (i = 0, 1, 2) and β _j ∈ ℝ (j = 0, 1, 2, 3, 4) will be obtained later. (15) $\begin{array}{l} u (ξ) = α_{0} + α_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + α_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ v (ξ) = β_{0} + β_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + β_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1} + β_{3} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{2} + β_{4} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 2} . \end{array}$ \left\{{\matrix{{u\left(\xi \right) = {\alpha_0} + {\alpha_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\alpha_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {v\left(\xi \right) = {\beta_0} + {\beta_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\beta_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}} + {\beta_3}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^2} + {\beta_4}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 2}}.} \hfill \cr}} \right.

The left side of Eq.(14) is transformed to polynomials in ${(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{j}$ {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^j} , (j = 0, ±1, ±2, …) by using Eq.(15). As a result, we get a system of equations for k₁, k₂, α_i(i = 0, 1, 2), and β_j(j = 0, 1, 2, 3, 4), by grouping the coefficients of these polynomials and setting them to zero. Maple software is used to generate three cases:

Case 1: $\begin{array}{l} α_{0} = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}, α_{1} = 0, α_{2} = 2 μ k_{1} \sqrt{β^{2} + γ}, k_{1} = k_{1}, k_{2} = k_{2}, \\ β_{0} = - 2 μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ)), β_{1} = 0, β_{2} = λ β_{0}, β_{3} = 0, β_{4} = μ β_{0} . \end{array}$ \left\{{\matrix{{{\alpha_0} = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}},{\alpha_1} = 0,{\alpha_2} = 2\mu {k_1}\sqrt {{\beta^2} + \gamma},{k_1} = {k_1},{k_2} = {k_2},} \hfill \cr {{\beta_0} = - 2\mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right),{\beta_1} = 0,{\beta_2} = \lambda {\beta_0},{\beta_3} = 0,{\beta_4} = \mu {\beta_0}.} \hfill \cr}} \right. hence (16) $\begin{array}{l} u_{1} (ξ) = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}} + 2 μ k_{1} \sqrt{β^{2} + γ} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ v_{1} (ξ) = - 2 μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ)) + λ β_{0} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1} + μ β_{0} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 2} . \end{array}$ \left\{{\matrix{{{u_1}\left(\xi \right) = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}} + 2\mu {k_1}\sqrt {{\beta^2} + \gamma} {{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {{v_1}\left(\xi \right) = - 2\mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right) + \lambda {\beta_0}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}} + \mu {\beta_0}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 2}}.} \hfill \cr}} \right.

Case 2: $\begin{array}{l} α_{0} = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}, α_{1} = 2 k_{1} \sqrt{β^{2} + γ}, α_{2} = 0, k_{1} = k_{1}, k_{2} = k_{2}, \\ β_{0} = μ β_{3}, β_{1} = λ β_{3}, β_{2} = 0, β_{3} = 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)), β_{4} = 0 . \end{array}$ \left\{{\matrix{{{\alpha_0} = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}},{\alpha_1} = 2{k_1}\sqrt {{\beta^2} + \gamma},{\alpha_2} = 0,{k_1} = {k_1},{k_2} = {k_2},} \hfill \cr {{\beta_0} = \mu {\beta_3},{\beta_1} = \lambda {\beta_3},{\beta_2} = 0,{\beta_3} = 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right),{\beta_4} = 0.} \hfill \cr}} \right. then (17) $\begin{array}{l} u_{2} (ξ) = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}} + 2 k_{1} \sqrt{β^{2} + γ} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}), \\ v_{2} (ξ) = μ β_{3} + λ β_{3} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)) {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{2} . \end{array}$ \left\{{\matrix{{{u_2}\left(\xi \right) = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}} + 2{k_1}\sqrt {{\beta^2} + \gamma} \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right),} \hfill \cr {{v_2}\left(\xi \right) = \mu {\beta_3} + \lambda {\beta_3}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right){{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^2}.} \hfill \cr}} \right.

Case 3: $\begin{array}{l} α_{0} = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}, α_{1} = 2 k_{1} \sqrt{β^{2} + γ}, α_{2} = μ α_{1}, k_{1} = k_{1}, k_{2} = k_{2}, β_{0} = 0, β_{1} = λ β_{3}, \\ β_{2} = - 2 λ μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ)), β_{3} = 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)), β_{4} = \frac{μ}{λ} β_{2} . \end{array}$ \left\{{\matrix{{{\alpha_0} = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}},{\alpha_1} = 2{k_1}\sqrt {{\beta^2} + \gamma},{\alpha_2} = \mu {\alpha_1},{k_1} = {k_1},{k_2} = {k_2},{\beta_0} = 0,{\beta_1} = \lambda {\beta_3},} \hfill \cr {{\beta_2} = - 2\lambda \mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right),{\beta_3} = 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right),{\beta_4} = {\mu \over \lambda}{\beta_2}.} \hfill \cr}} \right. therefore (18) $\begin{array}{l} u_{3} (ξ) = \frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}} + 2 k_{1} \sqrt{β^{2} + γ} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + μ α_{1} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ v_{3} (ξ) = λ β_{3} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) - 2 λ μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ)) {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1} + \\ 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)) {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{2} + \frac{μ}{λ} β_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 2} . \end{array}$ \left\{{\matrix{{{u_3}\left(\xi \right) = {{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}} + 2{k_1}\sqrt {{\beta^2} + \gamma} \left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + \mu {\alpha_1}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {{v_3}\left(\xi \right) = \lambda {\beta_3}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) - 2\lambda \mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right){{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}} +} \hfill \cr {2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right){{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^2} + {\mu \over \lambda}{\beta_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 2}}.} \hfill \cr}} \right.

Using the solutions given by Eq.(11), the Whitham-Broer-Kaup equation has the following new exact solutions:

Case (3.1). (λ² − 4μ) > 0, $\begin{matrix} u_{3.1} (ξ) = (\begin{array}{l} (\frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}) \\ + 2 k_{1} \sqrt{β^{2} + γ} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})})) \\ + μ α_{1} {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))}^{- 1} \end{array}), \\ v_{3.1} (ξ) = (\begin{array}{l} λ β_{3} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})})) \\ - \frac{2 λ μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ))}{(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))} \\ + 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)) \\ \times {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))}^{2} \\ + \frac{μ β_{2}}{λ {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))}^{2}} \end{array}), \\ ξ = (\frac{k_{1}}{α}) x^{α} + (\frac{k_{2}}{α}) t^{α} . \end{matrix}$ \matrix{{{u_{3.1}}\left(\xi \right) = \left({\matrix{{\left({{{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}}} \right)} \hfill \cr {+ 2{k_1}\sqrt {{\beta^2} + \gamma} \left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)} \hfill \cr {+ \mu {\alpha_1}{{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^{- 1}}} \hfill \cr}} \right),} \cr {{v_{3.1}}\left(\xi \right) = \left({\matrix{{\lambda {\beta_3}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)} \hfill \cr {- {{2\lambda \mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right)} \over {\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}}} \hfill \cr {+ 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right)} \hfill \cr {\times {{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^2}} \hfill \cr {+ {{\mu {\beta_2}} \over {\lambda {{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^2}}}} \hfill \cr}} \right),} \cr {\xi = \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{{k_2}} \over \alpha}} \right){t^\alpha}.} \cr}

Case (3.2). (λ² − 4μ) < 0, $\begin{matrix} u_{3.2} (ξ) = (\begin{array}{l} (\frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}) \\ + 2 k_{1} \sqrt{β^{2} + γ} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})})) \\ + μ α_{1} {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))}^{- 1} \end{array}), \\ v_{3.2} (ξ) = (\begin{array}{l} λ β_{3} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})})) \\ - \frac{2 λ μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ))}{(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))} \\ + 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)) \\ \times {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))}^{2} \\ + \frac{μ β_{2}}{λ {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))}^{2}} \end{array}), \\ ξ = (\frac{k_{1}}{α}) x^{α} + (\frac{k_{2}}{α}) t^{α} . \end{matrix}$ \matrix{{{u_{3.2}}\left(\xi \right) = \left({\matrix{{\left({{{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}}} \right)} \hfill \cr {+ 2{k_1}\sqrt {{\beta^2} + \gamma} \left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)} \hfill \cr {+ \mu {\alpha_1}{{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^{- 1}}} \hfill \cr}} \right),} \cr {{v_{3.2}}\left(\xi \right) = \left({\matrix{{\lambda {\beta_3}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)} \hfill \cr {- {{2\lambda \mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right)} \over {\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}}} \hfill \cr {+ 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right)} \hfill \cr {\times {{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^2}} \hfill \cr {+ {{\mu {\beta_2}} \over {\lambda {{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^2}}}} \hfill \cr}} \right),} \cr {\xi = \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{{k_2}} \over \alpha}} \right){t^\alpha}.} \cr}

Case (3.3). (λ² − 4μ) = 0, $\begin{array}{l} u_{3.3} (ξ) = (\frac{λ k_{1}^{2} (β^{2} + γ) - k_{2} \sqrt{β^{2} + γ}}{k_{1} \sqrt{β^{2} + γ}}) + 2 k_{1} \sqrt{β^{2} + γ} (- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}}) + μ α_{1} {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{- 1}, \\ v_{3.3} (ξ) = (\begin{array}{l} λ β_{3} (- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}}) - 2 λ μ k_{1}^{2} (β \sqrt{β^{2} + γ} + (β^{2} + γ)) {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{- 1} \\ + 2 k_{1}^{2} (β \sqrt{β^{2} + γ} - (β^{2} + γ)) {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{2} + \frac{1}{4} λ β_{2} {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{- 2} \end{array}), \\ ξ = (\frac{k_{1}}{α}) x^{α} + (\frac{k_{2}}{α}) t^{α} . \end{array}$ \left\{{\matrix{{{u_{3.3}}\left(\xi \right) = \left({{{\lambda k_1^2\left({{\beta^2} + \gamma} \right) - {k_2}\sqrt {{\beta^2} + \gamma}} \over {{k_1}\sqrt {{\beta^2} + \gamma}}}} \right) + 2{k_1}\sqrt {{\beta^2} + \gamma} \left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right) + \mu {\alpha_1}{{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^{- 1}},} \hfill \cr {{v_{3.3}}\left(\xi \right) = \left({\matrix{{\lambda {\beta_3}\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right) - 2\lambda \mu k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} + \left({{\beta^2} + \gamma} \right)} \right){{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^{- 1}}} \hfill \cr {+ 2k_1^2\left({\beta \sqrt {{\beta^2} + \gamma} - \left({{\beta^2} + \gamma} \right)} \right){{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^2} + {1 \over 4}\lambda {\beta_2}{{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^{- 2}}} \hfill \cr}} \right),} \hfill \cr {\xi = \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{{k_2}} \over \alpha}} \right){t^\alpha}.} \hfill \cr}} \right. The remaining solutions for Case 1 and Case 2 can be obtained analogously.

Figures (1,2,3) are the profiles of some solutions that are obtained as u_3.1 (ξ), v_3.1 (ξ), u_3.2 (ξ), v_3.2 (ξ), u_3.3 (ξ), and v_3.3 (ξ). These profiles were plotted using Maple 2022, providing a clearer understanding of the characteristics of each solution. The results show fascinating shapes of singular, periodic, and bright bell-shaped solitons for different parameters.

3.2

The coupled Burgers equations

The governing equations in the conformal derivative of the coupled Burgers equations are given by Eq.(19), where 0 < α ≤ 1, p and q are two constants, h = h(x,t) and g = g(x,t) are functions of x and t independent variables. (19) $\begin{array}{l} D_{t}^{α} h - D_{x}^{2 α} h + 2 h D_{x}^{α} h + p D_{x}^{α} (hg) = 0, \\ D_{t}^{α} g - D_{x}^{2 α} g + 2 g D_{x}^{α} g + q D_{x}^{α} (hg) = 0 . \end{array}$ \left\{{\matrix{{D_t^\alpha h - D_x^{2\alpha}h + 2hD_x^\alpha h + pD_x^\alpha \left({hg} \right) = 0,} \hfill \cr {D_t^\alpha g - D_x^{2\alpha}g + 2gD_x^\alpha g + qD_x^\alpha \left({hg} \right) = 0.} \hfill \cr}} \right. Using the following transformation: (20) $h (x, t) = h (ξ), g (x, t) = g (ξ), ξ = \frac{k_{1}}{α} x^{α} + \frac{k_{2}}{α} t^{α},$ h\left({x,t} \right) = h\left(\xi \right),\;\;\;\;g\left({x,t} \right) = g\left(\xi \right),\;\;\;\;\xi = {{{k_1}} \over \alpha}{x^\alpha} + {{{k_2}} \over \alpha}{t^\alpha}, where k₁ and k₂ are two constants, Eq.(19) becomes equivalent to Eq.(21). (21) $\begin{array}{l} k_{2} D_{ξ}^{α} h - k_{1}^{2} D_{ξ}^{2 α} h + 2 k_{1} h D_{ξ}^{α} h + p k_{1} g D_{ξ}^{α} h + p k_{1} h D_{ξ}^{α} g = 0, \\ k_{2} D_{ξ}^{α} g - k_{1}^{2} D_{ξ}^{2 α} g + 2 k_{1} g D_{ξ}^{α} g + q k_{1} g D_{ξ}^{α} h + q k_{1} h D_{ξ}^{α} g = 0 . \end{array}$ \left\{{\matrix{{{k_2}D_\xi^\alpha h - k_1^2D_\xi^{2\alpha}h + 2{k_1}hD_\xi^\alpha h + p{k_1}gD_\xi^\alpha h + p{k_1}hD_\xi^\alpha g = 0,} \hfill \cr {{k_2}D_\xi^\alpha g - k_1^2D_\xi^{2\alpha}g + 2{k_1}gD_\xi^\alpha g + q{k_1}gD_\xi^\alpha h + q{k_1}hD_\xi^\alpha g = 0.} \hfill \cr}} \right. If we balance both nonlinear terms and order of derivatives with high levels in both equations of Eq.(21), we find n = 1 and m = 1. Hence, the solutions of Eq.(21) are given by Eq.(22), where α_i ∈ ℝ (i = 0, 1, 2) and β_j ∈ ℝ (j = 0, 1, 2) will be obtained later. (22) $\begin{array}{l} h (ξ) = α_{0} + α_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + α_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ g (ξ) = β_{0} + β_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + β_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1} . \end{array}$ \left\{{\matrix{{h\left(\xi \right) = {\alpha_0} + {\alpha_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\alpha_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {g\left(\xi \right) = {\beta_0} + {\beta_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\beta_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}}.} \hfill \cr}} \right.

Substituting Eq.(22) into Eq.(21), the left side is transformed into polynomials in ${(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{j}$ {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^j} where j = 0, ±1, ±2, …. As a result, we get a system of equations for k₁, k₂, α_i(i = 0, 1, 2) and β_j(j = 0, 1, 2) by setting the coefficients of these polynomials to zero. The Maple software is used to generate two cases:

Case 1: $\begin{array}{l} α_{0} = α_{0}, α_{1} = α_{1}, α_{2} = 0, β_{0} = α_{0} (\frac{q - 1}{p - 1}), β_{1} = α_{1} (\frac{q - 1}{p - 1}), β_{2} = 0, \\ k_{1} = - α_{1} (\frac{p q - 1}{p - 1}), k_{2} = \frac{α_{1} {(p q - 1)}^{2} (2 α_{0} - λ α_{1})}{{(p - 1)}^{2}} . \end{array}$ \left\{{\matrix{{{\alpha_0} = {\alpha_0},{\alpha_1} = {\alpha_1},{\alpha_2} = 0,{\beta_0} = {\alpha_0}\left({{{q - 1} \over {p - 1}}} \right),{\beta_1} = {\alpha_1}\left({{{q - 1} \over {p - 1}}} \right),{\beta_2} = 0,} \hfill \cr {{k_1} = - {\alpha_1}\left({{{pq - 1} \over {p - 1}}} \right),{k_2} = {{{\alpha_1}{{\left({pq - 1} \right)}^2}\left({2{\alpha_0} - \lambda {\alpha_1}} \right)} \over {{{\left({p - 1} \right)}^2}}}.} \hfill \cr}} \right. hence (23) $\begin{array}{l} h_{1} (ξ) = α_{0} + α_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}), \\ g_{1} (ξ) = α_{0} (\frac{q - 1}{p - 1}) + α_{1} (\frac{q - 1}{p - 1}) (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}), \\ ξ = \frac{- α_{1} (\frac{p q - 1}{p - 1})}{α} x^{α} + \frac{\frac{α_{1} {(p q - 1)}^{2} (2 α_{0} - λ α_{1})}{{(p - 1)}^{2}}}{α} t^{α} . \end{array}$ \left\{{\matrix{{{h_1}\left(\xi \right) = {\alpha_0} + {\alpha_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right),} \hfill \cr {{g_1}\left(\xi \right) = {\alpha_0}\left({{{q - 1} \over {p - 1}}} \right) + {\alpha_1}\left({{{q - 1} \over {p - 1}}} \right)\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right),} \hfill \cr {\xi = {{- {\alpha_1}\left({{{pq - 1} \over {p - 1}}} \right)} \over \alpha}{x^\alpha} + {{{{{\alpha_1}{{\left({pq - 1} \right)}^2}\left({2{\alpha_0} - \lambda {\alpha_1}} \right)} \over {{{\left({p - 1} \right)}^2}}}} \over \alpha}{t^\alpha}.} \hfill \cr}} \right.

Case 2: $\begin{array}{l} α_{0} = β_{0} (\frac{p - 1}{q - 1}), α_{1} = 0, α_{2} = α_{2}, β_{0} = β_{0}, β_{1} = 0, β_{2} = α_{2} (\frac{q - 1}{p - 1}), \\ k_{1} = \frac{α_{2}}{μ} (\frac{p q - 1}{p - 1}), k_{2} = \frac{α_{2} {(p q - 1)}^{2} (λ α_{2} (q - 1) + 2 μ β_{0} (1 - p))}{μ^{2} {(p - 1)}^{2} (q - 1)} . \end{array}$ \left\{{\matrix{{{\alpha_0} = {\beta_0}\left({{{p - 1} \over {q - 1}}} \right),{\alpha_1} = 0,{\alpha_2} = {\alpha_2},{\beta_0} = {\beta_0},{\beta_1} = 0,{\beta_2} = {\alpha_2}\left({{{q - 1} \over {p - 1}}} \right),} \hfill \cr {{k_1} = {{{\alpha_2}} \over \mu}\left({{{pq - 1} \over {p - 1}}} \right),{k_2} = {{{\alpha_2}{{\left({pq - 1} \right)}^2}\left({\lambda {\alpha_2}\left({q - 1} \right) + 2\mu {\beta_0}\left({1 - p} \right)} \right)} \over {{\mu^2}{{\left({p - 1} \right)}^2}\left({q - 1} \right)}}.} \hfill \cr}} \right. then (24) $\begin{array}{l} h_{2} (ξ) = β_{0} (\frac{p - 1}{q - 1}) + α_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ g_{2} (ξ) = β_{0} + α_{2} (\frac{q - 1}{p - 1}) {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ ξ = \frac{\frac{α_{2}}{μ} (\frac{p q - 1}{p - 1})}{α} x^{α} + \frac{\frac{α_{2} {(p q - 1)}^{2} (λ α_{2} (q - 1) + 2 μ β_{0} (1 - p))}{μ^{2} {(p - 1)}^{2} (q - 1)}}{α} t^{α}, \end{array}$ \left\{{\matrix{{{h_2}\left(\xi \right) = {\beta_0}\left({{{p - 1} \over {q - 1}}} \right) + {\alpha_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {{g_2}\left(\xi \right) = {\beta_0} + {\alpha_2}\left({{{q - 1} \over {p - 1}}} \right){{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {\xi = {{{{{\alpha_2}} \over \mu}\left({{{pq - 1} \over {p - 1}}} \right)} \over \alpha}{x^\alpha} + {{{{{\alpha_2}{{\left({pq - 1} \right)}^2}\left({\lambda {\alpha_2}\left({q - 1} \right) + 2\mu {\beta_0}\left({1 - p} \right)} \right)} \over {{\mu^2}{{\left({p - 1} \right)}^2}\left({q - 1} \right)}}} \over \alpha}{t^\alpha},} \hfill \cr}} \right.

Using Case 2 and Eq.(17), the coupled Burgers equations have the following new exact solutions:

Case (2.1). (λ² − 4μ) < 0, $\begin{array}{l} h_{2.1} (ξ) = β_{0} (\frac{p - 1}{q - 1}) + α_{2} {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))}^{- 1}, \\ g_{2.1} (ξ) = β_{0} + α_{2} (\frac{q - 1}{p - 1}) {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}))}^{- 1}, \\ ξ = (\frac{\frac{α_{2}}{μ} (\frac{p q - 1}{p - 1})}{α}) x^{α} + (\frac{\frac{α_{2} {(p q - 1)}^{2} (λ α_{2} (q - 1) + 2 μ β_{0} (1 - p))}{μ^{2} {(p - 1)}^{2} (q - 1)}}{α}) t^{α} . \end{array}$ \left\{{\matrix{{{h_{2.1}}\left(\xi \right) = {\beta_0}\left({{{p - 1} \over {q - 1}}} \right) + {\alpha_2}{{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^{- 1}},} \hfill \cr {{g_{2.1}}\left(\xi \right) = {\beta_0} + {\alpha_2}\left({{{q - 1} \over {p - 1}}} \right){{\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right)}^{- 1}},} \hfill \cr {\xi = \left({{{{{{\alpha_2}} \over \mu}\left({{{pq - 1} \over {p - 1}}} \right)} \over \alpha}} \right){x^\alpha} + \left({{{{{{\alpha_2}{{\left({pq - 1} \right)}^2}\left({\lambda {\alpha_2}\left({q - 1} \right) + 2\mu {\beta_0}\left({1 - p} \right)} \right)} \over {{\mu^2}{{\left({p - 1} \right)}^2}\left({q - 1} \right)}}} \over \alpha}} \right){t^\alpha}.} \hfill \cr}} \right.

Case (2.2). (λ² − 4μ) = 0, $\begin{array}{l} h_{2.2} (ξ) = β_{0} (\frac{p - 1}{q - 1}) + α_{2} {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))}^{- 1}, \\ g_{2.2} (ξ) = β_{0} + α_{2} (\frac{q - 1}{p - 1}) {(- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}))}^{- 1}, \\ ξ = (\frac{\frac{α_{2}}{μ} (\frac{p q - 1}{p - 1})}{α}) x^{α} + (\frac{\frac{α_{2} {(p q - 1)}^{2} (λ α_{2} (q - 1) + 2 μ β_{0} (1 - p))}{μ^{2} {(p - 1)}^{2} (q - 1)}}{α}) t^{α} . \end{array}$ \left\{ {\matrix{{{h_{2.2}}\left(\xi \right) = {\beta_0}\left({{{p - 1} \over {q - 1}}} \right) + {\alpha_2}{{\left({ - {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right)}}} \right)} \right)}^{ - 1}},} \hfill \cr {{g_{2.2}}\left(\xi \right) = {\beta_0} + {\alpha_2}\left({{{q - 1} \over {p - 1}}} \right){{\left({ - {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha }} \right)}}} \right)} \right)}^{ - 1}},} \hfill \cr {\xi = \left({{{{{{\alpha_2}} \over \mu }\left({{{pq - 1} \over {p - 1}}} \right)} \over \alpha }} \right){x^\alpha } + \left({{{{{{\alpha_2}{{\left({pq - 1} \right)}^2}\left({\lambda {\alpha_2}\left({q - 1} \right) + 2\mu {\beta_0}\left({1 - p} \right)} \right)} \over {{\mu^2}{{\left({p - 1} \right)}^2}\left({q - 1} \right)}}} \over \alpha }} \right){t^\alpha }.} \hfill \cr } } \right.

Case (2.3). (λ² − 4μ) = 0, $\begin{array}{l} h_{2.3} (ξ) = β_{0} (\frac{p - 1}{q - 1}) + α_{2} {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{- 1}, \\ g_{2.3} (ξ) = β_{0} + α_{2} (\frac{q - 1}{p - 1}) {(- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}})}^{- 1}, \\ ξ = (\frac{\frac{α_{2}}{μ} (\frac{p q - 1}{p - 1})}{α}) x^{α} + (\frac{\frac{α_{2} {(p q - 1)}^{2} (λ α_{2} (q - 1) + 2 μ β_{0} (1 - p))}{μ^{2} {(p - 1)}^{2} (q - 1)}}{α}) t^{α} . \end{array}$ \left\{{\matrix{{{h_{2.3}}\left(\xi \right) = {\beta_0}\left({{{p - 1} \over {q - 1}}} \right) + {\alpha_2}{{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^{- 1}},} \hfill \cr {{g_{2.3}}\left(\xi \right) = {\beta_0} + {\alpha_2}\left({{{q - 1} \over {p - 1}}} \right){{\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right)}^{- 1}},} \hfill \cr {\xi = \left({{{{{{\alpha_2}} \over \mu}\left({{{pq - 1} \over {p - 1}}} \right)} \over \alpha}} \right){x^\alpha} + \left({{{{{{\alpha_2}{{\left({pq - 1} \right)}^2}\left({\lambda {\alpha_2}\left({q - 1} \right) + 2\mu {\beta_0}\left({1 - p} \right)} \right)} \over {{\mu^2}{{\left({p - 1} \right)}^2}\left({q - 1} \right)}}} \over \alpha}} \right){t^\alpha}.} \hfill \cr}} \right. Similarly, the remaining solutions for Case 1 can be obtained analogously.

Figures (4,5,6) are shown below: h_2.1 (ξ), h_2.2 (ξ), h_2.3 (ξ), g_2.1 (ξ), g_2.2 (ξ) and g_2.3 (ξ). This information was plotted using Maple 2022, which helped us see each solution's characteristics more clearly.

3.3

The coupled mKdV equation

The governing equations in the conformal derivative of the coupled mKdV equation are given by Eq.(25), (25) $\begin{array}{l} D_{t}^{α} ψ = \frac{1}{2} D_{x}^{3 α} ψ - 3 ψ^{2} D_{x}^{α} ψ + \frac{3}{2} D_{x}^{2 α} ψ + 3 D_{x}^{α} (ψ ϕ) - 3 δ D_{x}^{α} ψ, \\ D_{t}^{α} ϕ = - D_{x}^{3 α} ϕ - 3 ϕ D_{x}^{α} ϕ - 3 D_{x}^{α} ψ D_{x}^{α} ϕ + 3 ψ^{2} D_{x}^{α} ϕ + 3 δ D_{x}^{α} ϕ, \end{array}$ \left\{{\matrix{{D_t^\alpha \psi = {1 \over 2}D_x^{3\alpha}\psi - 3{\psi^2}D_x^\alpha \psi + {3 \over 2}D_x^{2\alpha}\psi + 3D_x^\alpha \left({\psi \phi} \right) - 3\delta D_x^\alpha \psi,} \hfill \cr {D_t^\alpha \phi = - D_x^{3\alpha}\phi - 3\phi D_x^\alpha \phi - 3D_x^\alpha \psi D_x^\alpha \phi + 3{\psi^2}D_x^\alpha \phi + 3\delta D_x^\alpha \phi,} \hfill \cr}} \right. where 0 < α ≤ 1, δ is a constant, ψ = ψ(x,t) and ϕ = ϕ(x,t) are functions of x and t independent variables. We consider the following traveling wave transformation where k₁, k₂ are two constants: (26) $ψ (x, t) = ψ (ξ), ϕ (x, t) = ϕ (ξ), ξ = \frac{k_{1}}{α} x^{α} + \frac{k_{2}}{α} t^{α},$ \psi \left({x,t} \right) = \psi \left(\xi \right),\;\;\phi \left({x,t} \right) = \phi \left(\xi \right),\;\;\xi = {{{k_1}} \over \alpha}{x^\alpha} + {{{k_2}} \over \alpha}{t^\alpha},

We find the following nonlinear conformal derivative ordinary differential equations by replacing Eq.(26) in Eq.(25): (27) $\begin{array}{l} k_{2} D_{ξ}^{α} ψ = \frac{1}{2} k_{1}^{3} D_{ξ}^{3 α} ψ - 3 k_{1} ψ^{2} D_{ξ}^{α} ψ + \frac{3}{2} k_{1}^{2} D_{ξ}^{2 α} ψ + 3 k_{1} ψ D_{ξ}^{α} ϕ + 3 k_{1} ϕ D_{ξ}^{α} ψ - 3 δ k_{1} D_{ξ}^{α} ψ, \\ k_{2} D_{ξ}^{α} ϕ = - k_{1}^{3} D_{ξ}^{3 α} ϕ - 3 k_{1} ϕ D_{ξ}^{α} ϕ - 3 k_{1}^{2} D_{ξ}^{α} ψ D_{ξ}^{α} ϕ + 3 k_{1} ψ^{2} D_{ξ}^{α} ϕ + 3 δ k_{1} D_{ξ}^{α} ϕ . \end{array}$ \left\{{\matrix{{{k_2}D_\xi^\alpha \psi = {1 \over 2}k_1^3D_\xi^{3\alpha}\psi - 3{k_1}{\psi^2}D_\xi^\alpha \psi + {3 \over 2}k_1^2D_\xi^{2\alpha}\psi + 3{k_1}\psi D_\xi^\alpha \phi + 3{k_1}\phi D_\xi^\alpha \psi - 3\delta {k_1}D_\xi^\alpha \psi,} \hfill \cr {{k_2}D_\xi^\alpha \phi = - k_1^3D_\xi^{3\alpha}\phi - 3{k_1}\phi D_\xi^\alpha \phi - 3k_1^2D_\xi^\alpha \psi D_\xi^\alpha \phi + 3{k_1}{{\psi}^2}D_\xi^\alpha \phi + 3\delta {k_1}D_\xi^\alpha \phi.} \hfill \cr}} \right.

We balance both nonlinear terms and order of derivatives with high levels in both equations of Eq.(27) to get n = 1 and m = 1. Thus, the solutions of Eq.(27) are given by Eq.(28), where α_i ∈ ℝ (i = 0, 1, 2) and β_j ∈ ℝ (j = 0, 1, 2) will be obtained later. (28) $\begin{array}{l} ψ (ξ) = α_{0} + α_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + α_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ ϕ (ξ) = β_{0} + β_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}) + β_{2} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1} . \end{array}$ \left\{{\matrix{{\psi \left(\xi \right) = {\alpha_0} + {\alpha_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\alpha_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {\phi \left(\xi \right) = {\beta_0} + {\beta_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right) + {\beta_2}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}}.} \hfill \cr}} \right.

By replacing Eq.(28) in Eq.(27), the left side of Eq.(27) is transformed into polynomials in ${(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{j}$ {\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)^j} where j = 0, ±1, ±2, …. As a result, we get a system of equations for k₁, k₂, α_i, i ∈ {0, 1, 2} and β_j, (j ∈ {0, 1, 2, 3, 4}) by setting the coefficients of these polynomials to zero. The Maple software is used to generate two cases:

Case 1: $\begin{array}{l} α_{0} = \frac{1}{2} (λ k_{1} + 1), α_{1} = k_{1}, α_{2} = 0, β_{0} = \frac{1}{2} (2 δ + λ k_{1}), \\ β_{1} = k_{1}, β_{2} = 0, k_{1} = k_{1}, k_{2} = \frac{1}{4} k_{1} (3 - k_{1}^{2} (λ^{2} - 4 μ)) . \end{array}$ \left\{{\matrix{{{\alpha_0} = {1 \over 2}\left({\lambda {k_1} + 1} \right),{\alpha_1} = {k_1},{\alpha_2} = 0,{\beta_0} = {1 \over 2}\left({2\delta + \lambda {k_1}} \right),} \hfill \cr {{\beta_1} = {k_1},{\beta_2} = 0,{k_1} = {k_1},{k_2} = {1 \over 4}{k_1}\left({3 - k_1^2\left({{\lambda^2} - 4\mu} \right)} \right).} \hfill \cr}} \right. then $\begin{array}{l} ψ_{1.1} (ξ) = \frac{1}{2} (λ k_{1} + 1) + k_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}), \\ ϕ_{1.1} (ξ) = \frac{1}{2} (2 δ + λ k_{1}) + k_{1} (\frac{D_{ξ}^{α} G (ξ)}{G (ξ)}), \\ ξ = k_{1} \frac{{(x)}^{α}}{α} + \frac{1}{4} k_{1} (3 - k_{1}^{2} (λ^{2} - 4 μ)) \frac{{(t)}^{α}}{α} . \end{array}$ \left\{{\matrix{{{\psi_{1.1}}\left(\xi \right) = {1 \over 2}\left({\lambda {k_1} + 1} \right) + {k_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right),} \hfill \cr {{\phi_{1.1}}\left(\xi \right) = {1 \over 2}\left({2\delta + \lambda {k_1}} \right) + {k_1}\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right),} \hfill \cr {\xi = {k_1}{{{{\left(x \right)}^\alpha}} \over \alpha} + {1 \over 4}{k_1}\left({3 - k_1^2\left({{\lambda^2} - 4\mu} \right)} \right){{{{\left(t \right)}^\alpha}} \over \alpha}.} \hfill \cr}} \right.

Case 2: $\begin{array}{l} α_{0} = \frac{1}{2} (1 - λ k_{1}), α_{1} = 0, α_{2} = - μ k_{1}, β_{0} = \frac{1}{2} (2 δ - λ k_{1}), \\ β_{1} = 0, β_{2} = - μ k_{1}, k_{1} = k_{1}, k_{2} = \frac{1}{4} k_{1} (3 - k_{1}^{2} (λ^{2} - 4 μ)) . \end{array}$ \left\{{\matrix{{{\alpha_0} = {1 \over 2}\left({1 - \lambda {k_1}} \right),{\alpha_1} = 0,{\alpha_2} = - \mu {k_1},{\beta_0} = {1 \over 2}\left({2\delta - \lambda {k_1}} \right),} \hfill \cr {{\beta_1} = 0,{\beta_2} = - \mu {k_1},{k_1} = {k_1},{k_2} = {1 \over 4}{k_1}\left({3 - k_1^2\left({{\lambda^2} - 4\mu} \right)} \right).} \hfill \cr}} \right. hence $\begin{array}{l} ψ_{2.1} (ξ) = \frac{1}{2} (1 - λ k_{1}) - μ k_{1} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ ϕ_{2.1} (ξ) = \frac{1}{2} (2 δ - λ k_{1}) - μ k_{1} {(\frac{D_{ξ}^{α} G (ξ)}{G (ξ)})}^{- 1}, \\ ξ = k_{1} \frac{{(x)}^{α}}{α} + \frac{1}{4} k_{1} (3 - k_{1}^{2} (λ^{2} - 4 μ)) \frac{{(t)}^{α}}{α} . \end{array}$ \left\{{\matrix{{{\psi_{2.1}}\left(\xi \right) = {1 \over 2}\left({1 - \lambda {k_1}} \right) - \mu {k_1}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {{\phi_{2.1}}\left(\xi \right) = {1 \over 2}\left({2\delta - \lambda {k_1}} \right) - \mu {k_1}{{\left({{{D_\xi^\alpha G\left(\xi \right)} \over {G\left(\xi \right)}}} \right)}^{- 1}},} \hfill \cr {\xi = {k_1}{{{{\left(x \right)}^\alpha}} \over \alpha} + {1 \over 4}{k_1}\left({3 - k_1^2\left({{\lambda^2} - 4\mu} \right)} \right){{{{\left(t \right)}^\alpha}} \over \alpha}.} \hfill \cr}} \right.

Using Case 1 and Eq.(17), the coupled mKdV equation has the following new exact solutions:

Case (1.1). (λ² − 4μ) > 0, $\begin{array}{l} ψ_{1.1} (ξ) = \frac{1}{2} (1 + λ k_{1}) + k_{1} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})})), \\ ϕ_{1.1} (ξ) = \frac{1}{2} (λ k_{1} + 2 δ) + k_{1} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} (\frac{C_{1} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})}{C_{1} sinh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α}) + C_{2} \cosh (\frac{1}{2} \sqrt{\frac{(λ^{2} - 4 μ)}{α^{2}}} ξ^{α})})), \\ = (\frac{k_{1}}{α}) x^{α} + (\frac{\frac{1}{4} k_{1} (k_{1}^{2} (4 μ - λ^{2}) + 3)}{α}) t^{α} . \end{array}$ \left\{{\matrix{{{\psi_{1.1}}\left(\xi \right) = {1 \over 2}\left({1 + \lambda {k_1}} \right) + {k_1}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right),} \hfill \cr {{\phi_{1.1}}\left(\xi \right) = {1 \over 2}\left({\lambda {k_1} + 2\delta} \right) + {k_1}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sinh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cosh \left({{1 \over 2}\sqrt {{{\left({{\lambda^2} - 4\mu} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right),} \hfill \cr {= \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{{1 \over 4}{k_1}\left({k_1^2\left({4\mu - {\lambda^2}} \right) + 3} \right)} \over \alpha}} \right){t^\alpha}.} \hfill \cr}} \right.

Case (1.2). (λ² − 4μ) < 0, $\begin{array}{l} ψ_{1.2} (ξ) = \frac{1}{2} (1 + λ k_{1}) + k_{1} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})})), \\ ϕ_{1.2} (ξ) = \frac{1}{2} (λ k_{1} + 2 δ) + k_{1} (- \frac{λ}{2} + \frac{α}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} (\frac{C_{1} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) - C_{2} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})}{C_{1} sin (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α}) + C_{2} \cos (\frac{1}{2} \sqrt{\frac{(4 μ - λ^{2})}{α^{2}}} ξ^{α})})), \\ ξ = (\frac{k_{1}}{α}) x^{α} + (\frac{\frac{1}{4} k_{1} (k_{1}^{2} (4 μ - λ^{2}) + 3)}{α}) t^{α} . \end{array}$ \left\{{\matrix{{{\psi_{1.2}}\left(\xi \right) = {1 \over 2}\left({1 + \lambda {k_1}} \right) + {k_1}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right),} \hfill \cr {{\phi_{1.2}}\left(\xi \right) = {1 \over 2}\left({\lambda {k_1} + 2\delta} \right) + {k_1}\left({- {\lambda \over 2} + {\alpha \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} \left({{{{C_1}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) - {C_2}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)} \over {{C_1}\sin \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right) + {C_2}\cos \left({{1 \over 2}\sqrt {{{\left({4\mu - {\lambda^2}} \right)} \over {{\alpha^2}}}} {\xi^\alpha}} \right)}}} \right)} \right),} \hfill \cr {\xi = \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{{1 \over 4}{k_1}\left({k_1^2\left({4\mu - {\lambda^2}} \right) + 3} \right)} \over \alpha}} \right){t^\alpha}.} \hfill \cr}} \right.

Case (1.3). (λ² − 4μ) = 0, $\begin{array}{l} ψ_{1.3} (ξ) = \frac{1}{2} (1 + λ k_{1}) + k_{1} (- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}}), \\ ϕ_{1.3} (ξ) = \frac{1}{2} (λ k_{1} + 2 δ) + k_{1} (- \frac{λ}{2} + \frac{C_{2} α}{C_{1} + C_{2} ξ^{α}}), \\ ξ = (\frac{k_{1}}{α}) x^{α} + (\frac{3 k_{1}}{4 α}) t^{α} . \end{array}$ \left\{{\matrix{{{\psi_{1.3}}\left(\xi \right) = {1 \over 2}\left({1 + \lambda {k_1}} \right) + {k_1}\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right),} \hfill \cr {{\phi_{1.3}}\left(\xi \right) = {1 \over 2}\left({\lambda {k_1} + 2\delta} \right) + {k_1}\left({- {\lambda \over 2} + {{{C_2}\alpha} \over {{C_1} + {C_2}{\xi^\alpha}}}} \right),} \hfill \cr {\xi = \left({{{{k_1}} \over \alpha}} \right){x^\alpha} + \left({{{3{k_1}} \over {4\alpha}}} \right){t^\alpha}.} \hfill \cr}} \right.

The remaining solutions for Case 2 can be obtained analogously.

We present the Figures (7,8,9) of some solutions obtained as ψ_3.1 (ξ), ϕ_3.1 (ξ), ψ_3.2 (ξ), ϕ_3.2 (ξ), ψ_3.3 (ξ), and ϕ_3.3 (ξ). These profiles were graphed using Maple 2022, aiding in a more transparent comprehension of their characteristics.

4

Conclusion

In this research paper, we presented the general properties of MCEM for solving various NPD equations. This method is based on the transforming the nonlinear equations into ordinary differential equations by utilizing Katugampola's derivative and the complex transform. The proposed method is effectively used to find exact solutions for three distinct space-time conformal derivative nonlinear equations: the Whitham-Broer-Kaup equations, the Burgers equations, and the modified Korteweg-de Vries equations. The exact solutions included singular kink wave solution, singular anti kink wave solution, multiple-soliton solution, and singular soliton were extracted. Under the suitable values of parameters in solutions, various wave simulations of results of this paper were simulated.

A new modified conformal expansion method to some nonlinear conformal partial differential equations

Full Article

Paradigm

My account