Integrability Properties of the Slepyan–Palmov Model Arising in the Slepyan–Palmov Medium

: This study investigates the Slepyan–Palmov (SP) model, which describes plane longitudinal waves propagating within a medium comprising a carrier medium and nonlinear oscillators. The primary objective is to analyze the integrability properties of this model. The research entails two key aspects. Firstly, the study explores the group invariant solution by utilizing reductions in symmetry subalgebras based on the optimal system. Secondly, the conservation laws are studied using the homotopy operator, which offers advantages over the conventional multiplier approach, especially when arbitrary functions are absent from both the equation and characteristics. This method proves advantageous in handling complex multipliers and yields signiﬁcant outcomes.


Introduction
The Cauchy stress principle is a fundamental assumption in classical continuum mechanics which establishes a key equivalence between the action of all internal forces applied to an elementary area and the action of their resultant force applied at the center of the area.This principle forms the backbone of classical continuum mechanics and plays a crucial role in understanding the behavior of continuous materials under various loads and deformations.However, as engineering and scientific investigations have advanced, it has become evident that the simple stress principle embodied by Cauchy's hypothesis might not fully capture the complexity of certain materials and deformation scenarios.In more general cases, the action of an arbitrary system of forces on a deformable solid cannot be solely represented by stresses alone.In addition to stresses, moment stresses emerge, leading to the formation of asymmetric tensors.These moment stresses significantly influence the mechanical response of materials, especially in situations where classical continuum mechanics fall short of describing real-world phenomena accurately.
To adequately account for these complex behaviors and to broaden the understanding of material responses, the introduction of additional degrees of freedom in the medium becomes necessary.This involves considering a physically infinitesimal volume (over which the properties of the medium are averaged) not as a simple material point, but as a more intricate object with new degrees of freedom.This recognition of additional degrees of freedom allows for the inclusion of a microstructure within the material, such as graininess, fibrousness, or cellular structures present in real-world materials.By embracing this expanded perspective, the range of properties that can be modeled within a continuous medium expands considerably.The presence of internal microstructure and the consideration of additional degrees of freedom enable a more comprehensive representation of the mechanical response of materials, particularly in cases where classical continuum mechanics might be limited in its accuracy.Both the classical continuum model and the generalized continuum model are extensively utilized in modern deformable solid mechanics [1][2][3].
In addition to the classical continuum model and the generalized continuum model, there are other models used in deformable solid mechanics.One example is the Cosserat continuum [4], which considers internal rotational degrees of freedom.This theory was originally published in Russian and it has been subject to comments and discussions [5].Gradient models, like the Leroux continuum [6], Jeremillo model, and Tupin model [7], also exist.These models incorporate gradients of certain fields to handle materials with spatial variations in properties, proving useful in describing complex materials.
Furthermore, there are models that account for media with oscillatory degrees of freedom.These models address the presence of oscillations or vibrations within the material.The works of Slepyan [8] in one-dimensional systems and Palmov [9] in three-dimensional systems have significantly contributed to this area of research.These diverse models enrich our understanding of material behaviors and offer valuable insights into various physical phenomena in deformable solids.
The SP model combines a linearly elastic carrier medium with non-interacting oscillators (elastic or viscoelastic) suspended at each point.It postulates that the dynamic behavior of the carrier model is described by the Lame equations and the oscillators fixed within it have continuously distributed eigenfrequencies.
The following equations govern the dynamics of the SP model where • ρ represents the mass density of the carrier medium.• and µ are the Lame elastic moduli characterizing the carrier medium's elasticity.• v denotes the displacement vector of points within the carrier medium.• K stands for the intensity of the external body force acting on the medium.• w q represents the absolute displacement vector of the oscillator mass with respect to its equilibrium position.• Q q denotes the external force applied to the mass of the oscillator.

