About this Journal Submit a Manuscript Table of Contents
Discrete Dynamics in Nature and Society
Volume 2013 (2013), Article ID 146976, 7 pages
Research Article

Invariant Solutions for Nonhomogeneous Discrete Diffusion Equation

1Department of Mathematics, University of Azad Jammu & Kashmir, Muzaffarabad 13100, Pakistan
2Department of Mathematics, COMSATS Institute, Abbottabad 22010, Pakistan

Received 29 April 2013; Accepted 18 August 2013

Academic Editor: R. Sahadevan

Copyright © 2013 M. N. Qureshi et al. 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.


One-dimensional optimal systems for nonhomogeneous discrete heat equation with different source terms are calculated. By utilizing these optimal systems invariant solutions are found. Also generating solutions are calculated, using the elements of the symmetry algebra.

1. Introduction

Mathematics provides models, for real life phenomena, to precisely understand the underlying laws governing them. Dynamic changes in a physical phenomenon are usually modeled by differential equations. This means that it is supposed that the changes are continuous and these are taking place in a continuous domain. But there are many real life situations where either the domain or the phenomenon itself or the both are not continuous. These situations occur but are not limited to the fields of biology, physics (classical and quantum), geometry, mathematical design, finance and so forth. Discrete situations are better modeled by difference equations, in contrast with the differential equations. Therefore, difference equations make their appearance in almost every branch of mathematics and their importance cannot be overemphasized.

Differences and their calculus are as old as is the differential calculus, but apart from its importance and usefulness the theory and methods of solving difference equations are not as developed as those of differential equations. Since differential equation is a limiting case of difference equation; therefore, it is natural to extend methods available for solving differential equations to the respective difference equations, and the same is done most of the times [13].

Solving differential equations by exploiting symmetry properties of their solution spaces is one of the standard methods, introduced by the Norwegian mathematician Sophus Lie (1842–1899). In contrast with adhoc techniques this method provides an algorithmic and unified procedure to solve almost all types of differential equations, linear or nonlinear [4, 5].

At present, the theory of difference equations is at the same stage as was the theory of differential equations at the time of S. Lie. A number of adhoc approaches are available to solve difference equations, for example, substitution, nonlinear functional relation, Schroder’s generation function, Maeda’s method, or the theory of integrable maps. But a number of authors have attempted, following Lie’s method for differential equations, to develop a unified integration procedure based on invariance properties of difference equation [6, 7]. It is worth noting that the symmetry properties are not only helpful to find solution but play an important role to understand the physical phenomenon more deeply. In this direction, a number of researchers have attempted to apply symmetries to analyze different discrete physical phenomena [619].

The discrete diffusion equation is widely used in many contexts [20, 21]. For instance, it has been applied to the area of population growth where one wishes to model geographic spread in addition to growth in number. In the area of physics it is used to model ionic diffusion on a lattice. It has also been used in digital filtering in the form of diffusion filtering [10]. Due to its wide utility, a number of researchers have analyzed different aspects of the discrete diffusion process. Recently, in [6], Levi et al. have attempted to generalize the Lie infinitesimal formulism for calculation of symmetries of difference equations and calculated generalized symmetries of diffusion type difference equations. The symmetry group, which they have found for the discrete equation, is a minimal extension of the Lie point symmetries of the corresponding differential equation. They also found that in the case of linear discrete equation the symmetry algebra is isomorphic to the continuous limit.

Symmetries are used not only to integrate but also to analyze the solution space for more insight into the physical system at hand. So, once symmetries are given it is natural to ask for further analysis and find solutions of the equation. There is limited literature on the analysis of solutions of diffusion type difference equations via symmetries. In this paper one-dimensional optimal systems are obtained for the symmetry algebra of diffusion difference equations with different source terms found in [6]. We have also calculated invariant solutions corresponding to each representative of the one-dimensional optimal system, wherever they exist. Using elements of the symmetry algebra generating solutions are also obtained. Moreover, all those cases in which no group invariant solution exist are mentioned.

The rest of the paper is arranged as follows. In the next two sections we give some preliminaries and summarize the results of [6], for completeness. In Section 3 one-dimensional optimal systems are obtained for all cases and corresponding group invariant solutions are given. Section 4 is devoted to obtain generating solutions corresponding to each symmetry generator. A brief conclusion is given in Section 5.

2. Preliminaries

Let be a continuous scalar function of independent variables evaluated at finitely many points on a uniformly distributed lattice with spacing . Then partial finite differences are defined by It is sometimes convenient to express finite difference operator in terms of shift operator defined as follows: Consequently, For example, for the factorial function (2) gives Also the action of on is given by

Consider a partial difference equation in the following form: where the operator operates on a function as follows: with , , and , are some fixed integers. The shift operators , are defined likewise.

A discrete vector field in its evolutionary form, with characteristic given by is an infinitesimal generalized symmetry generator of the difference equation (6) if it satisfies the infinitesimal symmetry criterion:

Here, is the discrete analog of the prolongation of the vector field . Prolonged vector field acts on the extended lattice and therefore has the following form: Here, the coefficient functions , are the discrete total variations of given by where the total variation acts on functions of as follows:

3. Symmetry Algebra of Nonhomogeneous Discrete Heat Equation

In this section we summarize the results of [6]. The difference equation is a discrete heat or diffusion equation. Here is a function of space variable , time and partial shifts and with respect to and , respectively. Let the following be the generalized symmetry generator in the evolutionary form: If (14) is admitted by (13) then the infinitesimal symmetry criterion gives where is the second prolongation of . Using (16) in (15) one gets the following set of determining equations:

To obtain the symmetry generator (14) one needs to solve the set of equations (17) to (20) for unknown functions , , and . In [6] these functions have been found for three special source terms . Here, we give the infinitesimal symmetry generators obtained in [6] and refer to the paper for details.

3.1. Free Heat Equation

In this case and (13) takes the form

The set of equations (17) to (20) are then solved to obtain the following symmetry algebra:

3.2. Nonzero Potential

In this case the potential is taken to be a nonzero function . Plugging in this assumption on in the set of equations (17) to (20) one gets where is an arbitrary function to be known. And where is an arbitrary function to be known. To find the unknown functions and one puts a condition on the potential function such that the commutator vanishes and (24) becomes Using (25) in (20) and simplifying one gets

The authors then consider the following two particular physically important models that obey the assumptions imposed on the function .

3.2.1. Discrete Harmonic Oscillator

The discrete harmonic oscillator follows the discrete heat equation with the following potential function: Solving the set of determining equations (17) to (20) for above value of one gets the following symmetry algebra:

3.2.2. Discrete Centrifugal Barrier

Likewise, the discrete centrifugal barrier with potential function has the following symmetry algebra:

We consider the following case where the potential function also depends on the variables , and yet satisfies the required condition.

3.2.3. Time-Dependent Discrete Harmonic Oscillator

Now, if the potential function considered in (27) depends on all the four variables , , , and and is given by the following

then the values of coefficient functions of infinitesimal generator (14) are given by: And one gets 3-dimensional subalgebra of the discrete harmonic oscillator obtained in ((28) and (29)). This is generated by the following infinitesimal generators:

4. One-Dimensional Optimal System and Invariant Solution

In this section one-dimensional optimal system for each of the above cases of discrete diffusion equation is calculated. It is done in analogy with the procedure to calculate optimal systems for differential equation, laid down in [4]. These optimal systems are then used to calculate invariant solutions.

4.1. Free Heat Equation

A calculation gives eight equivalence classes of one-dimensional optimal systems represented by the following generators: , , , , , , , and .

4.1.1. Invariant Solution

Employing the method given in [4] for calculating invariant solution for differential equations, the invariant solution corresponding to the representatives of the optimal systems given above are obtained below.

To obtain invariant solution corresponding to the subalgebra we find the group invariants and define new variables in terms of these invariants by , . This change of variables reduces the partial difference equation (21) to the ordinary 2nd order difference equation . This equation is then solved to get the invariant solution . Similarly, corresponding to , is the group invariant. Substituting , (21) reduces to the ordinary difference equation . A further change of variable reduces the equation to the following first-order difference equation with variable coefficient:

To solve (35) we use Laplace method given in [3]. Let , , and . On substituting these values (35) becomes

Now there are two possibilities either and this gives or ; therefore, or . So, , and one gets the following invariant solution: Similar calculations give invariant solutions corresponding to other representatives of optimal classes. The results are summarized in Table 1.

Table 1: Optimal systems for free heat equation and respective invariant solutions.
4.2. Discrete Harmonic Oscillators and Centrifugal Barrier

For the discrete harmonic oscillator and the discrete centrifugal barrier the representatives of their one-dimensional optimal systems and the corresponding invariant solutions are summarized in Tables 2 and 3.

Table 2: The optimal systems for discrete harmonic oscillator and respective invariant solutions.
Table 3: The optimal systems for discrete centrifugal barrier and respective invariant solutions.

Remark 1. The solution corresponding to is obtained by applying Laplace method in reduced form.

Since, in case of time-dependent harmonic oscillator, the symmetry algebra is a subalgebra of the time-independent case therefore one does not get new invariant solutions.

5. Generating Solutions

Since a symmetry of a given equation maps its solution space onto itself, so, one can find new solutions by applying symmetry transformations to a known solution of the equation. The solution thus obtained is called a generating solution. Here, these generating solutions are obtained by first obtaining the one parameter transformation corresponding to each infinitesimal generators obtained in Section 2. These solutions are obtained in analogy with the method given in [4] for differential equations. The results are summarized in the form of Tables 4, 5, and 6.

Table 4: Generating solutions for free heat equation.
Table 5: Generating solutions for discrete harmonic oscillator.
Table 6: Generating solutions for discrete centrifugal barrier.

6. Conclusion

In this paper, in analogy with differential equations, the symmetry analysis of a nonhomogeneous discrete heat equation (13) has been carried out. The one-dimensional optimal systems are obtained for the equation with different source terms by using symmetry algebra given in [6]. These optimal systems are then used to obtain invariant solutions if they exist. Moreover, one parameter group of transformations corresponding to all infinitesimal generators of the symmetry algebra has been calculated to obtain generating solutions of the equation. Main results are summarized in Tables 1, 2, 3, 4, 5, and 6. Working on the same lines one can calculate conservation laws which we will report soon.


The authors thank the main editor and anonymous referees for their valuable comments and suggestions leading to improvement of this paper. This work was supported by the Higher Education Commission of Pakistan.


  1. S. Goldberg, Introduction to Difference Equations, with Illustrative Examples from Economics, Psychology, and Sociology, John Wiley & Sons, New York, NY, USA, 1958. View at Zentralblatt MATH · View at MathSciNet
  2. C. Jordan, Calculus of Finite Differences, Chelsea Publishing CompanyBronx, New York, NY, USA, 3rd edition, 1965.
  3. R. Bellman and K. L. Cooke, Differential-Difference Equations, Academic Press, New York, NY, USA, 1963. View at Zentralblatt MATH · View at MathSciNet
  4. P. J. Olver, Applications of Lie Groups to Differential Equations, Springer, New York, NY ,USA, 1986. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. G. W. Bluman and S. Kumei, Symmetries and Differential Equations, Springer, New York, NY, USA, 1989. View at Zentralblatt MATH · View at MathSciNet
  6. D. Levi, L. Vinet, and P. Winternitz, “Lie group formalism for difference equations,” Journal of Physics A, vol. 30, no. 2, pp. 633–649, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. S. Maeda, “The similarity method for difference equations,” IMA Journal of Applied Mathematics, vol. 38, no. 2, pp. 129–134, 1987. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. G. R. W. Quispel and R. Sahadevan, “Lie symmetries and the integration of difference equations,” Physics Letters A, vol. 184, no. 1, pp. 64–70, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. Z. Jiang, “Lie symmetries and their local determinacy for a class of differential-difference equations,” Physics Letters A, vol. 240, no. 3, pp. 137–143, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. J. Campbell, “The SMM Model as a boundary value problem using the discrete heat equation,” NASA Langley Research Center Hampton Virginia, vol. 72, no. 4, pp. 539–546, 2007.
  11. R. Hernández Heredero, D. Levi, and P. Winternitz, “Symmetries of the discrete Burgers equation,” Journal of Physics A, vol. 32, no. 14, pp. 2685–2695, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. R. Hernandez Heredero, D. Levi, and P. Winternitz, “Point symmetries and generalized symmetries of nonlinear difference equations,” Preprint CRM 2568, CRM, Montréal, Canada, 1998.
  13. R. Floreanini and L. Vinet, “Lie symmetries of finite-difference equations,” Journal of Mathematical Physics, vol. 36, no. 12, pp. 7024–7042, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. P. A. Clarkson and F. W. Nijhoff, Eds., Symmetries and Integrability of Difference Equations, vol. 255 of London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, UK, 1999. View at Publisher · View at Google Scholar · View at MathSciNet
  15. J. M. Burgers, ‘The Nonlinear Diffusion Equation’ Asymptotic Solutions and Statistical Problems, Reidel, Dordrecht, The Netherlands, 1974.
  16. D. Levi and P. Winternitz, “Symmetries and conditional symmetries of differential-difference equations,” Journal of Mathematical Physics, vol. 34, no. 8, pp. 3713–3730, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. D. Levi and P. Winternitz, “Continuous symmetries of discrete equations,” Physics Letters A, vol. 152, no. 7, pp. 335–338, 1991. View at Publisher · View at Google Scholar · View at MathSciNet
  18. D. David, N. Kamran, D. Levi, and P. Winternitz, “Subalgebras of loop algebras and symmetries of the Kadomtsev-Petviashvili equation,” Physical Review Letters, vol. 55, no. 20, pp. 2111–2113, 1985. View at Publisher · View at Google Scholar · View at MathSciNet
  19. R. Hirota, “Nonlinear partial difference equations. V. Nonlinear equations reducible to linear equations,” Journal of Physical Society of Japan, vol. 46, no. 1, pp. 312–319, 1979.
  20. Z. Y. Lu and Y. Takeuchi, “Global asymptotic behavior in single-species discrete diffusion systems,” Journal of Mathematical Biology, vol. 32, no. 1, pp. 67–77, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. G. Fáth, “Propagation failure of traveling waves in a discrete bistable medium,” Physica D. Nonlinear Phenomena, vol. 116, no. 1-2, pp. 176–190, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet