Commit: 4adefb72974a84f21d707e200bff7d97c97b6c7c - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revision: 4adefb72974a84f21d707e200bff7d97c97b6c7c (tree)
• Zeit: 2007-07-04 19:12:13
• Autor: iselllo
• Commiter: iselllo

Log Message

I added the code correlations_deposition.py which calculates the
deposition velocity in a cylinder using different correlations and plots
them together.

Ändern Zusammenfassung

• add: Python-codes/correlations_deposition.py

Diff

diff -r 9d7126879836 -r 4adefb72974a Python-codes/correlations_deposition.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/correlations_deposition.py	Wed Jul 04 10:12:13 2007 +0000
@@ -0,0 +1,105 @@
+#! /usr/bin/env python
+from scipy import *
+import pylab  # used to read the .csv file
+
+nu_air=2.1e-5
+Re=linspace(3000.,30000.,100)
+Diff=2.e-10 # in the range 1e-6-->1e-11
+R=0.05
+
+def corr_vouitsis(nu_air,Diff,Re):
+	Sc=nu_air/Diff # definition of Schmidt number
+	f=0.0791*Re**(-0.25)
+	#Now I use a correlation to calculate the Sherwood number
+	Sh=0.042*Re*f**0.5*Sc**(1./3.) # 
+	V_d=Sh*Diff/(2.*R) # sh=Sherwood number; R= tube radius
+	return V_d
+
+v_dep_vouitsis=corr_vouitsis(nu_air,Diff,Re)
+
+def corr_housiadas(nu_air,Diff,Re):
+	Sc=nu_air/Diff # definition of Schmidt number
+	U=Re*nu_air/(2.*R)
+	print 'U is', U
+	f=0.0791*Re**(-0.25)
+	#Sh=0.042*Re*f**0.5*Sc**(1./3.)
+	u_star=U*sqrt(f/2.)
+	V_d=u_star*0.2*Sc**(-2./3.)*Re**(-1./8.)
+	return V_d
+
+def corr_housiadas_gliniesky(nu_air,Diff,Re):
+	Sc=nu_air/Diff # definition of Schmidt number
+	U=Re*nu_air/(2.*R)
+	print 'U is', U
+	factor=(1./(((8./Re)**10.+(Re/36500.)**20.))**0.5+(2.21*log(Re/7.))**10.)**(1./5.)
+	f=2./factor
+	#f=0.0791*Re**(-0.25)
+	#Sh=0.042*Re*f**0.5*Sc**(1./3.)
+	u_star=U*sqrt(f/2.)
+	V_d=u_star*0.2*Sc**(-2./3.)*Re**(-1./8.)
+	return V_d
+
+def corr_vouitsis_gliniesky(nu_air,Diff,Re):
+	Sc=nu_air/Diff # definition of Schmidt number
+	#f=0.0791*Re**(-0.25)
+	factor=(1./(((8./Re)**10.+(Re/36500.)**20.))**0.5+(2.21*log(Re/7.))**10.)**(1./5.)
+	f=2./factor
+	#Now I use a correlation to calculate the Sherwood number
+	Sh=0.042*Re*f**0.5*Sc**(1./3.) # 
+	V_d=Sh*Diff/(2.*R) # sh=Sherwood number; R= tube radius
+	return V_d
+
+
+v_dep_housiadas=corr_housiadas(nu_air,Diff,Re)
+
+pylab.plot(Re,v_dep_vouitsis,Re,v_dep_housiadas)
+pylab.xlabel('Reynolds number')
+pylab.ylabel('Deposition velocity')
+pylab.legend(('voutsis','housiadas'))
+pylab.title('V_dep')
+pylab.grid(True)
+pylab.savefig('v_dep_comparison')
+
+pylab.hold(False)
+
+print 'max_vouitsis,max_housiadas,max_diff', max(v_dep_vouitsis),max(v_dep_housiadas),\
+max(abs(v_dep_housiadas-v_dep_vouitsis))
+
+
+v_dep_housiadas_gli=corr_housiadas_gliniesky(nu_air,Diff,Re)
+v_dep_vouitsis_gli=corr_vouitsis_gliniesky(nu_air,Diff,Re)
+
+pylab.plot(Re,v_dep_vouitsis_gli,Re,v_dep_housiadas_gli)
+pylab.xlabel('Reynolds number')
+pylab.ylabel('Deposition velocity')
+pylab.legend(('voutsis gliniesky','housiadas gliniesky'))
+pylab.title('V_dep')
+pylab.grid(True)
+pylab.savefig('v_dep_comparison_gliniesky')
+
+pylab.hold(False)
+
+
+pylab.plot(Re,v_dep_housiadas,Re,v_dep_housiadas_gli)
+pylab.xlabel('Reynolds number')
+pylab.ylabel('Deposition velocity')
+pylab.legend(('housiadas no Gliniesky','housiadas Gliniesky'))
+pylab.title('V_dep')
+pylab.grid(True)
+pylab.savefig('v_dep_comparison_housiadas')
+
+pylab.hold(False)
+
+
+pylab.plot(Re,v_dep_vouitsis,Re,v_dep_vouitsis_gli)
+pylab.xlabel('Reynolds number')
+pylab.ylabel('Deposition velocity')
+pylab.legend(('vouitsis no Gliniesky','vouitsis Gliniesky'))
+pylab.title('V_dep')
+pylab.grid(True)
+pylab.savefig('v_dep_comparison_vouitsis')
+
+pylab.hold(False)
+
+
+print 'So far so good'
\ No newline at end of file