diff options
Diffstat (limited to 'Textbook_Of_Heat_Transfer/Chapter_4_Principles_of.ipynb')
| -rwxr-xr-x | Textbook_Of_Heat_Transfer/Chapter_4_Principles_of.ipynb | 333 |
1 files changed, 333 insertions, 0 deletions
diff --git a/Textbook_Of_Heat_Transfer/Chapter_4_Principles_of.ipynb b/Textbook_Of_Heat_Transfer/Chapter_4_Principles_of.ipynb new file mode 100755 index 00000000..ed17906e --- /dev/null +++ b/Textbook_Of_Heat_Transfer/Chapter_4_Principles_of.ipynb @@ -0,0 +1,333 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:c0b8814a6ac7ab401db03352df216c35ffd27c0ad0849fd68734e0e65b766ec4"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 4: Principles of Fluid Flow"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.1 , Page no:172"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "L = 3 ; #Length\n",
+ "D = 0.01 ; #ID\n",
+ "V = 0.2 ; #Average Velocity\n",
+ "rho =999.7 ; #kg /m^ 3\n",
+ "\n",
+ "#calculations\n",
+ "v =1.306 * 10** -6 ; #m^2/ s\n",
+ "ReD =0.2*0.01/(1.306*10** -6) ;\n",
+ "f = 16/ ReD ;\n",
+ "deltap = 4*f*(L/D)*( rho*V **2) /2;\n",
+ "Vmax = 2*V;\n",
+ "\n",
+ "#result\n",
+ "print\"Pressure drop is\",round(deltap,4),\"Pa\";\n",
+ "print\"Maximum velocity is\",round(Vmax,3),\"m/s\";"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Pressure drop is 250.6768 Pa\n",
+ "Maximum velocity is 0.4 m/s\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.2(a) , Page no:180"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "L = 3 ; #m\n",
+ "D = 0.01 ; #m\n",
+ "V = 0.2 ;#m/ s\n",
+ "rho = 971.8 ; #kg /m^ 3\n",
+ "\n",
+ "#calculations\n",
+ "v = 0.365 * 10** -6 ; #m^2/ s\n",
+ "ReD = D*V/v;\n",
+ "f =0.079*( ReD )**( -0.25) ;\n",
+ "deltap = (4*f*L* rho *V**2) /(D *2) ;\n",
+ "x = ((f/2) **0.5) *V ;\n",
+ "yplus = 0.005* x /(0.365*10** -6) ;\n",
+ "Vmax = x *(2.5* math.log ( yplus ) + 5.5) ;\n",
+ "ratio = Vmax /V;\n",
+ "\n",
+ "#result\n",
+ "print\"(a) If the temperature of water is increased to 80 degree C\";\n",
+ "print\"Pressure drop is\",round(deltap,4),\"Pa\";\n",
+ "print\"Vmax =\",round(Vmax,4),\"m/s\";\n",
+ "print\"Vmax/Vbar =\",round(ratio,4);"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "(a) If the temperature of water is increased to 80 degree C\n",
+ "Pressure drop is 214.1563 Pa\n",
+ "Vmax = 0.2515 m/s\n",
+ "Vmax/Vbar = 1.2575\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.2(b) , Page no:180"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "L = 3 ; #m\n",
+ "D = 0.01 ; #m\n",
+ "V = 0.2 ; #m/s\n",
+ "V1 =0.7;\n",
+ "v1 = 1.306 * 10** -6 ; #m^2/ s\n",
+ "V1 =0.7; #m/ s\n",
+ "\n",
+ "#calculations\n",
+ "ReD1 =V1*D /(1.306*10** -6) ;\n",
+ "f1 = 0.079*( ReD1 )**( -0.25) ;\n",
+ "deltap1 = (4* f1*L *999.7*0.7**2) /(0.01*2) ;\n",
+ "x1 = (( f1 /2)**0.5) *V1 ;\n",
+ "y1plus = 0.005* x1 /( v1);\n",
+ "Vmax1 = x1 *(2.5* math.log ( y1plus ) + 5.5) ;\n",
+ "ratio1 = Vmax1 /V1;\n",
+ "\n",
+ "#result\n",
+ "print\"(b) If the velocity is increased to 0.7 \";\n",
+ "print\"Reynolds no is\",round(ReD1,4);\n",
+ "print\"Pressure drop is\",round(deltap1,4),\"pa\";\n",
+ "print\"y+ at centre line =\",round(y1plus,4);\n",
+ "print\"Vmax is\",round(Vmax1,4),\"pa\";\n",
+ "print\"Vmax/Vbar =\",round(ratio1,4);"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "(b) If the velocity is increased to 0.7 \n",
+ "Reynolds no is 5359.8775\n",
+ "Pressure drop is 2713.6598 pa\n",
+ "y+ at centre line = 182.087\n",
+ "Vmax is 0.8804 pa\n",
+ "Vmax/Vbar = 1.2577\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.3 , Page no:181"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "P = 80 * 10**3 ; #Pa\n",
+ "L = 10 ; #m\n",
+ "Vbar = 1.9 ; #m/s\n",
+ "l = 0.25 ; #m\n",
+ "b = 0.15 ; #m\n",
+ "rho = 1.185 ; #kg /m^ 3\n",
+ "\n",
+ "#calculations\n",
+ "mew = 18.35 * 10** -6 ; #kg /m s\n",
+ "rho1 = rho *(80/101.3) ; #kg /m^ 3\n",
+ "r = b/l;\n",
+ "De = (4* l/2*b /2) /(l/2 + b /2) ;\n",
+ "Dl = (2/3 + 11/24*0.6*(2 -0.6)) * De ;\n",
+ "Re = rho1 * Dl * Vbar / mew ;\n",
+ "f = 0.079*( Re ** -0.25) ;\n",
+ "deltaP = 4*f*(L/ Dl )*( rho1 *( Vbar **2) /2) ;\n",
+ "power = deltaP *( Vbar *l*b)\n",
+ "\n",
+ "#result\n",
+ "print\"Reynolds no =\",round(Re,4);\n",
+ "print\"f =\",round(f,4);\n",
+ "print\"Pressure drop =\",round(deltaP,4),\"Pa\";\n",
+ "print\"Power required =\",round(power,4),\"W\";"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Reynolds no = 19107.1453\n",
+ "f = 0.0067\n",
+ "Pressure drop = 2.3024 Pa\n",
+ "Power required = 0.164 W\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.4 , Page no:189"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "l = 2 ; #m\n",
+ "b = 1 ; #m\n",
+ "V = 1 ; #m/s\n",
+ "rho = 1.060 ; #kg /m^ 3\n",
+ "x = 1.5 ; #m^2/s\n",
+ "\n",
+ "#calculations\n",
+ "v = 18.97 * 10** -6 ;\n",
+ "Re = V*x/v; #Reynolds number\n",
+ "d = 5*x/( Re **(1/2) ) *1000 ;\n",
+ "Rel = V*l/v;\n",
+ "cf = 1.328* Rel** -(1/2) ; #drag coefficient\n",
+ "Fd = 0.00409*(1/2) *rho *(2* l*b) *1**2;\n",
+ "\n",
+ "#result\n",
+ "print\"Thickness of Boundary layer at x =1.5 is \",round(d,4),\"mm\"\n",
+ "print\"Drag Coefficient cf =\",round(cf,5);\n",
+ "print\"Drag Force FD =\",round(Fd,6),\"N\";"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Thickness of Boundary layer at x =1.5 is 26.6716 mm\n",
+ "Drag Coefficient cf = 0.00409\n",
+ "Drag Force FD = 0.008671 N\n"
+ ]
+ }
+ ],
+ "prompt_number": 5
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 4.5 , Page no:195"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "from __future__ import division\n",
+ "\n",
+ "#Variable declaration\n",
+ "l = 2 ; #m\n",
+ "v = 4 ; #m/s\n",
+ "\n",
+ "#calculations\n",
+ "mew = 18.1*10** -6; #N s /m^ 2\n",
+ "rho = 1.205*1.5; #kg /m^ 3\n",
+ "Rel = rho*v*l/ mew ;\n",
+ "Cf = 0.074*(7.989*10**5) **( -0.2) - 1050/ Rel ;\n",
+ "Df = Cf *1/2* rho*l*v **2;\n",
+ "x = 3*10**5 * (18.1*10** -6) /(1.808*4) ;\n",
+ "\n",
+ "#result\n",
+ "print\"Drag coefficieent is\",round(Cf,6);\n",
+ "print\"Drag force per meter width =\",round(Df,6),\"N\";\n",
+ "print\"Value of xc is\",round(x,6),\"m\";"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Drag coefficieent is 0.003569\n",
+ "Drag force per meter width = 0.103221 N\n",
+ "Value of xc is 0.75083 m\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |