DifferentialGeometry
Flow
calculate the one parameter group of differeomorphisms (flow) of a vector field
Calling Sequence
Parameters
Description
Examples
Flow(X, var, options)
X
-
a vector field
var
an unassigned Maple name, the flow parameter
options
optional arguments to pass to the Maple command dsolve for solving the ordinary differential equations for the flow
The flow of a vector field X on a manifold M is a one parameter group of transformations Phi_t: M -> M such that for all p in M, diff(Phi_t(p), t) = X(Phi_t(p)) and Phi_0(p) = p. For each fixed t, Phi_t is a local diffeomorphism of M and Phi_t o Phi_s = Phi_(t + s).
The flow of X is calculated by solving a first order system of ordinary differential equations with the Maple dsolve command.
If dsolve fails to solve these odes, the Flow command returns NULL.
The command Flow returns a transformation whose domain and range coincide with the manifold on which X is defined.
With the option ode = true, the system of odes (with initial conditions) defining the flow is returned.
With the option initialpoint = [x1 = a, x2 = b, ...], the flow though the specific point [a, b, ...] is calculated.
With the option dsolvehints = [hints], the list of optional arguments hints is passed to dsolve.
A customized ode solver can be used in place of dsolve though the use of the Preference command.
This command is part of the DifferentialGeometry package, and so can be used in the form Flow(...) only after executing the command with(DifferentialGeometry). It can always be used in the long form DifferentialGeometry:-Flow.
with⁡DifferentialGeometry:
DGsetup⁡x,y,z,M:
Example 1.
Calculate the flow Phi_t for the vector field X.
X≔evalDG⁡−y⁢D_x+x⁢D_y+14⁢z⁢D_z
X≔−y⁢D_x+x⁢D_y+z⁢D_z4
Phi_t≔Flow⁡X,t
Phi_t≔x=−y⁢sin⁡t+x⁢cos⁡t,y=y⁢cos⁡t+x⁢sin⁡t,z=z⁢ⅇt4
Flow plots
Plot the flows for various initial conditions:
C≔map2⁡ApplyTransformation,Phi_t,0,1,12,1,0,0,1,0,12
C≔−sin⁡t,cos⁡t,ⅇt42,cos⁡t,sin⁡t,0,cos⁡t,sin⁡t,ⅇt42
plotsspacecurve⁡C,t=−π..π,axes=normal
We check that Phi_s o Phi_t is Phi_(t + s).
combine⁡ComposeTransformations⁡eval⁡Phi_t,t=s,Phi_t
x=−y⁢sin⁡t+s+x⁢cos⁡t+s,y=y⁢cos⁡t+s+x⁢sin⁡t+s,z=z⁢ⅇt4+s4
eval⁡Phi_t,t=t+s
We check that the derivative of the flow with respect to t coincides with the vector field evaluated along the flow:
C1≔ApplyTransformation⁡Phi_t,a,b,c
C1≔−b⁢sin⁡t+a⁢cos⁡t,b⁢cos⁡t+a⁢sin⁡t,c⁢ⅇt4
Differentiate the components of this curve with respect to t.
Y1≔DGzip⁡diff⁡C1,t,D_x,D_y,D_z,plus
Y1≔−b⁢cos⁡t+a⁢sin⁡t⁢D_x+−b⁢sin⁡t+a⁢cos⁡t⁢D_y+c⁢ⅇt4⁢D_z4
C2≔ApplyTransformation⁡Phi_t,x=a,y=b,z=c
C2≔x=−b⁢sin⁡t+a⁢cos⁡t,y=b⁢cos⁡t+a⁢sin⁡t,z=c⁢ⅇt4
Y2≔eval⁡X,C2
Y2≔−b⁢cos⁡t+a⁢sin⁡t⁢D_x+−b⁢sin⁡t+a⁢cos⁡t⁢D_y+c⁢ⅇt4⁢D_z4
Y1&minusY2
0⁢D_x
Example 2.
We find the flow of the vector X through the point (1, 0, 0).
Flow⁡X,t,initialpoint=x=1,y=0,z=0
x=cos⁡t,y=sin⁡t,z=0
Example 3.
We obtain the ode defining the flow for X. The result consists of a sequence of 3 sets: the ode, the initial conditions, and the dependent variables.
Flow⁡X,t,ode=true
−_z1⁡t+ⅆⅆt_z2⁡t,_z2⁡t+ⅆⅆt_z1⁡t,−_z3⁡t4+ⅆⅆt_z3⁡t,_z1⁡0=x,_z2⁡0=y,_z3⁡0=z,_z1⁡t,_z2⁡t,_z3⁡t
Example 4.
Flow⁡X,t,dsolvehints=method=laplace
x=−y⁢sin⁡t+x⁢cos⁡t,y=y⁢cos⁡t+x⁢sin⁡t,z=z⁢ⅇt4
See Also
ApplyTransformation
ComposeTransformations
InfinitesimalTransformation
Preferences
Transformation
Download Help Document