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/C0093.h5','r') C0093_T_uu_xz_inner =F['T_uu_xz_inner'][()] C0093_T_uu_xz_outer =F['T_uu_xz_outer'][()] C0093_T_uu_y_inner =F['T_uu_y_inner'][()] C0093_T_uu_y_outer =F['T_uu_y_outer'][()] C0093_yp_loc =F['yp_loc'][()] F.close()
In [3]:
F = h5.File('../data/C0220.h5','r') C0220_T_uu_xz_inner =F['T_uu_xz_inner'][()] C0220_T_uu_xz_outer =F['T_uu_xz_outer'][()] C0220_T_uu_y_inner =F['T_uu_y_inner'][()] C0220_T_uu_y_outer =F['T_uu_y_outer'][()] C0220_yp_loc =F['yp_loc'][()] F.close()
In [4]:
F = h5.File('../data/C0500.h5','r') C0500_T_uu_xz_inner =F['T_uu_xz_inner'][()] C0500_T_uu_xz_outer =F['T_uu_xz_outer'][()] C0500_T_uu_y_inner =F['T_uu_y_inner'][()] C0500_T_uu_y_outer =F['T_uu_y_outer'][()] C0500_yp_loc =F['yp_loc'][()] F.close()
In [5]:
F = h5.File('../data/P0550.h5','r') 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 [6]:
F = h5.File('../data/P1000.h5','r') 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 [7]:
F = h5.File('../data/P2000.h5','r') 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 [8]:
F = h5.File('../data/P5200.h5','r') 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 [9]:
plt.figure(figsize=(8,5)) plt.subplot(221) plt.plot(P0550_yp_loc, P0550_T_uu_y_inner, 'C0-' , label=r'$Re_\tau = 550$ ') plt.plot(P1000_yp_loc, P1000_T_uu_y_inner, 'C1-' , label=r'$Re_\tau = 1000$') plt.plot(P2000_yp_loc, P2000_T_uu_y_inner, 'C2-' , label=r'$Re_\tau = 2000$') plt.plot(P5200_yp_loc, P5200_T_uu_y_inner, 'C3-' , label=r'$Re_\tau = 5200$') plt.plot(P0550_yp_loc, P0550_T_uu_xz_inner, 'C0--' ) plt.plot(P1000_yp_loc, P1000_T_uu_xz_inner, 'C1--' ) plt.plot(P2000_yp_loc, P2000_T_uu_xz_inner, 'C2--' ) plt.plot(P5200_yp_loc, P5200_T_uu_xz_inner, 'C3--' ) plt.legend(frameon=False,loc='upper right') plt.xlabel(r'$y^+$'); plt.xlim([0,20]) plt.ylim([-0.01,0.04]) plt.ylabel(r'$T^+_{\lambda^+ \ge 1000}$') plt.title('Poiseuille') plt.text(-5,0.04,'(a)') plt.subplot(222) plt.plot(C0093_yp_loc, C0093_T_uu_y_inner, 'C4-' , label=r'$Re_\tau = 93$') plt.plot(C0220_yp_loc, C0220_T_uu_y_inner, 'C5-' , label=r'$Re_\tau = 220$') plt.plot(C0500_yp_loc, C0500_T_uu_y_inner, 'C6-' , label=r'$Re_\tau = 500$') plt.plot(C0093_yp_loc, C0093_T_uu_xz_inner, 'C4--' ) plt.plot(C0220_yp_loc, C0220_T_uu_xz_inner, 'C5--' ) plt.plot(C0500_yp_loc, C0500_T_uu_xz_inner, 'C6--' ) plt.text(-5,0.04,'(b)') plt.xlabel(r'$y^+$'); plt.xlim([0,20]) plt.ylim([-0.01,0.04]) plt.ylabel(r'$T^+_{\lambda^+ \ge 1000}$') plt.legend(frameon=False,loc='upper right') plt.title('Couette') plt.subplot(223) plt.plot(P0550_yp_loc, P0550_T_uu_y_outer, 'C0-') plt.plot(P1000_yp_loc, P1000_T_uu_y_outer, 'C1-') plt.plot(P2000_yp_loc, P2000_T_uu_y_outer, 'C2-') plt.plot(P5200_yp_loc, P5200_T_uu_y_outer, 'C3-') plt.plot(P0550_yp_loc, P0550_T_uu_xz_outer, 'C0--') plt.plot(P1000_yp_loc, P1000_T_uu_xz_outer, 'C1--') plt.plot(P2000_yp_loc, P2000_T_uu_xz_outer, 'C2--') plt.plot(P5200_yp_loc, P5200_T_uu_xz_outer, 'C3--') plt.xlabel(r'$y^+$'); plt.xlim([0,20]) plt.ylim([-0.005,0.015]) plt.ylabel(r'$T^+_{\lambda \ge \delta}$') plt.text(-5,0.015,'(c)') plt.subplot(224) plt.plot(C0093_yp_loc, C0093_T_uu_y_outer, 'C4-') plt.plot(C0220_yp_loc, C0220_T_uu_y_outer, 'C5-') plt.plot(C0500_yp_loc, C0500_T_uu_y_outer, 'C6-') plt.plot(C0093_yp_loc, C0093_T_uu_xz_outer, 'C4--') plt.plot(C0220_yp_loc, C0220_T_uu_xz_outer, 'C5--') plt.plot(C0500_yp_loc, C0500_T_uu_xz_outer, 'C6--') plt.xlabel(r'$y^+$'); plt.xlim([0,20]) plt.ylim([-0.08,0.10]) plt.ylabel(r'$T^+_{\lambda \ge \delta}$') plt.text(-5,0.1,'(d)') plt.subplots_adjust(hspace=0.4,wspace=0.3) plt.savefig('figure10.eps',dpi=300,bbox_inches='tight') plt.show() plt.close()
Out[9]:
