User Tools

Site Tools


gtspring2009:gibson:w03

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
gtspring2009:gibson:w03 [2009/03/26 09:30]
gibson
gtspring2009:gibson:w03 [2009/04/27 08:01]
predrag suggested a desymmetrized plot of Poincare sections
Line 132: Line 132:
 involved, so I'm glad I went ahead and did it. There is still tons to do. I will make notes tomorrow on  involved, so I'm glad I went ahead and did it. There is still tons to do. I will make notes tomorrow on 
 exactly how I did the calculations and provide you with the needed code. exactly how I did the calculations and provide you with the needed code.
-===== 2009-03-18 ​more slice and dice =====+===== 2009-03-18 ​EQ2 eigenvalues,​ full space  ​=====
  
-By the way, here are the eigenvalues for EQ2. ASA symmetry, for example, means+More slice and dice: the eigenvalues for EQ2. ASA symmetry, for example, means
 antisymmetry in s1 and s3 and symmetry in s2. antisymmetry in s1 and s3 and symmetry in s2.
  
Line 190: Line 190:
 These pictures are nice to look at but I'm not sure how much we can take away from them. You can't necessarily interpret an apparent fold in these sections as the classic nonlinear stretch and fold scenario, like in Rossler or Lorenz. It is true that an apparent fold, like the black part of the θ=0 section, shows that the unstable manifold has stopped stretching in the θ=0 direction and has returned back to lower values of (u, eθ). However, we can't attach too much meaning to the vertical axis. That "​fold"​ could be a very gentle curvature in the full space, and the portion of the unstable manifold that appears to be right above EQ2 (the origin) could be way, way far away from it, in directions we've projected out.  These pictures are nice to look at but I'm not sure how much we can take away from them. You can't necessarily interpret an apparent fold in these sections as the classic nonlinear stretch and fold scenario, like in Rossler or Lorenz. It is true that an apparent fold, like the black part of the θ=0 section, shows that the unstable manifold has stopped stretching in the θ=0 direction and has returned back to lower values of (u, eθ). However, we can't attach too much meaning to the vertical axis. That "​fold"​ could be a very gentle curvature in the full space, and the portion of the unstable manifold that appears to be right above EQ2 (the origin) could be way, way far away from it, in directions we've projected out. 
  
-===== 2009-03-26 ​ =====+===== 2009-03-26 ​EQ2 eigenvalues,​ S-invariant subspace ===== 
 + 
 +^ n ^ Re λ ^ Im λ ^ T-multip ^ symmetry ^ 
 +| 1,2 | 0.03252919 | 0.1070431 | 6.75 | SSS | 
 +| 3,4 | -0.01412161 | 0.05774779 | 0.447 | SSS | 
 +| 5,6 | -0.02429575 | 0.1479471 | 0.240 | SSS | 
 +| 7,8 | -0.02936259 | 0.1387545 | 0.178 | SSS | 
 +| 9,10 | -0.05782164 | 0.2838004 | 0.0336 | SSS | 
 +| 11,12 | -0.0620747 | 0.2420612 | 0.0262 | SSS | 
 +| 13,14 | -0.06761083 | 0.2166657 | 0.0189 | SSS | 
 +| 15,16 | -0.07723642 | 0.08462713 | 0.0107 | SSS | 
 +| 17,18 | -0.1008762 | 0.2608582 | 2.6e-03 | SSS | 
 +| 19,20 | -0.1023387 | 0.1500897 | 2.4e-03 | SSS | 
 +| 21,22 | -0.1068182 | 0.1773124 | 1.8e-03 | SSS | 
 +| 23,24 | -0.1090364 | 0.02625084 | 1.6e-03 | SSS | 
 +| 25,26 | -0.1114009 | 0.0739416 | 1.4e-03 | SSS | 
 +| 27 | -0.1118931 | 0 | 1.4e-03 | SSS | 
 +| 28,29 | -0.1134904 | 0.1401729 | 1.2e-03 | SSS | 
 + 
 +The 1,2 eigenvalues have a period of 58.7 and a multiplier of 6.75. The next SSS eigenvalue pair 3,4 has a multiplier of 
 +exp(-0.0141*58.7) = 0.437 in that period, which is not very strongly contracting. There'​s roughly a factor of three between 
 +the most expanding SSS direction and the next, most weakly stable. 
 + 
 +2009-04-06: added characteristic multiplier for period T=58.7 for each eigenvalue. %%mult = exp(57.8 Re λ)%% 
 + 
 +===== 2009-04-10 Poincare sections 3 ===== 
 + 
 +{{gtspring2009:​gibson:​eq2poincare_ubef2_1.png?​300}}{{gtspring2009:​gibson:​eq2poincare_ubef2_2.png?​300}}{{gtspring2009:​gibson:​eq2poincare_ubef2_3.png?​300}}{{gtspring2009:​gibson:​eq2poincare_ubef2_4.png?​300}} 
 + 
 +The figures show the unstable manifold of EQ2 (the Nagata upper branch) projected onto an orthogonal 
 +basis set spanning the leading three eigenfunctions in the S-symmetric subspace, whose eigenvalues  
 +are listed above. e1, e2 span the complex instability λ = 0.0325 + i 0.1070 (n=1,2 in the above table).  
 +e3 is the component of the real part of the most weakly contracting eigenfunction (n=3,4 above,  
 +eigenvalue λ = -0.0141 + i 0.0577) that is orthogonal to e1,e2. That sets an arbitrary choice of  
 +orientation within the plane of the n=3,4 oscillation;​ it would perhaps be better to make the third 
 +coordinate in the plots below be the norm of the projection in this plane, or even the distance from 
 +u(t) to the e1,e2 plane. I will try that later. 
 + 
 +The green trajectories are in the unstable manifold of EQ2, and the blue orbit is P47p81. The short, thick 
 +black lines are the real and imaginary parts of the n=1,2 eigenfunctions. The short think red line is the 
 +real part of the n=3,4 eigenfunction. The long thin black lines are intersections of the unstable manifold  
 +with Poincare sections placed at θ = 0, pi/4, pi/2, 3pi/4, etc. The θ = 0 section is marked with small 
 +circles and the θ = pi/4 with dots.  
 +===== 2009-04-13 Poincare sections 4 ===== 
 + 
 +{{gtspring2009:​gibson:​eq2poincare_ubef2_3.png?​400}}{{gtspring2009:​gibson:​eq2poincare_ubef2_section_zero.png?​400}}  
 + 
 +Here we have the poincare sections at θ=0 (black, and θ=π in green). The left plot 
 +repeats one shown immediately above, the right plot is the same but showing only the 
 +intersections with the Poincare section. The discontinuities in the unstable manifold  
 +are due to factoring out the 4th-order discrete translation symmetry. For any given state  
 +u(t), we apply whichever half-cell translation τ in {1,  τ<​sub>​x</​sub>,​ τ<​sub>​z</​sub>,​ τ<​sub>​xz</​sub>​} puts τ u(t) in the first quadrant of the [[gtspring2009:​gibson:​w03#​inctro|four-fold symmetric pictures above]], that is meets the conditions  
 + 
 +<​latex>​ $ \begin{align*} 
 +(1 + \tau_{x} - \tau_z - \tau_{xz} u_{EQ2}, \; \tau u(t)) & \geq 0 \\ 
 +(1 - \tau_{x} + \tau_z - \tau_{xz} u_{EQ2}, \; \tau u(t)) & \geq 0 
 +\end{align} $ </​latex>​ 
 + 
 +{{gtspring2009:​gibson:​eq2poincare_ubef2_2.png?​400}}{{gtspring2009:​gibson:​eq2poincare_ubef2_section_pi2.png?​400}} 
 + 
 +Same as above but with the Poincare section at θ=π/2 (black) and θ=3π/2 (green).  
 + 
 +The θ=0 section above looks potentially fruitful. It could be that the unstable manifold is 
 +folding back on itself due to the desymmetrization. Or it could be an artifact of the projection due to the arbitrary choice of the third (vertical) axis in the plots, e3. I will plot a few other projections of the same section to see. // 2009-04-13 11:40 EST// 
 + 
 + 
 +{{gtspring2009:​gibson:​eq2poincare_ubef2_section_zeroB.png?​400}}{{gtspring2009:​gibson:​eq2poincare_ubef2_section_zeroC.png?​400}} 
 + 
 +Ok, some Poincare sections at θ=0 orientation,​ with vertical axis (u, e3) on left and 
 +(u, e4) on right. e3, e4 are the orthogonal components of the real/​imaginary parts of 
 +the n=3,4 eigenfunction,​ respectively. ​ I have marked a few points on the unstable manifold with symbols so they can be compared between plots. Next is to assign arclengths along 
 +unstable manifold and plot a return map. This involves some more arbitrary decisions, namely, when to continue increasing arclength and when to project back onto previous portions of the unstable manifold. ​ // 2009-04-13 12:07 EST// 
 + 
 + 
 +{{gtspring2009:​gibson:​eq2poincare_ubef2_section_zeroD.png?​400}} 
 + 
 +The plot above shows a number of points on the unstable manifold and their iterates under the return map.  
 +A triangle symbol maps into the next triangle out along the unstable manifold, etc. The symbols are, in sequence, 
 +triangle, square, diamond, star, circle. The first iterates of the symbols are all clustered just to the right of 
 +EQ2 (at the origin). One return later they are spread from the triangle at (0.03, 0.03) to the circle at (-0.01, 0.08). 
 +The next return, they are all over the map, and it no longer makes sense to draw lines connecting them (though I have 
 +done so to show you that I am insane).  
 + 
 +This is bad news, I believe. In order for the approximate return map s^(n+1) = f(s^n) (where s is pseudo-arclength  
 +along the unstable manifold) to produce good guesses for periodic orbits, it needs to intersect the identity.  
 +has to intersect the identity. In this plot, that would correspond to some portion of the lower part of the fold 
 +to be mapped into the upper part. But the symbols show us this doesn'​t happen.  
 + 
 +Corrections?​ Suggestions?​  
 + 
 +//John Gibson 2009-04-13 14:52 EST// 
 + 
 +===== Desymmetrized Poincare sections ===== 
 + 
 +Can you plot these sections in desymmetrized space of Gibson //et al?// 
 + 
 +{{:​gtspring2009:​howto:​gtspring2009:​2008-spring-p7.png?​800}} 
 + 
 +It should be easier to visualize them without discontinuities... ​ --- //​[[predrag.cvitanovic@physics.gatech.edu|Predrag Cvitanovic]] 2009-04-27 07:53//
  
  
-Eigenvalues of EQ2 in the S-symmetric subspace 
  
gtspring2009/gibson/w03.txt · Last modified: 2010/02/02 07:55 (external edit)