Index of /raw_database/HST

[ICO]NameLast modifiedSizeDescription

[PARENTDIR]Parent Directory   -  
[DIR]re100_32/ 2016-06-23 17:41 -  
[DIR]re100_34/ 2016-06-24 05:02 -  
[DIR]re250_32/ 2017-04-30 16:28 -  

===========================================================================
          Statistically-stationary homogeneous shear turbulence 
            at several Reynolds numbers and box aspect ratios 
===========================================================================

The raw data is organized according to the Reynolds numbers based on Taylor
microscale and its box aspect ratios measured as (Lx/Lz,Ly/Lz).

---------------------------------------------------------------------------
		Nomenclature for the folders for each case
---------------------------------------------------------------------------

   Snapshots_hdf5:   HDF5 file containing {Omega_y, Phi=Lap(v)}.

The next files, when available, are written in binary:
  
   U_Snapshot : Streamwise velocity.
   V_Snapshot : Vertical velocity.
   W_Snapshot : Spanwise velocity.
   Q_Snapshot : The second invariant of the velocity gradient tensor.
   R_Snapshot : The third invariant of the velocity gradient tensor.
   Ox_Snapshot: Streamwise vorticity.
   Oy_Snapshot: Vertical vorticity.
   Oz_Snapshot: Spanwise vorticity.
   Uy_Snapshot: Vertical derivative of the streamwise velocity.

---------------------------------------------------------------------------
		STRUCTURE OF THE FILES:
---------------------------------------------------------------------------

/////////////////////////////////////////////////////////////// 
///////////////////////  Snapshots_hdf5 /////////////////////// 
/////////////////////////////////////////////////////////////// 

Contents of HDF5 file. For example,
# h5ls var.323.h5
  buff                     Dataset {4}
  dum                      Dataset {7}
  mpe                      Dataset {6}
  phi                      Dataset {2048, 1023, 2048}
  u00                      Dataset {2048}
  vor                      Dataset {2048, 1023, 2048}
  w00                      Dataset {2048}
  y                        Dataset {2048}

  mpe =(/mgalx,my,mgalz,numerop,mx,mz/)
  dum =(/time,Re,alp,Ly/pi,gam,s,chi/)
  buff=(/xoffb,xofft,zoffb,zofft/)

  mgalx: streamwise grid number after dealiasing
  my   : vertical grid number
  mgalz: spanwise grid number after dealiasing
  numerop: number of MPI tasks
  mx   : streamwise grid number before dealiasing, mx=2*(mgalx/3)
  mz   : spanwise grid number before dealiasing, mz=2*(mgalz/3)-1
  t    : time
  Re   : inverse of kinematic viscosity
  alp  : 2pi/Lx
  gam  : 2pi/Lz
  s    : constant mean shear, s=1.0
  chi  : rate of rotation, chi=0.0
  xoffb: S*Ly*t
  xofft:-S*Ly*t
  zoffb and zofft are zeros

  y   : vertical grid
  u00 : (kx,kz)=(0,0) of the streamwise velocity
  w00 : (kx,kz)=(0,0) of the spanwise velocity
  vor : omegay in Fourier-Fourier-Physical format
  phi : Laplacian of v in Fourier-Fourier-Physical format

  vor and phi are real*4, others are real*8.

  Dimensions:
  re100_32 (icase=161): (mgalx,my,mgalz)=(1152,512 ,384 )
  re100_34 (icase=169): (mgalx,my,mgalz)=(1152,1536,576 )
  re250_32 (icase=168): (mgalx,my,mgalz)=(3072,2048,1536)
---------------------------------------------------------------------------

//////////////////////////////////////////////////////////////// 
////////////////// {U,V,W,...}_Snapshots /////////////////////// 
//////////////////////////////////////////////////////////////// 
              Stored as XZY: Physical Fields
//////////////////////////////////////////////////////////////// 

The files are written in standard FORTRAN-UNFORMATTED format in
LITTLE endian instead of BIG as for channels and boundary layer.

After specifying dimensions {mgalx,my,mgalz} (see above), the file
can be read as following:

    open(1,file=filename,access='direct',form='unformatted',recl=mgalx*mgalz)
    read(1,rec=1) mgalx,my,mgalz,numerop
    read(1,rec=2) t,Re,alp,bet,gam,S,chi
    read(1,rec=3) y
    read(1,rec=4) xoffb,xofft,zoffb,zofft

    do j=1,my
       read(1,rec=j+4) ((phy(i,k,j),i=1,mgalx),k=1,mgalz)
    end do
    close(1)

where bet=Ly/pi.

-----------------------------------------------------------------------
                     CONTACT:
-----------------------------------------------------------------------
Please report any problems with the data, queries or doubts to
Siwei Dong, dongsiwei@torroja.dmt.upm.es