=========================================================================================================
Wall-resolved LES database of A-Airfoil near stall condition at 13.3 deg, Re_c = 2.1E6, M = 0.15

Originally uploaded on December  20, 2019
Updated             on September 24, 2020

 Copyright (c) Soshi Kawai. All Rights Reserved.
 Contact: kawai@tohoku.ac.jp
=========================================================================================================

Original paper:
  Kengo Asada & Soshi Kawai, 
  "Large-eddy simulation of airfoil flow near stall condition at Reynolds number 2.1 \times 10^6,"
  Physics of Fluids, 30 (8), 085103, (2018).

Related paper (Wall-modeled LES at the same condition):
  Yoshiharu Tamaki, Yuma Fukushima, Yuichi Kuya & Soshi Kawai,
  "Physics and modeling of trailing-edge stall phenomena for wall-modeled large-eddy simulation"
  Physical Review Fluids (5), 074602, 2020.
  (Note: minus of -<U'V'> in Fig 10 (b) is missing)

=========================================================================================================
This database is for the grid CM3 in the paper and includes

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

# profiles_{x/c coordinate}.txt
 Y/c (wall normal coordinate)
 <U> (mean wall-parallel velocity, normalized by freestrem velocity)
 <U'U'> (Reynolds normal stress in the streamwise direction, normalized by u_inf^2)
 <V'V'> (Reynolds normal stress in the wall-normal direction, normalized by u_inf^2)
 <W'W'> (Reynolds normal stress in the spanwise direction, normalized by u_inf^2)
-<U'V'> (Reynolds shear stress, normalized by u_tau^2)
 y+
 U+= <U>/u_tau
 <U'U'> (Reynolds normal stress in the streamwise direction, normalized by u_tau^2)
 <V'V'> (Reynolds normal stress in the wall-normal direction, normalized by u_tau^2)
 <W'W'> (Reynolds normal stress in the spanwise direction, normalized by u_tau^2)
-<U'V'> (Reynolds shear stress, normalized by u_tau^2)

# budget_{x/c coordinate}.txt (momentum budget, see Fig. 20 of the original paper)
 Y/c
 convection        (normalized by u_inf^2/c)
 viscous           (normalized by u_inf^2/c)
 pressure_gradient (normalized by u_inf^2/c)
 Reynolds_stress   (normalized by u_inf^2/c)

# 2D data (C binary, little-endian, single precision, near-field cropped)
- grid_2d.xyz
# (imax, jmax) = 4941 x 288
  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) = 4941 x 288
# (nvar)       = 8
  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, (u'u')/(u_inf^2)
  data(i,j,6): Reynolds normal stress in y direction, (v'v')/(u_inf^2)
  data(i,j,7): Reynolds normal stress in z direction, (w'w')/(u_inf^2)
  data(i,j,8): Reynolds shear stress, (-u'v')/(u_inf^2)
  (u, v, w) are the velocities in the Cartesian coordinates (x, y, z) (w being the spanwise velocity)
=========================================================================================================