Two approximate symmetry frameworks for nonlinear partial differential equations with a small parameter: Comparisons, relations, approximate solutions

Abstract The frameworks of Baikov–Gazizov–Ibragimov (BGI) and Fushchich–Shtelen (FS) approximate symmetries are used to study symmetry properties of partial differential equations with a small parameter. In general, it is shown that unlike the case of ordinary differential equations (ODEs), unstable BGI point symmetries of unperturbed partial differential equations (PDEs) do not necessarily yield local approximate symmetries for the perturbed model. While some relations between the BGI and FS approaches can be established, the two methods yield different approximate symmetry classifications. Detailed classifications are presented for two nonlinear PDE families. The second family includes a one-dimensional wave equation describing the wave motion in a hyperelastic material with a single family of fibers. For this model, approximate symmetries can be used to compute approximate closed-form solutions. Wave breaking times are found numerically and using the approximate solutions, which yield comparable results.


Introduction
Perturbed differential equations (DE) are differential equations of the form involving a small parameter and may be considered a perturbation of some given differential equation F 0 [u] = 0 that corresponds to = 0. Such small perturbations can disturb the Lie symmetry properties of the unperturbed model in the sense that some or all exact point or local symmetries of the unperturbed equation disappear from the symmetry classification of the perturbed PDE. In particular, the perturbed model cannot have more symmetries holding for all values of than for a given value of , including = 0. In order to find extra symmetry-like structures for models with a small parameter, several approximate symmetry methods have been developed. Baikov, Gazizov and Ibragimov [4][5][6] introduced the approximate symmetry transformation method that is based on expanding approximate symmetry generators with respect to the small parameter (the BGI method). A different approach to approximate symmetries, developed by Fushchich and Shtelen (the FS method) [21], combines a perturbation technique with the symmetry group method by expanding the dependent variables in a Taylor series in the small parameter and approximately replacing the original equations by a system of equations that are the coefficients at different powers of the small parameter; the classical Lie symmetry method is consequently applied to obtain symmetries of the new system. The BGI and FS approaches are not equivalent. They have been compared and used to obtain approximate symmetries and approximate solutions for several PDE models [23,46]. Burde [11] developed a new approach for approximate symmetries by constructing equations that could be reduced by exact transformations to an unperturbed equation and at the same time would coincide approximately with the perturbed equation.
In this paper, we consider the BGI and FS approximate symmetry frameworks for perturbed partial differential equations (PDEs), with the main goal to investigate and compare the BGI and FS approximate symmetry structures of perturbed PDEs in more detail, and study their relations, and relate them to the symmetry structure of the unperturbed equations.
The majority of differential equations involve arbitrary parameters or arbitrary functions. These parameters have physical meanings and assume values in some ranges or belong to certain classes. Thus, to study the symmetry properties of system of differential equations involving arbitrary elements, one needs to investigate what happens to symmetries as these parameters assume special values. Namely, one classifies admitted symmetries depending on the forms of the arbitrary elements. At the first step towards the classification of symmetries of a differential equation-based model involving arbitrary elements, it is essential to find equivalence transformations for this model. Equivalence transformations map the given DE system to another differential equation from the same general class [26,30]. Lie groups of equivalence transformations can be sought systematically using techniques similar to those used to determine admitted Lie symmetry groups. The first group classification was given by Sophus Lie for a class of linear second-order partial differential equations in two independent variables [29]. Later, Ovsiannikov [35] found a complete classification for the nonlinear diffusion equation. More recently, classification problems for PDEs and ODEs were widely considered by many researchers (see, e.g., [1,19,20,24,27,31,45]).
The main set of partial differential equations used in this paper as testbeds and physical examples of interest for approximate symmetry computations and comparisons are one-dimensional nonlinear wave equations. Extremely broad literature exists dedicated to various analytical methods and studies of such models, including the studies of symmetry-related properties. In [2], group properties of the nonlinear wave equation were studied (here and below, subscripts denote partial derivatives). Bluman and Cheviakov [8] extended the group classification of (1.1) through a systematic construction of nonlocal symmetries. Point symmetry classifications for the generalized PDE classes u tt = (f (x, u)u x ) x , u tt = (f (u)u x + g(x, u)) x (1.2) were considered in [42,43]. The complete group classification for the PDE family was implemented in [7]. Further classifications of different classes of one-dimensional wave equation can be found, for example, in [22,28,44,48].
The study of wave propagation in nonlinear elastic materials has numerous applications in the study of complex materials [38], medical imaging [37], and other areas [10]. Of particular interest are hyperelastic solids, a class of materials that act as ideal elastic solids, and more general classes of models of this type. In particular, the stress within a hyperelastic solid is nonlinearly related to the deformation through a strain energy density function. In many settings, finite displacements in hyperelastic materials in one dimension reduce to nonlinear wave equations of the form where R is the stored energy function [14,18]. Among others, our goal here is to study and compare exact and approximate symmetry properties and to construct approximate solutions for an important special form of the PDE (1.4) describing fully nonlinear shear waves in a fiber-reinforced material: where c > 0 is the linear wave speed and is a real parameter that controls the nonlinear term. The above perturbed model describes waves in anisotropic fiber-reinforced elastic solids, such as biological membranes, in general, and including the cases when fiber effects are relatively small (see, e.g., [18,36]). While BGI and FS approximate symmetry classifications have been found for some models (see, e.g., [27,50] and references therein), the exact and approximate (BGI and FS) symmetry classification for PDEs of the form (1.5) has not been considered and compared in detail in the literature. As hyperbolic systems, wave equations possess characteristic curves (or surfaces) along which the solutions to the equation are simplified. If the characteristic curves intersect, the solution may become multi-valued. This is referred to as a shock or break in the wave and can have the physical meaning of a discontinuous solution [25]. In elementary models, breaks often form when the leading edge of the wave moves slower than the trailing edge, leading to a steeper and eventually vertical wave front. There are multiple approaches to determine the breaking time of a wave. If one can explicitly determine the characteristic curves, then their intersections and therefore the breaking time can be analytically calculated. When characteristic curves are not known explicitly, a breaking criterion can be developed based on the shape or velocity of the wave (see [47] and references therein). Tissier et al. [41] used local energy dissipation in addition to a maximal slope criterion to determine the breaking time. In this paper, we use approximate solution and finite-difference numerical simulations to find the approximate and numerical wave breaking times for the perturbed one-dimensional wave equation (1.5) with T(u x ) = u 2 x . The paper is organized as follows. In Section 2, we review the framework of Lie point and local symmetries of perturbed and unperturbed PDEs in detail necessary for the further analysis. We also provide a minimal introduction to BGI [4][5][6] and FS [21] approximate transformations and approximate symmetries for PDE models with a small parameter.
In Section 3, we study and compare the BGI and FS approximate symmetry frameworks for a perturbed PDE. In particular, we discuss types of BGI and FS approximate symmetries that can arise, including trivial BGI and FS approximate symmetries that have trivial action on O(1) solution components. We classify stable and unstable point symmetries of a family of wave-type equation with respect to the forms of the perturbation defined by an arbitrary function of two variables. For this example, it is observed that FS symmetries are stable in more classification cases than BGI symmetries. It is also shown (Section 3.4) that it is possible to have an unstable point symmetry of an unperturbed PDE that yields a point or a higher-order approximate symmetry neither in BGI nor in FS framework. This makes the PDE case radically different from the ODE case [40]. In Section 3.5, we derive a relation between BGI and FS approximate symmetries, holding in the case of specific forms of symmetry components.
In Section 4, we classify exact and approximate (BGI and FS) point symmetries of the perturbed wave equation (1.5) and use the results to construct approximate solutions of (1.5) in the case of power T , The approximate solutions are compared with numerical solutions for the most physically relevant case s = 2 and show good agreement. Even for smooth initial conditions, the wave dynamics leads to finitetime wave breaking. We use the numerically computed characteristic curves to approximate the wave breaking time. It is shown that the behavior of the approximate solution arising from approximate symmetry approach, near the time of breaking, can be used to provide an alternative estimate of the wave breaking time. There appears to be an inverse relationship between the wave breaking time and the small parameter . An overview of the results of the current paper and related open problems that would require extensions of approximate symmetry methods are discussed in Section 6.
The current work uses GeM module [15][16][17] for the Maple symbolic suite to generate and simplify determining equations for equivalence transformations and approximate symmetries.

