PDEtools
InvariantEquation
computes the differential equation invariant under a given symmetry group
Calling Sequence
Parameters
Description
Examples
Compatibility
InvariantEquation(S, DepVars, 'options'='value')
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
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 (1⁢≤⁢j⁢≤⁢n) 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=Lambda⁡phi1,..,phik−1,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.
with⁡PDEtools:
Consider a PDE problem with two independent variables and one dependent variable, u⁡x,t, and consider the list of infinitesimals of a symmetry group.
S≔_ξx=x,_ξt=1,_ηu=u
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
G≔InfinitesimalGenerator⁡S,u⁡x,t
G≔f→x⁢∂∂x⁢f+∂∂t⁢f+u⁢∂∂u⁢f
The equation invariant under the symmetry transformation underlying the infinitesimal generator S is (you can equally use G instead of S)
PDE≔InvariantEquation⁡S,u⁡x,t,arbitraryfunctionname=Λ
PDE≔Λ⁡−ln⁡x+t,u⁡x,tx,∂∂xu⁡x,t,∂∂tu⁡x,tx
This is the explicit form of the PDE system as explained in the Description, in terms of two arbitrary functions f__1 and _F2
InvariantEquation⁡S,u⁡x,t,explicit
∂∂xu⁡x,t=f__1⁡−ln⁡x+t,u⁡x,tx,∂∂tu⁡x,t=f__2⁡−ln⁡x+t,u⁡x,tx⁢x
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.
SymmetryTest⁡S,PDE
0
A less abstract verification would be to explicitly construct the symmetry transformation related to G, for example, in terms of new variables r,s,v⁡r,s.
itr≔SymmetryTransformation⁡G,u⁡x,t,v⁡r,s
itr≔r=x⁢ⅇ_ε,s=_ε+t,v⁡r,s=ⅇ_ε⁢u⁡x,t
tr≔solve⁡itr,t,x,u⁡x,t
tr≔t=s−_ε,x=rⅇ_ε,u⁡x,t=v⁡r,sⅇ_ε
Now, change the variables in PDE using tr, and recall that the Lie group parameter _ε is real:
dchange⁡tr,PDE,known=Λ,r,s,v⁡r,s
Λ⁡−ln⁡rⅇ_ε+s−_ε,v⁡r,sr,∂∂rv⁡r,s,∂∂sv⁡r,sr
simplify⁡assuming_ε::real
Λ⁡−ln⁡r+s,v⁡r,sr,∂∂rv⁡r,s,∂∂sv⁡r,sr
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,u⁡x,t and partial derivatives of u⁡x,t of order 1 at most).
G≔InfinitesimalGenerator⁡G,u⁡x,t,prolongation=1,expanded
G≔f→x⁢∂∂x⁢f+∂∂t⁢f+u⁢∂∂u⁢f+ut⁢∂∂ut⁢f
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
jetPDE≔ToJet⁡PDE,u⁡x,t
jetPDE≔Λ⁡−ln⁡x+t,ux,ux,utx
normal⁡G⁡jetPDE
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)
Invariants⁡G,u⁡x,t
ux,−ln⁡x+t,ux,utx
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 ∂∂xu⁡x,t−1
dynamical_S≔_ξx=1,_ξt=0,_ηu=1ux
As mentioned, the differential invariants, for instance, of order 1, automatically come depending on higher derivatives
Invariants⁡dynamical_S,u⁡x,t,order=1
t,ux3+3⁢x⁢ux,x,2⁢u⁢ux,x+ux22⁢ux,x,ut⁢ux,x−ux⁢ux,tux,x
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
InvariantEquation⁡dynamical_S,u⁡x,t,order=2,name=Λ=0
* Partial match of 'name' against keyword 'arbitraryfunctionname'
Λ⁡t,∂∂xu⁡x,t,u⁡x,t⁢∂∂xu⁡x,t2−2⁢x⁢∂∂xu⁡x,t+u⁡x,t∂∂xu⁡x,t2+1,∂∂tu⁡x,t,∂∂xu⁡x,t3+2⁢∂2∂x2u⁡x,t⁢x+∂∂xu⁡x,t∂∂xu⁡x,t⁢∂2∂x2u⁡x,t⁢∂∂xu⁡x,t2+1,∂∂xu⁡x,t2+1⁢∂∂xu⁡x,t⁢∂2∂t∂xu⁡x,t∂2∂x2u⁡x,t,∂2∂x2u⁡x,t⁢∂2∂t2u⁡x,t−∂2∂t∂xu⁡x,t2∂2∂x2u⁡x,t=0
Note the φj depend on first and second derivatives of u⁡x,t, as it should be. In explicit form (use declare to have a more compact and readable output with derivatives displayed in jet notation)
declare⁡u⁡x,t
u⁡x,t⁢will now be displayed as⁢u
InvariantEquation⁡dynamical_S,u⁡x,t,order=2,explicit
ux,x=ux⁢ux2+1ux3⁢f__1⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut+ux⁢f__1⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut−2⁢x,ut,x=f__2⁡t,ux,ut,u⁢ux2−2⁢x⁢ux+uux2+1ux3⁢f__1⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut+ux⁢f__1⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut−2⁢x,ut,t=−f__2⁡t,ux,ut,u⁢ux2−2⁢x⁢ux+uux2+12−ux2⁢ux2+12⁢f__1⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut+2⁢ux⁢x⁢ux2+1+f__3⁡t,ux,u⁢ux2−2⁢x⁢ux+uux2+1,ut
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.
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
SymmetryTest
SymmetryTransformation
ToJet
Download Help Document