Approximate symmetries of a mathematical model describing one-dimensional motion in a medium with a small nonlinear viscosity are studied. In a physical application, the approximate solution is calculated making use of the approximate generator of the first-order approximate symmetry.
MSC: 35J25, 32A37, 43A15, 35A58, 42B20.
Keywords:existence; uniqueness and stability of solutions; partial differential equations; approximate method
We consider the third-order partial differential equation
Equation (1) can describe the behavior of a one-dimensional viscoelastic medium in which nonlinearities appear not only in the elastic part of the stress, but also in the viscoelastic one.
Some mathematical questions as the global existence, uniqueness and stability of solutions can be found in [1,2]. Moreover, shear wave solutions are found in , where some explicit examples of blow-up for boundary value problems with smooth initial data are shown. A symmetry analysis and some exact solutions are shown in [4-6], while when , with a positive constant, a symmetry analysis can be performed in [7-9].
It is well known that a small dissipation is able to prevent the breaking of the wave profile allowing to study the so called ‘far field’, and a technique widely used is the perturbation analysis performed by expanding the dependent variables in power series of a small parameter (may be a physical parameter or often artificially introduced). Having in mind to perform an ‘approximate symmetry analysis’, we introduce in (1) a small parameter ε, namely
The combination of the Lie group theory and the perturbation analysis gives rise to the so-called approximate symmetry theories. The first paper on this subject is due to Baikov, Gazizov and Ibragimov . Successively another method for finding approximate symmetries was proposed by Fushchich and Shtelen . In the method proposed by Baikov, Gazizov and Ibragimov, the Lie operator is expanded in a perturbation series so that an approximate operator can be found. But the approximate operator does not reflect well an approximation in the perturbation sense; in fact, even if one uses a first-order approximate operator, the corresponding approximate solution could contain higher-order terms.
In the method proposed by Fushchich and Shtelen the dependent variables are expanded in a perturbation series; equations are separated at each order of approximation and the approximate symmetries of the original equations are defined to be the exact symmetries of the system coming out from equating to zero the coefficients of the smallness parameter. This method is consistent with the perturbation theory and yields correct terms for the approximate solutions but a ‘drawback’ is present: it is impossible to work in hierarchy, i.e., in the search for symmetries, there is a coupled system between the equations at several orders of approximation, therefore the algebra can increase enormously.
In this paper we work in the framework of the approximate method proposed in [12,13], in which the expansions of the dependent variable are introduced also in the Lie group transformations so that one obtains an approximate Lie operator which permits to solve in hierarchy the invariance conditions starting from the classification of unperturbed equation (3). We obtain the symmetry classification of the functions and through which equation (2) is approximately invariant and search for approximate solutions.
The plan of the paper is the following. The approximate symmetry method is introduced in the next section; the group classification via approximate symmetries is performed in Section 3; in Section 4, in a physical application, the approximate solution is calculated by means of the approximate generator of the first-order approximate group of transformations.
2 Approximate symmetry method
where a is the group parameter.
Upon formal substitution of (7) in (2), equating to zero the coefficients of zero and first degree powers of ε, we arrive at the following system of PDEs:
where we have set
In order to have a one-parameter Lie group of infinitesimal transformations of the system (9)-(10), which is consistent with the expansions of the dependent variables (7) and (8), we introduce these expansions in the infinitesimal transformations (4)-(6). Upon formal substitution, equating to zero the coefficients of zero and first degree powers of ε, we get the following one-parameter Lie group of infinitesimal transformations in the -space
where we have set
We give the following definition: We call approximate symmetries of equation (2) the (exact) symmetries of the system (9)-(10) through the one-parameter Lie group of infinitesimal transformations (11)-(14). Consequently, the one-parameter Lie group of infinitesimal transformations (11)-(14), the associated Lie algebra and the corresponding infinitesimal operator
are called the approximate Lie group, the approximate Lie algebra and the approximate Lie operator of equation (2), respectively.
Moreover, after putting
the approximate Lie operator (18) can be rewritten as
It is worthwhile noticing that, thanks to the functional dependencies of the coordinates of the approximate Lie operator (18) (or (20)), now we are able to work in hierarchy in finding the invariance conditions of the system (9)-(10): firstly, by classifying unperturbed nonlinear wave equation (9) through the operator (19) and after by determining and from the invariance condition that follows by applying the operator (20) to the linear equation (10). In fact the invariance condition of the system (9)-(10) reads
3 Group classification via approximate symmetries
The classification of equation (9) is well known (see for details Ibragimov  and bibliography therein). From (21), we arrive at the following result:
Taking (23)-(27) into account, from (22) we obtain the following additional conditions:
For f and λ arbitrary functions, we obtain
The classification of and with the corresponding extensions of arising from (32)-(38) is reported in Table 1.
Table 1. Classification ofandwith the corresponding extensions of
4 A physical application
Let us consider a homogeneous viscoelastic bar of uniform cross-section and assume the material to be a nonlinear Kelvin solid. The classical equation of motion (the constant density is normalized to 1 and the mass forces are neglected)
assuming a stress-strain relation of the following form:
where τ is the stress, x the position of a cross-section in the homogeneous rest configuration of the bar, the displacement at time t of the section from its rest position, the elastic tension ( is the strain), and is the viscosity component of the stress, reduces to (2).
So, we fall in Case II of Table 1 with the following identifications:
and from the corresponding invariant surface conditions, we obtain the following representation for the different terms in the expansion of w:
which give the form of an invariant solution approximate at the first order in ε.
After solving (49)-(50) and taking (47)-(48) into account, we have
Therefore, the invariant solution up to the first order in ε is
We have an unperturbed state represented by a stretching modified by the viscosity effect.
In this paper we perform the group analysis of the nonlinear wave equation with a small dissipation (2) in the framework of the approximate method proposed in [12,13]. In order to remove the ‘drawback’ of the method proposed by Fushchich and Shtelen , we introduce, according to the perturbation theory, the expansions of the dependent variables in the one-parameter Lie group of infinitesimal transformations of equation (2). Equating to zero the coefficients of zero and first degree powers of ε, we obtain an approximate Lie operator which permits to solve in hierarchy the invariance condition of the system (9)-(10) starting from the classification of unperturbed nonlinear wave equation (3). The proposed strategy is consistent with the perturbation point of view and can be generalized in a simple way to the higher orders of approximation in ε.
The authors declare that they have no competing interests.
The authors wrote this paper in collaboration and with the same responsibility. All authors read and approved the final version of the manuscript.
MR acknowledges the support of GNFM through the project 2012 Metodologie di tipo analitico e numerico per lo studio di problemi iperbolici ed iperbolico-parabolici di natura ondosa.
Dafermos, CM: The mixed initial-boundary value problem for the equations of non-linear viscoelasticity. J. Differ. Equ.. 6, 71 (1969). Publisher Full Text
MacCamy, RC: Existence, uniqueness and stability of . Indiana Univ. Math. J.. 20, 231 (1970). Publisher Full Text
Rajagopal, KR, Saccomandi, G: Shear waves in a class of nonlinear viscoelastic solids. Q. J. Mech. Appl. Math.. 56, 311 (2003). Publisher Full Text
Ruggieri, M, Valenti, A: Symmetries and reduction techniques for a dissipative model. In: Monaco R, Mulone G, Rionero S, Ruggeri T (eds.) Proceedings of WASCOM 2005,p. 481. World Scientific, Singapore (2006)
Fushchich, WI, Shtelen, WM: On approximate symmetry and approximate solutions of the non-linear wave equation with a small parameter. J. Phys. A, Math. Gen.. 22, L887–L890 (1989). Publisher Full Text
Capriz, G: Non-linear dynamics of a taut elastic string. In: Anile, AM, Cattaneo, G, Patano, P (eds.) Atti delle Giornate di Lavoro su Onde e Stabilità nei Mezzi Continui, Cosenza (Italy), 6-11 June 1983. Quaderni del Consiglio Nazionale delle Ricerche (CNR): Gruppo Nazionale per la Fisica Matematica (GNFM). Tipografia dell’Università, Catania (1986)