PDEtools
SymmetrySolutions
computes the group invariant (point symmetry) solutions for a given PDE system
Calling Sequence
Parameters
Description
Examples
SymmetrySolutions(Solution, S, 'options'='value')
Solution
-
an equation, a set (list) of equations, or a set (list) of sets (lists) of equations, representing the solution to some system of equations (itself not required)
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
simplifier = ...
optional - indicates the simplifier to be used instead of the default simplify/size
Given the Solution to a problem, as an equation or a set (list) of equations, or given a set (list) of these solutions, and given the infinitesimal generators of a symmetry S, or a list of them, SymmetrySolutions computes other solutions obtained by exploring the given symmetries. This command complements Infinitesimals and InvariantSolutions, which respectively compute the symmetry infinitesimals and the group (symmetry) invariant solutions of a PDE system.
The dependent variables of the problem do not need to be indicated: SymmetrySolutions will automatically take the left-hand-sides of the given solutions as dependent variables.
What SymmetrySolutions does can be decomposed into two steps that can be reproduced interactively:
compute the finite form of the transformation associated to all the given symmetries (can be performed using SymmetryTransformation);
apply this transformation to each of the given solutions (can be performed using dchange).
You can optionally specify a simplifier to be applied to the result using the optional argument simplifier = ...
To avoid having to remember the optional keywords, if you type the keyword misspelled, or just a portion of it, a matching against the correct keywords is performed, and when there is only one match, the input is automatically corrected.
Consider the wave equation in four dimensions to avoid redundant typing on input and on the display use diff_table and declare
with⁡PDEtools,SymmetrySolutions,Infinitesimals:
Consider for instance the wave equation in two dimensions
pde1≔diff⁡u⁡x,y,t,x,x+diff⁡u⁡x,y,t,y,y=0
pde1≔∂2∂x2u⁡x,y,t+∂2∂y2u⁡x,y,t=0
A solution for pde[1] obtained separating the variables by product
sol1≔pdsolve⁡pde1,HINT=`*`,build
sol1≔u⁡x,y,t=f__4⁡t⁢c__1⁢ⅇ_c1⁢x⁢c__3⁢sin⁡_c1⁢y+f__4⁡t⁢c__1⁢ⅇ_c1⁢x⁢c__4⁢cos⁡_c1⁢y+f__4⁡t⁢c__2⁢c__3⁢sin⁡_c1⁢yⅇ_c1⁢x+f__4⁡t⁢c__2⁢c__4⁢cos⁡_c1⁢yⅇ_c1⁢x
Some Infinitesimals of polynomial type for pde[1] are
S1≔Infinitesimals⁡pde1,typeofsymmetry=polynomial,displayfunctionality=false−2..−1
S1≔_ξx=x22−y22,_ξy=x⁢y,_ξt=0,_ηu=0,_ξx=x⁢y,_ξy=−x22+y22,_ξt=0,_ηu=0
New solutions for pde[1] exploring these infinitesimals can be computed using SymmetrySolutions as follows
newsol1≔SymmetrySolutions⁡sol1,S1
newsol1≔u⁡x,y,t=f__4⁡t⁢ⅇ4⁢x2+y2⁢_ε+2⁢x⁢_c14+x2+y2⁢_ε2+4⁢_ε⁢x⁢c__1+c__2⁢cos⁡4⁢_c1⁢y4+x2+y2⁢_ε2+4⁢_ε⁢x⁢c__4+sin⁡4⁢_c1⁢y4+x2+y2⁢_ε2+4⁢_ε⁢x⁢c__3⁢ⅇ−2⁢x2+y2⁢_ε+2⁢x⁢_c14+x2+y2⁢_ε2+4⁢_ε⁢x,u⁡x,y,t=f__4⁡t⁢cos⁡2⁢_c1⁢x2+y2⁢_ε+2⁢y4+x2+y2⁢_ε2+4⁢y⁢_ε⁢c__4+sin⁡2⁢_c1⁢x2+y2⁢_ε+2⁢y4+x2+y2⁢_ε2+4⁢y⁢_ε⁢c__3⁢ⅇ8⁢_c1⁢x4+x2+y2⁢_ε2+4⁢y⁢_ε⁢c__1+c__2⁢ⅇ−4⁢_c1⁢x4+x2+y2⁢_ε2+4⁢y⁢_ε
Note the introduction of the Lie group parameter _epsilon; the new solutions are valid for any value (not depending on x, y or t) of it. These solutions can be verified using pdetest
map⁡pdetest,newsol1,pde1
0,0
See Also
dchange
Infinitesimals
InvariantSolutions
InvariantTransformation
pdetest
pdsolve
PolynomialSolutions
SimilaritySolutions
SimilarityTransformation
SymmetryTransformation
Download Help Document