Analysis of parametric effects in the wave profile of the variant Boussinesq equation through two analytical approaches

Shao-Wen Yao; Md. Ekramul Islam; Md. Ali Akbar; Mustafa Inc; Mohamed Adel; Mohamed S. Osman

doi:10.1515/phys-2022-0071

Article Open Access

Analysis of parametric effects in the wave profile of the variant Boussinesq equation through two analytical approaches

Shao-Wen Yao , Md. Ekramul Islam , Md. Ali Akbar , Mustafa Inc , Mohamed Adel and Mohamed S. Osman

Published/Copyright: August 14, 2022

Published by

Become an author with De Gruyter Brill

Submit Manuscript Author Information Explore this Subject

From the journal Open Physics Volume 20 Issue 1

Abstract

The variant Boussinesq equation has significant application in propagating long waves on the surface of the liquid layer under gravity action. In this article, the improved Bernoulli subequation function (IBSEF) method and the new auxiliary equation (NAE) technique are introduced to establish general solutions, some fundamental soliton solutions accessible in the literature, and some archetypal solitary wave solutions that are extracted from the broad-ranging solution to the variant Boussinesq wave equation. The established soliton solutions are knowledgeable and obtained as a combination of hyperbolic, exponential, rational, and trigonometric functions, and the physical significance of the attained solutions is speculated for the definite values of the included parameters by depicting the 3D profiles and interpreting the physical incidents. The wave profile represents different types of waves associated with the free parameters that are related to the wave number and velocity of the solutions. The obtained solutions and graphical representations visualize the dynamics of the phenomena and build up the mathematical foundation of the wave process in dissipative and dispersive media. It turns out that the IBSEF method and the NAE are powerful and might be used in further works to find novel solutions for other types of nonlinear evolution equations ascending in physical sciences and engineering.

Keywords: improved Bernoulli subequation function method; the new auxiliary equation method; the variant Boussinesq equation; exact solutions; soliton

1 Introduction

Since many processes in science, technology, and engineering are modeled through nonlinear evolution equations (NLEEs), the investigation of closed-form analytical solutions of NLEEs is very important. Closed-form solutions provide further physical information and help us to understand the processes of the associated physical systems. Consequently, their studies are of fundamental importance due to the effective application of the analytical solutions in various fields, such as plasma physics, solid-state physics, neural physics, chaos, diffusion process, reaction process, optical fibers, nonlinear optics, quantum mechanics, mathematical biology, propagation of shallow water waves, and electromagnetic theory [1,2,3,4,5,6,7,8,9,10,11]. On account of this, researchers established several techniques and tried to examine various types of NLEEs with the help of computer algebra like Maple, MATLAB, and Mathematica. However, not all models are solvable by a single method. Consequently, several methods [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] have been established by mathematicians, physicists, and engineers.

Among the approaches, the improved Bernoulli subequation function (IBSEF) method [32,33] and the new auxiliary equation (NAE) [34] method are effective, direct, and compatible algebraic methods for establishing exact soliton solutions to NLEEs. In 1871, the Boussinesq equation was derived to describe certain physical processes. Since then, several generalizations and some variants of this model have been established in the literature. To the best of our understanding, through different schemes, several researchers have investigated the variant Boussinesq equation, namely: Gao and Tian [35] examined it through the generalized tanh method; Yan and Zhang [36] used an improved sine–cosine method and Wu elimination method; Jabbari et al. [37] investigated it by the homotopy analysis method; Ul-Hassan [38] applied exp-function method; Manafian et al. [39] implemented the improved tanh( ϕ ( ξ ) / 2 ) method and the improved ( G ′ / G ) -expansion method; Osman et al. [40] searched based on the generalized unified method; and Alharbi and Almatrafi [41] exploited the improved exp( − ϕ ( η ) ) -expansion method. Mohapatra and Guedes Soares [42] used the hyperbolic tangent method to investigate the solitary wave solutions of the one-dimensional coupled nonlinear Boussinesq equations related to shallow water. They also use the first integral approach to search for exact solutions to the generalized nonlinear Boussinesq equations. Mohapatra et al. [43] studied the long nonlinear internal waves between analytical and numerical wave models in shallow water. The (G′/G)-expansion procedure was used to achieve the exact and explicit solitary wave solutions to the model in the presence of dispersive coefficient.

Nevertheless, the stated model has not been investigated through the IBSEF method and the NAE approach. Therefore, the purpose of this article is to accomplish broad-ranging and adequate standard soliton solutions to the variant Boussinesq equations by putting into use the proposed methods. Besides, we analyze various types of waves for different values of the free parameters of the obtained solutions illustrated in the 3D plot via MATLAB and highlight the significance of wave number, velocity, and speed of the solutions in changing the nature of the wave profile.

2 Description of two proposed methods

This section introduces and analyzes the IBSEF method and the NAE method in detail.

2.1 The IBSEF method

To illustrate the IBSEF method, we take into consideration the NLEE associated with two independent variables, x and t , of the form [39]:

(2.1.1) R u , ∂ u ∂ t , ∂ u ∂ x , ∂ 2 u ∂ x ∂ t , ∂ 2 u ∂ t 2 , ∂ 2 u ∂ x 2 , … = 0 ,

where u = u ( x , t ) is the elevation of waves; R is a polynomial of the wave variable u ( x , t ) and ∂ u ∂ x , ∂ u ∂ t , … ; t is the temporal variable; and x is the spatial variable.

Under the traveling wave variable

(2.1.2) u ( x , t ) = u ( η ) and η = α x − ct ,

where η is the traveling wave variable, α is the wave number, and c is the speed of the traveling wave, the NLEE (2.1.1) is reduced by the action of the wave variable to the subsequent nonlinear differential equation:

(2.1.3) L u , du d η , d 2 u d η 2 , … = 0 ,

where L is a polynomial of the variable u ( η ) and its derivatives. In accordance with the IBSEF method, the solution of Eq. (2.1.3) can be formulated as

(2.1.4) u ( η ) = ∑ i = 0 q a i H i ∑ j = 0 p b j H j = a 0 + a 1 H + a 2 H 2 + … + a q H q b 0 + b 1 H + b 2 H 2 + b 3 H 3 + … + b p H p ,

where a i ( i = 0 , 1 , 2 , 3 , … , q ) and b j ( j = 0 , 1 , 2 , 3 , … , p ) are indefinite constants to be determined and H = H ( η ) satisfies the improved Bernoulli equation

(2.1.5) H ′ = σ H + d H M , σ ≠ 0 , d ≠ 0 , M ∈ R − { 0 , 1 , 2 } .

The values of p and q of undefined parameters can be determined by considering the balancing principle of the highest order linear term with the maximum order nonlinear term. This technique provides the values of q , p , and M . Introducing solution (2.1.4) along with the improved Bernoulli equation into Eq. (2.1.3) yields a polynomial equation Ω ( H ) of H :

(2.1.6) Ω ( H ) = ρ s H s + … + ρ 1 H + ρ 0 = 0 .

Equalizing the coefficients of Ω ( H ) to zero will provide a system of algebraic equations:

(2.1.7) ρ k = 0 , k = 0 , …

With Mathematica software, we can unravel the system of algebraic equations to determine the values of a 0 , a 1 , … , a q and b 0 , b 1 , … , b p . The solutions of the improved nonlinear Bernoulli equation depend on the values of the parameters σ and d . The two types of solutions to the improved Bernoulli’s equation are as follows:

(2.1.8) H ( η ) = − d σ + τ e σ ( M − 1 ) η 1 1 − M , σ ≠ d ,

(2.1.9) H ( η ) = ( τ − 1 ) + ( τ + 1 ) tanh σ ( 1 − M ) η 2 1 − tanh σ ( 1 − M ) η 2 1 1 − M , σ = d , τ ∈ R .

Thus, by embedding the values of the parameters a i , ( i = 0 , 1 , 2 , … , p ) , b j , j = ( 0 , 1 , 2 , … , q ) , and the solutions of the improved Bernoulli equation (2.1.5) together with wave variable, further general and some acknowledged solutions accessible in the literature for the definite values of the subjective parameters to the NLEE (2.1.1) can be ascertained.

2.2 The NAE method

Suppose the NLEE [38]

(2.2.1) L ( q , q t , q x , q y , q z , q xx , … … … ) = 0 ,

where L is a function of the (3 + 1)-dimensional wave function q ( x , y , z , t ) and its derivatives.

Step 1: To format into an ordinary differential equation of partial differential equation (2.2.1), we need to choose a wave variable as

(2.2.2) q ( x , y , z , t ) = Q ( μ ) , μ = α x + β y + γ z − δ t ,

where α , β , and γ are wave numbers and δ is the wave speed. The wave variable transformed Eq. (2.2.1) into a nonlinear equation as follows:

(2.2.3) M ( Q , Q ′ , Q ″ , Q ″ ′ , … ) = 0 ,

where prime means the derivative with respect to μ .

Step 2: In harmony with the NAE method, the exact soliton solution of Eq. (2.2.3) is assigned to be

(2.2.4) Q ( μ ) = ∑ k = 0 L S k a kg ( μ ) ,

where S 0 , S 1 , S 2 , … , S L are unknown parameters to be calculated, wherein S L ≠ 0 , and g ( μ ) is the solution of the nonlinear equation

(2.2.5) g ′ ( μ ) ln a = { p a − g ( μ ) + q + r a g ( μ ) } .

Step 3: The balancing principle needs to be applied to find the value of the positive integer L in Eq. (2.2.4).

Step 4: Eq. (2.2.3) provides a polynomial of a kg ( μ ) , k = 0, 1, 2 , … , by substituting the solution (2.2.4) together with (2.2.5). Furthermore, we attain a system of algebraic equations by conveying each coefficient of the resulting polynomials to zero.

The values of S k , p , q , and r obtained from Step 4 and the solution g ( μ ) , which is attained in (2.2.5) by making use of various cases in Eq. (2.2.4), provide ample soliton solution to the evolution Eq. (2.2.1).

3 Solution analysis

In this section, we have established some standard and broad-ranging explicit wave solutions through the aforementioned methods of the variant Boussinesq equation, which has the following form [39]:

(3.1) u t + H x + u u x = 0 , H t + ( uH ) x + u xxx = 0 .

To convert the nonlinear model (3.1), we used the wave transformation as follows:

(3.2) u ( x , t ) = u ( η ) and H ( x , t ) = H ( η ) , where η = λ x − ct .

Using Eq. (3.2) into Eq. (3.1) and integrating, we obtain

(3.3) − cu + λ H + λ u 2 2 = 0 ,

(3.4) − cH + λ ( uH ) + λ 3 u ″ = 0 .

Eq. (3.3) gives

(3.5) H = cu λ − u 2 2 .

By replacing the value of H in (3.4), we accomplish the subsequent nonlinear equation as follows:

(3.6) 2 λ 4 u ″ − 2 c 2 u + 3 c λ u 2 − λ 2 u 3 = 0 .

3.1 Analysis of solutions through the IBSEF method

To apply the procedure of the IBSEF method, we gain the correspondence for q, p, and M using the balancing principle between u″ and u 3 in Eq. (3.6) as follows:

p + M = q + 1 .

Picking p = 1 and M = 3 provides q = 3 . Thus, the trial solution to Eq. (2.1.4) becomes

(3.1.1) u ( η ) = a 0 + a 1 H ( η ) + a 2 H 2 ( η ) + a 3 H 3 ( η ) b 0 + b 1 H ( η ) ,

where H ( ξ ) is the solution to the improved Bernoulli Eq. (2.1.5).

The substitution of Eq. (3.1.1) together with Eq. (2.1.5) in Eq. (3.6) yields a polynomial in H . Solving the system of the algebraic equations, we secure the subsequent set of solutions of the parameters as follows:

(3.1.2) Set 1 : c = ± 2 λ 2 σ , a 0 = ± 4 λ σ b 0 , a 1 = ± 4 λ σ b 1 , a 2 = ± 4 λ d b 0 , a 3 = ± 4 λ d b 1 ,

(3.1.3) Set 2 : c = ∓ 2 λ 2 σ , a 0 = 0 , a 1 = 0 , a 2 = ± 4 λ d b 0 , a 3 = ± 4 λ d b 1 .

For σ ≠ d ,

Case 1(a): By making use of the values stated in (3.1.2) along with (2.1.8) in Eq. (3.1.1), the subsequent exponential function solution to the variant Boussinesq equation is established as follows:

(3.1.4) u 1 ( x , t ) = ∓ 4 λ σ 2 τ de 2 σ η − σ τ ,

where σ ≠ d and η = ( λ x − ct ) .

Rewrite the solution function (3.1.4) into hyperbolic form

(3.1.5) u 1 0 ( x , t ) = ± 4 λ σ 2 τ ( sinh ψ − cosh ψ ) ( d − τ σ ) cosh ψ + ( d + τ σ ) sinh ψ ,

where ψ = σ ( λ x − ct ) .

Let us choose d = 2 σ τ , then the solution (3.1.5) becomes the rational function solution of the hyperbolic function as follows:

(3.1.6) u 1 1 ( x , t ) = ± 4 λ σ ( sinh ψ − cosh ψ ) cosh ψ + 3 sinh ψ ,

and we obtain a part of the solution of (3.1) by putting u 1 1 ( x , t ) in Eq. (3.5), which provides

(3.1.7) H 1 1 ( x , t ) = ± 8 σ ( c − 2 λ 2 σ ) ( cosh 2 ψ − sinh ψ cosh ψ ) + λ 2 σ − 3 c 2 ( 10 cosh 2 ψ + 6 sinh ψ cosh ψ − 9 ) .

For d = σ τ , we bring out the hyperbolic function solution to Eq. (3.1) of the mode:

(3.1.8) u 1 2 ( x , t ) = ± 2 λ σ ( 1 − coth ψ ) ,

and we undertake from (3.5)

(3.1.9) H 1 2 ( x , t ) = ± 2 σ ( 1 − coth ψ ) ( − σ ( 1 − coth ψ ) λ 2 + c ) ,

where ψ = σ ( λ x − ct ) .

On the other hand, by setting d = − σ τ in (3.1.5), we accomplish the subsequent solution of the hyperbolic function:

(3.1.10) u 1 3 ( x , t ) = ± 2 λ σ ( 1 − tanh ψ ) ,

and based on Eq. (3.5), we obtain

(3.1.11) H 1 3 ( x , t ) = ± 2 σ ( 1 − tanh ψ ) ( − σ ( 1 − tanh ψ ) λ 2 + c ) ,

where ψ = σ ( λ x − ct ) .

Case 2(a): Considering the value of the parameter arranged in (3.1.3) with the help of Eq. (2.1.8) in Eq. (3.1.1), we accomplish the following solution:

(3.1.12) u 2 ( x , t ) = ∓ 4 σ d λ ( e 2 σ η σ τ + e 4 σ η d ) ( d e 2 σ η − σ τ ) 2 ,

which can be converted to hyperbolic form as follows:

(3.1.13) u 2 0 ( x , t ) = ∓ 4 d σ λ ( σ τ + d ( cosh 2 ψ + sinh 2 ψ ) ) { ( d − σ τ ) cosh ψ + ( d + σ τ ) sinh ψ } 2 ,

where ψ = σ ( λ x − ct ) .

Setting d = 2 σ τ in Eq. (3.1.13) represents the solution of (3.1) as follows:

(3.1.14) u 2 1 ( x , t ) = ∓ 8 σ λ ( 1 + 2 ( cosh 2 ψ + sinh 2 ψ ) ) ( cosh ψ + 2 sinh ψ ) 2 ,

and on behalf of Eq. (3.5),

H 2 1 ( x , t ) = ∓ 8 c σ ( 1 + 2 sinh ( 2 ψ ) + 2 cosh ( 2 ψ ) ) ( cosh ψ + 2 sinh ψ ) 2 (3.1.17) − 32 λ 2 σ 2 ( 1 + 2 sinh ( 2 ψ ) + 2 cosh ( 2 ψ ) ) 2 ( cosh ψ + 2 sinh ψ ) 4 .

Choose d = σ τ , then Eq. (3.1.13) reshapes as follows:

(3.1.16) u 2 2 ( x , t ) = ∓ 2 σ λ ( coth ψ + coth 2 ψ )

and Eq. (3.5) yields

(3.1.17) H 2 2 ( x , t ) = ∓ 2 σ coth ψ ( coth ψ + 1 ) ( − λ 2 σ ( coth 2 ψ + coth ψ ) + c )

When d = − σ τ , the solution function (3.1.13) can be rewritten as:

(3.1.18) u 2 3 ( x , t ) = ∓ 2 σ λ ( tanh 2 ψ − tanh ψ )

and from Eq. (3.5),

(3.1.19) H 2 3 ( x , t ) = ∓ 2 σ tanh ψ ( tanh ψ − 1 ) × ( − λ 2 σ ( tanh 2 ψ − tanh ψ ) + c ) ,

where ψ = σ ( λ x − ct ) .

For σ = d ,

Case 1(b): By substituting the values of the parameters arranged in (3.1.2) along with Eqs. (2.1.2) and (2.1.9) in Eq. (3.1.1), we obtain the exponential function solution to the variant Boussinesq equation of the exponential form:

(3.1.20) u 3 ( x , t ) = ± 4 σ λ τ ( tanh ψ − 1 ) ( τ + 1 ) tanh ψ − ( τ − 1 ) ,

where ψ = σ ( λ x − ct ) .

Eq. (3.1.20) can be written as

(3.1.21) u 3 0 ( x , t ) = ± 4 σ λ τ ( sinh ψ − cosh ψ ) ( τ + 1 ) sinh ψ − ( τ − 1 ) cosh ψ .

Choosing τ = 2 , the solution (3.1.21) yields the solution of (3.1) as follows:

(3.1.22) u 3 1 ( x , t ) = ± 4 σ λ ( sinh ψ − cosh ψ ) 2 sinh ψ − cosh ψ ,

and another solution of Eq. (3.1) obtained by captivating (3.5) is

(3.1.23) H 3 1 ( x , t ) = ± 12 σ − 4 λ 2 σ 3 + c ( cosh 2 ψ − sinh ψ cosh ψ ) + 2 λ 2 σ 3 − 2 c 3 5 cosh 2 ψ − 4 cosh ψ sinh ψ − 4 .

If we pick τ = 1 , then Eq. (3.1.15) becomes

(3.1.24) u 3 2 ( x , t ) = ± 2 σ λ ( 1 − coth ψ ) ,

and Eq. (3.5) gives

(3.1.25) H 3 2 ( x , t ) = ± 2 σ ( 1 − coth ψ ) ( − λ ( 1 − coth ψ ) σ λ + c ) .

When τ = − 1 , then Eq. (3.1.15) becomes

(3.1.26) u 3 3 ( x , t ) = ± 2 σ λ ( 1 − tanh ψ ) ,

and using Eq. (3.5) gives

(3.1.27) H 3 3 ( x , t ) = ± 2 σ ( 1 − tanh ψ ) ( − λ ( 1 − tanh ψ ) σ λ + c ) ,

where ψ = σ ( λ x − ct ) .

Case 2(b): By using Eq. (3.1.3) with the help of Eqs. (2.1.2) and (2.1.9) from Eq. (3.1.1), we accomplish the following solution:

u 4 ( x , t ) = ∓ 4 σ λ ( 1 + tanh ψ ) ( τ + 1 ) tanh ψ − ( τ − 1 ) ,

where ψ = σ ( λ x − ct ) .

Calculating the result u 4 ( x , t ) , we then find

(3.1.28) u 4 0 ( x , t ) = ∓ 4 σ λ ( cosh ψ + sinh ψ ) ( τ + 1 ) sinh ψ − ( τ − 1 ) cosh ψ .

When τ = 1 , u 4 0 ( x , t ) can be rewritten as

(3.1.29) u 4 1 ( x , t ) = ∓ 2 σ λ ( coth ψ + 1 ) ,

and

(3.1.30) H 4 1 ( x , t ) = ∓ 2 σ ( 1 + coth ψ ) ( − λ ( 1 + coth ψ ) σ λ + c ) .

Applying the value τ = − 1 , the solution u 4 0 ( x , t ) turns out to be

(3.1.31) u 4 2 ( x , t ) = ∓ 2 σ λ ( tanh ψ + 1 ) ,

and

(3.1.32) H 4 1 ( x , t ) = ∓ 2 σ ( 1 + tanh ψ ) ( − λ ( 1 + tanh ψ ) σ λ + c ) ,

where ψ = σ ( λ x − ct ) .

It is important to note that the wave solutions of the variant Boussinesq equation found here are functional and useful and were not proven in the earlier research. The solutions derived earlier might be fruitful in investigating unidirectional wave propagation in nonlinear media and dispersive relativistic one-particle theory, etc.

3.2 Analysis of solutions through the NAE method

To find the solutions to the stated equation through the NAE method, we gain the value of L = 1 based on the balancing principle between u ′ ′ and u 3 in Eq. (3.6).

Using the NAE approach, substitute the value of L into Eq. (2.2.4), which is the auxiliary solution of Eq. (3.6) of the form:

(3.2.1) u ( x , t ) = a 0 + a 1 a g ( η ) ,

where g ( η ) is the root of the nonlinear Eq. (2.2.5).

Based on the solution (3.2.1) with help of (2.2.5) from Eq. (3.6), we assert

(3.2.2) ( 4 λ 4 n 2 A 1 − λ 2 A 1 3 ) { a g ( η ) } 3 + ( 6 λ 4 mn A 1 − 3 λ 2 A 0 A 1 2 + 3 c λ A 1 2 ) { a g ( η ) } 2 + ( 4 l λ 4 n A 1 + 2 λ 4 m 2 A 1 − 3 λ 2 A 0 2 A 1 + 6 c λ A 0 A 1 − 2 c 2 A 1 ) { a g ( η ) } + 2 l λ 4 m A 1 − λ 2 A 0 3 + 3 c λ A 0 2 − 2 c 2 A 0 = 0 .

Equalizing the like power of a g ( η ) from Eq. (3.2.2) yields a system of algebraic equations, which when solved provides the following results:

Set 1:

(3.2.3) c = ± λ 2 ( 2 ln m 2 − 4 ln − m 2 m 2 − 4 ln + 4 lmn − m 3 2 ln − m 2 − m m 2 − 4 ln , A 0 = ± ( m + m 2 − 4 ln ) λ , A 1 = 2 λ n .

Set 2:

(3.2.4) c = ± λ 2 ( 2 ln m 2 − 4 ln − m 2 m 2 − 4 ln − 4 lmn + m 3 m m 2 − 4 ln + 2 ln − m 2 , A 0 = ± ( − m + m 2 − 4 ln ) λ , A 1 = 2 λ n .

Based on the above values of the parameters in Eq. (3.2.4), as follows:

When m 2 − 4 ln < 0 and n ≠ 0 .

By substituting the values of the constants arranged in (3.2.3) and (2.2.5) into Eq. (3.2.1), we obtain the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.5) u 5 1 ( x , t ) = λ 4 l n − m 2 tan 4 l n − m 2 η 2 + − 4 l n + m 2 ,

and

(3.2.6) H 5 1 ( x , t ) = θ λ 2 tan 2 − θ η 2 2 + − θ ( − λ 2 θ + c ) × tan − θ η 2 + c θ − λ 2 θ 2 ,

where η = λ x − ct and θ = m 2 − 4 ln ,

(3.2.7) u 5 2 ( x , t ) = λ − 4 l n − m 2 cot 4 l n − m 2 η 2 + − 4 l n + m 2 ,

and

(3.2.8) H 5 2 ( x , t ) = θ λ 2 cot 2 − θ η 2 2 − − θ ( − λ 2 θ + c ) × cot − θ η 2 + c θ − λ 2 θ 2 ,

where η = λ x − ct and θ = m 2 − 4 ln .

By substituting the values of the constant arranged in (3.2.4) and (2.2.5) into Eq. (3.2.1), we obtain the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.9) u 6 1 ( x , t ) = λ − 4 l n − m 2 tan 4 l n − m 2 η 2 + − 4 l n + m 2 ,

and

(3.2.10) H 6 1 ( x , t ) = θ λ 2 tan 2 − θ η 2 2 − − θ ( − λ 2 θ + c ) × tan − θ η 2 + c θ − λ 2 θ 2 ,

(3.2.11) u 6 2 ( x , t ) = λ 4 l n − m 2 cot 4 l n − m 2 η 2 + − 4 l n + m 2 ,

and

(3.2.12) H 6 2 ( x , t ) = θ λ 2 cot 2 − θ η 2 2 + − θ ( − λ 2 θ + c ) × cot − θ η 2 + c θ − λ 2 θ 2 ,

where θ = m 2 − 4 ln and η = λ x − ct .

When m 2 − 4 ln > 0 and n ≠ 0 .

By combining (3.2.3) and (2.2.5) in Eq. (3.2.1), we accomplish the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.13) u 7 1 ( x , t ) = λ − 4 l n + m 2 1 − tanh − 4 l n + m 2 η 2 ,

and

(3.2.14) H 7 1 ( x , t ) = 1 − tanh θ η 2 θ λ 2 2 tanh θ η 2 + c θ − λ 2 θ 2 ,

(3.2.15) u 7 2 ( x , t ) = λ − 4 l n + m 2 1 − coth − 4 l n + m 2 η 2 ,

