[Nek5000-users] Different velocity values in different 'v' inlets

[nek5000-users at lists.mcs.anl.gov]

What is the variable e in userbc?

On 31 Dec 2017, at 19:35, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:

Hi all,

I'm trying to simulate a shear-layer which involves two different jets emerging at different velocities. The simulation is 2D for now. The width of the first jet extends from y=0.0 to y=1.0 and the second extends from y=1.0 to y=3.0

I generated a mesh using genbox with 'v' BCs for both inlets.

In the .usr file I have the following form for usrbc

 subroutine userbc (ix,iy,iz,iside,ieg)
c     NOTE ::: This subroutine MAY NOT be called by every process
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'
      if (ym1(i,j,k,e).lt.1) then
          ux=10.0
      elseif (ym1(i,j,k,e).gt.1) then 
          ux=1.0
      endif
      ux=1.0
      uy=0.0
      uz=0.0
      temp=0.0
      return
      end

As you can see I try assigning ux = 10.0 and 1.0 to the different zones. However the results show a uniform value of 1.0 across both the inlets.
What am I missing here?

Sincerely,
-- 
Amitvikram Dutta

Graduate Research Assistant

Fluid Mechanics Research Lab

Multi-Physics Interaction Lab

University of Waterloo


_______________________________________________

Nek5000-users mailing list

Nek5000-users at lists.mcs.anl.gov <mailto:Nek5000-users at lists.mcs.anl.gov> 

https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171231/3c53e7ad/attachment.html>