Home Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model
Article Open Access

Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model

  • Dexu Zhao , Dianchen Lu , Samir A. Salama and Mostafa M. A. Khater EMAIL logo
Published/Copyright: December 1, 2021

Abstract

This article investigates the dynamical and physical behavior of the second positive member in a new, utterly integrable hierarchy. This investigation depends on constructing novel analytical and approximate solutions to the Qiao model. The model’s name is after the researcher who derived the mathematical formula of it in 2007. This model possesses a Lax representation and bi-Hamiltonian structure. This study employs the unified and variational iteration (VI) method to obtain analytical and numerical solutions to the considered model. The obtained analytical solutions are used to calculate the necessary conditions for applying the suggested numerical method that makes checking the obtained solutions’ accuracy a valuable option. The obtained solutions are sketched in different techniques to explain more physical and dynamics details of the Qiao model and show the matching between obtained solutions.

1 Introduction

Recently, many distinct fields such as soil mechanics, nonlinear optics, plasma wave, condensed matter physics, wave propagation in plasma physics, fluid mechanics, population ecology, neural networks, infectious disease epidemiology, thermodynamics, solid-state physics, civil engineering, quantum mechanics, and so on have considered nonlinear partial differential equation (NLPDE) field as a keyword for representing and explaining many of complex phenomena [1,2,3, 4,5,6]. Significantly, this field is able to formulate many complex natural phenomena in mathematical formulas [7,8,9]. The formats of these mathematical models depend on the real experiments to determine the parameters and empirical functions [10,11]. Consequently, these complex phenomena can be explained and simplified, which allows using the models’ properties for improving their applications [12].

Subsequently, a lot of mathematicians and physicists have focused on studying these models by applying some novel and accurate analytical and approximate schemes for evaluating computational solutions [13,14]. Furthermore, some other researchers have been focusing on deriving these schemes, which depends on the following three main icons: ansatz techniques, reduction techniques, and Lie group techniques [15,16,17].

This article handles the second positive member in a new utterly integrable hierarchy using the unified and VI methods to construct novel solitary wave solutions. This model is given by ref. [18]

(1) 2 U t 1 U q x x x + 1 U q x = 0 .

Using the next wave transformation U ( x , t ) = u ( ζ ) , ζ = x + ω t and replacing q with r , convert equation (1) into

(2) 2 ω u ( u r ) + ( u r ) = 0 .

Replacing u with v 1 r , integrating equation (2) once with respect to ζ , and zero constant of integration, we obtain

(3) 2 ω v λ v + v = 0 ,

where λ = 1 r . Balancing the terms of equation (3), we obtain n = 2 λ 1 . Thus, we apply another transformation on equation (3), which is given by v = ν 2 λ 1 . This transformation converts equation (3) into the following form:

(4) 2 ν 4 ω + ν 2 + 2 ( λ 3 ) ( ν ) 2 2 ( λ 1 ) ν ν ( λ 1 ) 2 = 0 .

Balancing the highest order and nonlinear term in equation (4), we obtain n = 1 . Consequently, the general solutions of the Qiao model based on the unified method take the following form:

(5) ν ( ζ ) = a 0 + i = 1 n i 1 ( ζ ) [ a i ( ζ ) + b i G ( ζ ) ] = a 0 + a 1 ( ζ ) + b 1 G ( ζ ) ,

where a 0 , a 1 , b 1 are arbitrary constants to be determined through the method’s framework. While ( ζ ) , G ( ζ ) satisfies the following ordinary differential equations:

(6) ( ζ ) = β 1 ( ζ ) G ( ζ ) , G ( ζ ) = β 2 β 3 ( ζ ) + β 1 G ( ζ ) 2 , G ( ξ ) 2 = β 2 + ( β 3 2 1 ) ( ξ ) 2 β 2 2 β 3 ( ξ ) , when β 1 = 1 , β 4 = 1 , β 2 0 Case I . G ( ξ ) 2 = β 2 + ( β 3 2 + 1 ) ( ξ ) 2 β 2 2 β 3 ( ξ ) , when β 1 = 1 , β 4 = 1 , β 2 0 Case II . G ( ζ ) 2 = β 2 ( β 3 2 1 ) ( ζ ) 2 β 2 + 2 β 3 ( ζ ) , when β 1 = 1 , β 4 = 1 , β 2 0 Case III ,

where β 2 , β 3 are arbitrary constants to be calculated through the method’s framework.

The rest of the article is organized in the following order: Section 2 investigates analytical and approximate solutions and checks the analytical solutions’ stability characterization. Section 3 demonstrates the results’ novelty and accuracy and the article’s contributions. Section 4 explains the figures’ interpretations. Section 5 gives the conclusion of the whole article.

2 Accurate stable soliton wave solutions of the Qiao model

Here, we apply the unified method, the Hamiltonian system’s characterizations, and the VI method to the nonlinear Qiao model to construct novel, accurate, stable analytical solutions. This investigation is given in the following subsections.

2.1 Abundant wave solutions of Qiao model

Employing the unified method’s framework to equation (5) for constructing novel soliton wave solutions, we obtain the following values of the aforementioned parameters:

Case I: when [ β 1 = 1 & β 4 = 1 ].

a 0 = 0 , a 1 = β 3 2 1 4 2 ω , b 1 = 0 , λ = 3 , β 2 = 4 . a 0 = i 2 2 ω , a 1 = β 3 2 1 8 2 ω , b 1 = i 4 2 ω , λ = 3 , β 2 = 4 .

Case II: when [ β 1 = 1 & β 4 = 1 ].

a 0 = 0 , a 1 = β 3 2 1 4 2 ω , b 1 = 0 , λ = 3 , β 2 = 4 . a 0 = i 2 2 ω , a 1 = β 3 2 1 8 2 ω , b 1 = i 4 2 ω , λ = 3 , β 2 = 4 .

Case III: when [ β 1 = 1 & β 4 = 1 ].

a 0 = 0 , a 1 = 1 β 3 2 4 2 ω , b 1 = 0 , λ = 3 , β 2 = 4 . a 0 = i 2 2 ω , a 1 = 1 β 3 2 8 2 ω , b 1 = i 4 2 ω , λ = 3 , β 2 = 4 .

Consequently, the solitary wave solutions of the Qiao model are given in the following formulas:

(7) U 1 ( x , t ) = 1 8 2 3 / 4 β 2 β 3 2 1 ω ( β 3 + cosh ( β 2 ( t ω + x ) ) ) 3 / 2 .

(8) U 2 ( x , t ) = 1 32 2 4 β 3 2 1 β 2 + 2 i β 2 sinh ( β 2 ( t ω + x ) ) + 4 i ( β 3 + cosh ( β 2 ( t ω + x ) ) ) ω ( β 3 + cosh ( β 2 ( t ω + x ) ) ) 3 / 2 .

(9) U 3 ( x , t ) = 1 8 2 3 / 4 β 2 β 3 2 1 ω ( β 3 + sinh ( β 2 ( t ω + x ) ) ) 3 / 2 .

(10) U 4 ( x , t ) = 1 32 2 4 β 3 2 1 β 2 + 4 i ( β 3 + sinh ( β 2 ( t ω + x ) ) ) + 2 i β 2 cosh ( β 2 ( t ω + x ) ) ω ( β 3 + sinh ( β 2 ( t ω + x ) ) ) 3 / 2 .

(11) U 5 ( x , t ) = 1 8 2 3 / 4 β 2 1 β 3 2 ω ( β 3 + cos ( β 2 ( t ω + x ) ) ) 3 / 2 .

(12) U 6 ( x , t ) = 1 32 2 4 1 β 3 2 β 2 + 2 i β 2 sin ( β 2 ( t ω + x ) ) + 4 i ( β 3 + cos ( β 2 ( t ω + x ) ) ) ω ( β 3 + cos ( β 2 ( t ω + x ) ) ) 3 / 2 .

(13) U 7 ( x , t ) = 1 8 2 3 / 4 β 2 1 β 3 2 ω ( β 3 + sin ( β 2 ( t ω + x ) ) ) 3 / 2 .

(14) U 8 ( x , t ) = 1 32 2 4 1 β 3 2 β 2 + 4 i ( β 3 + sin ( β 2 ( t ω + x ) ) ) 2 i β 2 cos ( β 2 ( t ω + x ) ) ω ( β 3 + sin ( β 2 ( t ω + x ) ) ) 3 / 2 .

2.2 Stability

Studying the stability property of the Qiao model’s solution (1) through the Hamiltonian system’s characterizations with the following values of the constants [ β 2 = 4 , β 3 = 6 ] , we obtain:

T = 1 5,488,000 2 ω 5 / 2 [ i ( 3,430,000 ω log ( 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) ) + 630 37 log ( 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) ) + 264,600 log ( 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) ) 1,260 35 tanh 1 6 + e 20 ( ω + 1 ) 35 14,700 1,295 × tanh 1 6 + e 20 ( ω + 1 ) 35 1,260 35 tanh 1 6 + e 2 ω + 2 35 14,700 1,295 tanh 1 6 + e 2 ω + 2 35 5,040 35 ω × tanh 1 6 + e 2 ω + 20 35 19,892 1,295 ω tanh 1 6 + e 2 ω + 20 35 50,400 35 ω tanh 1 6 + e 20 ω + 2 35 198,920 1,295 ω tanh 1 6 + e 20 ω + 2 35 343,000 ω log ( 1 + 12 e 2 ω + 2 + e 4 ω + 4 ) + 264,600 log ( 1 + 12 e 2 ω + 2 + e 4 ω + 4 ) + 343,000 ω log ( 1 + 12 e 2 ω + 20 + e 4 ω + 40 ) 630 37 log ( 1 + 12 e 2 ω + 20 + e 4 ω + 40 ) 264,600 log ( 1 + 12 e 2 ω + 20 + e 4 ω + 40 ) + 3,430,000 ω log ( 1 + 12 e 20 ω + 2 + e 40 ω + 4 ) 630 37 log ( 1 + 12 e 20 ω + 2 + e 40 ω + 4 ) 264,600 log ( 1 + 12 e 20 ω + 2 + e 40 ω + 4 ) + 343,000 ω log 1 + e 2 ω + 2 6 35 2,520 35 ω log 1 + e 2 ω + 20 6 35 9,946 1,295 ω log 1 + e 2 ω + 20 6 35 343,000 ω log 1 + e 2 ω + 20 6 35 25,200 35 ω log 1 + e 20 ω + 2 6 35

99,460 1,295 ω log 1 + e 20 ω + 2 6 35 3,430,000 ω log 1 + e 20 ω + 2 6 35 + 3,430,000 ω log 1 e 20 ( ω + 1 ) 35 6 25,200 35 ω log 1 + e 20 ( ω + 1 ) 35 + 6 99,460 1,295 ω log 1 + e 20 ( ω + 1 ) 35 + 6 + 3,430,000 ω log 1 + e 20 ( ω + 1 ) 35 + 6 2,520 35 ω log 1 + e 2 ω + 2 35 + 6 9,946 1,295 ω log 1 + e 2 ω + 2 35 + 6 + 343,000 ω log 1 + e 2 ω + 2 35 + 6 343,000 ω log 1 + e 2 ω + 20 35 + 6 3,430,000 ω log 1 + e 20 ω + 2 35 + 6 + 171,500 Li 2 e 20 ( ω + 1 ) 35 6 + 171,500 × Li 2 e 2 ω + 2 35 6 1,260 35 Li 2 e 2 ω + 20 35 6 4,973 1,295 Li 2 e 2 ω + 20 35 6 171,500 Li 2 e 2 ω + 20 35 6 1,260 35 Li 2 e 20 ω + 2 35 6 4,973 1,295 Li 2 e 20 ω + 2 35 6 171,500 Li 2 e 20 ω + 2 35 6 1,260 35 × Li 2 e 20 ( ω + 1 ) 35 + 6 4,973 1,295 Li 2 e 20 ( ω + 1 ) 35 + 6 + 171,500 Li 2 e 20 ( ω + 1 ) 35 + 6 1,260 35 Li 2 e 2 ω + 2 35 + 6 4,973 1,295 Li 2 e 2 ω + 2 35 + 6 + 171,500 Li 2 e 2 ω + 2 35 + 6 171,500 Li 2 e 2 ω + 20 35 + 6 171,500

× Li 2 e 20 ω + 2 35 + 6 + 50,400 ω tanh 1 6 + e 20 ( ω + 1 ) 35 35 + 5,040 ω tanh 1 6 + e 2 ω + 2 35 35 + 1,260 × tanh 1 6 + e 2 ω + 20 35 35 + 1,260 tanh 1 6 + e 20 ω + 2 35 35 + 2,520 ω log 1 + e 2 ω + 2 6 35 35 + 25,200 ω log 1 e 20 ( ω + 1 ) 35 6 35 + 2,520 ω log 1 + e 2 ω + 20 35 + 6 35 + 25,200 ω log 1 + e 20 ω + 2 35 + 6 35 + 1,260 Li 2 e 20 ( ω + 1 ) 35 6 35 + 1,260 Li 2 e 2 ω + 2 35 6 35 + 1,260 Li 2 e 2 ω + 20 35 + 6 35 + 1,260 Li 2 e 20 ω + 2 35 + 6 × 35 + 630 log ( 1 + 12 e 2 ω + 2 + e 4 ω + 4 ) 37 + 173,950 e 2 ω + 20 37 1 + 12 e 2 ω + 20 + e 4 ω + 40 + 173,950 e 20 ω + 2 37 1 + 12 e 20 ω + 2 + e 40 ω + 4 + 198,920 ω × tanh 1 6 + e 20 ( ω + 1 ) 35 1,295 + 19,892 ω tanh 1 6 + e 2 ω + 2 35 1,295 + 14,700 tanh 1 6 + e 2 ω + 20 35 1,295 + 14,700 tanh 1 6 + e 20 ω + 2 35 1,295 + 9,946 ω log 1 + e 2 ω + 2 6 35 1,295 + 99,460 ω log 1 e 20 ( ω + 1 ) 35 6 1,295 + 9,946 ω log 1 + e 2 ω + 20 35 + 6 1,295 + 99,460 ω log 1 + e 20 ω + 2 35 + 6 1,295 + 4,973 Li 2 e 20 ( ω + 1 ) 35 6 1,295

(15) + 4,973 Li 2 e 2 ω + 2 35 6 1,295 + 4,973 Li 2 e 2 ω + 20 35 + 6 1,295 + 4,973 Li 2 e 20 ω + 2 35 + 6 1,295 173,950 37 e 20 ( ω + 1 ) 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) + 1,073,100 e 20 ( ω + 1 ) 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) + 178,850 1 + 12 e 20 ( ω + 1 ) + e 40 ( ω + 1 ) 173,950 37 e 2 ω + 2 1 + 12 e 2 ω + 2 + e 4 ω + 4 + 1,073,100 e 2 ω + 2 1 + 12 e 2 ω + 2 + e 4 ω + 4 + 178,850 1 + 12 e 2 ω + 2 + e 4 ω + 4 1,073,100 e 2 ω + 20 1 + 12 e 2 ω + 20 + e 4 ω + 40 178,850 1 + 12 e 2 ω + 20 + e 4 ω + 40 1,073,100 e 20 ω + 2 1 + 12 e 20 ω + 2 + e 40 ω + 4 178,850 1 + 12 e 20 ω + 2 + e 40 ω + 4 ,

where T is the Hamiltonian system’s momentum. Thus, we have

(16) T ω ω = 9 = 0.00142019356 > 0 .

Consequently, the investigated solution (8) is the stable solution in the following interval x [ 1 , 10 ] , t [ 1 , 10 ] . Employing the aforementioned steps for the other obtained solutions demonstrates their stability characterization.

2.3 Approximate solutions

Here, the approximate solution of the investigated model is constructed by applying the variational iterational (VI) method. The method’s head lines are given as follows. The NLPDE can be given in the following formula:

(17) U ( x , t ) + N U ( x , t ) = G ( x , t ) ,

where , N , G ( x , t ) represent, respectively, a linear operator, a nonlinear operator, and a known differential function. The correction functional in the VI method’s framework is constructed in the next formula:

(18) U n + 1 ( x , t ) = U n ( x , t ) + 0 t λ { L U n ( x , s ) + N U ˜ n ( x , s ) G ( x , s ) } d s ,

where λ , U n ( x , t ) , U ˜ n ( x , t ) are, respectively, a general Lagrange multiplier that can be identified optimally via variational theory, the n th approximate solution, and considered a restricted variation. On the other hand, this term 0 t λ { L U n ( x , s ) + N U ˜ n ( x , s ) G ( x , s ) } d s is called the correction.

