Next Article in Journal
Novel Composite Speed Control of Permanent Magnet Synchronous Motor Using Integral Sliding Mode Approach
Previous Article in Journal
Lump Waves in a Spatial Symmetric Nonlinear Dispersive Wave Model in (2+1)-Dimensions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Numerical and Analytical Solutions for Nonlinear Evolution Equations Using Updated Mathematical Methods

by
Abdulghani R. Alharbi
Department of Mathematics, College of Science, Taibah University, Al-Madinah Al-Munawarah 42353, Saudi Arabia
Mathematics 2023, 11(22), 4665; https://doi.org/10.3390/math11224665
Submission received: 27 September 2023 / Revised: 13 November 2023 / Accepted: 15 November 2023 / Published: 16 November 2023

Abstract

:
This study explores adapted mathematical methods to solve the couple-breaking soliton (BS) equations in two-dimensional spatial domains. Using these methods, we obtained analytical soliton solutions for the equations involving free parameters such as the wave number, phase component, nonlinear coefficient, and dispersion coefficient. The solutions are expressed as hyperbolic, rational, and trigonometric functions. We also examined the impact of wave phenomenon on two-dimensional diagrams and used composite two-dimensional and three-dimensional graphs to represent the solutions. We used the finite difference method to transform the proposed system into a numerical system to obtain numerical simulations for the Black–Scholes equations. Additionally, we discuss the stability and error analysis of numerical schemes. We compare the validity and accuracy of the numerical results with the exact solutions through analytical and graphical comparisons. The methodologies presented in this research can be applied to various forms of nonlinear evolutionary systems because they are appropriate and acceptable.
MSC:
65N06; 65N40; 65N45; 65N50; 35A24; 35B35; 35Q51; 35Q92

1. Introduction

One of the most challenging problems in mathematical physics is the quest to find numerical and soliton solutions for the nonlinear PDEs (NPDEs). NPDEs are utilized extensively for the purpose of doing analyses on nonlinear processes. Nonlinear partial differential equations (NPDEs) have practical applications in various scientific fields. Solutions are obtained numerically and analytically using methods such as ( P / Q ) -expansion, Hirota bilinear, and exp-function [1,2,3,4]. In order to acquire soliton solutions for NLPDEs, a variety of approaches were utilized in order to examine the equations that were presented. The authors of [5] analyze and discuss the regularized long-wave and adjusted equal-width (MEW) equations. Wazwaz studied the KP-MEW equations using the Tanh approach in 2004. Radha and Lakshmanan examined localized structures in a pair BS equation and the Painleve property in 1995. This equation is of some assistance. In 2004, Wazwaz conducted research on the KP-MEW equations by utilizing the Tanh approach, which utilizes the (MW) equation represented in the PK purpose and is described in [6,7]. In recent decades, Radha et al. [8] researched the presence of exponentially localized systems in a (2 + 1)-dimensional pair BS equation and the Painleve property. For the two-and-a-half-dimensional plus one-dimensional BS equation, Yan and Zhang [9] created families with a solution that is analogous to a soliton. After a period of one year had passed, Chen et al. [10] utilized the recently obtained generalized expansion approach of the Riccati issue in order to achieve a soliton solution to the (2 + 1)-dimensional pair BS equation. Peng and Krishnan [11] used the singular manifold approach to examine two new kinds of exact soliton solutions for the two-and-a-half-dimensional pair BS issue. Inan [12] conducted research on the two-and-a-half-dimensional pair BS equation using the generalized Jacobi elliptic function approach. As a result of his investigation, periodic and multiple soliton solutions were obtained for the above equation. In addition to the residual symmetry that was related to the truncated Painleve expansion, Cheng et al. [13] were able to construct the nonlocal symmetry that was gained from the Lax pair. On the basis of the generalized unified technique, Osman [14] was able to obtain the multiple-soliton solutions to the two-dimensional pair BS problem. These solutions were derived for the problem. In recent years, a substantial number of researchers have examined a wide variety of analytical methodologies for the (2 + 1)-dimensional pair BS model. These researchers come from a variety of academic fields. The 2-dimensional pair BS problem has been handled by utilizing several different methodologies, such as the Hirota bilinear technique [15], the resemblance modification process [16], the cos-Gordon expansion method [17], and a few more. These many approaches have been combined in order to produce a single answer. The numerical solution to the two-dimensional fourth-order parabolic thin film equation was derived by Alharbi [18] using the parabolic Monge–Ampere approach in more recent research. Ren and Chu [19] looked at the soliton molecules and explored the generalized BS equations in two dimensions and one dimension. According to the best of our knowledge, researchers have yet to previously document in the prior literature a solution to the ( 2 + 1 ) -dimensional Z K -MEW equations and a handful of BS equations utilizing the method that is offered here. In addition, the information provided in this paper can be utilized as a supplement for other publications connected to the same overall topic area. In addition, we have begun working on a solution that incorporates hyperbolic and trigonometric function solutions with free parameters. The solutions to these issues will be extremely useful in a wide variety of sectors, including, but not limited to, scientific research, maritime engineering, and many others. In more recent research, Alharbi [18] used the parabolic Monge–Ampere technique to derive the numerical solution to the two-dimensional fourth-order parabolic thin film equation. Ren and Chu [19] investigated the generalized BS equations in two dimensions and one dimension and looked at the soliton molecules. To the best of our knowledge, researchers have yet to previously document in the literature a solution to the ( 2 + 1 ) -dimensional Z K -MEW equations and a handful of BS equations utilizing the method that is provided here. In addition, the material presented in this paper can serve as a supplement to other publications that are related to the same overarching subject area. In addition, we have started working on a solution with hyperbolic and trigonometric function solutions that incorporate free parameters. The answers to these problems will find widespread use in scientific research, maritime engineering, and a variety of other fields. The Kudryashov technique employs the transformed rational function approach [20] and the coefficient of variation for generalized breaking solitons [21]. Kaplan and Kumar [20], Kaplan and associates [22], and Kaplan and Akbulut [23] investigated the generalized Schrodinger equation, the Jaulent–Miodek equation, and fractional PDEs, respectively, all in the same year. Alharbi et al. [24] successfully found a solution to the Kadomtsev–Petviashvili problem using the adaptive moving mesh technique. Ma [25] researched nonlocally integrable NLPDEs, calculating the soliton solutions to the nonlocally integrable modified Korteweg–de Vries problem, and investigating the inverse scattering changes and soliton solutions of nonlocal integrable equations utilizing the Riemann–Hilbert problems that were made available through one-group reduction. Researchers have also discovered accurate solutions to the BS system, including double solitons and other linked forms [26], and solutions supported by trigonometry and exponential development [27]. The BS system [13,14] is considered as follows:
Γ t +   4 β Γ Ψ x + 4 β Γ x Ψ + β Γ x x y = 0 , Ψ x   Γ y = 0 .
Bogoyovlenskii [28] initially introduced this equation. Calogero and Degasperis [29,30] employed the (2 + 1)-dimensional couple BS equation to express the interaction of a Riemann wave traveling down the y-axis with a long wave traveling along the x-axis in 1996. They did this by applying the equation to the Riemann wave’s matter. The wave events that are described in the equations are essential to maintaining the integrity of water wave events and maintaining tight contact. It may open the gate to further research into nonlinear wave phenomena, which are significant in the field of ocean engineering. Kazeykina and Klein [31] carried out an investigation in which they analyzed the stability of solutions for a particular equation and provided numerical solutions to the equation. In order to determine the numerical consequences of the proposed system, they also used the Kansa approach [32]. On the other hand, there has been no discussion addressing the system’s stability or an error analysis of the numerical method that was used (1). Obtaining various analytic solutions to a problem (1) by utilizing modified S expansion and generalized algebraic approaches was the primary purpose of this research. In order to accomplish this goal, I combined the numerical solution with finite differences in order to produce numerical results for the system that was researched. I made a substantial contribution to the comprehension of physical challenges in practice by analytically and visually comparing the solutions to traveling wave problems and the numerical outcomes of those problems. The remaining parts of this article are arranged as follows: Section 2 presents a general overview of the mathematical models. Section 3 discusses the methodology that we developed for extracting results from NLPDEs. The numerical schemes used to discover the numerical results of the proposed system (1) are presented in Section 4, together with details regarding the stability, accuracy, and convergence of those numerical schemes. After that, we derive a variety of solutions to the collection of differential equations by making use of those equations. The graphical and physical explanations of the solutions we identified will be discussed in the following section of this article. A comparison of the answers is displayed in Section 5, and the conclusion is presented in Section 6.

