Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 897912, 7 pages

http://dx.doi.org/10.1155/2013/897912

Research Article

Conservation Laws of Some Physical Models via Symbolic Package GeM

Rehana Naz,¹ Imran Naeem,² and M. Danish Khan²

¹Centre for Mathematics and Statistical Sciences, Lahore School of Economics, Lahore 53200, Pakistan

²Department of Mathematics, School of Science and Engineering, LUMS, Lahore Cantt 54792, Pakistan

Received 10 May 2013; Accepted 13 June 2013

Academic Editor: Chaudry Masood Khalique

Copyright © 2013 Rehana Naz 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.

Abstract.

We study the conservation laws of evolution equation, lubrication models, sinh-Poisson equation, Kaup-Kupershmidt equation, and modified Sawada-Kotera equation. The symbolic software GeM (Cheviakov (2007) and (2010)) is used to derive the multipliers and conservation law fluxes. Software GeM is Maple-based package, and it computes conservation laws by direct method and first homotopy and second homotopy formulas.

1. Introduction

The study of conservation laws plays a vital role in analysis, solution, and reductions of PDEs. For the PDEs, the conservation laws are used in wide variety of applications, for example, inverse scattering transform in soliton solutions [1], bi-Hamiltonian structures and recursion operators [2], Lax operators [3], and derivation of conserved quantities for jet flows [4].

Different methods have been developed so far for the construction of conservation laws and are well documented in [5–7]. In the last few decades, the researchers focused on the development of symbolic computational packages based on different approaches of conservation laws. These packages work with either Mathematica or Maple. The development of symbolic computational packages gives relief to perform complicated and tedious algebraic computation. Recently, several computational packages have been developed, for example, CONDENS.M by Göktaş and Hereman [8], RUDCE by Wolf et al. [9–11], TransPDEDensity.m by Adams and Hereman [12], GeM by Cheviakov [13, 14], Vessiot suite by Anderson and Cheb-Terrab [15], ConservationLawsMD.m by Poole and Hereman [16], and SADE by Rocha Filho and Figueiredo [17].

In this paper, we will use GeM package [13] to compute the conservation laws for partial differential equations (PDEs) arising in applications. GeM package works with Maple to obtain the symmetries and conservation laws of differential equations. In symmetry analysis, it first computes the overdetermined system of determining equations and then simplifies the system by Rif package routines. After simplification, a Maple command in GeM generates all symmetry generators of differential equation. In conservation laws analysis, GeM computes an overdetermined system of determining equation of conservation law multipliers, and then this system is simplified by Rif package which is solved by using the built-in Maple function pdsolve to get multipliers. After computing multipliers, the conservation laws fluxes are derived by one of the following four methods: direct method [18, 19], first homotopy formula [20], second homotopy formula [19], and scaling symmetry formula [21]. All these four methods have some limitations in their use. The direct method written in GeM [13] is a Maple implementation based on Wolf [11] program in REDUCE. For simple partial differential equation (PDE) systems and multipliers, direct method is used to calculate fluxes. It is also used if arbitrary functions are involved. The conservation laws fluxes for complicated PDEs or multipliers, not involving arbitrary functions, are established by using first and second homotopy formulas. The scaling symmetry method is used to compute fluxes for the scaling-homogeneous PDEs or/and multipliers. For the complicated scaling-homogeneous PDEs and/or multipliers involving arbitrary functions, this is only a systematic method for computing fluxes.

The evolution equations are important and arise in many applications. We compute the conservation laws of various nonlinear evolution equations using GeM Maple routines. This includes a -dimensional evolution equation [22], lubrication models [23], sinh-Poisson equation [24], Kaup-Kupershmidt equation [25], and modified Sawada-Kotera equation [26]. At last, we summarize and discuss our results.

2. Multipliers and Conservation Laws Using GeM Maple Routines

2.1. Evolution Equation

As a first example, consider the following evolution equation [22]:

where and ,  ,   are constants. We will explain this example in detail along with GeM Maple routines given in [13, 14]. The variables and partial differential equation (PDE) (1) are defined in GeM by the following Maple commands.

 With(GeM): gem_decl_vars(indeps=[t,x], deps=[u(t,x)]); gem_decl_eqs([diff(u(t,x),t,t)+a*diff(u(t,x),x,x)+b*u(t,x)+c*(t,x)=0], solve_for=[diff(u(t,x),t,t)]).

The option solve_for is used in the flux-computation routine, and actually it defines a set of leading derivatives the given PDE systems can be solved for.

Consider multipliers of the form . In GeM, we use the Maple routines,

 det_eqs:=gem_conslaw_det_eqs([t,x,u(t,x),diff(u(t,x),t),diff(u(t,x),x)]): CL_multipliers:=gem_conslaw_multipliers(); simplified_eqs:=DEtools[rifsimp](det_eqs, CL_multipliers, mindim=1),

to obtain the set of determining equations for the multipliers expressed in the simplified form as

To solve the system (2), we use the Maple command

 multipliers_sol:=pdsolve(simplified_eqs[Solved]),

and it yields

where ,  ,   are arbitrary constants. We obtain three linearly independent conservation laws, arising from the multipliers

Next step is the derivation of conservation laws associated with multipliers given in (4). The Maple command

 gem_get_CL_fluxes(multipliers_sol)

computes the flux expressions by the direct method. For the multipliers (4), we have the following conservation laws fluxes:

The multipliers given in (4) do not involve arbitrary functions, so homotopy formulas can be used to compute fluxes. We call the routine for first homotopy method

 gem_get_CL_fluxes(multipliers_sol, method=“Homotopy1”)

to get the following expressions for conservation law fluxes:

For second homotopy formula, the Maple command

 gem_get_CL_fluxes(multipliers_sol, method=“Homotopy2”)

yields divergence expressions in the same form as in (6).

The PDE (1) has no scaling symmetry; therefore, we cannot apply the scaling symmetry formula here for derivation of fluxes.

2.2. Lubrication Models

Now we will study two lubrication models for conservation laws point of view. Gandarias and Medina [23] performed the symmetry analysis of lubrication model

where is an arbitrary function. For and , this equation has some extra symmetry [23]. Without loss of generality, take in (7); we have

where is arbitrary constant. Consider the multipliers of form in GeM Maple routines, and then we obtain the following four multipliers:

The fluxes associated with the multipliers given in (7) are computed by homotopy first method and are given by

We will get the same fluxes for (7) if we define higher order multipliers in GeM Maple routines.

Another interesting lubrication model is

It is obtained by taking in (7). The GeM Maple routines yield the following four multipliers of form :

The corresponding fluxes obtained by homotopy first method are

The conservation laws fluxes derived here can be used to find the solution of lubrication models and will be considered in future work.

2.3. sinh-Poisson Equation

The -dimensional sinh-Poisson equation is [24]

where . The conservation laws for PDE (14) are derived here by using GeM routines. Consider the multipliers of form in GeM routines, then it will yield six multipliers not containing any arbitrary function. The expression for fluxes is computed by using first homotopy formula. The multipliers and associated conserved vectors computed by first homotopy formula are given in Table 1.

Table 1: Multipliers and conserved vectors for PDE (14).

Multiplier
Fluxes

2.4. Kaup-Kupershmidt Equation

Now, we will compute the conservation laws for the fifth order Kaup-Kupershmidt [25]:

The GeM Maple routines yield three multipliers of the form for PDE (14). The first homotopy formula is applied to derive the expressions for conservation laws fluxes. One more multiplier can be computed if we consider higher order multipliers of the form . All the multipliers and associated conserved vectors for PDE (14) computed by first homotopy formula are presented in Table 2.

Table 2: Multipliers and conserved vectors for PDE (15).

Multiplier
Fluxes

2.5. Modified Sawada-Kotera Equation

Consider the fifth order modified SK equation:

For PDE (16), two conserved densities were derived by first computing Lax pair (see [26]). The higher order conservation laws fluxes exist for higher order multipliers and are not reported in [26]. Consider the multipliers of form in GeM routines, then it will yield two simple and three higher order multipliers not containing any arbitrary function. The simple multipliers yield same fluxes as derived in [26], and three new fluxes corresponding to higher order multipliers are computed. The multipliers and associated conserved vectors computed by first homotopy formula are listed in Table 3.

Table 3: Multipliers and conserved vectors for PDE (16).

Multiplier
Fluxes

,  

3. Conclusions

The conservation laws for the evolution equation, Benjamin equation, lubrication models, sinh-Poisson equation, Kaup-Kupershmidt equation, and modified Sawada-Kotera equation were derived by using the symbolic software GeM. First of all, we considered the evolution equation, and the commands for all GeM Maple routines, were explicitly given. The first order multipliers were defined in GeM Maple routines and three multipliers were obtained. The expressions for fluxes were computed by direct method and first and second homotopy formulas and equivalent expressions for fluxes were obtained. The scaling symmetry method was not applicable here as no scaling symmetry exists for the nonlinear evolution equation. The conservation laws fluxes for the lubrication models, sinh-Poisson equation, Kaup-Kupershmidt equation, and modified Sawada-Kotera equation were derived by the first homotopy formula. For the modified Sawada-Kotera equation, three new fluxes were derived.

The fluxes derived here can be used in constructing the solutions of underlying PDEs and will be considered in the future work.

Conflict of Interests

The authors declare that there is no conflict of interests.

References

M. J. Ablowitz and P. A. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse Scattering, vol. 149 of London Mathematical Society Lecture Note Series, Cambridge University Press, Cambridge, 1991.

D. E. Baldwin and W. Hereman, “A symbolic algorithm for computing recursion operators of nonlinear partial differential equations,” International Journal of Computer Mathematics, vol. 87, no. 5, pp. 1094–1119, 2010.

V. Drinfel’d and V. Sokolov, “Lie algebras and equations of Korteweg-de vries type,” Journal of Soviet Mathematics, vol. 30, no. 2, pp. 1975–2036, 1985.

R. Naz, D. P. Mason, and F. M. Mahomed, “Conservation laws and conserved quantities for laminar two-dimensional and radial jets,” Nonlinear Analysis, vol. 10, no. 5, pp. 2641–2651, 2009.

G. W. Bluman, A. F. Cheviakov, and S. C. Anco, Applications of Symmetry Methods to Partial Differential Equations, vol. 168 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1st edition, 2010.

R. Naz, Symmetry solutions and conservation laws for some partial differ-ential equations in fluid mechanics [Ph.D. dissertation], University of the Witwatersrand, Johannesburg, South Africa, 2008.

R. Naz, F. M. Mahomed, and D. P. Mason, “Comparison of different approaches to conservation laws for some partial differential equations in fluid mechanics,” Applied Mathematics and Computation, vol. 205, no. 1, pp. 212–230, 2008.

Ü. Göktaş and W. Hereman, “Symbolic computation of conserved densities for systems of nonlinear evolution equations,” Journal of Symbolic Computation, vol. 24, no. 5, pp. 591–621, 1997.

T. Wolf, A. Brand, and M. Mohammadzadeh, “Computer algebra algorithms and routines for the computation of conservation laws and fixing of gauge in differential expressions,” Journal of Symbolic Computation, vol. 27, no. 2, pp. 221–238, 1999.

T. Wolf, “A comparison of four approaches to the calculation of conservation laws,” European Journal of Applied Mathematics, vol. 13, no. 2, pp. 129–152, 2002.

T. Wolf, “Crack, LiePDE, ApplySym and ConLaw. section 4.3.5 and computer programon CD-ROM,” in Computer Algebra Handbook, J. Grabmeier, E. Kaltofen, and V. Weispfenning, Eds., pp. 465–468, Springer, Berlin, Germany, 2002.

P. J. Adams and W. Hereman, TransPDEDensityFlux.m: Symbolic computation of conserved densities and fluxes for systems of partial differential equations with transcendental nonlinearities, Scientific Software, 2002, http://inside.mines.edu/

A. F. Cheviakov, “GeM software package for computation of symmetries and conservation laws of differential equations,” Computer Physics Communications, vol. 176, no. 1, pp. 48–61, 2007.

A. F. Cheviakov, “Computation of fluxes of conservation laws,” Journal of Engineering Mathematics, vol. 66, no. 1–3, pp. 153–173, 2010.

I. M. Anderson and E. S. Cheb-Terrab, Differential geometry package, Maple Online Help, 2009, http://www.maplesoft.com/support/help/view.aspx?sid=26040

L. D. Poole and W. Hereman, “ConservationLawsMD.m: a mathematica package for the symbolic computation of conservation laws of polynomial systems of nonlinear PDEs in multiple space dimensions,” scientific software, 2009, http://inside.mines.edu/~whereman/

T. M. Rocha Filho and A. Figueiredo, “[SADE] a Maple package for the symmetry analysis of differential equations,” Computer Physics Communications, vol. 182, no. 2, pp. 467–476, 2011.

S. C. Anco and G. Bluman, “Direct construction of conservation laws from field equations,” Physical Review Letters, vol. 78, no. 15, pp. 2869–2873, 1997.

S. C. Anco and G. Bluman, “Direct construction method for conservation laws of partial differential equations. II: general treatment,” European Journal of Applied Mathematics, vol. 13, no. 5, pp. 567–585, 2002.

W. Hereman, M. Colagrosso, R. Sayers et al., “Continuous and discrete homotopy operators and the computation of conservation laws,” in Differential Equations with Symbolic Computation, D. Wang and Z. Zheng, Eds., pp. 249–285, Birkhäuser, 2005.

S. C. Anco, “Conservation laws of scaling-invariant field equations,” Journal of Physics, vol. 36, no. 32, pp. 8623–8638, 2003.

C. Bai, “Exact solutions for nonlinear partial differential equation: a new approach,” Physics Letters A, vol. 288, no. 3-4, pp. 191–195, 2001.

M. L. Gandarias and E. Medina, “Analysis of a lubrication model through symmetry reductions,” Europhysics Letters, vol. 55, no. 2, pp. 143–149, 2001.

A. H. Khater, W. Malfliet, D. K. Callebaut, and E. S. Kamel, “Travelling wave solutions of some classes of nonlinear evolution equations in and dimensions,” Journal of Computational and Applied Mathematics, vol. 140, no. 1-2, pp. 469–477, 2002.

V. P. Gomes Neto, “Fifth-order evolution equations describing pseudospherical surfaces,” Journal of Differential Equations, vol. 249, no. 11, pp. 2822–2865, 2010.

H. Guoliang and G. Xian-Guo, “An extension of the modified Sawada-Kotera equation and conservation laws,” Chinese Physics B, vol. 21, no. 7, Article ID 070205, 2012.