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

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

Sorry -- I didn't see the code in the first email.


What I would do is the following:


 subroutine userbc (ix,iy,iz,iside,ieg)
c     NOTE ::: This subroutine MAY NOT be called by every process
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      ux=1.0
      uy=0.0
      uz=0.0
      temp=0.0

      if (y.lt.1) ux=10.   ! y is passed in through NEKUSE


      return
      end

________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Tuesday, January 2, 2018 8:21:26 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Different velocity values in different 'v' inlets

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<mailto: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<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
>
_______________________________________________
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
--

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/0c3ba223/attachment.html>