{
"metadata": {
"name": "",
"signature": "sha256:da0a57e05023d78d40f48e24ac486e01915e80a89fcf828d645acb89b5297ac2"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Chapter 13 : Miscellaneous Problems"
]
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.1 Page No : 282"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"W= 5000 \t#lb\n",
"vr= 6\n",
"e= 0.95\n",
"ep = 0.75\n",
"d= 9 \t#in\n",
"D= 45 \t#ft\n",
"t= 2 \t#min\n",
"v= 4.5 \t#ft/sec\n",
"\n",
"#CALCULATIONS\n",
"L= W*vr/(e*ep)\n",
"Pr= L/(math.pi*d**2/4)\n",
"s= D/vr\n",
"V= s*math.pi*ep**2/(4*t*60)\n",
"T= D/v\n",
"V1= s*math.pi*ep**2/4\n",
"V2= V*T\n",
"V3= V1-V2\n",
"\n",
"#RESULTS\n",
"print 'Pressure on ram = %.f ln/in**2 '%(Pr)\n",
"print ' Pump duty = %.4f cusec'%(V)\n",
"print ' Minimum capacity if accumulator = %.3f ft**3 '%(round(V3))\n",
"\n",
"# rounding off error"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Pressure on ram = 662 ln/in**2 \n",
" Pump duty = 0.0258 cusec\n",
" Minimum capacity if accumulator = 3.000 ft**3 \n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.2 Page No : 283"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"P1= 1100. \t#lb/in**2\n",
"P2= 85. \t#lb/in**2\n",
"f= 0.01\n",
"g= 32.2 \t#ft/sec**2\n",
"l= 1600. \t#ft\n",
"r= 1./8\n",
"W= 2500. \t#lb\n",
"d= 6. \t#in\n",
"\n",
"#CALCULATIONS\n",
"L= W*d\n",
"P= L*2.31/(math.pi*(d/2)**2)\n",
"s1= P1*2540/1100\n",
"s2= P2*196/85\n",
"vp= math.sqrt((s1-s2-P)/(4*f*l/(2*g*r)))\n",
"V= vp/16\n",
"Vl= V*d\n",
"Vp= math.sqrt((s1/3)/(4*f*l/(2*g*r)))\n",
"vl= Vp*d/16\n",
"Hr= s1-(s1/3)-s2\n",
"Lr= Hr*math.pi*(d/2)**2/(2.31*d)\n",
"\n",
"#RESULTS\n",
"print \"In case 1 velocity of load = %.2f ft/sec\"%(Vl)\n",
"print \" In case 2 velocity of load = %.2f ft/sec\"%(vl)\n",
"print ' Load to be lifted = %d lb '%(Lr)\n",
"\n",
"# note : roundin off error"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"In case 1 velocity of load = 4.45 ft/sec\n",
" In case 2 velocity of load = 3.87 ft/sec\n",
" Load to be lifted = 3054 lb \n"
]
}
],
"prompt_number": 5
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.3 Page No : 284"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math \n",
"\n",
"#initialisation of variables\n",
"bhp= 1500 \t#h.p\n",
"e= 0.86\n",
"h1= 300 \t#ft\n",
"h2= 15 \t#ft\n",
"w= 62.4 \t#lb/ft**3\n",
"t= 30 \t#days\n",
"t1= 10 \t#hr\n",
"t2= 3 \t#months\n",
"f= 0.005\n",
"l= 1000 \t#ft\n",
"\n",
"#CALCULATIONS\n",
"WHP= bhp/e\n",
"Ha= h1-h2\n",
"W= WHP*550\n",
"Q= W/(Ha*w)\n",
"Qt= Q*36009*t1*t*t2\n",
"Qp= Qt/(3600*t*45)\n",
"d= (f*l*(Q/2)**2/(t1*h2))**(1./5)\n",
"\n",
"#RESULTS\n",
"print \"Minimum size of basin required = %.1e cusecs\"%Qt\n",
"print \"Pump Discharge : %.f cusecs\"%(Qp/10)\n",
"print 'Diameter = %.2f ft '%(d)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Minimum size of basin required = 1.7e+09 cusecs\n",
"Pump Discharge : 36 cusecs\n",
"Diameter = 1.89 ft \n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.4 Page No : 285"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math \n",
"\n",
"#initialisation of variables\n",
"l= 140 \t#ft\n",
"P= 70 \t#percent\n",
"V= 3*10**8 \t#ft**3\n",
"w= 62.4 \t#lb/ft**3\n",
"SBD= 4.9*10**8 \t#ft**3\n",
"Q= 162 \t#cuses\n",
"s= 12.2*10**6 \t#ft**3/day\n",
"\n",
"#CALCULATIONS\n",
"O= Q*w*l*(P/1000.)/550.\n",
"\n",
"#RESULTS\n",
"print 'Size of reservoir= %.2e ft**3'%(SBD)\n",
"print ' output = %.f h.p '%(O)\n",
"print ' output = %d h.p '%(Q)\n",
"\n",
"# rounding off error"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Size of reservoir= 4.90e+08 ft**3\n",
" output = 180 h.p \n",
" output = 162 h.p \n"
]
}
],
"prompt_number": 7
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.5 Page No : 287"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"Q= 140 \t#cuses\n",
"w= 62.4 \t#lb/ft**3\n",
"l= 140 \t#ft\n",
"P= 70 \t#percent\n",
"k= 1.6\n",
"v= 3*10**8\n",
"\n",
"#CALCULATIONS\n",
"rv= k*v\n",
"HP= Q*l*w*(P/1000.)/550.\n",
"\n",
"#RESULTS\n",
"print 'Required size of reservoir = %.1e ft**3 '%(rv)\n",
"print ' horsepower = %.f h.p '%(HP)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Required size of reservoir = 4.8e+08 ft**3 \n",
" horsepower = 156 h.p \n"
]
}
],
"prompt_number": 6
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.6 Page No : 288"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math \n",
"\n",
"#initialisation of variables\n",
"P= 10. \t#lb/in**2\n",
"r1= 0.5 \t#ft\n",
"r= 0.25 \t#ft\n",
"f= 42.3 \t#ft/sec\n",
"b= 1./40\n",
"Tt= 1400. \t#lb\n",
"\n",
"#CALCULATIONS\n",
"Q= 2*math.pi*r*b*f\n",
"p1= 34+P\n",
"Fu= p1*math.pi*(r-(r/4))*144/2.3\n",
"Fr= Fu-Tt\n",
"\n",
"#RESULTS\n",
"print 'Quantity = %.2f cusecs '%(Q)\n",
"print ' Resultant force on the plate = %.f lb '%(round(Fr,-1))\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Quantity = 1.66 cusecs \n",
" Resultant force on the plate = 220 lb \n"
]
}
],
"prompt_number": 8
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.7 Page No : 289"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"r= 0.5 \t#ft\n",
"N= 300\n",
"w= 62.4 \t#lb/ft**3\n",
"g= 32.2 \t#ft/sec**2\n",
"\n",
"#CALCULATIONS\n",
"A= N*2*math.pi/60\n",
"Ft= math.pi*A**2*r**4*w/(4*g)\n",
"\n",
"#RESULTS\n",
"print 'total force = %.1f lb '%(Ft)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"total force = 93.9 lb \n"
]
}
],
"prompt_number": 8
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.8 Page No : 292"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"d= 4. \t#in\n",
"h= 12. \t#in\n",
"h1= 9. \t#in\n",
"g= 32. \t#ft/sec**2\n",
"\n",
"#CALCULATIONS\n",
"H= 2*(1-(h1/h))\n",
"A= math.sqrt((H*2*g/((d/24)**2)))\n",
"A1= math.sqrt((H*2*g*2/((d/24)**2)))\n",
"\n",
"#RESULTS\n",
"print 'speed when the axial is zero = %.f radn/sec '%(A)\n",
"print ' speed when the axial is zero = %.f radn/sec '%(A1)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"speed when the axial is zero = 34 radn/sec \n",
" speed when the axial is zero = 48 radn/sec \n"
]
}
],
"prompt_number": 10
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.10 Page No : 295"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"P= 14.7 \t#lb/in**2\n",
"T= 15. \t#C\n",
"v= 350. \t#ft/sec\n",
"R= 0.714\n",
"\n",
"#CALCULATIONS\n",
"P1= P*144\n",
"r= 3091*(273+T)\n",
"d1= P1/r\n",
"r1= r+(v**2/7)\n",
"P2= (r1*d1/(P1**R))**(1/(1-R))/144\n",
"dP= P2-P\n",
"T2= r1/3091\n",
"dT= T2-(273+T)\n",
"\n",
"#RESULTS\n",
"print 'rise in pressure = %.f lb/in**2 '%(dP)\n",
"print ' rise in temperature = %.1f C '%(dT)\n",
"\n",
"# rounding off error"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"rise in pressure = 1 lb/in**2 \n",
" rise in temperature = 5.7 C \n"
]
}
],
"prompt_number": 9
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.11 Page No : 297"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import math \n",
"\n",
"#initialisation of variables\n",
"T= 27. \t#C\n",
"P = 33. \t#lb/in**2\n",
"p1= 14.7 \t#lb/in**2\n",
"w= 250. \t#lb\n",
"g= 32.2 \t#ft/sec**2\n",
"Cd= 0.99\n",
"r= 1.4\n",
"\n",
"#CALCULATIONS\n",
"w1= P*144/(96*(273+T))\n",
"d= p1*144/(96*(273+T))\n",
"W= d*w/60\n",
"d= math.sqrt(W*4/(Cd*math.pi*math.sqrt(2*g*P*144*(r/(r-1))*w1*(0.528**(2/1.4)-0.528**(2.4/1.4)))))*12\n",
"\n",
"#RESULTS\n",
"print 'Diameter = %.3f in '%(d)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Diameter = 0.722 in \n"
]
}
],
"prompt_number": 13
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.12 Page No : 299"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math \n",
"\n",
"#initialisation of variables\n",
"sp= 13.6\n",
"hm= 800. \t#mm\n",
"d= 3. \t#in\n",
"r= 1.4\n",
"R= 1385. \t#ft-lb/lb/C\n",
"w= 62.4 \t#lb/ft**3\n",
"T= 15. \t#C\n",
"hm1= 765. \t#mm\n",
"r1= 9.\n",
"g= 32.2 \t#ft/sec**2\n",
"\n",
"#CALCULATIONS\n",
"p1= hm*sp*w/304.8\n",
"r2= (273+T)*R\n",
"w1= p1/r2\n",
"k= hm/hm1\n",
"v1= math.sqrt((2*g*r*r2*(1-k**0.286))/((1-r)*(r1**2*k**1.43-1)))\n",
"W= v1*w1*3600*(math.pi/64)\n",
"\n",
"#RESULTS\n",
"print 'Weight flowing = %.f lb/hr '%(W)\n",
"\n",
"#The answer is a bit different due to rounding off error in textbook\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Weight flowing = 115 lb/hr \n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 13.13 Page No : 301"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import math \n",
"\n",
"#initialisation of variables\n",
"p= 160. \t#lb/in**2\n",
"d= 1./3 \t#ft\n",
"T= 15. \t#C\n",
"R= 96. \n",
"V= 120. \t#ft**3\n",
"f= 0.004\n",
"a= 60*math.pi\n",
"l= 10560. \t#ft\n",
"g= 32.2 \t#ft/sec**2\n",
"\n",
"#CALCULATIONS\n",
"p1= p*144\n",
"w1= p*144/(R*(273+T))\n",
"v1= V*36/a\n",
"p2= math.sqrt(p1**2-((2*4*f*p1*w1*v1**2*l)/(2*g*d)))/144\n",
"v2= p*v1/p2\n",
"\n",
"#RESULTS\n",
"print ' pressure = %.1f lb/in**2 '%(p2)\n",
"print ' velocity = %.1f ft/sec '%(v2)\n",
"\n",
"#The answer is a bit different due to rounding off error in textbook\n",
"\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
" pressure = 134.0 lb/in**2 \n",
" velocity = 27.4 ft/sec \n"
]
}
],
"prompt_number": 12
}
],
"metadata": {}
}
]
}