This is an old revision of the document!
Let us see how this format works. For now I will just start writing. I hope that greater flexibility in the structure (multiple pages and links) will make it easy to reorganize as the material gathers.
I have been dissatisfied for quite a while with the projection basis sets I've used for plotting periodic orbits in the HKW cell. We don't have good equilibria there, so we can't easily do the translational-basis set we used for the narrow cell in the JFM08 paper. I have previously used a POD basis from the data of a single periodic orbit (P87p89). The plots were good, in that they showed turbulent trajectories lying in the same regions of state space as the periodic orbits, and forming similar pretty-well-defined shapes in whatever projection we looked at. But Predrag never liked this since it involved POD and averages, and to me, it was rhetorically weak to use a basis based on periodic orbits to argue that the periodic orbits live where turbulence is.
So, I decided to compare the orbit-POD basis to a translational basis based on a single velocity field (for Predrag) and a turbulent-data POD basis (for me). It turns out they're all remarkably similar.
Some key points:
The mean of the P87p89 orbit and its half-cell shift in x, , with and without the laminar flow.
Three basis sets
Main things to note:
The basis sets are very similar. That means that P87.89 explores the same region of state space as turbulence, in the crude let's-pretend-everything-is-a-Gaussian sense of the POD. Still, that is something, beyond the 1d plots showing that some orbits capture the mean and the Reynolds stresses within a few percent. Those 1d plots are averaged over x and z as well as time. Here we preserve the symmetry that breaks x,z homogeneity and find that the 3d average and multidimensional variances are very similar.
(Correction, 2009-01-05: The P87 POD basis elements shown above are corrected since the original posting. I think I had mistakenly copied the turbulent POD basis elements. The correspondence between the two sets is weaker but still strong. Note that the POD modes are ordered differently between the turbulent POD and the P87. E.g. turbulent POD mode e4 matches P87 POD mode e3. The sets are ordered by the variances of the imaginary Gaussian.)
It is not surprising that the translational basis based on the P87p89 mean is close to the P87p89 POD basis. The means of the P87p89 and P87p89 orbit are far from the origin and relatively close together. Therefor it is not surprising that the first two POD modes of the full orbit are close to of the means of the two orbits.
The P87p89 POD basis, however, continues on to e0, e1, e2, e3, … whereas the translational basis is just 2d. And the P87p89 POD basis is well-defined and easy to compute compared to the turbulent POD, because it is closed.
Note that for the POD bases, e0 and e2 have the same signs in the streaks but opposite signs in the rolls. That's in line with Wally's contention that the rolls and streaks are dynamically independent. Linear superposition of e0 and e2 can produce roll-streak combinations with different relative strengths of rolls and streaks. That such a superposition is present in the very lowest POD says that such variations in relative roll-streak strength happen frequently in the periodic orbit and turbulent flow.
Next thing is to compute state-space portraits of the orbits and turbulence for the periodic orbit paper. Also, it would be good to compare the e2 POD modes to the linear instability of the streaks a la self-sustaining process.