DNS 1-3 Statistical Data: Difference between revisions
Line 7: | Line 7: | ||
= Statistical data = | = Statistical data = | ||
==Volume data== | ==Volume data== | ||
Volumetric data on the statistics computed | Volumetric data on the statistics computed for the whole computational domain are provided in HDF5 format. This can be easily read through the [https://www.hdfgroup.org/solutions/hdf5/ HDF5 library] or python's [https://www.h5py.org/ h5py]. The tree scheme of the dataset is the following: | ||
* Statistics.h5 | * Statistics.h5 | ||
** 01_Info | ** 01_Info | ||
Line 103: | Line 103: | ||
# <math>{\overline{www}}</math> | # <math>{\overline{www}}</math> | ||
</div> | </div> | ||
[http://kbwiki-data.s3-eu-west-2.amazonaws.com/DNS-1/3/PressureVelocity.h5 Pressure Velocity (2.8 GB)] | |||
Revision as of 13:37, 15 February 2021
Statistical data
Volume data
Volumetric data on the statistics computed for the whole computational domain are provided in HDF5 format. This can be easily read through the HDF5 library or python's h5py. The tree scheme of the dataset is the following:
- Statistics.h5
- 01_Info
- Dimensions
- 02_Entries
- Inputs
- 01_Output
- AdditionalQuantities
- Convection
- Production
- TurbulentDiffusion01
- TurbulentDiffusion02
- PressureStrain
- Dissipation
- TripleCorrelation
- PressureVelocity
- 03_Nodes
- Nodes
- 01_Info
The available files are:
Statistics (9.1 kB) Is the h5 masterfile that contains the data tree and relates with the other files.
Nodes (2.8 GB) Contains the node positions where data is stored. An example of how to read this data in python would be:
import h5py
f = h5py.File('Statistics.h5','r')
xyz = np.array( f.get('03_Nodes').get('Nodes') )
f.close()
Inputs (26.0 GB) Contains the averaged pressure and velocity along with their gradients, the shear stress and the Reynolds stress tensor in the following order:
An example to read this file in python and recover the gradients and Reynolds stress tensor (indices are these of the list above) would be:
inp = np.array( f.get('02_Entries').get('Inputs') )
grad_velocity = inp[:,[17,21,25,18,22,26,19,23,27]].astype(np.double)
Rij = inp[:,[10,11,13,11,12,14,13,14,15]].astype(np.double)
AdditionalQuantities (3.7 GB) Contains additional quantities such as the Taylor and Kolmogorov scales and the pressure autocorrelation in the following order:
- Taylor microscale
- Kolmogorov length scale
- Kolmorogov time scale
TripleCorrelation (9.3 GB) Contains the components of the velocity triple correlation in the following order:
Finally, the components of the Reynolds stress budget equation come in the following order (for a generic budget component ):
Convection (5.6 GB) Contains the convection budget.
Production (5.6 GB) Contains the production budget.
Dissipation (5.6 GB) Contains the dissipation budget.
Molecular Diffusion (5.6 GB) Contains the molecular diffusion budget.
Pressure Strain (5.6 GB) Contains the pressure strain budget.
Turbulent Diffusion 01 (5.6 GB) Contains the turbulent diffusion 1 budget.
Turbulent Diffusion 02 (5.6 GB) Contains the turbulent diffusion 2 budget.
Additional data
No additional data yet.
Contributed by: Oriol Lehmkuhl, Arnau Miro — Barcelona Supercomputing Center (BSC)
© copyright ERCOFTAC 2024