Solitary wave solutions of the MRLW equation using a spatial five-point stencil of finite difference approximation

This paper proposes a finite difference scheme with a three-level time and a five-point stencil in space to solve an initial boundary value problem for the MRLW equation. The scheme is shown to be marginally stable and convergent with a fourth-order convergence in the space direction and a second-order convergence in the time variable direction with regard to the maximum norm. The conservation properties of the proposed scheme are assessed using the three motion invariants for mass, momentum, and energy. To validate the theoretical results, numerical experiments are given for both single and interaction of two and three solitary waves.


Introduction
Solitary waves, or solitons as they are also known, are nonlinear waves that have the ability to propagate through media over a prolonged period while maintaining key properties; for example, velocity and shape [1].As solitons remain stable when they collide with other solitons, they have found widespread application in a range of areas including optics, fluid mechanics, finance, biology, physics, engineering sciences, and neuroscience.Nonlinear partial differential equations are used to model solitons.For example, the regularized long wave (RLW) equation, which was originally presented by Peregrine [2] and Benjamin et al. [3], is as follows: where and are positive constants.It is considered with the homogeneous Dirichlet boundary conditions u → 0 as x → ±∞.
The RLW equation is particular case of the generalized long wave (GRLW) equation which has the form where p is a positive integer by setting p = 1 .In our paper, we consider another particular case of the GRLW equation called the modified regularized long wave (MRLW) equation when p = 2 and is given by The GRLW equation and its particulars cases: RLW equation and MRLW equation were initially suggested as a means of (1) modeling phenomena that exhibited dispersion waves in combination with weak nonlinearity; for example, pressure waves in a liquid gas bubble mixture, ion-acoustic and magnetohydrodynamic waves in plasma, phonon packets in nonlinear crystals, and nonlinear transverse waves in shallow water.The analytical solution that underpins the MRLW equation is limited to single solitary waves under the restricted initial and boundary conditions.Scholars have yet to develop formulae for other cases, for example, within the context of the Maxwellian initial condition and situations involving the interaction of more than one soliton.As such, the development of accurate numerical approximations remains critical such as finite difference methods [4][5][6][7][8][9], finite element method [10][11][12][13], mixed finite element method [14,15], collocation method [16][17][18], spectral method [19,20] and Adomian decomposition method [21,22].
Numerous finite difference schemes for the MRLW problem have been documented in the extant literature, with many scholars placing a specific focus on the finite difference methods and the MRLW equation.The MRLW problem, for instance, was solved using the finite difference approach by Khalifa et al. [6], who also investigated other aspects of the MRLW equation, including the interaction of solitary waves.Additionally, Fourier analysis was performed to demonstrate the stability of the scheme.The truncation error was also well controlled.
The generalized regularised long wave (GRLW) problem was solved by Akbari and Mokhtari using a new compact finite difference method (CFDM) [23].The stability analysis of the energy method was explored, and an error estimate was presented.The method was validated using two solitary waves interaction and the propagation of single solitons.To ascertain the method conservation properties, three motion invariants were assessed.
A fully implicit finite difference technique for the numerical solution of the MRLW equation was presented by Inan and Bahadir [24].The validity of the approach was tested using several MRLW equation examples.A comparison of the results with analytical and other numerical invariants demonstrated the accuracy and dependability of the outcomes achieved utilizing the fully implicit finite difference scheme.
In the current study, we suggest a finite difference approach to solve the MRLW equation that has three levels in time and a five-point stencil in space.The stability of the Fourier analysis-based method is considered, and the accuracy of the convergence rate of O(h 4 + k 2 ) is also proved.The remainder of this paper is structured as follows.The analytical solution of the MRLW equation and its conservative laws are reviewed in Sect. 2. The creation of the suggested scheme is the focus of Sect. 3. The stability and convergence rates of the scheme are examined in Sect. 4. To validate our theoretical findings, Sect. 5 presents various numerical experiments for single and interaction of solitons.Finally, our concluding remarks are contained in Sect.6.

