computes generalized integrating factors for a system of differential equations (DE) - Maple Help

Online Help

All Products    Maple    MapleSim


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

PDEtools[IntegratingFactors] - computes generalized integrating factors for a system of differential equations (DE)

PDEtools[IntegratingFactorTest] - tests whether a given list of expressions is a list of generalized integrating factors for the given system of differential equations

Calling Sequence

IntegratingFactors(PDESYS, DepVars, _mu = ..., displayfunctionality = ..., jetnotation = ..., simplifier = ..., split = ...)

IntegratingFactorTest(Mu, PDESYS, DepVars)

Parameters

PDESYS

-

a system consisting of an equation or a list of equations involving partial and/or ordinary (possibly not) differential equations

Mu

-

a generalized integrating factor returned by IntegratingFactors

DepVars

-

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

_mu = ...

-

optional - indicates the functional form of the generalized integrating factors

displayfunctionality = ...

-

optional - can be true (default) or false, to display the functionality on the left-hand-side of the mu[k] (generalized integrating factors) functions

jetnotation = ...

-

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

order

-

optional - indicates the maximum differential order of the derivatives entering the dependence of the integrating factors

simplifier = ...

-

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

split = ...

-

optional - can be true (default) or false, to split the DE system to be solved in order to compute the J[k] functions

typeofintegratingfactor = ...

-

optional - can be polynomial or functionfield.

Description

• 

Given a system Δ=0 consisting of N equations pden, n=1..N, where the independent variables are x1,x2,...=X, and the dependent variables are u1,u2,...=U, with dU denoting the set of partial derivatives of U, the generalized integrating factors are expressions μα,nX,U,dU such that Σμα,npden=DivergenceJα = 0, so Jα is a conserved current. These generalized integrating factors, also called characteristic functions of conserved currents (see reference [1]), coincide with the traditional integrating factors when there is only one independent variable, so that Δ is a system of ODEs.

• 

The command IntegratingFactors computes these generalized integrating factors. The command IntegratingFactorTest verifies the result for correctness.

• 

Given the system Δ=0 the output of IntegratingFactors is as a sequence of α lists, each one containing N μα,n, where n=1..N, satisfying Σμα,npden=DivergenceJα = 0.

• 

The μα,nX,U,dU are computed constructing the PDE system they satisfy by applying Euler's operator to Σμα,npden=0 , then solving this system for the μα,n using pdsolve.

• 

By default, the integrating factors are searched as functions depending on the derivatives of the unknowns of the system (specified as DepVars or automatically detected) up to the order d-1, where d is the highest order of derivatives entering PDESYS. This default can be changed by optionally passing the argument order = m, where m is a nonnegative integer.

• 

By default, the conserved currents are searched as functions with no pre-specified form, just with the dependency explained in the previous paragraph. This default can be changed with the option typeofconservedcurrent = ... where the right-hand-side can be polynomial or functionfield, respectively indicating a conserved current of polynomial type, or of a functionfield type with the meaning explained in FunctionField.

• 

By default, the functionality of μα,nX,U,dU, entering the left-hand-sides of each element in the returned lists, is displayed, the output is presented in functional notation instead of jet notation and is simplified with respect to its size. The PDE system solved to compute the μα,nX,U,dU is also split, when that is possible, before being tackled. All these defaults can be changed by passing the optional arguments displayfunctionality = ..., jetnotation = ..., simplifier = ..., split = false.

• 

It is also possible to directly specify the functionality expected for the μα,n using _mu = .... See the examples for a demonstration of the use of this parameter.

• 

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:

U:=diff_tableux,t:

declareU[]

ux,twill now be displayed asu

(1)

pde1:=Ut,t+Ux,x+UxU[]=0

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

(2)

Two generalized integrating factors are

μα:=IntegratingFactorspde1

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

(3)

Note that pde1 is already the divergence of a function, so that a constant (the number 1 in the result above) is an integrating factor. To verify for correctness these integrating factors use

mapIntegratingFactorTest,μα,pde1

0,0

(4)

The conserved currents are related to the generalized integrating factors via Sigma mu[alpha, n] pde[n] = Divergence J[alpha] = 0. These are the J[alpha] corresponding to the mu[alpha] computed above; they depend on arbitrary functions

Jα:=ConservedCurrentspde1

Jα:=_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

(5)

To verify these results use

mapConservedCurrentTest,Jα,pde1

0,0

(6)

An example where the integrating factor depends on an arbitrary function

pde2:=Ux,t+Ux,x,x+UxU[]=0

pde2:=uux+ut,x+ux,x,x=0

(7)

IntegratingFactorspde2

_μ1x,t,u,ux,ut,ux,x,ut,x,ut,t=_F1t,12u2+ut+ux,x

(8)

For this example, integrating factors up to order 1, that is, depending at most on first order derivatives, are

IntegratingFactorspde2,order=1

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

(9)

which is in agreement with the general result obtained first. This is a related conserved current of order 1

ConservedCurrentspde2,order=1

_Jxx,t,u,ux,ut=_F1x,t,uut+∫u_F1tⅆ_a+_F3x,t,_Jtx,t,u,ux,ut=_F1x,t,uux+∫u_F1xⅆ_a+∫_F3xⅆt+_F4x

(10)

Specifying directly the functionality expected also confirms that there is no non-trivial integrating factor depending only on ux and ut but there is one depending on an arbitrary function of u, ut and ux,x

IntegratingFactorspde2,_μ=fUx,Ut

_μ1ux,ut=1

(11)

IntegratingFactorspde2,_μ=fUx,x,Ut,U[]

_μ1ux,x,ut,u=_F1u2+2ut+2ux,x

(12)

In various cases it is simpler, or of more use, to compute integrating factors of polynomial type, or with a mathematical function dependency on the field of functions of the input system. For these purposes use the option typeofintegratingfactor = ... where the right-hand-side can be polynomial or functionfield. For example, for pde2, a polynomial integrating factor, presented without specializing the arbitrary constants (option split = false) is

IntegratingFactorspde2,type=polynomial,split=false

* Partial match of 'type' against keyword 'typeofintegratingfactor'

_μ1x,t,u,ux,ut,ux,x,ut,x,ut,t=u2+2ut+2ux,x_C3+_C4t2+_C2t+_C1

(13)

IntegratingFactorTest,pde2

0

(14)

The following application of Euler's operator to pde2 shows that pde2 is already a divergence of a function

Eulerpde2

0

(15)

This is a conserved current with the same functionality of the last integrating factor computed and a verification of the result

Jα:=ConservedCurrentspde2,_J=fUx,x,Ut,U[]

Jα:=_Jxux,x,ut,u=_F1u2+2ut+2ux,x,_Jtux,x,ut,u=1

(16)

ConservedCurrentTestJα,pde2

0

(17)

See Also

ConservedCurrents, ConservedCurrentTest, Euler, PDEtools

References

  

[1] Olver, P.J. Applications of Lie Groups to Differential Equations. Graduate Texts in Mathematics. Springer-Verlag, 1993.


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