Home Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering
Article Open Access

Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering

  • Oke Davies Adeyemo and Chaudry Masood Khalique EMAIL logo
Published/Copyright: March 26, 2024

Abstract

This research article analytically investigates a soliton equation of high dimensions, particularly with applications, and precisely in the fields of physical sciences and engineering. The soliton equation of high dimensions, particularly with applications, and precisely in the fields of physical sciences along with engineering, is examined with a view to securing various pertinent results of interest. For the first time, the conserved currents of an integrodifferential equation (especially those of higher dimensions) are calculated using a detailed optimal system of one-dimensional subalgebras. Infinitesimal generators of diverse structures ascribed to Lie point symmetries of the understudy model are first calculated via Lie group analysis technique. Additionally, we construct various commutations along Lie-adjoint representation tables connected to the nine-dimensional Lie algebra achieved. Further to that, detailed and comprehensive computation of the optimal system of one-dimensional subalgebras linked to the algebra is also unveiled for the under-investigated model. This, in consequence, engenders the calculation of abundant conserved currents for the soliton equation through Ibragimov’s conserved vector theorem by utilizing its formal Lagrangian. Later, the applications of our results are highlighted.

1 Introduction

Nonlinear equations with dispersive property picture a class of mathematical equations refereed usually to as partial differential equations (PDEs) [122]. These PDEs are key in delineating a number of physical models inclusive of waves residing in a shallow water channel, confinement of Bose–Einstein condensate, light propagation in an optical waveguide, and so forth.

Furthermore, it has been observed that these equations simply put are non-solvable explicitly, yet mathematical techniques with analysis, together with computational capabilities, are divulging the means of engendering these models as prognostic tools with a high level of efficiency. This provides a rich phenomenon, for instance, balancing nonlinearity alongside dispersion produces coherent structures like vortices, and solitons – metastable states that are long-lived, found out to be waves localized and propagates with little or no disfigurement. The applicability of such structures can be seen in optical communication, within which solitons are engaged in conveying information. Additionally, these structures also have physical interesting features as a result of their particle-like actions.

In the recent times, investigations have largely been turned to nonlinear partial differential equations (NLNPDEs) as well as exact travelling wave results associated with these NLNPDEs. As a result, diverse complex physical happenstances are depicted via these NLNPDEs. A few of these NLNPDEs including the Boussinesq–Burgers-type system recounting shallow water waves and also emerging near ocean beaches and lakes were given attention in this article [1]. Moreover, Adeyemo et al. [2] examined another generalized NLNPDEs called advection–diffusion equation with power law nonlinearity in fluid mechanics. This generalized equation characterized buoyancy-propelled plume movement embedded in a medium that is bent on nature. Additionally, the vector bright solitons, alongside their various interaction attributes related to the coupled Fokas–Lenells system [3], was studied in the given reference. The femtosecond optical pulses embedded in a double-refractive optical fiber, modeled into an NLNPDEs, were further investigated. Recently, Adeyemo et al. [4] examined a (3+1)-D nonlinear generalized type of potential Yu–Toda–Sasa–Fukuyama model existent in Physics alongside Engineering. Besides, Du et al. [5] investigated the modified as well as generalized Zakharov–Kuznetsov model, delineating the ion-acoustic meandering solitary waves resident in a magneto-plasma and possessive of electron–positron–ion observable in the autochthonous universe. This model was utilized in representing dust-magneto-acoustic, and ion-acoustic, together with dust-ion-acoustic waves in the laboratory dusty plasmas. Further to that, a generalized structure of the Korteweg–de Vries–Zakharov–Kuznetsov model was investigated by Khalique and Adeyemo [6]. The dilution of warm isentropic fluid alongside cold static framework species together with hot isothermal, applicable in fluid dynamics was recounted via the use of the model. The list continues unending, see more in previous studies [4,715].

Sophus Lie (1842–1899) with his quintessential work on Lie Algebras [1720] which is essentially a unified approach for the treatment of a wide class of differential equations (DEs). With the inspiration of Galois theory, Sophus Lie, a Norwegian mathematician, established symmetry methods and demonstrated that many of the known ad hoc methods of integration of DEs could be obtained in a systematic manner. The approach has evolved into a helpful tool for solving DEs, classifying them, and preserving the solution set.

Furthermore, it has been observed that conservation laws are established and entrenched natural laws that have been studied by many researchers in various scientific fields. Conservation laws that are commonly used in this context include conservation of linear momentum in an isolated system, conservation of electric charge, conservation of energy, conservation of mechanical energy in the absence of dissipative forces, and many others. Conservation laws are deliberated to be basic laws of nature, with extensive application in physics and numerous other fields. Some of the important criteria of conservation laws are as follows [21]:

  1. the stability analysis and the global behavior of solutions.

  2. the development of numerical methods and provide an essential starting point for finding potential variables and nonlocally related systems.

  3. the investigation of integrability and linearization mappings.

Securing soliton solutions to NLNPDEs, as a result of its pertinence, is thus becoming a crucial point of interest and active space of investigation to scientists. Consequently, in a bid to gain the soliton solutions, travelling wave solutions, and other interesting exact solutions to NLNPDEs, sturdy approaches have been developed in the literature by scientists. We have some of them as power series solution method [22], simplest equation method [23], Darboux transformation [24], multiple exponential function method [25], just to mention a few. Others include bifurcation technique [26], Painlevé expansion [27], homotopy perturbation technique [28], tanh–coth approach [29], extended homoclinic test approach [30], Cole–Hopf transformation technique [31], Adomian decomposition approach [32], Bäcklund transformation [33], F-expansion technique [34], rational expansion technique [35], extended simplest equation approach [36], Kudryashov’s technique [37], Hirota technique [38], Darboux transformation [39], tanh-function technique [40], G G -expansion technique [41], sine-Gordon equation expansion technique [42], generalized unified technique [43], exponential function technique [44], and so on. Since the inception of Kadomtsev and Petviashvili’s hierarchy of equations a little more than half a century ago, dozens of research papers have emerged, each exploring an aspect of this rich domain of equations, see for example previous studies [4551].

There have emerged NLNPDEs that have been solved by using already arisen mathematical techniques. Nevertheless, not all the emergent modeled equations are solvable. Examination of analytic explicit outcomes to soliton equations is already in the limelight and as such highly significant with influence in physics of mathematics among others [911]. Such soliton model includes the 3D soliton-modeled Jimbo–Miwa-type [52] given as

(1.1) 2 q y t + 3 q x x v y + 3 q x q y x 3 q z z + q x x x y = 0 ,

investigated under the variable-dependent-Cole–Hopf transformation that reads

u = 2 ( log e τ ) x , q = 2 ( log e ω ) x .

Additionally, mapping (1.1) and Hirota bilinear relation to each other results in

(1.2) ( 2 D t D y + D x 3 D y 3 D z 2 ) ω ω = 0 .

Thus, bilinear differential operators unveiled in (1.2), that is, D x , D y : Ω × Ω Ω are computed as

(1.3) D y n D x m [ h ( y , x ) s ( y , x ) ] = y y n x x m × [ h ( y , x ) s ( y , x ) ] y = y , x = x ,

where m , n 0 , h , s Ω : C C as well as C standing for differential functions space. Thus, we give Eq. (1.2) regarding ω as

(1.4) ω ω x x x y + 2 ω ω t y 3 ω ω z z 3 ω x x y ω x + 3 ω x x ω x y 2 ω t ω y ω x x x ω y + 3 ω z 2 = 0 .

Furthermore, Asaad [52] engaged the Pfaffian technique to achieve closed-form results to the soliton equation (1.4).

Next, the high dimensions soliton model divulged as

(1.5) 3 Ψ x z ( 2 Ψ t + Ψ x x x 2 Ψ Ψ x ) y + 2 ( Ψ x x 1 Ψ y ) x = 0 ,

where Ψ = Ψ ( t , x , y , z ) , was at first gleaned by Geng [53]. The author went ahead to decompose (1.5) via the infuse of bi-dimensional Ablowitz–Kaup–Newell–Segur equations into solvable ordinary differential equations system. Liu and Zhang [54] examined (1.5) and thus gained mixed lump strip results together with solitonic-lump outcomes satisfying the equation through popular Hirota bilinear approach. Besides, investigations on them reveal that solitonic lump outcomes localized rationally in every observed direction within the space. Not only that diverse copious solutions that are periodic to the model (1.5) were secured under three-wave solution approach alongside the Hirota bilinear technique, see also the study of Liu et al. [55]. Geng and Ma [56] engaged a non-linearized-Lax-pair-approach in fetching geometric algebraic results of (1.5) explicitly and in Riemann theta function structure. Moreover, the Wronskian approach alongside Hirota technique was utilized to gain N-solitonic-solutions of the equation by the authors. Based on Pfaffian derivative relation, Jian-Ping and Xian-Guo [57] found Grampian determinant solutions to (1.5). Meanwhile, on invoking the bilinear Bäcklund transformation, more viable results of the equation have been explicitly gained by Jian-Ping [58]. Additionally, Wang et al. [59] broke down (1.5) into three integrable bi-dimensional models. They achieved this on the basis of a quartet Lax pair condition. These three models are as follows: nonlinear Schrödinger model, Lakshmanan–Porsezian–Daniel model, and complex modified Korteweg–de Vries equation in dimensions that are dissimilar. General rational Nth-order outcome given in a compact structure way was achieved for (1.5) using the Darboux transformation technique together with the limit approach. Not only that, in the study of Wang and Wei [60], the decomposition technique was invoked to secure solution regarding N -anti-dark soliton of the equation under contemplation on a finite background via the application of the limit technique and the Darboux transformation. Meanwhile, for the N -anti-dark soliton solution, the asymptotic analysis was furnished. Additionally, the authors proved that elastic-nature collision existed between multiple antidark solitonic outcome.

Hence, in our study, we seek to explicitly examine the optimal solutions assented to the three-dimensional soliton equation (1.5) via robust Lie group theoretic technique. Consequently, the integral function emergent in (1.5) is first eliminated via the representation v = u y d x . Thus, (1.5) alters the equation system that reads

(1.6a) 3 u x z 2 u t y + 2 u x u y + 2 u u x y + 2 v u x x + 2 v x u x u x x x y = 0 ,

(1.6b) u y v x = 0 ,

denoted as HD-SOLeqn for short. We observed that recently Khalique and Adeyemo [61,62] bought into play Lie-symmetric approach to gain various abundant invariant solutions of system (1.6). In addition, copious solutions in terms of closed-formed travelling waves of the under-study model via the systematic polynomial complete discriminant alongside elementary integral approaches. Meanwhile, homotopy formula was employed in computing some conserved quantities of (1.6).

Nonetheless, this study invokes the optimal Lie algebraic systematic approach to compute various vectors via a nine-dimensional Lie subalgebras ascribed to (1.6) to gain more extensive conserved quantities of the system with various applications in sciences and engineering. We state for the purpose of emphasis that this research engage a detailed computation of one-dimensional optimal system of Lie subalgebras, obtained from a nine-dimensional Lie algebra, to generate abundant conservation laws to (1.6). Moreover, for the first time, the significance of the associated conserved quantities are highlighted within the fields of physical sciences. All these attest to the fact that the work is novel and original.

Now, we catalog the rest of the research article as follows. Section 2 supplies well-thought-out steps adhered-to in calculating the Lie symmetries ascribed to the soliton equation. Moreover, one parametric transformation groups alongside a Lie-sub-algebraic optimal system is computed for the gained Lie algebra. In addition, Section 3 achieves conserved current calculations for optimal system of Lie subalgebras of solitonic system (1.6) in conjunction with the formal Lagrangian using Ibragimov’s conserved quantities theorem after which conclusions are furnished.

2 Lie algebra and optimal system of the HD-SOLeqn (1.6)

Computations of Lie-point symmetries of HD-SOLeqn (1.6) are first done from where optimal systems of sub-algebras are contrived. In consequence, engagement of the gained symmetries is taken into consideration to attain diverse possible conserved quantities ascribed to system (1.6).

2.1 Conspectus of infinitesimal generators of (1.6)

We suppose first that Lie transformation group of infinitesimal generators be explicated in the following format:

(2.1) t ¯ t + ε ξ 1 ( t , x , y , z , u , v ) + O ( ε 2 ) , z ¯ z + ε ξ 4 ( t , x , y , z , u , v ) + O ( ε 2 ) , x ¯ x + ε ξ 2 ( t , x , y , z , u , v ) + O ( ε 2 ) , u ¯ u + ε ϕ 1 ( t , x , y , z , u , v ) + O ( ε 2 ) , y ¯ y + ε ξ 3 ( t , x , y , z , u , v ) + O ( ε 2 ) , v ¯ v + ε ϕ 2 ( t , x , y , z , u , v ) + O ( ε 2 ) .

We now contemplate an infinitesimal dimensional Lie algebra covered by vector fields

(2.2) w = ξ 1 x + ξ 2 y + ξ 3 z + ξ 4 t + ϕ 1 u + ϕ 2 v ,

where coefficient functions ξ 1 , ξ 2 , ξ 3 , ξ 4 , ϕ 1 , ϕ 2 , depending on t , x , y , z , u , and v . Hence, one examines an appurtenant theorem:

Theorem 2.1

Suppose vector w is assumed to be the infinitesimal generators ascribed to classical point symmetric group of HD-SOLeqn (1.6), where ξ i , i = 1 , 2 , 3 , 4 alongside ϕ 1 and ϕ 2 , regarded as smooth functions of variables t , x , y , z , u , and v . Thus, one engenders results which are formatted as:

(2.3) ξ 1 = c 1 + 3 2 ( c 2 + c 4 ) t , ξ 2 = 1 2 ( c 2 + c 4 ) x + H 1 ( z , t ) + Q ( t ) , ξ 3 = c 2 y + H ( z ) , ξ 4 = c 3 + c 4 z , ϕ 1 = ( c 2 c 4 ) u + 3 2 H ( z ) Q ( t ) H t 1 ( z , t ) , ϕ 2 = 1 2 v ( c 2 + c 4 ) + 3 2 H z 1 ( z , t ) ,

where constants c 1 , c 2 , c 3 , c 4 are arbitrary along with functions H ( z ) , H 1 ( z , t ) , and Q ( t ) being observed to be, respectively, depending on their respective arguments.

See the comprehensive proof of Theorem 2.1 in the study of Khalique and Adeyemo [61].

Suppose one assumes in gained solution (2.3) that arbitrary H 1 ( z , t ) = c 5 z + c 6 t + c 7 , H ( z ) = c 8 z + c 9 , as well as Q ( t ) = 0 , we institute more symmetries of (1.6). Hence, the arrival at the subsequent corollary:

Corollary 2.1

Lie algebra represented by g for infinitesimal symmetries ascribed to soliton system (1.6) given in three dimensions consists of vector fields alongside their respective names, viz,

(2.4) w 1 = t , t i m e t r a n s l a t i o n , w 2 = 3 t t + x x 2 y y 2 u u + v v , S c a l i n g , w 3 = z , z s p a c e t r a n s l a t i o n , w 4 = 3 t t + x x + 2 z z 2 u u v v , S c a l i n g , w 5 = 2 z x + 3 v , G a l i l e a n b o o s t , w 6 = t x u , G a l i l e a n b o o s t , w 7 = x , x s p a c e t r a n s l a t i o n , w 8 = 2 z y + 3 u , G a l i l e a n b o o s t , w 9 = y . y s p a c e t r a n s l a t i o n .

Thus, HD-SOLeqn system (1.6) admits Lie algebra of nine dimensions whose basis are formatted as { w 1 , w 2 , w 3 , , w 9 } .

One observes quickly here that the instituted infinitesimal generators can be explicated as a linearly combined vectors formatted as

(2.5) w = a 1 w 1 + a 2 w 2 + a 3 w 3 + a 4 w 4 + a 5 w 5 + a 6 w 6 + a 7 w 7 + a 8 w 8 + a 9 w 9 .

The next phase of the research engenders the Lie transformation-groups connected to Lie-algebra g .

2.2 One-parameter Lie group transformation of (2.4)

On involving the Lie equations in previous studies [17,63] alongside the related initial conditions in the calculation of one parameteric transformation group ascribed to the gained generators (2.4). Thus, we institute a theorem given shortly:

Theorem 2.2

Suppose that transformation group G ε i ( t , x , y , z , u , v ) , i = 1 , 2 , 3 , , 9 of one parameter is engendered by infinitesimal generators w 1 , w 2 , w 3 , w 9 in (2.4), thus, for each generator, we secure accordingly

G ε 1 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t + ε , x , y , z , u , v ) , G ε 2 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t e 3 ε , x e ε , y e 2 ε , z , u e 2 ε , v e ε ) , G ε 3 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x , y , z + ε , u , v ) , G ε 4 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t e 3 ε , x e ε , y , z e 2 ε , u e 2 ε , v e ε ) , G ε 5 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x + 2 ε z , y , z , u , v + 3 ε ) , G ε 6 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x + ε t , y , z , u ε , v ) , G ε 7 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x + ε , y , z , u , v ) , G ε 8 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x , y + 2 ε z , z , u + 3 ε , v ) , G ε 9 : ( t ¯ , x ¯ , y ¯ , z ¯ , u ¯ , v ¯ ) ( t , x , y + ε , z , u , v ) ,

with ε R admired as the parametric group.

We direct the reader to the designated references with a view to gaining a much better under-standing of the ascribed proof of Theorem 2.2. Now, by extension, the subsequent theorem suffices, that is:

Theorem 2.3

If u = h ( t , x , y , z ) alongside v = g ( t , x , y , z ) fulfills HD-SOLeqn (1.6), so are the functions explicated as

u 1 = h ( t ε , x , y , z ) , v 1 = g ( t ε , x , y , z ) , u 2 = u e 2 ε h ( t e 3 ε , x e ε , y e 2 ε , z ) , v 2 = v e ε g ( t e 3 ε , x e ε , y e 2 ε , z ) , u 3 = h ( t , x , y , z ε ) , v 3 = g ( t , x , y , z ε ) , u 4 = u e 2 ε h ( t e 3 ε , x e ε , y , z e 2 ε ) , v 4 = v e ε g ( t e 3 ε , x e ε , y , z e 2 ε ) , u 5 = h ( t , x 2 ε z , y , z ) , v 5 = g ( t , x 2 ε z , y , z ) 3 ε , u 6 = h ( t , x ε t , y , z ) + ε , v 6 = g ( t , x ε t , y , z ) , u 7 = h ( t , x ε , y , z ) , v 7 = g ( t , x ε , y , z ) , u 8 = h ( t , x , y 2 ε z , z ) 3 ε , v 8 = g ( t , x , y 2 ε z , z ) , u 9 = h ( t , x , y ε , z ) + ε , v 9 = g ( t , x , y ε , z ) ,

with u i ( t , x , y , z ) = G ε i h ( t , x , y , z ) , i = 1 , 2 , , 9 , where ε 1 is conjectured as any real-positive number.

2.3 Optimal system of Lie subalgebras of (1.6)

This section tends to utilize the accessible chance of a symmetry group in the computation of one-dimensional Lie optimal system ascribed to Lie subalgebra [17,64] concerning (1.6). Therefore, deciding on subgroups of a symmetry group allows one to gain various kinds of solutions with a well-standardized approach which is stressed. The piece of work involved in sub-algebraic classification that is single-dimensional is a connective factor compared to the signalized one with orbit codification of the Lie adjoint-representation [17]. Now, an subalgebraic optimal set is achieved through the choice of single representative of any given equivalent sub-algebraic class. Additionally, it is possible to decipher the involved issue in orbit classification via the engagement of a Lie algebra general member. Thereafter, simplification is done via diverse transformations of adjoint. So, grounded on a well-known algorithm outlined by Hu et al. [64], we institute an optimal system of one dimension for HD-SOLeqn (1.6). Now, we seek first the principal invariant function associated with (1.6), next, we calculate the transformation matrix and on the final analysis, classification of the Lie algebraic finite dimension [17,64] is computed for (1.6) under consideration.

2.3.1 Principal invariants of System (1.6)

In achieving the one dimensional sub-algebra optimal system of Lie algebra g associated with R 9 , it is obligatory to construct the basic invariant in order to aid the selection of representative elements. The table of commutation relations, explicated in Table 1, unveils various combinations of relations ascribed to Lie brackets concerning g . The ( i ; j ) th entry of Table 1 is then computed via Lie bracket occasioned by [ w i , w j ] = w i w j w j w i .

Table 1

Lie algebra commutator table for HD-SOLeqn (1.6)

[ w i , w j ] w 1 w 2 w 3 w 4 w 5 w 6 w 7 w 8 w 9
w 1 0 3 w 1 0 3 w 1 0 w 7 0 0 0
w 2 3 w 1 0 0 0 w 5 2 w 6 w 7 2 w 8 2 w 9
w 3 0 0 0 2 w 3 2 w 7 0 0 2 w 9 0
w 4 3 w 1 0 2 w 3 0 w 5 2 w 6 w 7 2 w 8 0
w 5 0 w 5 2 w 7 w 5 0 0 0 0 0
w 6 w 7 2 w 6 0 2 w 6 0 0 0 0 0
w 7 0 w 7 0 w 7 0 0 0 0 0
w 8 0 2 w 8 2 w 9 2 w 8 0 0 0 0 0
w 9 0 2 w 9 0 0 0 0 0 0 0

We observe that Table 1 is skew symmetric having got zero diagonal elements. The real function Φ is taken as the invariant and it satisfies the relation Φ ( w ) = Φ ( A d g ( w ) ) for all w g alongside any subgroup g . Thus, we let w = i = 1 9 as well as g = e v v = j = 1 9 b j w j , where for g , we have

(2.6) A d exp ( ε v ) ( w ) = e ε v w e ε v = w ε [ v , w ] + 1 2 ! ε 2 [ v , [ v , w ] ] = ( a 1 w 1 + + a 9 w 9 ) ε [ b 1 w 1 + + b 9 w 9 , a 1 w 1 + + a 9 w 9 ] + O ( ε 2 ) = ( a 1 w 1 + + a 9 w 9 ) ε ( Σ 1 w 1 + + Σ 9 w 9 ) + O ( ε 2 ) ,

with function Σ i Σ i ( a 1 , , a 9 , b 1 , , b 9 ) secured via basic computations with the aid of Table 1. In consequence, we have the values of Σ i as

(2.7) Σ 1 = 3 a 1 b 2 3 a 1 b 4 + 3 a 2 b 1 + 3 a 4 b 1 , Σ 2 = 0 , Σ 3 = 2 a 3 b 4 + 2 a 4 b 3 , Σ 4 = 0 , Σ 5 = a 2 b 5 a 4 b 5 a 5 b 2 + a 5 b 4 , Σ 6 = 2 a 2 b 6 2 a 4 b 6 + 2 a 6 b 2 + 2 a 6 b 4 , Σ 7 = a 1 b 6 + a 2 b 7 2 a 3 b 5 + a 4 b 7 + 2 a 5 b 3 + a 6 b 1 a 7 b 2 a 7 b 4 , Σ 8 = 2 a 2 b 8 2 a 4 b 8 + 2 a 8 b 2 + 2 a 8 b 4 , Σ 9 = 2 a 2 b 9 2 a 3 b 8 + 2 a 8 b 3 + 2 a 9 b 2 .

Now, for any b j , where 1 j 9 , it is required that we have

(2.8) Σ 1 Φ a 1 + Σ 2 Φ a 2 + Σ 3 Φ a 3 + Σ 4 Φ a 4 + Σ 5 Φ a 5 + Σ 6 Φ a 6 + Σ 7 Φ a 7 + Σ 8 Φ a 8 + Σ 9 Φ a 9 = 0 .

Thus, by equating the coefficients of all same powers of β j in (2.8), one achieves the needed nine DEs with regards to real-valued function invariant Φ ( a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 ) as

(2.9) 2 a 2 Φ a 9 = 0 , ( a 2 + a 4 ) Φ a 7 = 0 , ( a 2 a 4 ) Φ a 5 2 a 3 Φ a 7 = 0 ( 3 a 2 + 3 a 4 ) Φ a 1 + a 6 Φ a 7 = 0 , ( 2 a 2 2 a 4 ) Φ a 8 2 a 3 Φ a 9 = 0 , ( 2 a 2 2 a 4 ) Φ a 6 a 1 Φ a 7 = 0 , 2 a 4 Φ a 8 + 2 a 5 Φ a 7 + 2 a 8 Φ a 9 = 0 , 3 a 1 Φ a 1 a 5 Φ a 5 + 2 a 6 Φ a 6 a 7 Φ a 7 + 2 a 8 Φ a 8 + 2 a 9 Φ a 9 = 0 , 3 a 1 Φ a 1 2 a 3 Φ a 3 + a 5 Φ a 5 + 2 a 6 Φ a 6 a 7 Φ a 7 + 2 a 8 Φ a 8 = 0 .

On solving the system of equations displayed in (2.9), one secures the value of invariant Φ ( a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 ) = G ( a 2 , a 4 ) , which we call the principal invariant function of Lie algebra g associated with R 9 . In this case, function G is an arbitrary function depending on a 2 and a 4 . In consequence, HD-SOLeqn (1.6) has only two basic invariants representing the killing form as depicted by Olver [17].

2.3.2 Adjoint transformation matrix of (1.6)

Let g be the symmetry Lie algebra having basis { w 1 , w 2 , w 3 , , w 9 } of Section 2.1 and also identifying with R 9 as a vector space imploring the map w i e i with { e 1 , e 2 , e 3 , , e 9 } regarded as the standard basis of R 9 . Suppose g = exp ( ε i w i ) for real constants ε i , i = 1 , 2 , 3 , , 9 , by reconing Table 2 and calculating the exponential of matrices ε i A d ( w i ) , we gain adjoint matrix representations of A d ( w i ) . For instance, from Table 2, whose ( i ; j ) th entry is A d exp ε w i ( w j ) , we use the relation

A d exp ( ε 1 w 1 ) w = a 1 A d exp ( ε 1 w 1 ) w 1 + a 2 A d exp ( ε 1 w 1 ) w 2 + + a 9 A d exp ( ε 1 w 1 ) w 9 = ( a 1 3 ε 1 a 2 3 ε 1 a 4 ) w 1 + a 2 w 2 + + a 5 w 5 ε 1 a 6 w 6 + + a 9 w 9 , = ( a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 ) A 1 ε ( w 1 , w 2 , w 3 , w 4 , w 5 , w 6 , w 7 , w 8 , w 9 ) T ,

whereby we have A 1 ε to be

Table 2

Adjoint representation table of Lie algebra for HD-SOLeqn (1.6)

A d g w 1 w 2 w 3 w 4 w 5 w 6 w 7 w 8 w 9
w 1 w 1 Q 0 w 3 w 4 3 ε 1 w 1 w 5 Q 12 w 7 w 8 w 9
w 2 e 3 ε 2 w 1 w 2 w 3 w 4 e ε 2 w 5 e 2 ε 2 w 6 e ε 2 w 7 e 2 ε 2 w 8 e 2 ε 2 w 9
w 3 w 1 w 2 w 3 w 4 2 ε 3 w 3 Q 10 w 6 w 7 Q 11 w 9
w 4 e 3 ε 4 w 1 w 2 e 2 ε 4 w 3 w 4 e ε 4 w 5 e 2 ε 4 w 6 e ε 4 w 7 e 2 ε 4 w 8 w 9
w 5 w 1 Q 1 Q 8 w 4 + ε 5 w 5 w 5 w 6 w 7 w 8 w 9
w 6 Q 3 Q 4 w 3 w 4 + 2 ε 6 w 6 w 5 w 6 w 7 w 8 w 9
w 7 w 1 Q 5 w 3 w 4 ε 7 w 7 w 5 w 6 w 7 w 8 w 9
w 8 w 1 Q 6 Q 9 w 4 + 2 ε 8 w 8 w 5 w 6 w 7 w 8 w 9
w 9 w 1 Q 7 w 3 w 4 w 5 w 6 w 7 w 8 w 9

Q 0 = w 2 3 ε 1 w 1 , Q 1 = w 2 ε 5 w 5 , Q 3 = w 1 + ε 6 w 7 , Q 4 = w 2 + 2 ε 6 w 6 , Q 5 = w 2 ε 7 w 7 , Q 6 = w 2 + 2 ε 8 w 8 , Q 7 = w 2 + 2 ε 9 w 9 , Q 8 = w 3 + 2 ε 5 w 7 , Q 9 = w 3 + 2 ε 8 w 9 , Q 10 = w 5 2 ε 3 w 7 , Q 11 = w 8 2 ε 3 w 9 , Q 12 = w 6 ε 1 w 7 .

A 1 ε = 1 0 0 0 0 0 0 0 0 3 ε 1 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 3 ε 1 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 ε 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 .

Therefore, by adopting the same process, one secures the other eight adjoint transformation matrices as

A 2 ε = e 3 ε 2 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 e ε 2 0 0 0 0 0 0 0 0 0 e 2 ε 2 0 0 0 0 0 0 0 0 0 e ε 2 0 0 0 0 0 0 0 0 0 e 2 ε 2 0 0 0 0 0 0 0 0 0 e 2 ε 2 ,

A 3 ε = 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 2 ε 3 1 0 0 0 0 0 0 0 0 0 1 0 2 ε 3 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 2 ε 3 0 0 0 0 0 0 0 0 1 ,

A 4 ε = e 3 ε 4 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 e 2 ε 4 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 e ε 4 0 0 0 0 0 0 0 0 0 e 2 ε 4 0 0 0 0 0 0 0 0 0 e ε 4 0 0 0 0 0 0 0 0 0 e 2 ε 4 0 0 0 0 0 0 0 0 0 1 ,

A 5 ε = 1 0 0 0 0 0 0 0 0 0 1 0 0 ε 5 0 0 0 0 0 0 1 0 0 0 2 ε 5 0 0 0 0 0 1 ε 5 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 ,

A 6 ε = 1 0 0 0 0 0 ε 6 0 0 0 1 0 0 0 2 ε 6 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 2 ε 6 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 ,

A 7 ε = 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 ε 7 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 ε 7 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 ,

A 8 ε = 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 2 ε 8 0 0 0 1 0 0 0 0 0 2 ε 8 0 0 0 1 0 0 0 2 ε 8 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 ,

A 9 ε = 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 2 ε 9 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 .

In consequence, we secure the general adjoint transformation matrix as

A ε = e 3 ε 2 + 3 ε 4 0 0 0 0 0 e 3 ε 2 + 3 ε 4 ε 6 0 0 A 21 ε 1 0 0 ε 5 2 ε 6 A 27 ε 2 ε 8 2 ε 9 0 0 e 2 ε 4 0 0 0 2 e 2 ε 4 ε 5 0 2 e 2 ε 4 ε 8 A 41 ε 0 2 e 2 ε 4 ε 3 1 ε 5 2 ε 6 A 47 ε 2 ε 8 4 e 2 ε 4 ε 3 ε 8 0 0 0 0 e ε 2 ε 4 0 2 e ε 2 + ε 4 ε 3 0 0 0 0 0 0 0 e 2 ε 2 2 ε 4 e ε 2 + ε 4 ε 1 0 0 0 0 0 0 0 0 e ε 2 + ε 4 0 0 0 0 0 0 0 0 0 e 2 ε 2 2 ε 4 2 e 2 ε 2 ε 3 0 0 0 0 0 0 0 0 e 2 ε 2 ,

where A 21 ε = 3 ε 1 e 3 ε 2 + 3 ε 4 , A 41 ε = 3 ε 1 e 3 ε 2 + 3 ε 4 , A 27 ε = 3 ε 1 ε 6 e 3 ε 2 + 3 ε 4 ε 7 , A 47 ε = 4 ε 3 ε 5 e 2 ε 4 3 ε 1 ε 6 e 3 ε 2 + 3 ε 4 ε 7 .

2.3.3 Adjoint transformation equation and Lie subalgebras of (1.6)

Here, the adjoint transformation equation associated with parameters a 1 , a 2 , a 3 , , a 9 is presented as

(2.10) ( α ˜ 1 , α ˜ 2 , α ˜ 3 , α ˜ 4 , α ˜ 5 , α ˜ 6 , α ˜ 7 , α ˜ 8 , α ˜ 9 ) = ( a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 , a 9 ) A ε ,

where A ε is the global adjoint matrix. One can assert here, that v = i = 1 9 α ˜ i w i as well as w = i = 1 9 a i w i coincide under the adjoint action. So, the involved adjoint equations associated with (1.6) with regards to system (2.10) are computed as

(2.11) α ˜ 1 = a 1 e 3 ε 2 + 3 ε 4 3 a 2 e 3 ε 2 + 3 ε 4 ε 1 3 a 4 e 3 ε 2 + 3 ε 4 ε 1 , α ˜ 2 = a 2 , α ˜ 3 = a 3 e 2 ε 4 2 a 4 e 2 ε 4 ε 3 , α ˜ 4 = a 4 , α ˜ 5 = a 5 e ε 2 ε 4 a 2 ε 5 + a 4 ε 5 , α ˜ 6 = a 6 e 2 ε 2 2 ε 4 + 2 a 2 ε 6 + 2 a 4 ε 6 , α ˜ 7 = a 7 e ε 2 + ε 4 a 6 e ε 2 + ε 4 ε 1 2 a 5 e ε 2 + ε 4 ε 3 + 2 a 3 e 2 ε 4 ε 5 + a 1 e 3 ε 2 + 3 ε 4 ε 6 a 2 ( 3 e 3 ε 2 + 3 ε 4 ε 1 ε 6 + ε 7 ) a 4 ( 4 e 2 ε 4 ε 3 ε 5 + 3 e 3 ε 2 + 3 ε 4 ε 1 ε 6 + ε 7 ) , α ˜ 8 = a 8 e 2 ε 2 2 ε 4 + 2 a 2 ε 8 + 2 a 4 ε 8 , α ˜ 9 = a 9 e 2 ε 2 2 a 8 e 2 ε 2 ε 3 + 2 a 3 e 2 ε 4 ε 8 4 a 4 e 2 ε 4 ε 3 ε 8 + 2 a 2 ε 9 .

Remark 2.1

It is noteworthy to state clearly here that the existence of solution of system (2.10) with regards to real constants ε i , i = 1 , 2 , 3 , , 9 implies that the selected element that occasions the solution is an optimal system and so it is germane.

Next, we begin the optimal system computations proper in the subsequent part of the study by first contemplating the invariants earlier secured on the basis of its sign as presented in the algorithm adopted [64] in two stages as a 2 = 1 , a 4 = 1 , and a 2 a 4 = 0 .

Case 1. a 2 = 1 , a 4 = 1

We choose in this situation, the representative element v = w 2 + w 4 . On inserting the parametric values α ˜ i = 0 , i = 1 , 3 , 5 , , 9 alongside α ˜ i = 1 , i = 2 , 4 into adjoint system (2.11), we achieve the solution

(2.12) ε 1 = 1 6 a 1 , ε 3 = 1 2 a 3 , ε 6 = 1 4 a 6 e 2 ( ε 2 + ε 4 ) , ε 7 = 1 12 a 6 a 1 e ε 2 + ε 4 + 1 2 a 7 e ε 2 + ε 4 , ε 8 = 1 4 a 8 e 2 ( ε 2 + ε 4 ) , ε 9 = 1 2 a 8 e 2 ε 2 a 3 1 2 a 9 e 2 ε 2 .

Case 2. a 2 a 4 = 0

This case occasions three scenarios a 2 = 0 , a 4 0 ; a 2 0 , a 4 = 0 ; and a 2 = 0 , a 4 = 0 . Thus, without loss of generality, we contemplate a 4 = a 2 = 1 and investigate those cases shortly.

Case 2.1. a 2 = 0 , a 4 = 1

We select in this subcase, the representative element v = w 4 . On invoking the parameters α ˜ i = 0 , i = 1 , 2 , 3 , 5 , , 9 together with α ˜ i = 1 , i = 4 into adjoint system (2.11), we gain the result

(2.13) ε 1 = 1 3 a 1 , ε 3 = a 9 2 a 8 , ε 4 = ln a 5 ε 5 + ε 2 , ε 6 = a 6 ε 5 2 2 a 5 2 e 4 ε 2 , ε 7 = 1 3 ε 5 a 8 { e 2 ε 2 a 5 ( a 1 a 6 a 8 + 3 a 5 a 9 3 a 7 a 8 ) } , ε 8 = a 8 ε 5 2 2 a 5 2 e 4 ε 2 ,

Case 2.2. a 2 = 1 , a 4 = 0

We choose the representative element v = w 2 . Substituting the parametric values α ˜ i = 0 , i = 1 , 3 , 4 , 5 , , 9 as well as α ˜ i = 1 , i = 2 into adjoint system (2.11), we secure the outcome

(2.14) ε 1 = 1 3 a 1 , ε 4 = ln a 5 ε 5 + ε 2 , ε 6 = a 6 ε 5 2 2 a 5 2 e 4 ε 2 , ε 7 = 1 3 ε 5 { e 2 ε 2 a 5 ( a 1 a 6 + 6 a 5 ε 3 3 a 7 ) } , ε 8 = a 8 ε 5 2 2 a 5 2 e 4 ε 2 , ε 9 = a 8 ε 3 e 2 ε 2 1 2 a 9 e 2 ε 2 .

Case 2.3. a 2 = 0 , a 4 = 0

Here, we insert a 2 = 0 and a 4 = 0 into equation PDE system (2.9), solve the resultant equations, and thereby secure a new invariant with regards to a i , i = 1 , 3 , 5 , , 9 as

(2.15) Φ ( a 1 , a 3 , a 5 , a 6 , a 7 , a 8 , a 9 ) = G a 3 a 5 a 1 3 , a 6 a 1 2 3 , a 8 a 1 2 3 .

In consequence, we treat this invariants accordingly in the succeeding part of this work by contemplating three basic cases where the invariants are 0 , 1 , and 1 [64].

Case 2.3.1. a 3 a 5 a 1 3 = 1 , a 6 a 1 2 3 = 1 , a 8 a 1 2 3 = 0 .

In this situation, we have a 3 3 a 5 3 = a 6 3 2 . Hence, by taking a 8 = 0 , a 1 0 , precisely a 1 > 0 , we select the representative element v = w 3 + w 5 w 6 . As a result, we implore the parametric values α ˜ i = 0 , i = 2 , 4 , , 9 , α ˜ i = 1 , i = 1 , 3 , 5 with α 6 = 1 in adjoint system (2.11) and so gain the solution

(2.16) ε 2 = 1 2 ln 1 a 5 1 a 1 3 , ε 4 = ln ( a 5 ) + 1 2 ln 1 a 5 1 a 1 3 , ε 6 = 1 a 1 2 a 1 a 5 ε 3 1 a 1 3 a 1 a 7 1 a 1 3 2 a 1 ε 5 ε 1 , ε 8 = 1 2 a 1 a 5 a 9 1 a 1 2 3 .

Case 2.3.2. a 3 a 5 a 1 3 = 1 , a 6 a 1 2 3 = 0 , a 8 a 1 2 3 = 1 .

Now, just as earlier demonstrated, here, we also have a 3 3 a 5 3 = a 8 3 2 and by assuming that a 6 = 0 , a 1 0 , we choose the representative element v = w 3 + w 5 w 8 . Thus, we insert the parameters α ˜ i = 0 , i = 2 , 4 , 6 , 7 , 9 together with α ˜ i = 1 , i = 1 , 3 , 5 , and α 8 = 1 into adjoint system (2.11) and obtain the outcome

(2.17) ε 2 = 1 2 ln 1 a 5 1 a 1 3 , ε 4 = ln ( a 5 ) + 1 2 ln 1 a 5 1 a 1 3 , ε 6 = 2 a 5 1 a 1 3 ε 3 a 7 1 a 1 3 2 ε 5 , ε 8 = 1 2 a 5 1 a 1 3 a 1 a 9 1 a 1 1 3 + 2 ε 3 .

Case 2.3.3. a 3 a 5 a 1 3 = 0 , a 6 a 1 2 3 = 1 , a 8 a 1 2 3 = 1

In third subcase of Case 2.3., we note that a 6 = a 8 and by supposing that a 5 = 0 , a 1 = a 3 0 , we select the representative element v = w 6 + w 8 . Therefore, by engaging the parametric values α ˜ i = 0 , i = 2 , 4 , 5 , 7 , 9 together with α ˜ i = 1 , i = 1 , 3 , 6 , 8 in adjoint system (2.11), we secure

(2.18) ε 2 = 1 2 ln ( a 3 a 8 ) , ε 3 = 1 2 a 8 ( 2 a 3 a 8 ε 8 + a 9 ) , ε 4 = 1 2 ln 1 a 3 , ε 6 = a 8 a 3 a 8 1 a 3 1 3 ε 1 a 7 a 3 a 8 1 a 3 1 3 2 ε 5 .

Case 2.3.4. a 3 a 5 a 1 3 = 1 , a 6 a 1 2 3 = 1 , a 8 a 1 2 3 = 1

This scenario presents a 1 = a 3 3 a 5 3 and also gives a 6 = a 8 . Thus, we choose the representative element v = w 1 + w 3 + w 5 . On inserting the parameters α ˜ i = 0 , i = 2 , 4 , 7 , 9 along with α ˜ i = 1 , i = 1 , 3 , 5 , 6 , 8 into adjoint system (2.11), we achieve the solution

(2.19) ε 2 = 1 2 ln 1 a 3 a 5 2 , ε 3 = 1 2 a 3 2 a 5 2 a 9 + a 3 ε 8 , ε 4 = ln ( a 5 ) + 1 2 ln 1 a 3 a 5 2 , ε 6 = 1 a 3 3 a 5 3 { a 3 4 a 5 5 a 9 2 a 3 3 a 5 3 ε 5 + 2 a 3 3 a 5 3 ε 8 a 3 2 a 5 2 a 7 + ε 1 } .

Case 2.3.5. a 3 a 5 a 1 3 = 1 , a 6 a 1 2 3 = 1 , a 8 a 1 2 3 = 1

Here, a 1 = a 3 3 a 5 3 and so for a 1 = a 3 > 0 , a 5 = a 6 = a 8 < 0 , we choose the representative element v = w 1 + w 3 w 5 . On substituting the parametric values α ˜ i = 0 , i = 2 , 4 , 7 , 9 , α ˜ i = 1 , i = 1 , 3 with α ˜ i = 1 , i = 5 , 6 , 8 into adjoint system (2.11), we obtain the result

(2.20) ε 2 = 1 4 ln a 8 a 5 2 , ε 3 = 1 2 a 8 2 a 8 a 5 2 ε 8 + a 9 , ε 4 = ln ( a 5 ) + 1 4 ln a 8 a 5 2 , ε 6 = 1 a 8 a 5 a 8 2 a 8 a 5 2 ε 1 + a 5 2 a 9 a 8 a 5 2 a 5 a 7 a 8 a 8 a 5 2 + 2 a 8 ε 5 2 a 8 ε 8 .

Furthermore, for a 1 = a 5 > 0 , a 3 = a 6 = a 8 < 0 , we select the representative element v = w 1 w 3 + w 5 . Thus, invoking the relevant parameters of α ˜ i , i = 1 , , 9 into system (2.11), we obtain

(2.21) ε 2 = 1 2 ln 1 a 3 a 5 2 , ε 3 = 1 2 a 3 2 a 5 2 a 9 a 3 ε 8 , ε 4 = ln ( a 5 ) + 1 2 ln 1 a 3 a 5 2 , ε 6 = 1 a 3 3 a 5 3 { a 3 4 a 5 5 a 9 + 2 a 3 3 a 5 3 ε 5 + 2 a 3 3 a 5 3 ε 8 + a 3 2 a 5 2 a 7 + ε 1 } .

Remark 2.2

We note that other possible representatives from Case 2.3.4. and Case 2.3.5. have been obtained earlier thereby contributing no additional subalgebra to the optimal system.

Case 2.3.6. a 3 a 5 a 1 3 = 0 , a 6 a 1 2 3 = 0 , a 8 a 1 2 3 = 0

In this situation, we take a 3 = a 6 = a 8 = 0 with a 1 0 and substitute it back into PDE system (2.9). Hence, we secure another new invariant dependent on ( a 1 , a 5 , a 7 , a 9 ) given as Φ ( a 1 , a 5 , a 7 , a 9 ) = G ( a 5 , a 9 , a 1 3 ) . Thus, we explore the invariant to further scale down vector (2.5) shortly, bearing in mind the approach engaged in the study of Hu et al. [64] as demonstrated earlier.

Case 2.3.6.1. a 1 0 , a 5 0 , a 9 0

Now, we contemplate a 1 > 0 , a 5 > 0 , and a 9 > 0 . Then, we select the optimal representative v = w 1 + w 5 + w 9 . Substituting values of parameters, viz, α ˜ i = 0 , i = 2 , 3 , 4 , 6 , , 8 alongside α ˜ i = 1 , i = 1 , 5 , 9 into system (2.11) gives the outcome

(2.22) ε 2 = 1 6 ln 1 a 1 a 5 3 , ε 3 = 1 2 a 1 a 5 1 a 1 a 5 3 2 3 ε 6 + a 7 2 a 5 , ε 4 = 1 6 ln a 5 3 a 1 .

Conversely, we consider a 1 < 0 , a 5 < 0 , and a 9 < 0 . So, we choose the optimal representative v = w 1 w 5 w 9 . On invoking parameters α ˜ i = 0 , i = 2 , 3 , 4 , 6 , , 8 alongside α ˜ i = 1 , i = 1 , 5 , 9 into adjoint Eq. (2.11), we have the result as

(2.23) ε 2 = 1 6 ln 1 a 1 a 5 3 , ε 3 = 1 2 a 1 a 5 1 a 1 a 5 3 2 3 ε 6 + a 7 2 a 5 , ε 4 = 1 6 ln a 5 3 a 1 .

Moreover, contemplating a 1 > 0 , a 5 > 0 with a 9 < 0 which gives representative v = w 1 + w 5 w 9 and also a 1 < 0 , a 5 < 0 with a 9 > 0 that secures representative v = w 1 w 5 + w 9 and solving system (2.11) using appropriate values of α i , i = 1 , 2 , 3 , , 9 gives the same solution as earlier obtained. Furthermore, for a 1 > 0 , a 5 < 0 with a 9 < 0 whose representative is v = w 1 w 5 w 9 as well as a 1 < 0 , a 5 > 0 with a 9 > 0 which gives representative v = w 1 + w 5 + w 9 together with a 1 > 0 , a 5 < 0 with a 9 > 0 that occasions representative v = w 1 w 5 + w 9 , we solve system (2.11). On imploring adequate values α i , i = 1 , 2 , 3 , , 9 , one obtains

(2.24) ε 2 = 1 6 ln 1 a 1 a 5 3 , ε 3 = 1 2 a 1 a 5 1 a 1 a 5 3 2 3 ε 6 + a 7 2 a 5 , ε 4 = 1 6 ln a 5 3 a 1 .

Next, we contemplate the case a 5 a 9 a 1 3 = 0 whereby various other possibilities are explored.

Case 2.3.6.2. a 1 = 0 , a 5 0 , a 9 0

Here, we examine the case a 5 > 0 and a 9 > 0 . Therefore, we select the representative element v = w 5 + w 9 . In addition, for a 5 < 0 and a 9 > 0 whose representative is v = w 5 + w 9 , we substitute appropriate values of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11) as earlier manifested and gain

(2.25) ε 2 = 1 2 ln ( a 9 ) , ε 3 = a 7 2 a 5 , ε 4 = 1 2 ln ( a 5 2 a 9 ) .

Moreover, we study the case when a 5 < 0 and a 9 < 0 . We choose the representative element v = w 5 w 9 . Besides, for a 5 > 0 and a 9 < 0 whose representative element is v = w 5 w 9 , we insert adequate parametric values of α i , i = 1 , 2 , 3 , , 9 into adjoint system (2.11) as exercised before and obtain

(2.26) ε 2 = 1 2 ln 1 a 9 , ε 3 = a 7 2 a 5 , ε 4 = 1 2 ln 1 a 5 2 a 9 , a 9 < 0 .

Case 2.3.6.3. a 5 = 0 , a 1 0 , a 9 0

Reckoning a 1 > 0 and a 9 > 0 with a 5 = 0 , we select the representative element v = w 1 + w 9 . Not only that on making adequate choice of α i , i = 1 , 2 , 3 , , 9 and substituting same in adjoint system (2.11), we secure the solution

(2.27) ε 2 = 1 3 ln 1 a 1 1 3 ln 1 a 1 a 9 3 , ε 4 = 1 3 ln 1 a 1 a 9 3 , ε 6 = a 7 a 1 a 9 ( a 9 a 1 2 3 ) .

Besides, for a 1 < 0 and a 9 < 0 alongside a 5 = 0 , we choose the representative element v = w 1 w 9 . On invoking adequate parametric values of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11) as revealed earlier, one obtains

(2.28) ε 2 = 1 3 ln 1 a 1 1 3 ln 1 a 1 a 9 a 9 , ε 4 = 1 3 ln 1 a 1 a 9 a 9 , ε 6 = a 7 a 1 a 9 ( ( a 1 a 9 ) 2 3 ( a 9 ) 3 ) .