Analytical solution and conservation laws
The analytical solution of the MRLW equation Eq. ( 3) is given in the form [6] where 6c is the amplitude of the MRLW solitary wave.
The solitary wave is initially centred at x 0 and its speed and its width are represented by c and √ c (c+1) , respectively.The validity of the numerical methods can be determined using the three invariants of the motion that the MRLW equation has; that is, the mass, momentum, and energy conservative laws given as [6,25]: and The use of these invariants can be particularly pertinent in situations for which there are no analytic solutions or during soliton interactions [25].

Construction of the finite difference scheme
To construct the finite difference scheme of the MRLW equation Eq. ( 3) with a three-level scheme in time and a five-point stencil in space, the following notations for the derivatives are used: where k = Δt is the time step and h = Δx represents the spatial step size.The superscript n denotes a quantity associated with time level t n and subscript j denotes a quantity associated with space mesh point x j .The grid points are t n = nk, n = 0, 1, 2, … , N for time and x j = jh, j = 0, 1, 2, … , M for space, where M and N are positive integers.The finite difference scheme therefore becomes Then substituting Eq. ( 8) into Eq.( 9) yields ( 8) where Q n j = kh(1 + (u n j ) 2 ) .The scheme ( 10) is a tridiagonal system that can be easily simulated using MATLAB platform.From now in our work, we take = 1.

Convergence and linear stability analysis
Lemma 1 The finite difference scheme (10) is marginally stable.Proof In the case of applying the Von Neumann stability theory, the solution of Eq. ( 10) can be written as where l is a mode number.Now, set and then inserting Eqs. ( 11) and ( 12) into Eq.( 10) yields where g is the growth factor and with assuming that (u n j ) 2 in Eq. ( 10) is locally constant, and for simplicity, we write it as u 2 .To verify Eq. ( 14), we must show that To do this, we have (10) for small spatial step size h and small time step k, and is practically taken to be unity.Now, Eq. ( 13) yields that g 1 = −e i and g 2 = e −i , which implies that ‖g 1 ‖ = ‖g 2 ‖ = 1 , and therefore the finite dif- ference scheme (10)  Proof Let v n j = v(x j , t n ) represents the exact solution for the Eq. ( 3) with independent variables x and t.The local truncation error of Eq. ( 9) is thus as follows: Now, using Tylor expansion, it is easily shown that T n j at the point (x j , t n ) can be written as and hence we have T  Subdomain FEM [12] Petrov Galerkin method [

Numerical experiments
In this section, we present some numerical experiments to verify our theoretical results obtained in the previous section.The accuracy of the proposed scheme is measured using the L ∞ and L 2 errors at t = t N that are approximated by In addition, the invariants of mass, momentum, and energy for the MRLW equation are calculated for a single soliton and during the interaction of two and three solitary waves to measure the conservation properties of the proposed scheme.

Motion of single solitary waves
With h = 0.4 and k = 0.05 fixed, two experiments were performed to demonstrate the viability of our scheme in the situation of a single soliton motion.Calculations up to T = 20 were performed.In line with the work of [23], the model parameters were selected as x 0 = 0 , c = 0.1 with range [−40, 60] for the first example and x 0 = 0 , c = 0.03 with range [−80, 120] for the second example.Tables 1 and  2 show the values of L ∞ and L 2 errors and the invariants I 1 , I 2 , and I 3 .The invariant I 1 is changed by less than 10 −5 in both circumstances, whilst the changes for the invariants I 2 and I 3 approach zero, demonstrating the reasonable conservatism of our proposed scheme.At T = 20 , the errors for the first example are reasonably minimal, at 1.182044 × 10 −4 for L ∞ error and 2.93952 × 10 −4 for L 2 error.Given that our approach is highly accurate, similar results are obtained for the second example as 2.09146 × 10 −5 for L ∞ error and 5.92125 × 10 −5 for L 2 error.In Fig. 1 for the first example, the motion of the single wave is plotted at various time levels with an amplitude of 0.3, and in Fig. 2 for the second case, with an amplitude of 0.17.The fact ( 18) that the waves at t = 16 and t = 20 adequately agree with those at t = 4 demonstrates the reliability and accuracy of our scheme.
Furthermore, we compare our results in terms of maximum errors with the obtained results in [8,9,12,13] to examine the validity of our scheme.For the purpose of comparisons, the parameters are chosen as h = 0.2 , k = 0.025 , c = 1 , x 0 = 40 , and = 6 , with a range [0, 100].The computations are performed up to T = 10 and are listed in Table 3.It is clearly observed that L ∞ -errors obtained by our scheme are marginally smaller than those obtained by others, indicating that our scheme is more accurate.

Convergence rate
To calculate the convergence rates in space and in time, we use the following formula [23] for the convergence rate in space, and Table 4 The convergence rates in space and maximum errors at T = 20 for single solitary wave of MRLW equation on [−40, 60] with c = 0.1 and  for the convergence rate in time, with respect to the maximum norm errors.With the other parameters remaining the same as in Table 1, we started with h = k = 0.8 , and reduced the spatial and temporal variables by 2 and 4, respectively, to calculate the spatial convergence rates.The resultant L ∞ errors and the corresponding convergence rates of our scheme are recorded in Table 4, which reveals that the fourth order of convergence in the spatial direction was achieved.Table 5, in contrast, displays the temporal convergence rates, as we first set h = k = 0.8 , then scaled them back by a factor of 2 for k and 4 for h.Other parameters are chosen as in Table 1.
Based on the proposed scheme, the accuracy of order 2 in the temporal direction is obtained.These spatial and temporal rates are aligned with the theoretical conclusions in Lemma 2. Additionally, the log-log scale depicted in Fig. 3. shows the resultant L ∞ errors on the solitary wave solutions with regard to spatial variable h and temporal variable k.In Fig. 3a and b, lines of slope 4 and 2 are added as references.As can be observed in these figures, our scheme achieved an O(h 4 + k 2 ) accuracy, coinciding with the theoretical predictions.

Interaction of solitary waves
The interaction of two and three solitary waves travelling in the same direction is discussed in this section.The initial conditions in these scenarios can be described by a linear sum of two and three well-separated solitary waves of different amplitudes, as follows: where p = 2 and p = 3 for the interaction of two and three solitary waves, respectively.We performed the simulations up to T = 200 and on the range [0, 300], with fixed h = 0.4 and k = 0.05 to enable the interaction to take place.Two solitary waves with c 1 = 4 21 , c 2 = 9 91 , x 1 = 15 and x 2 = 35 interacted at various time levels, as shown Table 7 Invariants for interaction of three solitary waves on [0, 300] with h = 0.   7 lists the three invariants for this case.The outputs of the experiments reveal that the three invariants held relatively steady throughout the interaction process.Our approach effectively maintains the soliton properties because the waves interact and maintain their shape.

Conclusion
In this paper, we described the use of a finite difference method based on a three-level temporal scheme and a five-point space stencil to solve the initial boundary-value problem for the MRLW equation.Based on the maximum norm, the scheme is shown to be marginally stable and convergent with second-order accuracy in time and fourth-order accuracy in space.In order to demonstrate the effectiveness and accuracy of the proposed method, numerical experiments using single and solitary waves interaction were described.Investigations into the conservation quantities for mass, momentum, and energy were also conducted, and the results were deemed satisfactory.Comparisons with other previous results are given to show the accuracy and the efficiency of the proposed scheme.

Table 5
The convergence rates in time and maximum errors at T = 20 for single solitary wave of MRLW equation on [−40, 60] with c = 0.1 and x 0 = 0 The L ∞ errors with respect to h at time T = 20 with parameters chosen as in Table4.b The L ∞ errors with respect to k at time T = 20 with parameters chosen as in