InvariantEquation - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


PDEtools

  

InvariantEquation

  

computes the differential equation invariant under a given symmetry group

 

Calling Sequence

Parameters

Description

Examples

Compatibility

Calling Sequence

InvariantEquation(S, DepVars, 'options'='value')

Parameters

S

-

the infinitesimals of a symmetry generator or the corresponding infinitesimal generator operator, or a list of n of them representing an n-dimensional group of symmetries

DepVars

-

a function or a list of them indicating the dependent variables of the problem

arbitraryfunctionname

-

(optional) the name used to represent the arbitrary function that appears in the returned invariant equation

explicit

-

(optional) to return the invariant equation(s) with the highest derivative(s) isolated on the left-hand-side(s)

order

-

(optional) the order of the differential equation required. The default value is 1.

jetnotation = ...

-

(optional) can be true (same as jetvariables), false (default), jetvariables, jetvariableswithbrackets, jetnumbers or jetODE; to respectively return or not using the different jet notations available

simplifier = ...

-

(optional) indicates the simplifier to use instead of the default simplify/size simplifier

Description

• 

Given a list with the infinitesimals of a symmetry generator or the corresponding infinitesimal generator operator, or a list of n of these lists or operators, possibly representing an n-dimensional symmetry group, InvariantEquation returns an expression Λφ1,φ2,..., that is simultaneously invariant under all of the Gj symmetry generators (1jn) related to the given n infinitesimals. The φk form a complete set of m differential invariants, so they are functions of the independent and dependent variables of the problem (DepVars), as well as their partial derivatives up to order (if not indicated, its default value is 1), and satisfy Gjφk=0. The equation Λφ1,...,φm=0 is thus the most general partial differential equation simultaneously invariant under all the symmetries in S.

• 

InvariantEquation works also with dynamical symmetries, that is, symmetries where the infinitesimals themselves depend on derivatives of the dependent variables of the problem up to (at most) order - 1, where order is the parameter used to indicate the differential order of the desired invariant equation.

• 

Note that in the multidimensional case, depending on the form of the infinitesimals given, the problem may have no solution. Also the φk entering Λφ1,φ2,... are all independent of each other and their number is unique, but one can still rewrite the φk in different ways, by combining invariants to construct equivalent but algebraically different invariants of the same order.

• 

The representation Λφ1,φ2,...=0 for the invariant equation can be transformed into a system of equations by constructing the objects

phik=Lambdaphi1,..,phik1,phik+1,..,phim

  

so that when the φk, isolated on the left-hand-side, is linear in the highest derivatives, of differential order equal to the value of order used by InvariantEquation, the resulting system will be an explicit PDE system linear in the highest derivatives.

• 

You can request InvariantEquation to directly return a similar explicit PDE system of differential order equal to order by passing the optional argument explicit. In that case the system returned:

– 

will contain as many differential equations as the number - say N - of partial derivatives of order equal to order;

– 

all the partial derivatives of order = order will appear on the left-hand-sides isolated;

– 

the right-hand-sides will involve N arbitrary functions of the form _Fn... with n=1..N;

– 

the right-hand-sides will depend on all the partial derivatives of at most one order less than order, including the unknowns themselves and the independent variables of the problem.

• 

This invariant PDE system is returned as an ordered list of equations, where the ordering is dictated by the partial derivatives and is the same implemented in the SortDerivatives command of the PDEtools Library programming routines.

• 

The returned invariant equation is automatically simplified in size. To avoid that, pass the optional argument simplifier = none, or to use a different simplifier use simplifier = ....

• 

By default, the equation returned by InvariantEquation is in function notation. You can change this by specifying the optional argument jetnotation = ... where the value can be any of the jet notations available.

Examples

withPDEtools:

Consider a PDE problem with two independent variables and one dependent variable, ux,t, and consider the list of infinitesimals of a symmetry group.

S_ξx=x,_ξt=1,_ηu=u

S_ξx=x,_ξt=1,_ηu=u

(1)

In the input above you can also enter the symmetry S without labels for the infinitesimals, as in x,1,u, or use the corresponding infinitesimal generator

GInfinitesimalGeneratorS,ux,t

Gf→xxf+tf+uuf

(2)

The equation invariant under the symmetry transformation underlying the infinitesimal generator S is (you can equally use G instead of S)

PDEInvariantEquationS,ux,t,arbitraryfunctionname=Λ

PDEΛlnx+t,ux,tx,xux,t,tux,tx

(3)

This is the explicit form of the PDE system as explained in the Description, in terms of two arbitrary functions f__1 and _F2

InvariantEquationS,ux,t,explicit

xux,t=f__1lnx+t,ux,tx,tux,t=f__2lnx+t,ux,txx

(4)

The invariance of PDE can be verified in various ways, the simplest of which is perhaps to use SymmetryTest, which verifies that S is a symmetry of PDE.

SymmetryTestS,PDE

0

(5)

A less abstract verification would be to explicitly construct the symmetry transformation related to G, for example, in terms of new variables r,s,vr,s.

itrSymmetryTransformationG,ux,t,vr,s

itrr=xⅇ_ε,s=_ε+t,vr,s=ⅇ_εux,t

(6)

trsolveitr,t,x,ux,t

trt=s_ε,x=rⅇ_ε,ux,t=vr,sⅇ_ε

(7)

Now, change the variables in PDE using tr, and recall that the Lie group parameter _ε is real:

dchangetr,PDE,known=Λ,r,s,vr,s

Λlnrⅇ_ε+s_ε,vr,sr,rvr,s,svr,sr

(8)

simplifyassuming_ε::real

Λlnr+s,vr,sr,rvr,s,svr,sr

(9)

The equation above is identical to PDE, that is the meaning of invariant in the context of symmetries. Alternatively, the following is the prolongation to order 1 of the infinitesimal generator G (5.2) (that is, ready to act on functions depending on x,t,ux,t and partial derivatives of ux,t of order 1 at most).

GInfinitesimalGeneratorG,ux,t,prolongation=1,expanded

Gf→xxf+tf+uuf+ututf

(10)

If you apply this differential operator G to PDE, you obtain zero because of the invariance of the PDE under the transformation tr related to G. To apply G, an operator in jet notation, you need to rewrite PDE in the same notation

jetPDEToJetPDE,ux,t

jetPDEΛlnx+t,ux,ux,utx

(11)

normalGjetPDE

0

(12)

Underlying this zero, there is the way PDE is constructed: it is an arbitrary function of the Invariants associated to G, that is, the solutions of the differential operator (you can equally pass S instead of G)

InvariantsG,ux,t

ux,lnx+t,ux,utx

(13)

InvariantEquation also handles dynamical symmetries, that is, symmetries that depend on the partial derivatives of the dependent variables of the problem. In these cases it is however of no use to directly compute the invariants, as in the example above, because they will all be of higher order. The approach used is then different. Consider for instance these infinitesimals depending on xux,t−1

dynamical_S_ξx=1,_ξt=0,_ηu=1ux

dynamical_S_ξx=1,_ξt=0,_ηu=1ux

(14)

As mentioned, the differential invariants, for instance, of order 1, automatically come depending on higher derivatives

Invariantsdynamical_S,ux,t,order=1

t,ux3+3xux,x,2uux,x+ux22ux,x,utux,xuxux,tux,x

(15)

In fact, for order = n, and besides the invariant of order 0, all the other ones will always be of order n+1. To see the actual invariants of order 1 you would need to eliminate the higher derivatives using the related invariant equation, thus defeating the use of these invariants with the purposes of constructing invariant equations of order n. The approach used by InvariantEquation results instead in

InvariantEquationdynamical_S,ux,t,order=2,name=Λ=0

* Partial match of 'name' against keyword 'arbitraryfunctionname'

Λt,xux,t,ux,txux,t22xxux,t+ux,txux,t2+1,tux,t,xux,t3+22x2ux,tx+xux,txux,t2x2ux,txux,t2+1,xux,t2+1xux,t2txux,t2x2ux,t,2x2ux,t2t2ux,t2txux,t22x2ux,t=0

(16)

Note the φj depend on first and second derivatives of ux,t, as it should be. In explicit form (use declare to have a more compact and readable output with derivatives displayed in jet notation)

declareux,t

ux,twill now be displayed asu

(17)

InvariantEquationdynamical_S,ux,t,order=2,explicit

ux,x=uxux2+1ux3f__1t,ux,uux22xux+uux2+1,ut+uxf__1t,ux,uux22xux+uux2+1,ut2x,ut,x=f__2t,ux,ut,uux22xux+uux2+1ux3f__1t,ux,uux22xux+uux2+1,ut+uxf__1t,ux,uux22xux+uux2+1,ut2x,ut,t=f__2t,ux,ut,uux22xux+uux2+12ux2ux2+12f__1t,ux,uux22xux+uux2+1,ut+2uxxux2+1+f__3t,ux,uux22xux+uux2+1,ut

(18)

where the PDE system returned contains as many differential equations as partial derivatives of order 2, that is, 3 equations, and the same number of arbitrary functions _Fn.

Compatibility

• 

The explicit option was introduced in Maple 15.

• 

For more information on Maple 15 changes, see Updates in Maple 15.

See Also

dchange

InfinitesimalGenerator

Invariants

PDEtools

SymmetryTest

SymmetryTransformation

ToJet