=========================================================================================================
Wall-resolved LES database of a separating-reattaching turbulent boundary layer at Re_theta \approx 2E3, M = 0.2

Originally uploaded on March 20, 2023

 Copyright (c) Soshi Kawai. All Rights Reserved.
 Contact: kawai@tohoku.ac.jp
=========================================================================================================
Original paper:
  Ryo Kamogawa, Yoshiharu Tamaki, and Soshi Kawai, 
  "Ordinary-differntial-equation-based nonequilibrium wall modeling for large-eddy simulation," 
  Physical Review Fluids, 8, 064605, (2023).
  https://doi.org/10.2514/1.J062593

Note:
  <>: Reynolds averaged quantities
  {}: Favre-averaged quantities
  Since the Mach number is low, the favre-averaged quantities are essentially equal to the Reynolds-averaged quantities. 
=========================================================================================================
This database includes

# surface data.txt
 x/theta_in
 C_p
 C_f
 H (shape factor)
 Re_theta (Re based on the momentum thickness)
 Re_delta_star (Re based on the displacement thickness)


# profile_{x/c coordinate}.txt 
 Note: the quantities using u_tau is not available at the separated locations (x/theta_in=150,200,250)

 y/theta_in (wall normal coordinate)
 <rho> (density, normalized by freestream density)
 {u} (mean wall-parallel velocity, normalized by freestream velocity)
 {v} (mean wall-normal velocity, normalized by freestream velocity)
 <p> (pressure, normalized by freestream pressure)
 <rho>*{u''u''} (Reynolds normal stress in the streamwise direction, normalized by rho_inf*u_inf^2)
 <rho>*{v''v''} (Reynolds normal stress in wall-normal direction, normalized by rho_inf*u_inf^2)
 <rho>*{w''w''} (Reynolds normal stress in the spanwise direction, normalized by rho_inf*u_inf^2)
 -<rho>*{u''v''} (Reynolds shear stress, normalized by u_tau^2)
 y+ = <rho>y*u_tau/<mu>
 u+ = {U}/u_tau 
 <rho>*{u''u''} (Reynolds normal stress in the streamwise direction, normalized by tau_wall)
 <rho>*{v''v''} (Reynolds normal stress in wall-normal direction, normalized by tau_wall)
 <rho>*{w''w''} (Reynolds normal stress in the spanwise direction, normalized by tau_wall)
 -<rho>*{u''v''} (Reynolds shear stress, normalized by u_tau^2)


# momentum_budget_{x/c coordinate}.txt 
 momentum budget in the inner layer normalized by rho_inf*u_inf^2/theta_in (will appear in the reference paper as Eq. (4))
 y/theta_in (wall normal coordinate)
 Conv. (Convection)
 Pres. (Pressure gradient)
 Re-stress (Reynolds stress)
 Vis-stress (Viscous stress)
 Residual (Sum of above four terms)
 

# 2D data (C binary, little-endian, single precision)
- grid_2d.xyz (normalized by theta_in)
# (imax, jmax) = 2419 x 400
  read(iu) imax,jmax
  read(iu) ((x(i,j),i=1,imax),j=1,jmax),((y(i,j),i=1,imax),j=1,jmax)

- statistics_2d.fun 
# (imax, jmax) = 2419 x 400
# (nvar)       = 9
  read(iu) imax,jmax,nvar
  read(iu) (((data(i,j,n),i=1,imax),j=1,jmax),n=1,nvar)

  data(i,j,1): Mean density, <rho>/rho_inf
  data(i,j,2): Mean x-velocity, {u}/u_inf
  data(i,j,3): Mean y-velocity, {v}/u_inf
  data(i,j,4): Mean pressure, <p>/p_inf
  data(i,j,5): Reynolds normal stress in x direction, <rho>*{u''u''}/(rho_inf*u_inf^2)
  data(i,j,6): Reynolds normal stress in y direction, <rho>*{v''v''}/(rho_inf*u_inf^2)
  data(i,j,7): Reynolds normal stress in z direction, <rho>*{w''w''}/(rho_inf*u_inf^2)
  data(i,j,8): Reynolds shear stress, -<rho>{u''v''}/(rho_inf*u_inf^2)
  data(i,j,9): RMS pressure fluctuation, sqrt(<p'p'>)/(rho_inf*u_inf^2)
  (u, v, w) are the velocities in the Cartesian coordinates (x, y, z) (w being the spanwise velocity)

=========================================================================================================