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

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

Hi all,

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.

Sincerely,

On Tue, Jan 2, 2018 at 4:10 AM <nek5000-users at lists.mcs.anl.gov> wrote:

> You set ux=1.0 after your if-test...?
>
> On 12/31/2017 07:34 PM, 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
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> >
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-- 

*Amitvikram Dutta*

Graduate Research Assistant

Fluid Mechanics Research Lab

Multi-Physics Interaction Lab

University of Waterloo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180102/623699d3/attachment.html>