PDEtools - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : PDEtools : Symmetries : PDEtools/Euler

PDEtools

  

Euler

  

Euler's operator - return the conditions that a DE system must satisfy in order to be the divergence of a current

 

Calling Sequence

Parameters

Description

Examples

Calling Sequence

Euler(PDESYS, DepVars, integrabilityconditions = ..., jetnotation = ..., simplifier = ...)

Parameters

PDESYS

-

a system consisting of a single equation or a set or list of them, involving partial or ordinary (possibly not) differential equations

DepVars

-

(optional) a specification of the unknown(s) in PDESYS

jetnotation = ...

-

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

integrabilityconditions = ...

-

(optional) - can be true (default) or false; indicates whether or not to reduce the determining PDE output system taking into account its integrability conditions

simplifier = ...

-

optional - indicates the simplifier to be used instead of the default simplify/size

Description

• 

The Euler command represents Euler's operator: it returns zero when applied to an expression that is the divergence of another expression also called "current". When this divergence is also equal to zero this current is called conserved current.

• 

Denoting the unknowns of a DE system by U, their derivatives by dU and the independent variables by X, Euler's operator is constructed as (d/dU + (-1)^n d/dX d/ddU), where n is the differential order of dU, d/dX is the differentiation operator such that dU = d/dx U and in Euler's operator there is one term of the form (-1)^n d/dX d/ddU for each dU in the DE system.

• 

If DepVars is not given, Euler will consider all the differentiated unknown functions in PDESYS as unknowns of the problems. Specifying DepVars however permits not only restricting the unknowns in different ways but also specifying unknowns of the problems which do not appear differentiated in PDESYS.

• 

When applying Euler's operator to a DE system constituted by N equations in M unknowns, there will be N x M expressions; that is: each of the M unknowns generates an operator of the form just described and each of these operators generates one expression when applied to each of the N equations of the system.

• 

This set of expressions resulting from applying Euler's operator to a DE system, representing the conditions for the DE system to be the divergence of a current, is by default simplified with respect to its integrability conditions, then each expression in the set is simplified with respect to its size. This default behavior can be changed by providing one or more of the arguments integrabilityconditions = false, simplifier = none, or simplifier = any_desired_simplification.

• 

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.

Examples

Consider the following PDE "system" consisting of a single pde.

withPDEtools:

Udiff_tableux,t:

declareU[]

ux,twill now be displayed asu

(1)

pde1Ut,t+Ux,x+UxU[]=0

pde1:=uux+ut,t+ux,x=0

(2)

This system automatically satisfies the conditions for being a divergence mentioned in the Description:

Eulerpde1

0

(3)

Hence pde1 is the divergence of a current

ConservedCurrentspde1

_Jxx,t,u,ux,ut=ux+_F3x,t,uut+∫u_F3t+_aⅆ_a+_F5x,t,_Jtx,t,u,ux,ut=∫u_F3xⅆ_a+∫_F5xⅆt_F3x,t,uux+_F6x+ut,_Jxx,t,u,ux,ut=tux+_F3x,t,uut+∫u_at+_F3tⅆ_a+_F5x,t,_Jtx,t,u,ux,ut=∫u_F3xⅆ_a+∫_F5xⅆt_F3x,t,uux+_F6x+tutu

(4)

and admits a constant integrating factor:

IntegratingFactorspde1

_μ1x,t,u,ux,ut=1,_μ1x,t,u,ux,ut=t

(5)

When combined with the rest of the Maple library, the Euler operator can serve varied purposes. Consider for example deriving the most general form of the divergence of a current that is also a first order PDE in two variables. The starting point is a generic expression, Δx,t,u,ux,ut, so it depends only on the first order derivatives.

pde2Δx,t,U[],Ux,Ut=0

pde2:=Δx,t,u,ux,ut=0

(6)

The conditions that Delta must satisfy in order to be a divergence are:

conditions_for_divergenceEulerpde2

conditions_for_divergence:=Δu,ux=Δu,utut+ΔuΔt,utΔx,uxux,Δux,ux=0,Δut,ux=0,Δut,ut=0

(7)

These conditions can be integrated.

Divergence_1PDE_general_formpdsolveconditions_for_divergence

Divergence_1PDE_general_form:=Δx,t,u,ux,ut=_F3x,t,uux+_F4x,t,uut+∫_F4t+_F3xⅆu+_F6x,t

(8)

Verify that these conditions are sufficient by applying Euler's operator to this result. First convert the result from jet notation to function notation.

answerFromJetrhsDivergence_1PDE_general_form1,ux,t=0

answer:=_F3x,t,uux+_F4x,t,uut+∫u_F4t+_F3xⅆ_a+_F6x,t=0

(9)

So the above is the most general form of a divergence that is also a first order PDE. The following verifies that this form is correct.

Euleranswer,ux,t

0

(10)

The most general form of a second order linear PDE in two independent variables that is also a divergence of a current can be derived in a similar way, starting with the following definition.

declareAx,t

Ax,twill now be displayed asA

(11)

pde3A[0]x,tU[]+A[1]x,tUx+A[2]x,tUt+A[3]x,tUx,x+A[4]x,tUx,t+A[5]x,tUt,t=0

pde3:=uA0+A1ux+A2ut+A3ux,x+A4ut,x+A5ut,t=0

(12)

The conditions for divergence, in the form of equations satisfied by the A[j]( x, t ) are obtained by applying Euler's operator.

Eulerpde3,U[]

A0A2,tA1,x+A5,t,t+A4,t,x+A3,x,x

(13)

Note that in the above calculation, the dependent variable of the problem must be specified, otherwise A[j]( x, t ) for all j would be also picked up as dependent variables. The result above is a single expression from which A[j]( x, t ) for one j can be isolated; the simplest form is achieved by isolating A[0].

conditionisolate1,A[0]x,t

condition:=A0=A2,t+A1,xA5,t,tA4,t,xA3,x,x

(14)

Divergence_linear_2PDEsubscondition,pde3

Divergence_linear_2PDE:=A2,t+A1,xA5,t,tA4,t,xA3,x,xu+A1ux+A2ut+A3ux,x+A4ut,x+A5ut,t=0

(15)

This result can be verified applying Euler's operator to it.

EulerDivergence_linear_2PDE,U[]

0

(16)

See Also

ConservedCurrents

ConservedCurrentTest

IntegratingFactors

IntegratingFactorTest

PDEtools

pdsolve

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam