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 [2]:
F = h5.File('../data/P0550.h5','r') P0550_D_uu_inner =F['D_uu_inner'][()] P0550_D_uu_outer =F['D_uu_outer'][()] P0550_E_uu_inner =F['E_uu_inner'][()] P0550_E_uu_outer =F['E_uu_outer'][()] P0550_P_uu_inner =F['P_uu_inner'][()] P0550_P_uu_outer =F['P_uu_outer'][()] P0550_PI1_uu_inner =F['PI1_uu_inner'][()] P0550_PI1_uu_outer =F['PI1_uu_outer'][()] P0550_T_uu_xz_inner =F['T_uu_xz_inner'][()] P0550_T_uu_xz_outer =F['T_uu_xz_outer'][()] P0550_T_uu_y_inner =F['T_uu_y_inner'][()] P0550_T_uu_y_outer =F['T_uu_y_outer'][()] P0550_yp_loc =F['yp_loc'][()] F.close()
In [3]:
F = h5.File('../data/P1000.h5','r') P1000_D_uu_inner =F['D_uu_inner'][()] P1000_D_uu_outer =F['D_uu_outer'][()] P1000_E_uu_inner =F['E_uu_inner'][()] P1000_E_uu_outer =F['E_uu_outer'][()] P1000_P_uu_inner =F['P_uu_inner'][()] P1000_P_uu_outer =F['P_uu_outer'][()] P1000_PI1_uu_inner =F['PI1_uu_inner'][()] P1000_PI1_uu_outer =F['PI1_uu_outer'][()] P1000_T_uu_xz_inner =F['T_uu_xz_inner'][()] P1000_T_uu_xz_outer =F['T_uu_xz_outer'][()] P1000_T_uu_y_inner =F['T_uu_y_inner'][()] P1000_T_uu_y_outer =F['T_uu_y_outer'][()] P1000_yp_loc =F['yp_loc'][()] F.close()
In [4]:
F = h5.File('../data/P2000.h5','r') P2000_D_uu_inner =F['D_uu_inner'][()] P2000_D_uu_outer =F['D_uu_outer'][()] P2000_E_uu_inner =F['E_uu_inner'][()] P2000_E_uu_outer =F['E_uu_outer'][()] P2000_P_uu_inner =F['P_uu_inner'][()] P2000_P_uu_outer =F['P_uu_outer'][()] P2000_PI1_uu_inner =F['PI1_uu_inner'][()] P2000_PI1_uu_outer =F['PI1_uu_outer'][()] P2000_T_uu_xz_inner =F['T_uu_xz_inner'][()] P2000_T_uu_xz_outer =F['T_uu_xz_outer'][()] P2000_T_uu_y_inner =F['T_uu_y_inner'][()] P2000_T_uu_y_outer =F['T_uu_y_outer'][()] P2000_yp_loc =F['yp_loc'][()] F.close()
In [5]:
F = h5.File('../data/P5200.h5','r') P5200_D_uu_inner =F['D_uu_inner'][()] P5200_D_uu_outer =F['D_uu_outer'][()] P5200_E_uu_inner =F['E_uu_inner'][()] P5200_E_uu_outer =F['E_uu_outer'][()] P5200_P_uu_inner =F['P_uu_inner'][()] P5200_P_uu_outer =F['P_uu_outer'][()] P5200_PI1_uu_inner =F['PI1_uu_inner'][()] P5200_PI1_uu_outer =F['PI1_uu_outer'][()] P5200_T_uu_xz_inner =F['T_uu_xz_inner'][()] P5200_T_uu_xz_outer =F['T_uu_xz_outer'][()] P5200_T_uu_y_inner =F['T_uu_y_inner'][()] P5200_T_uu_y_outer =F['T_uu_y_outer'][()] P5200_yp_loc =F['yp_loc'][()] F.close()
In [6]:
plt.figure(figsize=(10,16)) plt.subplot(521) plt.plot(P0550_yp_loc, P0550_P_uu_inner, 'C0-') plt.plot(P1000_yp_loc, P1000_P_uu_inner, 'C1-') plt.plot(P2000_yp_loc, P2000_P_uu_inner, 'C2-') plt.plot(P5200_yp_loc, P5200_P_uu_inner, 'C3-') plt.ylabel(r'$P^+_{\lambda^+\ge 1000}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.10,0.10]) plt.text(-4,0.1,'(a)') plt.subplot(522) plt.plot(P0550_yp_loc, P0550_P_uu_outer, 'C0-', label=r'$Re_\tau = 550$') plt.plot(P1000_yp_loc, P1000_P_uu_outer, 'C1-', label=r'$Re_\tau = 1000$') plt.plot(P2000_yp_loc, P2000_P_uu_outer, 'C2-', label=r'$Re_\tau = 2000$') plt.plot(P5200_yp_loc, P5200_P_uu_outer, 'C3-', label=r'$Re_\tau = 5200$') plt.ylabel(r'$P^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.legend(frameon=False,ncol=1,loc='lower right') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.text(-4,0.1,'(b)') plt.subplot(525) plt.plot(P0550_yp_loc, P0550_T_uu_y_inner + P0550_T_uu_xz_inner, 'C0-') plt.plot(P1000_yp_loc, P1000_T_uu_y_inner + P1000_T_uu_xz_inner, 'C1-') plt.plot(P2000_yp_loc, P2000_T_uu_y_inner + P2000_T_uu_xz_inner, 'C2-') plt.plot(P5200_yp_loc, P5200_T_uu_y_inner + P5200_T_uu_xz_inner, 'C3-') plt.ylabel(r'$T^+_{\lambda^+ \ge 1000}$') plt.xlabel(r'$y^+$') plt.text(-4,0.1,'(e)') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.subplot(526) plt.plot(P0550_yp_loc, P0550_T_uu_y_outer + P0550_T_uu_xz_outer, 'C0-') plt.plot(P1000_yp_loc, P1000_T_uu_y_outer + P1000_T_uu_xz_outer, 'C1-') plt.plot(P2000_yp_loc, P2000_T_uu_y_outer + P2000_T_uu_xz_outer, 'C2-') plt.plot(P5200_yp_loc, P5200_T_uu_y_outer + P5200_T_uu_xz_outer, 'C3-') plt.ylabel(r'$T^+_{ \lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.text(-4,0.1,'(f)') plt.subplot(523) plt.plot(P0550_yp_loc, P0550_PI1_uu_inner, 'C0-') plt.plot(P1000_yp_loc, P1000_PI1_uu_inner, 'C1-') plt.plot(P2000_yp_loc, P2000_PI1_uu_inner, 'C2-') plt.plot(P5200_yp_loc, P5200_PI1_uu_inner, 'C3-') plt.text(-4,0.1,'(c)') plt.ylabel(r'$\Pi^+_{\lambda^+ \ge 1000}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.subplot(524) plt.plot(P0550_yp_loc, P0550_PI1_uu_outer, 'C0-') plt.plot(P1000_yp_loc, P1000_PI1_uu_outer, 'C1-') plt.plot(P2000_yp_loc, P2000_PI1_uu_outer, 'C2-') plt.plot(P5200_yp_loc, P5200_PI1_uu_outer, 'C3-') plt.ylabel(r'$\Pi^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.text(-4,0.1,'(d)') plt.subplot(527) plt.plot(P0550_yp_loc, P0550_D_uu_inner, 'C0-') plt.plot(P1000_yp_loc, P1000_D_uu_inner, 'C1-') plt.plot(P2000_yp_loc, P2000_D_uu_inner, 'C2-') plt.plot(P5200_yp_loc, P5200_D_uu_inner, 'C3-') plt.ylabel(r'$D^+_{\lambda^+ \ge 1000}$') plt.xlabel(r'$y^+$') plt.text(-4,0.1,'(g)') plt.xlim([0,20]) plt.ylim([-0.1,0.10]) plt.subplot(528) plt.plot(P0550_yp_loc, P0550_D_uu_outer, 'C0-') plt.plot(P1000_yp_loc, P1000_D_uu_outer, 'C1-') plt.plot(P2000_yp_loc, P2000_D_uu_outer, 'C2-') plt.plot(P5200_yp_loc, P5200_D_uu_outer, 'C3-') plt.ylabel(r'$D^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.text(-4,0.1,'(h)') plt.subplot(529) plt.plot(P0550_yp_loc, -P0550_E_uu_inner, 'C0-') plt.plot(P1000_yp_loc, -P1000_E_uu_inner, 'C1-') plt.plot(P2000_yp_loc, -P2000_E_uu_inner, 'C2-') plt.plot(P5200_yp_loc, -P5200_E_uu_inner, 'C3-') plt.ylabel(r'$-\epsilon^+_{\lambda^+ \ge 1000}$') plt.xlabel(r'$y^+$') plt.text(-4,0.1,'(i)') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.subplot(5,2,10) plt.plot(P0550_yp_loc, -P0550_E_uu_outer, 'C0-') plt.plot(P1000_yp_loc, -P1000_E_uu_outer, 'C1-') plt.plot(P2000_yp_loc, -P2000_E_uu_outer, 'C2-') plt.plot(P5200_yp_loc, -P5200_E_uu_outer, 'C3-') plt.ylabel(r'$-\epsilon^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.1,0.1]) plt.text(-4,0.1,'(j)') plt.subplots_adjust(hspace=0.4,wspace=0.3) plt.savefig('figure05.eps',dpi=300,bbox_inches='tight') plt.show() plt.close()
Out[6]:
