<div dir="auto">Could it be that he uses V instead of v in the rea file?<div dir="auto"><br></div><div dir="auto">++</div><div dir="auto">JC</div></div><div class="gmail_extra"><br><div class="gmail_quote">Le 14 sept. 2017 18:41,  <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> a écrit :<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
Not sure if I understand your question well, but I guess what happens during the first time step is that through the continuity equation and the incompressibility constraint you will have velocity about one everywhere. This is just to maintain mass conservation.<br>
<br>
Philipp<br>
<br>
On 2017-09-14 18:29, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Neks,<br>
<br>
I am trying to run a simulation, replicating the Nek5000 run from this youtube video:<br>
<a href="https://www.youtube.com/watch?v=RmmN-KFSvAY" rel="noreferrer" target="_blank">https://www.youtube.com/watch?<wbr>v=RmmN-KFSvAY</a><br>
<br>
The inlet BC is constant velocity (which is set to 1.0), and as can be seen from the video, the velocity magnitude reasonably away from the inlet is zero at starting times. However, in all of my simulations, I see the velocity magnitude as one at all the interior points away from the inlet from the beginning (the inlet BC ux=1.0 I am applying is somehow propagated to all interior elements). Am I doing something wrong somewhere, or is this the expected behavior? Please see my userbc() and useric() from my .usr file below:<br>
<br>
       subroutine userbc (ix,iy,iz,iside,eg)<br>
       include 'SIZE'<br>
       include 'TOTAL'<br>
       include 'NEKUSE'<br>
       integer ix,iy,iz,eg,e,i,n<br>
<br>
CC "My inlet is at x=0"<br>
<br>
       if (x.le.1e-4) then<br>
         ux=1.0<br>
       endif<br>
<br>
       return<br>
       end<br>
<br>
       subroutine useric (ix,iy,iz,eg)<br>
       include 'SIZE'<br>
       include 'TOTAL'<br>
       include 'NEKUSE'<br>
<br>
       ux=0<br>
       uy=0<br>
<br>
       return<br>
       end<br>
<br>
I believe userbc() is executed for only 'v' in my case. For the other constraints on the elements (like 'W', 'O', 'E'),  I think userbc() would not run (please correct me if I am wrong, I am learning here). Since the above functions wouldn't work as expected, I changed userchk() and removed the above if statements in userbc(). My userchk() is as follows:<br>
<br>
       subroutine userchk<br>
       include 'SIZE'<br>
       include 'TOTAL'<br>
       integer i,n<br>
<br>
       n = nx1*ny1*nz1*nelv<br>
<br>
       if (istep.eq.0) then<br>
         do i=1,n<br>
           if (xm1(i,1,1,1).le.1e-4) then<br>
             vx(i,1,1,1)=1<br>
           endif<br>
         enddo<br>
       endif<br>
<br>
       return<br>
       end<br>
<br>
This didn't work either. Can somebody point me in the right direction?<br>
<br>
Thanks,<br>
Ajay<br>
<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
</blockquote>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
</blockquote></div></div>