In this study, the exact solutions for the propagation of pulses in optical fibers are obtained. Special values are given in the model used, and two nonlinear differential equations are obtained. Nonlinear equations are reduced to ordinary differential equations with the help of wave transformations. Then, the obtained differential equations are solved by two different methods, namely the modified simplest equation and the modified Kudryashov procedures. The solutions are given by hyperbolic, trigonometric and rational functions and the results are useful for optics, engineering and other related areas. Finally three-dimensional, contour and two-dimensional shapes are given for some solutions. These figures are important for understanding the motion of the wave. The given methods are applied to the equations for the first time. To the best of the authors' knowledge, these results are new and have not been obtained in the literature. The results are useful for applied mathematics, physics and other related areas.
Arzu Akbulut, Waseem Razzaq, Filiz Tascan. Study for obtaining the exact solutions for the propagation of pulses in optical fibers[J]. Communications in Theoretical Physics, 2024, 76(10): 105003. DOI: 10.1088/1572-9494/ad526c
1. Introduction
In the literature, nonlinear partial differential equations (NLPDEs) have an important role; therefore there has been much research into obtaining exact solutions for NLPDEs. As a result of these studies, many new methods for finding exact solutions of NLPDEs have been discovered. For example, Wazwaz applied the tanh technique to the Dodd–Bullough–Mikhailov equation [1], Hashemi et al utilized the Nucci reduction method with the (2+1)-dimensional Korteweg–de Vries equation [2], Biswas et al handled the method of undetermined coefficients [3], Zhou et al studied the Riccati equation mapping approach [4], Biswas et al applied the trial equation procedure to the Radhakrishnan–Kundu–Lakshmanan equation with a couple of integration schemes [5], Karakoc et al discussed the modified regularized long wave equation [6], Raza et al studied the extended trial approach [7], Hosseini et al applied the unified method to the generalized third-order nonlinear Schrödinger equation [8], Zafar et al utilized the modified extended tanh expansion method in the Biswas–Arshed model with full nonlinear form [9], Kaplan et al applied the exponential rational function procedure to the wave packet envelope [10], Rabie et al applied the extended F-expansion technique to the generalized derivative nonlinear Schrödinger equation with quintic nonlinearity [11], Hashemi et al handled the Nucci reduction method [12], Kudryashov and Lavrova employed the simplest equation method with the Chavy–Waddy–Kolokolnikov model [13] and Niwas and Kumar used the inverse (G'/G)-expansion method [14]. Finally, Kumar and colleagues used a generalized exponential rational function method [15], similarity reduction [16], the extended Jacobian elliptic function expansion method [17] and Lie symmetry analysis [18, 19].
The propagation pulse in the equation describing the propagation of pulses in optical fibers serves as a cornerstone for understanding, predicting and optimizing various phenomena in optical communication, nonlinear optics and other related fields, driving advances in diverse technological areas. This equation is important for the understanding signal transmission and exploring nonlinear phenomena. The equation describing the propagation of pulses in optical fibers, often referred to as the nonlinear Schrödinger equation or its variants, holds immense importance in several fields, particularly optical communications and nonlinear optics. Therefore, we have concluded that it is crucial to work on this model. In this paper, the propagation pulse in an optical fiber is explained using the following new PDE:
where W(x, t) is the complex function representing an optical wave and α1, β1, γ1, δ1 and λ1 are the parameters. If we substitute m = 0 in equation (1) , the Schrödinger equation is obtained as follows:
The given models are handled in some studies in the literature, for example Kudryashov obtained the first integrals of the given models [20], Zayed and Alngar applied three different methods to obtain the optical soliton solutions [21] and Raza et al used the sine-Gordon expansion and modified auxiliary equation procedures [22].
In this paper we will obtain the exact solutions of equations (3) and (4) via the modified simplest equation and the modified Kudryashov procedures. For this purpose, a description of the procedures will be given in section 2. Then, the procedures will be applied in section 3. Some plots of the results will be given in section 4. Finally, conclusions will be given in section 5.
where bi (i = 1, 2, ..., m) are non-variables to be found and bm should be different from zero. m is a balance term and is calculated according to the principle of homogeneous balance. Let us determine the integer m by balancing the highest-order linear term with the highest-order nonlinear term in equation (7). The function σ(ϵ) fulfills some ODEs. In this research, we use the Riccati equations as the simplest equation
Substituting equations (8) into (7) with (9) and setting all the coefficients of function σi to zero, the system of algebraic equations in the terms bi, c1, k1 and w1 is obtained; then we solve these algebraic equations. Finally, the obtained results are substituted into the auxiliary solutions [23, 24].
2.2. The modified Kudryashov procedure
We transform equation (7) using the wave transformation below
where ${b}_{i}\left(\,{,}(,i=0,1,\,\ldots ,\,m\right)$ are constants that will be determined later, m is calculated by the homogeneous balance principle and the function $\sigma \left(\varepsilon \right)$ is given by
Without ignoring equation (23), equation (21) is substituted into equation (7) to produce a set of algebraic equations for bi, a, Γ, c1, k1 and w1. Eventually, after solving the resulting system, the exact solutions to equation (5) are determined [25, 26].
3. Application of the procedures
In this section, the given models will be reduced to an ODE by transformation (6). Then, two different procedures will be applied to the obtained reduced equations. Equations (3) and (4) can be reduced to the following nonlinear ODEs:
To obtain closed form solutions for equation (31), the constraint condition α1 = γ1 = δ1 = 0 is been imposed. Equation (31) becomes
By substituting equation (33) into (29) and using some mathematical operations, we obtain a nonlinear algebraic system. Solving this system for different values of b0, b1, λ1 and w1 results in Cases 1 and 2.
By substituting equation (33) into equation (32) and using some mathematical operations, we obtain a nonlinear algebraic system. Solving this system for different values of b0, b1, λ1 and ω1 results in Cases 1 and 2.
then substitute equation (74) into equation (29) and collect the coefficients of ${\sigma }^{i}\left(\varepsilon \right)$, we obtain the determining equation system. If we solve the obtained system, we get the following solution family (Cases 1 and 2).
If we substitute the obtained values of the coefficients from equations (75) into equation (3) without ignoring equations (26), (27) and (74), we get the following solution:
If we take the auxiliary solution for equation (32) as (74), then substitute equation (74) in (32), then collect the coefficients of ${\sigma }^{i}\left(\varepsilon \right),$ we obtain the determining equation system. If we solve the obtained system, we get the following solution family (Cases 1 and 2):
If we substitute the obtained values of the coefficients equations (79) in (4) without ignoring equations (26), (27) and (74), we get the following solution:
Plots of some of the results will be given in this section. First, we give the plots for equation (34) in figure 1. When we set the parameters ω1 = − 0.8, n = 2, k1 = 0.5, α1 = 0.2 and γ1 = 0.7 in equation (34), the three-dimensional (3D) and contour plots shown in figure 1 are obtained. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the two-dimensional (2D) plot in figure 1; this figure represents the kink shape solution.
Figure 1. Plots of equation (34), when ω1 = − 0.8, n = 2, k1 = 0.5, α1 = 0.2 and γ1 = 0.7.
We show the plots for equation (59) in figure 2. When we set the parameters ω1 = 0.8, n = 2, k1 = 0.5 and β1 = 0.2 in equation (59), the 3D and contour plots obtained are shown in figure 2. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the 2D plot in figure 2; this figure also represents the periodic shape solution.
Figure 2. Plots of equation (59), when ω1 = 0.8, n = 2, k1 = 0.5 and β1 = 0.2.
We give the plots for equation (76) in figure 3. When we set the parameters in equation (76) as a = 0.2, Γ = 0.1, n = 2, k1 = 0.8, γ1 = 1 and α1 = 0.3, the 3D and contour plots are obtained are shown in figure 3. When we set t = 0, t = 0.5 and t = 1, red, green and blue lines, respectively, are obtained in the 2D plot in figure 3; this figure represents the kink shape solution.
Figure 3. Plots of equation (76), when a = 0.2, Γ = 0.1, n = 2, k1 = 0.8, γ1 = 1 and α1 = 0.3.
We show the plots for equation (82) in figure 4. We set the parameters in equation (82) as a = 0.2, Γ = 0.1, n = 2, k1 = 0.8 and β1 = 1, and the 3D and contour plots obtained are shown in figure 4. When we set t = 0, t = 0 and t = 1, red, green and blue lines, respectively, are obtained in 2D plot shown in figure 4; this figure represents kink shape solution.
Figure 4. Plots of equation (82) when a = 0.2, Γ = 0.1, n = 2, k1 = 0.8 and β1 = 1.
5. Conclusion
In this study, a new form of propagation pulse in optical fibers is studied to obtain the exact solutions. We assumed different values for the arbitrary constants of the given model and obtained two different nonlinear differential equations. Then, we applied suitable wave transformations to the models. Thus, the given models were reduced to nonlinear ordinary differential equations. Then, we applied the two different procedures, the modified simplest equation and the modified Kudryashov procedures, to the reduced equations. The solutions are given by hyperbolic, trigonometric and rational functions. When we compared our methods with the other methods in the literature our methods gave simpler results: our solutions are expressed as different types of functions and our methods are easier to apply to equations. We concluded that the methods discussed in our study are effective and powerful. However, the methods have some disadvantages: while successful for specific types of equations, the modified Kudryashov method may not be applicable to all nonlinear equations. Its effectiveness depends heavily on the specific structure of the equation. We give the 3D, contour and 2D plots for some results. These figures are important for understanding the motion of the wave. The results can be used in various fields such as physics, engineering and optics. These methods can apply to the different equations or researchers could apply different methods to the given models.
WazwazA M2005 The tan h method: solitons and periodic solutions for the Dodd-Bullough-Mikhailov and the Tzitzeica-Dodd-Bullough equations Chaos Soliton Fract.25 55
HashemiM S, MirzazadehM, AhmadH2023 A reduction technique to solve the (2+1)-dimensional KdV equations with time local fractional derivatives Opt. Quant. Electron55 721
KarakocS B G, GeyikliT, BashanA2013 A numerical solution of the modified regularized Long Wave MRLW equation using quartic B splines Mathematics3 231
7
RazaN, AslamM R, RezazadehH2019 Analytical study of resonant optical solitons with variable coefficients in Kerr and non-Kerr law media Opt. Quant. Electron51 59
RabieW Bet al2023 Investigation of solitons and conservation laws in an inhomogeneous optical BER through a generalized derivative nonlinear Schrödinger equation with quintic nonlinearity Opt. Quant. Electron55 825
KudryashovN A, LavrovaS F2023 Painleve test, phase plane analysis and analytical solutions of the Chavy-Waddy-Kolokolnikov Model for the description of bacterial colonies Mathematics11 3203
NiwasM, KumarS2023 Multi-peakons, lumps, and other solitons solutions for the (2+1)-dimensional generalized Benjamin-Ono equation: an inverse (G'/G)-expansion method and real-world applications Nonlinear Dyn.111 22499
KumarSet al2021 Abundant different types of exact soliton solution to the (4+1)-dimensional Fokas and (2+1)-dimensional breaking soliton equations Commun. Theor. Phys.73 105007
KumarS, KumarD, WazwazA M2019 Group invariant solutions of (3+1)-dimensional generalized B-type Kadomstsev Petviashvili equation using optimal system of Lie subalgebra Phys. Scr.94 065204
KumarS, RaniS, MannN2022 Diverse analytical wave solutions and dynamical behaviors of the new (2+1)-dimensional Sakovich equation emerging in fluid dynamics Eur. Phys. J. Plus137 1226
KumarS, DhimanS K2024 Exploring cone-shaped solitons, breather, and lump-forms solutions using the lie symmetry method and unifed approach to a coupled breaking soliton model Phys. Scr.99 025243
KumarS, KumarD, KumarA2021 Lie symmetry analysis for obtaining the abundant exact solutions, optimal system and dynamics of solitons for a higher-dimensional Fokas equation Chaos Soliton Fract.142 110507
ZayedE M E, AlngarM E M2021 Optical soliton solutions for the generalized Kudryashov equation of propagation pulse in optical BER with power nonlinearities by three integration algorithms Math. Models Methods Appl. Sci.44 315
HosseiniK, MayeliP, AnsariR2017 Modified Kudryashov method for solving the conformable time-fractional Klein-Gordon equations with quadratic and cubic nonlinearities Optik130 737
AkbulutA2023 Obtaining the soliton type solutions of the conformable time-fractional complex Ginzburg-Landau Equation with Kerr law nonlinearity by using two kinds of Kudryashov methods J. Math.2023 4741219