Considering a 1 < 0 and a 9 > 0 with a 5 = 0 , we select the representative element v = w 1 + w 9 . Invoking adequate choice of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11), we secure the solution

(2.29) ε 2 = 1 3 ln 1 a 1 1 3 ln 1 a 1 a 9 3 , ε 4 = 1 3 ln 1 a 1 a 9 3 , ε 6 = a 7 a 1 ( a 1 2 3 ) .

In the same vein, for a 1 > 0 and a 9 < 0 alongside a 5 = 0 , we choose the representative element v = w 1 w 9 . On invoking adequate parametric values of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11) as revealed earlier, one obtains

(2.30) ε 2 = 1 3 ln 1 a 1 1 3 ln 1 a 1 a 9 a 9 , ε 4 = 1 3 ln 1 a 1 a 9 a 9 , ε 6 = a 7 a 1 a 9 ( ( a 1 a 9 ) 2 3 ( a 9 ) 3 ) .

Case 2.3.6.4. a 9 = 0 , a 1 0 , a 5 0

Regarding the case of a 1 > 0 and a 5 > 0 together with a 9 = 0 , we select representative element v = w 1 + w 5 . In addition, for a 1 < 0 , a 5 < 0 as well as a 9 = 0 , we choose representative element v = w 1 w 5 . On inserting relevant values of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11), we achieve the result

(2.31) ε 2 = 1 6 ln 1 a 1 a 5 3 , ε 3 = 1 2 a 1 a 5 1 a 1 a 5 3 2 3 ε 6 + a 7 2 a 5 , ε 4 = 1 6 ln a 5 3 a 1 .

On the contrary, if we reckon a 1 > 0 , a 5 < 0 as well as a 9 = 0 , we choose representative element v = w 1 w 5 . Besides, for a 1 < 0 , a 5 > 0 as well as a 9 = 0 , we choose representative element v = w 1 + w 5 . On invoking pertinent values of α i , i = 1 , 2 , 3 , , 9 in adjoint system (2.11), we secure the outcome

(2.32) ε 2 = 1 6 ln 1 a 1 a 5 3 , ε 3 = 1 2 a 1 a 5 1 a 1 a 5 3 2 3 ε 6 + a 7 2 a 5 , ε 4 = 1 6 ln a 5 3 a 1 .

Remark 2.3

It is noteworthy to state that some of the representative elements obtained with both positive and negative signs are equivalent. For instance, w 1 + w 5 + w 9 and w 1 w 5 w 9 are equivalent. The study reveals that their solutions ((2.22) and (2.23)) are exactly the same. Moreover, we observe the same occurrence for representative elements w 1 + w 5 and w 1 w 5 .

Finally, in view of the detailed calculations and analysis presented alongside remarks (2.2) and (2.3), we reduce the list of the representatives slightly by admitting that the discrete symmetry ( t , x , y , z , u ) ( t , x , y , z , u ) not in the connected component of the identity of the full symmetry group maps w 1 + w 5 + w 9 to w 1 w 5 w 9 , w 5 + w 9 to w 5 + w 9 , w 1 + w 9 to w 1 w 9 , and w 1 + w 5 to w 1 w 5 , and so on, thereby minimizing the number of inequivalent subalgebras [17]. Thus, we arrive at the theorem.

Theorem 2.4

The optimal system of one-dimensional Lie subalgebra of HD-SOLeqn (1.6) comprises the list: w 2 + w 4 ; w 4 ; w 2 ; w 3 + w 5 w 6 ; w 6 + w 8 ; w 3 + w 5 w 8 ; w 1 + w 3 + w 5 ; w 1 + w 5 + w 9 ; w 5 + w 9 ; w 1 + w 9 ; w 1 + w 5 ; w 1 w 5 + w 9 .

We note that various invariant solutions associated with the subalgebras presented in Theorem 2.4 have been copiously explored [61,62]. Thus, we compute the conserved vectors related to the vectors in the subsequent part of the research paper.

3 Conserved currents of system (1.6) with applications

This part of the article reveals the calculation of conserved vectors related to (1.6) by invoking the Ibragimov’s theorem [65,66] for determining conserved quantities. Hence, some salient information are divulged in order to understand the technique.

3.1 Preliminaries

A new theorem in [65] was introduced by Ibragimov for the computations of conserved vectors associated with a given DE. In addition, availability of classical Lagrangian is not demanded for the theorem to function. Ibragimov’s technique suggests fundamentally that infinitesimal generators are uniquely associated with their conserved current. Furthermore, the concept stands on theavailability of adjoint equation related to nonlinear DEs. Thus, a detailed outline of the theorem is furnished shortly.

Formal Lagrangian and adjoint equation

Theorem 3.1

[65] The system of adjoint relations given as

(3.1) Ξ σ * ( x , Ψ , Ω , , Ψ ( s ) , Ω ( s ) ) δ ( Ω β Ξ β ) δ Ψ σ = 0 , σ = 1 , , α , Ω = ( Ω 1 , , Ω α ) ,

which exists for a known system of α sets of relations presented as

(3.2) Ξ σ ( x , Ψ , Ψ ( 1 ) , , Ψ ( s ) ) = 0 , σ = 1 , , α ,

where Ω = Ω ( x ) with κ independent together with α dependent variables accordingly explicated as x = ( x 1 , x 2 , , x κ ) and Ψ = ( Ψ 1 , Ψ 2 , , Ψ α ) , alongside the variational derivative designated as Euler–Lagrange operator, expressed for each σ , by the formal sum defined as

(3.3) δ δ Ψ σ = Ψ σ + s = 1 ( 1 ) s D i 1 D i s × δ δ Ψ i 1 , i 2 , , i s σ , i = 1 , , κ ,

owns the symmetries inherited by set of relations (3.2). The complete derivative D i explicates as

(3.4) D i = x i + Ψ i σ Ψ σ + Ψ i j σ × Ψ j σ + , i = 1 , , κ , j = 1 , , κ .

Noteworthy, it is to declare that suppose set of relations (3.2) admit a point transformation group, having a generator convey as

(3.5) = ξ i x i + φ σ Ψ σ , ξ i = ξ i ( x , Ψ ) , φ σ = φ σ ( x , Ψ ) ,

then the system of adjoint relations (3.1) admit the operator (3.5) whose extension to a variable Ω σ is given by

(3.6) S = ξ i x i + φ σ Ψ σ + φ * σ Ω σ

with coefficient φ * σ = φ * σ ( x , Ψ , Ω , ) defined as

(3.7) φ * σ = [ λ β σ Ω β + Ω σ D i ( ξ i ) ] ,

which is selected approximately and the included variable Ω σ = Ω σ ( x ) . λ β σ is achievable by the formula:

(3.8) λ β σ Ξ β = X ( Ξ σ ) .

Hence, the formal Lagrangian is explicated as

(3.9) = Ω σ Ξ σ ( x , Ψ , Ψ ( 1 ) , , Ψ ( s ) ) ,

with the adjoint relation to (3.2) stated as

(3.10) Ξ σ * ( x , Ψ , Ω , , Ψ ( s ) , Ω ( s ) ) = 0 ,

together with criteria (3.1) holding.

Theorem 3.2

Every nonlocal symmetry (3.5), Lie–Bäcklund, as well as Lie point, admitted by the system of (3.2) produce a conserved current for relations (3.2) alongside the adjoint (3.10), with the conserved current T = ( T 1 , , T κ ) having components T i and decided by

(3.11) T i = ξ i + Π σ Ψ i σ D j Ψ i j σ + D j D k Ψ i j k σ + + D j ( Π σ ) Ψ i j σ D k Ψ i j k σ + + D j D k ( Π σ ) Ψ i j k + , i , j , k = 1 , , κ ,

where Π σ is the involved Lie characteristic function is engendered as

(3.12) Π σ = φ σ ξ j Ψ j σ , σ = 1 , , α , j = 1 , , κ .

Remark 3.1

We assert here that a given system of (3.2) is said to be self-adjoint if after replacing Ω = Ψ in the system of adjoint relations (3.10) gives that same system. A more detailed knowledge about the proof and various other copious information on the results made available here, can be accessed by the reader in the studies of Ibragimov [65,66].

3.2 Derivation of conservation laws via Ibragimov’s theorem

We give the Ibragimov’s conservation theorem [65] in this part of the article to find the conserved vectors of the HD-SOLeqn (1.6). Utilizing the highlighted information earlier presented, one achieves the theorem:

Theorem 3.3

The HD-SOLeqn (1.6) as well as its adjoint equation is accordingly explicated as

(3.13) G α 1 3 u x z 2 u t y + 2 u x u y + 2 u u x y + 2 v u x x + 2 v x u x u x x x y , G α 2 u y v x , G α 1 * 3 p x z 2 p t y + 2 p x v x + 2 u p x y + 2 p x x v q y p x x x y = 0 , G α 2 * q x 2 p x u x = 0 ,

with a second-order formal Lagrangian given as

(3.14) = p ( 2 u x x v + 2 u u x y 2 u t y + 2 u x v x + 2 u x u y + 3 u x z ) + q ( u y v x ) p x x u x y .

Obviously, one can frankly state that from the adjoint relation (3.13) and remark (3.1), HD-SOLeqn (1.6) is not self-adjoint. By applying of the earlier highlighted facts in Theorem 3.2, we calculate the conserved currents associated with the 12 elements of the achieved optimal system of one-dimensional subalgebras as well as other Lie point symmetries attained for (1.6). These conserved currents are

T 1 t = 2 u x v x p + 2 u x x p v + 2 u x u y p + 3 u x z p + 2 u x y p u u t y p + u y q v x q p y u t p x x u x y , T 1 x = 2 p x u t v 2 v t u x p 2 u t x p v + p y u t u u t u y p 3 2 u t z p u t y p u + v t q p t u x x y 1 2 u t p x x y + 1 2 p x x u t y + p t x u x y + 3 2 p z u t , T 1 y = p x u t u u t u x p u t x p u + u t t p u t q 1 2 p x x x u t + 1 2 p x x u t x p t u t , T 1 z = 3 2 p x u t 3 2 u t x p ;

T 2 t = u y z p p y u z , T 2 x = 2 p x u z v 2 u x v z p 2 u x z p v + p y u z u u y u z p 3 2 u z z p u y z p u + v z q 1 2 u z p x x y + p x z u x y + 1 2 p x x u y z p z u x x y + 3 2 p z u z , T 2 y = p x u z u u x u z p u x z p u + u t z p u z q p t u z 1 2 p x x x u z + 1 2 p x x u x z , T 2 z = 2 u x v x p + 2 u x x p v + 2 u x u y p + 3 2 u x z p + 2 u x y p u 2 u t y p + u y q v x q p x x u x y + 3 2 p x u z ;

T 3 t = 2 y u y p y 2 u p y x u x p y 3 t u t p y + 3 t q u y 2 y p u y y + 6 t p u y u x 3 t q v x + 6 t p u x v x + 9 t p u x z + x p u x y + 6 t p u u x y 3 t u x y p x x + 6 t p v u x x 3 t p u t y , T 3 x = 2 p y u 2 + 3 p z u 2 p u y u 2 y p y u y u + 2 y p u y y u + 4 v p x u + x p y u x u + x p u x y u p x x y u + 3 t p y u t u 3 t p u t y u + 2 y p u y 2 q v 3 p u z + x q u y 3 y p z u y 2 y q v y + 3 y p u y z 4 y v u y p x 4 p v u x + 3 2 x p z u x + x p u y u x + 4 y p v y u x + 2 x v p x u x + 3 2 x p u x z + 4 y p v u x y + p x u x y 2 y p x y u x y y u y y p x x + 9 2 t p z u t + 1 2 x u x y p x x + y u y p x x y 1 2 x u x p x x y + 2 y p y u x x y x p x u x x y 3 t u x x y p t 3 t p u y u t + 6 t v p x u t 3 2 t p x x y u t + 3 t q v t 6 t p u x v t 9 2 t p u t z 2 x p u t y + 3 2 t p x x u t y + 3 t u x y p t x 6 t p v u t x ,

T 3 y = 2 p x u 2 2 q u 2 y u y p x u 5 p u x u + x p x u x u 2 y p u x y u x p u x x u p x x x u 2 p t u + 3 t p x u t u 3 t p u t x u x p u x 2 x q u x 2 y p u y u x + 2 y q v x 4 y p u x v x 6 y p u x z + 3 2 u x p x x + y u x y p x x 4 y p v u x x + 1 2 x p x x u x x + y u y p x x x 1 2 x u x p x x x + 2 y u y p t x u x p t + 5 p u t 3 t q u t 3 t p u x u t 3 2 t p x x x u t 3 t p t u t + 2 y p u t y + x p u t x + 3 2 t p x x u t x + 3 t p u t t ,

T 3 z = 3 u p x 3 y u y p x + 3 2 x u x p x + 9 2 t u t p x 9 2 p u x + 3 y p u x y 3 2 x p u x x 9 2 t p u t x ;

T 4 t = 2 p u y 2 u p y 2 z u z p y x u x p y 3 t u t p y + 3 t q u y + 2 z p u y z + x p u x y + 6 t p u y u x 3 t q v x + 6 t p u x v x + 9 t p u x z + 6 t p u u x y 3 t u x y p x x + 6 t p v u x x 3 t p u t y , T 4 x = 2 p y u 2 + 3 p z u + 2 z u z p y u 4 p u y u 2 z p u y z u + 4 v p x u + x p y u x u + x p u x y u p x x y u + 3 t p y u t u 3 t p u t y u + q v 6 p u z + 3 z p z u z + 2 z q v z 3 z p u z z + x q u y 2 z p u z u y + 4 z v u z p x 8 p v u x + 3 2 x p z u x 4 z p v z u x + x p u y u x + 2 x v p x u x + 3 2 x p u x z 4 z p v u x z + p x u x y + 2 z p x z u x y + u y p x x + z u y z p x x + 1 2 x u x y p x x z u z p x x y 1 2 x u x p x x y 2 z p z u x x y x p x u x x y 3 t u x x y p t + 9 2 t p z u t 3 t p u y u t + 6 t v p x u t 3 2 t p x x y u t + 3 t q v t 6 t p u x v t 9 2 t p u t z 2 x p u t y + 3 2 t p x x u t y + 3 t u x y p t x 6 t p v u t x ,

T 4 y = 2 p x u 2 2 q u + 2 z u z p x u 5 p u x u + x p x u x u 2 z p u x z u x p u x x u p x x x u 2 p t u + 3 t p x u t u 3 t p u t x u x p u x 2 2 z q u z x q u x 2 z p u z u x + 3 2 u x p x x + z u x z p x x + 1 2 x p x x u x x z u z p x x x 1 2 x u x p x x x 2 z u z p t x u x p t + 5 p u t 3 t q u t 3 t p u x u t 3 2 t p x x x u t 3 t p t u t + 2 z p u t z + x p u t x + 3 2 t p x x u t x + 3 t p u t t ,

T 4 z = 2 z q u y + 4 z p u x u y + 3 u p x + 3 z u z p x 9 2 p u x + 3 2 x p x u x 2 z q v x + 4 z p u x v x + 3 z p u x z + 4 z p u u x y 2 z u x y p x x 3 2 x p u x x + 4 z p v u x x + 9 2 t p x u t 4 z p u t y 9 2 t p u t x ;

T 5 t = 2 z u x y p 2 z p y u x , T 5 x = 4 z p x u x v + 3 u x p + 2 z p y u x u + 2 z u x u y p + 3 z u x z p + 2 z u x y p u 1 2 z u x x x y p 4 z u t y p 2 z u y q + 3 q + 1 2 z p x x u x y + z u x x p x y 3 2 z u x p x x y z p x u x x y 1 2 z p y u x x x + 3 z p z u x , T 5 y = 2 z p x u x u 2 z u x 2 p 2 z u x x p u + 1 2 z u x x x x p + 2 z u t x p + 2 z u x q 2 z p t u x 1 2 z p x x x u x + 1 2 z p x x u x x 1 2 z p x u x x x ,

T 5 z = 3 z p x u x 3 z u x x p ; T 6 t = t u x y p t p y u x p y , T 6 x = 2 t p x u x v + p y u u y p + t p y u x u + t u x u y p + 3 2 t u x z p + t u x y p u 1 4 t u x x x y p 2 t u t y p + 2 p x v t u y q + 1 4 t p x x u x y + 1 2 t u x x p x y 3 4 t u x p x x y 1 2 t p x u x x y 1 4 t p y u x x x + 3 2 t p z u x 3 4 p x x y + 3 2 p z ,

