Hi All, I encounter some problems with Nektar++. I want to calculate the eigenmode of a period field using BiGlobal Floquet analysis. 1, I have gotten a series base flow file (xx_1.bse, xx_2.bse ......). The Residual error of base period flow (long enough) , L 2 error (variable u), is about 10, and others are o(10^-1). *Is the result **convergence*? (i also got the lift force of the bluff body, and it shows period smoothly ) 2, The total structural quadrangular element of the 2D mesh is about 14000, order of P is 8, the Re number is 4000. The external flow is similar to the example in User Guide 8.6.12 (version 4.0.1), but i simulate a plan-like bluff body instead of cylinder. when i run the "mpirun -np 4 IncNavierStokesSolver XX.xml " , the XX_0.chk folder is generate, but it stay all day long in that state (no error, no respones). Simultaneously, I run "IncNavierStokesSolver XX.xml" ,i.e., using single processor, it works well (although slowly). Anyone can tell me why and what should I do? 3, How can I get the corresponding eigenmode of eigenvalue in .evl file? 4, I am a newbie, would anyone gives me some advises in my simulation? I will appreciate! PS, the part of .xml content (base flow): <I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" /> <I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" /> <I PROPERTY="EvolutionOperator" VALUE="Nonlinear" /> <I PROPERTY="Projection" VALUE="Galerkin" /> <I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" /> </SOLVERINFO> <PARAMETERS> <P> TimeStep = 0.0001 </P> <P> NumSteps = 100000 </P> <P> IO_CheckSteps = 10000 </P> <P> IO_InfoSteps = 100 </P> <P> Re = 4000 </P> <P> Kinvis = 0.775/4000 </P> </PARAMETERS>