Abstract
This study's subject is a (3 + 1) dimensional new Hirota bilinear (NHB) equation that appears in the theory of shallow water waves. We investigate how particular dispersive waves behave in an NHB equation. In this regard, we first test the Painlavé integrability using the WTC-Kruskal method and second perform Lie point symmetry analysis on NHB. The algorithm outputs the Lie algebra, the symmetry reductions, and group invariant solutions. Using Lie point symmetries, NHB transforms into an ordinary differential equation. The integration architectures to solve this equation are the Bernoulli sub-ODE, 1/ G, and modified Kudryashov methods. We plot their graphs and observe how the solutions behaved in understanding the physical phenomenon. Additionally, we discuss the model utilizing nonlinear self-adjointness and Ibragimov's approach to generate conservation laws for each Lie symmetry generator.
1. Introduction
Materials science, ocean engineering, biology, solid state physics, chemistry, chemical physics, optical fibers, chemical kinetics, deep water waves, oceanography, signal processing, and all other nonlinear disciplines use nonlinear evolution equations (NLEE). A NLEE is a partial differential equation that explains how a physical system changes over time and under a set of given initial conditions. It affects physical processes such as system identification that appear beyond several scientific disciplines.
Energy is moved from one location to another through waves. The ocean's water level will rise and fall in intricate and unpredictable motions while transferring energy from a region of disturbance due to forces that include gravity and wind, as well as other factors like gravity or movement of the ocean floor. An ocean wave is the term used to describe this energy transfer via the ocean's surface water [Citation1]. The dynamic waves that characterize marine environments, examine sea movement, and display tropical tidal waves are the shallow water waves. These waves form as a result of wind activity on the sea and show long wavelength wonders where level length is more important than liquid depth. Sea foot influences these waves as they proliferate, aggravating the water's orbital movement. It might then result in unimaginable destruction to the coastal ecology, among other things [Citation2]. If the contact of the wave with the sea or ocean floor is having an impact on the waves at the surface, ocean waves are regarded as shallow water waves [Citation1].
As utilized in physical sciences, the exact solution for the most part alludes to an arrangement that captures the complete nature of a problem as contradicted to one that is approximate and perturbation solutions.
Setting up exact analytical solutions of NLEEs is imperative since it better portrays the physical phenomena. Exact analytical solutions are extremely beneficial for testing numerical algorithms used for the considered equation.
Group-invariant solutions, which is a special class of exact solutions, are the solution forms obtained by the Lie symmetry groups approach and remaining invariant through various geometric transformations of the considered model.
Conservation laws, which were closely related through Lie symmetry groups at the beginning of the twentieth century and which are quite well known in physical sciences, are the basic building blocks in the linearization of the equation, the use of numerical schemes and stability analysis, as well as the inference of the basic conservation laws of the model under consideration.
To obtain the exact solutions of NLEEs, various techniques, including Hirota bilinear theory [Citation3], Darboux transformation [Citation4], Painlavé analysis [Citation5], Lie symmetry analysis [Citation6], Bernoulli sub-ODE method [Citation7], -expansion [Citation8], the modified (G /G)-expansion method [Citation9], etc. were suggested in the literature. Some of the methods have even been used for fractional differential equations [Citation10, Citation11].
A self-reinforcing single wave within a wave container is referred to as a soliton if it continually maintains its form while moving at a fixed pace. Solitons are constructed by arranging a variety of dispersive NLEEs to represent physical structures [Citation12].
The (1+1)-dimensional Korteweg–de Vries (KdV) equation (1) (1) is a NLEE which is conventionally employed to model the behaviour of waves in shallow water.
Another NLEE that expresses the development of nonlinear, long waves with small amplitudes and gradual reliance on the transverse coordinate is the Kadomtsev–Petviashvili (KP) equation, which has two spatial and one temporal variables. The KP equation has two unique forms and may be expressed as follows in the normalized form: (2) (2) where may be a scalar function, x and y are separately the longitudinal and transverse spatial coordinates, and . The KPII equation is known under the case , and miniatures, for example, water waves with low surface pressure. The KPI equation, which may be used to demonstrate pulses in thin films with high surface pressure, is known under the case .
One of the such NLEEs (extended and high-dimensional version of the KdV equation) mentioned above lines is the following -dimensional new Hirota bilinear (NHB) equation [Citation13–15] (3) (3) Equation (Equation3(3) (3) ) is transformed into the KdV equation (4) (4) by using the transformations t = −P, x = −R, y = Q, z = Q and .
Equation (Equation4(4) (4) ) has N-soliton solutions via Hirota perturbation theory and is totally integrable. The literature on Equation (Equation3(3) (3) ) reveals a number of ongoing studies. Equation (Equation3(3) (3) ) was solved using the linear superposition approach in [Citation13], and two different resonant multiple wave solutions were discovered. Also, lump solutions to two different dimensional reductions with z = y and z = t in [Citation14]. Bilinear Bäcklund transformation is shown in [Citation16]. The N-solitary waves are further obtained in [Citation17] depending on the resulting bilinear form by applying Hirota's bilinear theory. In [Citation18], a periodic type II, rogue, bright, and dark wave evolution of Equation (Equation3(3) (3) ) with physical properties was developed.
The primary objective of this study is to carry out invariance analysis of the (3+1) dimensional NHB equation model, which emerged in shallow-water wave theory, through the technique of Lie-symmetry groups. In this way, variable changes corresponding to various geometric transformations and the formation of new reduced versions of the model and new exact (group invariant) solution forms will be investigated. However, the formation of local conservation vectors, which we cannot observe for this model, will also be demonstrated with the help of Lie symmetry groups and a coupled version of Noether's theorem developed by Ibragimov [Citation19].
This study consists of the following sections. In Section 2, we investigate the Painlavé integrability of Equation (Equation3(3) (3) ) using the WTC-Kruskal algorithm [Citation20–24]. In Sections 3 and 4, using Lie algebra [Citation25–30], we take into consideration the (3+1) dimensional Hirota equation. We investigate the Lie infinitesimals for Equation (Equation3(3) (3) ) using the fourth prolongation. We reduce Equation (Equation3(3) (3) ) to ODEs using some generators and their combinations. In Section 5 we have to reach at the solutions of the ODE's and, consequently, various solutions of Equation (Equation3(3) (3) ) using the Bernoulli Sub-Ode, and modified Kudryashov approaches mentioned in [Citation7]. In Section 6, as stated in [Citation31], we study and determine the quasi self-adjointness of Equation (Equation3(3) (3) ). Again in this section, we calculate the conservation vectors using Ibragimov's method [Citation19, Citation31]. In the Conclusion section we obtain 3D graphs of the solutions and analyse their behaviour over time, in order to better understand the physical phenomenon behind the model. In the last, we provide an overview of the study and conclusions.
2. Painlavé analysis
According to [Citation20], consider an NLEE, for example one dependent variable with respect to independent variables x, y, z, and t, (5) (5) where is the dependent variable, P is a polynomial about u and their derivatives. The Painlavé test [Citation20, Citation32, Citation33] is said to have been passed by Equation (Equation3(3) (3) ) if all solutions of Equation (Equation5(5) (5) ) are expressed as Laurent series, (6) (6) with enough arbitrary functions to be the order of (Equation5(5) (5) ), where u is analytic function and, . The WTC-Kruskal method consists of four phases [Citation20–23].
Phase 1. To find the leading order exponents α and the leading order coefficients , set , and put them into (Equation5(5) (5) ). The highest order derivative terms and the nonlinear terms can be balanced to find all possible . The process comes to a halt if the only α options are not integers.
Phase 2. Determine all potential truncated expansions of the form (7) (7) for each pair of (α, ) from Phase 1. can be determined by substituting (Equation7(7) (7) ) into Equation (Equation5(5) (5) ) and equating coefficients of same powers of ϕ. If cannot be determined, it is impossible to truncate the series (Equation6(6) (6) ) at constant terms.
Phase 3. Calculate all integer powers k, sometimes known as resonances, at which arbitrary functions may enter the series (Equation6(6) (6) ) if α is integer. is obtained by replacing the simplified version of (Equation6(6) (6) ) into Equation (Equation5(5) (5) ) and collecting the terms with the least ϕ powers. Given that Q is a matrix. Q's elements depends on k. The resonances of (Equation5(5) (5) ) are roots of . All resonances in a single PDE, whether positive or negative, should be distinct integers [Citation20, Citation33]. The method stops finding a non-integer root when trying to find k.
Phase 4. The coefficients at non-resonances are computed, and at each positive resonance, compatibility is tested, as the final phase. For each branch of (Equation5(5) (5) ), we insert the truncated expansions into Equation (Equation5(5) (5) ), where N is the largest resonance. If the compatibility criteria for one branch cannot be met, the next branch's compatibility conditions should be checked. If the above procedures can be carried out consistently and clearly, Equation (Equation5(5) (5) ) is considered to pass the Painlavé test. If the series (Equation6(6) (6) ) is self-consistent and single-valued, and the primary branch allows as many arbitrary functions as the order of (Equation5(5) (5) ), then (Equation5(5) (5) ) is said to have the Painlavé property. The singular manifold in the series (Equation6(6) (6) ) is defined in Phases 1 and 2 by using WTC [Citation20], , , whereas in the last two phases, Kruskal's gauge is used for the singular manifold, , to simplify the computations.
By referring to the above-mentioned steps, according to the leading order analysis, and (8) (8) The truncated expansion to Equation (Equation3(3) (3) ), (9) (9) an auto-Bäcklund transformation may be derived by assuming that is a given solution of Equation (Equation3(3) (3) ) and inserting (Equation9(9) (9) ) into Equation (Equation3(3) (3) ). Namely, since is a solution to Equations (Equation3(3) (3) ) and (Equation9(9) (9) ) must be as well.
Using this value in Equation (Equation6(6) (6) ), and separating the first term from the sum we get; (10) (10) Additionally, using Equations (Equation8(8) (8) ) and (Equation10(10) (10) ) in Equation (Equation3(3) (3) ), it has been possible to create the characteristic equation for resonances, which is then solved to get resonances for ,1,4, and 6. We know that, arbitrariness of singular manifold relates to the resonance at . Now we must calculate the coefficients and using the recursion relation and to satisfy the compatibility requirements for the existence of the free functions and . When necessary calculations are made, we get; The compatibility condition does not hold at resonance k=6 is where . As Equation (Equation3(3) (3) ) does not meet the conventional Painlavé property, we may deduce that it is not Painlavé integrable.
3. Lie symmetry analysis
Following similar procedures in [Citation30, Citation34–40] the Lie symmetry analysis of Equation (Equation3(3) (3) ) is studied in this section. First, we analyse the one-parameter Lie group of infinitesimal transformation; The Lie generator of Equation (Equation3(3) (3) ) is the following vector field (see, [Citation30, Citation34–40]). In this case, , , , , η are infinitesimal of the variables . , , , and η should meet the infinitesimal condition when the fourth prolongation of the infinitesimal generator is applied to Equation (Equation3(3) (3) ). Utilizing the prolongation and total derivative definitions to substitute from this condition, then matching the coefficients of the dependent variable and its derivatives, we encounter the following overdetermined system for infinitesimal functions: By solving it we get, where and are arbitrary constants and and are arbitrary functions. As a result, the infinitesimal symmetries of Equation (Equation3(3) (3) ) give rise to nine-dimensional Lie algebra, which is covered by the independent operators: (11) (11) To get group transformations, , we must solve the following initial value problems: (12) (12) Table is the commutator table of the Lie algebra for the NHB. It is clearly seen from the table that the structure whose bases are the generators in question is a Lie algebra. By virtue of (Equation12(12) (12) ) we build one parameter Lie groups of transformation; (13) (13) The right-hand elements show the transformed point .
If we assume that is a solution to Equation (Equation3(3) (3) ) and use the groups above, we write the new solutions as follows;
(14) (14)
4. The Lie symmetry reductions corresponding to some Lie vector fields
The Lie groups approach is an useful technique for the investigating the integrability aspects of the under study PDE (or systems of PDEs). A symmetry group of a system of differential equations is, in short, a group that transforms one solution of the system into another. In this section, symmetry reductions and group-invariant solutions corresponding to Lie point generators obtained in Section 3 will be systematically discussed.
For The characteristic equation will be By solving this equation we procure the following invariants: and thereby Equation (Equation3(3) (3) ) is reduced to (15) (15) This equation can be solved directly with Maple and the exact solution reads
(16) (16)
For . The characteristic equation is By solving this equation we deduce the following invariants: (17) (17) Hence, Equation (Equation3(3) (3) ) is rewritten in terms of the invariants of (Equation17(17) (17) ) as (18) (18) Equation (Equation18(18) (18) )'s exact solution can be constructed as shown below
For . The characteristic equation can be written as By solving the characteristic equation we have the following invariants: Then Equation (Equation3(3) (3) ) is reformulated in terms of the above invariants as (19) (19) By following the same procedure we reach to fourth-order ODE, (20) (20) Using various analytical approaches, we shall study the exact solutions to the resulting Equation (Equation20(20) (20) ).
5. Exact solutions through three distinct integration schemes
This section focuses on Equation (Equation20(20) (20) ). By through three well-automated integration architecture, we shall retrieve exact group-invariant solutions.
5.1. Bernoulli sub-ODE methods
Consider the nonlinear ordinary differential equation (NLODE) of the form [Citation7]: (21) (21) The following series expansion provides the exact solution of Equation (Equation21(21) (21) ): (22) (22) Here the constants () need to be computed, N is determined using the balancing principle and provides the following first-order ODE; (23) (23) where λ and μ are constants and (24) (24) An algebraic equation system is yielded by inserting Equation (Equation22(22) (22) ) into Equation (Equation21(21) (21) ) and accumulating the polynomial coefficients setting to zero. By using the Maple package program, the system is solved, and the results are then substituted in (Equation22(22) (22) ). Thus the exact solutions of Equation (Equation21(21) (21) ) are accomplished.
Now we use this procedure for Equation (Equation20(20) (20) ). We arrive at N = 1 using the balancing principle. The solution of (Equation20(20) (20) ) of the form (25) (25) Putting (Equation25(25) (25) ) into (Equation20(20) (20) ), we get the following system after collecting the coefficients of : (26) (26) The solutions of (Equation26(26) (26) ) are and Hence (27) (27) Going back sequentially, the exact group invariant solution of Equation (Equation3(3) (3) ) reads (28) (28)
5.2. method
Equation (Equation20(20) (20) ) has a solution of the type [Citation7] (29) (29) where satisfies the below second-order ODE. (30) (30) (31) (31) Substituting Equation (Equation29(29) (29) ) into Equation (Equation20(20) (20) ) and collecting the coefficients of we obtained the following algebraic system: (32) (32) Assisting via Maple we get; and . Thus, another exact solution of Equation(Equation20(20) (20) ) reads; (33) (33) and thereby the exact group-invariant solution of Equation (Equation3(3) (3) ) reads
(34) (34)
5.3. Modified Kudryashov method
The solution of (Equation20(20) (20) ) is of the form (35) (35) where satisfies the below equation: (36) (36) The exact solution of (Equation36(36) (36) ) is represented as follows: Inserting Equation (Equation35(35) (35) ) into Equation (Equation20(20) (20) ), we obtain the following algebraic system after gathering the coefficients of : Assisting via the Maple program we acquire, and
Hence, the exact solution of Equation (Equation20(20) (20) ) reads (37) (37) As a result, the non-travelling group invariant solution is given as (38) (38)
6. Conservation laws
This section shall serve two aims. In the first subsection we examine the quasi-self-adjointness of Equation (Equation3(3) (3) ) and subsequently will include local conservation laws via Ibragimov's nonlocal conservation theorem [Citation19, Citation31]. When we examine the studies on NHB, as far as we can see, we believe that conservation laws for NHB will make important contributions to the literature.
6.1. Quasi self adjoint
As detailed in [Citation31], we provide the fourth-order formal Lagrangian (39) (39) and adjoint equation (40) (40) for Equation (Equation3(3) (3) ). (41) (41) is obtained by substituting for in the adjoint equation and appears with it. Consequently, Equation (Equation3(3) (3) ) is not self-adjoint.
Substituting in adjoint equation (Equation40(40) (40) ) and using (42) (42) condition we get; When the system obtained through the matching of the derivatives of u 's coefficients is solved, v = c is found where c is a constant. As a result, Equation (Equation3(3) (3) ) is quasi-self-adjoint according to the definition given in [Citation31].
6.2. Local conservation laws via Ibragimov's nonlocal conservation method
Let us use the general theory for conservation laws [Citation19, Citation31] to construct the conserved vectors of Equation (Equation3(3) (3) ). The elements that make up the conservation vector are given as follows: (43) (43) where is the Lie characteristic function and is formal Lagrangian, η and in relation to the Lie point symmetry . The considered Equations (Equation3(3) (3) ) and (Equation43(43) (43) ) become
(44) (44) We have shown in the previous subsection that v = c (constant). We have taken v = 1 for convenience in calculations.
•Remark We calculated conservation laws for each generator, finding non-trivial conservation laws will be of importance to us. Since for each function with , conservation laws derived for generators and are trivial. In [Citation6], this is called the second type of triviality.
– | Case 1 For The associated Lie characteristic function is provided as Thus the conserved vector via Equation (Equation44(44) (44) ) has the following components: | ||||
– | Case 2 For = The associated Lie characteristic function is provided as Thus the conserved vector via Equation (Equation44(44) (44) ) has the following components: | ||||
– | Case 3 For The associated Lie characteristic function is provided as Thus the conserved vector via Equation (Equation44(44) (44) ) has the following components: | ||||
– | Case 4 For The associated Lie characteristic function is provided as Thus the conserved vector via Equation (Equation44(44) (44) ) has the following components: | ||||
– | Case 5 For The associated Lie characteristic function is provided as Thus the conserved vector has via Equation (Equation44(44) (44) ) the following components: | ||||
– | Case 6 For The associated Lie characteristic function is provided as Thus the conserved vector via Equation (Equation44(44) (44) ) has the following components: |
7. Conclusion
In Figure we examined how the kink wave profile for changed over time. At t = 2, the waves get smaller and then become stable t = 7.
In Figure as time passes, the annihilation of a multi-soliton for has been observed. After t = 7, the multisoliton shape is converted into a stable wave profile.
In Figure we observed the stabilization of the king-bright profile over time.
In Figure another kink wave profile was seen in a similar way and was stable at t = 2.5.
In this work, using the WTC Kruskal algorithm, we examined the Painlavé property of the (3+1) dimensional NHB equation. We have demonstrated that the compatibility condition does not hold at the k = 6 resonance and is not Painlevé integrable. We obtained 9-dimensional Lie symmetry algebra for Equation (Equation3(3) (3) ). We reduced the NHB equation to the ODE using some generators and their combination.Then we obtained the exact group invariant solutions of the NHB equation using Bernoulli sub-ODE, and modified Kudryashov methods. Lie symmetry analysis of many important equations has been studied in the literature, but as a result of our research, we have seen that the symmetry analysis of the NHB has not been studied before. For this reason, we believe that our results are new, useful and different from the solutions available in the literature. On the other hand, we have provided that Equation (Equation3(3) (3) ) is quasi-self- adjoint and for each Lie point generator we get local conservation laws via Ibragimov's nonlocal conservation method. In future works, we plan to discuss non-classical symmetries and the relationships between symmetries and conservation laws.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Correction Statement
This article has been corrected with minor changes. These changes do not impact the academic content of the article.
References
- https://study.com/academy/lesson/shallow-water-waves-definition-speed-calculation.html.
- Tanwar DV, Kumar M. On Lie symmetries and invariant solutions of Broer–Kaup–Kupershmidt equation in shallow water of uniform depth. J Ocean Eng Sci. 2022. DOI:10.1016/j.joes.2022.04.027
- Ma WX. Soliton solutions by means of Hirota bilinear forms. Partial Differ Equ Appl Math. 2022;5:100220.
- Guo B, Ling L, Liu QP. Nonlinear Schrödinger equation: generalized Darboux transformation and rogue wave solutions. Phys Rev E. 2012;85(2):026607.
- Musette M. Painlevé analysis for nonlinear partial differential equations. In: Conte R, editor. The Painlevé property, one century later, CRM Series in Mathematical Physics. New York (NY): Springer; 1999. p. 517–572.
- Olver PJ. Applications of Lie groups to differential equations. Vol. 107. New York: Springer Science & Business Media; 2000.
- Tariq KU, Bekir A, Zubair M. On some new travelling wave structures to the (3+1)-dimensional Boiti-Leon-Manna-Pempinelli model. J Ocean Eng Sci. 2022. DOI:10.1016/j.joes.2022.03.015
- Alam MN, Tunç C. New solitary wave structures to the (2+1)-dimensional KD and KP equations with spatio-temporal dispersion. J King Saud Univ Sci. 2020;32(8):3400–3409.
- Islam S, Alam M, Al-Asad M, et al. An analytical technique for solving new computational? Solutions of the modified Zakharov-Kuznetsov Equation arising in electrical engineering. J Appl Comput Mech. 2021;7(2):715–726.
- Alam MN, Tunç C. The new solitary wave structures for the (2+1)-dimensional time-fractional Schrodinger equation and the space-time nonlinear conformable fractional Bogoyavlenskii equations. Alexandria Eng J. 2020;59(4):2221–2232.
- Alam MN, Tunç C. Constructions of the optical solitons and other solitons to the conformable fractional Zakharov–Kuznetsov equation with power law nonlinearity. J Taibah Univ Sci. 2020;14(1):94–100.
- Ali A, Seadawy AR, Lu D. New solitary wave solutions of some nonlinear models and their applications. Adv Differ Equ. 2018;2018(1):1–12.
- Gao LN, Zhao XY, Zi YY, et al. Resonant behavior of multiple wave solutions to a Hirota bilinear equation. Comput Math Appl. 2016;72(5):1225–1229.
- Lü X, Ma WX. Study of lump dynamics based on a dimensionally reduced Hirota bilinear equation. Nonlinear Dyn. 2016;85(2):1217–1222.
- Wang C. Lump solution and integrability for the associated Hirota bilinear equation. Nonlinear Dyn. 2017;87(4):2635–2642.
- Gao LN, Zi YY, Yin YH, et al. Bäcklund transformation, multiple wave solutions and lump solutions to a (3+1)-dimensional nonlinear evolution equation. Nonlinear Dyn. 2017;89(3):2233–2240.
- Dong MJ, Tian SF, Yan XW, et al. Solitary waves, homoclinic breather waves and rogue waves of the (3+1)-dimensional Hirota bilinear equation. Comput Math Appl. 2018;75(3):957–964.
- Zeynel M, Yasar E. A new (3+1) dimensional Hirota bilinear equation: periodic, rogue, bright and dark wave solutions by bilinear neural network method. J Ocean Eng Sci. 2022. DOI:10.1016/j.joes.2022.04.017
- Ibragimov NH. A new conservation theorem. J Math Anal Appl. 2007;333(1):311–328.
- Xu GQ, Li ZB. Symbolic computation of the Painlevé test for nonlinear partial differential equations using Maple. Comput Phys Commun. 2004;161(1-2):65–75.
- Weiss J, Tabor M, Carnevale G. The Painlevé property for partial differential equations. J Math Phys. 1983;24(3):522–526.
- Hereman W, Göktaş Ü, Colagrosso MD, et al. Algorithmic integrability tests for nonlinear differential and lattice equations. Comput Phys Commun. 1998;115(2-3):428–446.
- Lou SY, Chen CL, Tang XY. (2+ 1)-dimensional (M+ N)-component AKNS system: Painlevé integrability, infinitely many symmetries, similarity reductions and exact solutions. J Math Phys. 2002;43(8):4078–4109.
- Zhang SL, Wu B, Lou SY. Painlevé analysis and special solutions of generalized Broer–Kaup equations. Phys Lett A. 2002;300(1):40–48.
- Kumar S, Ma WX, Kumar A. Lie symmetries, optimal system and group-invariant solutions of the (3+1)-dimensional generalized KP equation. Chin J Phys. 2021;69:1–23.
- Wang G, Vega-Guzman J, Biswas A, et al. (2+1)-dimensional Boiti–Leon–Pempinelli equation–domain walls, invariance properties and conservation laws. Phys Lett A. 2020;384(10):126255.
- Jadaun V, Kumar S. Symmetry analysis and invariant solutions of (3+1)-dimensional Kadomtsev–Petviashvili equation. Int J Geom Methods Mod Phys. 2018;15(08):1850125.
- Velan MS, Lakshmanan M. Lie symmetries and invariant solutions of the shallow-water equation. Int J Non Linear Mech. 1996;31(3):339–344.
- Sadat R, Kassem M. Explicit solutions for the (2+1)-dimensional Jaulent–Miodek equation using the integrating factors method in an unbounded domain. Math Comput Appl. 2018;23(1):15.
- Ali MR, Sadat R. Lie symmetry analysis, new group invariant for the (3+1)-dimensional and variable coefficients for liquids with gas bubbles models. Chin J Phys. 2021;71:539–547.
- Ibragimov NH. Nonlinear self-adjointness and conservation laws. J Phys A: Math Theor. 2011;44(43):432002.
- Steeb WH, Euler N. Nonlinear evolution equations and Painlevé test. Singapore: World Scientific; 1988.
- Roy-Chowdhury AK. Painlevé analysis and its applications. Vol. 105. Boca Raton (USA): CRC Press; 1999.
- Adem AR, Khalique CM, Biswas A. Solutions of Kadomtsev–Petviashvili equation with power law nonlinearity in 1+3 dimensions. Math Methods Appl Sci. 2011;34(5):532–543.
- Yang H, Liu W, Yang B, et al. Lie symmetry analysis and exact explicit solutions of three-dimensional Kudryashov–Sinelshchikov equation. Commun Nonlinear Sci Numer Simul. 2015;27(1–3):271–280.
- Yong X, Chen Y, Huang Y, et al. Lie symmetry analysis for a generalized Conde-Gordoa-Pickering equation via equivalence transformations. Chin J Phys. 2020;66:430–435.
- Liu YK, Li B. Nonlocal symmetry and exact solutions of the (2+1)-dimensional Gardner equation. Chin J Phys. 2016;54(5):718–723.
- Ali MR, Ma WX. New exact solutions of Bratu Gelfand model in two dimensions using Lie symmetry analysis. Chin J Phys. 2020;65:198–206.
- Kumar S, Kumar D, Wazwaz AM. Group invariant solutions of (3+1)-dimensional generalized B-type Kadomstsev Petviashvili equation using optimal system of Lie subalgebra. Phys Scr. 2019;94(6):065204.
- Kumar S, Nisar KS, Kumar A. A (2+1)-dimensional generalized Hirota–Satsuma–Ito equations: lie symmetry analysis, invariant solutions and dynamics of soliton solutions. Res Phys. 2021;28:104621.