Hello All !
I am trying to validate the Barkley JFM paper on the Biglobal floquet analysis of cylinder. I am doing the simulations at the following Re and Beta values to find the neutral modes and create a curve of neutral modes (as done in the paper) to identify the critical Re and Beta.

1. Re 190; Beta = 1.4 and 1.8.

2. Re 200, Beta = 1.2 and 2.0

For the above cases, the leading eigenvalue magnitude should be approximately 1. But the solver is producing eigenvalues with magnitude greater than 1.

For example, converged eigenvalues for Re=190 and Beta=1.4

Magnitude   Angle       Growth      Frequency
EV:  0     1.28858    0.886868    0.202835    0.709494
EV:  1     1.28858   -0.886868    0.202835   -0.709494
EV:  2     1.13585    0.342226    0.101906    0.273781
EV:  3     1.13585   -0.342226    0.101906   -0.273781

For the base flow, I am doing DNS for around 40 shedding periods (~200s) and then storing 32 snapshots for the 41st shedding period. DNS solution is converged and validated against CD, and St values. The DNS mesh has around (~2050 Quads and 7th order polynomial). I also made sure that the 1st snapshot is not equal to the 32nd snapshot as mentioned in the user-guide. Later, I am interpolating the DNS snapshots to a coarser mesh (~550 Quads and 7th order polynomial) using FieldConvert. These interpolated chk files are used as baseflow for the Floquet Analysis. 

I would be highly grateful, if anyone could help me out with this. I have attached my case setup files for Re=190 and Beta=1.4 case. 

Case Files: https://drive.google.com/drive/folders/1-0fbs1mxAuQVKfYpbSpj5PBod2aQjKKS?usp=sharing