Kernel: Python 3 (ipykernel)
In [1]:
import numpy as np import h5py as h5 import matplotlib.pyplot as plt from matplotlib import rcParams rcParams['font.family'] = 'serif' rcParams['font.serif'] = ['Computer Modern Roman'] rcParams['text.usetex'] = True rcParams['xtick.major.size']=7 rcParams['ytick.major.size']=7 rcParams['xtick.minor.size']=3.5 rcParams['ytick.minor.size']=3.5 rcParams['xtick.labelsize']=10 rcParams['ytick.labelsize']=10 rcParams['contour.negative_linestyle'] = 'solid'
In [6]:
F = h5.File('../data/P0550.h5','r') P0550_T_uu = F['T_uu'][()] P0550_T_uu_xz_neg = F['T_uu_xz_neg'][()] P0550_T_uu_xz_pos = F['T_uu_xz_pos'][()] P0550_T_uu_y_neg = F['T_uu_y_neg'][()] P0550_T_uu_y_pos = F['T_uu_y_pos'][()] P0550_dUdy = F['dUdy'][()] P0550_nu = F['nu'][()] P0550_ny = F['ny'][()] P0550_y_loc = F['y_loc'][()] F.close() P0550_u_tau = np.sqrt(P0550_dUdy*P0550_nu) P0550_Re_tau = P0550_u_tau / P0550_nu P0550_yp = P0550_y_loc*P0550_Re_tau P0550_scale = 1./(P0550_u_tau**2)/P0550_dUdy P0550_T_uu *= P0550_scale P0550_T_uu_xz_neg *= P0550_scale P0550_T_uu_xz_pos *= P0550_scale P0550_T_uu_y_neg *= P0550_scale P0550_T_uu_y_pos *= P0550_scale
In [7]:
F = h5.File('../data/P1000.h5','r') P1000_T_uu = F['T_uu'][()] P1000_T_uu_xz_neg = F['T_uu_xz_neg'][()] P1000_T_uu_xz_pos = F['T_uu_xz_pos'][()] P1000_T_uu_y_neg = F['T_uu_y_neg'][()] P1000_T_uu_y_pos = F['T_uu_y_pos'][()] P1000_dUdy = F['dUdy'][()] P1000_nu = F['nu'][()] P1000_ny = F['ny'][()] P1000_y_loc = F['y_loc'][()] F.close() P1000_u_tau = np.sqrt(P1000_dUdy*P1000_nu) P1000_Re_tau = P1000_u_tau / P1000_nu P1000_yp = P1000_y_loc*P1000_Re_tau P1000_scale = 1./(P1000_u_tau**2)/P1000_dUdy P1000_T_uu *= P1000_scale P1000_T_uu_xz_neg *= P1000_scale P1000_T_uu_xz_pos *= P1000_scale P1000_T_uu_y_neg *= P1000_scale P1000_T_uu_y_pos *= P1000_scale
In [8]:
F = h5.File('../data/P2000.h5','r') P2000_T_uu = F['T_uu'][()] P2000_T_uu_xz_neg = F['T_uu_xz_neg'][()] P2000_T_uu_xz_pos = F['T_uu_xz_pos'][()] P2000_T_uu_y_neg = F['T_uu_y_neg'][()] P2000_T_uu_y_pos = F['T_uu_y_pos'][()] P2000_dUdy = F['dUdy'][()] P2000_nu = F['nu'][()] P2000_ny = F['ny'][()] P2000_y_loc = F['y_loc'][()] F.close() P2000_u_tau = np.sqrt(P2000_dUdy*P2000_nu) P2000_Re_tau = P2000_u_tau / P2000_nu P2000_yp = P2000_y_loc*P2000_Re_tau P2000_scale = 1./(P2000_u_tau**2)/P2000_dUdy P2000_T_uu *= P2000_scale P2000_T_uu_xz_neg *= P2000_scale P2000_T_uu_xz_pos *= P2000_scale P2000_T_uu_y_neg *= P2000_scale P2000_T_uu_y_pos *= P2000_scale
In [10]:
F = h5.File('../data/P5200.h5','r') P5200_T_uu = F['T_uu'][()] P5200_T_uu_xz_neg = F['T_uu_xz_neg'][()] P5200_T_uu_xz_pos = F['T_uu_xz_pos'][()] P5200_T_uu_y_neg = F['T_uu_y_neg'][()] P5200_T_uu_y_pos = F['T_uu_y_pos'][()] P5200_dUdy = F['dUdy'][()] P5200_nu = F['nu'][()] P5200_ny = F['ny'][()] P5200_y_loc = F['y_loc'][()] F.close() P5200_u_tau = np.sqrt(P5200_dUdy*P5200_nu) P5200_Re_tau = P5200_u_tau / P5200_nu P5200_yp = P5200_y_loc*P5200_Re_tau P5200_scale = 1./(P5200_u_tau**2)/P5200_dUdy P5200_T_uu *= P5200_scale P5200_T_uu_xz_neg *= P5200_scale P5200_T_uu_xz_pos *= P5200_scale P5200_T_uu_y_neg *= P5200_scale P5200_T_uu_y_pos *= P5200_scale
In [11]:
plt.figure(figsize=(10,3.5)) plt.subplot(121) plt.plot(P0550_yp, P0550_T_uu_y_pos, 'C0-' ) plt.plot(P1000_yp, P1000_T_uu_y_pos, 'C1-' ) plt.plot(P2000_yp, P2000_T_uu_y_pos, 'C2-' ) plt.plot(P5200_yp, P5200_T_uu_y_pos, 'C3-' ) plt.plot(P0550_yp, P0550_T_uu_y_neg, 'C0--' ) plt.plot(P1000_yp, P1000_T_uu_y_neg, 'C1--' ) plt.plot(P2000_yp, P2000_T_uu_y_neg, 'C2--' ) plt.plot(P5200_yp, P5200_T_uu_y_neg, 'C3--' ) plt.ylim([-0.2,0.2]) plt.plot([-10,0],[0,0],'k-',label=r'$E_{u^{\prime 2}}^{T^\bot} \ge 0$') plt.plot([-10,0],[0,0],'k--',label=r'$E_{u^{\prime 2}}^{T^\bot} < 0$') plt.legend(frameon=False,loc='upper right') plt.xlabel(r'$y^+$'); plt.xlim([0,20])#; plt.xscale('log') #plt.ylim([0,14]) plt.text(-5,0.2,'(a)') plt.ylabel(r'$T^+_\bot$') #plt.title('Poiseuille') plt.subplot(122) plt.plot(P0550_yp, P0550_T_uu_xz_pos, 'C0-', label=r'$Re_\tau = 550$ ') plt.plot(P1000_yp, P1000_T_uu_xz_pos, 'C1-', label=r'$Re_\tau = 1000$') plt.plot(P2000_yp, P2000_T_uu_xz_pos, 'C2-', label=r'$Re_\tau = 2000$') plt.plot(P5200_yp, P5200_T_uu_xz_pos, 'C3-', label=r'$Re_\tau = 5200$') plt.plot(P0550_yp, P0550_T_uu_xz_neg, 'C0--' ) plt.plot(P1000_yp, P1000_T_uu_xz_neg, 'C1--' ) plt.plot(P2000_yp, P2000_T_uu_xz_neg, 'C2--' ) plt.plot(P5200_yp, P5200_T_uu_xz_neg, 'C3--' ) plt.plot([-10,0],[0,0],'k-',label=r'$E_{u^{\prime 2}}^{T^\|} \ge 0$') plt.plot([-10,0],[0,0],'k--',label=r'$E_{u^{\prime 2}}^{T^\|} < 0$') plt.text(-5,0.2,'(b)') plt.legend(frameon=False,ncol=1,loc='center right') plt.xlabel(r'$y^+$'); plt.xlim([0,20])#; plt.xscale('log') plt.ylabel(r'$T^+_\|$') plt.ylim([-0.2,0.2]) plt.subplots_adjust(hspace=0.4,wspace=0.3) plt.savefig('figure09.eps',dpi=300,bbox_inches='tight') plt.show() plt.close()
Out[11]:
In [ ]:
In [ ]: