[Nek5000-users] Time-dependant boundary condition

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Jun 11 02:50:52 CDT 2014


Thanks, I will try

Van

> 09.06.2014, 16:06:05 пользователь  (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
> _______________________________________________
> 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