<div>Hi Nek's,<br clear="all"></div><div><br></div><div>I am still struggling with the perturbation mode. </div><div>Cavity_20x20x10_6x6x6_R1000_25D0.f00001 contains my base flow. It is a 2D lid-driven cavity extruded in the third direction (z).<div>

Perturbation_20x20x10_6x6x6_25D0.f00001 contains my divergence-free perturbation. It is subject to no-slip conditions at the walls and periodic BC at zmin and zmax.</div></div><div><br></div><div>Following is how I set it up.</div>

<div><br></div><div>**********************************************</div><div><br></div><div>In my blah.rea<br></div><div><br></div><div style="margin-left: 40px">-1.                        p31 NPERT: #perturbation modes</div>

<div><br></div><div style="margin-left: 40px">[...]</div><div><div style="margin-left: 40px"><br></div><div style="margin-left: 40px">   2 PRESOLVE/RESTART OPTIONS  *****</div><div style="margin-left: 40px">Cavity_20x20x10_6x6x6_R1000_25D0.f00001</div>

<div style="margin-left: 40px">Perturbation_20x20x10_6x6x6_25D0.f00001</div></div><div style="margin-left: 40px"><div><br></div></div><div>**********************************************<br></div><div><br></div><div>In SIZE:</div>

<div><br></div><div style="margin-left: 40px"><div>parameter (lpelv=lelv,lpelt=lelt,lpert=1)  ! perturbation<br>      parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1)     ! array sizes<br>      parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)<br>

</div></div><div style="margin-left: 40px"><div><br></div></div><div>**********************************************<br></div><div><br></div><div>I have tried with various initial amplitude for the perturbation, and various time steps. However, I always have a very small CFL and after six timesteps I get:</div>

<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Failed in HMHOLTZ: VELX NAN<br></blockquote><div><br></div><div>Same for vely and velz. I may have missed a key point but being not particularly familiar with Nek numerics yet I'm not sure where I missed it.</div>

<div><br></div><div>Best regards,<br></div><div>-- <br></div>Jean-Christophe<br>