[Nek5000-users] turbulent 3d box

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

Hi Barak,

Just take the box file from turbChannel
and set nelx nely nelz to (say) -8 -8 -8
with periodic BCs  P  ,P  ,P  ,P  ,P  ,

Then, use a blank .usr file and set ffx fry ffz to be the
formulas you have below in userf() in the .usr file.

Finally, add the following to usrdat2()

      a=0
      b=2*pi
      call rescale_x(xm1,a,b)
      call rescale_x(ym1,a,b)
      call rescale_x(zm1,a,b)

which will put your box on [0,2pi]^3.

Paul


----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Sunday, December 15, 2013 4:14:22 PM
Subject: [Nek5000-users] turbulent 3d box

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