Next Article in Journal
Pattern Recognition of Different Window Size Control Charts Based on Convolutional Neural Network and Information Fusion
Next Article in Special Issue
Hyperbolic Equations with Unknown Coefficients
Previous Article in Journal
A New Operational Matrices-Based Spectral Method for Multi-Order Fractional Problems
Previous Article in Special Issue
On Solutions of the Initial Value Problem for the Three-Term Fractional Differential Equation with Caputo Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact and Numerical Solitary Wave Structures to the Variant Boussinesq System

by
Abdulghani Alharbi
and
Mohammed B. Almatrafi
*
Department of Mathematics, College of Science, Taibah University, P.O. Box 344, Al-Madinah Al-Munawarah 30002, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 19 August 2020 / Revised: 2 September 2020 / Accepted: 3 September 2020 / Published: 8 September 2020
(This article belongs to the Special Issue Ordinary and Partial Differential Equations: Theory and Applications)

Abstract

:
Solutions such as symmetric, periodic, and solitary wave solutions play a significant role in the field of partial differential equations (PDEs), and they can be utilized to explain several phenomena in physics and engineering. Therefore, constructing such solutions is significantly essential. This article concentrates on employing the improved exp ( ϕ ( η ) ) -expansion approach and the method of lines on the variant Boussinesq system to establish its exact and numerical solutions. Novel solutions based on the solitary wave structures are obtained. We present a comprehensible comparison between the accomplished exact and numerical results to testify the accuracy of the used numerical technique. Some 3D and 2D diagrams are sketched for some solutions. We also investigate the L 2 error and the CPU time of the used numerical method. The used mathematical tools can be comfortably invoked to handle more nonlinear evolution equations.
MSC:
35A20; 35A99; 83C15; 65L12; 65N06; 65N40; 65N45; 65N50; 35Q51

1. Introduction

Most natural phenomena arising in various nonlinear sciences and engineering are modelled by nonlinear evolution equations. More specifically, a dramatic increase in the use of PDEs has been recently seen in some areas such as physics, biology, chemistry, economics, and computer sciences. Equations that describe shallow water waves appear in various fields of physics. The search for finding the exact solutions of such equations has been considerably given more attention in recent years. Although several approaches have been effectively developed, the exact traveling wave solutions for a massive number of NPDEs cannot be obtained. We point out some proposed techniques such as the truncated Painleve expansion process [1], the improved exp ( ϕ ( η ) ) -expansion technique [2], the projective Riccati equation technique [3], the Weierstrass elliptic function method [4], the extended tanh-procedure [5,6], the exp ( f ( ζ ) ) -expansion process [7,8,9], the sine–cosine approach [10,11], the Adomian decomposition technique [12,13], the Hirota’s bilinear technique [14,15], the inverse scattering transform [16], etc. The availability of some mathematical software such as Matlab, Maple, and Mathematica stimulates mathematicians and scientists to deal with insolvable nonlinear PDEs. Among the available numerical approaches, we state the following: the finite element method, the finite differences, the adaptive moving mesh technique [17], and the Parabolic Monge–Ampere method [18]. For more information about analytical and numerical solutions of NPDEs, one can refer to [19,20,21,22,23,24,25].
The variant Boussinesq system [26,27] reads
Φ t γ Φ x x t + Φ Φ x + Ψ x = 0 , Ψ t Φ Ψ x Φ x Ψ λ Φ x x x = 0 .
This system describes shallow water waves, where Φ indicates the velocity of wave, Ψ presents the height of free wave surface, and x and t play the role of the spatial and the temporal derivatives. The constants γ and λ > 0 are non-zero. The dynamics of shallow water waves are described by the Korteweg-de Vries (KdV) equation. However, the Boussinesq equation perfectly governs shallow water waves and gives much better approximation to such waves [28]. The exact solutions of system (1) have been derived by various scientists. For instance, the authors in [29] presented some traveling wave solutions for system (1) using the exp ( f ( ζ ) ) -expansion method. Huiqun [26] applied the extended Jacobi elliptic function expansion process to construct some periodic solutions for system (1). The improved ( G / G )-expansion technique was used in [27] to extract the exact traveling wave solutions of system (1). Zheng [30] applied the generalized Bernoulli sub-ODE method to extract traveling wave solutions for system (1). Patel et al. [31] utilized nonlinear Boussinesq equations to model the shallow water waves. Then, the Adams–Bashfourth (AB) predictor–corrector approach was used to derive the numerical solution of a nonlinear Boussinesq equation.
The principal purpose of this work is to seek the exact traveling wave solutions and the numerical solutions for system (1) using the improved exp ( ϕ ( η ) ) -expansion method and the method of lines. Since the improved exp ( ϕ ( η ) ) -expansion process depends on Jacobi elliptic functions, the trigonometric and hyperbolic solutions from the obtained solutions can be simply generated. The procedure is achieved by utilizing a wave transformation to convert system (1) into a system of ordinary differential equations (ODEs) solved by using the proposed method. Regarding the numerical solutions, the spatial derivatives are replaced by the finite difference formulae, whereas the temporal derivatives are left continuous. We graphically develop appropriate boundary conditions. The behavior of the exact solutions at the end points of the domain is invoked to construct that Φ x = Ψ x = 0 , and Φ x x x = Ψ x x x = 0 , as x ± . Thus, the boundary conditions are given by
Φ x = Φ x x x = 0 , and Ψ x = 0 at x ± .

2. Analysis of the Improved exp ( ϕ ( η ) ) -Expansion Approach

This section is assigned to summarize the improved exp ( ϕ ( η ) ) -expansion technique, as illustrated in [2]. Let
P ( Φ , Ψ , Φ x , Ψ x , Φ t , Ψ x x , Φ x x , Ψ x x x , Φ x x t , ) = 0 , Q ( Φ , Ψ , Φ x , Ψ x , Φ x x , Ψ t , Φ x x x , Ψ x x x , ) = 0 ,
be a given system of PDEs on the unknown functions Φ = Φ ( x , t ) , and Ψ = Ψ ( x , t ) . P and Q are polynomials in Φ , Ψ and their partial derivatives. In order to change system (3) into ODEs, we use the following transformations:
Φ ( x , t ) = ϕ ( η ) , Ψ ( x , t ) = ψ ( η ) , η = h x w t .
Inserting the transformations (4) into system (3) yields
p ( ϕ , ψ , ψ η , ϕ η , ψ η η , ϕ η η , ψ η η η , ϕ η η η , ) = 0 , q ( ϕ , ψ , ψ η , ϕ η , ψ η η , ϕ η η , ψ η η η , ϕ η η η , ) = 0 .
System (5) is integrated, if possible, term by term. For the sake of simplicity, we equate the integral constants to zero. As claimed by the proposed method, the solutions of system (5) are given by
ϕ ( η ) = j = 0 N α j exp ( j g ( η ) ) , ψ ( η ) = j = 0 M β j exp ( j g ( η ) ) ,
where the function g ( η ) fulfills the following ODE:
g η 2 = k 2 ( r + a 0 ) exp ( g ( η ) ) 2 + k 2 a 1 .
The constants w , r , a 0 , a 1 , k , h , α j and β j , j = 0 , 1 , , n , are evaluated later. N and M are positive integers which can be easily evaluated using the homogeneous balance. Equation (7) has various Jacobi elliptic function solutions presented in Appendix A. Plugging Equations (6) and (7) into system (5) and equating the coefficients of exp ( g ( η ) ) to zero lead to a system of algebraic equations which can be solved using any mathematical software. The solutions of this system determine the above-mentioned constants. Substituting these constants into Equation (6) gives the exact traveling wave solutions.

3. Exact Solutions of the Variant Boussinesq System