T 6 y = p x u t u x 2 p + t p x u x u t u x x p u + 1 4 t u x x x x p + t u t x p + t u x q + q 1 4 t p x x x u x t p t u x + 1 4 t p x x u x x 1 4 t p x u x x x p t 1 4 p x x x , T 6 z = 3 2 t p x u x 3 2 t u x x p + 3 2 p x ; T 7 t = u x y p p y u x , T 7 x = 2 p x u x v + u x u y p + p y u x u + 3 2 u x z p + u x y p u 1 4 u x x x y p 2 u t y p u y q + 1 4 p x x u x y + 1 2 u x x p x y 3 4 u x p x x y 1 2 p x u x x y 1 4 p y u x x x + 3 2 p z u x ,

T 7 y = p x u x u u x 2 p u x x p u + 1 4 u x x x x p + u t x p + u x q p t u x 1 4 p x x x u x + 1 4 p x x u x x 1 4 p x u x x x , T 7 z = 3 2 p x u x 3 2 u x x p ; T 8 t = 2 z u y p y + 3 p y + 2 z p u y y , T 8 x = 3 z p z u y 2 z p u y 2 + 2 z u p y u y + 4 z v p x u y 3 2 z p x x y u y 9 2 p z 3 u p y 2 z q v y 3 z p u y z 2 z p u u y y 6 v p x 4 z p v y u x 4 z p v u x y + z p x y u x y z p x u x y y + 1 2 z u y y p x x + 9 4 p x x y 1 2 z p y u x x y + 3 2 z p u x x y y ,

T 8 y = 2 z v x q 3 q 3 u p x + 2 z u u y p x + 3 p u x + 2 z p u y u x + 4 z p u x v x + 6 z p u x z + 2 z p u u x y + 1 2 z u x y p x x + 4 z p v u x x 1 2 z p x u x x y 1 2 z u y p x x x + 3 4 p x x x 3 2 z p u x x x y 2 z u y p t + 3 p t 2 z p u t y , T 8 z = 3 z u y p x 9 2 p x 3 z p u x y ;

T 9 t = u y y p p y u y , T 9 x = 2 p x u y v 2 u x v y p 2 u x y p v u y 2 p + p y u y u 3 2 u y z p u y y p u + 3 4 u x x y y p v y q 3 4 u y p x x y + 1 2 p x y u x y 1 2 p x u x y y + 1 4 p x x u y y 1 4 p y u x x y + 3 2 p z u y , T 9 y = 2 u x v x p + 2 u x x p v + p x u y u + u x u y p + 3 u x z p + u x y p u 3 4 u x x x y p u t y p + v x q p t u y 1 4 p x u x x y + 1 4 p x x u x y 1 4 p x x x u y ,

T 9 z = 3 2 p x u y 3 2 u x y p ; T 10 t = 2 y u y p y 4 u p y 2 z u z p y 2 x u x p y 6 t u t p y + 2 p u y + 6 t q u y + 2 z p u y z 2 y p u y y + 12 t p u y u x 6 t q v x + 12 t p u x v x + 18 t p u x z + 2 x p u x y + 12 t p u u x y 6 t u x y p x x + 12 t p v u x x 6 t p u t y ,

T 10 x = 4 p y u 2 + 6 p z u + 2 z u z p y u 6 p u y u 2 y p y u y u 2 z p u y z u + 2 y p u y y u + 8 v p x u + 2 x p y u x u + 2 x p u x y u 2 p x x y u + 6 t p y u t u 6 t p u t y u + 2 y p u y 2 9 p u z + 3 z p z u z + 2 z q v z 3 z p u z z + 2 x q u y 3 y p z u y 2 z p u z u y 2 y q v y + 3 y p u y z + 4 z v u z p x 4 y v u y p x 12 p v u x + 3 x p z u x 4 z p v z u x + 2 x p u y u x + 4 y p v y u x + 4 x v p x u x + 3 x p u x z 4 z p v u x z + 4 y p v u x y + 2 p x u x y + 2 z p x z u x y 2 y p x y u x y + u y p x x + z u y z p x x y u y y p x x + x u x y p x x z u z p x x y + y u y p x x y x u x p x x y 2 z p z u x x y + 2 y p y u x x y 2 x p x u x x y 6 t u x x y p t + 9 t p z u t 6 t p u y u t + 6 t q v t + 12 t v p x u t 3 t p x x y u t 12 t p u x v t 9 t p u t z 4 x p u t y + 3 t p x x u t y + 6 t u x y p t x 12 t p v u t x ,

T 10 y = 4 p x u 2 4 q u + 2 z u z p x u 2 y u y p x u 10 p u x u + 2 x p x u x u 2 z p u x z u 2 y p u x y u 2 x p u x x u 2 p x x x u 4 p t u + 6 t p x u t u 6 t p u t x u 2 x p u x 2 2 z q u z 2 x q u x 2 z p u z u x 2 y p u y u x + 2 y q v x 4 y p u x v x 6 y p u x z + 3 u x p x x + z u x z p x x + y u x y p x x 4 y p v u x x + x p x x u x x z u z p x x x + y u y p x x x x u x p x x x 2 z u z p t + 2 y u y p t 2 x u x p t + 10 p u t 6 t q u t 6 t p u x u t 3 t p x x x u t 6 t p t u t + 2 z p u t z + 2 y p u t y + 2 x p u t x + 3 t p x x u t x + 6 t p u t t ,

T 10 z = 2 z q u y 3 y p x u y + 4 z p u x u y + 6 u p x + 3 z u z p x 9 p u x + 3 x p x u x 2 z q v x + 4 z p u x v x + 3 z p u x z + 3 y p u x y + 4 z p u u x y 2 z u x y p x x 3 x p u x x + 4 z p v u x x + 9 t p x u t 4 z p u t y 9 t p u t x ;

T 11 t = q u y + 2 p u x u y 2 z p y u x q v x + 2 p u x v x + 3 p u x z + 2 z p u x y + 2 p u u x y p y u t + 2 p v u x x u x y p x x p u t y ,

T 11 x = 2 z u y q + v t q 3 q + 3 p u x + 3 z p z u x + 2 z u p y u x + 2 z p u y u x + 4 z v p x u x + 3 z p u x z + 2 z p u u x y + z u x y p x x z u x p x x y 2 z p x u x x y u x x y p t + 3 2 p z u t + u p y u t p u y u t + 2 v p x u t 1 2 p x x y u t 2 p u x v t 3 2 p u t z 4 z p u t y p u u t y + 1 2 p x x u t y + u x y p t x 2 p v u t x ,

T 11 y = 2 z p u x 2 2 z q u x + 2 z u p x u x z p x x x u x 2 z p t u x p u t u x 2 z p u u x x + z p x x u x x q u t + u p x u t 1 2 p x x x u t p t u t + 2 z p u t x p u u t x + 1 2 p x x u t x + p u t t , T 11 z = 3 z p x u x 3 z p u x x + 3 2 p x u t 3 2 p u t x ;

T 12 t = t u x p y u z p y 2 z u x p y + p y + p u y z t p u x y + 2 z p u x y , T 12 x = v z q t u y q + 2 z u y q 3 q 3 2 p z + 3 2 p z u z 3 2 p u z z u p y + u u z p y + p u y p u u y z + 2 v u z p x + 3 p u x 3 2 t p z u x + 3 z p z u x 2 p v z u x t u p y u x + 2 z u p y u x 2 v p x t p u y u x + 2 z p u y u x 2 t v p x u x + 4 z v p x u x 3 2 t p u x z + 3 z p u x z 2 p v u x z t p u u x y + 2 z p u u x y + p x z u x y + 1 2 u y z p x x 1 2 t u x y p x x + z u x y p x x 1 2 u z p x x y + 1 2 t u x p x x y z u x p x x y + 1 2 p x x y p z u x x y + t p x u x x y 2 z p x u x x y + 2 t p u t y 4 z p u t y p u z u y ,

T 12 y = t p u x 2 2 z p u x 2 + t q u x 2 z q u x p u z u x t u p x u x + 2 z u p x u x + 1 2 t p x x x u x z p x x x u x + t p t u x 2 z p t u x + q q u z u p x + u u z p x p u u x z + 1 2 u x z p x x + t p u u x x 2 z p u u x x 1 2 t p x x u x x + z p x x u x x 1 2 u z p x x x + 1 2 p x x x u z p t + p t + p u t z t p u t x + 2 z p u t x , T 12 z = q u y + 2 p u x u y + 3 2 u z p x 3 2 p x 3 2 t p x u x + 3 z p x u x q v x + 2 p u x v x + 3 2 p u x z + 2 p u u x y u x y p x x + 3 2 t p u x x 3 z p u x x + 2 p v u x x 2 p u t y ;

T 13 t = t p u x y 2 z u y p y t u x p y + 2 p y + 2 z p u y y , T 13 x = 2 z u p y u y 2 z p u y 2 p u y + t q u y + 3 z p z u y + 4 z v p x u y + t p u x u y z p x x y u y 3 p z 2 u p y + 2 z q v y 3 z p u y z 2 z p u u y y 4 v p x + 3 2 t p z u x + t u p y u x 4 z p v y u x + 2 t v p x u x + 3 2 t p u x z + t p u u x y 4 z p v u x y + 2 z p x y u x y + z u y y p x x + 1 2 t u x y p x x 1 2 t u x p x x y + p x x y 2 z p y u x x y t p x u x x y 2 t p u t y ,

T 13 y = 3 p u x t p u x 2 t q u x + 2 z p u y u x + t u p x u x + 4 z p v x u x 1 2 t p x x x u x t p t u x + 2 q 2 u p x + 2 z u u y p x 2 z q v x + 6 z p u x z + 2 z p u u x y z u x y p x x t p u u x x + 4 z p v u x x + 1 2 t p x x u x x z u y p x x x + p x x x 2 z u y p t + 2 p t 2 z p u t y + t p u t x , T 13 z = 3 z u y p x + 3 2 t u x p x 3 p x 3 z p u x y 3 2 t p u x x ;

T 14 t = u z p y + 2 z u y p y 2 z u x p y 3 p y + p u y z 2 z p u y y + 2 z p u x y , T 14 x = 2 z p u y 2 + 2 z q u y 3 z p z u y p u z u y 2 z u p y u y 4 z v p x u y + 2 z p u x u y + z p x x y u y 3 q + 9 2 p z + 3 2 p z u z + q v z 3 2 p u z z + 3 u p y + u u z p y 2 z q v y + 3 z p u y z p u u y z + 2 z p u u y y + 6 v p x + 2 v u z p x + 3 p u x + 3 z p z u x 2 p v z u x + 2 z u p y u x + 4 z p v y u x + 4 z v p x u x + 3 z p u x z 2 p v u x z + 2 z p u u x y + 4 z p v u x y + p x z u x y 2 z p x y u x y + 1 2 u y z p x x z u y y p x x + z u x y p x x 1 2 u z p x x y z u x p x x y 3 2 p x x y p z u x x y + 2 z p y u x x y 2 z p x u x x y 4 z p u t y ,

T 14 y = 2 z u p x u x 2 z p u x 2 3 p u x 2 z q u x p u z u x 2 z p u y u x 4 z p v x u x z p x x x u x 2 z p t u x 3 q q u z + 3 u p x + u u z p x 2 z u u y p x + 2 z q v x 6 z p u x z p u u x z 2 z p u u x y + 1 2 u x z p x x + z u x y p x x 2 z p u u x x 4 z p v u x x + z p x x u x x 1 2 u z p x x x + z u y p x x x 3 2 p x x x u z p t + 2 z u y p t 3 p t + p u t z + 2 z p u t y + 2 z p u t x , T 14 z = q u y 3 z p x u y + 2 p u x u y + 3 2 u z p x + 9 2 p x + 3 z p x u x q v x + 2 p u x v x + 3 2 p u x z + 3 z p u x y + 2 p u u x y u x y p x x 3 z p u x x + 2 p v u x x 2 p u t y ;

T 15 t = q u y u z p y 2 z u x p y u t p y + p u y z + 2 p u y u x q v x + 2 p u x v x + 3 p u x z + 2 z p u x y + 2 p u u x y u x y p x x + 2 p v u x x p u t y , T 15 x = v z q + 2 z u y q + v t q 3 q + 3 2 p z u z 3 2 p u z z + u u z p y p u z u y p u u y z + 2 v u z p x + 3 p u x + 3 z p z u x 2 p v z u x + 2 z u p y u x + 2 z p u y u x + 4 z v p x u x + 3 z p u x z 2 p v u x z + 2 z p u u x y + p x z u x y + 1 2 u y z p x x + z u x y p x x 1 2 u z p x x y z u x p x x y p z u x x y 2 z p x u x x y u x x y p t + 3 2 p z u t + u p y u t p u y u t + 2 v p x u t 1 2 p x x y u t 2 p u x v t 3 2 p u t z 4 z p u t y p u u t y + 1 2 p x x u t y + u x y p t x 2 p v u t x ,

T 15 y = u u z p x 2 z p u x 2 2 z q u x p u z u x + 2 z u p x u x z p x x x u x 2 z p t u x p u t u x q u z p u u x z + 1 2 u x z p x x 2 z p u u x x + z p x x u x x 1 2 u z p x x x u z p t q u t + u p x u t 1 2 p x x x u t p t u t + p u t z + 2 z p u t x p u u t x + 1 2 p x x u t x + p u t t , T 15 z = q u y + 2 p u x u y + 3 2 u z p x + 3 z p x u x q v x + 2 p u x v x + 3 2 p u x z + 2 p u u x y u x y p x x 3 z p u x x + 2 p v u x x + 3 2 p x u t 2 p u t y 3 2 p u t x ;

T 16 t = q u y p y u y + 2 p u x u y + p u y y 2 z p y u x + 3 p u x z + 2 z p u x y + 2 p u u x y u x y p x x + 2 p v u x x p y u t p u t y , T 16 x = 2 z q u y p u y 2 + 3 2 p z u y + u p y u y + 2 v p x u y + 2 z p u x u y 1 2 p x x y u y p u t u y 3 q + q v y 3 2 p u y z p u u y y + 3 p u x + 3 z p z u x + 2 z u p y u x 2 p v y u x + 4 z v p x u x + 3 z p u x z + 2 z p u u x y 2 p v u x y + p x y u x y + 1 2 u y y p x x + z u x y p x x z u x p x x y p y u x x y 2 z p x u x x y u x x y p t + 3 2 p z u t + u p y u t + 2 v p x u t 1 2 p x x y u t + q v t 2 p u x v t 3 2 p u t z 4 z p u t y p u u t y + 1 2 p x x u t y + u x y p t x 2 p v u t x ,

T 16 y = p u y u x 2 z p u x 2 2 z q u x + 2 z u p x u x + 2 p v x u x z p x x x u x 2 z p t u x p u t u x + u u y p x q v x + 3 p u x z + p u u x y 1 2 u x y p x x 2 z p u u x x + 2 p v u x x + z p x x u x x 1 2 u y p x x x u y p t q u t + u p x u t 1 2 p x x x u t p t u t p u t y + 2 z p u t x p u u t x + 1 2 p x x u t x + p u t t ,

T 16 z = 3 2 u y p x + 3 z u x p x + 3 2 u t p x 3 2 p u x y 3 z p u x x 3 2 p u t x ; T 17 t = p u y y p y u y 2 z p y u x + 2 z p u x y , T 17 x = 2 z q u y p u y 2 + 3 2 p z u y + u p y u y + 2 v p x u y + 2 z p u x u y 1 2 p x x y u y 3 q + q v y 3 2 p u y z p u u y y + 3 p u x + 3 z p z u x + 2 z u p y u x 2 p v y u x + 4 z v p x u x + 3 z p u x z + 2 z p u u x y 2 p v u x y + p x y u x y + 1 2 u y y p x x + z u x y p x x z u x p x x y p y u x x y 2 z p x u x x y 4 z p u t y ,

T 17 y = p u y u x 2 z p u x 2 2 z q u x + 2 z u p x u x + 2 p v x u x z p x x x u x 2 z p t u x + u u y p x q v x + 3 p u x z + p u u x y 1 2 u x y p x x 2 z p u u x x + 2 p v u x x + z p x x u x x 1 2 u y p x x x u y p t p u t y + 2 z p u t x , T 17 z = 3 2 u y p x + 3 z u x p x 3 2 p u x y 3 z p u x x ; T 18 t = q u y p y u y + 2 p u x u y + p u y y q v x + 2 p u x v x + 3 p u x z + 2 p u u x y u x y p x x + 2 p v u x x p y u t p u t y ,

