Index of /turbdata/raw_database/Boundary_Layer/High-Reynolds
========================================================================
= STATISTICS FROM DNS OF ZERO PRESSURE GRADIENT =
= BOUNDARY LAYER =
= IN A VERY LARGE NUMERICAL BOX =
= =
= Re_theta = 2780-6650 (Production Simulations) =
= Re_theta = 1100-3000 (Auxiliary Simulations) =
========================================================================
This directory tree contain some of the restart files of a direct
numerical simulations of fully developed ZPG Boundary Layer flow
together with a sample matlab subroutine to read the files. Details
are offered below.
Folders: HR2-Reth-2800-6650: Main production Simulation. High Resolution.
LR1-Reth-1100-3000: Auxiliary Simulation. Lower Resolution.
HR2-d99_chunk : Boxes of d99(Re_theta) size of the Production
Simulation at several Reynolds numbers.
-----------------------------------------------------------------------
GENERAL INFORMATION:
-----------------------------------------------------------------------
Authors: Juan A. Sillero (sillero@torroja.dmt.upm.es)
Javier Jimenez (jimenez@torroja.dmt.upm.es)
Formulation: Primitive variables.
////////CASE: Production Simulation "HR2-Reth-2800-6650"
Reynolds number: Re_theta=2780-6650
Domain size: (Lx,Ly,Lz)/theta=(547,29,84)
Resolution: (dx+,dy_max+,dz+)=(6.5,10.2,3.8)
*Taken at the middle of the box.
Number of fields: 212 {u,v,w,p}. Single precission.
Number of files & size: 848. Each file of size 84 GB
Total Size: 70 TB.
Name: tbl2-HR-Prod.xxx.var.h5 (var=u,v,w,p)
Format: Binary HDF5 (stored as ZYX)
////////CASE: Auxiliary Simulation "LR1-Reth-1100-3000"
Reynolds number: Re_theta=1100-3000
Domain size: (Lx,Ly,Lz)/theta=(481,47,191)
Resolution: (dx+,dy_max+,dz+)=(13,10.2,7.3)
* Taken at the middle of the box.
Number of fields: 212 {u,v,w,p}. Single precission.
Number of files & size: 848. Each file of size 7.2 GB
Total Size: 6 TB.
Name: tbl1-HR-Prod.xxx.var.h5 (var=u,v,w,p)
Format: Binary HDF5 (stored as ZYX)
-----------------------------------------------------------------------
STRUCTURE OF THE FILES:
-----------------------------------------------------------------------
The files are written in HDF5 format, which 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:
(kz,ny,nx) ---> ( Fourier - Physical -Physical )
For instance, one file contains:
Re Dataset {1} : Inverse of the viscosity used in the code,
which defines the scaling in which the code works.
Variable Dataset {SCALAR} : Variable recorded "u", "v", "w" or "p".
procs Dataset {536} : Number of nodes used to run the simulation.
(first element of the vector!)
cfl Dataset {1} : Courant–Friedrichs–Lewy condition at which the
simulation was run.
dt Dataset {1} : Time step in code units for this flow realization.
tiempo Dataset {1} : Time stamp in code units for this flow realization.
timeinit Dataset {1} : Initial time stamp in code units (=0).
lx Dataset {1} : Defines the streamwise length Lx of the box, in code units.
ly Dataset {1} : Defines the wall-normal length Ly of the box, in code units.
lz Dataset {1} : Defines the spanwise length Lz of the box, in code units.
nx Dataset {1} : Number of collocation points in the streamwise direction.
ny Dataset {1} : Number of collocation points in the wall-normal direction.
nz2 Dataset {1} : Half the number of complex Fourier Gallerkin modes in z-> K_z[0:nz2]
We use the 2/3 dealiasing rule in z, so the number of collocation points in those directions are 3/2
times the number of Fourier Gallerkin modes, and the extra modes are padded with zeroes.
xout Dataset {1} : Location of the recycling plane.
um Dataset {536} : Mean streamwise velocity profile in code units at the recycling
station "xout".
y Dataset {537} : Non-uniform grid for the wall-normal distance; in code units.
Note that, because of the staggered grid: y=(-y1,0,y1,y2,....)
value Dataset {15361, 536, 2730} : Primitive variable in code units. The array is stored in single-precision
as [1:2*(nz2+1),1:Ny,1:Nx], equivalent to the complex array [Kz,Ny,Nx]
-----------------------------------------------------------------------
Utilities:
-----------------------------------------------------------------------
[1] The sketch of the grid lay-out is specified in the file
"sketch_grid_DNS.pdf" (If possible, use Adobe Reader).
[2] A sample matlab subroutine "convert_2_phys.m" to read and
transform a portion of a given field from Fourier to Physical
space.
-----------------------------------------------------------------------
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] G. Borrell, J.A. Sillero, J. Jimenez
"A code for direct numerical simulation of turbulent boundary
layers at high Reynolds numbers in BG/P supercomputers"
Computers & Fluids, Volume 80, 10 July 2013, Pages 37-43
[3] J.A. Sillero, J. Jimenez, R.D. Moser
"One-point statistics for turbulent wall-bounded flows at
Reynolds numbers up to \delta^+\approx2000"
Phys. Fluids 25, 105102, 3 October 2013
[4] J.A. Sillero, J. Jimenez, R.D. Moser
"Two-point statistics for turbulent boundary layers and channel flows
at Reynolds numbers up to \delta^+\approx2000"
Phys. Fluids 26, 105109, 28 October 2014
-----------------------------------------------------------------------
FUNDING:
-----------------------------------------------------------------------
This research was funded in part by the Spanish CICYT under Grant
No. TRA2009-11498, and by the European Research Council under Grant
No. ERC-2010.AdG-20100224. It used the computational resources of the
Argonne Leadership Computing Facility at Argonne National Laboratory,
which is supported by the Office of Science of the U.S. Department of
Energy under Contract No. DE-AC02-06CH11357. J. A. Sillero was
supported in part by an FPU fellowship from the Universidad
Politecnica de Madrid.
-----------------------------------------------------------------------
CONTACT:
-----------------------------------------------------------------------
Please report any problems with the data, queries or doubts to
J.A. Sillero, sillero@torroja.dmt.upm.es