<html><head></head><body><div dir="auto">How can i leave the mailing list?<br><br></div>
<div dir="auto"><!-- tmjah_g_1299s -->Inviato da <!-- tmjah_g_1299e --><a href="http://www.bluemail.me/r?b=10455"><!-- tmjah_g_1299s -->BlueMail<!-- tmjah_g_1299e --></a><!-- tmjah_g_1299s --> <!-- tmjah_g_1299e --></div>
<div class="gmail_quote" >Il giorno 14 set 2017, alle ore 18:41, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> ha scritto:<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<pre class="blue">Hi,<br>Not sure if I understand your question well, but I guess what happens <br>during the first time step is that through the continuity equation and <br>the incompressibility constraint you will have velocity about one <br>everywhere. This is just to maintain mass conservation.<br><br>Philipp<br><br>On 2017-09-14 18:29, nek5000-users@lists.mcs.anl.gov wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 1ex 0.8ex; border-left: 1px solid #729fcf; padding-left: 1ex;"> 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">https://www.youtube.com/watch?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 <br> inlet is zero at starting times. However, in all of my simulations, I <br> see the velocity magnitude as one at all the interior points away from <br> the inlet from the beginning (the inlet BC ux=1.0 I am applying is <br> somehow propagated to all interior elements). Am I doing something wrong <br> somewhere, or is this the expected behavior? Please see my userbc() and <br> 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 <br> constraints on the elements (like 'W', 'O', 'E'), I think userbc() <br> would not run (please correct me if I am wrong, I am learning here). <br> Since the above functions wouldn't work as expected, I changed userchk() <br> 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><hr><br> Nek5000-users mailing list<br> Nek5000-users@lists.mcs.anl.gov<br> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br> <br></blockquote><hr><br>Nek5000-users mailing list<br>Nek5000-users@lists.mcs.anl.gov<br><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br></pre></blockquote></div></body></html>