The data for the elliptic streamline flow simulations is on my web site: https://engineering.purdue.edu/AAE/Research/Groups/blaisdell/Elliptic%20Flows%20and%20Strained%20Vortices You will also find there copies of my 1994 and 1996 reports from the CTR summer workshops. The data files for the 1994 cases and the 1996 cases are somewhat different. For example, note that time is scaled differently in the two sets of data files (see below). The cases in the 1994 report in Table 1 on page 359 are labeled A, B, C, D and D'. The actual case names used in the data files are case7 (A), case 11 (B), case12 (C), case5 (D) and case10 (D'). The data for these cases are in the files case7.dat.tar, etc. Each tar file contains separate files for each statistical quantity. The files included are the following: d11.dat dissipation tensor d13.dat d22.dat d33.dat diss.dat dissipation rate iib.dat second invariant II of the Reynolds stress anisotropy tensor iid.dat second invariant II of the dissipation anisotropy tensor iiib.dat third invariant III of the Reynolds stress anisotropy tensor iiid.dat third invariant III of the dissipation anisotropy tensor iiiv.dat third invariant III of Vij defined by Bill Reynolds and Stavros Kassinos iiiy.dat third invariant III of Yij structure dimensionality tensor iiv.dat second invariant II of Vij iiy.dat second invariant II of Yij lam11.dat Taylor microscale, lambdaij lam12.dat lam13.dat lam21.dat lam22.dat lam23.dat lam31.dat lam32.dat lam33.dat lint11.dat Integral length scale, Lambdaij lint12.dat lint13.dat lint21.dat lint22.dat lint23.dat lint31.dat lint32.dat lint33.dat o11.dat Rotation term in the Reynolds stress equation o13.dat o22.dat o33.dat p11.dat Production term in the Reynolds stress equation p13.dat p22.dat p33.dat r11.dat Reynolds stress r11rhs.dat Sum of terms on the RHS of the Reynolds stress equation r13.dat r13rhs.dat r22.dat r22rhs.dat r31.dat r33.dat r33rhs.dat ret.dat Turbulent Reynolds number rotauomega.dat I don't remember rotauvort.dat I don't remember (Rossby number of sorts) skew.dat generalized skewness, defined in eqn. (5) on page 360 t11.dat Pressure strain term (t11 = tf11 + ts11) t13.dat t22.dat t33.dat tf11.dat Fast pressure strain tf13.dat tf22.dat tf33.dat tke.dat TKE tkeprod.dat TKE production term tkerhs.dat sum of the terms on the RHS of the TKE equation ts11.dat Slow pressure strain ts13.dat ts22.dat ts33.dat vratio.dat vorticity ratio, sqrt{enstrophy}/mean vorticity You can compute the Reynolds stress anisotropy tensor, bij, from the given values of the Reynolds stresses rij. The data are in ascii files. There are two columns. The first column contains time scaled by the mean flow period, t/T, where T = 2*pi/Omega and Omega = sqrt{gamma^2 - e^2} (see description below). The second column contains the value of the statistic. Some additional information about the runs is given in the table below (see the 1996 case description further down in this file for the meaning of the parameters): CTR'94 case : A B C D D' Actual case : 7 11 12 5 10 E : 1.1 1.1 1.1 2.0 2.0 \nu : 0.003 0.003 0.003 0.003 0.0015 \alpha_1 : -0.851468 -5.03297 -335.985 -9.15086 -17.9376 \alpha_2 : 0.703693 4.15948 277.674 2.2877 4.48440 \gamma : 0.7775805 4.596225 306.8295 5.71928 11.21100 e : 0.0738875 0.436745 29.1555 3.43158 6.72660 \beta : 0.095 0.095 0.095 0.6 0.6 ************************************************************************** The cases from the 1996 CTR report are labeled e2, e3, e4, and e5 in Table 1 on page 439. The actual case names used in the data files are 25b (e2), 22c (e3), 23h (e4), and 24a (e5). The statistics for these cases that I have extracted are the Reynolds stress anisotrpy, bij; the tke; and the dissipation rate, diss. These are all in one tar file ctr1996stats.tar. If additional statistics are needed, I can get them, but I will have to extract them from a larger data file. The data are in ascii files. There are two columns. The first column contains time scaled by the strain rate (e*t). The second column contains the value of the statistic. Some additional information about the runs is given in the table below: CTR'96 case : e2 e3 e4 e5 Actual case : 25b 22c 23h 24a E : 1.25 1.5 2.0 3.0 \nu : 0.001125 0.001125 0.001125 0.001125 \alpha_1 : -27.7982 -18.1149 -14.3170 -12.8494 \alpha_2 : 17.7908 8.05107 3.57924 1.42771 \gamma : 22.7945 13.0830 8.94811 7.13856 e : 5.00367 5.03193 5.36887 5.71085 \beta : 0.219512 0.384615 0.6 0.8 The mean velocity gradient tensor is of the form: - - | 0 0 \alpha_1 | | | {\partial U_i\over\partial x_j} = | 0 0 0 | | | | \alpha_2 0 0 | - - \alhpa_1 = - \gamma - e \alpha_2 = \gamma - e e = strain rate (negative) used to nondimensionalize time (for 1996 data) E = \sqrt{ -\alpha_1 / \alpha_2 } = Aspect ratio of the elliptic streamlines \beta = e/\gamma = (E^2 - 1)/(E^2 + 1) Let me know if you need any further information. Greg Blaisdell October 11, 2010