summaryrefslogtreecommitdiff
path: root/Introduction_To_Chemical_Engineering/ch8.ipynb
diff options
context:
space:
mode:
authordebashisdeb2014-06-20 15:42:42 +0530
committerdebashisdeb2014-06-20 15:42:42 +0530
commit83c1bfceb1b681b4bb7253b47491be2d8b2014a1 (patch)
treef54eab21dd3d725d64a495fcd47c00d37abed004 /Introduction_To_Chemical_Engineering/ch8.ipynb
parenta78126bbe4443e9526a64df9d8245c4af8843044 (diff)
downloadPython-Textbook-Companions-83c1bfceb1b681b4bb7253b47491be2d8b2014a1.tar.gz
Python-Textbook-Companions-83c1bfceb1b681b4bb7253b47491be2d8b2014a1.tar.bz2
Python-Textbook-Companions-83c1bfceb1b681b4bb7253b47491be2d8b2014a1.zip
removing problem statements
Diffstat (limited to 'Introduction_To_Chemical_Engineering/ch8.ipynb')
-rw-r--r--Introduction_To_Chemical_Engineering/ch8.ipynb40
1 files changed, 0 insertions, 40 deletions
diff --git a/Introduction_To_Chemical_Engineering/ch8.ipynb b/Introduction_To_Chemical_Engineering/ch8.ipynb
index a18b8152..d67175c5 100644
--- a/Introduction_To_Chemical_Engineering/ch8.ipynb
+++ b/Introduction_To_Chemical_Engineering/ch8.ipynb
@@ -30,13 +30,11 @@
"\n",
"import math \n",
"\n",
- "# Variables\n",
"pressure_difference = 3.4 #in mm water\n",
"pressure = 1.0133*10**5 #in pa\n",
"temperatue = 293. #in K\n",
"mass_of_air = 29. #in Kg\n",
"\n",
- "# Calculations and Results\n",
"density_air = pressure/(temperatue*8314)*mass_of_air #in kg/m3\n",
"print \"Density of air = %f kg/cu m\"%(density_air)\n",
"\n",
@@ -76,21 +74,16 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find viscosity of oil \n",
"\n",
"import math \n",
"\n",
- "# Variables\n",
"diameter=0.6; #in m\n",
"disk_distance=1.25*10**-3; #in m\n",
"speed=5.; #revolutions/min\n",
"torque=11.5; #in Joules\n",
"\n",
- "# Calculations\n",
- "#we know that torque= pi*omega*viscosity*radius**4/2*disc_distance\n",
"viscosity=(2*disk_distance*torque)/(3.14*(10*3.14)*(diameter/2)**4);\n",
"\n",
- "# Results\n",
"print \"viscosity = %f Pa-s\"%(viscosity)\n"
],
"language": "python",
@@ -118,22 +111,18 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the viscosity of solution using given parameters\n",
"\n",
"import math \n",
- "# Variables\n",
"diameter =10.; #in mm\n",
"density_of_solution = 1750.; #in kg/m3\n",
"density_of_air = 1.2; #in kg/m3\n",
"velocity = 0.9; #in mm/s\n",
"\n",
- "# Calculations and Results\n",
"viscosity = (density_of_solution-density_of_air)*9.8*(diameter*10**-3)**2/(18*velocity*10**-3); #expression for finding viscosity\n",
"\n",
"print \"viscosity of solution = %f Pa-s\"%(viscosity)\n",
"\n",
"\n",
- "#checking stoke's region validity\n",
"v=(0.2*viscosity)/(density_of_solution*diameter*10**-3);\n",
"if v>0.9 :\n",
" print \"system follows stokes law\"\n"
@@ -164,11 +153,9 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the flow rate in an orifice\n",
"\n",
"import math \n",
"\n",
- "# Variables\n",
"density_of_water = 1000.; #in kg/m3\n",
"viscosity = 1*10.**-3; #in Pa-s\n",
"pipe_diameter = 250.; #in mm\n",
@@ -176,14 +163,11 @@
"density_of_mercury = 13600.; # in mm\n",
"manometer_height = 242.; #in mm\n",
"\n",
- "# Calculations and Results\n",
"height_water_equivalent = (density_of_mercury-density_of_water)*(manometer_height*10**-3)/(density_of_water) #in m\n",
"\n",
- "#assuming Re>30000\n",
"Co = 0.61;\n",
"velocity = Co*(2*9.8*height_water_equivalent/(1-(orifice_diameter/pipe_diameter)**4))**0.5; #in m/s\n",
"\n",
- "#checking Reynold's number\n",
"Re = (orifice_diameter*10**-3*velocity*density_of_water)/viscosity;\n",
"print \"reynolds number = %f which is greater than 30000\"%(Re)\n",
"\n",
@@ -220,10 +204,8 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the coefficient of discharge for converging cone\n",
"\n",
"import math \n",
- "# Variables\n",
"pipe_diameter=0.15; #in m\n",
"venturi_diameter=0.05; #in m\n",
"pressure_drop=0.12; #m of water\n",
@@ -231,15 +213,12 @@
"density = 1000.; #in kg/m3\n",
"viscosity = 0.001 #in Pa-s\n",
"\n",
- "# Calculations and Results\n",
"velocity = ((4./3.14)*flow_rate)/(venturi_diameter**2*density);\n",
"print \"velociy = %f m/s\"%(velocity)\n",
"\n",
- "#calculating coefficient of discharge\n",
"Cv=velocity*((1-(venturi_diameter/pipe_diameter)**4)/(2*9.8*pressure_drop))**0.5;\n",
"print \"coefficient of discharge = %f\"%(Cv)\n",
"\n",
- "#calculating reynold's number\n",
"Re = velocity*(venturi_diameter/pipe_diameter)**2*pipe_diameter*density/viscosity;\n",
"print \"reynolds No = %f\"%(Re)\n"
],
@@ -270,10 +249,8 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find pA and pB\n",
"\n",
"import math \n",
- "# Variables\n",
"h1=0.66; #in m\n",
"h2=0.203; #in m\n",
"h3=0.305 #in m\n",
@@ -282,7 +259,6 @@
"s1=0.83;\n",
"s2=13.6;\n",
"\n",
- "# Calculations and Results\n",
"print (\"part 1\")\n",
"pA=pB+(h2*s2-(h1-h3)*s1)*density*9.81; #in Pa\n",
"print \"pressure at A = %f Pa\"%(pA)\n",
@@ -322,17 +298,14 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the rate of oil flow in l/s\n",
"\n",
"import math \n",
- "# Variables\n",
"density_oil=900.; #in kg/m3\n",
"viscosity_oil=38.8*10**-3; #in Pa-s\n",
"density_water = 1000.; #in kg/m3\n",
"diameter=0.102 #in m\n",
"manometer_reading=0.9; #m of water\n",
"\n",
- "# Calculations and Results\n",
"delta_H=manometer_reading*(density_water-density_oil)/density_oil;\n",
"print \"manometer reading as m of oil = %f m\"%(delta_H)\n",
"\n",
@@ -378,20 +351,16 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the maximum capacity of keroscene\n",
"\n",
"import math \n",
- "# Variables\n",
"flow_rate_steel=1.2; #l/s\n",
"density_steel=7.92;\n",
"density_kerosene=0.82;\n",
"density_water=1;\n",
"\n",
- "# Calculations\n",
"flow_rate_kerosene =(((density_steel-density_kerosene)/density_kerosene)/((density_steel-density_water)/density_water))**0.5*flow_rate_steel\n",
"\n",
"\n",
- "# Results\n",
"print \"maximum_flow rate of kerosene = %f litre/s\"%(flow_rate_kerosene)\n"
],
"language": "python",
@@ -419,28 +388,19 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "#to find the rate of flow of flue gas\n",
"\n",
"from scipy.optimize import fsolve \n",
"import math \n",
- "# Variables\n",
"initial_CO2 = 0.02; #weight fraction\n",
"flow_rate_CO2 = 22.5; #gm/s\n",
"final_CO2=0.031; #weight fraction\n",
"\n",
- "#flow rate of flue gas =x\n",
- "#amount of CO2 entering = 0.02*x\n",
- "#amount of CO2 leaving = 0.02x+0.0225\n",
- "#amount of gas leaving = x+0.0225\n",
- "#amount of CO2 leaving = 0.031*(x+0.0225)\n",
"\n",
- "# Calculations\n",
"def f(x): \n",
"\t return initial_CO2*x+0.0225 - 0.031*(x+0.0225)\n",
"\n",
"flow_rate_flue_gas=fsolve(f,0)\n",
"\n",
- "# Results\n",
"print \"flow rate of flue gas = %f kg/s\"%(flow_rate_flue_gas)\n",
"\n"
],