Journal of Applied Mathematics
Volume 2015 (2015), Article ID 601657, 8 pages
http://dx.doi.org/10.1155/2015/601657

Research Article

Derivation of Asymptotic Dynamical Systems with Partial Lie Symmetry Groups

Masatomo Iwasa

General Education Department, Aichi Institute of Technology, 1247 Yachigusa Yakusacho, Toyota 470-0392, Japan

Received 16 July 2015; Revised 12 September 2015; Accepted 14 September 2015

Academic Editor: Wan-Tong Li

Copyright © 2015 Masatomo Iwasa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Lie group analysis has been applied to singular perturbation problems in both ordinary differential and difference equations and has allowed us to find the reduced dynamics describing the asymptotic behavior of the dynamical system. The present study provides an extended method that is also applicable to partial differential equations. The main characteristic of the extended method is the restriction of the manifold by some constraint equations on which we search for a Lie symmetry group. This extension makes it possible to find a partial Lie symmetry group, which leads to a reduced dynamics describing the asymptotic behavior.

1. Introduction

Appropriate reduction of a dynamical system is one of the important approximation approaches when the system is too complicated to be solved analytically. The derivation of a reduced dynamical system which describes the asymptotic behavior of the solution has been achieved by using various singular perturbation methods developed for ordinary differential equations, partial differential equations, and difference equations [1–11]. On the other hand, the Lie group analysis, a systematic method to construct solutions with Lie symmetry groups that leave the manifold formed by the dynamical system invariant [12], has been extended in order to obtain approximation solutions or reduced dynamics by considering approximate symmetries [13–19], asymptotic symmetries [20–25], and renormalization group symmetries [26–30].

While the above singular perturbation methods and the Lie symmetry group analysis had been developed independently, relations between them have been studied recently. It has been shown that the asymptotic dynamics derived with the singular perturbation methods can also be derived with the Lie group analysis when it is applied to singular perturbation problems of ordinary differential equations or difference equations [31–33]. In other words, the Lie group analysis plays also a role of a singular perturbation method. While the method works for ordinary differential equations and difference equations, it has not been clear if it works or not also for partial differential equations. The straightforward application of the method does not work because, due to the increase of independent variables in the case of partial differential equations, systems have no Lie symmetry group suitable for the construction of asymptotic dynamics.

This paper reports that, by extending the method appropriately, we can reduce partial differential equations to asymptotic dynamics. The main characteristic of the extension is to add some constraint equations to the equations originally considered in order to find partial Lie symmetry groups. The extended method is applied to one ordinary differential equation and two partial differential equations in the paper. While examples taken up here are quite simple in order to make it easy to catch the essence of the procedure, this method works for all the same types of singular perturbation problems in which secular terms are included in the naïve expansion of the solution.

2. Outline of the Extended Method

The singular perturbation method with the Lie group analysis [31–33] is briefly reviewed. For the perturbed dynamical system under consideration here, referred to as the original system in what follows, we find an approximate Lie symmetry group which includes the transformation of the perturbation parameter that is a constant in terms of the dynamics. An approximate solution is generated from the nonperturbative solution by the transformation of the group; in other words, the approximate solution is given by the group-invariant solution. The differential equation, which the group-invariant solution satisfies, provides a reduced system describing the asymptotic behavior of the solution when we consider its long-time dynamics. It should be noted that another method presented by Gazizov and Khalique [19] also considers transformations of the perturbation parameter while the type of the transformation is different from that presented here. With their method, the regular expansion of the solution of the original system is successfully derived.

The method is extended in this study. The procedure is briefly summarized as follows. (1) Specify beforehand a nonperturbative solution on which we would like to see the effect of the perturbation. (2) Construct equations composed of the independent and dependent variables and their derivatives based on the nonperturbative solution specified in step (1). Those equations are referred to as constraint equations in what follows. (3) Find an approximate Lie symmetry group of the original equations on the manifold formed by the constraint equations constructed in the previous step. (4) For the approximate Lie symmetry group found in step (3), construct the differential equation that is satisfied by the group-invariant solution and make an appropriate approximation when necessary. As a consequence, this differential equation generates the reduced system, which describes the asymptotic behavior of the original system.

The procedure is almost the same as that presented before for ordinary differential equations and difference equations [31–33]. The main difference is in the second step, that is to say, the addition of the constraint equations for finding an approximate Lie symmetry group. In general, additional Lie symmetry groups can be found by the addition of such a constraint equation, because the groups do not need to be admitted by the entire manifold formed by the original system in the jet space but need to be admitted only by that part of the manifold restricted by the added equations. In this sense, a Lie symmetry group found by using this procedure is interpreted as a partial Lie symmetry group [24]. Because of the increase in the number of the admitted Lie symmetry groups, in some cases, this method allows us to find a Lie symmetry group when it is impossible to find such a group with the normal Lie group analysis.

3. Application to Ordinary Differential Equation

The method presented in this study is particularly meaningful when it is applied to partial differential equations. However, first, we apply the method to an ordinary differential equation, because the calculation is simpler and therefore serves to clarify the essence of the method.

Consider the following ordinary differential equation [34]:

Here, denotes the dependent variable, the overdot denotes the derivative with respect to the independent variable , and denotes a perturbation parameter supposed to be small sufficiently.

The two independent solutions of the nonperturbative equation, , are

where and are constants.

Focusing on one of these solutions, , we investigate the effect of the perturbation on this nonperturbative solution. For this solution,

is satisfied. From this relation, an equation is constructed,

which is used as a constraint equation in finding a Lie symmetry group.

We find a Lie symmetry group admitted by (1). Let

be the infinitesimal generator of a Lie symmetry group. Its second prolongation is given by

When we adopt the normal procedure followed by Lie group analysis [12], we find a Lie symmetry group by solving the differential equation for ,

which indicates the criterion for the invariance of the entire manifold formed by the original differential equation in the jet space. In the extended method presented in this study, instead, we find a Lie symmetry group by solving

Note that comes from (4). Then, (8) reads

Noting the homogeneity of the equation, we find an exact solution of this differential equation,

For this Lie symmetry group, the group-invariant solution, , satisfies

which reads

This equation, which can be analytically solved, has the following solution:

Here, we have used , which is the nonperturbative solution, as the boundary condition in the integration. This solution corresponds exactly to the solution derived with the renormalization group method [34].

The original differential equation (1) admits no exact Lie symmetry group when we adopt the normal procedure (see the Appendix). However, the addition of constraint equation (4) has allowed us to find a Lie symmetry group.

4. Application to Partial Differential Equation

In this section the present method is applied to two partial differential equations. In general, in the case of partial differential equations, because of the increase of independent variables, it often happens that a partial differential equation admits no approximate Lie symmetry group suitable for the derivation of a system describing the asymptotic dynamics. However, it is in some cases possible to find an approximate Lie symmetry group suitable for deriving the asymptotic dynamic properties by adding some constraint equations constructed from the nonperturbative solution, thereby restricting the manifold on which we find a Lie symmetry group. The reduced systems obtained here are consistent with those derived with the singular perturbation methods presented before.

4.1. First Example:

Consider a system of weakly nonlinear first-order partial differential equations,

Here, is a real vector that has components and depends on two independent variables and . is an nondegenerate matrix that depends on . is a perturbation parameter. The subscripts of denote the derivative with respect to the independent variables.

Similar to the approach followed in the previous section, we find a Lie symmetry group such that it is not admitted by the entire manifold formed by the original equation; rather, it is admitted by a part of the manifold determined from a solution of the nonperturbative system. In the following, we focus on such a nonperturbative solution, , that is an arbitrary function of where is an eigenvalue of ; namely, we set

This solution represents a traveling wave of which the velocity is . As shown by the direct substitution into the nonperturbative system, is an eigenvector of :

In addition,

is also satisfied. Based on these relations, we find a Lie symmetry group with constraints of

in the following.

Let

be the infinitesimal generator of a Lie symmetry group. Its first prolongation is given by

The criterion for the invariance of the original differential equation (14) on the manifold restricted by the constraint equations (18) is given by

This equation reads

Here, and are defined by matrices of which the components are given by, respectively,

where denotes the component of the matrix, .

Suppose and can be expanded as a power series of . Their expanded forms,

are substituted into (22), which is then rewritten as

We consider only the lowest order of , namely, (25), which is rewritten as

and reads

The vector on the right-hand side of this equation is separated into the following two components:

where denotes the unit vector with the same direction as in and is the unit vector parallel to the projection of on a dimensional subspace in spanned by all the eigenvectors of except . The solution corresponding to the second term is given by a linear combination of the eigenvectors of whose linear coefficients are a function of . On the other hand, a term proportional to appears in the solution corresponding to the first term, because every vector proportional to whose proportional coefficient is a function of is a zero eigenvector of the linear operator . For this reason, we can write a solution of (29) as follows:

The group-invariant solution, , for this approximate Lie symmetry group satisfies

which reads

The long-time behavior of the solution is obtained by considering the case in which and by neglecting the terms other than the first one. Then, (32) is reduced to

where is a long time-scale defined by . As a result, we have obtained a dynamical system which describes asymptotic behavior of the original system. This reduced dynamics corresponds to that presented before with the reductive perturbation method [3] or the renormalization group method [7].

4.2. Second Example:

We consider a perturbed Klein-Gordon equation represented by

Here, is a scalar function which depends on two independent variables, and , and is a perturbation parameter. The subscripts of denote the derivative with respect to the independent variables. By introducing new dependent variables as

we rewrite (34) in terms of four first-order partial differential equations for three dependent variables , , and as follows:

The last equation implies the integrable condition.

As in the previous examples, a specific nonperturbative solution is prespecified, on which we investigate the effect of the perturbation. Here we focus on a nonperturbative solution describing a traveling wave with a constant velocity, . The solution is represented by

where for and an arbitrary constant, . For this solution,

therefore,

are satisfied. By combining these relations with the nonperturbative system, , we obtain

From these relations, we construct the constraint equations,

which are used in the finding of Lie symmetry group.

Let

be the infinitesimal generator of a Lie symmetry group of the four partial differential equations (36). Its first prolongation is given by

where , , , and are defined as follows:

The criteria for the invariance of the original system on the manifold formed by the constraint equations (41) are given by

With the identities,

(45) read, in the lowest order of ,

Noting the homogeneity of the equation, we find a solution,

where is an arbitrary constant. For this infinitesimal generator of the approximate Lie symmetry group, the group-invariant solution satisfies

which reads

Asymptotic dynamic behavior can be obtained by considering the cases and and by neglecting the terms which are not proportional to or . By introducing and , we obtain

By eliminating , the equations are reduced to

This reduced equation describes the asymptotic behavior of the original system. This fact can be confirmed by comparing with the result derived before. By introducing a transformation of the variables,

the reduced equation (52) is transformed into

This is the nonlinear Schrödinger equation, which corresponds to the reduced equation derived with the renormalization group method [9].

5. Concluding Remarks

In summary, this paper presents an investigation of perturbation problems, which involves the application of a particular type of Lie group method to partial as well as ordinary differential equations. As a consequence, we have obtained the dynamics which describes the asymptotic behavior of the original system. The main characteristic of the present method is the addition of some constraint equations when searching for a Lie symmetry group, an approach which allows us to find a partial Lie symmetry group. The constraint equations are determined from the nonperturbative solution on which we try to investigate the effect of the perturbation. The reason why it is necessary to specify a nonperturbative solution beforehand for the partial differential equations is that, unlike for ordinary differential equations, the nonperturbative system has an infinite number of independent solutions. In terms of technique, those constraint equations are constructed in such a way that the arbitrary constants or functions included in the nonperturbative solution are eliminated as we have done in deriving (4), (18), and (41). In this study, we have investigated the effect of the perturbation of its lowest order; however, extending our work to higher-order approximations would require us to use higher-order constraint equations instead.

Examples to which the method is applied here are singular perturbation problems in which secular terms appear in the naïve expansion of the solution [3, 6, 8, 10, 11, 34]. Terms proportional to the independent variables emerge in the infinitesimal generator of the Lie symmetry group as seen in (10), (30), and (48) corresponding to the emergence of the secular terms in the naïve expansion. As inferred from this mathematical similarity, the present method is also applicable to other systems which exhibit the same type of singular perturbation problems. In the singular perturbation methods proposed thus far, it is only after the construction of the naïve expansion that it becomes clear whether the method should be used or not. However, the present method can be applied without making the expansion, enabling the appropriate reduced equations to be derived.

Future work should involve a clarification of the relations between the present method and other types of singular perturbation methods presented thus far, such as the WKB or boundary layer problem methods. In addition, there are existing methods that have enabled the successful derivation of the asymptotic behavior for various other partial differential equations [13–18, 20–30, 35–38]. Therefore, it should be clarified whether or not the present method can be appropriately extended and applied to the problems those methods have resolved.

Appendix

Searching for Lie Symmetry with the Normal Procedure

Let us search for the Lie symmetry group of (1) without constraint equations. The criterion for the invariance of the equation is given by (7); namely,

which reads

In the lowest order of , has to satisfy the following differential equations:

However, we can easily find out that there is no solution that satisfies all the above equations, which implies that there is no Lie symmetry group.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by JSPS KAKENHI Grant no. 26800208.

References

N. N. Bogoliubov and C. M. Place, Asymptotic Methods in the Theory of Non-Linear Oscillations, Gordon and Break, 1961.

G. Hori, “Theory of general perturbation with unspecified canonical variable,” Publications of the Astronomical Society of Japan, vol. 18, p. 287, 1966.

T. Taniuchi, “Reductive perturbation method and far fields of wave equations,” Progress of Theoretical Physics Supplement, vol. 55, pp. 1–35, 1974.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical systems, and Bifurcations of Vector Fields, Springer, New York, NY, USA, 1983.

A. H. Nayfeh, Method of Normal Forms, John Wiley & Sons, New York, NY, USA, 1993.

T. Kunihiro, “A geometrical formulation of the renormalization group method for global analysis,” Progress of Theoretical Physics, vol. 94, no. 4, pp. 503–514, 1995.

L.-Y. Chen, N. Goldenfeld, Y. Oono, and G. Paquette, “Selection, stability and renormalization,” Physica A: Statistical Mechanics and its Applications, vol. 204, no. 1–4, pp. 111–133, 1994.

L. Y. Chen, N. Goldenfeld, and Y. Oono, “Renormalization group and singular perturbations: multiple scales, boundary layers, and reductive perturbation theory,” Physical Review E, vol. 54, no. 1, pp. 376–394, 1996.

S.-I. Goto, Y. Masutomi, and K. Nozaki, “Lie-group approach to perturbative renormalization group method,” Progress of Theoretical Physics, vol. 102, no. 3, pp. 471–497, 1999.

R. E. L. DeVille, A. Harkin, M. Holzer, K. Josić, and T. J. Kaper, “Analysis of a renormalization group method and normal form theory for perturbed ordinary differential equations,” Physica D, vol. 237, no. 8, pp. 1029–1052, 2008.

H. Chiba, “ approximation of vector fields based on the renormalization group method,” SIAM Journal on Applied Dynamical Systems, vol. 7, no. 3, pp. 895–932, 2008.

P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, NY, USA, 1986.

V. A. Baĭkov, R. K. Gazizov, and N. H. Ibragimov, “Approximate symmetries,” Mathematics of the USSR-Sbornik, vol. 64, no. 2, p. 427, 1989.

W. I. Fushchich and W. M. Shtelen, “On approximate symmetry and approximate solutions of the nonlinear wave equation with a small parameter,” Journal of Physics A: Mathematical and General, vol. 22, no. 18, pp. L887–L890, 1989.

N. Eular, M. W. Shulga, and W.-H. Steeb, “Approximate symmetries and approximate solutions for a multidimensional Landau-Ginzburg equation,” Journal of Physics A: Mathematical and General, vol. 25, no. 18, Article ID L1095, 1992.

Y. Y. Bagderina and R. K. Gazizov, “Equivalence of ordinary differential equations ,” Differential Equations, vol. 43, no. 5, pp. 595–604, 2007.

V. N. Grebenev and M. Oberlack, “Approximate Lie symmetries of the Navier-Stokes equations,” Journal of Nonlinear Mathematical Physics, vol. 14, no. 2, pp. 157–163, 2007.

R. K. Gazizov and N. K. Ibragimov, “Approximate symmetries and solutions of the Kompaneets equation,” Journal of Applied Mechanics and Technical Physics, vol. 55, no. 2, pp. 220–224, 2014.

R. K. Gazizov and C. M. Khalique, “Approximate transformations for van der Pol-type equations,” Mathematical Problems in Engineering, vol. 2006, Article ID 68753, 11 pages, 2006.

G. Gaeta, “Asymptotic symmetries and asymptotically symmetric solutions of partial differential equations,” Journal of Physics A: Mathematical and General, vol. 27, no. 2, pp. 437–451, 1994.

G. Gaeta, “Asymptotic symmetries in an optical lattice,” Physical Review A, vol. 72, no. 3, Article ID 033419, 2005.

G. Gaeta, D. Levi, and R. Mancinelli, “Asymptotic symmetries of difference equations on a lattice,” Journal of Nonlinear Mathematical Physics, vol. 12, supplement 2, pp. 137–146, 2005.

G. Gaeta and R. Mancinelli, “Asymptotic scaling symmetries for nonlinear PDEs,” International Journal of Geometric Methods in Modern Physics, vol. 2, no. 6, pp. 1081–1114, 2005.

G. Gaeta and R. Mancinelli, “Asymptotic scaling in a model class of anomalous reaction-diffusion equations,” Journal of Nonlinear Mathematical Physics, vol. 12, no. 4, pp. 550–566, 2005.

D. Levi and M. A. Rodríguez, “Asymptotic symmetries and integrability: the KdV case,” Europhysics Letters, vol. 80, no. 6, Article ID 60005, 2007.

V. F. Kovalev, V. V. Pustovalov, and D. V. Shirkov, “Group analysis and renormgroup symmetries,” Journal of Mathematical Physics, vol. 39, no. 2, pp. 1170–1188, 1998.

V. F. Kovalev, “Renormalization group analysis for singularities in the wave beam self-focusing problem,” Theoretical and Mathematical Physics, vol. 119, no. 3, pp. 719–730, 1999.

V. F. Kovalev and D. V. Shirkov, “Functional self-similarity and renormalization group symmetry in mathematical physics,” Theoretical and Mathematical Physics, vol. 121, no. 1, pp. 1315–1332, 1999.

D. V. Shirkov and V. F. Kovalev, “The Bogoliubov renormalization group and solution symmetry in mathematical physics,” Physics Report, vol. 352, no. 4–6, pp. 219–249, 2001.

V. F. Kovalev and D. V. Shirkov, “The renormalization group symmetry for solution of integral equations,” Proceedings of the Institute of Mathematics of NAS of Ukraine, vol. 50, no. 2, pp. 850–861, 2004.

M. Iwasa and K. Nozaki, “A method to construct asymptotic solutions invariant under the renormalization group,” Progress of Theoretical Physics, vol. 116, no. 4, pp. 605–613, 2006.

M. Iwasa and K. Nozaki, “Renormalization group in difference systems,” Journal of Physics A: Mathematical and Theoretical, vol. 41, no. 8, Article ID 085204, 2008.

H. Chiba and M. Iwasa, “Lie equations for asymptotic solutions of perturbation problems of ordinary differential equations,” Journal of Mathematical Physics, vol. 50, no. 4, 2009.

Y. Y. Yamaguchi and Y. Nambu, “Renormalization group equations and integrability in hamiltonian systems,” Progress of Theoretical Physics, vol. 100, no. 1, pp. 199–204, 1998.

N. D. Goldenfeld, Lectures on Phase Transitions and the Renormalization Group, Addison-Wesley, 1992.

J. Bricmont, A. Kupiainen, and G. Lin, “Renormalization group and asymptotics of solutions of nonlinear parabolic equations,” Communications on Pure and Applied Mathematics, vol. 47, no. 6, pp. 893–922, 1994.

N. V. Antonov and J. Honkonen, “Field-theoretic renormalization group for a nonlinear diffusion equation,” Physical Review E, vol. 66, no. 4, Article ID 046105, 7 pages, 2002.

S. Yoshida and T. Fukui, “Exact renormalization group approach to a nonlinear diffusion equation,” Physical Review E, vol. 72, no. 4, Article ID 046136, 2005.