Index of /turbdata/raw_database/Boundary_Layer/Low-Reynolds

[ICO]NameLast modifiedSizeDescription
[PARENTDIR]Parent Directory   -  
[DIR]HR-Reth-620-2140/ 2015-03-24 14:05 -  
========================================================================
=                                                                      =
=            STATISTICS FROM DNS OF ZERO PRESSURE GRADIENT 	       =
=	             	   BOUNDARY LAYER			       =
=                   IN A VERY LARGE NUMERICAL BOX		       =
=                          Re_theta = 600-2000			       =
=								       =
========================================================================

This directory tree contain some of the restart files of a direct
numerical simulations of fully developed ZPG Boundary Layer flow
(Re_theta=600-2000). Details are offered below.


Authors:                Mark P. Simens (markphilsimens@gmail.com)
                        Sergio Hoyas (serhocal@mot.upv.es)
                        Javier Jimenez (jimenez@torroja.dmt.upm.es)

-----------------------------------------------------------------------
                     GENERAL INFORMATION:
-----------------------------------------------------------------------
Formulation:            Primitive variables.
Reynolds number:        Re_theta=600-2000
Domain size:            (Lx,Ly,Lz)/theta=(535,29,88)
Resolution:             (dx+,dy+,dz+)=(6.1,0.32,4.1)  
			*Taken at the middle of the box.

Number of fields:       298. Single precission.
Number of files & size: 298. Each file of size 34 GB
Total Size:             10 TB.
Name:                   tbl.2000.xxx
Format:                 Binary (stored as YXZ)


-----------------------------------------------------------------------
                     STRUCTURE OF THE FILES:
-----------------------------------------------------------------------
The files are originally written in binary format. We can provide with
a parallel code to convert those files to HDF5. HDF5 is supported by
most of the mathematical programs. For programming languages,
libraries can be downloaded for C, C++, Fortran 90, and Java
interfaces.

The primitive variables {u,v,w,p} are stored in Fourier space as:
(Ny,Nx,Kz) ---> ( Physical - Physical - Fourier )

For instance, one file is read (pseudo-code) as follow in Fortran:

**********************************************************
integer:: nx,ny,nz,k,i,j,l
real*8:: tiempo,cfl,Re,ax,ay,az
real*8:: x(0:nx+1),y(0;ny+1),um(1:ny+1)
real*4,dimension(1:ny  ,1:nx,2,1:nz/3-1):: v,p
real*4,dimension(1:ny+1,1:nx,2,1:nz/3-1):: u,w

open(100,file='tbl.2000.xxx',form='unformatted',status='old')
!Re=1/nu; ax=Lx/pi, ay=Ly/pi, az=Lz/2pi (L in code units)
read(100) tiempo,cfl,Re,ax,ay,az,nx,ny,nz,timeinit,dt

read(100) (x(i), i=0,nx+1)                   !Grid in X
read(100) (y(i), i=0,ny+1),(um(i), i=1,ny+1) !Grid in Y and <U> @ the Rescaling plane


!k_begin=1; k_end=[2*(Nz/3)]/2-1 !2/3 Dealiasing Rule
do k=1,nz/3-1  !(l=1,2 stand for Real,Imaginary parts)
   read(100) (((u(j,i,l,k), j=1,ny+1),i=1,nx),l=1,2), &
        &    (((v(j,i,l,k), j=1,ny  ),i=1,nx),l=1,2), &
        &    (((w(j,i,l,k), j=1,ny+1),i=1,nx),l=1,2), &
        &    (((p(j,i,l,k), j=1,ny)  ,i=1,nx),l=1,2)
enddo
**********************************************************


Detailed information about the code and the collocation of the
privitive variables on the grid can be consulted in,

   Simens, M.P., "The study and control of wall bounded flows"
   2008, (PhD Thesis).
   http://torroja.dmt.upm.es/pubs/2008/MarkThesis.pdf


-----------------------------------------------------------------------
                    RELATED PUBLICATIONS:
-----------------------------------------------------------------------
[1] M. P. Simens J. Jimenez, S. Hoyas, Y. Mizuno,
    "A high-resolution code for turbulent boundary layers"
    Journal of Computational Physics
    Volume 228, Issue 11, 20 June 2009, Pages 4218-4231

[2] J. Jimenez, S. Hoyas, M. P. Simens, Y. Mizuno
    "Turbulent boundary layers and channels at moderate
    Reynolds numbers"
    Journal of Fluid Mechanics 657, 2010, Pages 335-360


-----------------------------------------------------------------------
                     FUNDING:
-----------------------------------------------------------------------
This work was supported in part by the CICYT grant TRA2006-08226, by a
contract from Industria de Turbopropulsores S.A., and by the EU FP6
Wallturb Strep AST4-CT-2005-516008. The computations were made
possible by generous grants of computer time from the Barcelona
supercomputing centre and from the Cesvima centre of the UPM, as well
as by the equally generous collaboration of the Port d’Informacio
Cientifica (PIC), who lent their mass storage facilities to archive
raw data.  Special thanks are due to G. Hauët and S. Vaux for running
some of the early evaluation runs for this code. MPS and GH were
supported in part by the EU FP5 Training and Mobility Network
HPRN-CT-2002-00300, and SH, SV and YM by the Spanish Ministry of
Education and Science, under the Juan de la Cierva program.


-----------------------------------------------------------------------
                     CONTACT:
-----------------------------------------------------------------------
Please report any problems with the data, queries or doubts to
J.A. Sillero, sillero@torroja.dmt.upm.es