T 18 x = p u y 2 + 3 2 p z u y + u p y u y + 2 v p x u y 1 2 p x x y u y p u t u y + q v y 3 2 p u y z p u u y y 2 p v y u x 2 p v u x y + p x y u x y + 1 2 u y y p x x p y u x x y u x x y p t + 3 2 p z u t + u p y u t 1 2 p x x y u t + q v t 2 p u x v t 3 2 p u t z p u u t y + 1 2 p x x u t y + u x y p t x + 2 v p x u t 2 p v u t x ,

T 18 y = u u y p x + u u t p x + p u y u x q v x + 2 p u x v x + 3 p u x z + p u u x y 1 2 u x y p x x + 2 p v u x x 1 2 u y p x x x u y p t q u t p u x u t 1 2 p x x x u t p t u t p u t y p u u t x + 1 2 p x x u t x + p u t t , T 18 z = 3 2 u y p x + 3 2 u t p x 3 2 p u x y 3 2 p u t x ;

T 19 t = q u y p y u y + 2 p u x u y + p u y y + 2 z p y u x q v x + 2 p u x v x + 3 p u x z 2 z p u x y + 2 p u u x y u x y p x x + 2 p v u x x p y u t p u t y , T 19 x = p u y 2 2 z q u y + 3 2 p z u y + u p y u y + 2 v p x u y 2 z p u x u y 1 2 p x x y u y p u t u y + q v y 3 2 p u y z p u u y y 3 p u x 3 z p z u x 2 z u p y u x 2 p v y u x 4 z v p x u x 3 z p u x z 2 z p u u x y 2 p v u x y + p x y u x y + 1 2 u y y p x x z u x y p x x + z u x p x x y p y u x x y + 2 z p x u x x y u x x y p t + 3 2 p z u t + u p y u t + 2 v p x u t + 3 q 1 2 p x x y u t + q v t 2 p u x v t 3 2 p u t z + 4 z p u t y p u u t y + 1 2 p x x u t y + u x y p t x 2 p v u t x ,

T 19 y = 2 z p u x 2 + 2 z q u x + p u y u x 2 z u p x u x + 2 p v x u x + z p x x x u x + 2 z p t u x p u t u x + u u y p x q v x + 3 p u x z + p u u x y 1 2 u x y p x x + 2 z p u u x x + 2 p v u x x z p x x u x x 1 2 u y p x x x u y p t q u t + u p x u t 1 2 p x x x u t p t u t p u t y 2 z p u t x p u u t x + 1 2 p x x u t x + p u t t , T 19 z = 3 2 u y p x 3 z u x p x + 3 2 u t p x 3 2 p u x y + 3 z p u x x 3 2 p u t x .

3.2.1 Physical meaning of some of the obtained conservation laws in physical sciences and engineering

Local conservation laws for the higher dimensional integrodifferential soliton Eq. (1.6) observe a divergence criterion on the whole space ε solution of Eq. (1.6)

(3.15) ( D z T z + D y T y + D x T x + D t T t ) ε = 0 ,

where ( D z , D y , D x , D t ) are total differential operators with conserved density T as well as T z , T y , T x , which are spatial fluxes are functions depending on ( z , y , x , t ) . In physical sense, we observe that every conservation law generates a related conserved integral [67,68]

(3.16) P [ u , v ] = Γ C t d x ,

as can be obtained in this case, where Γ denotes the solution domain u ( t , x , y , z ) . It has been established that conservation laws with discontinuous coefficients, such as fluxes, densities as well as source terms, emerge in quite a large number of problems in the fields of physics and engineering. In this research, the resultant conservation laws have been observed to have relevance in physical sciences. Physically speaking, part of the obtained symmetries that delineate conserved quantities, viz, conservation of energy, momentum, and so on. More precisely, it is observed that time translation symmetries purvey energy, whereas in the case of space translation alongside boost symmetries, momenta are attained.

Theorem 3.4

The set of 19 non-trivial local conservation laws generated by HD-SOLeqn (1.6) via Ibragimov’s theorem comprise components ( T j t , T j x , T j y , T j z , j = 1 , 2 , , 19 ) , with density as the first component and fluxes constituting the other components.

Consequently, in the case of HD-SOLeqn (1.6), for example, conserved vectors ( T 1 t , T 1 x , T 1 y , T 1 z ) yield conservation of energy obtained as

(3.17) P 1 [ u , v ] = Γ ( 2 u x v x p + 2 u x x p v + 2 u x u y p + 3 u x z p + 2 u x y p u u t y p + u y q v x q p y u t p x x u x y ) d x ,

for solution u ( t , x , y , z ) , corresponding to time translation symmetry. Besides, conserved vectors ( T 2 t , T 2 x , T 2 y , T 2 z ) , accordingly give the conserved quantity of momentum explicated as

(3.18) P 2 [ u , v ] = Γ ( u y z p p y u z ) d x ,

which relates to the corresponding space translation symmetry. In addition, one obtains conserved quantity of dilation energy as

(3.19) P 3 [ u , v ] = Γ ( 2 y u y p y 2 u p y x u x p y 3 t u t p y + 3 t q u y 2 y p u y y + 6 t p u y u x 3 t q v x + 6 t p u x v x + 9 t p u x z + x p u x y + 6 t p u u x y 3 t u x y p x x + 6 t p v u x x 3 t p u t y ) d x ,

which correspond to ( T 3 t , T 3 x , T 3 y , T 3 z ) , which furnishes conservation laws associated with scaling symmetry. Moreover, another conserved quantity of dilation energy is obtainable as

P 4 [ u , v ] = Γ ( 2 p u y 2 u p y 2 z u z p y x u x p y 3 t u t p y + 3 t q u y + 2 z p u y z + x p u x y + 6 t p u y u x 3 t q v x + 6 t p u x v x + 9 t p u x z + 6 t p u u x y 3 t u x y p x x + 6 t p v u x x 3 t p u t y ) d x ,

corresponding to ( T 4 t , T 4 x , T 4 y , T 4 z ) . The rest could also be studied in the same way. For more interesting insights, see the study of Adeyemo and Khalique [68]. In particle physics, various other conservation laws apply to subatomic particles’ properties which are invariant in the course of interactions. In addition, an essential function of conserved vectors is the fact that they occasion the possibility of making prediction for the macroscopic character of a system without taking cognizance of the microscopic details in the course of a chemical reaction or physical process [69].

3.3 Application of conserved quantities in physical sciences and engineering

Frankly speaking, conservation law, also referred to as the law of conservation, in physics, depicts a principle that states that a certain physical property (i.e., a measurable quantity) remains unchanged with time within an isolated physical system. Thus, in physical sciences and engineering mathematics, conservation laws state that a particular measurable property associated with an isolated physical system remains constant (that is does not change) as the system evolves over time [70].

Conservation laws [7173] alternatively, are an area of applicable engrossment in engineering and physics, with the inclusion of theoretical vis-à-vis quantum-mechanics. This part of our research scrutinizes the conserved quantities of model HD-SOLeqn (1.6) with an observable trait that showcases results unveiling the availability of conservation of momentum alongside that of energy. Physical quantities resident in isolated systems, namely, mass, charge, angular momentum, energy, along with linear momentum are conserved. Meanwhile, it is unveiled that conserved quantities invoke an advantageous feature of DE ntegrability check. Further to that, conserved quantities significantly engender the establishment of existence-uniqueness characteristic for linearization mappings and solutions. Besides, they bring to play stability analysis and to ascertain global behavior of solutions.

In addition, conserved quantities play a leading role in the evolution of numerical techniques. They also furnish a crucial starting point in securing non-locally related systems and potential variables. In particular, a conserved quantity is fundamental in the investigation of a given DE, which implies that it holds for any posed data whether initial conditions and/or boundary conditions. Furthermore, the conservation law structure is such that it is not depending on co-ordinates since it involves a contact transformation mapping one to the other. Momentum is simply referred to as the resistance of a given object with regards to an alteration in the object’s velocity. Engineers employ this concept to make lives safer for people through the design of products in lengthening the time over which a deceleration happens.

Moreover, momentum is crucial in Physics due to the fact that it recounts the connection that is existing between mass, speed, as well as direction. Besides, it also delineates the force required to put an object to a halt and to further keep them in motion. It is to be noted that an apparently small object can deploy a large amount of force, provided it possesses sufficient momentum. Figures 1 and 2 present some practical applications of momentum.

Figure 1 
                  In an experiment recently conducted at Glasgow University, Heisenberg’s uncertainty principle is manifested in which a light beam that initially possesses no orbital angular momentum was observed. We define 
                        
                           
                           
                              Δ
                              ϕ
                           
                           \Delta \phi 
                        
                      as the range of azimuthal angular positions appearing for a photon in a cross-sectional part of the beam through a sector aperture. Moreover, at the upstream of the said aperture, the beam is observed to be in an 
                        
                           
                           
                              ℓ
                              =
                              0
                           
                           \ell =0
                        
                      eigen state of orbital angular momentum (i.e., upper panel). Nonetheless, the uncertainty principle states that the limitation in 
                        
                           
                           
                              ϕ
                           
                           \phi 
                        
                      causes a spread in orbital angular momentum represented as 
                        
                           
                           
                              L
                              =
                              ℓ
                              ℏ
                           
                           L=\ell \hslash 
                        
                      for every photon. Additionally, for strait apertures, the connection is portrayed as 
                        
                           
                           
                              Δ
                              ϕ
                              Δ
                              L
                              ≥
                              ℏ
                              ∕
                              2
                           
                           \Delta \phi \Delta L\ge \hslash /2
                        
                      [74].
Figure 1

In an experiment recently conducted at Glasgow University, Heisenberg’s uncertainty principle is manifested in which a light beam that initially possesses no orbital angular momentum was observed. We define Δ ϕ as the range of azimuthal angular positions appearing for a photon in a cross-sectional part of the beam through a sector aperture. Moreover, at the upstream of the said aperture, the beam is observed to be in an = 0 eigen state of orbital angular momentum (i.e., upper panel). Nonetheless, the uncertainty principle states that the limitation in ϕ causes a spread in orbital angular momentum represented as L = for every photon. Additionally, for strait apertures, the connection is portrayed as Δ ϕ Δ L 2 [74].

Figure 2 
                  In the models of atom introduced by Rutherford and Thomson. The later postulated a model that predicted that virtually at small angles, all of the incident alpha particles would be dispersed. Geiger and Rutherford discovered that almost none of the alpha particles were dispersed. However, via a very large angle, a few were found to deflect. Thus, the outcome of the Thomson model was in disagreement with Rutherford’s experiments. Rutherford engaged conservation of energy alongside momentum in the development of a new, and much better model of the atom, and that is the nuclear model [75].
Figure 2

In the models of atom introduced by Rutherford and Thomson. The later postulated a model that predicted that virtually at small angles, all of the incident alpha particles would be dispersed. Geiger and Rutherford discovered that almost none of the alpha particles were dispersed. However, via a very large angle, a few were found to deflect. Thus, the outcome of the Thomson model was in disagreement with Rutherford’s experiments. Rutherford engaged conservation of energy alongside momentum in the development of a new, and much better model of the atom, and that is the nuclear model [75].

Consequently, understanding momentum makes it possible for engineers to have an understanding of various kinds of collisions (Figure 3). As a result, having the knowledge can assist in ensuring that cars are much safer, predict the results observed when two objects bang into each other, or investigate the proof of a traffic accident. One typical example is the utilization of air bags in automobiles. Using air bags in automobiles makes it possible for the effect of the force exerted on an object that is involved in a collision to be minimized. Air bags attain this by lengthening the required time in stopping the momentum of the driver as well as passenger. It is to be acknowledged that whether it is achieved or otherwise, situations involving momentum and impulse can be experienced everywhere.

Figure 3 
                  We have a uni-dimensional inextensible collision in play between two objects. Momentum is conserved, whereas kinetic energy is found not to be conserved. In (a), two objects that are initially having equal mass, move directly toward each other at a uniform speed. In the case of (b), the objects cohere, making way for a perfectly inductile collision. Therefore, the consolidated objects stop and that is the instance revealed in this figure. However, this is untrue for all unyielding collisions [76].
Figure 3

We have a uni-dimensional inextensible collision in play between two objects. Momentum is conserved, whereas kinetic energy is found not to be conserved. In (a), two objects that are initially having equal mass, move directly toward each other at a uniform speed. In the case of (b), the objects cohere, making way for a perfectly inductile collision. Therefore, the consolidated objects stop and that is the instance revealed in this figure. However, this is untrue for all unyielding collisions [76].

Mass can be defined as a measure of the quantum of material present in an object, being directly associated with the type and number of atoms that is found in the object. Mass remains constant with regards to the position of the involved body, its movement or alteration of its shape, except material is removed or added.

In engineering science, mass is used to signal the size of something. Mass which is measured in grams or kilograms likewise referred to as drive is a measure of anxiety, which makes it important in human lives.

Moreover, weight and mass are essential in engineering due to the fact that the greater the mass of any given object is, the greater the force required in accomplishing the same change in motion needed. Thus, for a given object, a larger force initiates a larger change in motion. Figure 4 shows a set-up relaying the relationship between force and acceleration. Besides, Figure 5 portrays some connections among force, weight, and mass.

Figure 4 
                  A graphical depiction of experimental relationship between force and acceleration – Newton’s law [77].
Figure 4

A graphical depiction of experimental relationship between force and acceleration – Newton’s law [77].

Figure 5 
                  A diagrammatic representation of relationship between force, weight, and mass [78].
Figure 5

A diagrammatic representation of relationship between force, weight, and mass [78].

Besides, mass is crucial in science as a result of two major factors affecting the movement of things in space: gravity and inertia. The more mass, an object possesses, the more experience it has regarding both properties. That is the reason heavy things (things with large mass) are difficult to move. Next, energy quantity implies energy quantum found in a certain volume of natural gas expressed in kilowatt hour (kWh).

People walk, ride bicycles, move cars along roads as well as boats through water, cook food on stoves, make ice in freezers, light our homes together with offices, manufacture products, and also send astronauts into space using energy. There are various different forms of energy, inclusive of heat, mechanical, electrical, and so on.

Therefore, energy systems are utilized on daily basis by humans to make life easy. Some of these ways include washing clothes, watching television, taking a shower, heating and lighting the home, working from home on desktop computers or laptop, running appliances, cooking, and so on. Residential uses of energy on a global scale account for almost 40% of total energy utilized.

More applications of energy is found in the phosphagen system (Figure 6) that is active during all-out exercise lasting for about 5–10 s including a 100-m dash, diving, lifting a heavy weight, jumping, dashing up a flight of stairs, or any other scheme that engages a maximal, short burst of power [79].

Figure 6 
                  A diagrammatic representation revealing the Phosphagen system [80].
Figure 6

A diagrammatic representation revealing the Phosphagen system [80].

4 Conclusion

In this article, we clearly purveyed a conspectus investigation carried out on the HD-SOLeqn (1.6). We engage the universal technique, namely, Lie group analysis to which when engaged to solve a DE, it occasions the methodical procedures of generating Lie point symmetries of such equation. In the study, a demonstration of the robust usefulness of the aforementioned technique assisted us in performing a detailed and comprehensive construction of a one-dimensional optimal system of the Lie subalgebras for the nine-dimensional Lie algebra obtained for the equation, which affords us the chance to obtain various more general and robust combinations of the symmetries. Moreover, owing to the relevance of conserved quantities in the study of DEs, we attained diverse associated conservation laws to the HD-SOLeqn using the Lie subalgebras, where various quantities such as conservation of energy are derived. In consequence, this study clearly highlighted the importance of soliton solutions of higher-dimensional NLNPDEs in physics and engineering mathematics and the robust application of the Lie group theory of DEs in proffering solutions to them. Therefore, this research can be beneficial in various fields and in particular in the research area of physical sciences and engineering. Particularly, in an area where further analysis of the result could be of immense usefulness.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

References

[1] Gao XY, Guo YJ, Shan WR. Water-wave symbolic computation for the Earth, Enceladus and Titan: The higher-order Boussinesq–Burgers system, auto-and non-auto-Bäcklund transformations. Appl Math Lett. 2020;104:106170. 10.1016/j.aml.2019.106170Search in Google Scholar

[2] Adeyemo OD, Motsepa T, Khalique CM. A study of the generalized nonlinear advection-diffusion equation arising in engineering sciences. Alex Eng J. 2022;61:185–94. 10.1016/j.aej.2021.04.066Search in Google Scholar

[3] Zhang CR, Tian B, Qu QX, Liu L, Tian HY. Vector bright solitons and their interactions of the couple Fokas-Lenells system in a birefringent optical fiber. Z Angew Math Phys. 2020;71:1–19. 10.1007/s00033-019-1225-9Search in Google Scholar

[4] Adeyemo OD, Khalique CM, Gasimov YS, Villecco F. Variational and non-variational approaches with Lie algebra of a generalized (3+1)-dimensional nonlinear potential Yu-Toda-Sasa-Fukuyama equation in engineering and physics. Alex Eng J. 2023;63:17–43. 10.1016/j.aej.2022.07.024Search in Google Scholar

[5] Du XX, Tian B, Qu QX, Yuan YQ, Zhao XH. Lie group analysis, solitons, self-adjointness and conservation laws of the modified Zakharov-Kuznetsov equation in an electron-positron-ion magnetoplasma. Chaos Solitons Fract. 2020;134:109709. 10.1016/j.chaos.2020.109709Search in Google Scholar

[6] Khalique CM, Adeyemo OD. A study of (3+1)-dimensional generalized Korteweg-de Vries-Zakharov-Kuznetsov equation via Lie symmetry approach. Results Phys. 2020;18:103197. 10.1016/j.rinp.2020.103197Search in Google Scholar

[7] Younis M. Optical solitons in (n+1) dimensions with Kerr and power law nonlinearities. Mod Phys Lett B. 2017;31:1750186. 10.1142/S021798491750186XSearch in Google Scholar

[8] Liu MM, Yu JP, Ma WX, Khalique CM, Sun YL. Dynamic analysis of lump solutions based on the dimensionally reduced generalized Hirota bilinear KP-Boussinesq equation. Mod Phys Lett B. 2023;37:2250203. 10.1142/S0217984922502037Search in Google Scholar

[9] Bilal M, Seadawy AR, Younis M. Dispersive of propagation wave solutions to unidirectional shallow water wave Dullin-Gottwald-Holm system and modulation instability analysis. Math Methods Appl Sci. 2021;44:4094–104. 10.1002/mma.7013Search in Google Scholar

[10] Simbanefayi I, Gandarias ML, Khalique CM. Travelling wave solutions, symmetry reductions and conserved vectors of a generalized hyper-elastic rod wave equation. Partial Differ Equ Appl Math. 2023;7:10050110.1016/j.padiff.2023.100501Search in Google Scholar

[11] Younis M, Ali S, Rizvi STR, Tantawy M, Tariq KU. Investigation of solitons and mixed lump wave solutions with (3+1)-dimensional potential-YTSF equation. Commun Nonlinear Sci Numer Simul. 2021;94:105544. 10.1016/j.cnsns.2020.105544Search in Google Scholar

[12] Benoudina N, Zhang Y, Khalique CM. Lie symmetry analysis, optimal system, new solitary wave solutions and conservation laws of the Pavlov equation. Commun Nonlinear Sci Numer Simulat. 2021;94:105560. 10.1016/j.cnsns.2020.105560Search in Google Scholar

[13] Kudryashov NA. Analytical theory of nonlinear differential equations. Moskow - Igevsk: Institute of Computer Investigations; 2004. Search in Google Scholar

[14] Bluman GW, Cheviakov AF, Anco SC. Applications of symmetry methods to partial differential equations. New York: Springer; 2010. 10.1007/978-0-387-68028-6Search in Google Scholar

[15] Kumar S, Kumar D, Kumar A. Lie symmetry analysis for obtaining the abundant exact solutions, optimal system and dynamics of solitons for a higher-dimensional Fokas equation. Chaos Solit Fractals. 2021;142:110507. 10.1016/j.chaos.2020.110507Search in Google Scholar

[16] 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:065204. 10.1088/1402-4896/aafc13Search in Google Scholar

[17] Olver PJ. Applications of Lie groups to differential equations. In: Graduate Texts in Mathematics. Vol. 107. 2nd edition. Berlin: Springer-Verlag; 1993. 10.1007/978-1-4612-4350-2Search in Google Scholar

[18] Ovsiannikov LV. Group analysis of differential equations. New York, USA: Academic Press; 1982. 10.1016/B978-0-12-531680-4.50012-5Search in Google Scholar

[19] Zhang L, Kwizera S, Khalique CM. A study of a new generalized Burgers’ equation: symmetry solutions and conservation laws. Adv Math Models Appl. 2023;8:2. Search in Google Scholar

[20] Adeyemo OD, Zhang L, Khalique CM. Optimal solutions of Lie subalgebra, dynamical system, travelling wave solutions and conserved currents of (3+1)-dimensional generalized Zakharov-Kuznetsov equation type I. Eur Phys J Plus. 2022;137:95410.1140/epjp/s13360-022-03100-zSearch in Google Scholar

[21] Anco SC, Bluman GW. Direct construction method for conservation laws of partial differential equations. Part I: Examples of conservation law classifications. Eur J Appl Math. 2002;13:545–66. 10.1017/S095679250100465XSearch in Google Scholar

[22] Feng L, Tian S, Zhang T, Zhou J. Lie symmetries, conservation laws and analytical solutions for two-component integrable equations. Chinese J Phys. 2017;55:996–1010. 10.1016/j.cjph.2017.03.008Search in Google Scholar

[23] Yu J, Wang D, Sun Y, S, Wu. Modified method of simplest equation for obtaining exact solutions of the Zakharov-Kuznetsov equation, the modified Zakharov-Kuznetsov equation, and their generalized forms. Nonlinear Dyn. 2016;85:2449–65. 10.1007/s11071-016-2837-7Search in Google Scholar

[24] Zhang Y, Ye R, Ma WX. Binary Darboux transformation and soliton solutions for the coupled complex modified Korteweg-de Vries equations. Math Meth Appl Sci. 2020;43:613–27. 10.1002/mma.5914Search in Google Scholar

[25] Adem AR, Yildrim Y, Sar EY. Complexiton solutions and soliton solutions: (2+1)-dimensional Date-Jimbo-Kashiwara-Miwa equation. Pramana. 2019;92:1–12. 10.1007/s12043-018-1707-xSearch in Google Scholar

[26] Zhang L, Khalique CM. Classification and bifurcation of a class of second-order ODEs and its application to nonlinear PDEs. Discrete and Continuous dynamical systems Series S. 2018;11:777–90. 10.3934/dcdss.2018048Search in Google Scholar

[27] Weiss J, Tabor M, Carnevale G. The Painlevé property and a partial differential equations with an essential singularity. Phys Lett A. 1985;109:205–8. 10.1016/0375-9601(85)90303-2Search in Google Scholar

[28] Chun C, Sakthivel R. Homotopy perturbation technique for solving two point boundary value problems-comparison with other methods. Comput Phys Commun. 2010;181:1021–4. 10.1016/j.cpc.2010.02.007Search in Google Scholar

[29] Wazwaz AM. Traveling wave solution to (2+1)-dimensional nonlinear evolution equations. J Nat Sci Math. 2007;1:1–13. Search in Google Scholar

[30] Darvishi MT, Najafi M. A modification of extended homoclinic test approach to solve the (3+1)-dimensional potential-YTSF equation. Chin Phys Lett. 2011;28:040202. 10.1088/0256-307X/28/4/040202Search in Google Scholar

[31] Salas AH, Gomez CA. Application of the Cole–Hopf transformation for finding exact solutions to several forms of the seventh-order KdV equation. Math Probl Eng. 2010;2010. 10.1155/2010/194329Search in Google Scholar

[32] Wazwaz AM. Partial differential equations. Boca Raton, Florida, USA: CRC Press; 2002. Search in Google Scholar

[33] Gu CH. Soliton theory and its application. Zhejiang, China: Zhejiang Science and Technology Press; 1990. Search in Google Scholar

[34] Zhou Y, Wang M, Wang Y. Periodic wave solutions to a coupled KdV equations with variable coefficients. Phys Lett A. 2003;308:31–6. 10.1016/S0375-9601(02)01775-9Search in Google Scholar

[35] Zeng X, Wang DS. A generalized extended rational expansion method and its application to (1+1)-dimensional dispersive long wave equation. Appl Math Comput. 2009;212:296–304. 10.1016/j.amc.2009.02.020Search in Google Scholar

[36] Kudryashov NA, Loguinova NB, Extended simplest equation method for nonlinear differential equations. Appl Math Comput. 2008;205:396–402. 10.1016/j.amc.2008.08.019Search in Google Scholar

[37] Kudryashov NA. Simplest equation method to look for exact solutions of nonlinear differential equations. Chaos Solitons Fract. 2005;24:1217–3110.1016/j.chaos.2004.09.109Search in Google Scholar

[38] Hirota R. The direct method in soliton theory. Cambridge, UK: Cambridge University Press; 2004. 10.1017/CBO9780511543043Search in Google Scholar

[39] Matveev VB, Salle MA. Darboux transformations and solitons. New York, USA: Springer; 1991. 10.1007/978-3-662-00922-2Search in Google Scholar

[40] Wazwaz AM. The tanh method for generalized forms of nonlinear heat conduction and Burgers-Fisher equations. Appl Math Comput. 2005;169:321–38. 10.1016/j.amc.2004.09.054Search in Google Scholar

[41] Wang M, Li X, Zhang J. The (G′∕G)-expansion method and travelling wave solutions for linear evolution equations in mathematical physics. Phys Lett A. 2005;24:1257–68. Search in Google Scholar

[42] Chen Y, Z Yan. New exact solutions of (2+1)-dimensional Gardner equation via the new sine-Gordon equation expansion method. Chaos Solitons Fract. 2005;26:399–406. 10.1016/j.chaos.2005.01.004Search in Google Scholar

[43] Osman MS. One-soliton shaping and inelastic collision between double solitons in the fifth-order variable-coefficient Sawada-Kotera equation. Nonlinear Dynam. 2019;96:1491–6. 10.1007/s11071-019-04866-1Search in Google Scholar

[44] He JH, Wu XH. Exp-function method for nonlinear wave equations. Chaos Solitons Fract. 2006;30:700–8. 10.1016/j.chaos.2006.03.020Search in Google Scholar

[45] Date Y, Jimbo M, Kashiwara M, Miwa T. Operator approach of the Kadomtsev-Petviashvili equation - Transformation groups for soliton equations III. JPSJ. 1981;50:3806–12. 10.1143/JPSJ.50.3806Search in Google Scholar

[46] Kuo CK, Ma WX. An effective approach to constructing novel KP-like equations. Waves random complex media. 2020;32:629–40. 10.1080/17455030.2020.1792580Search in Google Scholar

[47] Ma WX, Fan E. Linear superposition principle applying to Hirota bilinear equations. Comput Math Appl. 2011;61:950–9. 10.1016/j.camwa.2010.12.043Search in Google Scholar

[48] Wazwaz AM. Multiple-soliton solutions for a (3+1)-dimensional generalized KP equation. Commun Nonlinear Sci Numer Simul. 2012;17:491–5. 10.1016/j.cnsns.2011.05.025Search in Google Scholar

[49] Ma WX. Lump solutions to the Kadomtsev-Petviashvili equation. Phys Lett A. 2015;379:1975–8. 10.1016/j.physleta.2015.06.061Search in Google Scholar

[50] Zhao Z, Han B. Lump solutions of a (3+1)-dimensional B-type KP equation and its dimensionally reduced equations. Anal Math Phys. 2017;9:119–30. 10.1007/s13324-017-0185-5Search in Google Scholar

[51] Simbanefayi I, Khalique CM. Group invariant solutions and conserved quantities of a (3+1)-dimensional generalized Kadomtsev-Petviashvili equation. Mathematics. 2020;8:1012. 10.3390/math8061012Search in Google Scholar

[52] Asaad MG. Soliton solution in (3+1)-dimensions. The 2015 International Academic Research Conference. September University of Nevada, USA. 2015. Vol. 59. p. 15–18. Search in Google Scholar

[53] Geng X. Algebraic-geometrical solutions of some multidimensional nonlinear evolution equations. J Phys Math Gen. 2003;36:2289. 10.1088/0305-4470/36/9/307Search in Google Scholar

[54] Liu J, Zhang Y. Construction of lump soliton and mixed lump stripe solutions of (3+1)-dimensional soliton equation. Results Phys. 2018;10:94–8. 10.1016/j.rinp.2018.05.022Search in Google Scholar

[55] Liu J, Wu P, Zhang Y. New periodic wave solutions of (3+1)-dimensional soliton equation. Therm Sci. 2017;21:169–76. 10.2298/TSCI17S1169LSearch in Google Scholar

[56] Geng X, Ma Y. N-soliton solution and its Wronskian form of a (3+1)-dimensional nonlinear evolution equation. Phys Lett A. 2007;369:285–9. 10.1016/j.physleta.2007.04.099Search in Google Scholar

[57] Jian-Ping W, Xian-Guo G. Grammian determinant solution and Pfaffianization for a (3+1)-dimensional soliton equation. Commun Theor Phys. 2009;52:791. 10.1088/0253-6102/52/5/05Search in Google Scholar

[58] Jian-Ping W. A bilinear Bäcklund transformation and explicit solutions for a (3+1)-dimensional soliton equation. Chinese Phys Lett. 2008;25:4192. 10.1088/0256-307X/25/12/002Search in Google Scholar

[59] Wang X, Wei J, Geng X. Rational solutions for a (3+1)-dimensional nonlinear evolution equation. Commun Nonlinear Sci Numer Simul. 2020;83:105116. 10.1016/j.cnsns.2019.105116Search in Google Scholar

[60] Wang X, Wei J. Antidark solitons and soliton molecules in a (3+1)-dimensional nonlinear evolution equation. Nonlinear Dyn. 2020;102:363–77. 10.1007/s11071-020-05926-7Search in Google Scholar

[61] Khalique CM, Adeyemo OD. Soliton solutions, travelling wave solutions and conserved quantities for a three-dimensional soliton equation in plasma physics. Commun Theor Phys. 2021;73:125003. 10.1088/1572-9494/ac27a1Search in Google Scholar

[62] Adeyemo OD, Khalique CM, Dynamical soliton wave structures of one-dimensional Lie subalgebras via group-invariant solutions of a higher-dimensional soliton equation in ocean physics and mechatronics engineering, Commun Appl Math Comput. 2022;4:1531–82. 10.1007/s42967-022-00195-0Search in Google Scholar

[63] Ibragimov NH, CRC handbook of lie group analysis of differential equations. Vol. 1–3. Boca Raton, Florida: CRC Press; 1994–1996. Search in Google Scholar

[64] Hu X, Li Y, Chen Y. A direct algorithm of one-dimensional optimal system for the group invariant solutions. J Math Phys. 2015;56:053504. 10.1063/1.4921229Search in Google Scholar

[65] Ibragimov NH. A new conservation theorem. J Math Anal Appl. 2007;333:311–28. 10.1016/j.jmaa.2006.10.078Search in Google Scholar

[66] Ibragimov NH. Integrating factors, adjoint equations and Lagrangians. J Math Anal Appl. 2006;318:742–57. 10.1016/j.jmaa.2005.11.012Search in Google Scholar

[67] Márquez AP, Bruzón MS. Lie point symmetries, traveling wave solutions and conservation laws of a non-linear viscoelastic wave equation. Mathematics. 2021;9:2131. 10.3390/math9172131Search in Google Scholar

[68] Adeyemo OD, Khalique CM. Lie group theory, stability analysis with dispersion property, new soliton solutions and conserved quantities of 3D generalized nonlinear wave equation in liquid containing gas bubbles with applications in mechanics of fluids, biomedical sciences and cell biology, Commun Nonlinear Sci Numer Simul. 2023;123:107261. 10.1016/j.cnsns.2023.107261Search in Google Scholar

[69] https://www.britannica.com/science/conservation-law. Search in Google Scholar

[70] https://en.wikipedia.org/wiki/Conservationlaw. Search in Google Scholar

[71] Khalique CM, Adeyemo OD. Closed-form solutions and conserved vectors of a generalized (3. 1)-dimensional breaking soliton equation of engineering and nonlinear science. Mathematics. 2020;8:1692. 10.3390/math8101692Search in Google Scholar

[72] Anco SC. Symmetry properties of conservation laws. Int J Modern Phys B. 2016;30:1640003. 10.1142/S0217979216400038Search in Google Scholar