•
The quantity m(q)dq corresponds to the mass of all oscillators with eigenfrequencies lying within the interval (q, q + dq) multiplied by a unit volume.• m = ∞ 0 m(q)dq is the total mass density of all oscillators fixed to the carrier medium.• c(q) = q 2 m(q) represents the static stiffness of the oscillator suspension.
The quantity R(∂/∂κ) characterizes the energy dissipation in the oscillator suspension.An interesting characteristic of this model is that even with the low damping of oscillators, the spatial attention of vibrations in the medium is finite [10].As a result, the model finds effective application in calculating the vibrations of aircraft, rockets, and space technology objects, as well as submarines.It proves to be a valuable tool in analyzing and predicting the behavior of these structures under various vibrational conditions.
In the context of an isolated environment with K = 0 and Q q = 0, we focus on a onedimensional version of the system (1).To account for nonlinearity attributed to the carrier medium, we introduce nonlinear terms into the system.The resulting one-dimensional system is expressed as where F is the coefficient characterizing the nonlinearity of the carrier medium.
Under the assumption that the absolute displacement of the oscillator does not depend on its eigenfrequency, we can express it as The system (2) is simplified to a single equation representing the longitudinal displacement of the carrier medium, denoted as v y given by We define the dimensionless parameters as By introducing v 0 as the maximum displacement, within which the deformation of the carrier medium remains elastic, Equation ( 3) is modified to take the following form Considering the nonlinearity of the medium in the absence of dissipation (R = 0, N 1 = 0), Equation (3) can be expressed as follows In this study, we conduct a Lie symmetry analysis of the SP model [11] in the given mathematical form Lie theorists like Ovsyannikov [12], Ibragimov [13], Bluman [14], Olver [15], Hydon [16], Stephani [17], and others [18,19] have played a vital role in connecting group structures to differential equations.This link has been crucial in discovering invariant solutions for nonlinear partial differential equations and gaining deeper insight into their behavior.A significant area of study has been the group analysis of the SP model, where Lie theory techniques have led to the identification of invariant solutions and conservation laws, highlighting essential underlying symmetries and dynamic properties.This study provides group invariant solutions that remain unchanged under a specific Lie symmetry generator, adding novelty to our research as such solutions are not present in the existing literature.Additionally, the conservation laws outlined in this study are reported for the first time, contributing to the enrichment of the SP model's applications from a divergence perspective.Our study is specifically centered on invariant solutions and nonlocal conservation laws, which represents a limitation in terms of the methods we have employed.References [20][21][22] represent valuable contributions to this field, shaping our understanding of the SP model and its applications in various scientific and engineering domains.
The article is organized as follows: Section 1 presents a detailed explanation of the SP model.In Section 2, the Lie group method and one-dimensional subalgebras for the SP model are discussed.Section 3 is devoted to exploring the invariant solutions of the SP model.Moving on to Section 4, the article examines the derivation of conservation laws using the homotopy operator.Section 5 focuses on the physical interpretation of the solutions.Lastly, in Section 6, the article concludes by highlighting potential future research directions.

Similarity Reductions and Invariant Solutions
Vector field J 6 = E 1 .
The characteristic equation associated with the vector field E 1 = ∂ ∂t is written as and provides a transformation θ(x, t) = k(r), r = x.With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, which gives, Hence, the solution of (4) in original variables becomes, Vector field J 7 = E 3 .
The characteristic equation associated with the vector field E 3 = ∂ ∂x is written as and provides a transformation θ(x, t) = k(r), r = x.With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, Hence, the solution of (4) in original variables becomes, Vector field The characteristic equation associated with the vector field and provides a transformation θ(x, t) = t 2 c 2 + k(r), r = x.With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, which gives, Hence, the solution of ( 4) in original variables becomes, The characteristic equation associated with the vector field and provides a transformation θ(x, t) = x c + k(r), r = t.With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, which gives, Hence, the solution of (4) in original variables becomes, The characteristic equation associated with the vector field and provides a transformation θ(x, t) = cxt + k(r), r = t.With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, Hence, the solution of (4) in original variables becomes, Vector field The characteristic equation associated with the vector field and provides a transformation θ(x, t) = k(r), r = t − x c .With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, We suggest solving the above equation numerically.
The characteristic equation associated with the vector field and provides a transformation θ(x, t) With the application of this transformation, we acquire the simplified version of Equation ( 4) presented as follows, We suggest solving the above equation numerically.

Local Conservation Laws via Homotopy Operator
When dealing with complex multipliers and/or equations in the context of inverting divergence operators, an effective approach involves the utilization of homotopy operators derived from differential geometry.These operators help simplify the task of determining fluxes, reducing it to a more manageable problem of integration in single-variable calculus.In this section, we introduce the initial series of such formulas, as outlined in reference [19].It is crucial to emphasize that when dealing with PDE systems and/or multipliers that are not excessively intricate but involve arbitrary constitutive functions, the direct method of flux computation is typically employed.
Let us consider a system of partial differential equations of order k, denoted as (x, θ (k) ) = 0, where x = (x 1 , x 2 , • • • , x n ) represents the independent variables and θ = (θ 1 , θ 2 , • • • , θ m ) represents the dependent variables.The system can be expressed as a collection of individual equations as follows Now, let us define a local conservation law for the system (48).This law is characterized by a divergence expression which holds true for all solutions of the system (48).Here, the operators D i represent the total derivatives with respect to the variable x i , while the quantities D i , i = 1, 2, • • • n, correspond to the fluxes associated with the conservation laws.The Euler operator is defined as follows Now, let us define the n-dimensional higher-order Euler operator concerning a function It is worth noting that when s 1 = 0 and s n = 0, we obtain the original Euler operator E (0,0) θ = E θ as defined in (50).
The n-dimensional homotopy operator [24] is introduced for an expression g[θ] = g(x, θ, ∂θ, • • • ), where θ = (θ 1 (x), θ 2 (x), • • • , θ m (x)) and x = (x 1 , x 2 , • • • , x n ).This operator is defined through its n components, each corresponding to a specific independent variable x i , i = 1, 2, • • • , n, and can be expressed as where ) is determined using the expression where j = 1, • • • , m.In essence, the n-dimensional homotopy operator enables us to deal with expressions g[θ] involving multiple variables θ and x and it provides a systematic way of handling each independent variable x i individually.The main theorem is presented below.
Proposition 2. Let g[θ] be a divergence expression given by [24] g and assume that g[0] = 0. Then the fluxes Ω i can be expressed as up to the corresponding fluxes of a trivial conservation law, provided that the integrals (53) converge.
By utilizing the direct (multiplier) method to solve Equation (4), we have derived the first-order multipliers as follows In the above equation, Ψ(x, t, θ, θ x ) represents the first-order multipliers obtained from the direct (multiplier) method applied to Equation (4).The constants C 1 to C 10 are coefficients determined during the solution process.Further details and the complete derivation process can be found in the reference [25,26].

Wave Nature of the Obtained Solutions
The graphical interpretation of a solution is of paramount importance as it offers a visual representation of complex mathematical relationships.It provides intuitive insights into the behavior, trends, and critical points of the solution.Graphs aid in understanding the sensitivity to varying parameters, validating results, and facilitating communication with a broader audience.They serve as valuable tools for exploration, model selection, and enhancing our understanding of intricate systems and phenomena.Figures 1 and 2 show the behavior of nonlinear longitudinal waves in an SP medium using the SP model.

Concluding Remarks
In this study, we successfully applied the Lie symmetry method to analyze the integrability properties of the SP model.By exploring the group invariant solutions via the reduction in symmetry subalgebras based on the optimal system, we gained valuable insights into the behavior of the system.Additionally, we listed the conservation laws using the homotopy operator, which proved advantageous over the traditional multiplier approach due to the absence of arbitrary functions in both the equation and characteristics.In the existing literature, Erofeev et al. [11] attempted to explore linear and non-linear plane longitudinal waves in the SP medium.However, their work did not present any soliton solutions or invariant ones.Our study focuses on invariant solutions, emphasizing solutions that remain unchanged under specific Lie symmetry transformations.This unique aspect adds novelty to our results.Furthermore, to illustrate the behavior of the solution at specific points in the model, we introduced nonlocal conservation laws reported for the first time.These findings contribute to the enhanced exploration of plane longitudinal waves in the SP medium via the application of the SP model.These findings shed light on the behavior of nonlinear longitudinal waves in the Slepyan-Palmov medium and the impact of various parameters on wave characteristics.These results motivate us to continue using the Lie symmetry method and the homotopy operator in tackling mathematical physics problems to further our understanding and contribute to this field of research.