2. Overview of Proposed Procedures

The following form represents the development equations with the physical fields Γ ( x , y , t ) and Ψ ( x , y , t ) in the variables x , y, and t.
G 1 ( Γ t , Γ x , Γ y , Γ x x y , Ψ , Ψ y , Ψ x ) = 0 , G 2 ( Ψ x , Γ y , Γ x x y , Ψ , Ψ y , Ψ x ) = 0 .
Step 1.
We obtain the System (1) traveling-wave solutions, which are created as follows:
Γ ( x , y , t ) = Φ ( ζ ) , ζ = x + y α t , Ψ ( x , y , t ) = Θ ( ζ ) .
where α is the wave’s speed.
Step 2.
The following ODE represents the nonlinear evolution (2):
G 3 ( Φ , Φ ζ , Φ ζ ζ ζ , Θ , Θ ζ , ) = 0 , G 4 ( Θ , Θ ζ , Φ , Φ ζ , ) = 0 .
where G 3 and G 4 are polynomials with respect to Θ ( ζ ) , Φ ( ζ ) , and their combined derivatives.
The generalized indirect algebraic approach suggests that the solution to (4) is
Φ ( ζ ) = ρ 0 + k = 1 N ρ k ψ ( ζ ) k + s k ψ ( ζ ) k ,
where ψ ( ζ ) is a solution of
ψ ( ζ ) = k = 0 4 γ k ψ k ( ζ ) 1 2 ,
where ρ k , and s k are to be determined, and N is a positive integer that balances the highest possible degree of the nonlinear terms and the highest order of the derivatives. The relationships between γ k and 0 k 4 are listed in Table 1 [33].
The Modified S—the expansion method yields the following solutions to Equation (4):
ψ ( ζ ) = λ 0 + k = 1 N λ k S ( ζ ) k + q k S ( ζ ) k ,
where S ( ζ ) is a solution of
S ( ζ ) = μ 0 + μ 1 S ( ζ ) + μ 2 S ( ζ ) 2
where μ k and 0 k 2 are listed in Table 2 [34] and q k and λ k are to be obtained later.

3. Methodology

