Difference between revisions of "Periodic orbits"

From openpipeflow.org
Jump to navigation Jump to search
Line 14: Line 14:
  
 
=== Periodic orbits in pipe flow ===
 
=== Periodic orbits in pipe flow ===
Here, by the 'state' of the system we mean the speed and direction of the flow of fluid at every position in the flow.  By listing the velocity at every point we form a state-vector, '''a'''=(a1,a2,a3,a4,...).  As there are infinitely many positions, the dimension (length) of '''a''' is formally infinite.  However, for dissipative systems, the dynamics typically lives on a space of much smaller dimension (like a 2D-sheet within 3D space), e.g. only a few dozen dimensions.  (Only realistic states occur as time progresses.)  To visualise this space, it is necessary to project onto a shorter, more meaningful vector, e.g. (x,y)=(energy,dissipation).
+
Here, by the 'state' of the system we mean the speed and direction of the flow of fluid at every position in the flow.  By listing the velocity at every point we form a state-vector, '''a'''=(a1,a2,a3,a4,...).  As there are infinitely many positions, the dimension (length) of '''a''' is formally infinite.  However, for dissipative systems, the dynamics typically lives on a space of much smaller dimension (like a 2D-sheet within 3D space; only realistic states, those on the sheet, occur as time progresses.)  To visualise this space, it is necessary to project onto a shorter, more meaningful vector, e.g. (x,y)=(energy,dissipation).
  
 
The following video shows a selection of periodic orbits in pipe flow, visualised in terms input energy (I), dissipation (D) and kinetic energy (E) (relative to their laminar (non-turbulent) values):
 
The following video shows a selection of periodic orbits in pipe flow, visualised in terms input energy (I), dissipation (D) and kinetic energy (E) (relative to their laminar (non-turbulent) values):

Revision as of 08:57, 27 June 2016

Shadowing of two periodic orbits (red and green) leads to a longer one (blue).

Apparently complex dynamics can originate from a small set of dynamically important recurring patterns. The image to the right shows two recurring patterns (green,red) that exist for the double pendulum. (A simple android application and screenshots of the double-pendulum can be found HERE;) As time progresses, the position of the lower pendulum traces out the lines shown, (x,y)=(position of lower pendulum). As the line is a closed loop, it traces out a repeating pattern - a 'periodic orbit'.

Observe that the longer orbit (blue) can be considered to be constructed from 'shadowing' of the two simpler orbits. To understand a lot about the dynamics it is sufficient to find the short periodic orbits. These we might call our 'alphabet', out of which infinitely many 'words', i.e. longer orbits, may be constructed.

An excellent (graduate level) free text on periodic orbits can be found at chaosbook.org.


Periodic orbits in pipe flow

Here, by the 'state' of the system we mean the speed and direction of the flow of fluid at every position in the flow. By listing the velocity at every point we form a state-vector, a=(a1,a2,a3,a4,...). As there are infinitely many positions, the dimension (length) of a is formally infinite. However, for dissipative systems, the dynamics typically lives on a space of much smaller dimension (like a 2D-sheet within 3D space; only realistic states, those on the sheet, occur as time progresses.) To visualise this space, it is necessary to project onto a shorter, more meaningful vector, e.g. (x,y)=(energy,dissipation).

The following video shows a selection of periodic orbits in pipe flow, visualised in terms input energy (I), dissipation (D) and kinetic energy (E) (relative to their laminar (non-turbulent) values):

The setup for the previous video was pipe flow with flow rate measured by Re=DU/nu=2500 (D=diameter, U=mean axial speed, nu=kinematic viscosity). The next video shows flow in the lab frame. Regions of flow slower than the mean flow are indicated in blue, called 'streaks'. Rotating 'vortex structures' are indicated in yellow. To simplify analysis a 4-fold rotational symmetry has been imposed. This does not significantly affect statistical properties of the flow, e.g. the turbulent friction factor.


When viewed in a moving frame it is possible to see recurring patterns, i.e. periodic orbits (next video). The speed at which the frame moves, however, is different for every periodic orbit. To determine this automatically we use the Method of slices.

Periodic orbits exist that encompass the broad range of observed flow patterns. The following video shows an example of a periodic orbit that is representative a spatially localised patch of turbulence: