New Solitary Wave Solutions for Variants of (3+1)-Dimensional Wazwaz-Benjamin-Bona-Mahony Equations

We solve distinct forms of (3+1)-Dimensional Wazwaz-Benjamin-Bona-Mahony [(3+1)-Dimensional WBBM] equations by employing the method of Sardar-subequation. When parameters involving this approach are taken to be special values, we can obtain the solitary wave solutions (sws) which is concluded from other approaches such as the functional variable method, the trail equation method, the first integral method and so on. We obtain new and general solitary wave solutions in terms of generalized hyperbolic and trigonometric functions. The results demonstrate the power of the proposed method for the determination of sws of non-linear evolution equations (NLEs).

In this paper, we utilize of an effective and efficient technique for manufacturing a range of solitary wave solutions for the variants of the (3+1)-dimensional WBBM equations.
The Benjamin-Bona-Mahony (BBM) equation proposed in [40] as a model to study the approximately the unidirectional propagation of smallamplitude long waves on certain non-linear dispersive systems as a good alternative to the KdV [41]. It is used in modeling surface waves of long wavelength in liquids; it covers hydro-magnetic waves in cold plasma, acoustic waves in anharmonic crystals, and acoustic gravity waves in compressible fluids. Many efforts have been offered to modified forms of this equation such as known as the modified Benjamin-Bona-Mahoney (mBBM) equation [42]. Furthermore, as higher dimensional samples convert to be more realistic; various modifications to Equation (2) have been proposed in the literature among which the (3+1)-dimensional mBBM equation given in Equations (3)(4)(5) by Wazwaz [43]. Using the tanh/sech method, Wazwaz obtained soliton, kink, and periodic solutions for the following three different types of Equation (2) and Based on these ideas, this paper is organized as follows. Section The Sardar-Subequation Method introduces a brief description of the Sardar-sub-equation method. Section The (3+1)-Dimensional WBBM Equation discusses the application of the Sardar-subequation method to variants of (3+1)dimensional WBBM equations represent by (3)-(5). The graphical presentation for the acquired solution is given in section Exact Solutions of the (3+1)-Dimensional WBBM Equation.
We complete the paper with conclusions part.

THE SARDAR-SUBEQUATION METHOD
We consider the following NLEs where u = u(x, t) is an unknown function and G is a polynomial in u and its partial derivatives.
To solve (Equation 6), we take the traveling wave transformation (twt) where c = 0 is constants to be determined later. By using (7), Equation (6) is turned into following ODE w.r.t. η in which U = U(η), U ′ = dU dη , U ′′ = d 2 U dη 2 , . . . . Suppose that the Equation (8) has a solution of the form where ̟ i , (i = 0, 1, . . . , s) are coefficients to be determined with (̟ s = 0) and ϕ(η) satisfies the ODE in the form where a and ρ are real constants. The solutions of ODE (10) are Case I: If a > 0 and ρ = 0 then Case II: If a < 0 and ρ = 0 then Case III: If a < 0 and ρ = a 2 4b then Case IV: If a > 0 and ρ = a 2 4 then The procedure starts by determining s by the assistance of the classical balance rule. When s is determined the predicted solution (8) is substituted into (Equation 7). Since we seek a non-zero solution that is ϕ(η) = 0, all the coefficient of power of ϕ(η) are equated to zero. Thus, the resultant algebraic system is solved for ̟ i s and c. One should note that ̟ s = 0 and c = 0 both have to be satisfied in the solution sets. Whenever ̟ i s and c are determined the solutions are constructed by using these parameters explicitly. The procedure ends by substitution of η = x − ct into the solutions satisfying (Equation 7). Remark 1. The Sardar-subequation method is considered among those general ones from which, under certain cases, various methods can be deduced such as the functional variable method [44], the first integral method [45], and so on.
The Third (3+1)-Dimensional WBBM Equation We next study the third (3+1)-dimensional WBBM equation (5). By making the traveling wave transformation Integrating (Equation 20) twice and setting the integration constants to zero yield Then, Equation (21) can be written as.

EXACT SOLUTIONS OF THE (3+1)-DIMENSIONAL WBBM EQUATION
Balancing U with U ′′ in Equations (13,17,22), we get s = 1. Therefore, the solutions form of Equations (13,17,22) has the following expression Substituting Equation (24) into (13) [or (17) or (22)] with along (Equation 10) and equating all the coefficients of ϕ(η) to zero, we obtain a highly complicated system of algebraic equations. A set of algebraic equations is obtained in ̟ 0 , ̟ 1 , and a as follows: Solving these algebraic equations with Maple, we acquire: The First (3+1)-Dimensional WBBM Equation

GRAPHICAL PRESENTATION
Graph is a strong tool for relationship describing clarity the solutions of the challenges. Therefore, some graphs of the

CONCLUSIONS
We have proposed a method, namely Sardar-subequation method to solve NLEs with the help of Maple Software. Distinct forms of (3+1)-dimensional WBBM equations are handled to display the effectiveness of the suggested method. From our results, some sws are obtained including the generalized hyperbolic and trigonometric function solutions. As far as we know, for the first time, we describe and introduce Sardarsubequation method which is a new method for solving NLEs.
Thus, all the solutions distinct forms of (3+1)-dimensional WBBM equations are new, which cannot be found in literature to our best knowledge. We can also see that the approach used in this letter is very effective, powerful and convenient and can be steadily applied to NLEs. We will extend the proposed method for some fractional models [46][47][48] in a future work.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.