[Nek5000-users] Time-dependant boundary condition

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

In userbc:

      ux = cos(pi*y)*cos(2.*pi*time)
      uy = 0.

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: Monday, June 09, 2014 2:30 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Time-dependant boundary condition

How should I impose time-dependant BC for velocity on the left vertical boundaty for a 2D box geometry.
The velocity profile should be
Vx=cos(pi*y)*cos(2*pi*t)
Vy=0

To illustrate the idea here are the pictures.
 Instead of constant velocity Vx=c on the left (inflow) boundary of a channel
 V __________________
  |
  |
  |
  |__________________

I would impose Vx=cos(pi*y)*cos(2*pi*t)
so for t = 0 velocity profile is
 ____________________
\
 \
  \
   \_________________
for t=1 velocity profile is
    _________________
   /
  /
 /
/____________________

and so on.

I know it should be possible in nek5000, but still not very familiar with syntax and fortran.

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