Next Article in Journal
Real Space Theory for Electron and Phonon Transport in Aperiodic Lattices via Renormalization
Previous Article in Journal
BRST Cohomologies of Mixed and Second Class Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fundamental Solutions for the Coupled KdV System and Its Stability

by
Mahmoud A. E. Abdelrahman
1,2,*,
M. B. Almatrafi
1 and
Abdulghani Alharbi
1
1
Department of Mathematics, College of Science, Taibah University, Al-Madinah Al-Munawarah 344, Saudi Arabia
2
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 2 February 2020 / Revised: 20 February 2020 / Accepted: 23 February 2020 / Published: 7 March 2020

Abstract

:
In this paper, we establish exact solutions for the non-linear coupled KdV equations. The exp-function method is used to construct the solitary travelling wave solutions for these equations. The numerical adaptive moving mesh PDEs (MMPDEs) method is also implemented in order to solve the proposed coupled KdV equations. The achieved results may be applicable to some plasma environments, such as ionosphere plasma. Some numerical simulations compared with the exact solutions are provided to illustrate the validity of the proposed methods. Furthermore, the modulational instability is analyzed based on the standard linear-stability analysis. The depiction of the techniques are straight, powerful, robust and can be applied to other nonlinear systems of partial differential equations.
MSC:
35A20, 35A99, 83C15, 65L12, 65N06, 65N40, 65N45, 65N50, 35Q51

1. Introduction

The nonlinear partial differential equations (NPDEs) are widely used to express a variety of physical circumstances in fields such as biology, economy, engineering, elastic media, meteorology, plasma physics, optics, fluid mechanics and chemical physics, see [1,2,3,4,5,6,7,8,9,10,11,12]. In the last three decades, most researchers have proposed and developed some analytical and numerical methods to find exact and numerical solutions for NPDEs, such as F-expansion method [13], ( G G ) expansion method [14], homogeneous balance method [15], Jacobi elliptic function method [16], exp-function method [17], sine-cosine method [18], tanh-sech method [19], Lie symmetry analysis [20], RB sub-ODE method [21,22], adaptive moving mesh method [23,24]. Indeed, there are recent developments of symmetry methods to obtain special reduced solutions of NPDEs [25].
The study of non-linear complex waves is very important in various applied sciences, especially in optics, biology, fluid mechanics, engineering, physics, higher-order symmetries, and solid state physics and chemical physics. The Kortweg-de Vries (KdV) equations have been part of an important class of non-linear evolution equations with numerous applications in physics, plasma and engineering fields. All these applications start from a more or less general physical model and end up in the KdV equation by considering a specific limit of the physical problem. In the theory of rogue waves, the nonlinear Schrödinger equation plays a dominant role for waves in deep water, while the KdV equation describes the effects in shallow water. In plasma physics, the KdV equations produce ion-acoustic solutions [26]. Several well-know researchers obtain symmetries because they can be used to obtain systematically exact solutions of the KdV equations [27,28].
This article is concerned with the following coupled KdV equations [29]:
V t α V x x x 6 α V V x β H H x = 0 , H t + H x x x + 3 V H x = 0 ,
subject to
V ( L , t ) = f l ( L , t ) , H ( L , t ) = g l ( L , t ) , V ( L , t ) = f R ( L , t ) , H ( L , t ) = g R ( L , t ) , t 0 ,
where α and β are nonzero parameters, x [ L , L ] , f l , f R , g l and g R are functions of t . System (1) was derived by Hirota and Satsuma [29] to describe iterations of water waves with different dispersion relations. Recently, many authors mainly had paid attention to study solutions of coupled equations by using various methods. Among these are the trigonometric function transform method [30], F-expansion method [13], the homotopy perturbation [31] and differential transformation Method [32]. For further methods used to solve system (1), see [33,34] and references therein.
The novelties and the significance of the results in this article are mainly exhibited in three aspects. First, we achieve some new solutions for the coupled KdV Equation (1), utilizing the exp-function method [17,35]. Second, we implement the adaptive moving mesh PDEs (MMPDEs) scheme to solve the coupled KdV Equation (1). We compare these numerical solutions with the exact solutions obtained by the exp-function method. To the best of our knowledge, the proposed methods have not been used to solve Equation (1). Third, the modulation instability is employed in order to study the stability of exact solutions.

2. Summary of the exp ( Γ ( ξ ) ) -Expansion Method

We give a brief explanation for the exp ( Γ ( ξ ) ) -expansion approach [17,35]. This technique is described in the following steps: Assume that
T ( φ , φ t , φ x , φ t t , φ x x , φ x t , ) = 0 ,
where φ ( x , t ) is an unknown function and T is a polynomial of φ is a PDE with two independent variables x and t. In order to reduce Equation (3) to an ordinary differential equation (ODE), we introduce the transformation
ξ = x ± w t , φ ( x , t ) = φ 0 ( ξ ) ,
where w is a constant. Expanding the derivatives in Equation (3) using the chain rule leads to the ODE
R ( φ 0 , φ 0 , ξ , φ 0 , ξ ξ , φ 0 , ξ ξ ξ , ) = 0 ,
where R is a polynomial in φ 0 ( ξ ) and
φ 0 , ξ = d φ 0 d ξ , φ 0 , ξ ξ = d 2 φ 0 d ξ 2 , φ 0 , ξ ξ ξ = d 3 φ 0 d ξ 3 ,
According to this technique, the solution of Equation (5) is given in the form:
φ 0 ( ξ ) = k = 0 n b k e k Γ ( ξ ) ,
where the constants b 0 , b 1 , b 2 , , b n will be estimated later. Furthermore, b n 0 , and Γ = Γ ( ξ ) satisfies the equation
Γ ξ = e Γ ( ξ ) + μ e Γ ( ξ ) + λ .
We now turn to present various cases for the solutions of Equation (8).
  • If λ 2 4 μ > 0 , μ 0 , then
    • Γ ( ξ ) = l n ( λ 2 4 μ ) tanh ( λ 2 4 μ ) 2 ( ξ + γ ) λ 2 μ ,
    • Γ ( ξ ) = l n ( λ 2 4 μ ) coth ( λ 2 4 μ ) 2 ( ξ + γ ) λ 2 μ .
  • If λ 2 4 μ < 0 , μ 0 , then
    • Γ ( ξ ) = l n ( 4 μ λ 2 ) tan ( 4 μ λ 2 ) 2 ( ξ + γ ) λ 2 μ ,
    • Γ ( ξ ) = l n ( 4 μ λ 2 ) cot ( 4 μ λ 2 ) 2 ( ξ + γ ) λ 2 μ .
  • If λ 2 4 μ > 0 , μ = 0 , λ 0 , then
    Γ ( ξ ) = l n λ exp ( λ ( ξ + γ ) ) 1 .
  • If λ 2 4 μ = 0 , μ 0 , λ 0 , then
    Γ ( ξ ) = l n 2 ( λ ( ξ + k ) + 2 λ 2 ( ξ + γ ) .
  • If λ 2 4 μ = 0 , μ = 0 , λ = 0 , then
    Γ ( ξ ) = l n ξ + γ ,
    where γ , μ , λ and w are arbitrary constants. When we insert Equation (7) into Equation (5) and equate the coefficients of the same order of e Γ ( ξ ) , we end up with an algebraic system that can be simply solved by using Mathematica to evaluate the values of b 0 , b 1 , b 2 , , b n .

Stability Analysis

We examine the stability of the achieved traveling wave solutions by using the Hamiltonian system form [36,37], which is
ρ 1 ( w ) = 1 2 φ 1 2 ( ξ ) d x , ρ 2 ( w ) = 1 2 1 2 φ 2 2 ( ξ ) d x ,
where w is the wave speed, ρ i ( w ) indicates the momentum and φ 1 and φ 2 are the obtained traveling wave solutions of the system (1). A sufficient condition for stability is given by
ρ i w > 0 , i = 1 , 2 .
Hence, if conditions (12) and (13) are met, then it is stable on unbounded domains. If conditions (12) and (13) are not satisfied, then it is unstable. The criteria for the stability of traveling wave solutions are discussed here for the system (1) at particular intervals. For more details about the stability of the exact solutions, we refer to [37,38,39] and references therein.

3. The Exact Solution of the Coupled KdV Equations

Consider that the traveling wave solutions of Equation (1) are given by
ξ = x + w t V ( x , t ) = v ( ξ ) and H ( x , t ) = h ( ξ ) .
to alter system (1) into the ODEs
w v ξ α v ξ ξ ξ 3 α ( v 2 ) ξ β ( h 2 ) ξ = 0 , w h ξ + h ξ ξ ξ + 3 v h ξ = 0 ,
Integrating Equation (15) concerning ξ once and neglecting the constant of integration, we have
w v α v ξ ξ 3 α v 2 β h 2 = 0 , w h ξ + h ξ ξ ξ + 3 v h ξ = 0 ,
Balancing the order of v ξ ξ with the exponent of v 2 from the first equation of the system (16), gives N = 2 and the order of h ξ ξ ξ with the exponent of v h ξ from the second equation of the system (16), gives M = 2 . Hence, the solutions take the form of
v ( ξ ) = r = 0 N a r e r Γ ( ξ ) , h ( ξ ) = r = 0 M b r e r Γ ( ξ ) .
And then
v ( ξ ) = a 0 + a 1 e Γ ( ξ ) + a 2 e 2 Γ ( ξ ) , h ( ξ ) = b 0 + b 1 e Γ ( ξ ) + b 2 e 2 Γ ( ξ ) ,
where the values of the constants a 0 , a 1 , a 2 , b 0 , b 1 and b 2 are computed later. Substituting Equation (18) into Equation (16) and solving the obtained algebraic system after equating the coefficients of e r Γ ( ξ ) , to zero, yields
  • Case I
    a 0 = 1 3 ( λ 2 8 μ w ) , a 1 = 4 λ , a 2 = 4 , b 0 = α λ 2 + 8 α μ + 2 α w + w 6 α β b 1 = 2 6 α λ β , b 2 = 2 6 α β ,
  • Case II
    a 0 = 1 3 ( λ 2 2 μ w ) , a 1 = 2 λ , a 2 = 2 , b 0 = λ α λ 2 + 4 α μ 2 α w w 2 β b 1 = 2 α λ 2 + 4 α μ 2 α w w β , b 2 = 0 ,
Putting λ and μ so that λ 2 4 μ > 0 , gives
  • From Case I, the exact solutions are given by
V ( x , t ) = 1 3 ( λ 2 8 μ w ) + 8 μ λ ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 16 μ 2 ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 2 , H ( x , t ) = α λ 2 + 8 α μ + 2 α w + w 6 α β + 4 μ 6 α λ β ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 8 μ 2 6 α β ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 2 ,
  • From Case II
V ( x , t ) = 1 3 ( λ 2 2 μ w ) + 4 μ λ ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 8 μ 2 ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ 2 , H ( x , t ) = λ α λ 2 + 4 α μ 2 α w w 2 β + μ 8 α λ 2 + 4 α μ 2 α w w β ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 + w t ) + λ ,
Figure 1a,b present the traveling wave solutions Equation (21) with the parameters given by α = 1.2 , β = 1.2 , x 0 = 30 , w = 1.0 , μ = 0.5 , λ = 2.5 , t = 0 10 and in the interval [ 0 , 40 ] . Figure 2a,b show the traveling wave solutions Equation (22) depicted under the parameters α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and in the interval [ 0 , 20 ] . All of the results are displayed at t = 0 10 . The sufficient conditions of stability (Equations (12) and (13)) are applied, and then we find that all of the obtained solutions are stable in the interval [ 0 , 20 ] .
The boundary conditions are initially generated from studying the obtained travelling wave solutions at the endpoints. Figure 1a illustrates that the traveling wave solution V ( x , t ) goes to 1.75 at the endpoints of the physical domain and Figure 1b shows that H ( x , t ) goes to 2.2113 at the endpoints. While Figure 2a illustrates that the traveling wave solution V ( x , t ) goes to 1.3667 at the endpoints of the physical domain and Figure 2b shows that H ( x , t ) goes to ± 2.8142 at the endpoints. Hence, we can conclude that
Case I : V ( x , t ) 1.75 , and H ( x , t ) 2.2113 , as x ± , Case II : V ( x , t ) 1.3667 , and H ( x , t ) ± 2.8142 , as x ± .

4. Numerical Results

The MMPDEs technique [40] is implemented to extract the numerical solutions of the considered equations. This method can be used to deal with one-dimensional and multi-dimensional PDEs as shown in [40,41,42], respectively. The solutions of the coupled KdV equations are approximated by employing new meshes named gradient flow equations. The main idea of utilizing this approach is to generate new meshes and diminish the error in the solutions by feeding the regions with a high error with more points. The coupled KdV Equation (1) can be simply rewritten by
V t α V x x x 6 α V V x β H H x = 0 , H t + H x x x + 3 V H x = 0 ,
subject to the given boundary values in Equation (23). The corresponding initial data is given by
V ( x , 0 ) = 1 3 ( λ 2 2 μ w ) + 4 μ λ ( λ 2 4 μ ) t a n h 1 2 λ 2 4 μ ( x + x 0 ) + λ 8 μ 2 ( λ 2 4 μ ) t a n h ( 1 2 λ 2 4 μ ( x + x 0 ) + λ 2 , H ( x , 0 ) = λ α λ 2 + 4 α μ 2 α w w 2 β + μ 8 α λ 2 + 4 α μ 2 α w w β ( λ 2 4 μ ) t a n h ( 1 2 λ 2 4 μ ( x + x 0 ) + λ ,
The used coordinate transform is expressed as x = x ( ξ , t ) : [ 0 , 1 ] [ L , L ] , where x [ L , L ] and ξ [ 0 , 1 ] denote the space and computational coordinates, respectively. We use the chain rule to descritise the spatial and temporal derivatives as follows:
V t V x x t α V x x x 6 α V V x β H H x = 0 , H t H x x t + H x x x + 3 V H x = 0 .
Now, we rewrite Equation (26) as follows:
V t V ξ x ξ x t α V x x x 3 α V 2 ξ x ξ β H 2 ξ x ξ = 0 , H t H ξ x ξ x t + H x x x + 3 V H ξ x ξ = 0 ,
where
V x x x = 1 x ξ 1 x ξ V ξ x ξ ξ ξ , H x x x = 1 x ξ 1 x ξ H ξ x ξ ξ ξ .
To generate the mesh x ( ξ , t ) , we use the one dimensional equation [43,44,45]
δ ( I λ ξ ξ ) x t = 1 Ψ Ψ x ξ ξ ,
is utilized with the following boundary conditions
x 0 = L , x N + 1 = L ,
where L is half of the length of the physical domain and Ψ is the monitor function. The initial data that we use is taken
x ( ξ , t = 0 ) = ( n 1 ) Δ x L , Δ x = 2 L N n = 1 , 2 , 3 , , N + 1 .
The relaxation parameter is presented by δ and δ ( I λ ξ ξ ) is the smoothing operator described by Ceniceros and Hou [46]. The fundamental purpose of Ψ ( V , H , x ) is to control the movement of the mesh. A suitable choice for the function Ψ leads to perfect and acceptable results [43,44,45]. A modified monitor function was developed in [47] as follows:
Ψ ( V , H , x ) = 1 + α 1 x ξ V ξ x ξ ξ 2 + 1 x ξ H ξ x ξ ξ 2 .
where α is a positive constant. The semi-discretisation of the spatial derivative converts the coupled KdV equations into ( N + 1 ) ODEs easily solved using line methods. The obtained system is numerically integrated using the MATLAB ODE solver (ode15i). The physical domain is firstly separated as follows:
x n < x n + 1 , n = 0 , 1 , 2 , , N .
The computational coordinates are illustrated by
ξ n = n Δ ξ , Δ ξ = 1 N , n = 0 , 1 , 2 , , N .
We now fix the boundary mesh by x 0 = L and x N + 1 = L . The interior mesh is completely approximated by resolving Equation (29) with the boundary conditions x t , 0 = x t , N = 0 . The discretisation of Equation (28) is simply developed by applying the central finite differences. That is
V t i V i + 1 V i 1 x i + 1 x i 1 x t i = 2 α V x x , i + 1 V x x , i x i + 1 x i 1 3 α V i + 1 2 V i 1 2 x i + 1 x i 1 + β H i + 1 2 H i 1 2 x i + 1 x i 1 , H t i H i + 1 H i 1 x i + 1 x i 1 x t i = 2 H x x , i + 1 H x x , i x i + 1 x i 1 3 V i + 1 2 H i + 1 H i 1 x i + 1 x i 1 ,
where
H x x , i = 2 x i + 1 x i 1 H i + 1 H i x i + 1 x i H i H i 1 x i x i 1 , V i + 1 2 = V i + 1 + V i 2 , V x x , i = 2 x i + 1 x i 1 V i + 1 V i x i + 1 x i V i V i 1 x i x i 1 .
Figure 3 compares the movement of one travelling wave solution for both the exact solution of V ( x , t ) (left) and the numerical solution of V ( x , t ) (right) as the time increases from t = 0 15 . In Figure 4, we illustrate the movement of the analytic solution of H ( x , t ) (left) compared with the movement of the numerical solution of H ( x , t ) (right) for t = 0 15 . The used parameters are given by α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 . The results appear to have approximately the same behaviour. This indicates that the used methods are reliable and can be utilized to solve more nonlinear PDEs. The conditions of the stability are employed and we discover that our results are stable inside the interval [ 0 , 20 ] .

5. Comparison

The MMPDEs method utterly confirms its validity to be used in dealing with most NLPDEs. In particular, the method gives excellent results for the coupled KdV equations. More specifically, Figure 3 and Figure 4 illustrate a close similarity between the obtained exact and numerical solutions. Figure 5 illustrate that the solution moves as the time increases and seems to be without variations. A satisfactory performance of this approach can be also seen in Figure 6a,b where the numerical solutions mostly agree with the exact ones. This method distributes the points adequately. Therefore, the error is often very small. As a result, this technique can be used to solve several other NLPDEs.

6. Conclusions

In this work, we applied two methods to solve the coupled KdV equations, specifically, the analytical exp-function method and numerical MMPDEs method. We introduced some 2D and 3D plots for some of the obtained solutions. The obtained solutions are stable on the interval [0,20]. The achieved solutions may be applicable for explaining important phenomena in science, namely, in plasma physics and crystal lattices. We are quit sure that these techniques can also be applied to other NPDEs and, therefore, we shall take it up in our future studies.
Now, we summarize the article as follows.
  • Constructing some new solutions for coupled KdV equations, using the exp-function method.
  • Studying the modulation instability to check the stability of the exact solutions.
  • Implementing (numerically) the MMPDEs scheme to solve coupled KdV equations.
  • Currently, work is in progress on the application of the proposed two schemes for further coupled systems of NPDEs.

Author Contributions

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

Funding

The authors received no direct funding for this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdelrahman, M.A.E.; Kunik, M. The ultra-relativistic Euler equations. Math. Meth. Appl. Sci. 2015, 38, 1247–1264. [Google Scholar] [CrossRef]
  2. Abdelrahman, M.A.E. On the shallow water equations. Z. Naturforsch. 2017, 72, 873–879. [Google Scholar] [CrossRef]
  3. Abdelrahman, M.A.E.; Sohaly, M.A. The development of the deterministic nonlinear PDEs in particle physics to stochastic case. Results Phys. 2018, 9, 344–350. [Google Scholar] [CrossRef]
  4. Abdelrahman, M.A.E.; Sohaly, M.A.; Alharbi, A.R. The new exact solutions for the deterministic and stochastic (2+1)-dimensional equations in natural sciences. J. Taibah. Univ. Sci. 2019, 13. [Google Scholar] [CrossRef] [Green Version]
  5. Abdelrahman, M.A.E.; Hassan, S.Z.; Inc, M. The coupled nonlinear Schrödinger-type equations. Mod. Phys. Lett. B 2020, 2050078. [Google Scholar] [CrossRef]
  6. Yaslan, H.C.; Girgin, E. New exact solutions for the conformable space-time fractional KdV, CDG, (2+1)-dimensional CBS and (2+1)-dimensional AKNS equations. J. Taibah Sci. 2019, 13, 1–8. [Google Scholar] [CrossRef] [Green Version]
  7. Wazwaz, A.M.; Kaur, L. Optical solitons for nonlinear Schrödinger equation in the normal dispersive regimes. Optik 2019, 184, 428–435. [Google Scholar] [CrossRef]
  8. 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]
  9. Hassan, S.Z.; Abdelrahman, M.A.E. Solitary wave solutions for some nonlinear time fractional partial differential equations. Pramana-J. Phys. 2018, 91, 67. [Google Scholar] [CrossRef]
  10. Abdelrahman, M.A.E.; Sohaly, M.A. On the new wave solutions to the MCH equation. Indian J. Phys. 2018, 93, 903–911. [Google Scholar] [CrossRef]
  11. Cattani, C.; Ciancio, A. Hybrid two scales mathematical tools for active particles modeling complex systems with learning hiding dynamics. Math. Mod. Meth. Appl. Sci. 2007, 17, 171–188. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Cattani, Y.C.; Yang, X.J. Local fractional homotopy perturbation method for solving nonhomogeneous heat conduction equations in fractal domains. Entropy 2015, 17, 6753–6764. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, J.L.; Wang, M.L.; Wang, Y.M.; Fang, Z.D. The improved F-expansion method and its applications. Phys. Lett. A 2006, 350, 103–109. [Google Scholar] [CrossRef]
  14. Zhang, S.; Tong, J.L.; Wang, W. A generalized ( G G ) -expansion method for the mKdv equation with variable coefficients. Phys. Lett. A 2008, 372, 2254–2257. [Google Scholar] [CrossRef]
  15. Fan, E.; Zhang, H. A note on the homogeneous balance method. Phys. Lett. A 1998, 246, 403–406. [Google Scholar] [CrossRef]
  16. Dai, C.Q.; Zhang, J.F. Jacobian elliptic function method for nonlinear differential difference equations. Chaos Solut. Fractals 2006, 27, 1042–1049. [Google Scholar] [CrossRef]
  17. Aminikhad, H.; Moosaei, H.; Hajipour, M. Exact solutions for nonlinear partial differential equations via Exp-function method. Numer. Methods Partial. Differ. Equations 2009, 26, 1427–1433. [Google Scholar]
  18. 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]
  19. Wazwaz, A.M. The tanh method for travelling wave solutions of nonlinear equations. Appl. Math. Comput. 2004, 154, 714–723. [Google Scholar]
  20. Hydon, P.E. Symmetry Methods for Differential Equations; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  21. Yang, X.F.; Deng, Z.C.; Wei, Y. A Riccati-Bernoulli sub-ODE method for nonlinear partial differential equations and its application. Adv. Diff. Equ. 2015, 1, 117–133. [Google Scholar] [CrossRef] [Green Version]
  22. Abdelrahman, M.A.E.; Sohaly, M.A. Solitary waves for the nonlinear Schrödinger problem with the probability distribution function in stochastic input case. Eur. Phys. J. Plus 2017, 132, 339. [Google Scholar] [CrossRef]
  23. Huang, W.; Russell, R.D. The Adaptive Moving Mesh Methods; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  24. Budd, C.J.; Huang, W.; Russell, R.D. Adaptivity with moving grids. Acta Numer. 2009, 18, 111–241. [Google Scholar] [CrossRef]
  25. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  26. Gao, Y.T.; Tian, B. Ion-acoustic shocks in space and laboratory dusty plasmas: Two-dimensional and non-traveling-wave observable effects. Phys. Plasmas 2001, 8, 3146. [Google Scholar] [CrossRef]
  27. Lou, S.Y. Symmetries of the KdV equation and four hierarchies of the integro-differential KdV equations. J. Math. Phys. 1994, 35, 2390–2396. [Google Scholar] [CrossRef]
  28. Wang, G.; Xu, T. Symmetry properties and explicit solutions of the nonlinear time fractional KDV equation. Bound. Value Probl. 2013, 1, 1–13. [Google Scholar] [CrossRef] [Green Version]
  29. Hirota, R.; Satsuma, J. Soliton solutions of a coupled Korteweg-de Vries equation. Phys. Lett. A 1981, 85, 407–408. [Google Scholar] [CrossRef]
  30. Caom, D.B.; Yan, J.R.; Zang, Y. Exact solutions for a new coupled MKdV equations and a coupled KdV equations. Phys. Lett. A 2002, 297, 68–74. [Google Scholar]
  31. Ganji, D.D.; Rafei, M. Solitary wave solutions for a generalized Hirota-Satsuma coupled KdV equation by homotopy perturbation method. Phys. Lett. A 2006, 356, 131–137. [Google Scholar] [CrossRef]
  32. Gokdogan, A.; Yildirim, A.; Merdan, M. Solving coupled KdV equations by differential transformation method. World Appl. Sci. J. 2012, 19, 1823–1828. [Google Scholar]
  33. Assas, L.M.B. Variational iteration method for solving coupled KdV equations. Chaos Solitons Fractals 2008, 38, 1225–1228. [Google Scholar] [CrossRef]
  34. Jaradat, H.M.; Syam, M.; Alquran, M. A two-mode coupled Korteweg–de Vries: multiple-soliton solutions and other exact solutions. Nonlinear Dyn. 2017, 90, 371–377. [Google Scholar] [CrossRef]
  35. Hassan, S.H.; Alyamani, N.A.; Abdelrahman, M.A.E. A construction of new traveling wave solutions for the 2D Ginzburg-Landau equation. Eur. Phys. J. Plus 2019, 134, 425. [Google Scholar] [CrossRef]
  36. Chandrasekhar, S. Hydrodynamic and Hydromagnetic Stability; Dover Publications Inc.: New York, NY, USA, 1981. [Google Scholar]
  37. Sandstede, B. Chapter 18—Stability of Travelling Waves. In Handbook of Dynamical Systems; Fiedler, Ed.; Elsevier: Amsterdam, The Netherlands, 2002; Volume 2, pp. 983–1055. [Google Scholar]
  38. Jafari, H.; Soltani, R.; Khalique, C.M.; Baleanu, D. Exact solutions of two nonlinear partial differential equations by using the first integral method. Bound. Value Probl. 2013, 2013, 117. [Google Scholar] [CrossRef] [Green Version]
  39. Alharbi, A.R.; Almatrafi, M.B.; Abdelrahman, M.A.E. Analytical and numerical investigation for Kadomtsev-Petviashvili equation arising in plasma physics. Phys. Scr. 2020, 95, 4. [Google Scholar] [CrossRef]
  40. Russell, R.D.; Huang, W.; Ren, Y. Moving mesh partial differential equations (MMPDES) based on the equidistribution principle. Siam. J. Numer. Anal. 1994, 31, 709–730. [Google Scholar]
  41. Huang, W.; Russell, R.D. A high dimensional moving mesh strategy. Appl. Numer. Math. 1997, 26, 63–76. [Google Scholar] [CrossRef]
  42. Budd, C.J.; Williams, J.F. Moving mesh generation using the parabolic Monge-Ampere equation. SIAM J. Sci. Comput. 2009, 31, 3438–3465. [Google Scholar] [CrossRef]
  43. Alharbi, A.R.; Naire, S. An adaptive moving mesh method for thin film flow equations with surface tension. J. Comput. Appl. Math. 2017, 319, 365–384. [Google Scholar] [CrossRef]
  44. Alharbi, A.R.; Naire, S. An adaptive moving mesh method for two-dimensional thin film flow equations with surface tension. J. Comput. Appl. Math. 2019, 356, 219–230. [Google Scholar] [CrossRef] [Green Version]
  45. Alharbi, A.R. Numerical Solution of Thin-Film Flow Equations Using Adaptive Moving MESH methods. Ph.D. Thesis, Keele University, Keele, UK, 2016. [Google Scholar]
  46. Ceniceros, H.D.; Hou, T.Y. An efficient dynamically adaptive mesh for potentially singular solutions. J. Comput. Phys. 2001, 172, 609–639. [Google Scholar] [CrossRef] [Green Version]
  47. Cook, S.P. Adaptive Mesh Methods for Numerical Weather Prediction. Ph.D. Thesis, University of Bath, Bath, UK, 2016. [Google Scholar]