Lie groups of exact and approximate point and local symmetries
We now briefly review the framework of Lie point and local symmetries in comparison with the BGI [4][5][6] and Fushchich-Shtelen [21] approximate symmetry frameworks for PDE models involving a small parameter.
For simplicity of presentation and analysis, in this work, we consider single PDEs rather than PDE systems.
Let x = (x 1 , x 2 , ..., x n ), n > 1, and u = u(x) denote respectively the independent variables and the dependent variable of a given problem. We also denote partial derivatives by subscripts: ∂u/∂x j ≡ u j , etc., and the set of all partial derivatives of u of order k by ∂ k u. A general k th -order scalar PDE on u has the form We assume that the PDE (2.1) as it stands, or after a point transformation, is of generalized Kovalevskaya type [33], that is can be written in a solved form with respect to the highest pure derivative of u by one of the independent variables. In (2.1) and below, F 0 [u] and similar notation denotes differential functions (i.e., functions depending on x, u, and derivatives of u up to some prescribed order s), defined in a domain of the jet space J s (x|u). (The latter is viewed as a multi-dimensional space with coordinates x, u, ∂u, . . . , ∂ s u. ) We also note that repeated indices, where appropriate, assume summation.
The solution set S of the PDE (2.1) in J s (x|u), s ≥ k, is a hypersurface defined by the relations F 0 [u] = 0 and its differential consequences ∂F 0 [u] = 0, . . ., solved for the corresponding differential consequences of the leading derivative, up to the highest order s. Any differential function f [u] can be evaluated on the solution set of (2.1) by substituting the expressions of the leading derivative and its differential consequences into f [u], and the result is denoted by

Local symmetries of a PDE
A one-parameter (a) Lie group of transformations in the space of the problem variables (x,u), with the corresponding infinitesimal generator defines a point symmetry of the PDE (2.1) when the solution set of (2.1) is invariant under the action of the mapping (2.2); that is, each solution u(x) is mapped to a solution u * (x * ) of the same PDE. The local condition of this invariance is given by the determining equation in terms of the prolonged generator X 0(k) . The determining equation splits into an overdetermined PDE system on the unknown symmetry components ξ i 0 , η 0 (see, e.g., [9,33]). For a point symmetry (2.2) of the PDE (2.1), the evolutionary (characteristic) form providing the same mapping between solutions has the form (x * ) i = x i , i = 1, 2, . . . , n, with the evolutionary component ζ 0 [u] = η 0 (x, u) − u i ξ i 0 (x, u), and the generator is given bŷ Local (point, contact and higher-order) transformations and the related point, contact and higher-order local symmetries of the PDE (2.1) generalize (2.5), (2.6) by allowing the evolutionary infinitesimal components ζ 0 = ζ 0 [u] to be general differential functions of u, depending on first-and/or higherorder derivatives of u (see, e.g., [9,33] and references therein). The invariance condition (2.4) is replaced byX The prolongation ofX 0 is defined by [33] X 0(k) = ζ 0 ∂ ∂u where the higher-order components are computed using A point or local symmetry of a PDE is trivial if its components vanish on the solution set of the PDE (2.1). Trivial symmetries provide identity transformations (2.2), (2.5). Symmetry generators are commonly simplified modulo trivial symmetries.
The following elementary example will serve as a basis of further examples involving PDEs with a small parameter, and their exact and approximate symmetries.
Example 2.1. Consider a nonlinear wave-type equation [34] u tt = u x u xx , u = u(x, t). (2.9) The exact symmetry generator for the PDE (2.9) is given by (2.10) The determining equations (2.4) yield the solution Consequently, the PDE (2.9) admits a six-dimensional Lie algebra of point symmetries spanned by (2.12)

Exact and approximate symmetries of a PDE with a small parameter
A general first-order perturbation of a PDE (2.1) is a partial differential equation involving a small parameter . We assume that the perturbation F 1 [u] is regular, in the sense that the Kovalevskaya forms of the unperturbed PDE (2.1) and its perturbation (2.13) have the same leading derivatives.

Exact local symmetries
Exact point and local symmetry generators of (2.13) have the forms (2.14) To find the exact symmetries of (2.13) holding for an arbitrary , one solves the determining equations (2.4) or (2.7). Since (2.13) is a PDE family that includes an arbitrary element , the dimension of Lie algebra of point or local symmetries holding for a general cannot exceed that for some fixed , including = 0. Therefore, the family (2.13) of perturbed PDEs will admit the same or smaller number of local symmetries than its unperturbed version (2.1).

Example 2.2.
We compute exact point symmetries of a perturbed version of the PDE (2.9), u tt + uu t = u x u xx (2.15) holding for an arbitrary . The leading derivative u tt can be chosen for both (2.9) and (2.15). We obtain a Lie algebra of point symmetries spanned by a three-dimensional subalgebra of the six-dimensional Lie algebra of point symmetries (2.12) of the unperturbed wave equation (2.9).

Baikov-Gazizov-Ibragimov approximate point and local symmetries
Approximate symmetries can be useful for finding additional symmetry-like structures for the perturbed equation (2.13) that are not its exact symmetries but rather preserve (2.13) approximately, up to o( ) [27]. A one-parameter Lie group of Baikov-Gazizov-Ibragimov (BGI) approximate point transformations with the group parameter a, acting on the (x,u)-space, is given by .., n, (u * ) = g(x, u; a, ) = g 0 (x, u; a) + g 1 (x, u; a) + o( ), (2.17) where, as in (2. 2), f i j , g j are sufficiently smooth functions. The generator of the approximate group (2.17), up to o( ), is  Below in this paper, Definition 2.2 will be used not only for BGI approximate point symmetries, but more generally, for BGI and FS approximate point and local symmetries.
It is straightforward to show that O( ) term of the determining equation (2.19) leads to the PDEs where H is obtained from the coefficients of in the expression  . This leads to the symmetry generated by X 0 being unstable (or respectively, stable). If all symmetries of the equations (2.1) are stable, the perturbed equations (2.13) are said to inherit the symmetry structure of the unperturbed equations [27].

Remark 2.2.
Similarly to exact local transformations with generators of the form (2.3), one can define more general local approximate BGI transformations with generators in evolutionary form given bŷ Approximate local (including point, contact and higher-order) BGI symmetries of the perturbed PDE (2.13) can be found using the same procedure as described above for point approximate BGI symmetries.
In particular, the analog of the first-order condition (2.20) takes the form

Fushchich-Shtelen approximate point and local symmetries
Unlike the BGI approach where the symmetry generator is expanded in a power series in terms of the small parameter, the Fushchich-Shtelen method [21] applies the perturbation technique to the solution u(x) and the given PDE. In particular, the solution is written as Substituting (6.3) into the PDE (2.13) with a small parameter, expanding the result and setting to zero the O(1) and O( ) terms independently, one obtains a coupled system of equations on v(x) and w(x) without the small parameter, given by

Definition 2.3. A Lie group of point transformations with the group parameter a
with the generator

defines a FS approximate point symmetry of the PDE (2.13) if it is an exact Lie point symmetry group of the Fushchich-Shtelen system (2.25).
In a similar manner, a generalized local (point or higher-order) transformation group in the evolutionary form could contain additional conditions on the v-component ψ 1 . (2.30) Such FS symmetries do not originate from stable local symmetries of the unperturbed PDE u tt = u xx and cannot arise in the BGI framework.

Trivial BGI approximate symmetries
Consider a local BGI approximate transformation with the evolutionary generator (2.22): which means that suchX 1 must be a local symmetry generator of the unperturbed equation (2.1). The opposite is also true: ifX 0 is a local symmetry generator of the unperturbed equation (2.1), then is a BGI approximate symmetry generator of the perturbed PDE (2.13). In the light of the above, we call a BGI approximate symmetry that has a generator with vanishing O(1) part a trivial BGI approximate symmetry. This triviality relates not to the trivial action of such symmetries but rather to the fact that every local symmetryX 0 of the unperturbed equation ( with the first Taylor term of the transformation having the order ∼ a = o(a, ).

Trivial FS approximate symmetries
is an evolutionary component of the local symmetry of the unperturbed equation (2.1) generated bŷ Indeed, the action of (3.5) on the first PDE (2.25a) of Fushchich-Shtelen system is trivial, and the action on the linear equation (2.25b) is equivalent to the local symmetry determining equation (2.7) of the unperturbed PDE (2.1). The converse also holds: every local symmetry generator (3.6) of the unperturbed equation yields a trivial FS approximate symmetry generator (3.5).

Types of approximate symmetries arising in classifications
In the computation of BGI approximate symmetries of a PDE (2.13) with a small parameter, the following three types of symmetries can arise.
For FS approximate symmetries, the following three types can arise. Genuine BGI and FS approximate symmetries are the main focus of the approximate symmetry study.

A computational example: Exact and approximate point symmetry classification for a second-order nonlinear PDE with a small parameter
In this section, we compare exact point symmetries (2.12) of the (1+1)-dimensional wave-type equation (2.9) u tt = u x u xx with BGI and FS approximate point symmetry classifications for the family of perturbed equations where F 1 (u, u t ) is an arbitrary function. The computation is limited to point symmetries for simplicity of presentation. The cases for the arbitrary function F 1 (u, u t ) are given modulo equivalence transformations where C 1 and C 2 are arbitrary constants. The PDE (2.9) admits six exact point symmetries given by (2.12). The BGI approximate point symmetries are computed and classified following the procedure described in Section 2.2.2, and the FS approximate local symmetries are obtained following Section 2.2.3. In particular, the Fushchich-Shtelen system (2.25) for (3.7) is given by where u(x, t) = v(x, t) + w(x, t). For the Fushchich-Shtelen system (3.9), point equivalence transformations are given by (3.8) for x, t and F 1 , with dependent variable transformations v =ṽ, The resulting classification is presented in Table 1. In the table, Q i denote arbitrary functions of their arguments and a j arbitrary constants. The table is organized as follows. The first column lists evolutionary forms of the six point symmetry generators X 0 k (2.12), k = 1, . . . , 6, of the unperturbed PDE (2.9). The second column lists the forms of the arbitrary function F 1 (u, u t ) for which the corresponding X 0 k is stable in the BGI sense, and the corresponding BGI approximate point symmetry of the perturbed wave equation (3.7). The third column contains the same information for the FS approximate point symmetries of the perturbed equation (3.7). Trivial approximate symmetries are not listed (see Section 3.1).
In Table 1, the generators of BGI and FS symmetries were sought in the evolutionary forms, depending on at most first derivatives of the fields. As a result, approximate BGI and FS point symmetries were obtained, as well as some local first-order approximate FS symmetries. We also note that for the point symmetry defined byX 0 5 to be stable in the FS sense, the arbitrary function F 1 = F 1 (v, v t ) must satisfy three coupled nonlinear PDEs. Because of the complexity of the latter, a general form of F 1 could not be obtained explicitly; however, it was verified that the form of F 1 as listed satisfies those equations. The same takes place for the point symmetry defined byX 0 6 . Table 1 illustrates differences between BGI and FS frameworks as there are stable symmetries in one framework and unstable in the other framework. Some specific examples are listed below. The following observations can be made.
• When F 1 (u, u t ) = u t or F 1 = const, all point symmetries (2.12) of the unperturbed PDE (2.9) are stable as BGI and FS symmetries. • For all forms F 1 (u, u t ), the PDE (3.7) and the FS system (3.9) are invariant under t-and x-translations.
• As an example, as seen in the first row corresponding toX 0 1 , for F 1 = uu t , this symmetry is stable as a point BGI approximate symmetry and as a local (but not point) first-order FS approximate symmetry.
• It can be separately proven that for every form of F 1 where a symmetryX 0 i is stable as a BGI point approximate symmetry, it is also stable as a local (point or first-order) FS approximate symmetry.
• The converse is not true; for example, the scaling symmetryX 0 6 is stable as a local FS approximate symmetry for F 1 = v a 3 Q 6 (v t ) + a 2 v t + a 1 , but is only stable as a BGI approximate symmetry when a 3 = −1.

Instability of local symmetries of unperturbed PDEs in terms of higher-order approximate symmetries: An example
For an ordinary differential equation (ODE), all local symmetries are stable in the BGI sense: each local symmetry of a given ODE corresponds to an approximate BGI local, often higher-order, symmetry of its perturbed version [40]. For a PDE, in general, this is not the case. For the BGI framework, differential functions ζ 1 (p) ..ip in the determining equation (2.23) contain derivatives of u of orders higher than those in the differential function ζ 1 . It follows that the left-hand side of equation (2.23) splits into a system of linear PDEs in ζ 1 . On the other side, the function H may contain derivatives of u with respect to other variables different than those in the left-hand side of equation (2.23). This can lead to some constraints on the unperturbed symmetry component ζ 0 ; in that case, an exact local symmetry of the unperturbed PDE (2.1) may not correspond to a local approximate BGI symmetry of the perturbed PDE (2.13). A similar argument holds for FS approximate symmetries. The main reason, as it can be seen in the example below, is the existence of multiple kinds of derivatives in PDEs and thus more restrictive conditions that arise for ζ 0 when the determining equations are being split with respect to higher-order derivatives.
As an illustration, consider the PDE (3.7) with F 1 (u, u t ) = uu t : and the related Fushchich-Shtelen system (3.9) From Table 1, one can see thatX 0 2 = t ∂/∂u is unstable as a point symmetry in both BGI and FS frameworks; that is, the point symmetryX 0 2 admitted by the PDE (3.11) with = 0 corresponds to no approximate point or first-order symmetry arising from BGI or FS approaches. First, we examine whether or not it is possible to construct a local, possibly higher-order, BGI approximate symmetry of (3.11) that would correspond toX 0 2 . The generator of such a symmetry would have the form where C 2 = const = 0. The determining equation (2.23) for BGI local symmetries reads where one readily finds One can show by a direct computation that whatever the dependence of ζ 1 on partial derivatives of u is chosen to be, higher-order derivatives of u that arise in (3.14) lead to constraints on C 2 that result in C 2 = 0, which means that no nontrivial BGI point symmetry (3.13) corresponding toX 0 2 exists. Second, we seek a local, possibly higher-order, approximate FS symmetry of the PDE (3.11) corresponding toX 0 2 . Such a symmetry would arise as an exact local symmetry of Fushchich-Shtelen system (3.12). The corresponding evolutionary generator (2.29) has the form As noted in Section 2.2.3, the determining equation for the first equation of the system (3.12) is satisfied when ψ 1 = C 2 t as inX 0 2 . Now the determining equation for the second PDE of (3.12) leads to It can be shown that for any dependence ψ 2 [v, w], constraints on C 2 exist, leading to C 2 = 0. Consequently, there is no higher-order FS symmetry corresponding to the unstable point symmetrŷ X 0 2 admitted by the wave equation (3.11) with = 0.

A relation between BGI and FS approximate symmetries
The computational example of Section 3.3 above illustrated the fact that BGI and FS frameworks can yield rather different approximate point symmetry classifications for the same PDE with a small parameter. However, in certain situations, the two approaches can lead to related results. We now show that for a specific class of (1+1)-dimensional PDEs, a stable BGI approximate point symmetry always correspond to a stable FS approximate local symmetry. Consider the following class of PDEs on u(x,t), written in the Kovalevskaya form with respect to an independent variable t: and its perturbed version with a small parameter :  (3.19) and split the orders of to get the Fushchich-Shtelen system (3.22) The evolutionary generator of a local FS approximate symmetry has the form (2.29). The determining equations (2.4) for exact local symmetries of (3.22) arê holding on solutions of (3.22).

Theorem 3.1. If (3.20) is a BGI approximate local symmetry generator of a PDE (3.19) having the specific formX
= (ζ 0 (x, t) + ζ 1 (x, t, u, u x , u t )) ∂ ∂u (3.24) and additionally, F 0 [v] in (3.19) satisfies the following system of equations is a FS approximate local symmetry of the perturbed PDE (3.19) to the point symmetry generatorX 0 = ζ 0 ∂/∂v of the unperturbed PDE (3.18 The second FS determining equation (3.23b) with ψ 1 = ζ 0 , ψ 2 = ζ 1 can be rewritten as As ζ 0 and F 0 [v] satisfy (3.25), G reduces to which is equivalent to the left-hand side of (3.27). Now, the right-hand side of (2.20), the function H, is the coefficient of in

(3.29)
Since ζ 0 = ζ 0 (x, t), none of the terms in (3.29) contains ∂ n u/∂t n . Hence, the coefficient of in (3.29) is The latter is equivalent to G (3.28). It follows that the determining equation ( The above theorem states that when a point symmetry of an unperturbed PDE yields a BGI approximate point symmetry but not a FS approximate point symmetry of the perturbed PDE, under the conditions of the theorem, there exists a corresponding higher-order FS approximate symmetry of the perturbed PDE instead. This is illustrated by the following example. Example 3.1 Consider again the PDE (3.11) u tt + uu t = u x u xx . Using Table 1, we observe that X 0 1 = ∂/∂u is unstable as a FS approximate point symmetry but it is a stable point symmetry in the sense of BGI; the corresponding BGI generator is given bŷ with H (3.30) given by (3.15) with C 1 = 1. Consider now the Fushchich-Shtelen system (3.12) for the PDE (3.11). Using determining equation (3.5) for exact local symmetries of (3.12), one can find that ) ∂ ∂w is a higher-order FS approximate symmetry generator of the PDE (3.11).  admits an approximate point symmetry with the evolutionary form (3.32) corresponding to the stable point symmetryX 0 = ζ 0 ∂/∂u = (2 − xu x ) ∂/∂u. Here, ζ 0 = 2 − xu x does not satisfy the conditions of Theorem 3.1 since it involves u x . It turns out thatX 0 is unstable as a FS approximate point symmetry of (3.31). Indeed, it is easy to check that is not a local symmetry of the Fushchich-Shtelen system of the PDE (3.31) given by However, a direct computation shows that a slightly modified version of (3.33) given bŷ ∂ ∂w yields the corresponding first-order FS approximate symmetry.

Exact and approximate point symmetry classification of a one-dimensional perturbed wave model in a fiber-reinforced solid
One-dimensional nonlinear wave equations on the unknown u(x,t) and various forms of K(u x ) arise in multiple physical contexts, in particular, in nonlinear mechanics [32]. The point symmetry classification of the PDE family (4.1) has been performed by Oron and Rosenau [34]. If K(u x ) = c 2 = const, the PDE (4.1) becomes linear: and consequently, its Lie symmetry algebra is infinite-dimensional and consists of the vector fields of the form parameterized by an arbitrary constant C 1 and four arbitrary functions α 1 (x + ct), β 1 (x + ct), α 2 (x − ct), and β 2 (x − ct).
In the current section, we consider a special form of the arbitrary function K(u x ) = c 2 + Q(u x ) in (4.1), which yields a PDE family with a small parameter . It is assumed that Q(u x ) = const. Such models arise, for example, in the analysis of wave propagation in fiber-reinforced elastic solids [3,13] with small fiber strengths. The PDEs (ũx), (4.5) involving arbitrary constants C i , i = 1, . . . , 8. It follows that by taking C 1 = 1/c, C 3 = C 5 = 1, and other constants zero, upon dropping tildes, one obtains the PDE (4.4) with c 2 = 1: which will be considered below.
The results below are presented modulo the equivalence transformations (4.5), usually without obvious trivial approximate symmetries (see Section 3.1); some trivial approximate symmetries will be pointed out.

Exact point symmetries of the wave equation family (4.6)
The exact symmetry generator for the PDE (4.6) has the form The following cases arise, holding for an arbitrary and non-constant Q.
1. In the general case of arbitrary Q(u x ), one has the five-dimensional Lie group of point symmetries generated by describing respectively translations in t, x, u, the Galilei transformation in the direction of the displacement u, and a homogeneous space-time scaling. 2. In the case when Q(u x ) = u x , the Lie algebra (4.8) is extended by a point symmetry generator (4.9)

BGI approximate point symmetry classification of the wave equation family (4.6)
The BGI approximate point symmetry generator for the PDE (4.6) has the form where, according to Theorem 2.1, the freedom in X 0 does not exceed that in X 0 ∞ (4.3). From the determining equations (2.19) for BGI point symmetries, the following cases arise.
1. Q(u x ) arbitrary: the O(1) and O( ) components of the generator (4.10) are given by where C i = const, λ 1 and λ 3 are arbitrary functions of x + t, and λ 2 and λ 4 are arbitrary functions of x − t. Consequently, the nonlinear wave equation (4.6) for an arbitrary Q(u x ) admits the approximate point symmetries which are, respectively, the re-numbered exact point symmetries (4.8), and a trivial approximate symmetry X ∞ corresponding to the infinite symmetry set (4.3) of the linear wave equation (4.2). The difference between the freedom in (4.3) and the space spanned by X 1 , . . . , X 5 in (4.12) corresponds to unstable point symmetries of the linear wave equation.
2. Q(u x ) = A ln (u x + B) + C, where A, B and C are arbitrary constants: here, the nonlinear wave equation (4.6) admits an additional genuine approximate symmetry given by (4.13) 3. Q(u x ) = u x : the exact symmetry generator (4.3) of the linear wave equation (4.2) reduces to 14) and the O( ) approximate symmetry components have the form H(x, t), where H(x,t) is an arbitrary solution of the PDEs: In this third case, the point symmetries of the unperturbed linear wave equation (4.2) with arbitrary β 1 (x + t) and β 2 (x − t) remain stable and yield genuine approximate symmetries with O( ) components given by the terms in (4.15) that contain β 1 and β 2 . In particular, a genuine approximate symmetry similar to (4.13) is admitted: (4.16)

FS approximate point symmetry classification of the wave equation family (4.6)
For the perturbed PDE (4.6) with the solution form the Fuschich-Shtelen system (2.25) reads The basic group of point equivalence transformations of the PDE system (4.18) is given by in terms of arbitrary constants C 1 , . . . , C 11 . We now find exact point symmetries of the system (4.18) that correspond to FS approximate point symmetries of the PDE (4.6). The infinitesimal generator of such symmetries has the form (4.20) The solution of the determining equations (2.4) leads to the following classification.

Q(v x )
= v s x , s = 0: (4.23) The symmetries given by Z 7 and Z 7 are genuine FS approximate point symmetries of the perturbed PDE (4.6). In particular, Z 7 corresponds to a scaling which is different for O (1) and O( ) components of the approximate solution (4.17).

Summary
For an arbitrary Q, the perturbed one-dimensional wave equation (4.6) admits five exact symmetries given by (4.8) and it has these five symmetries and a trivial approximate symmetry as BGI approximate symmetries (4.12). For Q = u x , the equation (4.6) admits (4.8) and an additional exact symmetry given by (4.9). For BGI classification with Q(u x ) = u x , the PDE (4.6) has an infinite set of BGI approximate symmetries with approximate symmetry components given by (4.15). Note that the exact symmetry generator Y 6 in (4.9) can be obtained from the BGI approximate components (4.15) by taking It follows that the BGI approximate symmetry classification of the wave equation (4.6) includes the exact symmetry classification of (4.6) but corresponds to a subset of exact point symmetries (4.3) of the unperturbed (linear) wave equation (4.2). An additional case appears in the BGI approximate symmetry classification when Q = A ln (u x + B) + C, with a corresponding additional approximate symmetry given by (4.13). This case does not arise in the FS symmetry classification.
For an arbitrary Q, the PDE (4.6) admits exact point symmetry generators (4.8) and trivial approximate symmetries given by (4.21). In comparison with the exact and BGI symmetry classifications of (4.6), two different cases appear in FS approximate symmetry classification: Q = v s x , s = 0 and Q = e vx . For Q = v s x , the PDE (4.6) admits an additional FS approximate symmetry given by (4.22). (By contrast, in the exact and BGI symmetry classifications of (4.6), this case appears only when s = 1.) For Q = e vx , a stable exact symmetry x ∂/∂u of the linear wave equation (4.2) yields a genuine FS approximate symmetry of (4.6) given by (4.23).

Approximate and numerical solutions modeling breaking waves in fiber-reinforced materials
The displacements in shear waves propagating in an incompressible hyperelastic material with a single family of fibers directed along the wave propagation are governed by a nonlinear one dimensional wave equation where the constants α, β > 0 are the material parameters [13]. In this section, we consider wave equations (4.4) with Q(u x ) = Bu s x , B > 0, s = 0, which include the model (5.1). By a re-scaling of x, t and u, these PDEs can be brought into a simpler form

Approximate solutions of the PDE (5.2) in the FS framework
Here we use Fuschich-Shtelen approximate symmetries to construct an approximate solution for the PDE (5.2) in the usual FS form In the first-order of precision in , the equation (5.2) is equivalent to the Fuschich-Shtelen system (4.18) with which admits the symmetry generator (4.22). The corresponding characteristic equations are given by Substituting (5.6) and w = g s+1 φ into the second PDE of (5.4), one gets to the PDE in φ When s = −1, the PDE (5.7) has a general solution where h = h(x, t) satisfies h tt = h xx . Similarly, when s = −1, the solution form changes to In light of the above results, the higher-order solution part w has the form Finally, when s = −1, the perturbed equation (5.2) has the approximate solution (5.3) given by When s = −1, the approximate solution takes the form Example 5.1. As a specific example that will be used below, we consider the PDE (5.1) describing shear waves in a fiber-reinforced solid, re-scaled to the form (5.2) with s = 2: and assume in (5.1) that β/α ∼ 1, which corresponds to weak fiber effects. We also choose v(x, t) = exp (−(x − t) 2 ). Then, the solution (5.9a) of the PDE (5.2) with h = 0 reduces to Note that for any fixed t, the approximate solution (5.11) approaches zero as x → ∞. Also, for any x ∈ (−∞, ∞), the solution (5.11) is bounded as t → ∞. The solution (5.11) is not a purely right-traveling wave solution but describes an evolving wave form (see Section 5.2 below). In particular, the PDE (5.10) is known to have breaking wave-type solutions [13].

Numerical simulations of (5.10) and breaking waves
We now compute numerical solutions of the wave equation (5.10) in Example 5.1 in order to model its breaking wave behavior (see Section 5.3 below for details) and provide a reference for comparison with the approximate solutions developed in Section 5.1. Gaussian-type initial conditions corresponding to a right-traveling wave and periodic boundary conditions (5.12) are posed in the space-time domain x ∈ [−L, L], t ≥ 0, and the equation (5.10) is solved using an explicit finite-difference cross-stencil scheme with constant spatial and temporal stepsh,τ . Following [12], we use a conservative finite-difference scheme developed for the PDE (5.10): with U = U n m approximating the value of u(x, t) at the mesh node (x m , t n ). Here, U tť and U xx represent the second-order central differences, U x the first-order forward difference and Ux the first-order backward difference: 14) The numerical solutions provide a good agreement with approximate solutions (5.11), for a broad range of values of , from the initial dimensionless time to the time close to the wave breaking. The breaking time increases approximately as −1 as decreases (see Section 5.3 below). Here we present sample computation results for a relatively large value of the small parameter, = 0.5. The computation is performed from t = 0 to t = 4 close to the wave breaking time. A comparison of the numerical solution and the approximate solution (5.11) of PDE (5.10) with initial and boundary conditions (5.12) at several time snapshots is presented in Figure 1(a). The relative difference at the time step t n between the approximate and numerical solutions is calculated using 2-norms according to the formula and is shown in Figure 1(b).

Estimates of wave breaking times using approximate and numerical solutions
As discussed in [13], the nonlinearity (1 + (u x ) 2 ) leads to greater characteristic speed values at points where |u x | is larger. This can lead to the intersection of characteristic curves, which corresponds to the wave breaking and the formation of a shock. This behavior can be studied using the method of characteristics. While for linear hyperbolic PDEs, such as the constant-coefficient wave equation u tt = c 2 u xx in the simplest case, characteristic curves can be found in terms of explicit formulas such as x = x 0 ± ct and lead to explicit exact solutions, the situation is significantly more complex for nonlinear hyperbolic PDEs. Using the method described in [25,49], one can show that the equation (5.10) can be reduced to the first-order characteristic form on the characteristic curves dx dt In physical terms, the part of the wave with a time derivative given by (5.16) moves at a finite velocity given by (5.17). The integration of (5.17) yields a constant of integration x 0 that corresponds to the point on the characteristic curve where t = t 0 is some initial time. Thus, x = x(x 0 , t), from which u t (x(x 0 , t), t) in ( For example, when the numerical computation has reached the time layer t = t n , the linearized characteristics are launched from each spatial grid point x m , m = 1, . . . , M − 1. In particular, the right-traveling characteristics are approximated by the lines where u x = U x is the first-order forward finite difference (5.14). To approximate the breaking time t = τ numerically, we solve (5.18) for the time τ when two different characteristics intersect. Given two starting points, x m 1 and x m 2 , we get the system Solving for τ yields where are the slopes of the characteristic lines. We choose x m 1 and x m 2 to be adjacent grid points, x m 1 +h = x m 2 . The numerator of (5.20) is constant, so the approximate breaking time corresponds to the largest denominator of (5.20). We determine x m 1 corresponding to the largest difference between the slopes m 1 and m 2 , then solve for τ . The meaning of τ is thus the estimated time from t n to the wave breaking time T b ; one consequently has an estimate As the wave evolves, the slopes of the linearized characteristics will change and therefore so will τ . To account for this, at each time step, we repeat the calculation for τ . We use the first-order forward finite-difference approximation to compute u x in (5.21) at each time step. Figure 2   spatial step sizesh = 0.01 andh = 0.005 and temporal step sizeτ = 0.00125, we found the break times for each in the limit ash → 0 as shown in Figure 2(a). This approach uses the fact that as the wave approaches the breaking time, the second spatial derivative u xx tends to negative infinity (see Figure 2(b)). We determined the actual, non-linearized characteristic curves by numerically integrating (5.17). The curves are shown in Figure 3.
It is interesting that one can also approximately determine the wave breaking time from the approximate solution (5.11) by finding the time when the second spatial derivative u xx of the approximate solution develops an additional root, as shown in Figure 4. This corresponds to an additional inflection point in the wave itself, which can be observed in Figure 1(a) closer to the wave breaking.
The numerically determined break time values (τ num ) and the approximate-determined break time values (τ approx ) are given in Table 2. Both sets of data are also plotted in Figure 5. The results suggest that the wave breaking time behaves like τ ∼ −1 .

Discussion
In this paper, we considered exact and approximate symmetries of partial differential equations with small parameter, in the Baikov-Gazizov-Ibragimov (BGI) [4][5][6] and the Fushchich-Shtelen (FS) [21] frameworks. The goals of this paper were to study relations, properties and applications of approximate symmetries of perturbed PDEs in comparison with exact symmetries of the unperturbed and perturbed PDEs and to use the approximate symmetries to construct approximate solutions of nonlinear PDEs. For ordinary differential equations with a small parameter, all BGI approximate symmetries are stable: each symmetry of an unperturbed equation reappears as a point or higher-order approximate BGI symmetry [40]. For PDEs, the situation is different: we found that the point symmetry of the unperturbed PDE, in general, does not correspond to a point or a higher-order approximate symmetry of the perturbed model. The main reason of this symmetry instability is that the determining equations (2.23) for BGI local symmetries of the perturbed PDE (2.13), whatever the dependence of the approximate symmetry components ζ 1 (2.22), always contains derivatives of u higher than those in ζ 1 . This leads to a split system of PDEs in ζ 1 with some restrictions on the unperturbed symmetry components ζ 0 . A similar argument holds in FS framework. As an illustration, in Section 3.4, we showed that there was no higher-order (BGI and FS) approximate symmetry for the perturbed PDE (3.11) corresponding to the point symmetry t∂/∂u of the unperturbed wave equation (2.9). A detailed stability classification of point symmetries of the nonlinear wave equation (2.9) in terms of approximate BGI and FS frameworks was performed in Section 3.3. The results illustrate that the two approximate symmetry are significantly different. Yet as shown in Section 3.5, there exists a connection between the BGI and FS approximate symmetries; in particular, each stable BGI point symmetry of the form (ζ 0 (x, t) + ζ 1 (x, t, u, u x , u t ))∂/∂u yields a higher-order approximate FS symmetry in the form ζ 0 (x, t)∂/∂v + ζ 1 (x, t, v, v x , v t )∂/∂w (Theorem 3.1).
In the FS framework, we found that there exist FS approximate point symmetries of perturbed PDEs that do not correspond to the stable point symmetries of the unperturbed models and also cannot appear in the BGI framework (Remark 2.3). In such FS symmetries, the O(1) and O( ) solution and symmetry components 'mix' together; in particular, the O(1) symmetry component depends on O( ) solution component. An example is given by (2.30).
The classifications of exact and approximate (BGI and FS) symmetries of a perturbed PDE not only illustrates differences between these approximate symmetry approaches but also leads to different types of approximate symmetries that can be used to construct approximate solutions for the given PDE. In Section 4, exact and approximate point symmetries for the perturbed one-dimensional wave equation (4.6) arising in the context of nonlinear waves in elastic materials were classified. It was shown that the exact symmetry classification of the perturbed PDE (4.6) is included in the BGI approximate symmetry classification of (4.6). Moreover, the classifications of BGI and Fushchich-Shtelen approximate symmetries of (4.6) yielded different classification cases in terms of the arbitrary constitutive function. While the BGI classification included a logarithmic case that did not appear in the FS classification, two new cases appeared in the FS classification: Q(v x ) = v s x and Q(v x ) = e vx that did not arise in the BGI framework except for s = 1. All these cases yielded genuine BGI and FS approximate point symmetries for the nonlinear wave equation (4.6). In Section 5, for the nonlinear wave PDE family (5.2) with a small power nonlinearity, an FS approximate symmetry generator was used to derive explicit families of approximately invariant solutions (5.9). For the case of quadratic nonlinearity (5.10), these approximate solutions were compared with finitedifference numerical solutions based on the method of characteristics. The latter indicate the appearance of breaking waves from smooth initial conditions in finite time. Though the approximate solutions remain smooth, their behavior can be used to predict the breaking times; these breaking times showed a good agreement with those estimated from the numerical computations. Both approaches indicate that the breaking time is inversely proportional to the magnitude of the small parameter that controls the nonlinearity.
The current work employed GeM module [15][16][17] for Maple symbolic manipulation software to generate and simplify determining equations for approximate symmetry and equivalence transformation components, the determining equations which were subsequently solved with standard Maple routines.
The approximate symmetry area still contains many open questions, including the general question of whether or not there exists some 'unified' approach that would include, and perhaps extend, both BGI and FS approximate symmetry frameworks. Another related open question is that of systematic computation of approximate conservation laws and approximately conserved quantities of dynamic ODEs and PDEs, and relations between such conservation laws and approximate symmetries. An obvious area of application of approximate conservation laws is the development and validation of specialized numerical methods for nonlinear DE models with small parameters.
A further extension of interest is the computation of approximate symmetries of DEs that include higher-order perturbation terms: Another open area in the domain of approximate symmetries involves DE systems with more than one independent small parameters. Such models arise in multiple contexts. For instance, the dimensionless Serre-Su-Gardner-Green-Naghdi equations [39] u t + εuu x + h x = δ 2 3h (h) 3 u xt + εuu xx − ε(u x ) 2 x , h t + (hu) x = 0 involve two independent parameters ε and δ; a dimensionless viscoelastic wave model [13] u tt = α + 3βu 2 x u xx + ηu x 8u 2 x + 2 u xx u tx + 2u 2 x + 1 u x u txx + ζ u 3 x 24u 2 x + 4 u xx u tx + 4u 2 x + 1 u x u txx includes four real parameters, three of which may be small: β, η, ζ α.