In this section, we attempt to introduce some new solitary wave solutions for the variant Boussinesq system [26,27], which is given by
Φ t γ Φ x x t + Φ Φ x + Ψ x = 0 , Ψ t Φ Ψ x Φ x Ψ λ Φ x x x = 0 ,
where γ and λ are arbitrary constants. Inserting Equation (4) into system (8) leads to
w ϕ η + w h 2 γ ϕ η η η + h ϕ ϕ η + h ψ η = 0 , w ψ η h ϕ ψ η h ϕ η ψ h 3 λ ϕ η η η = 0 .
We now integrate each equation in system (9) once with respect to η . Achieving this, we have
w ϕ + w h 2 γ ϕ η η + h 2 ϕ 2 + h ψ = 0 , w ψ h ϕ ψ h 3 λ ϕ η η = 0 .
In the first equation in system (10), we balance the highest order ϕ η η with the nonlinear term ϕ 2 while, in the second equation, we consider the homogeneous balance between ϕ η η and ϕ ψ . Consequently, we have N = 2 and M = 2 . Thus, the solutions are expressed by
ϕ ( η ) = α 0 + α 1 exp ( g ( η ) ) + α 2 exp ( 2 g ( η ) ) , ψ ( η ) = β 0 + β 1 exp ( g ( η ) ) + β 2 exp ( 2 g ( η ) ) .
The values of the constants α j and β j are evaluated later. Plugging system (11) into system (10) and equating the coefficients of exp ( n g ( η ) ) , n = 0 , 1 , 2 , 3 , 4 , to zero lead to some algebraic equations whose solutions are shown in various cases as follows:
  • Case 1
    α 0 = λ a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 3 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 3 λ r 2 γ a 1 2 3 a 0 r , β 0 = λ γ 2 k 4 a 1 2 3 a 0 r + a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 3 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = k 2 λ γ a 1 2 3 a 0 r 2 6 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = 1 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 2
    α 0 = λ a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 3 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 3 λ r 2 γ a 1 2 3 a 0 r , β 0 = λ γ 2 k 4 a 1 2 3 a 0 r a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 3 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = ± i k 2 λ γ a 1 2 3 a 0 r 2 6 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = i 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 3
    α 0 = λ a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 3 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 3 λ r 2 γ a 1 2 3 a 0 r , β 0 = λ γ 2 k 4 a 1 2 3 a 0 r a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 3 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = i k 2 λ γ a 1 2 3 a 0 r 2 6 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = i 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 4
    α 0 = λ a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 3 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 3 λ r 2 γ a 1 2 3 a 0 r , β 0 = λ γ 2 k 4 a 1 2 3 a 0 r + a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 3 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = ± k 2 λ γ a 1 2 3 a 0 r 2 6 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = 1 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 5
    α 0 = λ 3 a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 9 λ r 2 γ a 1 2 3 a 0 r , β 0 = 3 λ γ 2 k 4 a 1 2 3 a 0 r + a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 9 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = 3 2 k 2 λ γ a 1 2 3 a 0 r 2 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = 3 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 6
    α 0 = λ 3 a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 9 λ r 2 γ a 1 2 3 a 0 r , β 0 = 3 λ γ 2 k 4 a 1 2 3 a 0 r a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 9 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = ± i 3 2 k 2 λ γ a 1 2 3 a 0 r 2 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = i 3 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 7
    α 0 = λ 3 a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 9 λ r 2 γ a 1 2 3 a 0 r , β 0 = 3 λ γ 2 k 4 a 1 2 3 a 0 r a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 9 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = i 3 2 k 2 λ γ a 1 2 3 a 0 r 2 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = i 3 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
  • Case 8
    α 0 = λ 3 a 1 γ 2 k 4 a 1 2 3 a 0 r 3 a 0 γ k 2 r + a 1 2 γ k 2 2 γ a 1 2 3 a 0 r γ 2 k 4 a 1 2 3 a 0 r , α 1 = 0 , α 2 = 9 λ r 2 γ a 1 2 3 a 0 r , β 0 = 3 λ γ 2 k 4 a 1 2 3 a 0 r + a 1 γ k 2 4 γ γ 2 k 4 a 1 2 3 a 0 r , β 1 = 0 , β 2 = 9 k 2 λ r 4 γ 2 k 4 a 1 2 3 a 0 r , w = ± 3 2 k 2 λ γ a 1 2 3 a 0 r 2 γ 2 k 4 a 1 2 3 a 0 r 3 / 4 , h = 3 2 2 γ 2 k 4 a 1 2 3 a 0 r 4 .
