LieAlgebrasOfVectorFields
MapDE
algorithmically checks if an input system of polynomially nonlinear ODE or PDE can be mapped into a target class of linear DEs by an invertible Map of the variables of the input system
Calling Sequence
Parameters
Description
Examples
References
Compatibility
MapDE(InputDE, TargetClass, Map, Options)
MapDE([DE,[x,u],[xi,eta]], ToLinearDE, [psi,phi], Options)
MapDE([DE,[x,u],[xi,eta]], [ToLinearDE, [xh,uh]], [psi,phi], Options)
InputDE
-
A list of lists [DE, vars, infinitesimals] for an input (source) DE system. DE is a list of polynomially nonlinear partial or ordinary differential equations and inequations, and vars is a list x,u where x=x1,x2,…,xn is a list of independent variable names and u=u1,u2,…,um is a list of dependent variable names for the DE. infinitesimals is a list ξ,η where ξ=ξ1,ξ2,…,ξn and η=η1,…,ηm are lists of the components of vectorfields of the Lie algebra corresponding to continuous Lie point symmetry transformations of the DE in the space of independent and dependent variables of the DE.
TargetClass
Properties of the class of target DE. Currently TargetClass is either ToLinearDE or has form [ToLinearDE, xh,uh] where xh=xh1,…,xhn is a list of independent variable names and uh=uh1,…,uhm is a list of dependent variable names for the target DE.
Map
ψ,φ where ψ and φ are lists of mapping functions corresponding to the independent variables xh and dependent variables uh. Thus Map transforms the InputDE to a member of the TargetClass where xh,uh=Map⁡x,u=ψ⁡x,u,φ⁡x,u at points x,u where the Jacobian of Map is nonzero.
Options
OutputDetails when included in Options returns additional details of the computation (including structure and dimensions of various objects computed in intermediate computations).
The command MapDE determines the invertible linearizability of nonlinear differential polynomial systems and attempts to construct the map if it exists. Thus it determines the existence of an invertible mapping xh,uh=ψ⁡x,u,φ⁡x,u on points x,u in the space of independent and dependent variables of the InputDE to points xh,uh in the space of independent and dependent variables of a target linear DE.
The output of MapDE is a table. The output table includes the label MapExists = true or MapExists = false, if a mapping of the InputDE to a linear DE exists or does not exist. Additional outputs include the reduced Map DE (RedMapDE) which are differential equations satisfied by the functions ψ,φ of the mapping, together with inequality (pivot) conditions on ψ,φ given by the output table entry labeled by MapDEPivs. If pdsolve is successful in explicitly integrating the RedMapDE system, then the output also includes the integrated form of the transformations labeled by IntegratedMap, and an explicit form of the target DE, labeled by TargetDE.
A user can get more information about the mapping system, including the structure and dimensions of various objects (solution spaces, derived algebra, Lie algebra structure, etc) by including OutputDetails or by assigning infolevel[MapDE] to a non-negative integer value.
with⁡LieAlgebrasOfVectorFields:
DE[1] below is obviously linearizable by the transformation uh=u2,xh=x:
DE1≔diff⁡u⁡x2,x,x,x+u⁡x2=0
DE1≔6⁢ⅆⅆxu⁡x⁢ⅆ2ⅆx2u⁡x+2⁢u⁡x⁢ⅆ3ⅆx3u⁡x+u⁡x2=0
MDE1≔MapDE⁡DE1,x,u,ξ,η,ToLinearDE,ψ,φ
MDE1≔table⁡RedMapDE=∂2∂u2φ⁡x,u=∂∂uφ⁡x,uu,∂∂uψ⁡x,u=0,MapDEPivs=∂∂uφ⁡x,u≠0,∂∂xψ⁡x,u≠0,TargetDE=ⅆ3ⅆxh3uh⁡xh=−uh⁡xhwhereuh⁡xh≠0,IntegratedMap=xh=x,uh=u2,MapExists=true
DE2≔u⁡x⁢diff⁡diff⁡diff⁡u⁡x,x,x,x+3⁢diff⁡u⁡x,x⁢diff⁡diff⁡u⁡x,x,x−3⁢diff⁡u⁡x,x2+u⁡x⁢diff⁡diff⁡u⁡x,x,x+12u⁡x⁢diff⁡u⁡x,x+x−8⁢x⁢u⁡x⁢diff⁡u⁡x,x+x4⁢1+u⁡x2+x2u⁡x2+x2=0
DE2≔u⁡x⁢ⅆ3ⅆx3u⁡x+3⁢ⅆⅆxu⁡x⁢ⅆ2ⅆx2u⁡x−3⁢ⅆⅆxu⁡x2+u⁡x⁢ⅆ2ⅆx2u⁡x+12u⁡x⁢ⅆⅆxu⁡x+x−8⁢x⁢u⁡x⁢ⅆⅆxu⁡x+x4⁢1+u⁡x2+x2u⁡x2+x2=0
MDE2≔MapDE⁡DE2,x,u,ξ,η,ToLinearDE,xh,uh,ψ,φ
MDE2≔table⁡RedMapDE=∂2∂x2φ⁡x,u=−−2⁢∂2∂u∂xφ⁡x,u⁢u2⁢x+∂2∂u2φ⁡x,u⁢u⁢x2−∂∂uφ⁡x,u⁢u2−∂∂uφ⁡x,u⁢x2u3,∂∂xψ⁡x,u=∂∂uψ⁡x,u⁢xu,MapDEPivs=∂∂uψ⁡x,u≠0,∂∂xφ⁡x,u⁢u−∂∂uφ⁡x,u⁢x≠0,TargetDE=ⅆ3ⅆxh3uh⁡xh=−xh+1⁢uh⁡xhxhwhereuh⁡xh⁢ⅆⅆxhuh⁡xh≠0,−uh⁡xh2+xh≠0,uh⁡xh≠0,IntegratedMap=xh=u2+x2,uh=−x,MapExists=true
For the example above, the transformation xh=u2+x2,uh=−x maps DE[2] to the given TargetLinearDE, which can be independently checked by applying PDEtools[dchange].
Burgers' PDE below is the simplest nonlinear model combining convection and diffusion:
DE3≔diff⁡u⁡x,t,x,x=diff⁡u⁡x,t,t−u⁡x,t⁢diff⁡u⁡x,t,x
DE3≔∂2∂x2u⁡x,t=∂∂tu⁡x,t−u⁡x,t⁢∂∂xu⁡x,t
MDE3≔MapDE⁡DE3,x,t,u,ξ,τ,η,ToLinearDE,x,t,u,ψ,ϒ,φ
MDE3≔table⁡MapExists=false
MapDE proves above that there is no invertible change of variables in x,t,u-space linearizing Burgers' equation. Bluman et al (2010) nonlocally embed Burgers' PDE in the system below
DE4≔diff⁡v⁡x,t,x=u⁡x,t,diff⁡v⁡x,t,t=diff⁡u⁡x,t,x+u⁡x,t22
DE4≔∂∂xv⁡x,t=u⁡x,t,∂∂tv⁡x,t=∂∂xu⁡x,t+u⁡x,t22
The reader can check that Burgers' PDE results from computing cross partials. Application of MapDE below reveals that Burgers' DE[4] is linearizable in (x,t,u,v)-space and recovers the famous Hopf-Cole linearizing transformation (also see Bluman et al (2010)).
MDE4≔MapDE⁡DE4,x,t,u,v,ξ,τ,η,β,ToLinearDE,xh,th,uh,vh,ψ,ϒ,ϕ,φ
MDE4≔table⁡RedMapDE=∂2∂u2φ⁡x,t,u,v=0,∂2∂u2ϕ⁡x,t,u,v=0,∂2∂u∂vφ⁡x,t,u,v=∂∂uφ⁡x,t,u,v2,∂2∂u∂vϕ⁡x,t,u,v=∂∂uϕ⁡x,t,u,v2,∂2∂v2φ⁡x,t,u,v=∂∂vφ⁡x,t,u,v2,∂2∂v2ϕ⁡x,t,u,v=∂∂vϕ⁡x,t,u,v2,∂∂uϒ⁡x,t,u,v=0,∂∂uψ⁡x,t,u,v=0,∂∂vϒ⁡x,t,u,v=0,∂∂vψ⁡x,t,u,v=0,MapDEPivs=∂∂vϕ⁡x,t,u,v⁢∂∂uφ⁡x,t,u,v−∂∂uϕ⁡x,t,u,v⁢∂∂vφ⁡x,t,u,v≠0,∂∂xψ⁡x,t,u,v⁢∂∂tϒ⁡x,t,u,v−∂∂tψ⁡x,t,u,v⁢∂∂xϒ⁡x,t,u,v≠0,∂∂xψ⁡x,t,u,v≠0,TargetDE=∂∂xhvh⁡xh,th=uh⁡xh,th2,∂∂thvh⁡xh,th=∂∂xhuh⁡xh,th2,∂2∂xh2uh⁡xh,th=∂∂thuh⁡xh,thwherevh⁡xh,th≠0,IntegratedMap=xh=x,th=t,uh=u⁢ⅇv2,vh=ⅇv2,MapExists=true
The above apparently paradoxical result is resolved, by noting that the two spaces x,t,u-space and x,t,u,v-space are noninvertibly related, since non-unique v yields u by integration. Systematically exploring such nonlocal embeddings is discussed in Bluman et al (2010).
Z. Mohammadi, G. Reid, and S.-L.T Huang. "Introduction of the MapDE algorithm for determination of mappings relating differential equations." Proceedings of ISSAC 2019, pp.331-338. ACM Press, 2019.
Z. Mohammadi, G. Reid, and S.-L.T. Huang. "Symmetry-based algorithms for invertible mappings of polynomially nonlinear PDE to linear PDE". Mathematics in Computer Science, 1–24, 2020.
G. W. Bluman, A. F. Cheviakov, and S. C. Anco. Applications of Symmetry Methods to Partial Differential Equations, Springer, Chapter 2, 2010.
D. Lyakhov, V. Gerdt, and D. Michels. "Algorithmic verification of linearizability for ordinary differential equations." Proceedings of ISSAC 2017, pp. 285–292. ACM Press, 2017.
The LieAlgebrasOfVectorFields[MapDE] command was introduced in Maple 2021.
For more information on Maple 2021 changes, see Updates in Maple 2021.
See Also
PDEtools
Download Help Document