Abstract

In this paper, we present new numerical results for the dispersive optical soliton solutions of the nonlinear Schrödinger-Hirota equation. The spatio-temporal dispersion term is included, in addition to group velocity dispersion Kerr law of nonlinearity are studied. A general recursive numerical scheme for the equation is devised via the Improved Adomian Decomposition Method (IADM) and further sought for some analytical results for validation. The scheme is shown to be efficient and possessed high level of accuracy as demonstrated.

1. Introduction

Nonlinear Schrodinger equations (NLSE) play a significant role in the fiber optics, optical soliton, prorogation of pulses in metamaterials and network engineering applications among others [1, 2]. Thus, different forms of NLSE exist and serve for a variety of purposes mostly in communication and networking engineering including, for example, the trans-continental and trans-oceanic data transfer requiring phase modulation [3] and beyond; see [416] for different mathematical studies on these equations. Furthermore, due to their immense application, many analytical methods and few numerical methods have been proposed in the past decades. Such methods include the Backlund transformation method [17], Hirota's direct method [18, 19], tanh-sech method [20, 21], extended tanh method [22, 23], sine-cosine method [24], and homogeneous balance method [25].

However, we will propose in this paper a numerical method for solving a class of NLSE called the nonlinear Schrodinger-Hirota equation. The method will heavily depend on the Improved Adomian Decomposition Method (IADM) [2631]. The IADM is an efficient numerical method for functional and integral solutions based on the Adomian decomposition method [32]. Further, a recent analytical study on the nonlinear Schrodinger-Hirota equation will be sought for to validate the proposed scheme. This form of solution is being reported for the first time in this paper.

2. Governing Equation

The nonlinear Schrodinger equation used in modeling propagation of solitons through optical fibers with Third-Order Dispersion (TOD) is given bywhere is a complex-valued function of (space) and (time); is the coefficient of TOD. The first term on the left hand side is the linear temporal evolution, the second term is the group velocity dispersion term, and the third term accounts for Kerr law nonlinearity. Also when the group velocity is low, the TOD is justified.

We now introduced the Lie symmetry concept to be able to study (1). With , let which transforms (1) toIn (3), higher order terms are neglected [11, 33, 34] with Kerr law nonlinearity that models transmission of dispersive optical solitons through nonlinear fibres. We express (3) with general coefficients as where physically represents nonlinear dispersion. Further, (4) happens to be ill-posed due to the group velocity dispersion term. However, with the addition of the spatiotemporal dispersion (STD) term, (4) possesses a well-posedness [35, 36] given by where the coefficient represents STD. Finally, in presence of perturbation terms (5) becomes [3739] However, we consider an analytical 1-soliton solution of (6) as special case presented by Al Qarni et al. [40] for the numerical simulation sake in Section 4 given by where is the frequencies of the solitons, is the wave number, is the phase velocity, and and are the amplitude and width given, respectively, byThe constraint conditions, for the existence of bright solitons from (7)-(9), are as follows:

3. Numerical Method

In [2631], authors introduced the IADM to convert a special case of the complex-valued system into a real-valued system by splitting as where are real functions. By substituting equation (12) into (6) we obtain the following system:whereThe decomposition method [32] decomposes the solution into infinite sums of components defined bywhere the components will be determined recursively. In an operator form with , (14) and (15) becomeTaking the inverse operator to both sides of (18) and (19) gives Assuming that the nonlinear terms in (20) and (21) are represented by the series, are the so-called Adomian polynomials that can be constructed for all forms of nonlinearity according to specific algorithms set by Adomian [32]. Substituting the nonlinear terms into (22) and (23) and the solution from (16) and (17) into (20) and (21) givesFollowing the decomposition analysis, the following recursive relations are introduced:Thus, we determine and as follows: and the overall approximate solution for (6) is obtained by substituting the above into (12) coupled to (26)-(29) to get

4. Numerical Results

In this section, we consider three different cases for the Schrödinger-Hirota equation with spatiotemporal dispersion given in (6) to illustrate the application of the IADM scheme we presented in the above section. We also consider the bright soliton solution given in (7)-(11) for numerical simulation with the following fixed parameters: , , , , . We also represent the graphs for absolute value for each of exact solution and approximate solution .

Example 1. We consider Schrödinger-Hirota equation with spatiotemporal dispersion equation (6) with and the following three special cases.

Case 1. Let

Case 2. Let .

Case 3. Let

The results and the profiles of this example are presented in Table 1 and Figure 1.

Example 2. We consider Schrödinger-Hirota equation with spatiotemporal dispersion equation (6) with and the following three special cases.

Case 1. Let .

Case 2. Let .

Case 3. Let .

The result and the profiles of this example are presented in Table 2 and Figure 2.

Example 3. We consider Schrödinger-Hirota equation with spatiotemporal dispersion equation (6) with and the following three special cases.

Case 1. Let .

Case 2. Let .

Case 3. Let .

The result and the profiles of this example are presented in Table 3 and Figure 3.

5. Conclusion

In this work, an Improved Adomian Decomposition Method (IADM) has been proposed to solve the nonlinear Schrodinger-Hirota equation in presence of several Hamiltonian type perturbation terms. The obtained results possess high precision and converged to the exact solution with less computational efforts. It was noticed that a high accuracy of results is obtained when is the coefficient of 3OD and within the range , ( is the coefficient of spatio-temporal dispersion term) as commented also in [40]. It is also clear from Figures 13 that the complete correspondence of the two solutions is indeed remarkable. Thus, the IADM is an efficient method for nonlinear Schrodinger equations since it shows a high level of accuracy with lesser computational efforts as compared to other numerical methods. This topic is still open for further researches for different values of parameters and other types of soliton solutions.

Data Availability

All the data used for the numerical simulations and comparison purpose have been reported in the tables included and visualized in the graphical illustrations and nothing is left.

Conflicts of Interest

The authors declare that they have no conflicts of interest.