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>