Consider the couple-breaking soliton (BS) Equation (1). By applying the appropriate transformations, the equation denoted as Equation (1) can be simplified into a system of ordinary differential equations (ODEs) in the following manner:
Γ ( x , y , t ) = Φ ( ζ ) , ζ = x + y α t , Ψ ( x , y , t ) = Θ ( ζ ) .
As a result, Equation (1) is presented as
α Φ ζ + 4 β Φ Θ ζ + 4 β Φ ζ Θ + β Φ ζ ζ ζ = 0 , Θ ζ Φ ζ = 0 .
Integrating the second equation in (10) with respect to ζ yields Θ = Φ . Thus, the first equation in (10) is given by
α Φ + 4 β Φ 2 + β Φ ζ ζ = 0 .
In (11), the value of N = 2 is determined by balancing Φ ζ ζ with Φ 2 .
Based on the modified S expansion method with N = 2 , the solution to (11) is
ψ ( ζ ) = λ 0 + λ 1 S ( ζ ) + q 1 S ( ζ ) + λ 2 S ( ζ ) 2 + q 2 S ( ζ ) 2 ,
where S ( ζ ) is a solution of
F ( ζ ) = μ 0 + μ 1 S ( ζ ) + μ 2 S ( ζ ) 2 .
In addition, μ i , i = 0 , 1 , 2 , are listed in Table 2. In order to examine the analytical solutions to (11), the steps outlined below should be followed.
(1)
A system of equations for λ 0 , λ m , q m , where m = 1 , 2 , can be obtained by incorporating (12) and (13) into Equation (11) and setting all of the coefficients of S ( ζ ) m , 4 m 4 to zero;
(2)
Mathematical software, such as Mathematica 13.2 or Maple 2023.1, should be employed to solve the resulting system;
(3)
Using the values of μ 0 , μ 1 and μ 2 , and the function S ( ζ ) from Table 2, and substituting them along with λ 0 , λ k , q k , k = 1 , 2 into Equation (12), several trigonometric functions and rational solutions to Equation (11) are obtained.
Using the preceding steps, I determined the following values for λ 0 , λ 1 , λ 2 , q 1 , q 2 , and α :
1.
When μ 0 = 0 , μ 1 = ± 1 , and μ 2 = 1 , the solutions of System (1) Γ are provided using
Γ 1 ( x , y , t ) = 1 4 + 3 4 1 + tanh ( ( x + y + β t ) / 2 ) 3 8 1 + tanh ( ( x + y + β t ) / 2 ) 2 , Ψ 1 ( x , y , t ) = 1 8 1 3 tanh 2 1 2 ( β t + x + y ) , Γ 2 ( x , y , t ) = 1 4 + 3 4 1 + coth ( ( x + y + β t ) / 2 ) 1 1 2 1 + coth ( ( x + y + β t ) / 2 ) , Ψ 2 ( x , y , t ) = 1 8 1 3 coth 2 1 2 ( β t + x + y ) , Γ 3 ( x , y , t ) = 3 4 1 + tanh ( ( x + y β t ) / 2 ) 3 8 1 + tanh ( ( x + y β t ) / 2 ) 2 , Ψ 3 ( x , y , t ) = 3 8 tanh 2 1 2 ( β t + x + y ) 1 , Γ 4 ( x , y , t ) = 3 4 1 + coth ( ( x + y + β t ) / 2 ) 1 1 2 1 + coth ( ( x + y + β t ) / 2 ) , Ψ 4 ( x , y , t ) = 1 8 ( 3 ) coth 2 1 2 ( β t + x + y ) 1 .
2.
When μ 0 = 1 , μ 1 = 0 , and μ 2 = 1 , two cases exist. The solutions are
Γ 5 ( x , y , t ) = 3 2 3 2 tanh 2 ( ( x + y 4 β t ) ) , Ψ 5 ( x , y , t ) = 3 2 sech 2 ( 4 β t + x + y ) , Γ 6 ( x , y , t ) = 3 2 3 2 coth 2 ( ( x + y 4 β t ) ) , Ψ 6 ( x , y , t ) = 3 2 csch 2 ( 4 β t + x + y ) , Γ 7 ( x , y , t ) = 1 3 2 tanh 2 ( ( x + y + 16 β t ) ) + coth 2 ( ( x + y + 16 β t ) ) , Ψ 7 ( x , y , t ) = 3 2 tanh 2 ( 16 β t + x + y ) + coth 2 ( 16 β t + x + y ) 1 , Γ 8 ( x , y , t ) = 1 2 3 2 tanh 2 ( 4 β t + x + y ) , Ψ 8 ( x , y , t ) = 1 2 1 3 tanh 2 ( 4 β t + x + y ) , Γ 9 ( x , y , t ) = 1 2 3 2 coth 2 ( 4 β t + x + y ) , Ψ 9 ( x , y , t ) = 1 2 1 3 coth 2 ( 4 β t + x + y ) , Γ 10 ( x , y , t ) = 3 3 2 tanh 2 ( 16 β t x y ) 3 2 coth 2 ( 16 β t x y ) , Ψ 10 ( x , y , t ) = 6 csch 2 ( 2 ( 16 β t + x + y ) ) .
3.
When μ 0 = 1 2 , μ 1 = 0 , and μ 2 = 1 2 , the solutions are
Γ 11 ( x , y , t ) = 3 8 1 + sec x + y + β t ± tan x + y + β t 2 , Ψ 11 ( x , y , t ) = 3 8 ( sec ( β t + x + y ) ± tan ( β t + x + y ) ) 2 3 8 Γ 12 ( x , y , t ) = 1 4 3 8 sec x + y + 4 β t + tan x + y + 4 β t 2 + sec x + y + 4 β t + tan x + y + 4 β t 2 , Ψ 12 ( x , y , t ) = 1 8 3 ( tan ( 4 β t + x + y ) + sec ( 4 β t + x + y ) ) 2 3 ( tan ( 4 β t + x + y ) + sec ( 4 β t + x + y ) ) 2 + 2 .
4.
When μ 0 = 1 2 , μ 1 = 0 , and μ 2 = 1 2 , the solutions are
Γ 13 ( x , y , t ) = 1 8 3 tanh 2 1 2 ( 4 β t + x + y ) 3 coth 2 1 2 ( 4 β t + x + y ) 2 , Ψ 13 ( x , y , t ) = 3 8 tanh 2 1 2 ( 4 β t + x + y ) 1 8 3 coth 2 1 2 ( 4 β t + x + y ) + 1 4 , Γ 14 ( x , y , t ) = 1 8 3 8 ( coth ( β t + x + y ) + csch ( β t + x + y ) ) 2 , Ψ 14 ( x , y , t ) = 1 8 ( cosh ( β t + x + y ) + 2 ) csch 2 1 2 ( β t + x + y ) , Γ 15 ( x , y , t ) = 1 8 3 8 coth ( β t + x + y ) + csch ( β t + x + y ) 2 , Ψ 15 ( x , y , t ) = 1 8 1 3 coth ( β t + x + y ) + csch ( β t + x + y ) 2 , Γ 16 ( x , y , t ) = 1 8 1 3 coth 2 1 2 ( β t + x + y ) , Ψ 16 ( x , y , t ) = 1 8 3 8 coth 2 1 2 ( β t + x + y ) , Γ 17 ( x , y , t ) = 3 8 sech 2 1 2 ( β t x y ) , Ψ 17 ( x , y , t ) = 3 4 1 + cosh ( x + y β t ) 1 .
Based on the generalized algebraic method with N = 2 , the solution to (11) is
Φ ( ζ ) = ρ 0 + ρ 1 ψ ( ζ ) + ρ 2 ψ ( ζ ) 2 + s 1 ψ ( ζ ) + s 2 ψ ( ζ ) 2 ,
where ψ ( ζ ) is a solution of
ψ ( ζ ) = γ 0 + γ 1 ψ ( ζ ) + γ 2 ψ 2 ( ζ ) + γ 3 ψ 3 ( ζ ) + γ 4 ψ 4 ( ζ ) .
Table 1 contains a list of all possible values for γ j , j = 0 , 1 , 2 , 3 , 4 . To determine the values of the coefficients ρ 0 , ρ 1 , ρ 2 , s 1 , s 2 , and α in all of the scenarios mentioned earlier and their subsequent solutions, I used mathematical software, specifically Mathematica 13.2, to obtain the values. As a result, the analytic solutions to Equation (13) will be shown in this section using the generalized algebraic method with various alternative values for γ k , k = 0 , 1 , 3 , 4 .
4.
If γ 0 = 0 , γ 1 = 0 , γ 2 > 0 , γ 3 = 0 and γ 4 < 0 , then
Γ 18 ( x , y , t ) = 3 2 γ 2 sech 2 γ 2 x + y 4 β γ 2 t , Γ 19 ( x , y , t ) = γ 2 + 3 2 γ 2 sech 2 γ 2 x + y 4 β γ 2 t .
Ψ 18 ( x , y , t ) = 3 2 γ 2 sech 2 γ 2 x + y 4 β γ 2 t , Ψ 19 ( x , y , t ) = γ 2 + 3 2 γ 2 sech 2 γ 2 x + y 4 β γ 2 t .
5.
If γ 0 = γ 2 2 4 γ 4 4 , γ 1 = 0 , γ 2 < 0 , γ 3 = 0 , and γ 4 > 0 , then
Γ 20 ( x , y , t ) = 3 γ 2 csch 2 2 γ 2 8 β γ 2 t + x + y , Γ 21 ( x , y , t ) = 3 4 γ 2 tanh 2 γ 2 2 β γ 2 t + x + y 2 1 , Γ 22 ( x , y , t ) = 1 4 γ 2 3 tanh 2 γ 2 ( 2 β γ 2 t + x + y ) 2 1 Γ 23 ( x , y , t ) = 3 4 γ 2 csch 2 γ 2 2 β γ 2 t + x + y 2 , Γ 24 ( x , y , t ) = 1 4 γ 2 3 tanh 2 γ 2 ( 8 β γ 2 t + x + y ) 2 + 3 coth 2 γ 2 ( 8 β γ 2 t + x + y ) 2 + 2 .
Ψ 20 ( x , y , t ) = 3 γ 2 csch 2 2 γ 2 8 β γ 2 t + x + y , Ψ 21 ( x , y , t ) = 3 4 γ 2 tanh 2 γ 2 2 β γ 2 t + x + y 2 1 , Ψ 22 ( x , y , t ) = 1 4 γ 2 3 tanh 2 γ 2 ( 2 β γ 2 t + x + y ) 2 1 Ψ 23 ( x , y , t ) = 3 4 γ 2 csch 2 γ 2 2 β γ 2 t + x + y 2 , Ψ 24 ( x , y , t ) = 1 4 γ 2 3 tanh 2 γ 2 ( 8 β γ 2 t + x + y ) 2 + 3 coth 2 γ 2 ( 8 β γ 2 t + x + y ) 2 + 2 .
6.
If γ 0 = 0 , γ 1 = 0 , γ 2 > 0 , γ 3 0 , and γ 4 = 0 , then
Γ 25 ( x , y , t ) = 3 8 γ 2 sech 2 1 2 γ 2 ( β γ 2 t + x + y ) , Γ 26 ( x , y , t ) = γ 2 4 + 3 8 γ 2 sech 2 1 2 γ 2 ( β γ 2 t + x + y ) .
Ψ 25 ( x , y , t ) = 3 8 γ 2 sech 2 1 2 γ 2 ( β γ 2 t + x + y ) , Ψ 26 ( x , y , t ) = γ 2 4 + 3 8 γ 2 sech 2 1 2 γ 2 ( β γ 2 t + x + y ) .
7.
If γ 0 = 0 , γ 1 = 0 , γ 2 < 0 , γ 3 = 0 , γ 4 > 0 , then
Γ 27 ( x , y , t ) = 3 2 γ 2 sec 2 γ 2 ( 4 β γ 2 t + x + y ) , Γ 28 ( x , y , t ) = 3 2 γ 2 sec 2 γ 2 ( 4 β γ 2 t + x + y ) γ 2 .
Ψ 27 ( x , y , t ) = 3 2 γ 2 sec 2 γ 2 ( 4 β γ 2 t + x + y ) , Ψ 28 ( x , y , t ) = 3 2 γ 2 sec 2 γ 2 ( 4 β γ 2 t + x + y ) γ 2 .
8.
If γ 0 = γ 2 2 4 γ 4 , γ 1 = 0 , γ 2 > 0 , γ 3 = 0 , γ 4 > 0 , then
Γ 29 ( x , y , t ) = 3 4 γ 2 tan 2 γ 2 ( 8 β γ 2 t + x + y ) 2 1 4 3 γ 2 cot 2 γ 2 ( 8 β γ 2 t + x + y ) 2 3 γ 2 2 , Γ 30 ( x , y , t ) = 1 4 3 γ 2 tan 2 γ 2 ( 2 β γ 2 t + x + y ) 2 3 γ 2 4 , Γ 31 ( x , y , t ) = 1 4 3 γ 2 cot 2 γ 2 ( 2 β γ 2 t + x + y ) 2 3 γ 2 4 , Γ 32 ( x , y , t ) = 1 4 3 γ 2 tan 2 γ 2 ( 2 β γ 2 t + x + y ) 2 γ 2 4 , Γ 33 ( x , y , t ) = 1 4 3 γ 2 cot 2 γ 2 ( 2 β γ 2 t + x + y ) 2 γ 2 4 ,
Ψ 29 ( x , y , t ) = 3 4 γ 2 tan 2 γ 2 ( 8 β γ 2 t + x + y ) 2 1 4 3 γ 2 cot 2 γ 2 ( 8 β γ 2 t + x + y ) 2 3 γ 2 2 , Ψ 30 ( x , y , t ) = 1 4 3 γ 2 tan 2 γ 2 ( 2 β γ 2 t + x + y ) 2 3 γ 2 4 , Ψ 31 ( x , y , t ) = 1 4 3 γ 2 cot 2 γ 2 ( 2 β γ 2 t + x + y ) 2 3 γ 2 4 , Ψ 32 ( x , y , t ) = 1 4 3 γ 2 tan 2 γ 2 ( 2 β γ 2 t + x + y ) 2 γ 2 4 , Ψ 33 ( x , y , t ) = 1 4 3 γ 2 cot 2 γ 2 ( 2 β γ 2 t + x + y ) 2 γ 2 4 ,
where the solution Ψ is equal to Γ in the corresponding cases and β is an arbitrary constant. In this discussion, we will cover the BS equation, which involves a constant known as β . These equations describe the interaction between a long wave along the x axis and a Riemann wave that moves along the y axis. α is used to represent the speed of the singularity. Figure 1 shows surface plots of the soliton solutions (a) Γ 1 ( x , y , t ) and (b) Γ 3 for t = 0 50 . Additionally, Figure 2 depicts the 3D wave profiles for the soliton solutions (a) Γ 5 and (b) Γ 17 with β = 1.2 and t = 0 10 . The surface plots in the solution represent a brilliant-type wave profile.

4. Numerical Solution

We employed various numerical techniques to ensure accurate analytical results when solving the ODEs (11). Using Γ 7 as an example, we set Φ and Θ to 0 at the endpoint where ζ is 0 and then guessed the starting values of Φ ζ and Θ ζ . At t = 0 , we applied the nonlinear shooting and BVP methods to find the second Φ = 0 boundary condition at a specific left endpoint of the domain. Once we obtain the numerical data, we compared them to the analytical solution Γ 7 at t = 0 . We generated the numerical solution using FSOLVE and ODE15s in MATLAB [35]. The discretized version of the resulting ODE (11), which is
α Φ + 2 β Θ 2 + 2 β Φ 2 + β δ ζ 2 Φ i + 1 2 Φ i + Φ i 1 = 0 ,
was subjected to the shooting method. Figure 3 displays a comparison of the analytic and numerical solutions obtained using the aforementioned numerical approaches, all of which were identical. The accuracy of the analytical solution could, therefore, be checked using this method. The acquired numerical solution can serve as an initial condition for the numerical scheme discussed in the following section.
To calculate the numerical outcomes of System (1) within the rectangular domain [ a , b ] × [ 0 , 1 ] , I will be using the finite differences method. The starting point and ending point of the rectangle along the x axis are represented by a and b , respectively. The rectangular area enclosed by [ a , b ] × [ 0 , 1 ] is divided into ( N + 1 ) × ( J + 1 ) mesh points. These points are arranged as
x i = a + i Δ x , i = 0 , 1 , 2 , , N , y j = j Δ y , j = 0 , 1 , 2 , , J ,
where Δ x and Δ y are the x a n d y domain step sizes, respectively. System (1) is turned into an ODE system by discretizing the space derivatives while keeping the time derivative continuous. Completing this process yields the desired results.
Γ t | i , j n = 2 β Δ x Γ i + 1 2 , j n + 1 ( Ψ i + 1 , j n + 1 Ψ i 1 , j n + 1 ) + Ψ i + 1 2 , j n + 1 ( Γ i + 1 , j n + 1 Γ i 1 , j n + 1 ) β 2 Δ y δ x 2 ( Γ i , j + 1 n + 1 Γ i , j 1 n + 1 ) , 0 = 1 2 Δ x ( Ψ i + 1 , j n + 1 Ψ i 1 , j n + 1 ) 1 2 Δ y ( Γ i , j + 1 n + 1 Γ i , j 1 n + 1 ) ,
where
δ x 2 Γ i , j n + 1 = Γ i + 1 , j n + 1 2 Γ i , j n + 1 + Γ i 1 , j n + 1 .
According to the boundary conditions
Γ x ( a , y , t ) = Γ x ( b , y , t ) = 0 , y , Γ y ( x , 0 , t ) = Γ y ( x , 1 , t ) = 0 , x ,
The boundary constraints enabled us to evaluate the space derivatives at the domain’s endpoints using fictitious points. The following factors contributed to the initial conditions:
Γ 10 ( x , y , t = 0 ) = 3 2 γ 2 sech 2 γ 2 x + y + x 0 ,
where γ 2 > 0 is the setting chosen by the user. The parameter values are held constant throughout all of the numerical results presented in this section as β = 0.1 , x 0 = 5.0 , y = 0 1 , x = 0 20 and t = 0 25 . It is important to note that we utilized the DDASPK solver [36], a FORTRAN software 95 that effectively solves ODEs, to solve System (31). The solver employed a regressive differentiation formula that necessitated an LU-factorization approximation of the Jacobian matrix of the linearized system, as there was no initial condition for the space derivatives. The numerical outcomes we acquired are satisfactory, and the figures below demonstrate their comprehension.

5. Stability of the Numerical Scheme

The stability of the numerical solution was investigated using the Fourier stability technique. This technique employs Fourier analysis, also known as von Neumann analysis, to assess the stability of the scheme (31). It is important to note that this method is only suitable for situations that involve linear equations. Therefore, the goal was to transform the given equations linearly.
Γ t +   4 β Γ Ψ x + 4 β Γ x Ψ + β Γ x x y = 0 , Ψ x   Γ y = 0 .
By examining the second equation of System (34), it can be observed that Γ is equal to Ψ . Hence, the first equation in (34) can be expressed as
Γ t + 8 β Γ Γ x + β Γ x x y = 0 .
where β is a constant. In order to employ the Fourier stability approach, it is important to possess linear equations. Hence, in order to align Equation (35) with this approach, it is necessary to convert it into a linear equation.
Γ t + d 0 Γ x + β Γ x x y = 0 ,
where d 0 = 8 β Γ is quantized and held constant by
d 0 = max 1 i N 1 j J 8 β Γ i , j n .
Equation (36) can be expressed using the finite difference method.
Γ i , j n = Γ i , j n + 1 + D 1 Γ i + 1 , j n + 1 D 1 Γ i 1 , j n + 1 + D 2 Γ i + 1 , j + 1 n + 1 2 Γ i , j + 1 n + 1 + Γ i 1 , j + 1 n + 1 D 2 Γ i + 1 , j 1 n + 1 2 Γ i , j 1 n + 1 + Γ i 1 , j 1 n + 1 ,
where D 1 = 0.5 d 0 Δ t Δ x , D 2 = 0.5 β Δ t Δ y . It is essential to consider the boundary conditions as neglecting them can lead to significant consequences. Assume that x i = i Δ x ,   y j = j Δ y and t n = n Δ t ; then,
Γ i , j n = μ n e ι ξ π i Δ x e ι η π j Δ y and then Γ i , j n + 1 = μ Γ i , j n , i = 1 , 2 , 3 , , N , j = 1 , 2 , 3 , , M and n .
By substituting (38) into (36), we obtain the subsequent results:
Γ i , j n = μ Γ i , j n + 2 D 1 sin ( Δ x ξ π ) μ Γ i , j n 8 i D 2 sin 2 Δ x ξ π 2 sin η π Δ y μ Γ i , j n .
Thus,
μ = 1 1 + 2 D 1 sin ( ξ π Δ x ) 8 i D 2 sin 2 Δ x ξ π 2 sin Δ y η π .
Let
v = 2 D 1 sin ( Δ x ξ π ) 8 i D 2 sin 2 Δ x ξ π 2 sin Δ y η π .
Equation (40) can be rewritten as:
μ = 1 1 + v 1 .
Based on our comprehensive study, we can assert with confidence that the scheme under consideration exhibits stability and adheres to the established norms of Fourier stability. Maintaining a value of mu below one is crucial, and we can verify that our system successfully meets this criterion. The equation shown in (40) provides clear evidence that the absolute value of μ is less than one, ensuring the stability of our numerical system under all circumstances.

6. Error Analysis

To determine the accuracy order of Scheme (31), I utilized the Taylor series. The order was calculated based on the evaluation of the truncation error. Here, we assume that
e i , j n + 1 = Γ i , j n + 1 Γ x i , y j , t n + 1 .
In order to study the error analysis at a specific position and time ( x i , y j , t n + 1 ) , we utilize e i , j n + 1 , which denotes the error. The estimated solution is referred to as Γ i , j n , while the analytical solution is Γ x i , y j , t n + 1 . The ultimate outcome can be achieved by incorporating (43) into scheme (37).
e i , j n = e i , j n + 1 + D 0 e i + 1 , j n + 1 D 0 e i 1 , j n + 1 + D 1 e i + 1 , j + 1 n + 1 2 e i , j + 1 n + 1 + e i 1 , j + 1 n + 1 D 1 e i + 1 , j 1 n + 1 2 e i , j 1 n + 1 + e i 1 , j 1 n + 1 Δ t T i , j n ,
where D 0 = d 0 Δ t 2 Δ x and D 1 = β Δ t 2 Δ y . In addition, T i , j n presents a way to express the truncation error, which is as follows:
T i , j n = Γ x i , y j , t n + 1 + D 0 Γ x i + 1 , y j , t n + 1 D 0 Γ x i 1 , y j , t n + 1 + D 1 Γ x i + 1 , y j + 1 , t n + 1 2 Γ x i , y j + 1 , t n + 1 + Γ x i 1 , y j + 1 , t n + 1 D 1 Γ x i + 1 , y j 1 , t n + 1 2 Γ x i , y j 1 , t n + 1 + Γ x i 1 , y j 1 , t n + 1 ,
Hence,
T i , j n Δ t 2 Γ t t x i , y j , η 0 n + 1 Δ x Δ y 4 Γ x x x y y η 1 i , η 2 j , t n + 1 Δ x 2 6 Γ x x x η 1 i , y j , t n + 1 ,
where η 0 n + 1 t n , t n + 1 , η 1 i x i 1 , x i + 1 , and η 2 j y j 1 , y j + 1 . Thus, the numerical scheme’s truncation error is T i , j n = O ( Δ t ) + O ( Δ y 2 ) + O Δ x 2 . Consequently, every step’s truncation error is indicated by
T i , j n = O Δ t , Δ y 2 , Δ x 2 .
Performing a series of computations using the initial data and refining it with distinct meshes, gives ( Δ x , Δ y , Δ t 0 ) .

7. Convergence of the Numerical Schemes

In order to determine whether a numerical scheme is convergent within a specific domain D , the first step is to evaluate each fixed point ( x * , y * , t * ) D . If ( x i , y j , t n ) ( x * , y * , t * ) , then Γ i , j n = Γ ( x * , y * , t * ) and the numerical scheme is convergent. Earlier, I demonstrated that the implicit scheme guarantees the unconditional stability of the system. In order to confirm that its convergence is unconditional, let us assume that the mistake is represented by e, which is made possible by
e i , j n = Γ i , j n Γ x i , y j , t n .
Despite the truncation error indication, Γ i , j n removes it and successfully fulfilled the strategy outlined in Equation (31). Thus,
e i , j n = e i , j n + 1 + D 1 e i + 1 , j n + 1 D 1 e i 1 , j n + 1 + D 2 e i + 1 , j + 1 n + 1 2 e i , j + 1 n + 1 + e i 1 , j + 1 n + 1 D 2 e i + 1 , j 1 n + 1 2 e i , j 1 n + 1 + e i 1 , j 1 n + 1 Δ t T i , j n ,
where D 1 = 0.5 d 0 Δ t Δ x , D 2 = 0.5 β Δ t Δ y . Assuming that the following is the correct way to define the maximum error for each time step:
E n = max e i , j n , i , j and n 0 .
Consequently, the result of implementing Equation (49) within Equation (48) is
e i , j n + 1 E n + Δ t T i , j n ,
Subsequently, for every value of i = 1 , , N , j = 1 , , M gives
E n + 1 E n + Δ t T i , j n .
Based on the original data, it is possible to determine that E 0 = 0 . Therefore, the inequality is expressed as
E 1 E 0 + Δ t T i , j n = Δ t T i , j n E 2 E 1 + Δ t T i , j n 2 Δ t T i , j n .
Thus, we have
E n n × Δ t T i , j n .
When Δ x , Δ y , Δ t 0 , the value of T i , j n 0 ; hence,
E k k × Δ t T n , m k 0 as Δ t 0 .
Therefore, it may be concluded that the scheme (Equation (31)) exhibits convergence when Δ x , Δ y , Δ t 0 .

8. Result and Discussion

This study focuses on the significant outcomes and methods employed. The exact solutions for Equation (1) can be obtained using the generalized indirect algebraic and modified S expansion methods. Figure 1, Figure 2 and Figure 4 display some solitary wave solutions for the parameters β = 1.20 , N = 100 , and M = 3200 , respectively. I also used a numerical scheme (31) to confirm these solutions and obtain numerical results. Figure 5 and Table 3 present the main findings, facilitating the comparison between the analytical and numerical solutions. Based on the above analysis, it is evident that the solutions exhibit a high degree of similarity, and the error significantly diminishes as Δ x and Δ y approach zero. Any value of parameter β must remain unchanged for the numerical schemes to maintain stability. Any alteration to this parameter may lead to instability and should be avoided at all costs. Figure 1a,b present the time evolution of the analytic solutions (a) Γ 1 with x = 70 10 and (b) Γ 3 with x = 10 70 . These figures are plotted at t = 0 50 and β = 1.2 . The firing method efficiently addresses the stated problem. Figure 2a,b present the time evolution of the analytic solutions (a) Γ 5 with x = 10 60 and (b) Γ 8 with x = 10 60 . These figures are plotted at t = 0 10 and β = 1.2 .  Figure 4a,b present the time evolution of the analytic solutions (a) Γ 11 with x = 1 1 and (b) Γ 17 with x = 10 60 . These figures are plotted at t = 0 10 , and β = 1.2 . To create the numerical results of Equation (31), the solution was used as an initial condition for the numerical scheme. Figure 6 shows the exact solution for a single traveling wave and the shooting method solutions. These solutions are mutually conducive and behave as a unit, indicating that the methods employed are valid and successful. The Taylor series expansion was also used to examine the scheme’s precision. The precision was determined beginning with the second order. The proposed numerical technique is only approximately accurate to the second order, as indicated in Figure 5. The Von Neumann analysis demonstrates the unconditional stability of the numerical system.
Figure 7a,b depict the time progression of the exact and numerical results with β = 0.1 , x 0 = 5.0 , y = 0 1 , x = 0 20 , t = 0 25 , J = 100 , and N = 3200 . It can be seen that the constructed accurate traveling solutions match with the numerical solutions. The resulting accurate solitary traveling wave and numerical solutions at t = 25 are displayed in 3 D in Figure 8a,b. These graphs make it clear that the exact solutions exhibit behaviors similar to those of the numerical solutions. Table 3 and Figure 5 also demonstrate the efficacy of the techniques employed. The L 2 error and the amount of CPU time needed to arrive to time step 5 using the numerical scheme are shown in Table 3. Increasing the number of points dramatically reduces inaccuracy and increases the CPU time required to process them. The finite difference method is more suitable for solving nonlinear partial differential equations computationally because it provides the high-error locations with sufficient data points, quickly decreasing the error.

9. Conclusions

In this study, we used the generalized indirect algebraic and modified S-expansion methods to derive various types of solitary wave solutions for the BS equations. These solutions were presented in the form of hyperbolic and trigonometric functions. To solve the suggested problem, we used the shooting method to construct a solution, which we then used as an initial condition for a numerical scheme. The numerical solution to the problem was achieved successfully without any issues. We conducted an in-depth analysis of the differences and similarities between the analytical and numerical results to validate the constructed solutions. The solutions behaved almost identically to one another. The finite differences approach achieved a significantly higher accuracy, contributing to the method’s significant improvement. We determined that the numerical scheme is unconditionally stable, and there is a correlation between the total number of points and the relative error. From a computational standpoint, the methods utilized are effective.

Funding

The authors extend their appreciation to the Deanship of Research and Innovation at the Ministry of Education in Saudi Arabia for funding this research through project number 445-9-535.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The author declares that he has no potential conflicts of interest.