Figure 1. (a,b) depict the solutions (21) with α = 1.2 , β = 1.2 , x 0 = 30 , w = 1.0 , μ = 0.5 , λ = 2.5 , t = 0 10 and x = 0 40 .
Figure 1. (a,b) depict the solutions (21) with α = 1.2 , β = 1.2 , x 0 = 30 , w = 1.0 , μ = 0.5 , λ = 2.5 , t = 0 10 and x = 0 40 .
Symmetry 12 00429 g001
Figure 2. (a,b) depict the solutions (22) with α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Figure 2. (a,b) depict the solutions (22) with α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Symmetry 12 00429 g002
Figure 3. Time evolution ( t = 0 15 ) for the exact solution of V ( x , t ) is shown in (a) and the numerical solution of V ( x , t ) is illustrated in (b). The parameter values are α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Figure 3. Time evolution ( t = 0 15 ) for the exact solution of V ( x , t ) is shown in (a) and the numerical solution of V ( x , t ) is illustrated in (b). The parameter values are α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Symmetry 12 00429 g003
Figure 4. Time evolution ( t = 0 15 ) for the exact solution of H ( x , t ) is presented in (a) and the numerical solution of H ( x , t ) is depicted in (b). The parameter values are α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Figure 4. Time evolution ( t = 0 15 ) for the exact solution of H ( x , t ) is presented in (a) and the numerical solution of H ( x , t ) is depicted in (b). The parameter values are α = 2.5 , β = 2.5 , x 0 = 15 , w = 0.5 , μ = 0.10 , λ = 2.0 , t = 0 10 and x = 0 20 .
Symmetry 12 00429 g004
Figure 5. (a,b) depict the time development for one wave solution for V ( x , t ) and H ( x , t ) , respectively, at 0 t 20 . (c) presents the movement of the mesh x = x ( ξ , t ) .
Figure 5. (a,b) depict the time development for one wave solution for V ( x , t ) and H ( x , t ) , respectively, at 0 t 20 . (c) presents the movement of the mesh x = x ( ξ , t ) .
Symmetry 12 00429 g005
Figure 6. The movement of the exact and numerical solutions for V ( x , t ) and H ( x , t ) is plotted in (a) and (b), respectively, over the domain 0 x 20 . (c) exhibits the behaviour of the adaptive and uniform mesh.
Figure 6. The movement of the exact and numerical solutions for V ( x , t ) and H ( x , t ) is plotted in (a) and (b), respectively, over the domain 0 x 20 . (c) exhibits the behaviour of the adaptive and uniform mesh.
Symmetry 12 00429 g006

Share and Cite

MDPI and ACS Style

Abdelrahman, M.A.E.; Almatrafi, M.B.; Alharbi, A. Fundamental Solutions for the Coupled KdV System and Its Stability. Symmetry 2020, 12, 429. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030429

AMA Style

Abdelrahman MAE, Almatrafi MB, Alharbi A. Fundamental Solutions for the Coupled KdV System and Its Stability. Symmetry. 2020; 12(3):429. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030429

Chicago/Turabian Style

Abdelrahman, Mahmoud A. E., M. B. Almatrafi, and Abdulghani Alharbi. 2020. "Fundamental Solutions for the Coupled KdV System and Its Stability" Symmetry 12, no. 3: 429. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030429

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