Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 22 March 2023
Sec. Mathematical Physics
This article is part of the Research Topic Symmetry and Exact Solutions of Nonlinear Mathematical Physics Equations View all 20 articles

Lie symmetries and reductions via invariant solutions of general short pulse equation

Muhammad Mobeen Munir
Muhammad Mobeen Munir1*Hajra BashirHajra Bashir2Muhammad AtharMuhammad Athar2
  • 1Department of Mathematics, University of Punjab, Lahore, Pakistan
  • 2Department of Mathematics, Division of Science and Technology, University of Education, Lahore, Pakistan

Around 1880, Lie introduced an idea of invariance of the partial differential equation (PDE) under one-parameter Lie group of transformation to find the invariant, similarity, or auto-model solutions. Lie symmetry analysis (LSA) provides us an algorithm to search for point symmetries for solving related linear systems for infinitesimal generators. Actually, point symmetries lead us to one-parameter family of solutions from a known solution. LSA is a program that provides us the exact solutions for the non-linear differential equations (DEs) in analogy of the program designed by Galois for algebraic polynomial equations. In this paper, we have carried out the LSA for computing the similarity solutions (symmetries) of the non-linear short pulse equation (SPE) for the cases when h(u) = eu, k(u) = uxx, h(u)=eun, and k(u) = uxx. In addition, an optimal system of one-dimensional sub-algebra has been set up. The reductions and invariant solutions for the generalized SPE are calculated corresponding to this optimal system as well. Reductions reduce the non-linear PDE or system of PDEs into non-linear reduced ordered ODE or system of PDEs. This helps to solve these systems of PDEs to reduced form. Graphical behavior of the transformed points of the 1-parameter solution functions have drawn.

1 Introduction

Galois used the group theory to discuss the solvability of algebraic polynomial equations. Sophus Lie used the same idea foe differential equations and devised a comprehensive program now known as Lie symmetry analysis (LSA). In his attempt, he also developed the theory of Lie groups with broad applications in many areas of mathematics, physics, and in applied sciences [1, 2]. [3] have explained the procedure of symmetry reductions and exact solutions for the non-linear PDEs using three different methods that are direct, classical, and non-classical. [4] used LSA for systems of non-linear PDEs including the solutions, for system of non-linear coupled PDEs in real physical application, for the unsteady liquid and gas flow in long pipelines, for approximated long wave equations in shallow water and for the general dispersive long-wave equation.

Non-linear short pulse equation (SPE) represents the propagation of ultra-short pulses (light pulses) in optical fibers of silica [5]. Propagation of pulses in optical fibers was first depicted by the cubic non-linear Schrodinger equation (NLSE) which are used to provide the actual fiber connections and refer new systems of fiber communication to attain very high data transmission [6, 7]. Research studies on a large scale have been performed for the propagation of ultra-short pulses (very narrow pulses) that permit high quality fast data transmission [6, 8]. In case of short pulses (or ultra-short pulses), the rationality of NLSE lacks due to the breakdown, [9]. Also, the higher order terms that are involved in cubic NLSE cause difficulties, see Figure 1, for the behavior of NLSE as an output [10]. Therefore, determined the SPE which provides more accurate approximation of the solution of Maxwell’s equation in non-linear media rather than the NLSE [6]. The SPE has vast applications in many applied fields such as systems of impulse, mechanics, neural networks, and in many other fields of sciences. Determined the symmetries of SPE and travelling wave solution by parametric representation and power series process, respectively, [11]. Evaluated the symmetry reductions and conservation laws by using the direct method for SPE, [12]. Authors also determined the Lie symmetries for SPE through the non-local system. Established the results for well-posedness of solutions which are bounded for homogenous IBVP and Cauchy problem connected with SPE, [5]. Matsuno constructed multiple exact solutions by using the direct method for three novel PDEs related with generalizations of SPE that are integrable, [13]. He gave the parametric representation of multi-soliton solutions of generalized SPE. LSA has been used by many mathematicians to explore the results related to the exact solutions of non-linear PDEs which depict physical phenomena [14]. Discussed the class of non-linear PDEs having an arbitrary order [15]. Authors estimated the determining equations for non-classical symmetries by using compatibility of original equations with invariant surface conditions.

FIGURE 1
www.frontiersin.org

FIGURE 1. Pulse propagation in NL-dispersive optics.

