Datenpaket: Dataset for "A framework for assessing the Reynolds analogy in turbulent forced convection over rough walls"

Direct numerical simulations (DNS) are used to assess the similarity, in rough-wall turbulent flow, between the time-averaged transport of mechanical energy (i.e. the sum of kienetic energy and pressure work) and of the squared passive-scalar field, theta^2/2. The database comprises three simulations:

1. k00re180, smooth wall at a friction Reynolds number Re_tau=180;
2. k15re180, rough wall at a friction Reynolds number Re_tau=180 and roughness Reynolds number k^+=15 (k^+=k u_tau/nu, k is the roughness characteristics size, u_tau is the friction velocity, and nu is the kinematic viscosity);
3. k90re540, rough-wall at a friction Reynolds number Re_tau=540 and roughness Reynolds number k^+=90. The roughness topography is taken from: 'Thakkar, M. (2017) Investigation of turbulent flow over irregular rough surfaces using direct numerical simulations (PhD Thesis), University of Southampton, Southampton, UK'. The DNS are performed using the spectral element flow solver Nek5000 (Version 19.0. 2019, Argonne National Laboratory, Illinois. https://nek5000.mcs.anl.gov). The dataset contains the following files:
4. Mean wall-normal profiles (ascii),
   • k00re180.dat
   • k15re180.dat
   • k90re540.dat
5. Volume-averaged quantities (i.e. volume integrated quantities per unit volume) (ascii)
   • BULK_k00re180.dat
   • BULK_k15re180.dat
   • BULK_k90re540.dat
6. Time-averaged dispersive fluctuations distributions on the plane y/delta=1.2 (y is the spanwise coordinate and delta is the half channel height) (binary)
   • sfA_y120_k15re180.vtk
   • sfA_y120_k90re540.vtk In each file, the variables are made dimensionless using the friction velocity u_tau and the channel half-height delta. For all the cases, a unit Prandtl number is used. Therefore, the friction Peclet number is equal to the friction Reynolds number, Pe_tau=Re_tau. Furthermore, the flow is always considered incompressible. In the following, Ui indicates the i-velocity component; C indicates the scalar field; P indicates the pressure; K=UiUi/2 is the kinetic energy; G=C^2/2 is the squared scalar field; ,i denotes differentiation in the i-direction; Oi indicates the i-vorticity component; Sij indicates the (i,j)-component of the strain-rate tensor; < . >t indicates a time-averaged quantity; < . >ts indicates a space- and time-averaged quantity; a prime superscript, (.)', indicates temporal fluctuations; a double prime superscript, (.)'', indicates dispersive fluctuations of the time-averaged field induced by the roughness. Note that the coordinate directions are: Streamwise, x, i=1; Spanwise, y, i=2; Wall-normal, z, i=3. The dataset files contain:
7. kXXreXXX.dat zts<C,i>ts<C,i>ts/Pe_tau<C,i'C,i'>/Pe_tautsts/Re_tau<Oi'Oi'>/Re_tau 2tsts/Re_tau 2<Sij'Sij'>ts/Re_tau<(Ui P),i>ts<(Ui K),i>ts<P,ii>ts/Re_tau<K,ii>ts/Re_tauts<(Ui G),i>ts<G,ii>ts/Pe_tau
8. BULK_kXXreXXX.dattt<Oi'Oi'>tt''t''tststtstst''<C,i>t<C,i>t<C,i'C,i'>t<C,i>t''<C,i>t''<C,i>ts<C,i>ts<C,i>t<C,i>ts<C,i>ts<C,i>t'' AER[tt] AER[<Oi'Oi'>t] AER[t''t''] AER[tsts] AERtts AER[tst'' AER[<C,i>t<C,i>t] AER[<C,i'C,i'>t] AER[<C,i>t''<C,i>t''] AER[<C,i>ts<C,i>ts] AER[<C,i>t<C,i>ts] AER[<C,i>ts<C,i>t''] WMR[tt] WMR[<Oi'Oi'>t] WMR[t''t''] WMR[tsts ] WMR[tts] WMR[tst''] WMR[<C,i>t<C,i>t] WMR[<C,i'C,i'>t] WMR[<C,i>t''<C,i>t''] WMR[<C,i>ts<C,i>ts] WMR[<C,i>t<C,i>ts] WMR[<C,i>ts<C,i>t''] Note: AER[ ] and WMR[ ] indicate conditional integration over "attached, exposed" and "well-mixed" region. Further information is provided in the header of each file.
9. sfA_y120_kXXreXXX.vtk VTK NAME: corresponding quantity p:

   t s00: t s01: t'' s02: t'' s03: t'' s04:

   t'' s05: t'' s06: t''t'' s07: t''t'' s08: t''t'' s09: <C,x>t''<C,x>t'' s10: <C,y>t''<C,y>t'' s11: <C,z>t''<C,z>t'' s12: t'' s13: t'' s14: t'' s15: <C,x>t'' s16: <C,y>t'' s17: <C,z>t'' s18: t s19: ts s20: <C,z>t s21: <C,z>ts s22: ts s23: ts vx: t vy: t vz: t zwb: wall-height of the bottom rough wall as a function of x.