<div dir="ltr"><div>Thanks a lot Paul.<br><br>That's pretty much the idea I had, but I just wanted to make sure.<br><br></div>Regards,<br>JC<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-01-30  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span>:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear JP,<br>
<br>
In perturb.f, the driver reads:<br>
<br>
      do jp=1,npert<br>
<br>
         if (nid.eq.0.and.igeom.eq.2) write(6,1) istep,time,jp<br>
   1     format(i9,1pe14.7,' Perturbation Solve:',i5)<br>
<br>
         call perturbv (igeom)<br>
<br>
      enddo<br>
<br>
      jp=0   ! set jp to zero, for baseline flow<br>
<br>
<br>
So, the appropriate discriminator in userbc would be:<br>
<br>
      if (jp.eq.0) then<br>
<br>
         ux=1<br>
<br>
      else<br>
<br>
         ux=0<br>
<br>
      end<br>
<br>
Best regards,<br>
<br>
Paul<div><div class="h5"><br>
<br>
<br>
<br>
<br>
On Thu, 30 Jan 2014, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Hi Nek's,<br>
<br>
I have a question regarding the linearised Navier-Stokes solver while the<br>
underlying base flow is time-dependant. As far as I understand, one should<br>
simply have to set param(31) = 1 to achieve this (and p31 = -1 for a steady<br>
state base flow). What about the boundary conditions however?<br>
<br>
For instant, I want to investigate the linear dynamics of perturbations<br>
evolving onto a lid-driven cavity flow exhibiting a periodic limit cycle I<br>
have computed. How could I impose in userbc:<br>
<br>
ux = 1<br>
uy = 0<br>
uz = 0<br>
<br>
for the periodic base flow time-evolution, and in the meantime:<br>
<br>
ux = 0<br>
uy = 0<br>
uz = 0<br>
<br>
for the perturbation?<br>
<br>
Cheers,<br>
JC<br>
<br>
-- <br>
Jean-Christophe Loiseau<br></div></div>
Homepage <<a href="https://sites.google.com/site/loiseaujc/" target="_blank">https://sites.google.com/<u></u>site/loiseaujc/</a>><br>
<br>
</blockquote>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">Jean-Christophe Loiseau<br><a href="https://sites.google.com/site/loiseaujc/" target="_blank">Homepage</a><br></div>
</div>