In this article, we have discussed the LSA for one of the modified form of SPE and see graphical behavior of the functions depending upon 1-parameter (ϵ) Lie groups. The non-linear SPE is as follows:

uxt=αu+13βu3xx,(1)

where u(x, t) is the magnitude of electric field. α and β are the real parameters. Considering the SPE of the following form

uxt=αhu+13βk3u,(2)

where we let the general functions h(u) and k(u) as:

h(u) = eu and k(u) = uxx,

h(u)=eun for nN, (n > 1) and k(u) = uxx.

It is worth mentioning that the case h(u) = un and k(u) = uxx for Eq. 2 has been recently discussed in the article [16]. We will find Lie point symmetries corresponding to the aforementioned cases and the optimal system with reductions and see their graphical behavior corresponding to the Lie symmetries.

2 Results

In the present section, we give our main results with computations.

2.1 Lie symmetries of SPE for the case of h(u) = eu and k(u) = uxx

Eq. 2 becomes

uxt=αeu+13βuxx3,(3)

Consider the one parameter Lie group of point transformations for Eq. 3.

x*=x+ϵλx,t,u+Oϵ2,t*=t+ϵμx,t,u+Oϵ2,u*=u+ϵνx,t,u+Oϵ2,(4)

where ϵR is the group parameter.The group generator of (4) is defined in the following vector form as:

W=λx,t,ux+μx,t,ut+νx,t,uu,(5)

where λ, μ and ν are infinitesimal functions of group variables. The second prolongation of the infinitesimal generator along with coefficients has the following form:

Pr2W=W+νxux+νtut+νxxuxx+νxtuxt+νttutt,νxx=DxDxνλuxμut+λuxxx+μuxxt,νxt=DxDtνλuxμut+λuxtx+μuxtt.(6)

where Dx and Dt are the total derivatives.

Apply the second prolongation of the infinitesimal generator Eq. 5 onto Eq. 3. Then, in order to calculate symmetry of Eq. 3, we have the equation of the following form:

Pr2Wuxtαeu13βuxx3|uxt=αeu+13βuxx3=0.(7)

Solving this equation

ανeuβνxxu2xx+νxt|uxt=αeu+13βuxx3=0.(8)

Substitute the values of νxx, νxt and Eq. 3 which leads to an under-determined system of equations given as:

μxx=0,μxu=0,μuu=0,λuu=0,μu=0,λu=0,λtu=0,μx=0,νuu=0,νxu=0,νxx=0,νtu=0,λt=0,λxx=0,ανeu+νxt+ανuλxμteu=0,23νu+53λx13μt=0.(9)

The solution of the aforementioned determining equations gives the coefficient functions in the form of

λx=15c1x+c3,μt=c1t+c2,νx,t,u=65c1.(10)

c1,c2 and c3 are arbitrary constants. Thus, the Lie algebra of the infinitesimal symmetries for the case n = 1 is

W1=15xx+tt65u,W2=t,W3=x.(11)

Theorem 3.1 The set of these generators is closed under the one parameter Lie groups Hiϵ which are generated by infinitesimal generators Wi for i = 1, 2, and 3 are given in the following table. The entries give the transformed points exp(ϵWi)(x, t, u) = (x*, t*, u*).

H1ϵ:x,t,ue15ϵx,eϵt,u65ϵ,H2ϵ:x,t,ux,t+ϵ,u,H3ϵ:x,t,ux+ϵ,t,u.(12)

where ϵR is the group parameter.Theorem 3.2 If u=B(x,t) satisfies Eq. 3, then, u(i)(i = 1, 2, and 3) are solutions of Eq. 3:

u1=65ϵBe15ϵx,eϵt,u2=Bx,tϵ,u3=Bxϵ,t.(13)

where ui=Hiϵ.Bi(x,t), (i = 1,2, and 3), ϵ ≪ 1 is any positive number.

The Eq. 13 provides a class of solutions for Eq. 3 depending upon the parameter ϵ and general function Bi where (i = 1, 2, 3). The Figures 25 show the graphical view of the functions ui, (i = 1, 2, 3) where ui attained from Lie symmetry groups Wi. These graphs are formed by letting different general functions in place of Bi in Eq. 13. The graphs are constructed from the maple.

FIGURE 2
www.frontiersin.org

FIGURE 2. For u(1)=65ϵ[cos(e15ϵx+eϵt)] and ϵ =0.000005.

FIGURE 3
www.frontiersin.org

