[Nek5000-users] userbc to set BC from velocity field computed in userchk

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

Hi Arnold,

Correct syntax would be:

   subroutine userbc(i,j,k,f,eg)
   include 'SIZE'
   include 'TOTAL'
   YOUR COMMON BLOCK HERE (to pass boundx...)

   integer e,eg,f

   e = gllel(eg)   !  To get the local element number

   ux = boundx(i,j,k,e)
   uy = boundy(i,j,k,e)
   uz = boundz(i,j,k,e)

   return
   end


Paul



   
________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [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: Thursday, October 06, 2016 10:37 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] userbc to set BC from velocity field computed in       userchk

Good afternoon Neks,

I am trying to set boundary conditions from a velocity field computed in
the userchk subroutine. The nature of the velocity field is defined as
follows :

real boundx(lx1,ly1,lz1,lelt)
real boundy(lx1,ly1,lz1,lelt)
real boundz(lx1,ly1,lz1,lelt)

I was wondering what was the proper syntax to set the user defined
Dirichlet velocity; is the following syntax in the userbc correct ?

ux=boundx(ix,iy,iz,iside,ie)
uy=boundy(ix,iy,iz,iside,ie)
uz=boundz(ix,iy,iz,iside,ie)

Thank you for your time,

Best wishes,

Arnold
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users