diff options
Diffstat (limited to 'Engineering_Thermodynamics/ch1.ipynb')
| -rwxr-xr-x | Engineering_Thermodynamics/ch1.ipynb | 249 |
1 files changed, 249 insertions, 0 deletions
diff --git a/Engineering_Thermodynamics/ch1.ipynb b/Engineering_Thermodynamics/ch1.ipynb new file mode 100755 index 00000000..fff2e7be --- /dev/null +++ b/Engineering_Thermodynamics/ch1.ipynb @@ -0,0 +1,249 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:0acf769d31306f5bc291e2574308a74de8651b4a4723b52ee59419306e6b7112" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 1 : Introduction" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.1 Page No : 20" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "d_r = 13640.; \t\t\t# Density of mercury in kg/m3\n", + "g = 9.79; \t\t\t # Acceleration due to gravity in m/s2\n", + "z = 562e-03; \t\t\t# Difference in height in m\n", + "z0 = 761e-03; \t\t\t# Reading of barometer in m\n", + "\n", + "# Calculation\n", + "P = (d_r*g*(z+z0))*(0.987/1e05); \t\t\t# Gas Pressure in bar\n", + "\n", + "# Results\n", + "print \"Gas Pressure is %.3f bar\"%P\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Gas Pressure is 1.744 bar\n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.2 Page No : 21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "d_r = 13.6e03; \t\t\t# Density of mercury in kg/m3\n", + "g = 9.81; \t\t\t# Acceleration due to gravity in m/s2\n", + "z = 710e-03; \t\t\t# Stean flow pressure in m\n", + "z0 = 772e-03; \t\t\t# Reading of barometer in m\n", + "P = 1.4e06; \t\t\t# Gauge pressure of applied steam in Pa\n", + "\n", + "# Calculation\n", + "P0 = d_r*g*z0; \t\t\t# Atmospheric pressure in Pa\n", + "Pi = P+P0 ; \t\t\t# Inlet steam pressure in Pa\n", + "Pc = d_r*g*(z0-z); \t\t\t# Condenser pressure in Pa\n", + "\n", + "# Results\n", + "print \"Inlet steam pressure is %.3e Pa\"%Pi\n", + "print \"Condenser pressure is %.2e Pa\"%Pc\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Inlet steam pressure is 1.503e+06 Pa\n", + "Condenser pressure is 8.27e+03 Pa\n" + ] + } + ], + "prompt_number": 6 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.3 Page No : 21" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "z = 0.760; \t\t\t# Barometer reading in m\n", + "# Part (a)\n", + "h1 = 40e-02; \t\t\t# Mercury height in vaccume in m\n", + "d_r = 13.6e03; \t\t\t# Density of mercury in kg/m3\n", + "g = 9.80; \t\t\t# Acceleration due to gravity in m/s2\n", + "\n", + "# Calculation and Results\n", + "Patm = z*d_r*g; \t\t\t# Atmospheric pressure in Pas\n", + "Pv = h1*d_r*g; \t\t\t# Pressue in vaccume in Pa\n", + "Pabst = Patm-Pv; \t\t\t# Absolute pressure in Pa\n", + "Pabs = 101.325 - Pv/1000\n", + "\n", + "print \"Pvaccum = %.2f kPa\"%(Pv/1000)\n", + "print \"Pabsolute = %.f\"%(Pabst/1000),\"kPa\"\n", + "print \"40 cmHg vacuum : %.3f kPa\"%Pabs\n", + "\n", + "# Part (b)\n", + "h2 = 90e-02; \t\t\t# Mercury height in gauge in m\n", + "Pg = h2*d_r*g; \t\t\t# Gauge Pressure in Pa\n", + "Pabs1 = Patm + Pg ; \t\t\t# Absolute pressure in Pa\n", + "\n", + "print \"\\nPgauge = %.f kPa\"%(Pg/1000)\n", + "print \"90cmHg gauge is %.3f\"%(Pabs1/1000),\"kPa\"\n", + "\n", + "# Part(c)\n", + "d_w = 1e03 ; \t\t\t# Density of water in kg/m3\n", + "h3 = 1.2 ; \t\t\t# Gauge Pressure water height in m\n", + "Pga = d_w*h3*g; \t\t\t# Gauge Pressure in Pa\n", + "Pabs3 = Patm + Pga ; \t\t\t# Absolute pressure in Pa\n", + "print \"\\n1.2 m H2O gauge is %.3f\"%(Pabs3/1000),\"kPa\"\n", + "\n", + "# rounding off error" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Pvaccum = 53.31 kPa\n", + "Pabsolute = 48 kPa\n", + "40 cmHg vacuum : 48.013 kPa\n", + "\n", + "Pgauge = 120 kPa\n", + "90cmHg gauge is 221.245 kPa\n", + "\n", + "1.2 m H2O gauge is 113.053 kPa\n" + ] + } + ], + "prompt_number": 3 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.4 Page No : 22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from scipy.integrate import quad\n", + "\n", + "# Variables\n", + "Pr = 1.033e05; \t\t\t# Required Pressure in bar\n", + "\n", + "# Calculation\n", + "def pressure(p):\n", + " return p**(-0.714);\n", + "g = 9.81; \t\t\t# Acceleration due to gravity in m/s2\n", + "H = ((2.5e05**0.714)/g)* quad(pressure,0,Pr)[0]; \t\t\t# Depth of atmosphere required in m\n", + "\n", + "# Results\n", + "print \"The depth of atmosphere required is %.3f Km\"%(H/1000)\n", + " \n", + "# note : there will be rounding off error because of quad function" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The depth of atmosphere required is 69.203 Km\n" + ] + } + ], + "prompt_number": 23 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 1.5 Page No : 22" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "\n", + "# Variables\n", + "m = 68. ; \t\t\t# Astronaut mass in Kg\n", + "g = 9.806; \t\t\t# Acceleration due to gravity in m/s2\n", + "\n", + "# Calculation\n", + "a = 10*g ; \t\t\t# Lift off acceleration in m/s2\n", + "F = m*a; \t\t\t# Net vertical force in N\n", + "\n", + "# Results\n", + "print \"Net vertical force experienced by astronaut is %.0f N\"%F\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Net vertical force experienced by astronaut is 6668 N\n" + ] + } + ], + "prompt_number": 6 + } + ], + "metadata": {} + } + ] +}
\ No newline at end of file |