<div dir="ltr"><div><div>Hi JC, I double-checked and am using v (rather than V) only.<br><br></div>Hi Philipp, Thanks for your answer. Can someone else confirm this too? Also, is there a way to have zero velocities in the interior, just like in the youtube simulation?<br><br></div>Thanks, <br>Ajay<br><div><br><br><div><div>Could it be that he uses V instead of v in the rea file?<br>
<br>++<br>JC<br>
<br>Le 14 sept. 2017 18:41, <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.<wbr>gov</a>> a ?crit :<br>
<br>> Hi,<br>> Not sure if I understand your question well, but I guess what happens<br>> during the first tim<span class="gmail-aBn" tabindex="0"><span class="gmail-aQJ">e step is that t</span></span>hrough the continuity equation and the<br>> incompressibility constraint you will have velocity about one everywhere.<br>> 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">nek5000-users@lists.mcs.anl.<wbr>gov</a> wrote:<br>><br>>> Dear Neks,<br>>><br>>> I am trying to run a simulation, replicating the Nek5000 run from this<br>>> 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<br>>> seen from the video, the velocity magnitude reasonably away from the inlet<br>>> is zero at starting times. However, in all of my simulations, I see the<br>>> velocity magnitude as one at all the interior points away from the inlet<br>>> from the beginning (the inlet BC ux=1.0 I am applying is somehow propagated<br>>> to all interior elements). Am I doing something wrong somewhere, or is this<br>>> the expected behavior? Please see my userbc() and useric() from my .usr<br>>> 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<br>>> constraints on the elements (like 'W', 'O', 'E'),  I think userbc() would<br>>> not run (please correct me if I am wrong, I am learning here). Since the<br>>> above functions wouldn't work as expected, I changed userchk() and removed<br>>> 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</div></div></div></div>