Corrected Python scripts to plot Rayleigh collapse results

examples/rayleighcollapse/plot_result_emissions.py

@@ -14,7 +14,11 @@
 plt.rcParams['font.size']=10
 
 cm = 1/2.54
-tc = 0.09136264691 # Rayleigh collapse time
+R0 = 1 # Initial radius
+pinf = 1.0e5 # Ambient pressure
+pG0 = 1.0e3 # Initial gas pressure
+tc = 0.915*R0*np.sqrt(997/(pinf-pG0)) # Rayleigh collapse time
+uc = R0/tc # Characteristic velocity
 
 # Load the computed results.
 Bubble = np.genfromtxt(path + "Gilmore_R1.000e+00.txt", delimiter=" ")
@@ -59,8 +63,8 @@
 ax1.set_ylim(ymin=0,ymax=1)
 ax1.set_yticks([0,0.5,1])
 ax1.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax1.plot(MbR[:, 0]/tc, MbR[:, 2], linestyle='solid',linewidth=2,  color='mediumseagreen')
-ax1.plot(Bubble[:, 1]/tc, Bubble[:, 3], linestyle='solid',linewidth=0.75,  color='navy')
+ax1.plot(MbR[:, 0]/tc, MbR[:, 2]/R0, linestyle='solid',linewidth=2,  color='mediumseagreen')
+ax1.plot(Bubble[:, 1]/tc, Bubble[:, 3]/R0, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax3.set_yscale('log')
 ax3.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$p_\mathrm{G}(t)/p_\mathrm{G,0}$')
@@ -69,19 +73,19 @@
 ax3.set_ylim(ymin=1,ymax=1e6)
 ax3.set_yticks([1,1e2,1e4,1e6])
 ax3.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax3.plot(Mbpg[:, 0]/tc, Mbpg[:, 1]*1e-3, linestyle='solid',linewidth=2,  color='mediumseagreen', label=r'Navier-Stokes')
-ax3.plot(Bubble[:, 1]/tc, Bubble[:, 5]*1e-3, linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
+ax3.plot(Mbpg[:, 0]/tc, Mbpg[:, 1]/pG0, linestyle='solid',linewidth=2,  color='mediumseagreen', label=r'Navier-Stokes')
+ax3.plot(Bubble[:, 1]/tc, Bubble[:, 5]/pG0, linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
 
 ax3.legend(ncol=1,labelspacing=0.2,markerfirst=True,loc=(1.1,0.7),fontsize='small',facecolor='None',edgecolor='None',framealpha=1,frameon=True)
 
 ax4.set_yscale('log')
 ax4.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$\Delta p(t)/p_\infty$')
-ax4.set_xlim(xmin=1.004,xmax=1.0325)
+ax4.set_xlim(xmin=1.0,xmax=1.0275)
 ax4.set_ylim(ymin=5e-1,ymax=2e3)
 ax4.set_yticks([1,1e1,1e2,1e3])
 ax4.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax4.plot(Mp[:, 0]/tc, (Mp[:, 3]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax4.plot(Ep2[:, 0]/tc, (Ep2[:, 1]-1e5)*1e-5, linestyle='solid',linewidth=0.75,  color='navy')
+ax4.plot(Mp[:, 0]/tc, (Mp[:, 3]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax4.plot(Ep2[:, 0]/tc, (Ep2[:, 1]-1e5)/pinf, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax5.set_yscale('log')
 ax5.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$\Delta p(t)/p_\infty$')
@@ -90,8 +94,8 @@
 ax5.set_ylim(ymin=5e-1,ymax=2e3)
 ax5.set_yticks([1,1e1,1e2,1e3])
 ax5.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax5.plot(Mp[:, 0]/tc, (Mp[:, 4]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax5.plot(Ep5[:, 0]/tc, (Ep5[:, 1]-1e5)*1e-5, linestyle='solid',linewidth=0.75,  color='navy')
+ax5.plot(Mp[:, 0]/tc, (Mp[:, 4]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax5.plot(Ep5[:, 0]/tc, (Ep5[:, 1]-1e5)/pinf, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax6.set_yscale('log')
 ax6.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$\Delta p(t)/p_\infty$')
@@ -100,34 +104,34 @@
 ax6.set_ylim(ymin=5e-1,ymax=2e3)
 ax6.set_yticks([1,1e1,1e2,1e3])
 ax6.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax6.plot(Mp[:, 0]/tc, (Mp[:, 5]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
-ax6.plot(EpX[:, 0]/tc, (EpX[:, 1]-1e5)*1e-5,linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
+ax6.plot(Mp[:, 0]/tc, (Mp[:, 5]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
+ax6.plot(EpX[:, 0]/tc, (EpX[:, 1]-1e5)/pinf,linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
 
 ax7.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
-ax7.set_xlim(xmin=1.004,xmax=1.0325)
+ax7.set_xlim(xmin=1.0,xmax=1.0275)
 ax7.set_ylim(ymin=-8,ymax=8)
 ax7.set_yticks([-8,0,8])
 ax7.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax7.plot(Mu[:, 0]/tc, (Mu[:, 3]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax7.plot(Ep2[:, 0]/tc, (Ep2[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax7.plot(Mu[:, 0]/tc, (Mu[:, 3]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax7.plot(Ep2[:, 0]/tc, (Ep2[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax8.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
 ax8.set_xlim(xmin=0.92,xmax=1.12)
 ax8.set_xticks([0.92,1.02,1.12])
 ax8.set_ylim(ymin=-2,ymax=2)
 ax8.set_yticks([-2,0,2])
 ax8.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax8.plot(Mu[:, 0]/tc, (Mu[:, 4]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax8.plot(Ep5[:, 0]/tc, (Ep5[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax8.plot(Mu[:, 0]/tc, (Mu[:, 4]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax8.plot(Ep5[:, 0]/tc, (Ep5[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax9.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
 ax9.set_xlim(xmin=0.92,xmax=1.12)
 ax9.set_xticks([0.92,1.02,1.12])
 ax9.set_ylim(ymin=-0.6,ymax=0.6)
 ax9.set_yticks([-0.6,0,0.6])
 ax9.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax9.plot(Mu[:, 0]/tc, (Mu[:, 5]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
-ax9.plot(EpX[:, 0]/tc, (EpX[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax9.plot(Mu[:, 0]/tc, (Mu[:, 5]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
+ax9.plot(EpX[:, 0]/tc, (EpX[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax1.xaxis.set_label_coords(0.5,-0.24)
 ax3.xaxis.set_label_coords(0.5,-0.24)

examples/rayleighcollapse/plot_result_shock.py

@@ -14,7 +14,11 @@
 plt.rcParams['font.size']=10
 
 cm = 1/2.54
-tc = 0.09136264691 # Rayleigh collapse time
+R0 = 1 # Initial radius
+pinf = 1.0e5 # Ambient pressure
+pG0 = 1.0e2 # Initial gas pressure
+tc = 0.915*R0*np.sqrt(997/(pinf-pG0)) # Rayleigh collapse time
+uc = R0/tc # Characteristic velocity
 
 # Load the computed results.
 Bubble = np.genfromtxt(path + "Gilmore_R1.000e+00.txt", delimiter=" ")
@@ -45,8 +49,8 @@
 ax1.set_ylim(ymin=0,ymax=1)
 ax1.set_yticks([0,0.5,1])
 ax1.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax1.plot(MbR[:, 0]/tc, MbR[:, 2], linestyle='solid',linewidth=2,  color='mediumseagreen')
-ax1.plot(Bubble[:, 1]/tc, Bubble[:, 3], linestyle='solid',linewidth=0.75,  color='navy')
+ax1.plot(MbR[:, 0]/tc, MbR[:, 2]/R0, linestyle='solid',linewidth=2,  color='mediumseagreen')
+ax1.plot(Bubble[:, 1]/tc, Bubble[:, 3]/R0, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax3.set_yscale('log')
 ax3.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$p_\mathrm{G}(t)/p_\mathrm{G,0}$')
@@ -55,8 +59,8 @@
 ax3.set_ylim(ymin=1,ymax=1e8)
 ax3.set_yticks([1,1e2,1e4,1e6,1e8])
 ax3.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax3.plot(Mbpg[:, 0]/tc, Mbpg[:, 1]*1e-2, linestyle='solid',linewidth=2,  color='mediumseagreen', label=r'Navier-Stokes')
-ax3.plot(Bubble[:, 1]/tc, Bubble[:, 5]*1e-2, linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
+ax3.plot(Mbpg[:, 0]/tc, Mbpg[:, 1]/pG0, linestyle='solid',linewidth=2,  color='mediumseagreen', label=r'Navier-Stokes')
+ax3.plot(Bubble[:, 1]/tc, Bubble[:, 5]/pG0, linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
 
 ax3.legend(ncol=1,labelspacing=0.2,markerfirst=True,loc=(1.1,0.7),fontsize='small',facecolor='None',edgecolor='None',framealpha=1,frameon=True)
 
@@ -66,8 +70,8 @@
 ax4.set_ylim(ymin=1,ymax=5e3)
 ax4.set_yticks([1,1e1,1e2,1e3])
 ax4.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax4.plot(Mp[:, 0]/tc, (Mp[:, 3]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax4.plot(Ep2[:, 0]/tc, (Ep2[:, 1]-1e5)*1e-5, linestyle='solid',linewidth=0.75,  color='navy')
+ax4.plot(Mp[:, 0]/tc, (Mp[:, 3]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax4.plot(Ep2[:, 0]/tc, (Ep2[:, 1]-1e5)/pinf, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax5.set_yscale('log')
 ax5.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$\Delta p(t)/p_\infty$')
@@ -76,8 +80,8 @@
 ax5.set_ylim(ymin=1,ymax=5e3)
 ax5.set_yticks([1,1e1,1e2,1e3])
 ax5.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax5.plot(Mp[:, 0]/tc, (Mp[:, 4]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax5.plot(Ep5[:, 0]/tc, (Ep5[:, 1]-1e5)*1e-5, linestyle='solid',linewidth=0.75,  color='navy')
+ax5.plot(Mp[:, 0]/tc, (Mp[:, 4]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax5.plot(Ep5[:, 0]/tc, (Ep5[:, 1]-1e5)/pinf, linestyle='solid',linewidth=0.75,  color='navy')
 
 ax6.set_yscale('log')
 ax6.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$\Delta p(t)/p_\infty$')
@@ -86,34 +90,34 @@
 ax6.set_ylim(ymin=1,ymax=5e3)
 ax6.set_yticks([1,1e1,1e2,1e3])
 ax6.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax6.plot(Mp[:, 0]/tc, (Mp[:, 5]-1e5)*1e-5, linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
-ax6.plot(EpX[:, 0]/tc, (EpX[:, 1]-1e5)*1e-5,linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
+ax6.plot(Mp[:, 0]/tc, (Mp[:, 5]-1e5)/pinf, linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
+ax6.plot(EpX[:, 0]/tc, (EpX[:, 1]-1e5)/pinf,linestyle='solid',linewidth=0.75,  color='navy', label=r'APECSS')
 
 ax7.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
 ax7.set_xlim(xmin=0.996,xmax=1.02)
 ax7.set_ylim(ymin=-10,ymax=10)
 ax7.set_yticks([-10,0,10])
 ax7.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax7.plot(Mu[:, 0]/tc, (Mu[:, 3]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax7.plot(Ep2[:, 0]/tc, (Ep2[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax7.plot(Mu[:, 0]/tc, (Mu[:, 3]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax7.plot(Ep2[:, 0]/tc, (Ep2[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax8.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
 ax8.set_xlim(xmin=0.92,xmax=1.12)
 ax8.set_xticks([0.92,1.02,1.12])
 ax8.set_ylim(ymin=-4,ymax=4)
 ax8.set_yticks([-4,0,4])
 ax8.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax8.plot(Mu[:, 0]/tc, (Mu[:, 4]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen')
-ax8.plot(Ep5[:, 0]/tc, (Ep5[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax8.plot(Mu[:, 0]/tc, (Mu[:, 4]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen')
+ax8.plot(Ep5[:, 0]/tc, (Ep5[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax9.set(xlabel=r'$t/t_\mathrm{c}$',ylabel=r'$u(t)/u_\mathrm{c}$')
 ax9.set_xlim(xmin=0.92,xmax=1.12)
 ax9.set_xticks([0.92,1.02,1.12])
 ax9.set_ylim(ymin=-2,ymax=2)
 ax9.set_yticks([-2,0,2])
 ax9.grid(color='gainsboro', linestyle='-', linewidth=0.5)
-ax9.plot(Mu[:, 0]/tc, (Mu[:, 5]*tc), linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
-ax9.plot(EpX[:, 0]/tc, (EpX[:, 2]*tc), linestyle='solid',linewidth=0.75,  color='navy')
+ax9.plot(Mu[:, 0]/tc, (Mu[:, 5]/uc), linestyle='solid',linewidth=3,  color='mediumseagreen', label=r'Navier-Stokes')
+ax9.plot(EpX[:, 0]/tc, (EpX[:, 2]/uc), linestyle='solid',linewidth=0.75,  color='navy')
 
 ax1.xaxis.set_label_coords(0.5,-0.24)
 ax3.xaxis.set_label_coords(0.5,-0.24)