FIGURE 3. For u(1*)=65ϵ[cos(ϵ15ϵx)+sin(eϵt)] and ϵ = 0.000005.

FIGURE 4
www.frontiersin.org

FIGURE 4. For u(2) = x3+2(tϵ) and ϵ = 0.000005.

FIGURE 5
www.frontiersin.org

FIGURE 5. For u(3)=2t ln (xϵ).

For first equation in Eq. 13, letting the general trigonometric function in place of B(x,t)

u1=65ϵcose15ϵx+eϵt,(14)

along-with ϵ = 0.000005 and abscissa x = −5 to 5, ordinate t = −5 to 5.

u1=0.000006cose0.000001x+e0.000005t.(15)

Figure 2 shows the graphical behavior of Eq. 15.letting another general value of function B(x,t)=cos(ϵ15ϵx)+sin(eϵt). The function becomes

u1*=0.000006cose0.000001x+sine0.000005t.(16)

Figure 3 shows the graphical view of Eq. 16.

For second equation of Eq. 13, considering a general function B(x,t)=x3+2(tϵ) for the same values of ϵ = 0.000005 and aforementioned coordinates for Eq. 14.

u2=x3+2t0.000005,(17)

Figure 4 shows its graphical view.

For last equation of Eq. 13, we let a general logarithmic function B=2tln(xϵ) and for similar values of ϵ, x, and t coordinates.

u3=2tlnx0.000005.(18)

Its graph is in Figure 5.

2.2 Optimal system of subalgebras

In this part, we will find the optimal system of one dimensional Lie subalgebras for Eq. 3 by using the adjoint representation. The corresponding commutator table and the adjoint table are as follows:Commutator Table: Adjoint Table:

Let us take a generator

W=β1W1+β2W2+β3W3,(19)

Case No.1 For β1 ≠ 0, the generator turns to

W=W1+β2W2+β3W3.(20)

Applying Adjeβ2W2 on Y gives

W=W1+β3W3,(21)

furthermore, proceeding in the same way

W=Adje5β3W3W3=W1,(22)

which successively makes the coefficients β2 and β3 equal to 0and implies that WW1.Case No.2 Without loss of generality, here we take β1 = 0 and β2 = 1, the generator becomes

W=W2+β3W3,(23)

Now, act Adjeβ3W3 on the aforementioned W,

W=W2+β3W3,(24)

Subcase No.2.1 If β3 < 0, then

W=W2W3.(25)

Subcase No.2.2 If β3 > 0, then

W=W2+W3.(26)

Case No.3 For β1 = β2 = 0 and β3 = 1. Thus, in the meanwhile we have WW3.Case No.4 Let consider β1 = 0 = β3 and β2 ≠ 0. In this case, the generator is WW2.

The optimal system of one-dimensional subalgebras admitted by Eq. 3) is as follows:

W=W1,W2,W3,W2±W3.(27)

2.3 Reductions and invariant solutions

2.3.1 Reduction by W2

The invariants for corresponding characteristic equation are as follows:

x=a,u=b,(28)

where a and b are arbitrary constants.

The invariant solution can be written in the form of b = f(a), implies that

u=fa,(29)

substituting this value in Eq. 3), we obtain

βfa3+3αefa=0.(30)

The solution of this reduced equation for β = 1 is given in the form of solution set as.

fa163αebk1dbxk2=0,
6k1arctan3I3563313αefa13+1k12k1xk2=0,
6k1arctan3I3563313αefa131k12k1xk2=0.

2.3.2 Reduction by W3

The corresponding characteristic equation to this generator is as follows:

dx1=dt0=du0,(31)

this gives two invariants

t=a1,u=b1,(32)

where a1 and b1 are arbitrary constants. It implies

u=ft,(33)

putting this in Eq. 3, we obtain

αeu=0,(34)

which gives a trivial solution for u = f(x).

2.3.3 Reduction by W1

The characteristic equation is

5dxx=dtt=56du,(35)

solving this, we obtain corresponding invariants of the form

r=tx5,s=eux6,(36)

from this

u=lnx6ftx5,(37)

where we obtain

ux=1xfr5rfr+6fr,uxx=fr25r2fr+6fr+30rfr25r2f2rx2f2r,uxt=fr5rfr5fr+5rf2rx6f2r.(38)

substituting these derivatives into Eq. 3, we obtain

3f5r5rfr5rfr+15rf4rf2r3αf7r+βfr6fr+25r2fr+30rfr25r2f2r3=0.(39)

Thus, non-linear PDE (3) reduces to a non-linear ODE.

2.3.4 Reduction by W2 + W3

The invariants that we gain by solving characteristic equation are as follows:

a3=xt,u=b3,(40)

a3 and b3 are arbitrary constants. The invariant solution corresponding to them is u = f (a3). Inserting this solution into Eq. 3 will give us a non-linear ODE of the form

3fa3+βfa33+3αefa3=0.(41)

2.3.5 Reduction by W2 W3

The invariants corresponding to characteristic equation for this case are a4 = x + t and b4 = u. Furthermore, its invariant solution is given as u = f (a4). Therefore, the Eq. 3 will be converted into an ODE of the form

3fa4βfa433αefa4=0.(42)

2.4 Determining lie symmetry of SPE for the case h(u)=eun and k(u) = uxx (n > 1)

The equation becomes

uxt=αeun+13βuxx3.(43)

The one-parameter Lie group of transformations and the second prolongation with coefficients are given in Eqs 4, 6, respectively for Eq. 43. Let the generator be

Z=λx,t,ux+μx,t,ut+νx,t,uu,(44)

Therefore, we have

Pr2Zuxtαeun13βuxx3|uxt=αeun+13βuxx3=0,(45)

simplification gives the following equation:

αnνeunun1βνxxuxx2+νxt=0,(46)

which is solved for the values of νxx and νxt, will give us the equation involving derivatives of infinitesimals with respect to dependent and independent variables and also the derivatives of dependent variable w.r.to independent variables. Substituting Eq. 43 and comparing the values of coefficients on both sides gives an under-determined system of equations

μxx=0μxu=0,μuu=0,λuu=0,μu=0,λu=0,λtu=0,μx=0,νuu=0,νxu=0,νxx=0,νtu=0,λt=0,λxx=0,23νu+53λx13μt=0,αnνeunun1+νxt+ανuλxμteun=0.(47)

To solve this system, we consider ν as:

ν=Ltx+Mu+Nt,(48)

which satisfies the aforementioned equations and then by solving the aforementioned system, we obtain

λx=c2,μt=c1,νx,t,u=0.(49)

c1 and c2 are any arbitrary constants. The infinitesimal generators for the one-parameter of Lie groups of transformations admitted in Eq. 43) are given by

Z1=t,Z2=x.(50)

These symmetry generators give us the symmetry groups Qiϵ for i = 1, 2:

Q1ϵ=x,t+ϵ,u,Q2ϵ=x+ϵ,t,u.(51)

If u = R (x,t) is a solution of Eq. 43), then ui for i = 1, 2, and 3 and ϵ ≪ 1 also satisfies Eq. 43,

u1=Rx,tϵ,u2=Rxϵ,t.(52)

Commutator Table:also,

Adjoint Table.

Proposition 5.1: The generators Z1 = t and Z2 = x form a two-dimensional abelian Lie symmetry algebra.

2.5 Optimal system, reductions and invariant solutions

Considering a generator Z = b1Z1 + b2Z2. This generator will established a set of optimal system comprising of Lie algebra

Z=Z1,Z2,b1Z1+b2Z2(53)

where b1, and b2 are arbitrary constants. The reduction of PDE Eq. 39 by using the generator Z1 leads to an invariant solution u = f (c1). The reduced non-linear ODE will be

3αeun+βu3=0(54)

The reduction through Z2 generates a trivial case for Eq. 39.

3 Conclusion

In this paper, we have carried out the LSA for computing the similarity solutions (symmetries) of the non-linear SPE for the cases when h(u) = eu and k(u) = uxx and h(u)=eun and k(u) = uxx in SPE (2). In addition, an optimal system of one-dimensional subalgebra has been set up. The reductions and invariant solutions for the generalized SPE are calculated corresponding to this optimal system as well. The graphs are formed by the maple for the functions obtained from the transformed points of one-parameter Lie groups.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

Article has been concieved by MM, drafted by HB and MA. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Oliveri F. Lie symmetries of differential equations: Classical results and recent contributions. Symmetry (2010) 2(2):658–706. doi:10.3390/sym2020658

CrossRef Full Text | Google Scholar