References

  1. Akbulut, A.; Islam, S.; Rezazadeh, H.; Tascan, F. Obtaining exact solutions of nonlinear partial differential equations via two different methods. Int. J. Mod. Phys. B 2022, 36, 2250041. [Google Scholar] [CrossRef]
  2. Wen-Xiu, M.; Yong, X.; Zhang, H. Diversity of interaction solutions to the (2+1)-dimensional Ito equation. Comput. Math. Appl. 2018, 75, 289–295. [Google Scholar] [CrossRef]
  3. Ozkam, Y.; Yasar, E.; Osman, M. Novel multiple soliton and front wave solutions for the 3D-Vakhnenko–Parkes equation. Mod. Phys. Lett. B 2022, 36, 2250003. [Google Scholar] [CrossRef]
  4. Bashar, M.; Islam, S. Exact solutions to the (2 + 1)-Dimensional Heisenberg ferromagnetic spin chain equation by using modified simple equation and improve F-expansion methods. Phys. Open 2020, 5, 100027. [Google Scholar] [CrossRef]
  5. Zaki, S. Solitary wave interactions for the modified equal width equation. Comput. Phys Commun. 2000, 126, 219–231. [Google Scholar] [CrossRef]
  6. Wazwaz, A. The tanh method and the sine–cosine method for solving the KP-MEW equation. Int. J. Comput. Math. 2005, 82, 235–246. [Google Scholar] [CrossRef]
  7. Alharbi, A.; Faisal, M.; Shah, F.; Waseem, M.; Ullah, R.; Sherbaz, S. Higher Order Numerical Approaches for Nonlinear Equations by Decomposition Technique. IEEE Access 2019, 7, 44329–44337. [Google Scholar] [CrossRef]
  8. Radha, R.; Lakshmanan, M. Dromion like structures in the (2 + 1)-dimensional breaking soliton equation. Phys. Lett. A 1995, 197, 7–12. [Google Scholar] [CrossRef]
  9. Yan, Z.; Zhang, H. Constructing families of soliton-like solutions to a (2+1)-dimensional breaking soliton equation using symbolic computation. Int. J. Comput. Math. Appls. 2002, 44, 1439–1444. [Google Scholar] [CrossRef]
  10. Chen, Y.; Biao, L.; Zhang, H. Symbolic Computation and Construction of Soliton-Like Solutions to the (2+1)-Dimensional Breaking Soliton Equation. Commun. Theor. Phys. 2003, 40, 137–142. [Google Scholar] [CrossRef]
  11. Peng, Y.; Krishnan, E. Two Classes of New Exact Solutions to (2+1)-Dimensional Breaking Soliton Equation. Commun. Theor. Phys. 2005, 44, 807–809. [Google Scholar] [CrossRef]
  12. Inan, I. Generalized Jacobi Elliptic Function Method for Traveling Wave Solutions of (2+1)-Dimensional Breaking Soliton Equation. Cankaya Univ. J. Sci. Eng. 2010, 7, 39–50. [Google Scholar]
  13. Cheng, W.; Chen, Y. Nonlocal symmetry and exact solutions of the (2+1)- dimensional breaking soliton equation. Commun. Nonlinear Sci. Numer. Simulat. 2015, 29, 198–207. [Google Scholar] [CrossRef]
  14. Osman, M.S. On multi-soliton solutions for the (2+1)-dimensional breaking soliton equation with variable coefficients in a graded-index waveguide. Comput. Math. Appls. 2018, 75, 1–6. [Google Scholar] [CrossRef]
  15. Manafian, J.; Behnam, M.; Abapour, M. Lump-type solutions and interaction phenomenon to the (2+1)-dimensional Breaking Soliton equation. Appl. Math. Comput. 2019, 356, 13–41. [Google Scholar] [CrossRef]
  16. Kumar, M.; Tanwar, D. Lie symmetries and invariant solutions of (2 + 1)-dimensional breaking soliton equation. Pranama J. Phys. 2020, 94, 23. [Google Scholar] [CrossRef]
  17. Baskonus, H.; Kumar, A.; Wei, G. Deeper investigations of the (4 + 1)-dimensional Fokas and (2 + 1)-dimensional Breaking soliton equations. Int. J. Mod. Phys. B 2020, 34, 2050152. [Google Scholar] [CrossRef]
  18. Alharbi, A. Numerical solutions to two-dimensional fourth order parabolic thin film equations using the Parabolic Monge-Ampere method. AIMS Math. 2023, 8, 16463–16478. [Google Scholar] [CrossRef]
  19. Ren, B.; Chu, P. Dynamics of D’Alembert wave and soliton molecule for a (2+1)-dimensional generalized breaking soliton equation. Chin. J. Phys. 2021, 74, 296–301. [Google Scholar] [CrossRef]
  20. Kaplan, M.; Akbulut, A. The analysis of the soliton-type solutions of conformable equations by using generalized Kudryashov method. Opt. Quantum. Electron. 2021, 53, 498. [Google Scholar] [CrossRef]
  21. Qin, Y.; Gao, Y.; Shen, Y.; Sun, Y.; Meng, G.; Yu, X. Solitonic interaction of a variable coefficient (2 + 1)-dimensional generalized Breaking Soliton equation. Phys. Scr. 2013, 88, 1–7. [Google Scholar] [CrossRef]
  22. Mirzazadeh, M.; Hosseini, K.; Dehingia, K.; Salahshour, S. A second-order nonlinear Schrödinger equation with weakly nonlocal and parabolic laws and its optical solitons. Optic 2021, 242, 166911. [Google Scholar] [CrossRef]
  23. Xia, T.; Xiong, S. Exact solutions of (2 + 1)-dimensional Bogoyavlenskii’s Breaking Soliton equation with symbolic computation. Comput. Math. Appl. 2010, 60, 919–923. [Google Scholar] [CrossRef]
  24. Alharbi, A.; Almatrafi, M.; Abdelrahman, M. Analytical and numerical investigation for Kadomtsev–Petviashvili equation arising in plasma physics. Phys. Scr. 2020, 95, 045215. [Google Scholar] [CrossRef]
  25. Alharbi, A. A Study of Traveling Wave Structures and Numerical Investigation of Two-Dimensional Riemann Problems with Their Stability and Accuracy. Comput. Model. Eng. Sci. 2023, 134, 2193–2209. [Google Scholar] [CrossRef]
  26. Cao, L.; Wang, D.; Chen, L. Symbolic computation and q-deformed function solutions of (2 + 1)-dimensional Breaking Soliton equation. Commun. Theor. Phys. 2007, 47, 270–274. [Google Scholar] [CrossRef]
  27. Zhang, S. A generalized new auxiliary equation method and its application to the (2 + 1)-dimensional Breaking Soliton equations. Appl. Math. Comput. 2007, 190, 510–516. [Google Scholar] [CrossRef]
  28. Bogoyavlensky, O. Overturning solitons in new two-dimensional integrable equations. Izv. Akad. Nauk SSSR Ser. Mat. 1989, 53, 243–257. [Google Scholar] [CrossRef]
  29. Calogero, F.; Degasperis, A. Nonlinear evolution equations solvable by the inverse spectral transform—I. Il Nuovo Cimento B 1976, 32, 201–242. [Google Scholar] [CrossRef]
  30. Calogero, F.; Degasperis, A. Nonlinear evolution equations solvable by the inverse spectral transform— II. Il Nuovo Cimento B 1977, 39, 1–54. [Google Scholar] [CrossRef]
  31. Kazeykina, A.; Klein, C. Numerical study of blow-up and stability of line solitons for the Novikov- Veselov equation. Nonlinearity 2017, 30, 2566. [Google Scholar] [CrossRef]
  32. Sagar, B.; Saha, S. Numerical soliton solutions of fractional (2+1)-dimensional Nizhnik-Novikov-Veselov equations in nonlinear optics. Int. J. Mod. Phys. B 2021, 35, 2150090. [Google Scholar] [CrossRef]
  33. Bai, C.; Bai, C.; Zhao, H. A new generalized algebraic method and its application in nonlinear evolution equations with variable coefficients. Z. Naturforsch. A 2005, 60, 211–220. [Google Scholar] [CrossRef]
  34. Aasaraai, A. The application of modified F-expansion method solving the Maccari’s system. J. Adv. Math. Comput. Sci. 2015, 11, 1–14. [Google Scholar] [CrossRef]
  35. Shampine, L.; Reichelt, M. The matlab ode suite. SIAM J. Sci. Comput. 1994, 18, 1–22. [Google Scholar] [CrossRef]
  36. Brown, P.; Hindmarsh, A.; Petzold, L.E. Using Krylov methods in the solution of large-scale differential-algebraic systems. SIAM J. Sci. Comput. 1994, 15, 1467–1488. [Google Scholar] [CrossRef]
Figure 1. The time evolution of the exact solutions (a) Γ 1 , and (b) Γ 3 . The parameter is supplied by β = 1.20 , with t = 0 : 10 : 50 .
Figure 1. The time evolution of the exact solutions (a) Γ 1 , and (b) Γ 3 . The parameter is supplied by β = 1.20 , with t = 0 : 10 : 50 .
Mathematics 11 04665 g001
Figure 2. The time evolution of the exact solutions (a) Γ 5 and (b) Γ 17 . The parameter is supplied by β = 1.20 , with t = 0 : 2 : 10 .
Figure 2. The time evolution of the exact solutions (a) Γ 5 and (b) Γ 17 . The parameter is supplied by β = 1.20 , with t = 0 : 2 : 10 .
Mathematics 11 04665 g002
Figure 3. Comparison of Γ 7 at t = 0 via analytical and numerical methods. With parameters β = 1.2 , x = 0 20 , and N = 3200 .
Figure 3. Comparison of Γ 7 at t = 0 via analytical and numerical methods. With parameters β = 1.2 , x = 0 20 , and N = 3200 .
Mathematics 11 04665 g003
Figure 4. The time evolution of the exact solutions (a) Γ 5 and (b) Γ 17 . The parameter is supplied by β = 1.20 , with t = 0 : 2 : 10 .
Figure 4. The time evolution of the exact solutions (a) Γ 5 and (b) Γ 17 . The parameter is supplied by β = 1.20 , with t = 0 : 2 : 10 .
Mathematics 11 04665 g004
Figure 5. The convergence records based on Δ x utilizing the l 2 norm relative error from Table 3). As for the variable y, the value of 0.5 at t = 25 and x = 0 20 was specifically selected.
Figure 5. The convergence records based on Δ x utilizing the l 2 norm relative error from Table 3). As for the variable y, the value of 0.5 at t = 25 and x = 0 20 was specifically selected.
Mathematics 11 04665 g005
Figure 6. Change in the time frame for the numerical results. Maintaining the constant values of y = 0.5 at t = 25 is crucial. Upon close observation of the wave, it is clear that there is a close resemblance between the numerical and analytic solutions.The parameter values are determined by β = 0.1 , x 0 = 5.0 , x = 0 20 , t = 0 25 , J = 100 , and N = 3200 .
Figure 6. Change in the time frame for the numerical results. Maintaining the constant values of y = 0.5 at t = 25 is crucial. Upon close observation of the wave, it is clear that there is a close resemblance between the numerical and analytic solutions.The parameter values are determined by β = 0.1 , x 0 = 5.0 , x = 0 20 , t = 0 25 , J = 100 , and N = 3200 .
Mathematics 11 04665 g006
Figure 7. The analytical (a) and numerical (b) solutions for Γ 10 are shown in 3D. The parameter values are determined by β = 0.1 , x 0 = 5.0 , y = 0 1 , x = 0 20 , t = 0 25 , J = 100 , and N = 3200 .
Figure 7. The analytical (a) and numerical (b) solutions for Γ 10 are shown in 3D. The parameter values are determined by β = 0.1 , x 0 = 5.0 , y = 0 1 , x = 0 20 , t = 0 25 , J = 100 , and N = 3200 .
Mathematics 11 04665 g007
Figure 8. The resulting (a) accurate solitary traveling wave and (b) numerical solutions at t = 25 with β = 1.2 , x = 0 20 , y = 0 1 , N = 3200 , and J = 100 .
Figure 8. The resulting (a) accurate solitary traveling wave and (b) numerical solutions at t = 25 with β = 1.2 , x = 0 20 , y = 0 1 , N = 3200 , and J = 100 .
Mathematics 11 04665 g008
Table 1. The relationships between γ k , 0 k 4 , and the function ψ ( ζ ) .
Table 1. The relationships between γ k , 0 k 4 , and the function ψ ( ζ ) .
γ 0 , γ 1 , γ 2 , γ 3 , γ 4 ψ ( ζ )
γ 0 = γ 2 2 4 γ 4 , γ 1 = 0 , γ 2 < 0 , γ 3 = 0 , γ 4 > 0 ψ ( ζ ) = γ 2 2 γ 4 tanh γ 2 2 ζ .
γ 0 = 0 , γ 1 = 0 , γ 2 > 0 , γ 3 = 0 , γ 4 < 0 ψ ( ζ ) = γ 2 γ 4 sech γ 2 ζ .
γ 0 = γ 2 2 4 γ 4 , γ 1 = 0 , γ 2 > 0 , γ 3 = 0 , γ 4 > 0 ψ ( ζ ) = γ 2 2 γ 4 tan γ 2 2 ζ .
γ 0 = 0 , γ 1 = 0 , γ 2 = 0 , γ 3 = 0 , γ 4 > 0 , ψ ( ζ ) = 1 γ 4 ζ .
γ 0 = 0 , γ 1 = 0 , γ 2 < 0 , γ 3 = 0 , γ 4 > 0 ψ ( ζ ) = γ 2 γ 4 sec γ 2 ζ .
γ 0 = 0 , γ 1 = 0 , γ 2 > 0 , γ 3 0 , γ 4 = 0 ψ ( ζ ) = γ 2 γ 3 sech 2 γ 2 2 ζ .
Table 2. The relationships between μ k , 0 k 2 , and the function S ( ζ ) .
Table 2. The relationships between μ k , 0 k 2 , and the function S ( ζ ) .
μ 0 , μ 1 , μ 2 S ( ζ )
μ 0 = 0.5 , μ 1 = 0.0 , μ 2 = 0.5 S ( ζ ) = sec ( ζ ) + tan ( ζ ) , csc ( ζ ) cot ( ζ ) .
μ 0 = 0.5 , μ 1 = 0.0 , μ 2 = 0.5 S ( ζ ) = sec ( ζ ) tan ( ζ ) , csc ( ζ ) + cot ( ζ ) .
μ 0 = ± 1 , μ 1 = 0.0 , μ 2 = ± 1 S ( ζ ) = tan ( ζ ) , cot ( ζ ) .
μ 0 = 0.0 , μ 1 = 1.0 , μ 2 = 1.0 S ( ζ ) = 1 2 1 + tanh 1 2 ζ .
μ 0 = 1.0 , μ 1 = 0.0 , μ 2 = 1.0 S ( ζ ) = tanh ( ζ ) , coth ( ζ ) .
μ 0 = 1 / 2 , μ 1 = 0.0 , μ 2 = 1 / 2 S ( ζ ) = tanh ( ζ ) ± sech ( ζ ) , coth ( ζ ) ± csch ( ζ ) .
Table 3. The relative error formula RLTV is used here to measure the error with fixing Δ t = 1 × 10 2 at a specific time t = 25 .
Table 3. The relative error formula RLTV is used here to measure the error with fixing Δ t = 1 × 10 2 at a specific time t = 25 .
Δ x RLTVCPU
1.0 × 10 1 1.5 × 10 3 56.0 s
5.0 × 10 2 3.2 × 10 4 133.36 s
2.5 × 10 2 1.6 × 10 4 316.20 s
1.25 × 10 2 4.4 × 10 5 619.80 s
6.3 × 10 3 1.2 × 10 5 1492.86 s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alharbi, A.R. New Numerical and Analytical Solutions for Nonlinear Evolution Equations Using Updated Mathematical Methods. Mathematics 2023, 11, 4665. https://doi.org/10.3390/math11224665

AMA Style

Alharbi AR. New Numerical and Analytical Solutions for Nonlinear Evolution Equations Using Updated Mathematical Methods. Mathematics. 2023; 11(22):4665. https://doi.org/10.3390/math11224665

Chicago/Turabian Style

Alharbi, Abdulghani R. 2023. "New Numerical and Analytical Solutions for Nonlinear Evolution Equations Using Updated Mathematical Methods" Mathematics 11, no. 22: 4665. https://doi.org/10.3390/math11224665

APA Style

Alharbi, A. R. (2023). New Numerical and Analytical Solutions for Nonlinear Evolution Equations Using Updated Mathematical Methods. Mathematics, 11(22), 4665. https://doi.org/10.3390/math11224665

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