[Nek5000-users] Time-dependant boundary condition
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jun 9 02:30:07 CDT 2014
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.
More information about the Nek5000-users
mailing list