<div dir="ltr">Hi all,<div><br></div><div>There's an error in the above code with ux = 1.0 having been left after the loop. However, the real issue was, I had to set the loop conditions in both usrbc and usric before starting the sim. That resolved the issue for me.</div><div><br></div><div>Sincerely,</div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Jan 2, 2018 at 4:10 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You set ux=1.0 after your if-test...?<br>
<br>
On 12/31/2017 07:34 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
> Hi all,<br>
><br>
> I'm trying to simulate a shear-layer which involves two different jets<br>
> emerging at different velocities. The simulation is 2D for now. The<br>
> width of the first jet extends from y=0.0 to y=1.0 and the second<br>
> extends from y=1.0 to y=3.0<br>
><br>
> I generated a mesh using genbox with 'v' BCs for both inlets.<br>
><br>
> In the .usr file I have the following form for usrbc<br>
><br>
>  subroutine userbc (ix,iy,iz,iside,ieg)<br>
> c     NOTE ::: This subroutine MAY NOT be called by every process<br>
>       include 'SIZE'<br>
>       include 'TOTAL'<br>
>       include 'NEKUSE'<br>
>       if (ym1(i,j,k,e).lt.1) then<br>
>           ux=10.0<br>
>       elseif (ym1(i,j,k,e).gt.1) then<br>
>           ux=1.0<br>
>       endif<br>
>       ux=1.0<br>
>       uy=0.0<br>
>       uz=0.0<br>
>       temp=0.0<br>
>       return<br>
>       end<br>
><br>
> As you can see I try assigning ux = 10.0 and 1.0 to the different zones.<br>
> However the results show a uniform value of 1.0 across both the inlets.<br>
> What am I missing here?<br>
><br>
> Sincerely,<br>
> --<br>
><br>
> *Amitvikram Dutta*<br>
><br>
> Graduate Research Assistant<br>
><br>
> Fluid Mechanics Research Lab<br>
><br>
> Multi-Physics Interaction Lab<br>
><br>
> University of Waterloo<br>
><br>
><br>
><br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> <a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><div dir="ltr">-- <br></div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p><p>Multi-Physics Interaction Lab</p>
<p dir="ltr">University of Waterloo</p>
</div></div>