[Nek5000-users] problem with applying boundary conditions
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Thu Sep 14 12:35:51 CDT 2017
Hi,
if you have unit velocity at the inflow, then due to mass conservation
each cross section in the channel needs to have the same mass flux. I am
not sure what boundary conditions were used in this youtube movie; it
appears it was not simply one everywhere, but perhaps slowly ramped up
from zero. In this way your flow would slowly accelerate.
Best,
Philipp
On 2017-09-14 19:21, nek5000-users at lists.mcs.anl.gov wrote:
> Hi JC, I double-checked and am using v (rather than V) only.
>
> Hi Philipp, Thanks for your answer. Can someone else confirm this too?
> Also, is there a way to have zero velocities in the interior, just like
> in the youtube simulation?
>
> Thanks,
> Ajay
>
>
> Could it be that he uses V instead of v in the rea file?
>
> ++
> JC
>
> Le 14 sept. 2017 18:41, <nek5000-users at lists.mcs.anl.gov
> <mailto:nek5000-users at lists.mcs.anl.gov>> a ?crit :
>
> > Hi,
> > Not sure if I understand your question well, but I guess what happens
> > during the first time step is that through the continuity equation
> and the
> > incompressibility constraint you will have velocity about one everywhere.
> > This is just to maintain mass conservation.
> >
> > Philipp
> >
> > On 2017-09-14 18:29, nek5000-users at lists.mcs.anl.gov
> <mailto:nek5000-users at lists.mcs.anl.gov> wrote:
> >
> >> Dear Neks,
> >>
> >> I am trying to run a simulation, replicating the Nek5000 run from this
> >> youtube video:
> >> https://www.youtube.com/watch?v=RmmN-KFSvAY
> <https://www.youtube.com/watch?v=RmmN-KFSvAY>
> >>
> >> The inlet BC is constant velocity (which is set to 1.0), and as can be
> >> seen from the video, the velocity magnitude reasonably away from the
> inlet
> >> is zero at starting times. However, in all of my simulations, I see the
> >> velocity magnitude as one at all the interior points away from the inlet
> >> from the beginning (the inlet BC ux=1.0 I am applying is somehow
> propagated
> >> to all interior elements). Am I doing something wrong somewhere, or
> is this
> >> the expected behavior? Please see my userbc() and useric() from my .usr
> >> file below:
> >>
> >> subroutine userbc (ix,iy,iz,iside,eg)
> >> include 'SIZE'
> >> include 'TOTAL'
> >> include 'NEKUSE'
> >> integer ix,iy,iz,eg,e,i,n
> >>
> >> CC "My inlet is at x=0"
> >>
> >> if (x.le.1e-4) then
> >> ux=1.0
> >> endif
> >>
> >> return
> >> end
> >>
> >> subroutine useric (ix,iy,iz,eg)
> >> include 'SIZE'
> >> include 'TOTAL'
> >> include 'NEKUSE'
> >>
> >> ux=0
> >> uy=0
> >>
> >> return
> >> end
> >>
> >> I believe userbc() is executed for only 'v' in my case. For the other
> >> constraints on the elements (like 'W', 'O', 'E'), I think userbc()
> would
> >> not run (please correct me if I am wrong, I am learning here). Since the
> >> above functions wouldn't work as expected, I changed userchk() and
> removed
> >> the above if statements in userbc(). My userchk() is as follows:
> >>
> >> subroutine userchk
> >> include 'SIZE'
> >> include 'TOTAL'
> >> integer i,n
> >>
> >> n = nx1*ny1*nz1*nelv
> >>
> >> if (istep.eq.0) then
> >> do i=1,n
> >> if (xm1(i,1,1,1).le.1e-4) then
> >> vx(i,1,1,1)=1
> >> endif
> >> enddo
> >> endif
> >>
> >> return
> >> end
> >>
> >> This didn't work either. Can somebody point me in the right direction?
> >>
> >> Thanks,
> >> Ajay
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
More information about the Nek5000-users
mailing list