Sequentially, we construct the following several cases for the traveling wave solutions of the approached problem when m = 1 , (see Appendix A).
ϕ 1 ( η ) = ± λ γ 2 k 4 2 γ k 2 2 3 γ 3 / 2 k 2 3 λ 2 γ tanh 2 ( k η ) , ψ 1 ( η ) = 3 λ γ 2 k 4 2 γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 tanh 2 ( k η ) .
ϕ 2 ( η ) = ± λ γ 2 k 4 2 γ k 2 2 3 γ 3 / 2 k 2 3 λ 2 γ coth 2 ( k η ) , ψ 2 ( η ) = 3 λ γ 2 k 4 2 γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 coth 2 ( k η ) .
ϕ 3 ( η ) = ± λ γ k 2 3 γ 2 k 4 2 γ γ 2 k 4 ± 9 λ 2 γ cosh 2 ( k η ) , ψ 3 ( η ) = 3 λ γ 2 k 4 + γ k 2 4 γ γ 2 k 4 + 9 k 2 λ 4 γ 2 k 4 cosh 2 ( k η ) .
ϕ 4 ( η ) = ± λ γ k 2 3 γ 2 k 4 2 γ γ 2 k 4 9 λ 2 γ sinh 2 ( k η ) , ψ 4 ( η ) = 3 λ γ 2 k 4 + γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 sinh 2 ( k η ) .
Thus, the exact solutions of system (8) are
ϕ 1 ( x , t ) = ± λ γ 2 k 4 2 γ k 2 2 3 γ 3 / 2 k 2 3 λ 2 γ tanh 2 ( k ( 3 2 2 γ 2 k 4 4 x ± 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) , ψ 1 ( x , t ) = 3 λ γ 2 k 4 2 γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 tanh 2 ( k ( 3 2 2 γ 2 k 4 4 x ± 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) .
ϕ 2 ( x , t ) = ± λ γ 2 k 4 2 γ k 2 2 3 γ 3 / 2 k 2 3 λ 2 γ coth 2 ( k ( 3 2 2 γ 2 k 4 4 x ± 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) , ψ 2 ( x , t ) = 3 λ γ 2 k 4 2 γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 coth 2 ( k ( 3 2 2 γ 2 k 4 4 x ± 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) .
ϕ 3 ( x , t ) = ± λ γ k 2 3 γ 2 k 4 2 γ γ 2 k 4 ± 9 λ 2 γ cosh 2 ( k ( 3 2 2 γ 2 k 4 4 x 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) , ψ 3 ( x , t ) = 3 λ γ 2 k 4 + γ k 2 4 γ γ 2 k 4 + 9 k 2 λ 4 γ 2 k 4 cosh 2 ( k ( 3 2 2 γ 2 k 4 4 x 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) .
ϕ 4 ( x , t ) = ± λ γ k 2 3 γ 2 k 4 2 γ γ 2 k 4 9 λ 2 γ sinh 2 ( k ( 3 2 2 γ 2 k 4 4 x 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) , ψ 4 ( x , t ) = 3 λ γ 2 k 4 + γ k 2 4 γ γ 2 k 4 9 k 2 λ 4 γ 2 k 4 sinh 2 ( k ( 3 2 2 γ 2 k 4 4 x 3 2 γ k 2 λ 2 γ 2 k 4 3 / 4 t ) ) .

4. Numerical Solutions of the Variant Boussinesq System

The numerical solutions of system (8) are discussed in this section using the method of lines. The domain on which we work is given by [ 0 , L x ] . We start by writing the variable U on the form:
U = Φ γ Φ x x .
Thus, system (8) is reformed as
U t + Φ Φ x + Ψ x = 0 , Ψ t Φ Ψ x λ Φ x x x = 0 .
The related boundary conditions are presented by Equation (2). For much better numerical results for system (29), we apply a uniform mesh technique on the domain [ 0 , L x ] . Then, we divide the domain into N x sub-intervals [ x n , x n + 1 ] with fixed step size h x = L x / N x such that
x n = ( n 1 ) h x , n = 1 , 2 , , N x + 1 ,
where h x plays the role of a uniform width of each sub-interval. Note that the spatial derivatives in system (29) are replaced with finite differences, whereas the temporal differentiation is left continuous. As a result, system (29) is discretized as
U t | n = Φ n Φ n + 1 / 2 Φ n 1 / 2 h x Ψ n + 1 / 2 Ψ n 1 / 2 h x , Ψ t | n = ( Φ Ψ ) n + 1 / 2 ( Φ Ψ ) n 1 / 2 h x + β Φ x x x | n , U n = Φ n γ Φ x x | n ,
where n = 2 , , N x . The space discretization of Φ x x x | n , ( Φ Ψ ) n + 1 / 2 , Φ x x | n and Φ n + 1 / 2 , Ψ n + 1 / 2 is presented in Appendix B. The relevant boundary conditions for Equation (2) are provided by
U t , 1 = U t , N x + 1 = 0 , Ψ t , 1 = Ψ t , N x + 1 = 0 .
The initial condition is derived from Equation (24) by taking t = 0 , as shown in Figure 1. In Figure 2, we present the time evolution of the exact and numerical solutions of Φ ( x , t ) for 0 t 24 . Moreover, Figure 3 illustrates the time evolution of the exact and numerical solution of Ψ ( x , t ) for 0 t 24 . In order to evaluate Φ x x x , Ψ x and Φ x at x = 0 and x = L x , we use some fictitious points given by
Φ 0 = Φ 2 , Φ N x + 2 = Φ N x , t [ 0 , T f ] , Ψ 0 = Ψ 2 , Φ x x | 0 = Φ x x | 2 , Ψ N x + 2 = Ψ N x , Φ x x | N x + 2 = Φ x x | N x .
It is notable mentioning that MATLAB ODE solver (ode15i) is employed to solve the numerical system. This solver is a variable order implicit time-stepping approach depended on the numerical differentiation formulas.

5. Results and Discussion

The improved exp ( ϕ ( η ) ) -expansion method is greatly applied on system (1) to derive several exact traveling wave solutions. Although this method is based on the Jacobi elliptic functions, its solutions can be converted into trigonometric and hyperbolic functions. We effectively extract many solutions expressed on the form of trigonometric and hyperbolic functions. The validity of the solutions is verified by substituting the obtained solutions into the leading equations. The presented exact solutions are more general than those obtained in [29]. Zheng [30] employed the generalized Bernoulli sub-ODE approach on system (1) and introduced two rational solutions. On the contrary, by using the improved exp ( ϕ ( η ) ) -expansion approach in this article, we obtain several solutions for system (1).
The execution of the method of lines leads to efficient and adequate results. For example, an appropriate coincidence between the numerical solutions is depicted in Figure 4 and Figure 5. The solutions almost have the same behavior. Moreover, Figure 6 illustrates the accuracy of the method of lines employed in this paper. As can be seen from Figure 6, we have used a small value for h x = 0.1 . However, the error is high. This error has been successfully reduced by taking smaller values of h x . When we consider h x = 0.01 , the numerical solutions (green sold lines) nicely converge to exact solutions. For h x = 0.0035 , the numerical results nearly meet the exact results. Note that the above presented figures are sketched under the values λ = 0.5 , γ = 0.7 , a 1 = 2 , a 0 = 1 , r = 1 , k = 1 , x 0 = 20 , x = 0 35 and t = 0 24 . The error is digitally shown in Table 1. L 2 has rapidly decreased for small values of h x . The L 2 error in the numerical solutions of Φ ( x , t ) and Ψ ( x , t ) has reached 9.30 × 10 3 and 2.30 × 10 3 , respectively, during 3.43 × 10 2 seconds when h x = 1 × 10 1 . Nevertheless, the method works well enough when h x = 3.5 × 10 3 . Here, the error stands at 4.42 × 10 7 and 2.62 × 10 6 for Φ ( x , t ) and Ψ ( x , t ) , respectively. The CPU time has slowly increased to hit 1.29 × 10 + 1 s.

6. Conclusions

This article has focused on developing the exact and numerical solutions of system (1) by taking advantage of the improved exp ( ϕ ( η ) ) -expansion approach and the method of lines, respectively. The improved exp ( ϕ ( η ) ) -expansion method depends on Jacobi elliptic functions which have been used to degenerated trigonometric functions. Numerous exact solutions have been well introduced. The obtained numerical solutions roughly match the exact solutions for a small value of h x . In other words, the curves of the numerical solutions differ substantially for huge h x and quickly converge together for small h x . The method of lines performs adequately well when we take the step size smaller. The utilized techniques are practical and effective to be employed on more sophisticated nonlinear PDEs.

Author Contributions

Both authors made an equal contribution to prepare the manuscript. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Jacobi Elliptic Function Solutions

Here, we mention some significant solutions for Equation (7).
  • If a 0 = 1 , a 1 = ( 1 + m 2 ) , and r = m 2 . Then, g 1 = l n ( s n ( k η , m ) ) .
  • If a 0 = 1 m 2 , a 1 = 1 + 2 m 2 , and r = m 2 . Then, g 2 = l n ( n c ( k η , m ) ) .
  • If a 0 = m 2 , a 1 = ( 1 + m 2 ) , and r = 1 . Then, g 3 = l n ( n s ( k η , m ) ) .
  • If a 0 = m 2 , a 1 = 2 m 2 1 , and r = 1 m 2 . Then, g 4 = l n ( c n ( k η , m ) ) .
  • If a 0 = 1 , a 1 = 2 m 2 , and r = 1 m 2 . Then, g 5 = l n ( c s ( k η , m ) ) .
  • If a 0 = 1 m 2 , a 1 = 2 m 2 , and r = 1 . Then, g 6 = l n ( s c ( k η , m ) ) .
Here, m indicates the modulus of the Jacobi elliptic function.

The Modulus of the Jacobi Elliptic Function

It should be noted that, if the modulus m 1 , then s n ( ξ , m ) t a n h ( ξ ) , c n ( ξ , m ) s e c h ( ξ ) , n s c o t h ( ξ ) , c s ( ξ , m ) c s c h ( ξ ) . However, if the modulus m 0 , then n s ( ξ , m ) c s c ( ξ ) , n c ( ξ , m ) s e c ( ξ ) , s c ( ξ , m ) t a n ( ξ ) , c s ( ξ , m ) c o t ( ξ ) .

Appendix B. Space Discretization

The high-order terms of Equation (30) are discretized as follows:
Φ n + 1 / 2 = Φ n + 1 + Φ n 2 , Ψ n + 1 / 2 = Ψ n + 1 + Ψ n 2 , ( Φ Ψ ) n + 1 / 2 = Φ n + 1 Ψ n + 1 + Φ n Ψ n 2 , Φ x x x | n = 1 2 h x 3 Φ n + 2 2 Φ n + 1 + 2 Φ n 1 Φ n 2 , Φ x x | n = 1 h x 2 Φ n + 1 2 Φ n + Φ n 1 .

References

  1. Wang, M.L.; Li, X.Z. Extended F-expansion and periodic wave solutions for the generalized Zakharov equations. Phys. Lett. A 2005, 343, 48–54. [Google Scholar] [CrossRef]
  2. Chen, G.; Xin, X.; Liu, H. The improved exp(-ϕ(η))-expansion method and new exact solutions of nonlinear evolution equations in mathematical physics. Adv. Math. Phys. 2019, 2019. [Google Scholar] [CrossRef] [Green Version]
  3. Conte, R.; Musette, M. Link between solitary waves and projective Riccati equations. J. Phys. A Math. 1992, 25, 5609–5623. [Google Scholar] [CrossRef]
  4. Chow, K.W. A class of exact periodic solutions of nonlinear envelope equation. J. Math. Phys. 1995, 36, 4125–4137. [Google Scholar] [CrossRef] [Green Version]
  5. Fan, E. Extended tanh-function method and its applications to nonlinear equations. Phys. Lett. A 2000, 277, 212–218. [Google Scholar] [CrossRef]
  6. Wazwaz, A.M. The extended tanh method for abundant solitary wave solutions of nonlinear wave equations. Appl. Math. Comput. 2007, 187, 1131–1142. [Google Scholar] [CrossRef]
  7. Alam, M.N.; Tunc, C. An analytical method for solving exact solutions of the nonlinear Bogoyavlenskii equation and the nonlinear diffusive predator-prey system. Alex Eng. J. 2016, 55, 1855–1865. [Google Scholar] [CrossRef] [Green Version]
  8. Khan, K.; Akbar, M.A. Application of exp(-φ(ζ))-expansion method to find the exact solutions of modified Benjamin-Bona-Mahony equation. World Appl. Sci. J. 2013, 24, 1373–1377. [Google Scholar]
  9. Khan, K.; Akbar, M.A. Exact traveling wave solutions of Kadomtsev-Petviashvili equation. J. Egypt Math. Soc. 2015, 23, 278–281. [Google Scholar] [CrossRef] [Green Version]
  10. Wazwaz, A.M. Exact solutions to the double sinh-Gordon equation by the tanh method and a variable separated ODE. method. Comput. Math. Appl. 2005, 50, 1685–1696. [Google Scholar] [CrossRef] [Green Version]
  11. Wazwaz, A.M. A sine–cosine method for handling nonlinear wave equations. Math. Comput. Model. 2004, 40, 499–508. [Google Scholar] [CrossRef]
  12. Adomain, G. Solving Frontier Problems of Physics: The Decomposition Method; Kluwer Academic Publishers: Boston, MA, USA, 1994. [Google Scholar]
  13. Wazwaz, A.M. Partial Differential Equations: Method and Applications; Taylor and Francis: Abingdon, UK, 2002. [Google Scholar]
  14. Hirota, R. Exact envelope soliton solutions of a nonlinear wave equation. J. Math. Phys. 1973, 14, 805–810. [Google Scholar] [CrossRef]
  15. Hirota, R.; Satsuma, J. Soliton solution of a coupled KdV equation. Phys. Lett. A 1981, 85, 407–408. [Google Scholar] [CrossRef]
  16. Ablowitz, M.J.; Clarkson, P.A. Solitons, Non-Linear Evolution Equations and Inverse Scattering Transform; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  17. Huang, W.; Russell, R.D. The Adaptive Moving Mesh Methods; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  18. Budd, C.J.; Huang, W. Russell, R.D. Adaptivity with moving grids. Acta Numer. 2009, 18, 111–241. [Google Scholar] [CrossRef]
  19. Alharbi, A.R.; Almatrafi, M.B.; Abdelrahman, M.A.E. The Extended Jacobian Elliptic Function Expansion Approach to the Generalized Fifth Order KdV Equation. J. Phy. Math. 2019, 10, 310. [Google Scholar]
  20. Alharbi, A.R.; Almatrafi, M.B. Numerical investigation of the Dispersive Long Wave Equation using an adaptive moving mesh method and its stability. Results Phys. 2020, 16, 102870. [Google Scholar] [CrossRef]
  21. Alharbi, A.R.; Almatrafi, M.B. Riccati-Bernoulli Sub-ODE approach on the partial differential equations and applications. Int. J. Math. Comput. Sci. 2020, 15, 367–388. [Google Scholar]
  22. Abdelrahman, M.A.E.; Almatrafi, M.B.; Alharbi, A. Fundamental solutions for the coupled KdV system and its stability. Symmetry 2020, 12, 429. [Google Scholar] [CrossRef] [Green Version]
  23. Alharbi, A.; Abdelrahman, M.A.E.; Almatrafi, M.B. Analytical and numerical investigation for the DMBBM equation. Comput. Model. Eng. Sci. 2020, 122, 743–756. [Google Scholar] [CrossRef]
  24. Alam, M.N.; Tunç, C. Constructions of the optical solitons and other solitons to the conformable fractional Zakharov-Kuznetsov equation with power law nonlinearity. J. Taibah Univ. Sci. 2019, 14, 94–100. [Google Scholar] [CrossRef] [Green Version]
  25. Shahida, N. Tunç, C. Resolution of coincident factors in altering the flow dynamics of an MHD elastoviscous fluid past an unbounded upright channel. J. Taibah Univ. Sci. 2019, 13, 1022–1034. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, H.Q. Extended Jacobi elliptic function expansion method and its applications. Commun. Nonlinear Sci. Numer. Simul. 2007, 12, 627–635. [Google Scholar] [CrossRef]
  27. Zayed, E.M.E.; Al-Joudi, S. An Improved (G/G)-expansion Method for Solving Nonlinear PDEs in Mathematical Physics. ICNAAM AIP Conf. Proc. 2010, 1281, 2220–2224. [Google Scholar]
  28. Ayub, K.; Saeed, M.; Ashraf, M.; Yaqub, M.; Hassan, M. Soliton solutions of Variant Boussinesq equations through Exp-function method. UW J. Sci. Technol. 2017, 1, 24–30. [Google Scholar]
  29. Alharbi, A.R.; Almatrafi, M.B. Analytical and numerical solutions for the variant Boussinseq equations. J. Taibah Univ. Sci. 2020, 14, 454–462. [Google Scholar] [CrossRef]
  30. Zheng, B. Application Of A generalized Bernoulli Sub-ODE method for finding traveling solutions of some nonlinear equations. Wseas Trans. Math. 2012, 11, 618–626. [Google Scholar]
  31. Patel, P.; Kumar, P.; Rajni. The Numerical Solution of Boussinesq Equation for Shallow Water Waves. AIP Conf. Proc. 2020, 2214, 020019. [Google Scholar] [CrossRef]
Figure 1. Figure (a) illustrates the behavior of the solution Φ ( x , 0 ) (Equation (24)) while Figure (b) shows the behavior of the solution Ψ ( x , 0 ) (Equation (24)). The used parameters are λ = 0.5 , γ = 0.7 , a 1 = 2 , a 0 = 1 , r = 1 , k = 1 , x 0 = 20 , x = 0 35 a n d t = 0 24 .
Figure 1. Figure (a) illustrates the behavior of the solution Φ ( x , 0 ) (Equation (24)) while Figure (b) shows the behavior of the solution Ψ ( x , 0 ) (Equation (24)). The used parameters are λ = 0.5 , γ = 0.7 , a 1 = 2 , a 0 = 1 , r = 1 , k = 1 , x 0 = 20 , x = 0 35 a n d t = 0 24 .
Symmetry 12 01473 g001
Figure 2. Figure (a) shows the time evolution of the exact solution of Φ ( x , t ) while Figure (b) presents the time evolution of the numerical solution of Φ ( x , t ) . We consider 0 t 24 .
Figure 2. Figure (a) shows the time evolution of the exact solution of Φ ( x , t ) while Figure (b) presents the time evolution of the numerical solution of Φ ( x , t ) . We consider 0 t 24 .
Symmetry 12 01473 g002
Figure 3. Figure (a) shows the time evolution of the exact solution of Ψ ( x , t ) while Figure (b) presents the time evolution of the numerical solution of Ψ ( x , t ) . We consider 0 t 24 .
Figure 3. Figure (a) shows the time evolution of the exact solution of Ψ ( x , t ) while Figure (b) presents the time evolution of the numerical solution of Ψ ( x , t ) . We consider 0 t 24 .
Symmetry 12 01473 g003
Figure 4. 3D figures comparing the performance of the numerical method with the exact solution. The exact solution of Φ ( x , t ) (Figure (a)) and the numerical solution of Φ ( x , t ) (Figure (b)) are graphically compared in these plots.
Figure 4. 3D figures comparing the performance of the numerical method with the exact solution. The exact solution of Φ ( x , t ) (Figure (a)) and the numerical solution of Φ ( x , t ) (Figure (b)) are graphically compared in these plots.
Symmetry 12 01473 g004
Figure 5. Figure (a) presents a 3D surface for the exact solution of Ψ ( x , t ) while Figure (b) illustrates a 3D surface for the numerical solution of Ψ ( x , t ) . The numerical graph seems to be identical with the exact one.
Figure 5. Figure (a) presents a 3D surface for the exact solution of Ψ ( x , t ) while Figure (b) illustrates a 3D surface for the numerical solution of Ψ ( x , t ) . The numerical graph seems to be identical with the exact one.
Symmetry 12 01473 g005
Figure 6. Figure (a) compares some numerical solutions of Φ ( x , t ) with the exact solution of Φ ( x , t ) for various values of h x . In Figure (b), we present a comparison between the exact and numerical solutions of Ψ ( x , t ) . The numerical solutions of Φ ( x , t ) and Ψ ( x , t ) approach the exact solutions if h x is very small, as can be observed in these figures.
Figure 6. Figure (a) compares some numerical solutions of Φ ( x , t ) with the exact solution of Φ ( x , t ) for various values of h x . In Figure (b), we present a comparison between the exact and numerical solutions of Ψ ( x , t ) . The numerical solutions of Φ ( x , t ) and Ψ ( x , t ) approach the exact solutions if h x is very small, as can be observed in these figures.
Symmetry 12 01473 g006
Table 1. L 2 error and CPU time consumed to arrive at t = 10 for the numerical process.
Table 1. L 2 error and CPU time consumed to arrive at t = 10 for the numerical process.
h x L 2 Error for Φ L 2 Error for Ψ CPU
1 × 10 1 9.30 × 10 3 2.30 × 10 3 3.43 × 10 2 s
5 × 10 2 5.80 × 10 4 1.63 × 10 4 6.32 × 10 2 s
2 × 10 2 1.53 × 10 5 7.04 × 10 6 2.92 × 10 1 s
1 × 10 2 1.36 × 10 6 2.80 × 10 6 8.30 × 10 1 s
5 × 10 3 4.86 × 10 7 2.66 × 10 6 0.33 × 10 + 1 s
3.5 × 10 3 4.42 × 10 7 2.62 × 10 6 1.29 × 10 + 1 s

Share and Cite

MDPI and ACS Style

Alharbi, A.; Almatrafi, M.B. Exact and Numerical Solitary Wave Structures to the Variant Boussinesq System. Symmetry 2020, 12, 1473. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12091473

AMA Style

Alharbi A, Almatrafi MB. Exact and Numerical Solitary Wave Structures to the Variant Boussinesq System. Symmetry. 2020; 12(9):1473. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12091473

Chicago/Turabian Style

Alharbi, Abdulghani, and Mohammed B. Almatrafi. 2020. "Exact and Numerical Solitary Wave Structures to the Variant Boussinesq System" Symmetry 12, no. 9: 1473. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12091473

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop