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

Online Help

All Products    Maple    MapleSim


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

Calling Sequence

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)

Parameters

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=Mapx,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).

Description

• 

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.

Examples

withLieAlgebrasOfVectorFields:

 

DE[1] below is obviously linearizable by the transformation uh=u2,xh=x:

DE1diffux2,x,x,x+ux2=0

DE16ⅆⅆxuxⅆ2ⅆx2ux+2uxⅆ3ⅆx3ux+ux2=0

(1)

MDE1MapDEDE1,x,u,ξ,η,ToLinearDE,ψ,φ

MDE1tableRedMapDE=2u2φx,u=uφx,uu,uψx,u=0,MapDEPivs=uφx,u0,xψx,u0,TargetDE=ⅆ3ⅆxh3uhxh=uhxhwhereuhxh0,IntegratedMap=xh=x,uh=u2,MapExists=true

(2)

 

DE2uxdiffdiffdiffux,x,x,x+3diffux,xdiffdiffux,x,x3diffux,x2+uxdiffdiffux,x,x+12uxdiffux,x+x8xuxdiffux,x+x41+ux2+x2ux2+x2=0

DE2uxⅆ3ⅆx3ux+3ⅆⅆxuxⅆ2ⅆx2ux3ⅆⅆxux2+uxⅆ2ⅆx2ux+12uxⅆⅆxux+x8xuxⅆⅆxux+x41+ux2+x2ux2+x2=0

(3)

MDE2MapDEDE2,x,u,ξ,η,ToLinearDE,xh,uh,ψ,φ

MDE2tableRedMapDE=2x2φx,u=22uxφx,uu2x+2u2φx,uux2uφx,uu2uφx,ux2u3,xψx,u=uψx,uxu,MapDEPivs=uψx,u0,xφx,uuuφx,ux0,TargetDE=ⅆ3ⅆxh3uhxh=xh+1uhxhxhwhereuhxhⅆⅆxhuhxh0,uhxh2+xh0,uhxh0,IntegratedMap=xh=u2+x2,uh=x,MapExists=true

(4)

 

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:

DE3diffux,t,x,x=diffux,t,tux,tdiffux,t,x

DE32x2ux,t=tux,tux,txux,t

(5)

MDE3MapDEDE3,x,t,u,ξ,τ,η,ToLinearDE,x,t,u,ψ,ϒ,φ

MDE3tableMapExists=false

(6)

 

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

DE4diffvx,t,x=ux,t,diffvx,t,t=diffux,t,x+ux,t22

DE4xvx,t=ux,t,tvx,t=xux,t+ux,t22

(7)

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)).

MDE4MapDEDE4,x,t,u,v,ξ,τ,η,β,ToLinearDE,xh,th,uh,vh,ψ,ϒ,ϕ,φ

MDE4tableRedMapDE=2u2φx,t,u,v=0,2u2ϕx,t,u,v=0,2uvφx,t,u,v=uφx,t,u,v2,2uvϕx,t,u,v=uϕx,t,u,v2,2v2φx,t,u,v=vφx,t,u,v2,2v2ϕ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,vuφx,t,u,vuϕx,t,u,vvφx,t,u,v0,xψx,t,u,vtϒx,t,u,vtψx,t,u,vxϒx,t,u,v0,xψx,t,u,v0,TargetDE=xhvhxh,th=uhxh,th2,thvhxh,th=xhuhxh,th2,2xh2uhxh,th=thuhxh,thwherevhxh,th0,IntegratedMap=xh=x,th=t,uh=uⅇv2,vh=ⅇv2,MapExists=true

(8)

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).

References

  

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.

Compatibility

• 

The LieAlgebrasOfVectorFields[MapDE] command was introduced in Maple 2021.

• 

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

See Also

PDEtools