diff options
Diffstat (limited to 'A_Heat_Transfer_Text_Book/Chapter1.ipynb')
| -rwxr-xr-x | A_Heat_Transfer_Text_Book/Chapter1.ipynb | 332 |
1 files changed, 332 insertions, 0 deletions
diff --git a/A_Heat_Transfer_Text_Book/Chapter1.ipynb b/A_Heat_Transfer_Text_Book/Chapter1.ipynb new file mode 100755 index 00000000..9c3e159f --- /dev/null +++ b/A_Heat_Transfer_Text_Book/Chapter1.ipynb @@ -0,0 +1,332 @@ +{
+ "metadata": {
+ "name": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter 1 - \"Introduction\""
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.1, Page number: 13"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "\n",
+ "#Variables\n",
+ "k=35; #Thermal Conductivity, [W/m*K]\n",
+ "T1=110 # Temperature of front[C]\n",
+ "T2=50; # Temperature of back,[C]\n",
+ "A=0.4 #area of slab,[m**2]\n",
+ "x=0.03; #Thickness of slab,[m]\n",
+ "\n",
+ "#Calculations\n",
+ "q=-k*(T2-T1)/(1000*x); #formula for heat flux[KW/m^2]\n",
+ "Q=q*A; #formula for heat transfer rate[KW]\n",
+ "\n",
+ "#Results\n",
+ "print \"Heat flux is:\",q,\"KW/m^2\\n\"\n",
+ "print \"Heat transfer rate is:\",Q,\"KW \\n\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Heat flux is: 70.0 KW/m^2\n",
+ "\n",
+ "Heat transfer rate is: 28.0 KW \n",
+ "\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.2, Page number: 16"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "from sympy import solve,symbols\n",
+ "\n",
+ "#Variables\n",
+ "x=symbols('x');\n",
+ "k1=372; # Thermal Conductivity of slab,W/m*K\n",
+ "x1=0.003; # Thickness of slab,m\n",
+ "x2=0.002 # Thickness of steel,m\n",
+ "k2=17; # Thermal Conductivity of steel,W/m*K\n",
+ "T1=400; # Temperature on one side,C\n",
+ "T2=100 #Temperature on other side,C\n",
+ "\n",
+ "#Calculations\n",
+ "Tcu=solve(x+2*x*(k1/x1)*(x2/k2)-(T1-T2),x);\n",
+ "#q=k1*(Tcu/x1)=k2*(Tss/x2);\n",
+ "Tss = Tcu[0]*(k1/x1)*(x2/k2); # formula for temperature gradient in steel side\n",
+ "Tcul=T1-Tss;\n",
+ "Tcur=T2+Tss;\n",
+ "q=k2*Tss/(1000*x2); # formula for heat conducted, kW\\m^2\n",
+ "\n",
+ "#Results\n",
+ "print \"Temperature on left copper side is :\",round(Tcul,3),\"C\\n\"\n",
+ "print \"Temperature on right copper side is :\",round(Tcur,3),\"C\\n\"\n",
+ "print \"Heat conducted through the wall is :\",round(q,3),\"kW\\m^2\\n\"\n",
+ "print \"Our initial approximation was accurate within a few percent.\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Temperature on left copper side is : 254.971 C\n",
+ "\n",
+ "Temperature on right copper side is : 245.029 C\n",
+ "\n",
+ "Heat conducted through the wall is : 1232.749 kW\\m^2\n",
+ "\n",
+ "Our initial approximation was accurate within a few percent.\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.3, Page number: 22"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "\n",
+ "#Variables\n",
+ "q1=6000; #Heat flux, W*m**-2\n",
+ "T1=120; #Heater Temperature, C\n",
+ "T2=70; #final Temperature of Heater, C\n",
+ "q2=2000; #final heat flux, W*m**-2\n",
+ "\n",
+ "#Calculations\n",
+ "h=q1/(T1-T2) #formula for average heat transfer cofficient\n",
+ "Tnew=T2+q2/h; #formula for new Heater temperature, C\n",
+ "\n",
+ "#Results\n",
+ "print \"Average Heat transfer coefficient is:\",h,\"W/(m^2*K)\\n\"\n",
+ "print \"New Heater Temperature is:\",round(Tnew,3),\"C\\n\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Average Heat transfer coefficient is: 120.0 W/(m^2*K)\n",
+ "\n",
+ "New Heater Temperature is: 86.667 C\n",
+ "\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.4, Page number: 25"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "from numpy import array\n",
+ "from numpy import linspace\n",
+ "import matplotlib.pyplot as plt\n",
+ "from pylab import *\n",
+ "%pylab inline\n",
+ "\n",
+ "#Variables\n",
+ "h=250; #Heat Transfer Coefficient, W/(m**2*K)\n",
+ "k=45; #Thermal Conductivity, W/(m*K)\n",
+ "c=180; #Heat Capacity, J/(kg*K)\n",
+ "a=9300; #density, kg/m**3\n",
+ "T1=200; #temperature, C\n",
+ "D=0.001; #diameter of bead, m\n",
+ "t1=linspace(0,5,50); #defining time interval of 0.1 seconds\n",
+ "T=linspace(0,5,50);\n",
+ "i=0;\n",
+ "\n",
+ "#Calculations\n",
+ "while i<50:\n",
+ " T[i]=T1-c*math.exp(-t1[i]/((a*c*D)/(6*h))); #Calculating temperature at each time in degree C\n",
+ " i=i+1;\n",
+ "\n",
+ "plt.plot(t1,T);\n",
+ "plt.xlabel(\"Time(in sec)\");\n",
+ "plt.ylabel(\"Temperature(in degree C)\");\n",
+ "plt.title(\"Thermocouple response to a hot gas flow\");\n",
+ "plt.show();\n",
+ "\n",
+ "Bi = h*(D/2)/k; #biot no.\n",
+ "\n",
+ "#Results\n",
+ "print \"The value of Biot no for this thermocouple is\",round(Bi,5);\n",
+ "print \"Bi is <0.1 and hence the thermocouple could be considered as a lumped heat capacity system and the assumption taken is valid.\\n\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "ename": "ImportError",
+ "evalue": "DLL load failed: %1 is not a valid Win32 application.",
+ "output_type": "pyerr",
+ "traceback": [
+ "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mImportError\u001b[0m Traceback (most recent call last)",
+ "\u001b[1;32m<ipython-input-4-48a98dda89ad>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 1\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0m__future__\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mdivision\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mmath\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 3\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0marray\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 4\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mlinspace\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mmatplotlib\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mpyplot\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mplt\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;32mC:\\Users\\jaidev\\Anaconda\\lib\\site-packages\\numpy\\__init__.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 151\u001b[0m \u001b[1;32mreturn\u001b[0m \u001b[0mloader\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0mpackages\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0moptions\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 152\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 153\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0madd_newdocs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 154\u001b[0m \u001b[0m__all__\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;34m'add_newdocs'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'ModuleDeprecationWarning'\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 155\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;32mC:\\Users\\jaidev\\Anaconda\\lib\\site-packages\\numpy\\add_newdocs.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 11\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0m__future__\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mdivision\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mabsolute_import\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mprint_function\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 12\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 13\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlib\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0madd_newdoc\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 14\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 15\u001b[0m \u001b[1;31m###############################################################################\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;32mC:\\Users\\jaidev\\Anaconda\\lib\\site-packages\\numpy\\lib\\__init__.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 6\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mversion\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mversion\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0m__version__\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 7\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 8\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m\u001b[0mtype_check\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[1;33m*\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 9\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m\u001b[0mindex_tricks\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[1;33m*\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 10\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m\u001b[0mfunction_base\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[1;33m*\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;32mC:\\Users\\jaidev\\Anaconda\\lib\\site-packages\\numpy\\lib\\type_check.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 9\u001b[0m 'common_type']\n\u001b[0;32m 10\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 11\u001b[1;33m \u001b[1;32mimport\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcore\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mnumeric\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0m_nx\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 12\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcore\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mnumeric\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0masarray\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0masanyarray\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0marray\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0misnan\u001b[0m\u001b[1;33m,\u001b[0m\u001b[0;31m \u001b[0m\u001b[0;31m\\\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 13\u001b[0m \u001b[0mobj2sctype\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mzeros\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;32mC:\\Users\\jaidev\\Anaconda\\lib\\site-packages\\numpy\\core\\__init__.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mversion\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mversion\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0m__version__\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 6\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mmultiarray\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 7\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mumath\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 8\u001b[0m \u001b[1;32mfrom\u001b[0m \u001b[1;33m.\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0m_internal\u001b[0m \u001b[1;31m# for freeze programs\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
+ "\u001b[1;31mImportError\u001b[0m: DLL load failed: %1 is not a valid Win32 application."
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.5, Page number: 32"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "from sympy import solve,symbols\n",
+ "\n",
+ "#Variables\n",
+ "x=symbols('x');\n",
+ "T1=293; #Temperature of air around thermocouple, K\n",
+ "T2=373; #Wall temperature, K\n",
+ "h=75; #Average Heat Transfer Coefficient, W/(m**2*K)\n",
+ "s=5.67*10**-8; #stefan Boltzman constant, W/(m**2*K**4)\n",
+ "\n",
+ "#Calculations\n",
+ "x=solve((h*(x-T1)+s*(x**4-T2**4)),x);\t #Calculating Thermocouple Temperature, K\n",
+ "y=x[1]-273;\t\t\t\t #Thermocouple Temperature, C\n",
+ "\n",
+ "#Results\n",
+ "print \"Thermocouple Temperature is :\",round(y,3),\"C\\n\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Thermocouple Temperature is : 28.395 C\n",
+ "\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 3,
+ "metadata": {},
+ "source": [
+ "Example 1.6, Page number: 34"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "import math\n",
+ "from sympy import solve,symbols\n",
+ "\n",
+ "#Variables\n",
+ "x=symbols('x');\n",
+ "e=0.4; #emissivity\n",
+ "T1=293; #Temperature of air around Thermocouple, K\n",
+ "T2=273; #wall Temperature, K\n",
+ "h=75; #Average Heat Transfer Coefficient, W/(m**2*K)\n",
+ "s=5.6704*10**-8; #stefan Boltzman constant, W/(m**2*K**4)\n",
+ "\n",
+ "#Calculations\n",
+ "z=solve(((s*e*((373)**4 - (x)**4)) - h*(x-293)),x);\t#Calculating Thermocouple Temperature, K\n",
+ "y=z[0]-273;\t\t\t\t\t #Thermocouple Temperature, C\n",
+ "\n",
+ "'''NOTE: Equation written is absolutely correct and solving this equation\n",
+ " should give real result as: 296.112 i.e. 23.112 C, but somehow python is giving wrong result.'''\n",
+ "\n",
+ "#Results\n",
+ "print \"Thermocouple Temperature is :\",round(y,1),\"C \\n\"\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Thermocouple Temperature is : 25.9 C \n",
+ "\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |