[Nek5000-users] how to define boundary conditions with a specified function.

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

Wei,

I would just do

          ux=0
          uy=0
          uz=0
          eps = 1.e-4
          if (abs(y-1.).le.eps) ux=2*x

That will work.

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: Tuesday, June 17, 2014 10:33 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] how to define boundary conditions with a specified     function.

Hello everyone,
Now I have a problem.  that how to define boundary conditions in the .usr file.
I know that I can write some code in the usrbc()

what I want do is like this

if (y .eq. 1) then
        u=2*x
endif
where y is the value of boundary coodinates .

what I wrote is like
if (ym1(1,1,11) .eq. 1) then
        ux=2*x
        uy=0
        uz=0
endif

it can be compiled, but not imposed. so what should I do? How can I know the relation between ux(i,j,k,e) and index informations? thanks.

Wei

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