[73] Anco SC, Kara A. Symmetry invariance of conservation laws. Euro J Appl Math. 2018;29:78–117. 10.1017/S0956792517000055Search in Google Scholar

[74] Physics today: Light’s Orbital Angular Momentum. [Online]. Available: https://physicstoday.scitation.org/doi/10.1063/1.1768672. Search in Google Scholar

[75] University Physics : Types of Collisions. [Online]. Available: https://courses.lumenlearning.com/suny-osuniversityphysics/chapter/9-4-types-of-collisions/. Search in Google Scholar

[76] Elastic and Inelastic Collisions. [Online]. Available: https://www.texasgateway.org/resource/83-elastic-and-inelastic-collisions. Search in Google Scholar

[77] Investigate the Relationship Between Force and Acceleration. [Online]. Available. https://classnotes.gidemy.com/topic/investigate-the-relationship-between-force-and-acceleration/. Search in Google Scholar

[78] Understanding the Relationship Between Mass and Weight. [Online]. Available: https://www.arborsci.com/blogs/cool/understanding-the-relationship-between-mass-and-weight. Search in Google Scholar

[79] https://fittrakker.com/scale-of-energy-pools-used-in-energy-physiology-zones/. Search in Google Scholar

[80] Kiika D. “The Basics of Energy Production”: The Phosphagen System. [Online]. Available: https://thesportsedu.com/the-phosphagen-system/. [ Accessed 27 March 2023]. Search in Google Scholar

Received: 2023-10-07
Revised: 2023-10-23
Accepted: 2023-11-12
Published Online: 2024-03-26

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Numerical study of flow and heat transfer in the channel of panel-type radiator with semi-detached inclined trapezoidal wing vortex generators
  3. Homogeneous–heterogeneous reactions in the colloidal investigation of Casson fluid
  4. High-speed mid-infrared Mach–Zehnder electro-optical modulators in lithium niobate thin film on sapphire
  5. Numerical analysis of dengue transmission model using Caputo–Fabrizio fractional derivative
  6. Mononuclear nanofluids undergoing convective heating across a stretching sheet and undergoing MHD flow in three dimensions: Potential industrial applications
  7. Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface
  8. The electrically conducting water-based nanofluid flow containing titanium and aluminum alloys over a rotating disk surface with nonlinear thermal radiation: A numerical analysis
  9. Growth, characterization, and anti-bacterial activity of l-methionine supplemented with sulphamic acid single crystals
  10. A numerical analysis of the blood-based Casson hybrid nanofluid flow past a convectively heated surface embedded in a porous medium
  11. Optoelectronic–thermomagnetic effect of a microelongated non-local rotating semiconductor heated by pulsed laser with varying thermal conductivity
  12. Thermal proficiency of magnetized and radiative cross-ternary hybrid nanofluid flow induced by a vertical cylinder
  13. Enhanced heat transfer and fluid motion in 3D nanofluid with anisotropic slip and magnetic field
  14. Numerical analysis of thermophoretic particle deposition on 3D Casson nanofluid: Artificial neural networks-based Levenberg–Marquardt algorithm
  15. Analyzing fuzzy fractional Degasperis–Procesi and Camassa–Holm equations with the Atangana–Baleanu operator
  16. Bayesian estimation of equipment reliability with normal-type life distribution based on multiple batch tests
  17. Chaotic control problem of BEC system based on Hartree–Fock mean field theory
  18. Optimized framework numerical solution for swirling hybrid nanofluid flow with silver/gold nanoparticles on a stretching cylinder with heat source/sink and reactive agents
  19. Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations
  20. Convective flow of a magnetohydrodynamic second-grade fluid past a stretching surface with Cattaneo–Christov heat and mass flux model
  21. Analysis of the heat transfer enhancement in water-based micropolar hybrid nanofluid flow over a vertical flat surface
  22. Microscopic seepage simulation of gas and water in shale pores and slits based on VOF
  23. Model of conversion of flow from confined to unconfined aquifers with stochastic approach
  24. Study of fractional variable-order lymphatic filariasis infection model
  25. Soliton, quasi-soliton, and their interaction solutions of a nonlinear (2 + 1)-dimensional ZK–mZK–BBM equation for gravity waves
  26. Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering
  27. Nonlinear fractional-order differential equations: New closed-form traveling-wave solutions
  28. Sixth-kind Chebyshev polynomials technique to numerically treat the dissipative viscoelastic fluid flow in the rheology of Cattaneo–Christov model
  29. Some transforms, Riemann–Liouville fractional operators, and applications of newly extended M–L (p, s, k) function
  30. Magnetohydrodynamic water-based hybrid nanofluid flow comprising diamond and copper nanoparticles on a stretching sheet with slips constraints
  31. Super-resolution reconstruction method of the optical synthetic aperture image using generative adversarial network
  32. A two-stage framework for predicting the remaining useful life of bearings
  33. Influence of variable fluid properties on mixed convective Darcy–Forchheimer flow relation over a surface with Soret and Dufour spectacle
  34. Inclined surface mixed convection flow of viscous fluid with porous medium and Soret effects
  35. Exact solutions to vorticity of the fractional nonuniform Poiseuille flows
  36. In silico modified UV spectrophotometric approaches to resolve overlapped spectra for quality control of rosuvastatin and teneligliptin formulation
  37. Numerical simulations for fractional Hirota–Satsuma coupled Korteweg–de Vries systems
  38. Substituent effect on the electronic and optical properties of newly designed pyrrole derivatives using density functional theory
  39. A comparative analysis of shielding effectiveness in glass and concrete containers
  40. Numerical analysis of the MHD Williamson nanofluid flow over a nonlinear stretching sheet through a Darcy porous medium: Modeling and simulation
  41. Analytical and numerical investigation for viscoelastic fluid with heat transfer analysis during rollover-web coating phenomena
  42. Influence of variable viscosity on existing sheet thickness in the calendering of non-isothermal viscoelastic materials
  43. Analysis of nonlinear fractional-order Fisher equation using two reliable techniques
  44. Comparison of plan quality and robustness using VMAT and IMRT for breast cancer
  45. Radiative nanofluid flow over a slender stretching Riga plate under the impact of exponential heat source/sink
  46. Numerical investigation of acoustic streaming vortices in cylindrical tube arrays
  47. Numerical study of blood-based MHD tangent hyperbolic hybrid nanofluid flow over a permeable stretching sheet with variable thermal conductivity and cross-diffusion
  48. Fractional view analytical analysis of generalized regularized long wave equation
  49. Dynamic simulation of non-Newtonian boundary layer flow: An enhanced exponential time integrator approach with spatially and temporally variable heat sources
  50. Inclined magnetized infinite shear rate viscosity of non-Newtonian tetra hybrid nanofluid in stenosed artery with non-uniform heat sink/source
  51. Estimation of monotone α-quantile of past lifetime function with application
  52. Numerical simulation for the slip impacts on the radiative nanofluid flow over a stretched surface with nonuniform heat generation and viscous dissipation
  53. Study of fractional telegraph equation via Shehu homotopy perturbation method
  54. An investigation into the impact of thermal radiation and chemical reactions on the flow through porous media of a Casson hybrid nanofluid including unstable mixed convection with stretched sheet in the presence of thermophoresis and Brownian motion
  55. Establishing breather and N-soliton solutions for conformable Klein–Gordon equation
  56. An electro-optic half subtractor from a silicon-based hybrid surface plasmon polariton waveguide
  57. CFD analysis of particle shape and Reynolds number on heat transfer characteristics of nanofluid in heated tube
  58. Abundant exact traveling wave solutions and modulation instability analysis to the generalized Hirota–Satsuma–Ito equation
  59. A short report on a probability-based interpretation of quantum mechanics
  60. Study on cavitation and pulsation characteristics of a novel rotor-radial groove hydrodynamic cavitation reactor
  61. Optimizing heat transport in a permeable cavity with an isothermal solid block: Influence of nanoparticles volume fraction and wall velocity ratio
  62. Linear instability of the vertical throughflow in a porous layer saturated by a power-law fluid with variable gravity effect
  63. Thermal analysis of generalized Cattaneo–Christov theories in Burgers nanofluid in the presence of thermo-diffusion effects and variable thermal conductivity
  64. A new benchmark for camouflaged object detection: RGB-D camouflaged object detection dataset
  65. Effect of electron temperature and concentration on production of hydroxyl radical and nitric oxide in atmospheric pressure low-temperature helium plasma jet: Swarm analysis and global model investigation
  66. Double diffusion convection of Maxwell–Cattaneo fluids in a vertical slot
  67. Thermal analysis of extended surfaces using deep neural networks
  68. Steady-state thermodynamic process in multilayered heterogeneous cylinder
  69. Multiresponse optimisation and process capability analysis of chemical vapour jet machining for the acrylonitrile butadiene styrene polymer: Unveiling the morphology
  70. Modeling monkeypox virus transmission: Stability analysis and comparison of analytical techniques
  71. Fourier spectral method for the fractional-in-space coupled Whitham–Broer–Kaup equations on unbounded domain
  72. The chaotic behavior and traveling wave solutions of the conformable extended Korteweg–de-Vries model
  73. Research on optimization of combustor liner structure based on arc-shaped slot hole
  74. Construction of M-shaped solitons for a modified regularized long-wave equation via Hirota's bilinear method
  75. Effectiveness of microwave ablation using two simultaneous antennas for liver malignancy treatment
  76. Discussion on optical solitons, sensitivity and qualitative analysis to a fractional model of ion sound and Langmuir waves with Atangana Baleanu derivatives
  77. Reliability of two-dimensional steady magnetized Jeffery fluid over shrinking sheet with chemical effect
  78. Generalized model of thermoelasticity associated with fractional time-derivative operators and its applications to non-simple elastic materials
  79. Migration of two rigid spheres translating within an infinite couple stress fluid under the impact of magnetic field
  80. A comparative investigation of neutron and gamma radiation interaction properties of zircaloy-2 and zircaloy-4 with consideration of mechanical properties
  81. New optical stochastic solutions for the Schrödinger equation with multiplicative Wiener process/random variable coefficients using two different methods
  82. Physical aspects of quantile residual lifetime sequence
  83. Synthesis, structure, IV characteristics, and optical properties of chromium oxide thin films for optoelectronic applications
  84. Smart mathematically filtered UV spectroscopic methods for quality assurance of rosuvastatin and valsartan from formulation
  85. A novel investigation into time-fractional multi-dimensional Navier–Stokes equations within Aboodh transform
  86. Homotopic dynamic solution of hydrodynamic nonlinear natural convection containing superhydrophobicity and isothermally heated parallel plate with hybrid nanoparticles
  87. A novel tetra hybrid bio-nanofluid model with stenosed artery
  88. Propagation of traveling wave solution of the strain wave equation in microcrystalline materials
  89. Innovative analysis to the time-fractional q-deformed tanh-Gordon equation via modified double Laplace transform method
  90. A new investigation of the extended Sakovich equation for abundant soliton solution in industrial engineering via two efficient techniques
  91. New soliton solutions of the conformable time fractional Drinfel'd–Sokolov–Wilson equation based on the complete discriminant system method
  92. Irradiation of hydrophilic acrylic intraocular lenses by a 365 nm UV lamp
  93. Inflation and the principle of equivalence
  94. The use of a supercontinuum light source for the characterization of passive fiber optic components
  95. Optical solitons to the fractional Kundu–Mukherjee–Naskar equation with time-dependent coefficients
  96. A promising photocathode for green hydrogen generation from sanitation water without external sacrificing agent: silver-silver oxide/poly(1H-pyrrole) dendritic nanocomposite seeded on poly-1H pyrrole film
  97. Photon balance in the fiber laser model
  98. Propagation of optical spatial solitons in nematic liquid crystals with quadruple power law of nonlinearity appears in fluid mechanics
  99. Theoretical investigation and sensitivity analysis of non-Newtonian fluid during roll coating process by response surface methodology
  100. Utilizing slip conditions on transport phenomena of heat energy with dust and tiny nanoparticles over a wedge
  101. Bismuthyl chloride/poly(m-toluidine) nanocomposite seeded on poly-1H pyrrole: Photocathode for green hydrogen generation
  102. Infrared thermography based fault diagnosis of diesel engines using convolutional neural network and image enhancement
  103. On some solitary wave solutions of the Estevez--Mansfield--Clarkson equation with conformable fractional derivatives in time
  104. Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres
  105. Review Article
  106. Transformer-based intelligent fault diagnosis methods of mechanical equipment: A survey
  107. Special Issue on Predicting pattern alterations in nature - Part II
  108. A comparative study of Bagley–Torvik equation under nonsingular kernel derivatives using Weeks method
  109. On the existence and numerical simulation of Cholera epidemic model
  110. Numerical solutions of generalized Atangana–Baleanu time-fractional FitzHugh–Nagumo equation using cubic B-spline functions
  111. Dynamic properties of the multimalware attacks in wireless sensor networks: Fractional derivative analysis of wireless sensor networks
  112. Prediction of COVID-19 spread with models in different patterns: A case study of Russia
  113. Study of chronic myeloid leukemia with T-cell under fractal-fractional order model
  114. Accumulation process in the environment for a generalized mass transport system
  115. Analysis of a generalized proportional fractional stochastic differential equation incorporating Carathéodory's approximation and applications
  116. Special Issue on Nanomaterial utilization and structural optimization - Part II
  117. Numerical study on flow and heat transfer performance of a spiral-wound heat exchanger for natural gas
  118. Study of ultrasonic influence on heat transfer and resistance performance of round tube with twisted belt
  119. Numerical study on bionic airfoil fins used in printed circuit plate heat exchanger
  120. Improving heat transfer efficiency via optimization and sensitivity assessment in hybrid nanofluid flow with variable magnetism using the Yamada–Ota model
  121. Special Issue on Nanofluids: Synthesis, Characterization, and Applications
  122. Exact solutions of a class of generalized nanofluidic models
  123. Stability enhancement of Al2O3, ZnO, and TiO2 binary nanofluids for heat transfer applications
  124. Thermal transport energy performance on tangent hyperbolic hybrid nanofluids and their implementation in concentrated solar aircraft wings
  125. Studying nonlinear vibration analysis of nanoelectro-mechanical resonators via analytical computational method
  126. Numerical analysis of non-linear radiative Casson fluids containing CNTs having length and radius over permeable moving plate
  127. Two-phase numerical simulation of thermal and solutal transport exploration of a non-Newtonian nanomaterial flow past a stretching surface with chemical reaction
  128. Natural convection and flow patterns of Cu–water nanofluids in hexagonal cavity: A novel thermal case study
  129. Solitonic solutions and study of nonlinear wave dynamics in a Murnaghan hyperelastic circular pipe
  130. Comparative study of couple stress fluid flow using OHAM and NIM
  131. Utilization of OHAM to investigate entropy generation with a temperature-dependent thermal conductivity model in hybrid nanofluid using the radiation phenomenon
  132. Slip effects on magnetized radiatively hybridized ferrofluid flow with acute magnetic force over shrinking/stretching surface
  133. Significance of 3D rectangular closed domain filled with charged particles and nanoparticles engaging finite element methodology
  134. Robustness and dynamical features of fractional difference spacecraft model with Mittag–Leffler stability
  135. Characterizing magnetohydrodynamic effects on developed nanofluid flow in an obstructed vertical duct under constant pressure gradient
  136. Study on dynamic and static tensile and puncture-resistant mechanical properties of impregnated STF multi-dimensional structure Kevlar fiber reinforced composites
  137. Thermosolutal Marangoni convective flow of MHD tangent hyperbolic hybrid nanofluids with elastic deformation and heat source
  138. Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks
  139. Single-channel cooling system design by using perforated porous insert and modeling with POD for double conductive panel
  140. Special Issue on Fundamental Physics from Atoms to Cosmos - Part I
  141. Pulsed excitation of a quantum oscillator: A model accounting for damping
  142. Review of recent analytical advances in the spectroscopy of hydrogenic lines in plasmas
  143. Heavy mesons mass spectroscopy under a spin-dependent Cornell potential within the framework of the spinless Salpeter equation
  144. Coherent manipulation of bright and dark solitons of reflection and transmission pulses through sodium atomic medium
  145. Effect of the gravitational field strength on the rate of chemical reactions
  146. The kinetic relativity theory – hiding in plain sight
  147. Special Issue on Advanced Energy Materials - Part III
  148. Eco-friendly graphitic carbon nitride–poly(1H pyrrole) nanocomposite: A photocathode for green hydrogen production, paving the way for commercial applications
Downloaded on 10.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2023-0155/html
Scroll to top button