2. Knight JC, Broeng J, Birks TA, Russell PSJ. Photonic band gap guidance in optical fibers. Science (1998) 282(5393):1476–8. doi:10.1126/science.282.5393.1476

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Clarkson PA, Ludlow DK, Priestley TJ. The classical, direct, and nonclassical methods for symmetry reductions of nonlinear partial differential equations. Methods Appl Anal (1997) 4(2):173–95. doi:10.4310/maa.1997.v4.n2.a7

CrossRef Full Text | Google Scholar

4. Liu H, Zhang L. Symmetry reductions and exact solutions to the systems of nonlinear partial differential equations. Physica Scripta (2018) 94(1):015202. doi:10.1088/1402-4896/aaeeff

CrossRef Full Text | Google Scholar

5. Coclite GM, di Ruvo L. Well-posedness results for the short pulse equation. Z für Angew Mathematik Physik (2015) 66(4):1529–57. doi:10.1007/s00033-014-0478-6

CrossRef Full Text | Google Scholar

6. Schäfer T, Wayne CE. Propagation of ultra-short optical pulses in cubic nonlinear media. Physica D: Nonlinear Phenomena (2004) 196(1-2):90–105. doi:10.1016/j.physd.2004.04.007

CrossRef Full Text | Google Scholar

7. Agrawal GP. Nonlinear fiber optics. In: Nonlinear science at the dawn of the 21st century. Berlin, Heidelberg: Springer (2000). p. 195–211.

Google Scholar

8. Karasawa N, Nakamura S, Nakagawa N, Shibata M, Morita R, Shigekawa H, et al. Comparison between theory and experiment of nonlinear propagation for a-few-cycle and ultrabroadband optical pulses in a fused-silica fiber. IEEE J Quan Electron (2001) 37(3):398–404. doi:10.1109/3.910449

CrossRef Full Text | Google Scholar

9. Moloney JV, Newell AC. Nonlinear optics. Physica D: Nonlinear Phenomena (1990) 44(1-2):1–37. doi:10.1016/0167-2789(90)90045-q

CrossRef Full Text | Google Scholar

10. Xiao Y, Maywar DN, Agrawal GP. New approach to pulse propagation in nonlinear dispersive optical media. JOSA B (2012) 29(10):2958–63. doi:10.1364/josab.29.002958

CrossRef Full Text | Google Scholar

11. Liu H, Li J. Lie symmetry analysis and exact solutions for the short pulse equation. Nonlinear Anal Theor Methods Appl (2009) 71(5-6):2126–33. doi:10.1016/j.na.2009.01.075

CrossRef Full Text | Google Scholar

12. Fakhar K, Wang G, Kara AH. Symmetry reductions and conservation laws of the short pulse equation. Optik (2016) 127(21):10201–7. doi:10.1016/j.ijleo.2016.08.013

CrossRef Full Text | Google Scholar

13. Matsuno Y. Parametric solutions of the generalized short pulse equations. J Phys A: Math Theor (2020) 53(10):105202. doi:10.1088/1751-8121/ab6f18

CrossRef Full Text | Google Scholar

14. Ali AT. New generalized Jacobi elliptic function rational expansion method. J Comput Appl Math (2011) 235(14):4117–27. doi:10.1016/j.cam.2011.03.002

CrossRef Full Text | Google Scholar

15. El-Sabbagh MF, Ali AT. New generalized Jacobi elliptic function expansion method. Commun Nonlinear Sci Numer Simulation (2008) 13(9):1758–66. doi:10.1016/j.cnsns.2007.04.014

CrossRef Full Text | Google Scholar

16. Zhao W, Munir MM, Bashir H, Ahmad D, Athar M. Lie symmetry analysis for generalized short pulse equation. Open Phys (2022) 20(1):1185–93. doi:10.1515/phys-2022-0212

CrossRef Full Text | Google Scholar

Keywords: short pulse equation, Lie point symmetry analysis, optimal system for lie subalgebras, reductions, invariant solutions

Citation: Mobeen Munir M, Bashir H and Athar M (2023) Lie symmetries and reductions via invariant solutions of general short pulse equation. Front. Phys. 11:1149019. doi: 10.3389/fphy.2023.1149019

Received: 20 January 2023; Accepted: 03 March 2023;
Published: 22 March 2023.

Edited by:

Gangwei Wang, Hebei University of Economics and Business, China

Reviewed by:

Xinyue Li, Shandong University of Science and Technology, China
Zenggui Wang, Liaocheng University, China

Copyright © 2023 Mobeen Munir, Bashir and Athar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Muhammad Mobeen Munir, mmunir.math@pu.edu.pk

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.