summaryrefslogtreecommitdiff
path: root/sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb
diff options
context:
space:
mode:
authorkinitrupti2017-05-12 18:53:46 +0530
committerkinitrupti2017-05-12 18:53:46 +0530
commit6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d (patch)
tree22789c9dbe468dae6697dcd12d8e97de4bcf94a2 /sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb
parentd36fc3b8f88cc3108ffff6151e376b619b9abb01 (diff)
downloadPython-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.gz
Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.tar.bz2
Python-Textbook-Companions-6279fa19ac6e2a4087df2e6fe985430ecc2c2d5d.zip
Removed duplicates
Diffstat (limited to 'sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb')
-rw-r--r--sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb1396
1 files changed, 0 insertions, 1396 deletions
diff --git a/sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb b/sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb
deleted file mode 100644
index 0cef27c9..00000000
--- a/sample_notebooks/RahulJoshi/Chapter_1_An_Overview_of_Heat.ipynb
+++ /dev/null
@@ -1,1396 +0,0 @@
-{
- "cells": [
- {
- "cell_type": "markdown",
- "metadata": {
- "collapsed": true
- },
- "source": [
- "# Example 1.1 Page Number 2"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 3,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Heat Flow through the surface is 17850.0 W\n",
- "Temprature Gradient in flow direction -700.0 C/m\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T = 100 # temperature of wall 1 in deg celcius\n",
- "\n",
- "t = 30 # temperature of wall 2 in deg celcius\n",
- "\n",
- "L = 0.1 # distance between the walls in meters\n",
- "\n",
- "k = 8.5 # thermal conductivity in W/mK\n",
- "\n",
- "A = 3 # area is meters square\n",
- "\n",
- "#calculation\n",
- "\n",
- "Q = (T-t)/(L/(k*A)) # heat flow rate in (W)\n",
- "\n",
- "tempgrad = (-1*Q)/(k*A) # temperature gradient in celcius/meter\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Heat Flow through the surface is\",Q,\"W\")\n",
- "\n",
- "print(\"Temprature Gradient in flow direction\",tempgrad,\"C/m\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Example 1.2 Page Number 6"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 6,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Convective heat transfer rate 1500.0 W\n",
- "Resistance 0.08 C/W\n",
- "Temprature Gradient along y direction -3000.0 C/m\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T = 160 # temperature of wall 1 in deg celcius\n",
- "\n",
- "t = 40 # temperature of wall 2 in deg celcius\n",
- "\n",
- "k = 1 # thermal conductivity in W/mK\n",
- "\n",
- "h = 25 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "A = 0.5 # area is meters square\n",
- "\n",
- "#calculation\n",
- "\n",
- "Q = h*A*(T-t) # heat tranfer by convection (W)\n",
- "\n",
- "r = 1/(h*A) # resistance (C/W)\n",
- "\n",
- "tempgrad = (-1*Q)/(k*A) # temperature gradient in celcius/meter along y\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Convective heat transfer rate \",Q,\"W\")\n",
- "\n",
- "print(\"Resistance\",r,\"C/W\")\n",
- "\n",
- "print(\"Temprature Gradient along y direction\",tempgrad,\"C/m\")\n"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Example 1.3 Page number 7"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 11,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Heat Flow through the surface is 2171.37 W\n",
- "Resistance 0.0783 K/W\n",
- "Equivalent thermal coefficient 6.3864 W/m2K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T = 473 # temperature of wall 1 kelvin\n",
- "\n",
- "t = 303 # temperature of wall 2 in kelvin\n",
- "\n",
- "sigma = 5.67*10**-8 # Stefen-Boltzmann constant\n",
- "\n",
- "F = 0.46 # emmissivity \n",
- "\n",
- "A = 2 # area is meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "Q = F*sigma*A*(T**4-t**4) # heat exchange in (W)\n",
- "\n",
- "R = (T-t)/Q # Resistance in (K/W)\n",
- "\n",
- "hr = 1/(R*A) # equivalent thermal coefficient W/m2K\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Heat Flow through the surface is\",round(Q,2),\"W\")\n",
- "\n",
- "print(\"Resistance\",round(R,4),\"K/W\")\n",
- "\n",
- "print(\"Equivalent thermal coefficient\",round(hr,4),\"W/m2K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Example 1.4 Page number 8"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 17,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Heat Received 7092.23 W\n",
- "T2 = 368.479 K\n",
- "Temprature on other side of the wall 263.3 K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "Tinf1 = 500 # temperature of wall 1 Kelvin\n",
- "\n",
- "T1 = 400 # temperature of wall 2 in Kelvin\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 50 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "k = 45 # thermal conductivity in W/mK\n",
- "\n",
- "L = 0.2 # slab thickness in meters\n",
- "\n",
- "#calculation\n",
- "\n",
- "Q = sigma*((Tinf1/100)**4 - (T1/100)**4)+ h*(Tinf1-T1) # heat received (W)\n",
- "\n",
- "dT = Q*(L/k) # temp gradient (K)\n",
- "\n",
- "T2 = T1-dT #\n",
- "\n",
- "Tinf2 = 263.3 # temperature on the other side of the wall using trial and error\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Heat Received \",round(Q,3),\"W\")\n",
- "\n",
- "print(\"T2 = \",round(T2,3),\"K\")\n",
- "\n",
- "print(\"Temprature on other side of the wall\",round(Tinf2,3),\"K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Example 1.5 Page number "
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 26,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Rate of change of temperature 0.03984 C/s\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T1 = 400 # temperature of wall 1 Kelvin\n",
- "\n",
- "T2 = 100 # temperature of wall 2 in Kelvin\n",
- "\n",
- "sigma = 5.67*10**-8 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 200 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "q = 1.5*10**6 # heat generated in W/m3\n",
- "\n",
- "H = 0.3 # height in meters\n",
- "\n",
- "r = 0.15 # radius in meters\n",
- "\n",
- "rho = 19000 # density in kg/m3\n",
- "\n",
- "cp = 118 # specific heat capacity in kJ/kgK\n",
- "\n",
- "#calculation\n",
- "\n",
- "Sa = 2*3.14*r*H+2*3.14*r**2 # Surface area in meters sq\n",
- "\n",
- "Hc = 3.14*r**2*H*rho*cp # heat capacity in J/deg C\n",
- "\n",
- "Hg = 3.14*r**2*H*q # Heat generated in W\n",
- "\n",
- "Hcon = h*Sa*(T1-T2) # convective heat transfer in W\n",
- "\n",
- "Hrad = sigma*Sa*((T1+273)**4 - (T2+273)**4)\n",
- "\n",
- "Th = Hg-Hcon-Hrad\n",
- "\n",
- "dTbydt = Th/Hc\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Rate of change of temperature \",round(dTbydt,5),\"C/s\")\n"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Example 1.6 Page number 11"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 40,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "BTU/hrftF = 1.7322 W/mC\n",
- "BTU/hrft2F = 5.6831 W/m2C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "J = 9.47*10**-4 # Joule to BTU conversion\n",
- "\n",
- "m = 39.37 # meter to inch conversion\n",
- "\n",
- "kg = 2.2046 # kg to lb conversion\n",
- "\n",
- "C = 9/5 # Celcius to Farhenight\n",
- "\n",
- "# Calculation\n",
- "\n",
- "BTU = 1/J # in Joule\n",
- "\n",
- "ft = 12/m # in feet\n",
- "\n",
- "a = (BTU/(3600*ft*(5/9))) # in BTU/hrftF\n",
- "\n",
- "print(\"BTU/hrftF = \",round(a,4),\"W/mC\")\n",
- "\n",
- "b = (BTU/(3600*ft**2*(5/9))) # in BTU/hrftF\n",
- "\n",
- "print(\"BTU/hrft2F = \",round(b,4),\"W/m2C\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 1 page Number 11"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 43,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Temprature gradient along surface -1111.7 C/m\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T1 = 200 # temperature of wall 1 Kelvin\n",
- "\n",
- "T2 = 60 # temperature of wall 2 in Kelvin\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 80 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "k = 12 # thermal conductivity in W/mK\n",
- "\n",
- "#L = 0.2 # slab thickness in meters\n",
- "\n",
- "#calculation\n",
- "\n",
- "Q = sigma*(((T1+273)/100)**4 - ((T2+273)/100)**4)+ h*(T1-T2) # heat received (W)\n",
- "\n",
- "dTbydx = Q/(-1*k) # temp gradient (K)\n",
- "\n",
- "print(\"Temprature gradient along surface\",round(dTbydx,1),\"C/m\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 2 Page number 12"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 46,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Temprature only conduction and convection 682.174 C\n",
- "Temprature only conduction and radiation 1139.148 C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "dTbydx = -9000 # temperature gradient \n",
- "\n",
- "T2 = 30 # temperature of wall 2 in C\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "k = -25 # Convective heat transfer coefficient W/mK\n",
- "\n",
- "h = 345 # thermal conductivity in W/m2K\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "#only conduction and convection\n",
- "\n",
- "T11 = k*A*dTbydx/(h*A) + T2\n",
- "\n",
- "#only conduction and radiation\n",
- "\n",
- "T12 = (((k*A*dTbydx/(sigma)) + ((T2+273)/100)**4)*100**4)**(1/4)-273\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Temprature only conduction and convection \",round(T11,3),\"C\")\n",
- "\n",
- "print(\"Temprature only conduction and radiation \",round(T12,3),\"C\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem number 3 Page number 12"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 47,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Wall surface temperature 330.4 K\n",
- "Heat Generated 2252.765 W/m2\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "Qc = 2250 # heat conducted in W/m2\n",
- "\n",
- "T1 = 303 # temperature of wall 2 in C\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 75 # thermal conductivity in W/m2K\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "# taking the approximate value from table 330.4\n",
- "\n",
- "Tapprox = 330.4\n",
- "\n",
- "Q = h*(Tapprox-T1)+sigma*((Tapprox/100)**4-(T1/100)**4)\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Wall surface temperature \",round(Tapprox,3),\"K\")\n",
- "\n",
- "print(\"Heat Generated \",round(Q,3),\"W/m2\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 4 page number 13"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 48,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Wall surface temperature 277.75 K\n",
- "Heat Generated 65.479 W\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "Hc = 65.5 # heat conducted in W/m\n",
- "\n",
- "T1 = 263 # temperature of wall in K\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 4.35 # thermal conductivity in W/m2K\n",
- "\n",
- "r = 0.08 # area in meters \n",
- "\n",
- "#calculation\n",
- "\n",
- "# taking the approximate value from table 277.75 K\n",
- "\n",
- "Tapprox = 277.75\n",
- "\n",
- "Q = h*3.14*r*2*(Tapprox-T1)+sigma*2*3.14*r*((Tapprox/100)**4-(T1/100)**4)\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Wall surface temperature \",round(Tapprox,3),\"K\")\n",
- "\n",
- "print(\"Heat Generated \",round(Q,3),\"W\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 5 Page number 14"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 51,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Wall surface temperature 386.1 K\n",
- "Heat Generated 449.65 W\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "Hc = 450 # heat conducted in W/m\n",
- "\n",
- "T1 = 396.4 # temperature of wall in K\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 1.5 # thermal conductivity in W/m2K\n",
- "\n",
- "r = 0.08 # area in meters \n",
- "\n",
- "A = 4*3.14*0.48**2 # area in meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "# taking the approximate value from table 386.1 K\n",
- "\n",
- "Tapprox = 386.1\n",
- "\n",
- "Q = h*A*(T1-Tapprox)+sigma*A*((T1/100)**4-(Tapprox/100)**4)\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Wall surface temperature \",round(Tapprox,3),\"K\")\n",
- "\n",
- "print(\"Heat Generated \",round(Q,3),\"W\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 6 Page number 14"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 52,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Heat Capacity 1000.0 J/C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "dTbydt = 0.5 # Temperature transition in C/s\n",
- "\n",
- "Qr = 4000 # Heat Received in J/s\n",
- "\n",
- "Qc = 5200 # Heat Convection in J/s\n",
- "\n",
- "qdot = 1700 # Heat generated in J/s\n",
- "\n",
- "#calculation\n",
- "\n",
- "HeatCapacity = (Qr-Qc+qdot)/dTbydt\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Heat Capacity \",round(HeatCapacity,3),\"J/C\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 7 page number 15"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 54,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Time rate of temperature change 0.1 C/s\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "Q = 240 # Heat Received in J/s\n",
- "\n",
- "qdot = 100000 # Heat generated in J/m3/s\n",
- "\n",
- "rho = 2500 # density in kg/m3\n",
- "\n",
- "cp = 0.52*10**3 # heat capacity in kJ/KgK\n",
- "\n",
- "a = 0.2 # side of the cube in meters\n",
- "\n",
- "#calculation\n",
- "\n",
- "V = a**3\n",
- "\n",
- "dTbydt = (Q+qdot*V)/(rho*V*cp)\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Time rate of temperature change \",round(dTbydt,3),\"C/s\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 8 Page Number 15 "
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 55,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Heat Convected 1309.5 W/m2\n",
- "Heat Received 1303.428 W/m2\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T1 = 160 # heat conducted in W/m\n",
- "\n",
- "T2 = 30 # temperature of wall in K\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 45 # thermal conductivity in W/m2K\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "# taking the approximate value from table 332 K\n",
- "\n",
- "Tapprox = 332.1 \n",
- "\n",
- "Hc = h*A*(Tapprox-(273+T2)) # Heat Convected in W/m2\n",
- "\n",
- "Hr = sigma*A*(((T1+273)/100)**4-(Tapprox/100)**4) # Heat received in W/m2\n",
- "\n",
- "# Result\n",
- "\n",
- "print(\"Heat Convected \",round(Hc,3),\"W/m2\")\n",
- "\n",
- "print(\"Heat Received \",round(Hr,3),\"W/m2\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 9 Page number 16"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 60,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Total Heat loss case 1 150.6 W\n",
- "Total Heat loss case 2 81.9 W\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T1 = 37 # temperature of body in C\n",
- "\n",
- "T21 = 26 # temperature of air in C\n",
- "\n",
- "T22 = 5 # temperature of walls in room in C\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 6 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "A = 0.6 # area in meters sq\n",
- "\n",
- "#calculation\n",
- "\n",
- "Hc = h*A*(T1-T21) # Heat Convected in W/m2\n",
- "\n",
- "Hr1 = sigma*A*(((T1+273)/100)**4-((T22+273)/100)**4) # Heat received in W/m2\n",
- "\n",
- "Ht1 = Hr1+Hc\n",
- "\n",
- "# calculate when temperature is 26C\n",
- "\n",
- "Hr2 = sigma*A*(((T1+273)/100)**4 - ((T21+273)/100)**4) # Heat received in W/m2\n",
- "\n",
- "Ht2 = Hc+Hr2\n",
- "\n",
- "print(\"Total Heat loss case 1\",round(Ht1,1),\"W\")\n",
- "\n",
- "print(\"Total Heat loss case 2\",round(Ht2,1),\"W\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 10 Page number 16"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 63,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Net heat Gain 291.1 W\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "#variable declaration\n",
- "\n",
- "T1 = 37 # temperature of body in C\n",
- "\n",
- "T21 = 650 # temperature of air in C\n",
- "\n",
- "T22 = 5 # temperature of walls in room in C\n",
- "\n",
- "sigma = 5.67 # Stefen-Boltzmann constant\n",
- "\n",
- "h = 6 # Convective heat transfer coefficient W/m2K\n",
- "\n",
- "A = 0.6 # area in meters sq\n",
- "\n",
- "F = 0.01 # fraction of radiation \n",
- "\n",
- "#calculation\n",
- "\n",
- "# Heat loss by convection in W\n",
- "\n",
- "Hc = h*A*(T1-T22)\n",
- "\n",
- "# Heat Gain by radiation\n",
- "\n",
- "Hr = sigma*F*(((T21+273)/100)**4 - ((T1+273)/100)**4) # Heat received in W/m2\n",
- "\n",
- "Hnet = Hr-Hc\n",
- "\n",
- "print(\"Net heat Gain\",round(Hnet,1),\"W\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 11 Page number 16"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 66,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Equilibrium Temperature = 960.01 K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "Q = 1500 # heat dissipation in W\n",
- "\n",
- "sigma = 5.67 # stefan-Boltzmann constant\n",
- "\n",
- "T2 = 288 # temperature in K\n",
- "\n",
- "r = 0.04 # radius in meters\n",
- "\n",
- "H = 0.25 # height in meters\n",
- "\n",
- "T1 = ((Q/(sigma*3.14*r*H)+(288/100)**4)*100**4)**(1/4)\n",
- "\n",
- "print(\"Equilibrium Temperature = \",round(T1,2),\"K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem number 12 Page number 17"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 68,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Equilibrium Temperature = 62.0 C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable declaration\n",
- "\n",
- "Hr = 800 # Heat Rate in W/m2\n",
- "\n",
- "h1 = 10 # convective heat transfer rate on back of plate in W/m2K\n",
- "\n",
- "h2 = 15 # convective heat transfer rate on front of plate in W/m2K\n",
- "\n",
- "T2 = 30 # temperature on both sides of the plate\n",
- "\n",
- "T = (Hr+h1*30+h2*30)/25\n",
- "\n",
- "print(\"Equilibrium Temperature = \",round(T,2),\"C\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem number 13 Page number 17"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 73,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "temperature of the plate 784.57 K\n",
- "heat transfer with sheet 19668.31 W\n",
- "heat transfer without sheet 39336.61 W\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "T1 = 650 # temperature from one side of the source in C\n",
- "\n",
- "T2 = 150 # temperature on other side of the surface in C\n",
- "\n",
- "sigma = 5.67 # stefan-boltzmann constant\n",
- "\n",
- "T = (((((T1+273)/100)**4 + ((T2+273)/100)**4)/2)*100**4)**(1/4)\n",
- "\n",
- "print(\"temperature of the plate\",round(T,2),\"K\")\n",
- "\n",
- "Q1 = sigma*(((T1+273)/100)**4 - (T/100)**4)\n",
- "\n",
- "print(\"heat transfer with sheet\",round(Q1,2),\"W\")\n",
- "\n",
- "Q2 = sigma*(((T1+273)/100)**4 - ((T2+273)/100)**4)\n",
- "\n",
- "print(\"heat transfer without sheet\",round(Q2,2),\"W\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 14 Page number 18"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 83,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "COnvective heat trasnfer rate 375.0 W/m2K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "Tair = 120 # temperature of air in C\n",
- "\n",
- "T1 = 42 # temperature of plate 1 in C\n",
- "\n",
- "T2 = 30 # temperature of plate 2 in C\n",
- "\n",
- "L = 0.01 # length of the slab in meters\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "k = 22.5 # thermal conductivity in W/mK\n",
- "\n",
- "# Calculation\n",
- "\n",
- "Q = (T1-T2)/(L/(k*A))\n",
- "\n",
- "Tnew = T1+6\n",
- "\n",
- "h = Q/(A*(Tair-Tnew))\n",
- "\n",
- "print(\"COnvective heat trasnfer rate\",round(h,2),\"W/m2K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 15 Page number 19"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 87,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Efficiency of the collector when temperature is 32 deg C 47.5 %\n",
- "Efficiency of the collector when temperature is 45 deg C 75.62 %\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable declaration\n",
- "\n",
- "T1 = 60 # temperature of the tube in C\n",
- "\n",
- "T2 = 32 # temperature of air in C\n",
- "\n",
- "h = 15 # convective heat transfer in W/m2K\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "Qf = 800 # heat flux in W/m2\n",
- "\n",
- "Tnew = 45 # new temperature in C\n",
- "\n",
- "# Calculation\n",
- "\n",
- "Q = h*A*(T1-T2) # heat transfer in W\n",
- "\n",
- "eff = ((Qf-Q)/Qf)*100\n",
- "\n",
- "print(\"Efficiency of the collector when temperature is 32 deg C\",round(eff,2),\"%\")\n",
- "\n",
- "# Heat lost by convection when T = 45 C\n",
- "\n",
- "Q2 = h*A*(Tnew-T2) # heat transfer in W\n",
- "\n",
- "eff1 = ((Qf-Q2)/Qf)*100 \n",
- "\n",
- "print(\"Efficiency of the collector when temperature is 45 deg C\",round(eff1,2),\"%\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 16 Page Number 19"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 88,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Temperature of the air 428.89 C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable declaration\n",
- "\n",
- "Tg = 40 # temperature of the glass plate in C\n",
- "\n",
- "dT = 5 # temperature graditent in C\n",
- "\n",
- "L = 0.001 # length in meters\n",
- "\n",
- "k = 1.4 # conductive heat transfer coefficient in W/mK\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "h = 18 # convective heat trasnfer coefficient in W/m2K\n",
- "\n",
- "# Calculation\n",
- "\n",
- "Q = dT/(L/(k*A))\n",
- "\n",
- "Tair = (Q/h)+ Tg\n",
- "\n",
- "print(\"Temperature of the air\",round(Tair,2),\"C\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 17 Page number 20"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 89,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Temperature gradient in the solid -631.58 C/m\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "h = 30 # Convective heat transfer coefficient in W/m2K\n",
- "\n",
- "k = 9.5 # conductive heat trasnfer coefficient in W/mK\n",
- "\n",
- "T1 = 260 # temperature of the surface in C\n",
- "\n",
- "T2 = 60 # temperature of the air in C\n",
- "\n",
- "tempgrad = (h/k)*(T2-T1)\n",
- "\n",
- "print(\"Temperature gradient in the solid\",round(tempgrad,2),\"C/m\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 18 Page number 20"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 90,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Steady state temperature of plate 313.33 K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "h1 = 100 # convective heat transfer coeffcient in W/m2K\n",
- "\n",
- "h2 = 15 # convective heat transfer coeffcient in W/m2K\n",
- "\n",
- "# solving by trial and error we get T1 = 313.33 K\n",
- "\n",
- "T1 = 313.33 \n",
- "\n",
- "print(\"Steady state temperature of plate\",round(T1,2),\"K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 19 Page number 21"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 91,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Surface temperature 674.39 K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "k = -22.5 # conductive heat trasnfer coefficient in W/mK\n",
- "\n",
- "tempgrad = -500 # temperature gradient in C/m\n",
- "\n",
- "sigma = 5.67 # stefan-boltzmann constant\n",
- "\n",
- "Ts = 303 # temperatre of surroundings in K\n",
- "\n",
- "# Calculation\n",
- "\n",
- "T2 = ((((k*tempgrad)/sigma)+(Ts/100)**4)*100**4)**(1/4)\n",
- "\n",
- "print(\"Surface temperature\",round(T2,2),\"K\")\n",
- "\n",
- "\n"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 20 Page number 21"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 92,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Surface temperatrue 230.2 K\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# variable declaration \n",
- "\n",
- "Hc = 2000 # heat generated in W\n",
- "\n",
- "r = 1 # radius in meters\n",
- "\n",
- "sigma = 5.67 # stefan boltzmann constant\n",
- "\n",
- "T2 = 0 # temperate of space in K\n",
- "\n",
- "# Calculation \n",
- "\n",
- "T = ((Hc/(4*3.14*r**2*sigma))*100**4)**(1/4)\n",
- "\n",
- "print(\"Surface temperatrue\",round(T,2),\"K\")"
- ]
- },
- {
- "cell_type": "markdown",
- "metadata": {},
- "source": [
- "# Problem 21 page Number 22"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 93,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Temperature drop through the wall 1.33 C\n"
- ]
- }
- ],
- "source": [
- "import math\n",
- "\n",
- "# Variable Declaration\n",
- "\n",
- "Q = 10 # heat flux in W/m2\n",
- "\n",
- "A = 1 # area in meters sq\n",
- "\n",
- "k = 1.5 # thermal cnductivity in W/mK\n",
- "\n",
- "t = 0.2 # wall thockness in m\n",
- "\n",
- "# Calculation\n",
- "\n",
- "dT = (Q*t)/(k*A)\n",
- "\n",
- "print(\"Temperature drop through the wall\",round(dT,2),\"C\")"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": null,
- "metadata": {
- "collapsed": true
- },
- "outputs": [],
- "source": []
- }
- ],
- "metadata": {
- "anaconda-cloud": {},
- "kernelspec": {
- "display_name": "Python [default]",
- "language": "python",
- "name": "python3"
- },
- "language_info": {
- "codemirror_mode": {
- "name": "ipython",
- "version": 3
- },
- "file_extension": ".py",
- "mimetype": "text/x-python",
- "name": "python",
- "nbconvert_exporter": "python",
- "pygments_lexer": "ipython3",
- "version": "3.5.2"
- }
- },
- "nbformat": 4,
- "nbformat_minor": 1
-}