Handling Qiao model by the VIM, we obtain the following approximate solution:

(19) U Approx. ( x , t ) = 1 32 2 4 [ 2 Ξ sech 2 ( x ) + ( sinh ( x ) Ξ cosh ( x ) Ξ ) sech ( x ) + ( cosh ( x ) Ξ + sinh ( x ) Ξ ) sech ( x ) + 2 Ξ + 1 ( Ξ ) 5 [ sech 7 ( x ) ( 768 t cosh ( 5 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 + 768 t sinh ( 5 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ) ] + sech 7 ( x ) ( Ξ ) 5 [ 576 t cosh ( 5 x ) 2 Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 + 576 t sinh ( 5 x ) 2 Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ] + sech 7 ( x ) ( Ξ ) 5 [ 768 t cosh ( 3 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 768 t sinh ( 3 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ] + sech 7 ( x ) ( Ξ ) 5 [ 576 2 t sinh ( 3 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 576 t cosh ( 3 x ) 2 Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ] + sech 6 ( x ) ( Ξ ) 5 [ 1,920 t cosh ( 4 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 1,920 t sinh ( 4 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ] + sech 6 ( x ) ( Ξ ) 5 [ 768 2 t sinh ( 4 x ) Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 768 t cosh ( 4 x ) 2 Ξ ( ( Ξ ) 3 / 2 ) 2 / 3 ] ] + ,

where Ξ = 2 sech 2 ( x ) + 2 tanh ( x ) + 2 .

3 Results and discussion

Here, we discuss the obtained solutions in detail to show the article’s contributions. Also, we compare our obtained results with recently published articles that have investigated the same studied model. Finally, we show our obtained solutions’ accuracy by studying the above-shown tables. We have obtained some solutions in different formats such as hyperbolic and trigonometric. These solutions represent bright soliton (7), dark soliton (8), formal soliton (9), singular dark soliton (10), and periodic wave solutions (11), (12), (13), and (14).

Second, the result’s novelty can be explained by comparing our solutions with those that have been published in ref. [19], where Pan and Liu have derived the traveling wave systems of the Qiao model. In that article, they have obtained just two solutions (31) and (32). Both solutions [(31), (32) in [19]] are similar to our solutions (7) and (9) when p + q = 3 2 , ω = 1 1,296 c 2 p 2 ( p + q ) 2 . All our other solutions are novel solutions.

Finally, the results, accuracy can be tested by studying Tables 1, and 2, where these tables show the value of approximate, analytical, and absolute errors between both of them. These two tables show the comparison of both solutions, which refers to the solutions’ accuracy. Furthermore, Figure 5 shows the absolute error between analytical and approximate solutions when t { 2 , 4 , 6 , 8 } in two different zooms for more explanation of the table’s values.

Table 1

Analytical and approximate solutions of the investigated model through different values of t { 2 , 4 , 6 , 8 , 10 }

t = 2 t = 4 t = 6 t = 8 t = 10
Ex. Approx. Ex. Approx. Ex. Approx. Ex. Approx. Ex. Approx.
0 2.98 × 1 0 11 1.9826002 1.12 × 1 0 21 4.130979 0 6.2793579 0 8.4277367 0 10.576116
2 1.2 × 1 0 8 0.2314726 4.53 × 1 0 19 0.2505121 0 0.2695517 0 0.2885913 0 0.3076309
4 4.84 × 1 0 6 0.2116478 1.83 × 1 0 16 0.2130278 0 0.2144078 0 0.2157878 0 0.2171677
6 0.001872 0.2102506 7.38 × 1 0 14 0.2102763 2.94 × 1 0 24 0.2103019 0 0.2103276 0 0.2103533
8 0.165779 0.2102246 2.98 × 1 0 11 0.2102251 1.12 × 1 0 21 0.2102255 0 0.210226 0 0.2102265
10 0.212433 0.2102241 1.2 × 1 0 8 0.2102241 4.53 × 1 0 19 0.2102241 0 0.2102241 0 0.2102241
12 0.210268 0.2102241 4.84 × 1 0 6 0.2102241 1.83 × 1 0 16 0.2102241 0 0.2102241 0 0.2102241
14 0.210225 0.2102241 0.001872 0.2102241 7.38 × 1 0 14 0.2102241 2.94 × 1 0 24 0.2102241 0 0.2102241
16 0.210224 0.2102241 0.165779 0.2102241 2.98 × 1 0 11 0.2102241 1.12 × 1 0 21 0.2102241 0 0.2102241
18 0.210224 0.2102241 0.212433 0.2102241 1.2 × 1 0 8 0.2102241 4.53 × 1 0 19 0.2102241 0 0.2102241
20 0.210224 0.2102241 0.210268 0.2102241 4.84 × 1 0 6 0.2102241 1.83 × 1 0 16 0.2102241 0 0.2102241
22 0.210224 0.2102241 0.210225 0.2102241 0.001872 0.2102241 7.38 × 1 0 14 0.2102241 2.94 × 1 0 24 0.2102241
24 0.210224 0.2102241 0.210224 0.2102241 0.165779 0.2102241 2.98 × 1 0 11 0.2102241 1.12 × 1 0 21 0.2102241
26 0.210224 0.2102241 0.210224 0.2102241 0.212433 0.2102241 1.2 × 1 0 8 0.2102241 4.53 × 1 0 19 0.2102241
28 0.210224 0.2102241 0.210224 0.2102241 0.210268 0.2102241 4.84 × 1 0 6 0.2102241 1.83 × 1 0 16 0.2102241
30 0.210224 0.2102241 0.210224 0.2102241 0.210225 0.2102241 0.001872 0.2102241 7.38 × 1 0 14 0.2102241
32 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.165779 0.2102241 2.98 × 1 0 11 0.2102241
34 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.212433 0.2102241 1.2 × 1 0 8 0.2102241
36 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210268 0.2102241 4.84 × 1 0 6 0.2102241
38 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210225 0.2102241 0.001872 0.2102241
40 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.165779 0.2102241
42 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.212433 0.2102241
44 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210268 0.2102241
46 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210225 0.2102241
48 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
50 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
52 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
54 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
56 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
58 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
60 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241 0.210224 0.2102241
Table 2

Absolute values of error between analytical and approximate solutions that have been obtained by unified method and VI method through different values of t { 2 , 4 , 6 , 8 , 10 }

t = 2 t = 4 t = 6 t = 8 t = 10
0 1.982600202 4.130979035 6.279357868 8.4277367 10.57611553
2 0.231472548 0.250512137 0.269551713 0.28859129 0.307630867
4 0.211642993 0.213027812 0.214407789 0.215787767 0.217167744
6 0.208378366 0.210276254 0.210301927 0.210327601 0.210353275
8 0.044445958 0.210225059 0.21022553 0.210226 0.21022647
10 0.00220887 0.210224109 0.21022413 0.210224139 0.210224147
12 4.37535 × 1 0 5 0.210219263 0.210224104 0.210224104 0.210224105
14 8.02511 × 1 0 7 0.20835189 0.210224104 0.210224104 0.210224104
16 1.46989 × 1 0 8 0.044445473 0.210224104 0.210224104 0.210224104
18 2.6922 × 1 0 10 0.002208879 0.210224092 0.210224104 0.210224104
20 4.93092 × 1 0 12 4.37536 × 1 0 5 0.210219262 0.210224104 0.210224104
22 9.02611 × 1 0 14 8.02514 × 1 0 7 0.20835189 0.210224104 0.210224104
24 1.63758 × 1 0 15 1.46989 × 1 0 8 0.044445473 0.210224104 0.210224104
26 5.55112 × 1 0 17 2.69221 × 1 0 10 0.002208879 0.210224092 0.210224104
28 2.77556 × 1 0 17 4.93092 × 1 0 12 4.37536 × 1 0 5 0.210219262 0.210224104
30 2.77556 × 1 0 17 9.02611 × 1 0 14 8.02514 × 1 0 7 0.20835189 0.210224104
32 2.77556 × 1 0 17 1.63758 × 1 0 15 1.46989 × 1 0 8 0.044445473 0.210224104
34 2.77556 × 1 0 17 5.55112 × 1 0 17 2.69221 × 1 0 10 0.002208879 0.210224092
36 5.55112 × 1 0 17 2.77556 × 1 0 17 4.93092 × 1 0 12 4.37536 × 1 0 5 0.210219262
38 5.55112 × 1 0 17 2.77556 × 1 0 17 9.02611 × 1 0 14 8.02514 × 1 0 7 0.20835189
40 2.77556 × 1 0 17 2.77556 × 1 0 17 1.63758 × 1 0 15 1.46989 × 1 0 8 0.044445473
42 2.77556 × 1 0 17 2.77556 × 1 0 17 5.55112 × 1 0 17 2.69221 × 1 0 10 0.002208879
44 5.55112 × 1 0 17 5.55112 × 1 0 17 2.77556 × 1 0 17 4.93092 × 1 0 12 4.37536 × 1 0 5
46 2.77556 × 1 0 17 5.55112 × 1 0 17 2.77556 × 1 0 17 9.02611 × 1 0 14 8.02514 × 1 0 7
48 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 1.63758 × 1 0 15 1.46989 × 1 0 8
50 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 5.55112 × 1 0 17 2.69221 × 1 0 10
52 2.77556 × 1 0 17 5.55112 × 1 0 17 5.55112 × 1 0 17 2.77556 × 1 0 17 4.93092 × 1 0 12
54 2.77556 × 1 0 17 2.77556 × 1 0 17 5.55112 × 1 0 17 2.77556 × 1 0 17 9.02611 × 1 0 14
56 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 1.63758 × 1 0 15
58 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 2.77556 × 1 0 17 5.55112 × 1 0 17
60 2.77556 × 1 0 17 2.77556 × 1 0 17 5.55112 × 1 0 17 5.55112 × 1 0 17 2.77556 × 1 0 17

4 Figures’ interpretation

This part is responsible for explaining the above-shown sketches and their physical interpretation of each. Generally, these figures illustrate the perspective view, the wave propagation pattern along x -axis, and the overhead view of the tested solution. This explanation is given in the following form.

  • Figure 1 explains bright soliton wave of equation (7) when β 2 = 4 , β 3 = 2 in three-dimensional (3D), two-dimensional (2D), and contour plots.

  • Figure 2 clarifies dark soliton wave of equation (8) when β 2 = 4 , β 3 = 6 in 3D, 2D, and contour plots.

  • Figure 3 shows combined formal soliton wave of equation (19) in 3-D, 2-D, and contour plots.

  • Figure 4 demonstrates the matching between analytical and approximate solutions based on the data shown in Table 1.

  • Figure 5 represents the absolute error through the values shown in Table 2 in 2D plot.

Figure 1 
               3D (a), 2D (b), and density (c) elucidations of equation (7).
Figure 1

3D (a), 2D (b), and density (c) elucidations of equation (7).

Figure 2 
               3D (a), 2D (b), and density (c) elucidations of equation (8).
Figure 2

3D (a), 2D (b), and density (c) elucidations of equation (8).

Figure 3 
               3D (a), 2D (b), and density (c) elucidations of equation (19).
Figure 3

3D (a), 2D (b), and density (c) elucidations of equation (19).

Figure 4 
               (a) 2D, (b) distribution plots for the obtained analytical wave solutions with different values of 
                     
                        
                        
                           t
                        
                        t
                     
                  , (c) 2D, (d) distribution/Rug plots for the obtained approximate wave solutions with different values of 
                     
                        
                        
                           t
                        
                        t
                     
                  , (e), (f) 2D plots of analytical and approximate in two different zoom when 
                     
                        
                        
                           t
                           ∈
                           
                              {
                              
                                 2
                                 ,
                                 4
                                 ,
                                 6
                                 ,
                                 8
                                 ,
                                 10
                              
                              }
                           
                        
                        t\in \left\{2,4,6,8,10\right\}
                     
                  .
Figure 4

(a) 2D, (b) distribution plots for the obtained analytical wave solutions with different values of t , (c) 2D, (d) distribution/Rug plots for the obtained approximate wave solutions with different values of t , (e), (f) 2D plots of analytical and approximate in two different zoom when t { 2 , 4 , 6 , 8 , 10 } .

Figure 5 
               (a) Zoom 2D plot for the absolute error and (b) regular 2D plot for the absolute error based on the values shown in Table 2.
Figure 5

(a) Zoom 2D plot for the absolute error and (b) regular 2D plot for the absolute error based on the values shown in Table 2.

5 Conclusion

This article has successfully investigated the analytical and approximate solutions of the Qiao model. The unified and VI methods have been applied, and many novel analytical solutions have been obtained. Additionally, the approximate solution has been obtained. The obtained solutions have been demonstrated through some different figures in 2D, 3D, and contour plots. The stability property of the results has been checked by using the Hamiltonian system’s properties. Finally, all obtained solutions have been checked by putting them back into the original model using Mathematica 12.

  1. Funding information: The authors greatly thank Taif University for providing fund for this work through Taif University Researchers Supporting Project number (TURSP-2020/52), Taif University, Taif, Saudi Arabia.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: This work does not have any conflicts of interest.

  4. Data availability statement: The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

[1] Rezazadeh H , Kumar D , Neirameh A , Eslami M , Mirzazadeh M . Applications of three methods for obtaining optical soliton solutions for the Lakshmanan-Porsezian-Daniel model with kerr law nonlinearity. Pramana. 2020;94(1):39. 10.1007/s12043-019-1881-5Search in Google Scholar

[2] Mirhosseini-Alizamini SM , Rezazadeh H , Eslami M , Mirzazadeh M , Korkmaz A . New extended direct algebraic method for the Tzitzica type evolution equations arising in nonlinear optics. Comput Meth Differ Equ. 2020;8(1):28–53. Search in Google Scholar

[3] Souleymanou A , Korkmaz A , Rezazadeh H , Mukam SPT , Bekir A . Soliton solutions in different classes for the Kaup-Newell model equation. Modern Phys Lett B. 2020;34(3):2050038. 10.1142/S0217984920500384Search in Google Scholar

[4] Attia RA , Lu D , Ak T , Khater MM . Optical wave solutions of the higher-order nonlinear Schrödinger equation with the non-kerr nonlinear term via modified Khater method. Modern Phys Lett B. 2020;34(5):2050044. 10.1142/S021798492050044XSearch in Google Scholar

[5] Khater MM , Park C , Abdel-Aty A-H , Attia RA , Lu D . On new computational and numerical solutions of the modified Zakharov-Kuznetsov equation arising in electrical engineering. Alex Eng J. 2020;59(3):1099–105. 10.1016/j.aej.2019.12.043Search in Google Scholar

[6] Wang K-J . Generalized variational principle and periodic wave solution to the modified equal width-burgers equation in nonlinear dispersion media. Phys Lett A. 2021;419:127723. 10.1016/j.physleta.2021.127723Search in Google Scholar

[7] Khater MM , Attia RA , Abdel-Aty A-H , Abdou M , Eleuch H , Lu D . Analytical and semi-analytical ample solutions of the higher-order nonlinear Schrödinger equation with the non-Kerr nonlinear term. Result Phys. 2020;16:103000. 10.1016/j.rinp.2020.103000Search in Google Scholar

[8] Park C , Khater MM , Attia RA , Alharbi W , Alodhaibi SS . An explicit plethora of solution for the fractional nonlinear model of the low-pass electrical transmission lines via Atangana-Baleanu derivative operator. Alex Eng J. 2020;59(3):1205–14. 10.1016/j.aej.2020.01.044Search in Google Scholar

[9] Wang K-J , Wang G-D . Variational theory and new abundant solutions to the (1.2)-dimensional chiral nonlinear Schrödinger equation in optics. Phys Lett A. 2021;412:127588. 10.1016/j.physleta.2021.127588Search in Google Scholar

[10] Khater MM , Attia RA , Baleanu D . Abundant new solutions of the transmission of nerve impulses of an excitable system. Eur Phys J Plus. 2020;135(2):1–12. 10.1140/epjp/s13360-020-00261-7Search in Google Scholar

[11] Yue C , Khater MM , Attia RA , Lu D . The plethora of explicit solutions of the fractional KS equation through liquid-gas bubbles mix under the thermodynamic conditions via Atangana-Baleanu derivative operator. Adv Differ Equ. 2020;2020(1):1–12. 10.1186/s13662-020-2540-3Search in Google Scholar

[12] Khater MM , Park C , Lu D , Attia RA . Analytical, semi-analytical, and numerical solutions for the Cahn-Allen equation. Adv Differ Equ. 2020;2020(1):1–12. 10.1186/s13662-019-2475-8Search in Google Scholar

[13] Qin H , Khater M , Attia RA , Lu D . Approximate Simulations for the Non-linear Long-Short Wave Interaction System. Front Phys. 2020;7:230. 10.3389/fphy.2019.00230Search in Google Scholar

[14] Khater M , Attia RA , Lu D . Computational and numerical simulations for the nonlinear fractional Kolmogorov-Petrovskii-Piskunov (FKPP) equation. Phys Scripta. 2020;95(5):055213. 10.1088/1402-4896/ab76f8Search in Google Scholar

[15] Qin H , Attia RA , Khater M , Lu D . Ample soliton waves for the crystal lattice formation of the conformable time-fractional (n+1) Sinh-Gordon equation by the modified Khater method and the Painlevé property. J Intell Fuzzy Syst. 2020;38(3):2745–52. 10.3233/JIFS-179560Search in Google Scholar

[16] Khater MM , Lu D , Attia RA . Dispersive long wave of nonlinear fractional Wu-Zhang system via a modified auxiliary equation method. AIP Adv. 2019;9(2):025003. 10.1063/1.5087647Search in Google Scholar

[17] Sun Y , Jayaraman AS , Chirikjian GS . Lie group solutions of advection-diffusion equations. Phys Fluids. 2021;33(4):046604. 10.1063/5.0048467Search in Google Scholar

[18] Qiao Z . New integrable hierarchy, its parametric solutions, cuspons, one-peak solitons, and M/W-shape peak solitons. J Math Phys. 2007;48(8):082701. 10.1063/1.2759830Search in Google Scholar

[19] Pan C , Liu Z . Further results on the traveling wave solutions for an integrable equation. J Appl Math. 2013;2013:681383. 10.1155/2013/681383Search in Google Scholar

Received: 2021-10-14
Revised: 2021-10-22
Accepted: 2021-10-24
Published Online: 2021-12-01

© 2021 Dexu Zhao et al., published by De Gruyter

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

Articles in the same Issue

  1. Regular Articles
  2. Circular Rydberg states of helium atoms or helium-like ions in a high-frequency laser field
  3. Closed-form solutions and conservation laws of a generalized Hirota–Satsuma coupled KdV system of fluid mechanics
  4. W-Chirped optical solitons and modulation instability analysis of Chen–Lee–Liu equation in optical monomode fibres
  5. The problem of a hydrogen atom in a cavity: Oscillator representation solution versus analytic solution
  6. An analytical model for the Maxwell radiation field in an axially symmetric galaxy
  7. Utilization of updated version of heat flux model for the radiative flow of a non-Newtonian material under Joule heating: OHAM application
  8. Verification of the accommodative responses in viewing an on-axis analog reflection hologram
  9. Irreversibility as thermodynamic time
  10. A self-adaptive prescription dose optimization algorithm for radiotherapy
  11. Algebraic computational methods for solving three nonlinear vital models fractional in mathematical physics
  12. The diffusion mechanism of the application of intelligent manufacturing in SMEs model based on cellular automata
  13. Numerical analysis of free convection from a spinning cone with variable wall temperature and pressure work effect using MD-BSQLM
  14. Numerical simulation of hydrodynamic oscillation of side-by-side double-floating-system with a narrow gap in waves
  15. Closed-form solutions for the Schrödinger wave equation with non-solvable potentials: A perturbation approach
  16. Study of dynamic pressure on the packer for deep-water perforation
  17. Ultrafast dephasing in hydrogen-bonded pyridine–water mixtures
  18. Crystallization law of karst water in tunnel drainage system based on DBL theory
  19. Position-dependent finite symmetric mass harmonic like oscillator: Classical and quantum mechanical study
  20. Application of Fibonacci heap to fast marching method
  21. An analytical investigation of the mixed convective Casson fluid flow past a yawed cylinder with heat transfer analysis
  22. Considering the effect of optical attenuation on photon-enhanced thermionic emission converter of the practical structure
  23. Fractal calculation method of friction parameters: Surface morphology and load of galvanized sheet
  24. Charge identification of fragments with the emulsion spectrometer of the FOOT experiment
  25. Quantization of fractional harmonic oscillator using creation and annihilation operators
  26. Scaling law for velocity of domino toppling motion in curved paths
  27. Frequency synchronization detection method based on adaptive frequency standard tracking
  28. Application of common reflection surface (CRS) to velocity variation with azimuth (VVAz) inversion of the relatively narrow azimuth 3D seismic land data
  29. Study on the adaptability of binary flooding in a certain oil field
  30. CompVision: An open-source five-compartmental software for biokinetic simulations
  31. An electrically switchable wideband metamaterial absorber based on graphene at P band
  32. Effect of annealing temperature on the interface state density of n-ZnO nanorod/p-Si heterojunction diodes
  33. A facile fabrication of superhydrophobic and superoleophilic adsorption material 5A zeolite for oil–water separation with potential use in floating oil
  34. Shannon entropy for Feinberg–Horodecki equation and thermal properties of improved Wei potential model
  35. Hopf bifurcation analysis for liquid-filled Gyrostat chaotic system and design of a novel technique to control slosh in spacecrafts
  36. Optical properties of two-dimensional two-electron quantum dot in parabolic confinement
  37. Optical solitons via the collective variable method for the classical and perturbed Chen–Lee–Liu equations
  38. Stratified heat transfer of magneto-tangent hyperbolic bio-nanofluid flow with gyrotactic microorganisms: Keller-Box solution technique
  39. Analysis of the structure and properties of triangular composite light-screen targets
  40. Magnetic charged particles of optical spherical antiferromagnetic model with fractional system
  41. Study on acoustic radiation response characteristics of sound barriers
  42. The tribological properties of single-layer hybrid PTFE/Nomex fabric/phenolic resin composites underwater
  43. Research on maintenance spare parts requirement prediction based on LSTM recurrent neural network
  44. Quantum computing simulation of the hydrogen molecular ground-state energies with limited resources
  45. A DFT study on the molecular properties of synthetic ester under the electric field
  46. Construction of abundant novel analytical solutions of the space–time fractional nonlinear generalized equal width model via Riemann–Liouville derivative with application of mathematical methods
  47. Some common and dynamic properties of logarithmic Pareto distribution with applications
  48. Soliton structures in optical fiber communications with Kundu–Mukherjee–Naskar model
  49. Fractional modeling of COVID-19 epidemic model with harmonic mean type incidence rate
  50. Liquid metal-based metamaterial with high-temperature sensitivity: Design and computational study
  51. Biosynthesis and characterization of Saudi propolis-mediated silver nanoparticles and their biological properties
  52. New trigonometric B-spline approximation for numerical investigation of the regularized long-wave equation
  53. Modal characteristics of harmonic gear transmission flexspline based on orthogonal design method
  54. Revisiting the Reynolds-averaged Navier–Stokes equations
  55. Time-periodic pulse electroosmotic flow of Jeffreys fluids through a microannulus
  56. Exact wave solutions of the nonlinear Rosenau equation using an analytical method
  57. Computational examination of Jeffrey nanofluid through a stretchable surface employing Tiwari and Das model
  58. Numerical analysis of a single-mode microring resonator on a YAG-on-insulator
  59. Review Articles
  60. Double-layer coating using MHD flow of third-grade fluid with Hall current and heat source/sink
  61. Analysis of aeromagnetic filtering techniques in locating the primary target in sedimentary terrain: A review
  62. Rapid Communications
  63. Nonlinear fitting of multi-compartmental data using Hooke and Jeeves direct search method
  64. Effect of buried depth on thermal performance of a vertical U-tube underground heat exchanger
  65. Knocking characteristics of a high pressure direct injection natural gas engine operating in stratified combustion mode
  66. What dominates heat transfer performance of a double-pipe heat exchanger
  67. Special Issue on Future challenges of advanced computational modeling on nonlinear physical phenomena - Part II
  68. Lump, lump-one stripe, multiwave and breather solutions for the Hunter–Saxton equation
  69. New quantum integral inequalities for some new classes of generalized ψ-convex functions and their scope in physical systems
  70. Computational fluid dynamic simulations and heat transfer characteristic comparisons of various arc-baffled channels
  71. Gaussian radial basis functions method for linear and nonlinear convection–diffusion models in physical phenomena
  72. Investigation of interactional phenomena and multi wave solutions of the quantum hydrodynamic Zakharov–Kuznetsov model
  73. On the optical solutions to nonlinear Schrödinger equation with second-order spatiotemporal dispersion
  74. Analysis of couple stress fluid flow with variable viscosity using two homotopy-based methods
  75. Quantum estimates in two variable forms for Simpson-type inequalities considering generalized Ψ-convex functions with applications
  76. Series solution to fractional contact problem using Caputo’s derivative
  77. Solitary wave solutions of the ionic currents along microtubule dynamical equations via analytical mathematical method
  78. Thermo-viscoelastic orthotropic constraint cylindrical cavity with variable thermal properties heated by laser pulse via the MGT thermoelasticity model
  79. Theoretical and experimental clues to a flux of Doppler transformation energies during processes with energy conservation
  80. On solitons: Propagation of shallow water waves for the fifth-order KdV hierarchy integrable equation
  81. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part II
  82. Numerical study on heat transfer and flow characteristics of nanofluids in a circular tube with trapezoid ribs
  83. Experimental and numerical study of heat transfer and flow characteristics with different placement of the multi-deck display cabinet in supermarket
  84. Thermal-hydraulic performance prediction of two new heat exchangers using RBF based on different DOE
  85. Diesel engine waste heat recovery system comprehensive optimization based on system and heat exchanger simulation
  86. Load forecasting of refrigerated display cabinet based on CEEMD–IPSO–LSTM combined model
  87. Investigation on subcooled flow boiling heat transfer characteristics in ICE-like conditions
  88. Research on materials of solar selective absorption coating based on the first principle
  89. Experimental study on enhancement characteristics of steam/nitrogen condensation inside horizontal multi-start helical channels
  90. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part I
  91. Numerical exploration of thin film flow of MHD pseudo-plastic fluid in fractional space: Utilization of fractional calculus approach
  92. A Haar wavelet-based scheme for finding the control parameter in nonlinear inverse heat conduction equation
  93. Stable novel and accurate solitary wave solutions of an integrable equation: Qiao model
  94. Novel soliton solutions to the Atangana–Baleanu fractional system of equations for the ISALWs
  95. On the oscillation of nonlinear delay differential equations and their applications
  96. Abundant stable novel solutions of fractional-order epidemic model along with saturated treatment and disease transmission
  97. Fully Legendre spectral collocation technique for stochastic heat equations
  98. Special Issue on 5th International Conference on Mechanics, Mathematics and Applied Physics (2021)
  99. Residual service life of erbium-modified AM50 magnesium alloy under corrosion and stress environment
  100. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part I
  101. Diverse wave propagation in shallow water waves with the Kadomtsev–Petviashvili–Benjamin–Bona–Mahony and Benney–Luke integrable models
  102. Intensification of thermal stratification on dissipative chemically heating fluid with cross-diffusion and magnetic field over a wedge
Downloaded on 16.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2021-0078/html
Scroll to top button