Structure of analytical ion-acoustic solitary wave solutions for the dynamical system of nonlinear wave propagation

Hanadi Zahed; Aly R. Seadawy; Mujahid Iqbal

doi:10.1515/phys-2022-0030

Article Open Access

Structure of analytical ion-acoustic solitary wave solutions for the dynamical system of nonlinear wave propagation

Hanadi Zahed , Aly R. Seadawy and Mujahid Iqbal

Published/Copyright: May 6, 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

In the present study, the ion-acoustic solitary wave solutions for Kadomtsev–Petviashvili (KP) equation, potential KP equation, and Gardner KP equation are constructed. The nonlinear KP equations are studying the nonlinear process of waves without collisions plasma and having non-isothermal electrons and cold ions. Two-dimensional ion-acoustic solitary waves (IASWs) in magnetized plasma are consisting of electrons and ions. We obtained the ion-acoustic solitary wave solutions same as dark and bright, kink and anti-kink wave solitons. The physical phenomena of various structures for IASWs are represented graphically with symbolic computations. These results are more helpful in the development of soliton dynamics, quantum plasma, dynamic of adiabatic parameters, fluid dynamics, and industrial phenomena.

Keywords: modified mathematical method; Kadomtsev–Petviashvili equation; potential KP equation; Gardner KP equation; ion-acoustic solitary wave solutions; solitons

1 Introduction

Kadomtsev and Petviashvili in 1970, first time introduced the important nonlinear Kadomtsev–Petviashvili (KP) equation, which is the generalized Korteweg–de Vries (KdV) equation of two space variables [1]. The KP equation explains the phenomena for weakly dispersive waves. The KP equation can be used as a model of long wavelength water waves having weak nonlinear restoring forces and frequency dispersion. From the last five decades, a much analytical and numerical research has been carried out on various forms of KP equation. In recent years, KP equation has large deal of interest because of its explicit solutions including multi-solitons, periodic, and rational solutions in variables x and y [2]. The solutions and explicit solutions of KP equations are obtained by using the elementary techniques, see ref. [2,3]. The numerical solutions of KP equation by numerical approach are found in ref. [4]. By adomian decomposition method, the numerical solutions including solitons are derived in ref. [5].

The dusty plasma is ionized gas, which contains a little particle of matter in solid form, has size range from tens of nanometers to hundreds of microns. The solitons are solitary waves which represent an important feature for nonlinearity marvels in a system of space. The nonlinear partial differential equations (NLPDEs) have a specified type of contained results and possess different important features [6,7,8, 9,10,11, 12,13,14, 15,16].

The study of dust plasmas containing huge dust particles is very significant to understand various behaviors for industrial physical applications, space, and astrophysical [17]. The dust plasma has potential applications in the research for space medium, e.g., zone of asteroid, cometary tails, planet rings, earth environments, medium of interstellar, and astrophysical [18,19,20, 21,22]. The grains in the form of dust have negative charge because the number of various processes of charging, for example, ultra violet radiation, field emission, and current of plasma [23, 24,25]. The effects of nonthermal supply for dust fluid and ions temperature on random amplitude applied in solitary construction of electrostatic and effective models, and exist in hot nonthermal dusty plasmas models containing hot nonthermal distributed ions and dust fluid [26,27].

The Sagdeev pseudo-potential and solitons of rarefactive are finding Bolltzmannian electrons and ions for dust-acoustic solitary waves under dust plasma [28]. Sagdeev reported that ion-acoustic waves under unmagnetized plasma have isothermal electrons, hot and cold ions [29]. He studied oscillation of scientific particles for nonlinear wave in good potential such as equation of integral energy to find the families of basic equations for plasma dynamic. He derived the Sagdeev potential, which is useful to measure the presence of various characteristics of soliton in plasma. In a dynamical system, the number of trapped ions is ignored because of a very small number of ions. The number of trap electrons is more, but the electron distribution is less under an ion-acoustic wave. So the distribution effect of trap electron is not Maxwellian, and these are not counted under the observation.

In the recent years, many researchers constructed the solutions for NLPDEs. Therefore, to construct the solutions for these NLPDEs many scholars and mathematicians developed many techniques. The list of some techniques are the scheme of Backlund transformed, the Hirota bilinear technique, the scheme of Darboux transformed, Exp-function method, the scheme of Jacobian technique, the scheme of trial equation, simple equation technique, the extension of fan-sub equation technique, the extend mapping technique, the technique of sinh-cosh, the direct algebraic technique, the scheme F-expansion technique, the modified tanh technique, the Riccati equation mapping technique, the extension of auxiliary equation technique, and Seadawy techniques [30,31,32, 33,34,35, 36,37,38, 39,40,41, 42,43,44, 45,46,47, 48,49,50, 51,52,53]. We constructed the bright-dark solitons as in form ion-acoustic solitary wave solutions for three nonlinear PDEs by modified technique [54,55,56, 57,58,59, 60,61,62, 63,64].

This article is arranged as follows: Introduction is given in Section 1. In Section 2, the proposed technique is described. Formulation for Sagdeev potential by basic equations is explained in Section 3. In Section 4, we construct solitonic solutions for equations by the described technique. In Section 5, the obtained results are discussed. Conclusion of this study is given in Section 6.

2 Algorithm of proposed method

Consider the NLPDEs in general form as:

(1) S ( ϕ t , ϕ x , ϕ y , ϕ t t , ϕ x x , ϕ y y , ϕ x t , … ) = 0 ,

where S denotes polynomial function for ϕ ( x , y , t ) and its all derivatives. The main features for proposed technique are explained as follows:

Step 1: Consider the traveling wave transform as:

(2) V ( ζ ) = ϕ ( x , y , t ) , ζ = κ x + λ y + μ t .

By using Eq. (2), the ODE (ordinary differential equation) for Eq. (1), we obtain

(3) R ( V ′ , V ″ , V ‴ , … ) = 0 ,

where R denotes polynomial function in V ( ζ ) and their derivatives.

Step 2: The trial solution for Eq. (3) is considered as follows:

(4) V ( ζ ) = ∑ l = 0 m a l Ψ ( ζ ) l + ∑ l = − 1 − m b − l Ψ ( ζ ) l + ∑ l = 2 m c l Ψ ( ζ ) l − 2 Ψ ′ ( ζ ) + ∑ l = 1 m d l Ψ ′ ( ζ ) Ψ ( ζ ) l ,

where a l , b l , c l , d l , ( l = 0 , 1 , 2 , 3 , … ) denote constant parameters to be determined later, while Ψ ( ζ ) satisfy the following auxiliary equation:

(5) d Ψ d ζ = β 1 Ψ 2 ( ζ ) + β 2 Ψ 3 ( ζ ) + β 3 Ψ 4 ( ζ ) .

Step 3: Apply the homogeneous rule in Eq. (3), to balance the highest order derivative and nonlinear term to determine m for Eq. (4).

Step 4: Substitute Eq. (4) in Eq. (5) and Eq. (3), and combine each coefficient of Ψ ′ k ( ζ ) Ψ l ( ζ ) ( k = 0 , 1 ; l = 1 , 2 , 3 , … m ) then each coefficient equates to zero. We obtained a system of algebraic equations and solved them by symbolic computation and obtained parameter values. Substituting parameter values and Ψ ( ζ ) in Eq. (4), solutions for Eq. (1) are obtained.

3 Formulation of Sagdeev potential by set of basic equations

Consider the collisionless plasma containing nonisothermal cold ions and electrons. The plasma model considered the free electrons temperature ( T e f ) , with trap electrons temperature ( T e t ) , move with good potential, continuously lose energy, and a consequences upper electronsmove front and back with good potential and become trap in good potential under plasma. The electron density and effect of trapped electrons [65,66] are defined as:

(6) n e ( ϕ ) = ∫ f e ( x , r ) d x = exp ( ϕ ) e r f e ( ϕ ) + 1 σ { exp ( σ ϕ ) e r f ( σ ϕ ) , where σ > 0 ; 2 π exp − ( − σ ϕ ) 2 ∫ 0 − σ ϕ exp ( X 2 d X ) , where σ < 0 .

The function of distribution of electron is f e ( x , r ) and ratio between temperature of free and trapped electrons is σ = T e f T e t . By applying Taylor series in Eq. (6), for ϕ < 1 , electron density n e is obtained as:

(7) n e = 1 + ϕ − 4 3 β ϕ 3 2 + 1 2 ϕ 2 − 8 15 γ ϕ 5 2 + 1 6 ϕ 3 + … ,

where n e denotes electron density and ϕ is the potential normalized by unperturbed electron density n 0 and κ B T e t e , respectively. κ B is the Boltzmann constant, β = 1 − σ π , and γ = 1 − σ 2 π , σ = T e f T e t . The plasma has flat top and Maxwellian distribution for σ = 0 and σ = 1 . If β = 0 and γ = 0 , determine electron densities isothermal plasmas. The conditions are 0 < β < 1 π , 0 < γ < 1 π , for the non-isothermal plasma. Consider the families of basic equations represent dynamic plasma under describe fluid, followed by Das and Sen (1994) as:

(8) ∂ n ∂ t + ∂ n ∂ x u = 0 ,

(9) − ∂ ϕ ∂ x = u ∂ n ∂ x + ∂ u ∂ t .

Poisson’s equation is accompanied as:

(10) ∂ 2 ϕ ∂ x 2 = n e − n ,

where n is the ion density, normal by the n 0 move by velocity, u normal with ion acoustic speed c s = κ B T e f m i 1 2 . The space x has Debye length λ D = κ B T e f 4 π n 0 e 2 1 2 and time t − 1 has ion-frequency 4 π n 0 e 2 κ B T e f 1 2 . For checking the stable solutions for Eqs. (8)–(10), suppose the physical parameters depend on space and time x , t with relation ( ξ = x − μ t ) , and the wave Mach number is μ . Consider ξ ⟶ ± ∞ for bounded solution, the plasma parameters u and ϕ are set to be zero. By applying these conditions the families of basic Eqs. (8)–(10) are changed as:

(11) 1 2 d ϕ d ξ 2 + V ( ϕ , μ ) = 0 .

Sagdeev potential [24] is taken as:

(12) V ( ϕ , μ ) = 1 − exp ( ϕ ) + 8 15 β ϕ 5 2 + 16 105 γ ϕ 7 2 − μ 2 1 − 2 ϕ μ 2 1 2 − 1 .

With parameter selection μ , β , and γ , various types of solutions for Eq. (12) are obtained.

4 Applications of described method for nonlinear dynamical equations

Here, we applied the proposed technique for the construction of solitary wave solutions for three nonlinear dynamical equations.

4.1 KP equation

Consider the KP equation as:

(13) ( ϕ t + 6 ϕ ϕ x + ϕ x x x ) x + 3 ϕ y y = 0 ,

and the transformation of wave as:

(14) V ( ζ ) = ϕ ( x , y , t ) , ζ = κ x + λ y + μ t ,

putting Eq. (14) in Eq. (13), we obtain:

(15) μ κ V ″ + 6 κ 2 ( V ′ ) 2 + 6 κ 2 V V ″ + κ 4 V ⁗ + 3 λ 2 V ″ = 0 ,

and integrate Eq. (15), once with respect to ζ with zero constant of integration,which becomes

(16) μ κ V ′ + 6 κ 2 V V ′ + κ 4 V ‴ + 3 λ 2 V ′ = 0 ,

and balancing higher order derivative and nonlinear term in Eq. (16), we obtain m = 2 . The general solution for Eq. (16) is obtained as:

(17) V ( ζ ) = a 0 + a 1 Ψ ( ζ ) + a 2 Ψ ( ζ ) 2 + b 1 Ψ ( ζ ) + b 2 Ψ ( ζ ) 2 + c 2 Ψ ′ ( ζ ) + d 1 Ψ ′ ( ζ ) Ψ ( ζ ) + d 2 Ψ ′ ( ζ ) Ψ ( ζ ) 2 .

Substitute Eq. (17) in Eq. (16), and combine each coefficient of Ψ ′ k ( ζ ) Ψ l ( ζ ) ( k = 0 , 1 ; i = 1 , 2 , 3 , 4 , 5 , … m ), and equate each coefficient to zero. We obtain system of equations. After solving the system of equations by Mathematica, the values of constants are obtained as:

Case-I

(18) a 0 = a 0 , a 1 = − β 2 ( d 2 + κ 2 ) , a 2 = − 2 β 3 κ 2 − β 3 d 2 , b 1 = b 2 = c 2 = d 1 = 0 , β 1 = β 2 2 4 β 3 , d 2 = d 2 , μ = − 24 a 0 β 3 κ 2 − β 2 2 κ 4 − 12 β 3 λ 2 − 6 β 2 2 d 2 κ 2 4 β 3 κ .

Substituting Eq. (18) in Eq. (17), then solutions for Eq. (13) are obtained as:

(19) ϕ 1 ( x , y , t ) = a 0 + β 1 ( d 2 + κ 2 ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 × β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − β 1 2 β 3 ( d 2 + 2 κ 2 ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 .

(20) ϕ 2 ( x , y , t ) = a 0 + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − 1 2 β 2 β 1 β 3 ( − d 2 − κ 2 ) ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 − β 1 ( d 2 + 2 κ 2 ) ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 4 ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 .

(21) ϕ 3 ( x , y , t ) = a 0 + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 3 ( d 2 + 2 κ 2 ) ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 + β 2 ( − d 2 − κ 2 ) − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 .

(Figure 1).

$Figure 1 Three dimensional, two dimensional, and contour plots for Eq. (19) represent dark soliton while ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 1

Three dimensional, two dimensional, and contour plots for Eq. (19) represent dark soliton while ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-II

(22) a 0 = − β 1 κ 4 − 6 β 1 d 2 κ 2 − κ μ − 3 λ 2 6 κ 2 , a 1 = ± β 1 β 3 ( 2 κ 2 + d 2 ) , a 2 = − 2 β 3 κ 2 − β 3 d 2 , β 2 = ∓ 2 β 1 β 3 , d 2 = d 2 , b 1 = b 2 = c 2 = d 1 = 0 .

Putting Eq. (22), in Eq. (17), only positive value for a 1 , solutions for Eq. (13) are obtained as:

(23) ϕ 4 ( x , y , t ) = β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − β 1 ( 6 d 2 κ 2 + κ 4 ) + κ μ + 3 λ 2 6 κ 2 − 2 β 1 3 / 2 β 3 ( d 2 + κ 2 ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 β 2 − β 1 2 β 3 ( d 2 + 2 κ 2 ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 .

(24) ϕ 5 ( x , y , t ) = β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 ( d 2 + κ 2 ) ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) η + cosh [ β 1 ( ζ + ζ 0 ) ] − β 1 ( d 2 + 2 κ 2 ) ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 4 ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 − β 1 ( 6 d 2 κ 2 + κ 4 ) + κ μ + 3 λ 2 6 κ 2 .

(25) ϕ 6 ( x , y , t ) = β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 ( 6 d 2 κ 2 + κ 4 ) + κ μ + 3 λ 2 6 κ 2 + 2 β 1 β 3 ( d 2 + κ 2 ) − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 − β 3 ( d 2 + 2 κ 2 ) ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

(Figure 2)

$Figure 2 Three dimensional, two dimensional, and contour plots for Eq. (20) representing bright soliton while ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 2

Three dimensional, two dimensional, and contour plots for Eq. (20) representing bright soliton while ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-III

(26) a 0 = − β 1 κ 4 − 6 β 1 d 2 κ 2 − κ μ − 3 λ 2 6 κ 2 , a 1 = ± 2 d 2 β 1 β 3 , a 2 = − β 3 d 2 , β 2 = ∓ 2 β 1 β 3 , b 1 = b 2 = c 2 = d 1 = 0 , d 2 = d 2 .

Substituting Eq. (26), in Eq. (17), only positive value for a 1 , solutions for Eq. (13) are obtained as:

(27) ϕ 7 ( x , y , t ) = − κ μ + 3 λ 2 6 κ 2 + β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − 1 − κ 2 6 − 2 β 1 3 / 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 β 2 − β 1 2 β 3 d 2 ε coth 1 2 β 1 ( κ x + λ y + μ t + ζ 0 ) + 1 2 β 2 2 ,

(28) ϕ 8 ( x , y , t ) = β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 − 1 4 β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 − β 1 ( 6 d 2 κ 2 + κ 4 ) + κ μ + 3 λ 2 6 κ 2 .

(Figure 3).

$Figure 3 Three dimensional, two dimensional, and contour plots for Eq. (21) representing anti-kink wave soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 2.5 \mu =2.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , p = 0.3 p=0.3 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 3

Three dimensional, two dimensional, and contour plots for Eq. (21) representing anti-kink wave soliton when ζ 0 = 0.3 , κ = 2 , μ = 2.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , p = 0.3 , a 0 = 1 , d 2 = 1 , y = 1 .

(29) ϕ 9 ( x , y , t ) = β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 ( 6 d 2 κ 2 + κ 4 ) + κ μ + 3 λ 2 6 κ 2 + 2 β 1 β 3 d 2 − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ) ] η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 ,

where ζ = κ x + λ y + μ t (Figure 4).

$Figure 4 Three dimensional, two dimensional, and contour plots for Eq. (23) representing dark soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 4

Three dimensional, two dimensional, and contour plots for Eq. (23) representing dark soliton when ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , a 0 = 1 , d 2 = 1 , y = 1 .

4.2 Potential KP equation

Consider the potential KP equation as:

(30) ϕ x t + 3 2 ϕ x ϕ x x + 1 4 ϕ x x x x + 3 4 ϕ y y = 0 .

We apply the wave transformation:

(31) V ( ζ ) = ϕ ( x , y , t ) , ζ = κ x + λ y + μ t .

Putting Eq. (31) in Eq. (30), we obtain:

(32) μ κ V ″ + 3 2 κ 3 V ′ V ″ + 1 4 κ 4 V ⁗ + 3 4 λ 2 V ″ = 0 .

Integrating Eq. (32), once with respect to ζ and integration constant equate to zero, we obtain:

(33) 3 4 λ 2 + μ κ V ′ + 3 4 κ 3 ( V ′ ) 2 + 1 4 κ 4 V ‴ = 0 .

Balancing highest order derivative and nonlinear term in Eq. (33), we obtain m = 2 . The general solution for Eq. (33) is given as:

(34) V ( ζ ) = a 0 + a 1 Ψ ( ζ ) + a 2 Ψ ( ζ ) 2 + b 1 Ψ ( ζ ) + b 2 Ψ ( ζ ) 2 + c 2 Ψ ′ ( ζ ) + d 1 Ψ ′ ( ζ ) Ψ ( ζ ) + d 2 Ψ ′ ( ζ ) Ψ ( ζ ) 2 .

Putting Eq. (34) in Eq. (33), and after solving, we obtain:

Case-I

(35) a 1 = ± 2 d 2 β 1 β 3 , a 2 = − β 3 d 2 , d 1 = − 2 κ , d 2 = d 2 , μ = − β 1 κ 4 − 3 λ 2 4 κ , β 2 = ∓ 2 β 1 β 3 , b 1 = b 2 = c 2 = 0 .

Substitute Eq. (35) in Eq. (34), only positive value for a 1 , then solutions for Eq. (30) are obtained as:

(36) ϕ 1 ( x , y , t ) = a 0 + ( β 1 ε csch 1 2 β 1 ( ζ + ζ 0 ) 2 β 2 2 ( β 1 d 2 ε + 2 κ ( cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] − 1 ) ) − 8 β 1 β 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 3 − 4 β 1 3 / 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 4 4 β 2 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 ,

(37) ϕ 2 ( x , y , t ) = a 0 − 2 β 1 κ ε ( η cosh [ β 1 ( ζ + ζ 0 ) ] + 1 ) ∕ ( η + cosh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 − 1 4 β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 ,

(38) ϕ 3 ( x , y , t ) = a 0 − 2 β 1 κ ε ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) / ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) + 2 β 1 β 3 d 2 − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

(Figure 5).

$Figure 5 Three dimensional, two dimensional, and contour plots for Eq. (36) representing solitary wave when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 5

Three dimensional, two dimensional, and contour plots for Eq. (36) representing solitary wave when ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-II

(39) a 1 = ± β 3 ( κ + 2 β 1 d 2 ) , a 2 = − β 3 d 2 , d 1 = − κ , μ = − β 1 κ 4 − 3 λ 2 4 κ , d 2 = d 2 , β 2 = − 2 β 1 β 3 , b 1 = b 2 = c 2 = 0 .

Substituting Eq. (39) in Eq. (34), only positive value for a 1 , then solutions for Eq. (30) are obtained as:

(40) ϕ 4 ( x , y , t ) = a 0 + β 1 κ ε csch 1 2 β 1 ( ζ + ζ 0 ) 2 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 2 + β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − β 1 β 3 ( 2 β 1 d 2 + κ ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 β 2 − β 1 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 ,

(41) ϕ 5 ( x , y , t ) = a 0 − β 1 κ ε ( η cosh [ β 1 ( ζ + ζ 0 ) ] + 1 ) ∕ ( η + cosh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) − 1 2 β 1 ( 2 β 1 d 2 + κ ) ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − 1 4 β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 ,

(42) ϕ 6 ( x , y , t ) = a 0 + β 1 κ ε ( η ( − p 2 + 1 ) cosh [ β 1 ( ζ + ζ 0 ) ] + p sinh [ β 1 ( ζ + ζ 0 ) ] − 1 ) / ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] × η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) + β 3 ( 2 β 1 d 2 + κ ) − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

(Figure 6).

$Figure 6 Three dimensional, two dimensional, and contour plots for Eq. (37) representing kink wave soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 1.5 \mu =1.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 6

Three dimensional, two dimensional, and contour plots for Eq. (37) representing kink wave soliton when ζ 0 = 0.3 , κ = 2 , μ = 1.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-III

(43) a 1 = − β 2 d 2 , a 2 = − β 2 2 d 2 4 β 1 , d 1 = − 2 κ , μ = − β 1 κ 4 − 3 λ 2 4 κ , d 2 = d 2 , β 3 = β 2 2 4 β 1 , b 1 = b 2 = c 2 = 0 .

Substituting Eq. (43) in Eq. (34), then solutions for Eq. (30), we obtain:

(44) ϕ 7 ( x , y , t ) = a 0 + β 1 κ ε csch 1 2 β 1 ( ζ + ζ 0 ) 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 + 1 4 3 + 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) − ε 2 coth 2 1 2 β 1 ( ζ + ζ 0 ) + ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 d 2 β 1 .

(45) ϕ 8 ( x , y , t ) = a 0 − 2 β 1 κ ε ( η cosh [ β 1 ( ζ + ζ 0 ) ] + 1 ) ∕ ( η + cosh [ β 1 ( ζ + ζ 0 ) ] × η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) + 1 2 β 2 β 1 β 3 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 − β 2 2 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 16 β 3 .

(46) ϕ 9 ( x , y , t ) = a 0 − 2 β 1 κ ε ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) / × ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 2 d 2 − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 − β 2 2 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 4 β 1 .

(Figure 7).

$Figure 7 Three dimensional, two dimensional, and contour plots for Eq. (38) representing bright soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 1.5 \mu =1.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , p = 0.8 p=0.8 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 7

Three dimensional, two dimensional, and contour plots for Eq. (38) representing bright soliton when ζ 0 = 0.3 , κ = 2 , μ = 1.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , p = 0.8 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-IV

(47) a 1 = β 2 κ − 2 β 1 β 2 d 2 2 β 1 , a 2 = − β 2 2 d 2 4 β 1 , b 1 = b 2 = c 2 = 0 , d 1 = − κ , d 2 = d 2 , μ = − β 1 κ 4 − 3 λ 2 4 κ , β 3 = β 2 2 4 β 1 .

Substituting Eq. (47) in Eq. (34), then solutions for Eq. (30), we obtain:

(48) ϕ 10 ( x , y , t ) = 1 4 4 a 0 + 2 β 1 κ − ε coth 1 2 β 1 ( ζ + ζ 0 ) + ε csch 1 2 β 1 ( ζ + ζ 0 ) 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 − 1 + 3 + 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) − ε 2 coth 1 2 β 1 ( ζ + ζ 0 ) 2 + ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 d 2 β 1 .

(49) ϕ 11 ( x , y , t ) = a 0 − β 1 κ ε ( η cosh [ β 1 ( ζ + ζ 0 ) ] + 1 ) ∕ ( η + cosh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 2 β 1 β 3 ( κ − 2 β 1 d 2 ) ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 4 β 1 − β 2 2 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 16 β 3 .

ϕ 12 ( x , y , t ) = a 0 + β 1 κ ε ( η ( − p 2 + 1 ) cosh [ β 1 ( ζ + ζ 0 ) ] + p sinh [ β 1 ( ζ + ζ 0 ) ] − 1 ) / ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] × η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) + β 2 ( κ − 2 β 1 d 2 ) − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 2 β 1 − β 2 2 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 4 β 1 , (50)

where ζ = κ x + λ y + μ t (Figure 8).

$Figure 8 Three dimensional, two dimensional, and contour plots for Eq. (57) representing solitary wave when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = − 1 \eta =-1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , α = 1 \alpha =1 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 8

Three dimensional, two dimensional, and contour plots for Eq. (57) representing solitary wave when ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = 1 , η = − 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , α = 1 , a 0 = 1 , d 2 = 1 , y = 1 .

4.3 Gardner KP equation

Here we take the Gardner KP equation as:

(51) ( ϕ t + 6 ϕ ϕ x + 6 α ϕ 2 ϕ x + ϕ x x x ) x + ϕ y y = 0 .

We apply the wave transformation as:

(52) V ( ζ ) = ϕ ( x , y , t ) , ζ = κ x + λ y + μ t .

Substituting Eq. (52) in Eq. (51), ODE of Eq. (51), we obtain:

(53) ( μ κ V ′ + 6 κ 2 V V ′ + 6 α κ 2 V 2 V ′ + κ 4 V ‴ ) ′ + λ 2 V ″ = 0 .

Integrating Eq. (53), once with respect to ζ and integration constant equal to zero, Eq. (53) becomes

(54) μ κ V ′ + 6 κ 2 V V ′ + 6 α κ 2 V 2 V ′ + κ 4 V ‴ + λ 2 V ′ = 0 .

Balancing the highest order derivative and nonlinear term in Eq. (54), we obtain m = 2 . The general solution for Eq. (54) is as follows:

(55) V ( ζ ) = a 0 + a 1 Ψ ( ζ ) + a 2 Ψ ( ζ ) 2 + b 1 Ψ ( ζ ) + b 2 Ψ ( ζ ) 2 + c 2 Ψ ′ ( ζ ) + d 1 Ψ ′ ( ζ ) Ψ ( ζ ) + d 2 Ψ ′ ( ζ ) Ψ ( ζ ) 2 .

Substituting Eq. (55) in Eq. (54), and after solving, we obtain:

Case-I

(56) a 0 = − α β 2 κ − 2 β 3 − 4 α β 3 β 1 d 2 4 α β 3 , a 1 = − β 3 κ − α β 2 d 2 α , a 2 = β 3 ( − d 2 ) , d 2 = d 2 , μ = 3 α β 2 2 κ 4 − 8 α β 1 β 3 κ 4 − 8 α β 3 λ 2 + 12 β 3 κ 2 8 α β 3 κ , b 1 = b 2 = c 2 = d 1 = 0 ,

where α < 0 , and substitute Eq. (56) in Eq. (55), then we obtain solutions for Eq. (51) as:

(57) ϕ 1 ( x , y , t ) = 1 4 β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − α β 2 κ β 3 + 4 α β 1 d 2 + 2 α + 4 β 1 ( β 3 κ + α β 2 d 2 ) ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 α β 2 − 4 β 1 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 .

(Figure 9).

$Figure 9 Three dimensional, two dimensional, and contour plots for Eq. (58) representing kink wave soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 1.5 \mu =1.5 , λ = 1.4 \lambda =1.4 , ε = 1 \varepsilon =1 , η = 1 \eta =1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , α = 1 \alpha =1 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 9

Three dimensional, two dimensional, and contour plots for Eq. (58) representing kink wave soliton when ζ 0 = 0.3 , κ = 2 , μ = 1.5 , λ = 1.4 , ε = 1 , η = 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , α = 1 , a 0 = 1 , d 2 = 1 , y = 1 .

(58) ϕ 2 ( x , y , t ) = 1 4 4 β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 − α β 2 κ β 3 + 4 α β 1 d 2 + 2 α − 2 β 1 β 3 ( − β 3 κ − α β 2 d 2 ) ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 α .

(59) ϕ 3 ( x , y , t ) = β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − α β 2 κ β 3 + 4 α β 1 d 2 + 2 4 α + ( − β 3 κ − α β 2 d 2 ) − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 α − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

Case-II

(60) a 0 = − 2 α β 1 d 2 − 1 2 α , a 1 = − β 3 κ α , a 2 = − β 3 d 2 , d 2 = d 2 , μ = − 2 α β 1 κ 4 − 2 α λ 2 + 3 κ 2 2 α κ , b 1 = b 2 = c 2 = d 1 = 0 ,

where α < 0 , substituting Eq. (60) in Eq. (55), then solutions for Eq. (51) become:

(61) ϕ 4 ( x , y , t ) = 1 4 4 α β 1 β 3 κ ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 − 2 β 2 α β 2 + × d 2 β 1 − 4 + ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − 4 β 1 β 3 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 .

(62) ϕ 5 ( x , y , t ) = 1 4 − 2 α − 4 β 1 d 2 + 4 β 1 d 2 4 β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 / ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 + 2 β 1 β 3 β 3 κ ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 α .

(63) ϕ 6 ( x , y , t ) = − 1 2 α − β 1 d 2 − β 3 κ − ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] − 1 α + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

(Figure 10).

$Figure 10 Three dimensional, two dimensional, and contour plots for Eq. (59) representing dark soliton when ζ 0 = 0.3 {\zeta }_{0}=0.3 , κ = 2 \kappa =2 , μ = 0.5 \mu =0.5 , λ = 1.4 \lambda =1.4 , ε = − 1 \varepsilon =-1 , η = − 1 \eta =-1 , β 1 = 0.2 {\beta }_{1}=0.2 , β 2 = 0.4 {\beta }_{2}=0.4 , β 3 = 0.2 {\beta }_{3}=0.2 , p = 0.8 p=0.8 , α = 1 \alpha =1 , a 0 = 1 {a}_{0}=1 , d 2 = 1 {d}_{2}=1 , y = 1 y=1 .$

Figure 10

Three dimensional, two dimensional, and contour plots for Eq. (59) representing dark soliton when ζ 0 = 0.3 , κ = 2 , μ = 0.5 , λ = 1.4 , ε = − 1 , η = − 1 , β 1 = 0.2 , β 2 = 0.4 , β 3 = 0.2 , p = 0.8 , α = 1 , a 0 = 1 , d 2 = 1 , y = 1 .

Case-III

(64) a 0 = − 2 α β 1 d 2 − 1 2 α , a 2 = − β 3 d 2 , d 1 = ± κ α , d 2 = d 2 , μ = 4 α β 1 κ 4 − 2 α λ 2 + 3 κ 2 2 α κ , a 1 = b 1 = b 2 = c 2 = 0 ,

where α < 0 , substituting Eq. (64) only positive value for d 1 in Eq. (55), then solutions for Eq. (51) become:

(65) ϕ 7 ( x , y , t ) = 1 4 − 2 α − 2 β 1 κ ε csch 1 2 β 1 ( ζ + ζ 0 ) 2 α ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 + β 1 d 2 ε 2 csch 1 2 β 1 ( ζ + ζ 0 ) 4 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 − 4 − 4 β 1 2 β 3 d 2 ε coth 1 2 β 1 ( ζ + ζ 0 ) + 1 2 β 2 2 .

(66) ϕ 8 ( x , y , t ) = − 1 2 α + β 1 κ ε ( η cosh [ β 1 ( ζ + ζ 0 ) ] + 1 ) ∕ ( α ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) × η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) − β 1 d 2 + β 1 d 2 ( η ε cosh [ β 1 ( ζ + ζ 0 ) ] + ε ) 2 ∕ ( ( η + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − 1 4 β 1 d 2 ε sinh [ β 1 ( ζ + ζ 0 ) ] η + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 .

(67) ϕ 9 ( x , y , t ) = − 1 2 α + β 1 κ ε ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) / ( α ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) − β 1 d 2 + β 1 d 2 ε 2 ( η p 2 + 1 cosh [ β 1 ( ζ + ζ 0 ) ] − p sinh [ β 1 ( ζ + ζ 0 ) ] + 1 ) 2 / ( ( η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] ) 2 × ( η p 2 + 1 + p ε + cosh [ β 1 ( ζ + ζ 0 ) ] + ε sinh [ β 1 ( ζ + ζ 0 ) ] ) 2 ) − β 3 d 2 ε ( p + sinh [ β 1 ( ζ + ζ 0 ) ] ) η p 2 + 1 + cosh [ β 1 ( ζ + ζ 0 ) ] + 1 2 ,

where ζ = κ x + λ y + μ t .

5 Results and discussion

We successfully constructed the new solitary wave solutions including bright and dark, kink and anti-kink wave solitons for three nonlinear PDEs by the proposed new modified technique. Here we discuss the similarity and difference between our new determined results with previous literature which have been found by other methods. The scholar found a numerical solutions including solitons by using the adomian decomposition method which is mentioned in ref. [5]. Our constructed solutions are new, more general, and different from that obtained in the existing literature. These new solutions are obtained in rational, elliptic, trigonometric, and hyperbolic form including singular bright and dark solitons, kink and anti-kink wave solitions.

From the above complete discussion we conclude that our constructed solutions are different and new, which have not been determined in the existing literature. The all calculation shows that the proposed technique is more powerful, effective, and reliable for analytically investigating other various kinds of nonlinear PDEs.

6 Conclusion

In the present work, we successfully proposed modified technique on three NLPDEs. We have determined new ion acoustic solitary wave solutions for KP equation, potential KP equation, and Gardner KP equation by the modified mathematical method. We obtained the ion-acoustic solitary wave solutions including singular dark solitons, singular bright solitons, kink and anti-kink wave solitons. We also represent physical interpretation of some obtained results by two-dimensional, three-dimensional, and contour graphs by Mathematica to understand the physical structure for ion-acoustic solitary wave. The obtained solutions may have potential applications in quantum plasma, optical fibers, soliton dynamics, laser optics, dynamic of adiabatic parameters, fluid dynamics, and various other fields. The complete calculations show that the proposed technique is more efficient, powerful, straightforward, and fruitful to investigate analytically other various kinds of NLPDEs.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work with the project number (141/442). Also, the authors would like to extend their appreciation to Taibah University for its supervision support.

Funding information: Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work with the project number (141/442).
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] Kadomtsev BB, Petviashvili VI. On the stability of solitary waves in weakly dispersing media. Sov Phys Dokl. 1970;15:539–41. Search in Google Scholar

[2] Grünbaum FA. The Kadomtsev-Petviashvili equation: an alternative approach to the rank two solutions of Krichever and Notikov. Phys Lett A. 1989;139:146–50. 10.1016/0375-9601(89)90349-6Search in Google Scholar

[3] Latham GA. Solutions of KP equation associated to rank-three commuting differential operators over a singular curve. Phys D. 1990;41:55–66. 10.1016/0167-2789(90)90027-MSearch in Google Scholar

[4] Bratsos AG, Twizell EH. An explicit fine difference scheme for the solutions of Kadomtsev-Petviashvili. Int J Comput Math. 1998;68:175–87. 10.1080/00207169808804685Search in Google Scholar

[5] Wazwaz AM. A computational approach to soliton solutions of the Kadomtsev-Petviashili equation. Appl Math Comput. 2001;123(2):205–17. 10.1016/S0096-3003(00)00065-5Search in Google Scholar

[6] Misra AP, Wang Y. Dust-acoustic solitary waves in a magnetized dusty plasma with nonthermal electrons and trapped ions. Commun Nonlinear Sci Numer Simul. 2015;22:1360–9. 10.1016/j.cnsns.2014.07.017Search in Google Scholar

[7] Rani M, Ahmed N, Dragomir SS, Mohyud-Din ST, Khan I, Nisar KS. Some newly explored exact solitary wave solutions to nonlinear inhomogeneous Murnaghan’s rod equation of fractional order. J TAIBAH Univ Sci. 2021;15(1):97–110. 10.1080/16583655.2020.1841472Search in Google Scholar

[8] Ahmed N, Rafiq M, Baleanu D, Rehman MA, Khan I, Ali M, et al. Structure preserving algorithms for mathematical model of auto-catalytic glycolysis chemical reaction and numerical simulations. Europ Phys J Plus. 2020;135:522. 10.1140/epjp/s13360-020-00539-wSearch in Google Scholar

[9] Ghaffar A, Ali A, Ahmed S, Akram S, Junjua M-D, Baleanu D, et al. A novel analytical technique to obtain the solitary solutions for nonlinear evolution equation of fractional order. Adv Differ Equ. 2020;2020:308. 10.1186/s13662-020-02751-5Search in Google Scholar

[10] Mustafa G, TehminaEjaz S, Baleanu D, Ghaffar A, SooppyNisar K. A subdivision-based approach for singularly perturbed boundary value problem. Adv Differ Equ. 2020;2020:282. 10.1186/s13662-020-02732-8Search in Google Scholar

[11] Subashini R, Jothimani K, Nisar KS, Ravichandrand C. New results on nonlocal functional integro-differential equations via Hilfer fractional derivative. Alexandr Eng J. 2020;59(5):2891–9. 10.1016/j.aej.2020.01.055Search in Google Scholar

[12] Jajarmi A, Baleanu D, Zarghami Vahid K, Mohammadi Pirouz H, Asad JH. A new and general fractional Lagrangian approach: a capacitor microphone case study. Results Phys. 2021;31:104950. 10.1016/j.rinp.2021.104950Search in Google Scholar

[13] Erturk VS, Godwe E, Baleanu D, Kumar P, Asad J, Jajarmi A. Novel fractional-order Lagrangian to describe motion of beam on nanowire. Acta Physica Polonica A. 2021;140(3):265–72. 10.12693/APhysPolA.140.265Search in Google Scholar

[14] Baleanu D, Hassan Abadi M, Jajarmi A, Zarghami Vahid K, Nieto JJ. A new comparative study on the general fractional model of COVID-19 with isolation and quarantine effects. Alexandr Eng J. 2022;61(6):4779–91. 10.1016/j.aej.2021.10.030Search in Google Scholar

[15] Jajarmi A, Baleanu D, Zarghami Vahid K, Mobayen S. A general fractional formulation and tracking control for immunogenic tumor dynamics. Math Meth Appl Sci. 2022;45(2):667–80. 10.1002/mma.7804Search in Google Scholar

[16] Baleanu D, Zibaei S, Namjoo M, Jajarmi A. A nonstandard finite difference scheme for the modelling and nonidentical synchronization of a novel fractional chaotic system. Adv Differ Equ. 2021;2021:308. 10.1186/s13662-021-03454-1Search in Google Scholar

[17] Shukla PK, Mamun AA. Introduct Dusty Plasma Phys. London, UK: IOP; 2002. 10.1887/075030653XSearch in Google Scholar

[18] Horanyi M, Mendis DA. Trajectories of charged dust grains in the cometary environment. Astrophys J. 1985;294:357–68. 10.1086/163303Search in Google Scholar

[19] Khater AH, Callebaut DK, Seadawy AR. General soliton solutions for nonlinear dispersive waves in convective type instabilities. Phys Scr. 2006;74:384–93. 10.1088/0031-8949/74/3/015Search in Google Scholar

[20] Çlik N, Seadawy AR, Özkan YS, Yasar E. A model of solitary waves in a nonlinear elastic circular rod: Abundant different type exact solutions and conservation laws. Chaos Solitons Fractal. 2021;143:110486. 10.1016/j.chaos.2020.110486Search in Google Scholar

[21] Mendis DA, Rosenberg M. Some aspects of dust-plasma interactions in the cosmic environment. IEEE Trans Plasma Sci. 1992;20(6):929–34. 10.1109/27.199553Search in Google Scholar

[22] Verheest F. Waves and instabilities in dusty space plasmas. Space Sci Rev. 1996;77(3-4):267–302. 10.1007/BF00226225Search in Google Scholar

[23] Feuerbacher B, Willis RF, Fitton B. Electrostatic potential of interstellar grains. Astrophys J. 1973;181:101–14. 10.1086/152033Search in Google Scholar

[24] Fechtig H, Grün E, Morfill G. Micrometeoroids within ten earth radii. Planetary Space Sci. 1979;27(4):511–31. 10.1016/0032-0633(79)90128-4Search in Google Scholar

[25] Havnes O, Goertz CK, Morfill GE, Grün E, Ip W. Dust charges, cloud potential, and instabilities in a dust cloud embedded in a plasma. J Geophys Res Space Phys. 1987;92(A3):2281–7. 10.1029/SP027p0346Search in Google Scholar

[26] Rao NN, Shukla PK, Yu MY. Dust-acoustic waves in dusty plasmas. Planetary Space Sci. 1990;38(4):543–6. 10.1016/0032-0633(90)90147-ISearch in Google Scholar

[27] Mamun AA, Cairns RA, Shukla PK. Effects of vortex-like and non-thermal ion distributions on nonlinear dust-acoustic waves. Phys Plasmas. 1996;3(7):2610–4. 10.1063/1.871973Search in Google Scholar

[28] Akhtar N, Mahmood S, Saleem H. Dust acoustic solitary waves in the presence of hot and cold dust. Phys Lett A. 2007;361(1):126–32. 10.1016/j.physleta.2006.09.017Search in Google Scholar

[29] Sagdeev RZ. Cooperative phenomena and shock waves in collisionless plasmas. Rev Plasma Phys. 1966;4:23. Search in Google Scholar

[30] Juan L, Xu T, Meng XH, Zhang YX, Zhang HQ. Lax pair, Backlund transformation and N-soliton-like solution for a variable-coefficient Gardner equation from nonlinear lattice, plasma physics and ocean dynamics with symbolic computation. J Math Anal Appl. 2007;336(2):1443–55. 10.1016/j.jmaa.2007.03.064Search in Google Scholar

[31] Hirota R. Exact solution of the Korteweg–de Vries equation for multiple collisions of solitons. Phys Rev Lett. 1971;27(18):1456–8. 10.1103/PhysRevLett.27.1192Search in Google Scholar

[32] Zheng Yi M, Hua GS, Zheng CL. Multisoliton excitations for the Kadomtsev-Petviashvili equation. Zeitschrift Fur Naturforschung A. 2006;61(1–2):32–38. 10.1515/zna-2006-1-205Search in Google Scholar

[33] Naher H, Abdullah FA, Akbar MA. New traveling wave solutions of the higher dimensional nonlinear partial differential equation by the exp-function method. J Appl Math. 2012;2012(2):1–7. 10.1155/2012/575387Search in Google Scholar

[34] Seadawy AR, Manafian J. New soliton solution to the longitudinal wave equation in a magneto-electro-elastic circular rod. Results Phys. 2018;8:1158–67. 10.1016/j.rinp.2018.01.062Search in Google Scholar

[35] Zhang S, Xia TC. A further improved extended Fan sub-equation method and its application to the (3+1)-dimensional Kadomstev-Petviashvili equation. Phys Lett A. 2006;356(2):119–23. 10.1016/j.physleta.2006.03.027Search in Google Scholar

[36] Seadawy AR, El-Rashidy K. Traveling wave solutions for some coupled nonlinear evolution equations. Math Comput Modell. 2013;57(5–6):1371–9. 10.1016/j.mcm.2012.11.026Search in Google Scholar

[37] Seadawy AR, Lu D. Ion acoustic solitary wave solutions of three-dimensional nonlinear extended Zakharov-Kuznetsov dynamical equation in a magnetized two-ion-temperature dusty plasma. Results Phys. 2016;6:590–3. 10.1016/j.rinp.2016.08.023Search in Google Scholar

[38] Seadawy AR. The solutions of the Boussinesq and generalized fifth-order KdV equations by using the direct algebraic method. Appl Math Sci. 2012;6(82):4081–90. Search in Google Scholar

[39] Iqbal M, Seadawy AR, Lu D. Construction of solitary wave solutions to the nonlinear modified Kortewege-de Vries dynamical equation in unmagnetized plasma via mathematical methods. Modern Phys Lett A. 2018;33:1850183, 1–13. 10.1142/S0217732318501833Search in Google Scholar

[40] Iqbal M, Seadawy AR, Lu D, Dispersive solitary wave solutions of nonlinear further modified Kortewege-de Vries dynamical equation in a unmagnetized dusty plasma via mathematical methods, Modern Phys Lett A. 2018;33:1850217, 1–19. 10.1142/S0217732318502176Search in Google Scholar

[41] Seadawy AR, Lu D, Iqbal M. Application of mathematical methods on the system of dynamical equations for the ion sound and Langmuir waves. Pramana. 2019;93(1):10. 10.1007/s12043-019-1771-xSearch in Google Scholar

[42] Iqbal M, Seadawy AR, Lu D. Applications of nonlinear longitudinal wave equation in a magneto-electro-elastic circular rod and new solitary wave solutions. Modern Phys Lett B. 2019;33(18):1950210. 10.1142/S0217984919502105Search in Google Scholar

[43] Seadawy AR, Iqbal M, Lu D. Propagation of kink and anti-kink wave solitons for the nonlinear damped modified Korteweg-de Vries equation arising in ion-acoustic wave in an unmagnetized collisional dusty plasma. Physica A Statist Mech Appl. 2020;544:123560. 10.1016/j.physa.2019.123560Search in Google Scholar

[44] Seadawy AR. Stability analysis for Zakharov-Kuznetsov equation of weakly nonlinear ion-acoustic waves in a plasma. Comput Math Appl. 2014;67(1):172–80. 10.1016/j.camwa.2013.11.001Search in Google Scholar

[45] Seadawy AR. Stability analysis for two-dimensional ion-acoustic waves in quantum plasmas. Phys Plasmas. 2014;21(5):052107. 10.1063/1.4875987Search in Google Scholar

[46] Seadawy AR. Approximation solutions of derivative nonlinear Schrodinger equation with computational applications by variational method. Europ Phys J Plus. 2015;130(9):182. 10.1140/epjp/i2015-15182-5Search in Google Scholar

[47] Seadawy AR. Nonlinear wave solutions of the three dimensional Zakharov Kuznetsov Burgers equation in dusty plasma. Phys A Statist Mechanic Appl. 2015;439:124–31, S0378437115006354. 10.1016/j.physa.2015.07.025Search in Google Scholar

[48] Seadawy AR. Ion acoustic solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili-Burgers equation in quantum plasma. Math Meth Appl Sci. 2017;40(5):1598–607. 10.1002/mma.4081Search in Google Scholar

[49] Seadawy AR. Three dimensional nonlinear modified Zakharov Kuznetsov equation of ion acoustic waves in a magnetized plasma. Comput Math Appl. 2016;71(1):201–12. 10.1016/j.camwa.2015.11.006Search in Google Scholar

[50] Seadawy AR. Solitary wave solutions of two-dimensional nonlinear Kadomtsev-Petviashvili dynamic equation in dust-acoustic plasmas. Pramana. 2017;89(3):49. 10.1007/s12043-017-1446-4Search in Google Scholar

[51] Seadawy AR, Wang J. Three-dimensional nonlinear extended Zakharov-Kuznetsov dynamical equation in a magnetized dusty plasma via acoustic solitary wave solutions. Brazilian J Phys. 2019;49(1):67–78. 10.1007/s13538-018-0617-1Search in Google Scholar

[52] Ahmed I, Seadawy AR, Lu D. Kinky breathers, W-shaped and multi-peak solitons interaction in (2+1)-dimensional nonlinear Schrödinger equation with Kerr law of nonlinearity. Europ Phys J Plus. 2019;134(3):120. 10.1140/epjp/i2019-12482-8Search in Google Scholar

[53] Seadawy AR, El-Rashidy K. Water wave solutions of the coupled system Zakharov-Kuznetsov and generalized coupled KdV equations. Scientific World J. 2014;2014:1–6. 10.1155/2014/724759Search in Google Scholar PubMed PubMed Central

[54] Lu D, Seadawy AR, Iqbal M. Mathematical method via construction of traveling and solitary wave solutions of three coupled system of nonlinear partial differential equations and their applications. Results Phys. 2018;11:1161–71. 10.1016/j.rinp.2018.11.014Search in Google Scholar

[55] Lu D, Seadawy AR, Iqbal M. Construction of new solitary wave solutions of generalized Zakharov-Kuznetsov-Benjamin-Bona-Mahony and simplified modified form of Camassa-Holm equations. Open Phys. 2018;16(1):896–909. 10.1515/phys-2018-0111Search in Google Scholar

[56] Seadawy AR, Iqbal M, Lu D. Propagation of long-wave with dissipation and dispersion in nonlinear media via generalized Kadomtsive-Petviashvili modified equal width-Burgers equation. Indian J Phys. 2020;94(5):675–87. 10.1007/s12648-019-01500-zSearch in Google Scholar

[57] El-Rashidy K., Seadawy AR, Althobaiti S, Makhlouf MM. Multiwave, Kinky breathers and multi-peak soliton solutions for the nonlinear Hirota dynamical system. Results Phys. 2020;19:103678. 10.1016/j.rinp.2020.103678Search in Google Scholar

[58] Seadawy AR, Iqbal M, Lu D. Analytical methods via bright-dark solitons and solitary wave solutions of the higher-order nonlinear Schrödinger equation with fourth-order dispersion. Modern Phys Lett B. 2019;33(35):1950443. 10.1142/S0217984919504438Search in Google Scholar

[59] Seadawy AR, Lu D, Yue C. Travelling wave solutions of the generalized nonlinear fifth-order KdV water wave equations and its stability. J Taibah Univ Sci. 2017;11:623–33. 10.1016/j.jtusci.2016.06.002Search in Google Scholar

[60] Seadawy AR, Iqbal M, Lu D. Applications of propagation of long-wave with dissipation and dispersion in nonlinear media via solitary wave solutions of generalized Kadomtsev-Petviashvili modified equal width dynamical equation. Comput Math Appl. 2019;78:3620–32. 10.1016/j.camwa.2019.06.013Search in Google Scholar

[61] Iqbal M, Seadawy AR, Lu D, Xia X. Construction of bright-dark solitons and ion-acoustic solitary wave solutions of dynamical system of nonlinear wave propagation. Modern Phys Lett A. 2019;34:1950309. 10.1142/S0217732319503097Search in Google Scholar

[62] Seadawy AR, Iqbal M, Lu D. Nonlinear wave solutions of the Kudryashov-Sinelshchikov dynamical equation in mixtures liquid-gas bubbles under the consideration of heat transfer and viscosity. J Taibah Univ Sci. 2019;13(1):1060–72. 10.1080/16583655.2019.1680170Search in Google Scholar

[63] Iqbal M, Seadawy AR, Khalil OH, Lu D. Propagation of long internal waves in density stratified ocean for the (2+1)-dimensional nonlinear Nizhnik-Novikov-Vesselov dynamical equation. Results Phys. 2020;16:102838. 10.1016/j.rinp.2019.102838Search in Google Scholar

[64] Khater, AH, Helal MA, Seadawy AR. General soliton solutions of n-dimensional nonlinear Schrödinger equation. IL Nuovo Cimento. 2000;115B:1303–12. Search in Google Scholar

[65] Schamel H. Stationary solitary, snoidal and sinusoidal ion acoustic waves. Plasma Phys. 1972;14(10):905. 10.1088/0032-1028/14/10/002Search in Google Scholar

[66] Schamel H. A modified Korteweg-de Vries equation for ion acoustic wavess due to resonant electrons. J Plasma Phys. 1973;9(3):377–87. 10.1017/S002237780000756XSearch in Google Scholar

Received: 2021-12-01

Revised: 2022-03-13

Accepted: 2022-03-29

Published Online: 2022-05-06

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-0030

Keywords for this article

modified mathematical method; Kadomtsev–Petviashvili equation; potential KP equation; Gardner KP equation; ion-acoustic solitary wave solutions; solitons

Creative Commons

BY 4.0