and

(3.2.16) H 7 2 ( x , t ) = 1 − coth θ η 2 θ λ 2 2 coth θ η 2 + c θ − λ 2 θ 2 ,

where θ = m 2 − 4 ln and η = λ x − ct .

Using (3.2.4) and (2.2.5) into Eq. (3.2.1), we achieve the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.17) u 8 1 ( x , t ) = λ − 4 l n + m 2 1 + tanh − 4 l n + m 2 η 2 ,

and

(3.2.18) H 8 1 ( x , t ) = 1 + tanh θ η 2 × − θ λ 2 2 tanh θ η 2 + c θ − λ 2 θ 2 ,

(3.2.19) u 8 2 x , t ) = λ − 4 l n + m 2 1 + coth − 4 l n + m 2 η 2 ,

and

(3.2.20) H 8 2 ( x , t ) = 1 + coth θ η 2 × − θ λ 2 2 coth θ η 2 + c θ − λ 2 θ 2 ,

where θ = m 2 − 4 ln and η = λ x − ct .

When m 2 + 4 l 2 < 0 , n ≠ 0 , and n = − l .

Employing (3.2.3) together with (2.2.5) into Eq. (3.2.1), we achieve the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.21) u 9 1 ( x , t ) = λ − 4 l 2 − m 2 tan − 4 l 2 − m 2 η 2 + 4 l 2 + m 2 ,

and

(3.2.22) H 9 1 ( x , t ) = λ 2 ϕ 2 tan 2 − ϕ η 2 − 1 + − ϕ ( c − ϕ ) tan − ϕ η 2 + c ϕ ,

(3.2.23) u 9 2 ( x , t ) = λ − − 4 l 2 − m 2 cot − 4 l 2 − m 2 η 2 + 4 l 2 + m 2 ,

and

(3.2.24) H 9 2 ( x , t ) = λ 2 ϕ 2 cot 2 − ϕ η 2 − 1 − − ϕ ( c − ϕ ) cot − ϕ η 2 + c ϕ ,

where ϕ = 4 l 2 + m 2 and η = λ x − ct .

Considering (3.2.4) together with (2.2.5) into Eq. (3.2.1), we attain the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) in the subsequent form

(3.2.25) u 10 1 ( x , t ) = λ − − 4 l 2 − m 2 tan − 4 l 2 − m 2 η 2 + 4 l 2 + m 2 ,

and

(3.2.26) H 10 1 ( x , t ) = λ 2 ϕ 2 tan 2 − ϕ η 2 − 1 − − ϕ ( c − ϕ ) tan − ϕ η 2 + c ϕ ,

(3.2.27) u 10 2 ( x , t ) = λ − 4 l 2 − m 2 cot − 4 l 2 − m 2 η 2 + 4 l 2 + m 2 ,

and

(3.2.28) H 10 2 ( x , t ) = λ 2 ϕ 2 cot 2 − ϕ η 2 − 1 + − ϕ ( c − ϕ ) cot − ϕ η 2 + c ϕ ,

where ϕ = 4 l 2 + m 2 and η = λ x − ct .

When m 2 + 4 l 2 > 0 , n ≠ 0 , and n = − l .

Inserting (3.2.3) together with (2.2.5) into Eq. (3.2.1), we secure the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.29) u 11 1 ( x , t ) = λ 4 l 2 + m 2 1 − tanh 4 l 2 + m 2 η 2 ,

and

(3.2.30) H 11 1 ( x , t ) = c ϕ 1 − tanh ϕ η 2 − λ 2 ϕ 2 1 − tanh ϕ η 2 2 ,

(3.2.31) u 11 2 ( x , t ) = λ 4 l 2 + m 2 1 − coth 4 l 2 + m 2 η 2 ,

and

(3.2.32) H 11 2 ( x , t ) = c ϕ 1 − coth ϕ η 2 − λ 2 ϕ 2 1 − coth ϕ η 2 2 ,

where ϕ = 4 l 2 + m 2 and η = λ x − ct .

By involving (3.2.4) along with (2.2.5) into Eq. (3.2.1), we derive the solution of stated Eqs. (3.1)–(3.2.1) and (3.5) as follows:

(3.2.33) u 12 1 ( x , t ) = λ 4 l 2 + m 2 1 + tanh 4 l 2 + m 2 η 2 ,

and

(3.2.34) H 12 1 ( x , t ) = c ϕ 1 + tanh ϕ η 2 − λ 2 ϕ 2 1 + tanh ϕ η 2 2 ,

(3.2.35) u 12 2 ( x , t ) = λ 4 l 2 + m 2 1 + coth 4 l 2 + m 2 η 2 ,

and

(3.2.36) H 12 2 ( x , t ) = c ϕ 1 + coth ϕ η 2 − λ 2 ϕ 2 1 + coth ϕ η 2 2 ,

where ϕ = m 2 + 4 l 2 and η = λ x − ct .

When m 2 − 4 l 2 < 0 and n = l .

Appling (3.2.3) along with (2.2.5) into Eqs. (3.2.1) and (3.5) gives

(3.2.37) u 13 1 ( x , t ) = λ 4 l 2 − m 2 tan 4 l 2 − m 2 η 2 + − 4 l 2 + m 2 ,

and

(3.2.38) H 13 1 ( x , t ) = c ω 1 + i tan − ω η 2 − λ 2 ω 2 1 − tan − ω η 2 ,

(3.2.39) u 13 2 ( x , t ) = λ − 4 l 2 − m 2 cot 4 l 2 − m 2 η 2 + − 4 l 2 + m 2 ,

and

(3.2.40) H 13 2 ( x , t ) = c ω 1 − i cot − ω η 2 − λ 2 ω 2 1 + cot − ω η 2 .

Substituting (3.2.4) along with (2.2.5) into Eqs. (3.2.1) and (3.5) gives

(3.2.41) u 14 1 ( x , t ) = λ − 4 l 2 − m 2 tan 4 l 2 − m 2 η 2 + − 4 l 2 + m 2 ,

and

(3.2.42) H 14 1 ( x , t ) = c ω 1 − i tan − ω η 2 − λ 2 ω 2 1 + tan − ω η 2 ,

(3.2.43) u 14 2 ( x , t ) = λ 4 l 2 − m 2 cot 4 l 2 − m 2 η 2 + − 4 l 2 + m 2 ,

and

(3.2.44) H 14 2 ( x , t ) = c ω 1 + i cot − ω η 2 − λ 2 ω 2 1 − cot − ω η 2 ,

where m 2 − 4 l 2 = ω and η = λ x − ct .

When m 2 − 4 l 2 > 0 , n ≠ 0 and n = l .

Setting (3.2.3) together with (2.2.5) into Eq. (3.2.1), we derive the solutions of (3.1) from Eqs. (3.2.1) and (3.5), which provides

(3.2.45) u 15 1 ( x , t ) = λ − 4 l 2 + m 2 1 − tanh − 4 l 2 + m 2 η 2 ,

and

(3.2.46) H 15 1 ( x , t ) = 1 − tanh ω η 2 × λ 2 ω 2 tanh ω η 2 − 1 + c ω ,

(3.2.47) u 15 2 ( x , t ) = λ − 4 l 2 + m 2 1 − coth − 4 l 2 + m 2 η 2 ,

and

(3.2.48) H 15 2 ( x , t ) = 1 − coth ω η 2 × λ 2 ω 2 coth ω η 2 − 1 + c ω ,

By means of (3.2.4) together with (2.2.5) from Eq. (3.2.1), we derive the solutions of (3.1) in the subsequent form

(3.2.49) u 16 1 ( x , t ) = λ − 4 l 2 + m 2 1 + tanh − 4 l 2 + m 2 η 2 ,

and

(3.2.50) H 16 1 ( x , t ) = 1 + tanh ω η 2 × c ω − λ 2 ω 2 tanh ω η 2 + 1 ,

(3.2.51) u 16 2 ( x , t ) = λ − 4 l 2 + m 2 1 + coth − 4 l 2 + m 2 η 2 ,

and

(3.2.52) H 16 2 ( x , t ) = 1 − coth ω η 2 × c ω − λ 2 ω 2 coth ω η 2 + 1 ,

where m 2 − 4 l 2 = ω and η = λ x − ct . When

m 2 = 4 ln .

By means of (3.2.3) together with (2.2.5) into Eq. (3.2.1), we gain the solutions of (3.1) from Eqs. (3.2.1) and (3.5), which gives

(3.2.53) u 17 1 ( x , t ) = − 2 λ η ,

and

(3.2.54) H 17 1 ( x , t ) = − 2 c η − 2 λ 2 η 2 .

Using (3.2.4) along with (2.2.5) into Eq. (3.2.1), we derive the solutions of (3.1) from Eqs. (3.2.1) and (3.5), which provides

(3.2.55) u 18 1 ( x , t ) = 2 λ η ,

and

(3.2.56) H 18 1 ( x , t ) = 2 c η − 2 λ 2 η 2 ,

where η = λ x − ct .

When m = 0 and l = − n , we secure the solutions of (3.1) by putting (3.2.3) and (2.2.5) into Eqs. (3.2.1) and (3.5) as follows:

(3.2.57) u 19 1 ( x , t ) = 2 ( 1 + e − 2 n η ) λ n ( − 1 + e − 2 n η ) + 4 n 2 λ , = 2 n λ ( 1 − coth ( n η ) ) ,

and

(3.2.58) H 19 1 ( x , t ) = 2 n ( 1 − coth ( n η ) ) ( c − λ 2 n ( 1 − coth ( n η ) ) .

Providing (3.2.4) along with (2.2.5) into Eq. (3.2.1), we attain the solution functions (3.2.1) and (3.5) as shown below

(3.2.59) u 20 1 = − 2 ( 1 + e − 2 n η ) λ n ( − 1 + e − 2 n η ) + 4 n 2 λ = 2 n λ ( coth ( n η ) − 1 ) ,

and

(3.2.60) H 20 1 ( x , t ) = 2 n ( coth ( n η ) − 1 ) ( c − λ 2 n ( coth ( n η ) − 1 ) ,

where η = λ x − ct .

When l = m = K , n = 0 .

Applying the parametric values stated in (3.2.3) into Eqs. (3.2.1) and (3.5) represent the solutions

(3.2.61) u 21 1 ( x , t ) = 2 K λ ,

and

(3.2.62) H 21 1 ( x , t ) = − 2 K 2 λ 2 + 2 c K .

When m = − ( l + n ) .

Substituting the value assigned in (3.2.3) and Eq. (2.2.5) into Eq. (3.2.1) yields the solution

(3.2.63) u 24 1 ( x , t ) = 2 ( l − e ( l − n ) η ) λ n ( n − e ( ( l − n ) η ) ) − 2 l λ .

Using the value u 24 1 ( x , t ) in Eq. (3.5), we secure

(3.2.64) H 24 1 ( x , t ) = 2 ( l − n ) e ( l − n ) η ( c n + ( − l λ 2 + λ 2 n − c ) e ( l − n ) η ) ( n − e ( l − n ) η ) 2 ,

where η = λ x − ct .

Applying (3.2.4) and Eq. (2.2.5) to Eqs. (3.2.1) and (3.5) gives

(3.2.65) u 25 1 ( x , t ) = − 2 ( l − e ( l − n ) η ) λ n ( n − e ( ( l − n ) η ) ) + 2 l λ ,

and

(3.2.66) H 25 1 ( x , t ) = − 2 ( l − n ) e ( l − n ) η ( c n + ( l λ 2 − λ 2 n − c ) e ( l − n ) η ) ( n − e ( l − n ) η ) 2 ,

respectively, where η = λ x − ct .

When m = l + n .

Selecting (3.2.3) and (2.2.5) for Eq. (3.2.1) and after simplification, we obtain

(3.2.67) u 22 ( x , t ) = − 2 λ n ( 1 − l ) cosh ( l − n ) η 2 − ( 1 + l ) sinh ( l − n ) η 2 ( 1 − n ) cosh ( l − n ) η 2 − ( 1 + n ) sinh ( l − n ) η 2 + 2 l λ

Choosing n = 1 in u 22 ( x , t ) gives

(3.2.68) u 22 1 ( x , t ) = λ ( 1 − l ) coth ( l − 1 ) η 2 − ( 1 + l ) + 2 l λ ,

and from Eq. (3.5),

(3.2.69) H 22 1 ( x , t ) = λ 2 2 ( l − 1 ) 2 1 − coth ( l − 1 ) η 2 × coth ( l − 1 ) η 2 − 1 + c λ 2 ( l − 1 ) ,

where η = λ x − ct .

If n = − 1 , u 22 ( x , t ) pays

(3.2.70) u 22 2 ( x , t ) = λ ( 1 − l ) − ( 1 + l ) tanh ( l + 1 ) η 2 + 2 l λ ,

and (3.5) gives

(3.2.71) H 22 2 ( x , t ) = λ 2 2 ( l + 1 ) 2 1 − tanh ( l + 1 ) η 2 × tanh ( l + 1 ) η 2 − 1 + c λ 2 ( l + 1 ) ,

where η = λ x − ct .

Based on (3.2.4) with the help of (2.2.5), after computation, the solution function (3.2.1) becomes

(3.2.72) u 23 ( x , t ) = 2 λ n ( 1 − l ) cosh ( l − n ) η 2 − ( 1 + l ) sinh ( l − n ) η 2 ( 1 − n ) cosh ( l − n ) η 2 − ( 1 + n ) sinh ( l − n ) η 2 − 2 l λ .

For n = 1 , the solution function u 23 ( x , t ) provides

(3.2.73) u 23 1 ( x , t ) = − λ ( 1 − l ) coth ( l − 1 ) η 2 − ( 1 + l ) − 2 l λ ,

and from Eq. (3.5), we obtain

(3.2.74) H 23 1 ( x , t ) = λ 2 2 ( l − 1 ) 2 1 − coth ( l − 1 ) η 2 × coth ( l − 1 ) η 2 − 1 − c λ 2 ( l − 1 ) ,

where η = λ x − ct .

Selecting n = − 1 , we obtain from u 23 ( x , t ) ,

(3.2.75) u 23 2 ( x , t ) = − λ ( 1 − l ) − ( 1 + l ) tanh ( l + 1 ) η 2 − 2 l λ ,

and obtain from Eq. (3.5) as

(3.2.76) H 23 2 ( x , t ) = λ 2 2 ( l + 1 ) 2 1 − tanh ( l + 1 ) η 2 × tanh ( l + 1 ) η 2 − 1 − c λ 2 ( l + 1 ) ,

where η = λ x − ct .

When n = m = l ≠ 0 .

We attain the required solution of (3.1) by considering (3.2.3) and (2.2.5) into (3.2.1) and (3.5) as follows:

(3.2.77) u 26 1 ( x , t ) = 3 λ l tan 3 l η 2 − 1 ,

and

(3.2.78) H 26 1 ( x , t ) = − l 2 1 + tanh 3 l η 2 × 3 l λ 2 tanh 3 l η 2 + 3 λ 2 l − 2 3 c ,

where η = λ x − ct .

Again, by selecting the values of (3.2.4) and (2.2.5) and putting into Eqs. (3.2.1) and (3.5), we obtain

(3.2.79) u 27 1 ( x , t ) = − 3 λ l tan 3 l η 2 + 1 ,

and

(3.2.80) H 27 1 ( x , t ) = − l 2 1 + tanh 3 l η 2 × 3 l λ 2 tanh 3 l η 2 + 3 λ 2 l + 2 3 c ,

where η = λ x − ct .

When = l , m = 0 .

Proceeding in this way, for (3.2.3) we preserve

(3.2.81) u 28 1 ( x , t ) = 2 λ l ( tan ( l η ) − 1 ) ,

and

(3.2.82) H 28 1 ( x , t ) = 2 l ( tan ( l η ) − 1 ) ( − l λ 2 tan ( l η ) + l λ 2 + c ) ,

where η = λ x − ct .

Also for (3.2.4), we gain the required solution of (3.1) as follows:

(3.2.83) u 29 1 ( x , t ) = 2 λ l ( − tan ( l η ) − 1 ) ,

and

(3.2.84) H 29 1 ( x , t ) = − 2 l ( tan ( l η ) + 1 ) ( l λ 2 tan ( l η ) + l λ 2 + c ) ,

where η = λ x − ct .

When n = 0 , depending on (3.2.3), the solution functions (3.2.1) and (3.5) provide the results as

(3.2.85) u 30 1 ( x , t ) = 2 m λ ,

and

(3.2.86) H 30 1 ( x , t ) = − 2 m 2 λ 2 + 2 c m .

From the above investigation, we have found various types of solutions of the variant Boussinesq equation, such as hyperbolic function, trigonometric solution, exponential function solution, and rational solution.

4 Result discussions and wave profile description

In this section, the obtained traveling wave solutions of the stated equations are represented in the figures and reviewed the natures of these waves for dissimilar values of the free parameters with the aid of software MATLAB.

4.1 Wave profile interpretations of attained solutions using the IBSEF method

In this article, we have found 16 solutions (whereas Manafian et al. [39] established 15 solutions) using the IBSEF approach, which are illustrative and have rich structured wave profiles that might be helpful in describing the unidirectional wave propagation in nonlinear media with dispersion relativistic one-particle theory. It is important to note that the wave solutions of the variant Boussinesq equation found here are functional and useful and were not proven in the earlier research.

For λ = 0.001 and σ = 1 , the solution u 1 2 ( x , t ) represents a spike-type singular soliton as shown in Figure 1a, where the spike’s amplitude is infinite near the origin and gradually decreases. Increasing the values of the free parameters proportionally, likewise λ = 0.5 and σ = 0.5 , as well as descending the velocity speed c , the same solution is converted to breather soliton with singularity asserted as shown in Figure 1b. For the values λ = 0.5 and σ = 1 , the exact solution delivers a singular soliton presented in Figure 1c. Furthermore, u 1 2 ( x , t ) demonstrates kink shape soliton with a singularity at ( 0,0 ) traced in Figure 1(d) for λ , σ ∈ R − [ − 3 , 3 ] . The other part H 1 2 ( x , t ) of the solution provides a four petal-type wave graphed in Figure 1e, where two spikes are outer and the other two spikes are inner for λ = 0.001 and σ = 1 . By varying the values of the parameters equally, as λ = 25 and σ = 0.2 , the same solution portrays double figures than the previous one that has been outlined in Figure 1(f). Proceeding in this way, by increasing the value of those parameters, we achieve the same result.

$Figure 1 Plot of u 1 2 ( x , t ) {u}_{{1}_{2}}(x,t) and H 1 2 ( x , t ) {H}_{{1}_{2}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = 0.001 and σ = 1 \lambda =0.001\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =1 ; (b) λ = 0.5 and σ = 0.5 \lambda =0.5\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =0.5 ; (c) λ = 0.5 and σ = 1 \lambda =0.5{\ }{\rm{and}}{\ }\sigma =1 ; (d) λ = 4 and σ = 4 \lambda =4\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =4 ; (e) λ = 0.001 and σ = 1 \lambda =0.001\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =1 ; and (f) λ = 0.2 and σ = 0.2 \lambda =0.2\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =0.2 .$

Figure 1

Plot of u 1 2 ( x , t ) and H 1 2 ( x , t ) for various values of λ and σ . (a) λ = 0.001 and σ = 1 ; (b) λ = 0.5 and σ = 0.5 ; (c) λ = 0.5 and σ = 1 ; (d) λ = 4 and σ = 4 ; (e) λ = 0.001 and σ = 1 ; and (f) λ = 0.2 and σ = 0.2 .

The profile of the solution u 1 3 ( x , t ) shows a smooth kink shape soliton for the proportional value of λ and σ as λ , σ ∈ R − ( − 0.7 , 0.7 ) . Furthermore, both solutions u 1 3 ( x , t ) and H 1 3 ( x , t ) represent steep shape kink for − 0.7 ≤ λ and σ ≤ 0.7 , but H 1 3 ( x , t ) is a compacton soliton shown in Figure 2a for λ = 0.8 and σ = 0.8 . Alternatively, for the lower value of λ and/or the upper value of σ , as earlier interval, the same solution displays the kink shape soliton asserted in Figure 2b, since the velocity c depends much on λ than σ .

$Figure 2 Plot of u 1 3 ( x , t ) {u}_{{1}_{3}}(x,t) and H 1 3 ( x , t ) {H}_{{1}_{3}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = 0.8 and σ = 0.8 \lambda =0.8\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =0.8 and (b) λ = 0.5 and σ = 10 \lambda =0.5\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =10 .$

Figure 2

Plot of u 1 3 ( x , t ) and H 1 3 ( x , t ) for various values of λ and σ . (a) λ = 0.8 and σ = 0.8 and (b) λ = 0.5 and σ = 10 .

For very small values of λ and σ , as λ = σ = 0.001 , the exact solution u 2 2 ( x , t ) portrays a single spike-type singular soliton illustrated in Figure 3a. By ascending the values of the free parameters as well as increasing the speed of wave c , the same solution turns out to be breather-type soliton as displayed in Figure 3b for λ = 0.5 and σ = 0.5 . Moreover, for the negative values of σ , the wave just overturns for λ = 0.5 and σ = − 0.5 as shown in Figure 1c. Moreover, by increasing the values gradually, like λ = 3.5 and σ = 3.5 , the same solution changed into ideal kink shape soliton as sketched in Figure 3d.

$Figure 3 Plot of u 2 2 ( x , t ) {u}_{{2}_{2}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = 0.001 and σ = 0.001 \lambda =0.001\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =0.001 ; (b) λ = 0.5 and σ = 0.5 \lambda =0.5{\ }{\rm{and}}{\ }\sigma =0.5 ; (c) λ = 0.5 and σ = − 0.5 \lambda =0.5\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =-0.5 ; and (d) λ = 3.5 and σ = 3.5 . \lambda =3.5{\ }{\rm{and}}{\ }\sigma =3.5.$

Figure 3

Plot of u 2 2 ( x , t ) for various values of λ and σ . (a) λ = 0.001 and σ = 0.001 ; (b) λ = 0.5 and σ = 0.5 ; (c) λ = 0.5 and σ = − 0.5 ; and (d) λ = 3.5 and σ = 3.5 .

The solution u 3 1 ( x , t ) represents an oblique plane for − 0.4 < λ , σ < 0.4 , and when increasing the values of the parameters, it demonstrates singular soliton as traced in Figure 4a–c for λ = σ = 0.5 ; λ = 1 and σ = 0.5 ; and λ = 3 and σ = 3 , respectively. On the other hand, solution H 3 1 ( x , t ) portrays steep soliton as shown in Figure 4d for λ = 0.3 and σ = 0.3 and a special soliton as marked in Figure 4e for λ = 10 and σ = 1 .

$Figure 4 Plot of u 3 1 ( x , t ) {u}_{{3}_{1}}(x,t) and H 3 1 ( x , t ) {H}_{{3}_{1}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = 0.5 and σ = 0.5 \lambda =0.5{\ }{\rm{and}}{\ }\sigma =0.5 ; (b) λ = 1 and σ = 0.5 \lambda =1{\ }{\rm{and}}{\ }\sigma =0.5 ; (c) λ = 3 and σ = 3 \lambda =3\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =3 ; (d) λ = 0.3 and σ = 0.3 \lambda =0.3\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =0.3 ; and (e) λ = 10 and σ = 1 . \lambda =10\hspace{0.25em}{\rm{and}}\hspace{0.25em}\sigma =1.$

Figure 4

Plot of u 3 1 ( x , t ) and H 3 1 ( x , t ) for various values of λ and σ . (a) λ = 0.5 and σ = 0.5 ; (b) λ = 1 and σ = 0.5 ; (c) λ = 3 and σ = 3 ; (d) λ = 0.3 and σ = 0.3 ; and (e) λ = 10 and σ = 1 .

4.2 Wave profile description through the NAE method

Manafian et al. [39] investigated the Boussinesq equation and established only 15 solutions. On the other hand, we have derived 30 solutions in this article, which is twice as many as the solutions developed in ref. [39]. Of these 30 solutions, some of our solutions coincide with those found by Manafian et al., indicating that the solutions we have ascertained are valid, while the others are new.

The exact solutions u 5 1 ( x , t ) and H 5 1 ( x , t ) represent two symmetrical solitons directed in opposite directions as sketched in Figure 5a for λ = σ = 0.5; l = m = n = 1, respectively. Considering larger values of λ and σ as well as speed of the wave c , like λ = σ = l = m = n = 1 , the solutions yield solitons as shown in Figure 5b. For decreasing the negative value of ( m 2 − 4 ln ) and λ = σ = m = 1 and l = n = 3 , the solution u 5 1 ( x , t ) shows irregular periodic soliton and H 5 1 ( x , t ) as portrayed in Figure 5c. On the other hand, another solution u 5 2 ( x , t ) and H 5 2 ( x , t ) display unequal spike-type singular wave asserted in Figure 5d for λ = σ = l = m = n = 1 .

$Figure 5 Plot of u 5 1 ( x , t ) {u}_{{5}_{1}}(x,t) , H 5 1 ( x , t ) {H}_{{5}_{1}}(x,t) , and u 5 2 ( x , t ) {u}_{{5}_{2}}(x,t) , H 5 2 ( x , t ) {H}_{{5}_{2}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = σ = 0 . 5 and l = m = n = 1 \lambda =\sigma =0.5\hspace{0.25em}{\rm{and}}\hspace{0.25em}l=m=n=1 ; (b) λ = σ = l = m = n = 1 \lambda =\sigma =l=m=n=1 ; (c) λ = σ = m = 1 and l = n = 3 \lambda =\sigma =m=1\hspace{0.25em}{\rm{and}}\hspace{0.25em}l=n=3 ; and (d) λ = σ = l = m = n = 1 \lambda =\sigma =l=m=n=1 .$

Figure 5

Plot of u 5 1 ( x , t ) , H 5 1 ( x , t ) , and u 5 2 ( x , t ) , H 5 2 ( x , t ) for various values of λ and σ . (a) λ = σ = 0 . 5 and l = m = n = 1 ; (b) λ = σ = l = m = n = 1 ; (c) λ = σ = m = 1 and l = n = 3 ; and (d) λ = σ = l = m = n = 1 .

For the values of free parameters λ = σ = l = n = 1 and m = 4 , solutions u 7 1 ( x , t ) and H 7 1 ( x , t ) depict standard kink soliton and compacton as illustrated in Figure 6a. By increasing the positive value of m 2 − 4 ln , as λ = σ = l = n = 1 and m = 8 , the same solutions yield kink and smooth kink soliton as marked out in Figure 6b, respectively.

$Figure 6 Plot of u 7 1 ( x , t ) {u}_{{7}_{1}}(x,t) and H 7 1 ( x , t ) {H}_{{7}_{1}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = σ = l = n = 1 and m = 4 \lambda =\sigma =l=n=1{\ }{\rm{and}}{\ }m=4 and (b) λ = σ = l = n = 1 and m = 8 \lambda =\sigma =l=n=1\hspace{0.25em}{\rm{and}}\hspace{0.25em}m=8 .$

Figure 6

Plot of u 7 1 ( x , t ) and H 7 1 ( x , t ) for various values of λ and σ . (a) λ = σ = l = n = 1 and m = 4 and (b) λ = σ = l = n = 1 and m = 8 .

Choosing the infinitesimal value of the parameters, for instance λ = σ = 0.001 , the rational function solution u 17 1 ( x , t ) represents singular soliton with descending spike and H 17 1 ( x , t ) depicts singular soliton with two equal spikes as traced in Figure 7a. Moreover, for the higher value of the velocity c corresponds to λ = σ = 5 , both solutions are representing same shape of wave as portrayed in Figure 7b.

$Figure 7 Plot of u 17 1 ( x , t ) {u}_{{17}_{1}}(x,t) and H 17 1 ( x , t ) {H}_{{17}_{1}}(x,t) for various values of λ \lambda and σ \sigma . (a) λ = σ = 0.001 \lambda =\sigma =0.001 and (b) λ = σ = 5 \lambda =\sigma =5 .$

Figure 7

Plot of u 17 1 ( x , t ) and H 17 1 ( x , t ) for various values of λ and σ . (a) λ = σ = 0.001 and (b) λ = σ = 5 .

From the earlier discussion, it is clear that the obtained solutions of the variant Boussinesq equation through the IBSEF and the NAE methods provide a variety of functional waves, such as kink, singular kink, breather-type, symmetric, compacton, bell-shape soliton, and irregular periodic wave. The wave profiles depend on the positive and negative values of the free parameters and disclose the dynamic behavior of the studied phenomenon. It has been established that wave profiles are strongly interrelated with free parameters, especially with wave velocity c and wave number λ compared to other parameters.

5 Conclusion

In this study, we successfully conducted two significant approaches, namely the IBSEF method and the new auxiliary method to the variant Boussinesq equation, and achieved some standard and advanced solutions and some broad-ranging exact solutions. Setting certain values of the associated parameters, broad-spectrum solutions are able to recover the known solutions and provide the optimal solutions for other values from which the standard functional waves, such as kink, singular kink, breather-type, symmetric, compacton, bell-shape soliton, and irregular periodic wave can be extracted. Furthermore, through some of the documented solutions, we have shown how wave numbers and velocities affect wave profiles. The solutions attained in this study are illustrative and have rich structured wave profiles (Figures 1–7) that might help describe the unidirectional wave propagation in nonlinear media with dispersion relativistic one-particle theory, the long waves on the surface of the liquid layer under gravity action, etc.

adel@sci.cu.edu.eg

Funding information: This research was supported by National Natural Science Foundation of China (No. 71601072), Key Scientific Research Project of Higher Education Institutions in Henan Province of China (No. 20B110006), and the Fundamental Research Funds for the Universities of Henan Province.
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.
Conflict of interest: The authors state no conflict of interest.

References

[1] Kudryashov NA. Exact solutions of the equation for surface waves in a convecting fluid. Appl Math Comput. 2019;344(345):97–106.10.1016/j.amc.2018.10.005Search in Google Scholar

[2] Barman HK, Akbar MA, Osman MS, Nisar KS, Zakarya M, Abdel-Aty AH, et al. Solutions to the Konopelchenko-Dubrovsky equation and the Landau-Ginzburg-Higgs equation via the generalized Kudryashov technique. Results Phys. 2021;24:104092.10.1016/j.rinp.2021.104092Search in Google Scholar

[3] Dhawan S, Machado JAT, Brzeziński DW, Osman MS. A Chebyshev wavelet collocation method for some types of differential problems. Symmetry. 2021;13(4):536.10.3390/sym13040536Search in Google Scholar

[4] Bayones FS, Nisar KS, Khan KA, Raza N, Hussien NS, Osman MS, et al. Magneto-hydrodynamics (MHD) flow analysis with mixed convection moves through a stretching surface. AIP Adv. 2021;11(4):045001.10.1063/5.0047213Search in Google Scholar

[5] Aljahdali M, El-Sherif AA. Equilibrium studies of binary and mixed-ligand complexes of zinc (II) involving 2-(aminomethyl)-benzimidazole and some bio-relevant ligands. J Solution Chem. 2012;41(10):1759–76.10.1007/s10953-012-9908-2Search in Google Scholar

[6] Mohamed MM, El-Sherif AA. Complex formation equilibria between zinc (II), nitrilo-tris (methyl phosphonic acid) and some bio-relevant ligands. the kinetics and mechanism for zinc (II) ion promoted hydrolysis of glycine methyl ester. J Solution Chem. 2010;39(5):639–53.10.1007/s10953-010-9535-8Search in Google Scholar

[7] Chowdhury MA, Miah MM, Ali HS, Chu YM, Osman MS. An investigation to the nonlinear (2 + 1)-dimensional soliton equation for discovering explicit and periodic wave solutions. Results Phys. 2021;23:104013.10.1016/j.rinp.2021.104013Search in Google Scholar

[8] Malik S, Almusawa H, Kumar S, Wazwaz AM, Osman MS. A (2 + 1)-dimensional Kadomtsev–Petviashvili equation with competing dispersion effect: Painlevé analysis, dynamical behavior and invariant solutions. Results Phys. 2021;23:104043.10.1016/j.rinp.2021.104043Search in Google Scholar

[9] Yang M, Osman MS, Liu JG. Abundant lump-type solutions for the extended (3 + 1)-dimensional Jimbo–Miwa equation. Results Phys. 2021;23:104009.10.1016/j.rinp.2021.104009Search in Google Scholar

[10] El-Sherif AA, Shoukry MM. Equilibrium investigation of complex formation reactions involving copper (II), nitrilo-tris (methyl phosphonic acid) and amino acids, peptides or DNA constitutents. The kinetics, mechanism and correlation of rates with complex stability for metal ion promoted hydrolysis of glycine methyl ester. J Coord Chem. 2006;59(14):1541–56.10.1080/00958970600561399Search in Google Scholar

[11] El-Sherif AA, Shoukry MM. Coordination properties of tridentate (N, O, O) heterocyclic alcohol (PDC) with Cu (II): Mixed ligand complex formation reactions of Cu (II) with PDC and some bio-relevant ligands. Spectrochim Acta A Mol Biomol Spectrosc. 2007;66(3):691–700.10.1016/j.saa.2006.04.013Search in Google Scholar PubMed

[12] Khan K, Akbar MA. Exact and solitary wave solutions for the Tzitzeica-Dodd-Bullough and the modified KdV-Zakharov-Kuznetsov equations using the modified simple equation method. Ain Shams Eng J. 2013;4:903–9.10.1016/j.asej.2013.01.010Search in Google Scholar

[13] Zhang S, Li J, Zhang L. A direct algorithm of exp-function method for non-linear evolution equations in fluids. Therm Sci. 2016;20(3):881–4.10.2298/TSCI1603881ZSearch in Google Scholar

[14] Osman MS, Rezazadeh H, Eslami M, Neirameh A, Mirzazadeh M. Analytical study of solitons to Benjamin-Bona-Mahony-Peregrine equation with power law nonlinearity by using three methods. U Politeh Buch Ser A. 2018;80(4):267–78.Search in Google Scholar

[15] Mohyud-Din ST, Irshad A, Ahmed N, Khan U. Exact solutions of (3 + 1)-dimensional generalized KP equation arising in physics. Results Phys. 2017;127(3):67–85.10.1016/j.rinp.2017.10.007Search in Google Scholar

[16] Akbar MA, Akinyemi L, Yao SW, Jhangeer A, Rezazadeh H, Khater MMA, et al. Soliton solutions to the Boussinesq equation through sine-Gordon method and Kudryashov method. Results Phys. 2021;25:1–10.10.1016/j.rinp.2021.104228Search in Google Scholar

[17] Mirzazadeh M, Eslami M, Zerrad E, Mahmood MF, Biswas A, Belic M. Optical solitons in nonlinear directional couplers by sine-cosine function method and Bernoulli's equation approach. Nonlinear Dyn. 2015;81(4):1933–49.10.1007/s11071-015-2117-ySearch in Google Scholar

[18] Inc M, Kilic B, Baleanu D. Optical soliton solutions of the pulse propagation generalized equation in parabolic-law media with space-modulated coefficients. Optik. 2016;127:1056–8.10.1016/j.ijleo.2015.10.020Search in Google Scholar

[19] Zhang LD, Tian SF, Peng WQ, Zhang TT, Yan XJ. The dynamics of lump, lumpoff and rogue wave solutions of (2 + 1)-dimensional Hirota-Satsuma-Ito equations; East Asian. J Appl Math. 2020;10(2):243–55.10.4208/eajam.130219.290819Search in Google Scholar

[20] Khan K, Akbar MA. Travelling wave solutions of some coupled nonlinear evolution equations. ISRN Math Phys. 2013;2013:8.10.1155/2013/685736Search in Google Scholar

[21] Inc M, Kilic B. The first integral method for the perturbed Wadati-Segur-Ablowitz equation with time dependent coefficient. Kuwait J Sci. 2016;43:81–7.Search in Google Scholar

[22] Hosseini K, Aligoli M, Mirzazadeh M, Eslami M, Gómez-Aguilar JF. Dynamics of rational solutions in a new generalized Kadomtsev-Petviashvili equation. Modern Phys Lett B. 2019;33(35):1950437.10.1142/S0217984919504372Search in Google Scholar

[23] Yépez-Martíneza H, Gómez-Aguilarb JF, Baleanu D. Beta-derivative and sub-equation method applied to the optical solitons in medium with parabolic law nonlinearity and higher order dispersion. Optik. 2018;155:357–65.10.1016/j.ijleo.2017.10.104Search in Google Scholar

[24] Osman MS. Multi-soliton rational solutions for quantum Zakharov-Kuznetsov equation in quantum magnetoplasmas. Waves Random Complex. 2016;26(4):434–43.10.1080/17455030.2016.1166288Search in Google Scholar

[25] Ahmad H, Seadawy AR, Khan TA, Thounthong P. Analytic approximate solutions for some nonlinear parabolic dynamical wave equations. J Taibah Uni Sci. 2020;14(1):346–58.10.1080/16583655.2020.1741943Search in Google Scholar

[26] Kundu PR, Almusawa H, Fahim MRA, Islam ME, Akbar MA, Osman MS. Linear and nonlinear effects analysis on wave profiles in optics and quantum physics. Results Phys. 2021;23:103995.10.1016/j.rinp.2021.103995Search in Google Scholar

[27] Helal MA, Seadawy AR, Zekry MH. Stability analysis of solitary wave solutions for the fourth-order nonlinear Boussinesq water wave equation. Appl Math Comput. 2014;232:1094–103.10.1016/j.amc.2014.01.066Search in Google Scholar

[28] Kumar S, Kour B, Yao SW, Inc M, Osman MS. Invariance analysis, exact solution and conservation laws of (2 + 1) dim fractional Kadomtsev-Petviashvili (KP) system. Symmetry. 2021;13(3):477.10.3390/sym13030477Search in Google Scholar

[29] Khater MMA, Lu D, Attia RAM. Erratum: Dispersive long wave of nonlinear fractional Wu-zhang system via a modified auxiliary equation method [AIP Adv. 9, 025003 (2019)]. AIP Adv. 2019;9(4):049902.10.1063/1.5096005Search in Google Scholar

[30] Akbulut A, Almusawa H, Kaplan M, Osman MS. On the conservation laws and exact Solutions to the (3 + 1)-dimensional modified KdV-Zakharov-Kuznetsov equation. Symmetry. 2021;13(5):765.10.3390/sym13050765Search in Google Scholar

[31] Li J, Qui Y, Lu D, Attia RAM, Khater MMA. Study on the solitary wave solutions of the ionic currents on microtubules equation by using the modified Khater method. Thermal Sci. 2019;23(6):2053–62.10.2298/TSCI190722370LSearch in Google Scholar

[32] Islam ME, Akbar MA. Stable wave solutions to the Landau-Ginzburg-Higgs equation and the modified equal width wave equation using the IBSEF method. Arab J Basic Appl Sci. 2020;27(1):270–8.10.1080/25765299.2020.1791466Search in Google Scholar

[33] Islam ME, Kundu PR, Akbar MA, Gepreel KA, Alotaibi H. Study of the parametric effect of self-control waves of the Nizhnik-Novikov-Veselov equation by the analytical solutions. Results Phys. 2021;22:103887.10.1016/j.rinp.2021.103887Search in Google Scholar

[34] Akbar MA, Ali NHM, Tanjim T. Outset of multiple soliton solutions to the nonlinear Schrödinger equation and the coupled burgers equation. J Phys Commun. 2019;3:095013.10.1088/2399-6528/ab3615Search in Google Scholar

[35] Gao YT, Tian B. On the variant Boussinesq equations. Z Naturforsch. 1997;52a:335–6.10.1515/zna-1997-0406Search in Google Scholar

[36] Yan Z, Zhang H. New explicit and exact travelling wave solutions for a system of variant Boussinesq equations in mathematical physics. Phys Lett A. 1999;252:291–6.10.1016/S0375-9601(98)00956-6Search in Google Scholar

[37] Jabbari A, Kheiri H, Bekir A. Analytical solution of variant Boussinesq equations. Math Meth Appl Sci. 2014;37(37):931–6.10.1002/mma.2853Search in Google Scholar

[38] Ul-Hasssan QM. Soliton solutions of the variant Boussinesq equation through the exp-function method. Univ Wah J Sci Tech. 2017;1:24–30.Search in Google Scholar

[39] Manafian J, Jalali J, Alizadehdiz A. Some new analytical solutions of the variant Boussinesq equations. Opt Quant Electron. 2018;50:80.10.1007/s11082-018-1345-zSearch in Google Scholar

[40] Osman MS, Baleanu D, Adem AR, Hosseini K, Mirzazadeh M, Eslami M. Double-wave solutions and Lie symmetry analysis to the (2 + 1)-dimensional coupled burgers equations. Chin J Phys. 2020;63:122–9.10.1016/j.cjph.2019.11.005Search in Google Scholar

[41] Alharbi A, Almatrafi MB. Exact and numerical solitary wave structures to the variant Boussinesq system. Symmetry. 2020;12:1473.10.3390/sym12091473Search in Google Scholar

[42] Mohapatra SC, Guedes Soares C. Comparing solutions of the coupled Boussinesq equations in shallow water. In: Guedes Soares C, Santos TA, editors, Maritime technology and engineering. London, UK: Taylor & Francis Group; 2015. p. 947–54. ISBN: 978-1-138-02727-5.Search in Google Scholar

[43] Mohapatra SC, Fonseca RB, Guedes Soares C. Comparison of analytical and numerical simulations of long nonlinear internal waves in shallow water. J Coastal Res. 2018;34(4):928–38.10.2112/JCOASTRES-D-16-00193.1Search in Google Scholar

Received: 2021-03-08

Revised: 2021-05-21

Accepted: 2021-11-29

Published Online: 2022-08-14

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

https://doi.org/10.1515/phys-2022-0071

Keywords for this article

improved Bernoulli subequation function method; the new auxiliary equation method; the variant Boussinesq equation; exact solutions; soliton

Creative Commons

BY 4.0