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_D_uu_inner =F['D_uu_inner'][()] C0093_D_uu_outer =F['D_uu_outer'][()] C0093_E_uu_inner =F['E_uu_inner'][()] C0093_E_uu_outer =F['E_uu_outer'][()] C0093_P_uu_inner =F['P_uu_inner'][()] C0093_P_uu_outer =F['P_uu_outer'][()] C0093_PI1_uu_inner =F['PI1_uu_inner'][()] C0093_PI1_uu_outer =F['PI1_uu_outer'][()] 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_D_uu_inner =F['D_uu_inner'][()] C0220_D_uu_outer =F['D_uu_outer'][()] C0220_E_uu_inner =F['E_uu_inner'][()] C0220_E_uu_outer =F['E_uu_outer'][()] C0220_P_uu_inner =F['P_uu_inner'][()] C0220_P_uu_outer =F['P_uu_outer'][()] C0220_PI1_uu_inner =F['PI1_uu_inner'][()] C0220_PI1_uu_outer =F['PI1_uu_outer'][()] 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_D_uu_inner =F['D_uu_inner'][()] C0500_D_uu_outer =F['D_uu_outer'][()] C0500_E_uu_inner =F['E_uu_inner'][()] C0500_E_uu_outer =F['E_uu_outer'][()] C0500_P_uu_inner =F['P_uu_inner'][()] C0500_P_uu_outer =F['P_uu_outer'][()] C0500_PI1_uu_inner =F['PI1_uu_inner'][()] C0500_PI1_uu_outer =F['PI1_uu_outer'][()] 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]:
plt.figure(figsize=(10,16)) plt.subplot(521) plt.plot(C0093_yp_loc, C0093_P_uu_inner, 'C4-') plt.plot(C0220_yp_loc, C0220_P_uu_inner, 'C5-') plt.plot(C0500_yp_loc, C0500_P_uu_inner, 'C6-') 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(C0093_yp_loc, C0093_P_uu_outer, 'C4-', label=r'$Re_\tau = 93$') plt.plot(C0220_yp_loc, C0220_P_uu_outer, 'C5-', label=r'$Re_\tau = 220$') plt.plot(C0500_yp_loc, C0500_P_uu_outer, 'C6-', label=r'$Re_\tau = 500$') 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.25,0.25]) plt.text(-4,0.25,'(b)') plt.subplot(525) plt.plot(C0093_yp_loc, C0093_T_uu_y_inner + C0093_T_uu_xz_inner, 'C4-') plt.plot(C0220_yp_loc, C0220_T_uu_y_inner + C0220_T_uu_xz_inner, 'C5-') plt.plot(C0500_yp_loc, C0500_T_uu_y_inner + C0500_T_uu_xz_inner, 'C6-') 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(C0093_yp_loc, C0093_T_uu_y_outer + C0093_T_uu_xz_outer, 'C4-') plt.plot(C0220_yp_loc, C0220_T_uu_y_outer + C0220_T_uu_xz_outer, 'C5-') plt.plot(C0500_yp_loc, C0500_T_uu_y_outer + C0500_T_uu_xz_outer, 'C6-') plt.ylabel(r'$T^+_{ \lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.25,0.25]) plt.text(-4,0.25,'(f)') plt.subplot(523) plt.plot(C0093_yp_loc, C0093_PI1_uu_inner, 'C4-') plt.plot(C0220_yp_loc, C0220_PI1_uu_inner, 'C5-') plt.plot(C0500_yp_loc, C0500_PI1_uu_inner, 'C6-') 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(C0093_yp_loc, C0093_PI1_uu_outer, 'C4-') plt.plot(C0220_yp_loc, C0220_PI1_uu_outer, 'C5-') plt.plot(C0500_yp_loc, C0500_PI1_uu_outer, 'C6-') plt.ylabel(r'$\Pi^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.25,0.25]) plt.text(-4,0.25,'(d)') plt.subplot(527) plt.plot(C0093_yp_loc, C0093_D_uu_inner, 'C4-') plt.plot(C0220_yp_loc, C0220_D_uu_inner, 'C5-') plt.plot(C0500_yp_loc, C0500_D_uu_inner, 'C6-') 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(C0093_yp_loc, C0093_D_uu_outer, 'C4-') plt.plot(C0220_yp_loc, C0220_D_uu_outer, 'C5-') plt.plot(C0500_yp_loc, C0500_D_uu_outer, 'C6-') plt.ylabel(r'$D^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.25,0.25]) plt.text(-4,0.25,'(h)') plt.subplot(529) plt.plot(C0093_yp_loc, -C0093_E_uu_inner, 'C4-') plt.plot(C0220_yp_loc, -C0220_E_uu_inner, 'C5-') plt.plot(C0500_yp_loc, -C0500_E_uu_inner, 'C6-') 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(C0093_yp_loc, -C0093_E_uu_outer, 'C4-') plt.plot(C0220_yp_loc, -C0220_E_uu_outer, 'C5-') plt.plot(C0500_yp_loc, -C0500_E_uu_outer, 'C6-') plt.ylabel(r'$-\epsilon^+_{\lambda \ge \delta}$') plt.xlabel(r'$y^+$') plt.xlim([0,20]) plt.ylim([-0.25,0.25]) plt.text(-4,0.25,'(j)') plt.subplots_adjust(hspace=0.4,wspace=0.3) plt.savefig('figure06.eps',dpi=300,bbox_inches='tight') plt.show() plt.close()
Out[5]:
