[Nek5000-users] turbulent 3d box

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

Dear Paul,

Thanks for the example. However, it runs but something seems to be 
wrong: the value of the velocity magnitude is too high to my experience. 
(it should be about |f| where f is the forcing  (the correct max value 
is sqrt(3)).
Two questions:
1. In b3d.rea  the parameter p2 (viscos) is 10935.  Is't too high? I 
used values such as 1/12 and smaller ones.
2. in b3d.f there is a routine userbc(..)
where the parameter ux=uy=uz=0. does it mean that the velocity on the 
boundaries is set to zero?
This is not the case which I solved (in previous DNS) and also it 
conflicts with the forcing. I think that this condition should be removed.

Thanks very much for your help.
Barak

On 12/18/2013 03:29 AM, nek5000-users at lists.mcs.anl.gov wrote:
>
> Hi Barak,
>
> I added a 3d periodic box under "3dbox" in the examples
> directory.
>
> Paul
>
>
> Dear nek's,
> I would like to run the following problem:
> a 3D periodic box, small viscosity and the following forcing:
>
> fx=cos(y)+sin(z)
> fy=cos(z)+sin(x)
> fz=cos(x)+sin(y)
>
> Do you have a similar example?
>
> Many thanks,
> Barak
>
>
> _______________________________________________
> 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
>