From df60071cf1d1c18822d34f943ab8f412a8946b69 Mon Sep 17 00:00:00 2001 From: hardythe1 Date: Wed, 3 Jun 2015 15:27:17 +0530 Subject: add books --- .../Chapter05_2.ipynb | 618 +++++++++++++++++++++ 1 file changed, 618 insertions(+) create mode 100755 Fluid_Mechanics-Fundamentals_&_Applications/Chapter05_2.ipynb (limited to 'Fluid_Mechanics-Fundamentals_&_Applications/Chapter05_2.ipynb') diff --git a/Fluid_Mechanics-Fundamentals_&_Applications/Chapter05_2.ipynb b/Fluid_Mechanics-Fundamentals_&_Applications/Chapter05_2.ipynb new file mode 100755 index 00000000..405f5ec7 --- /dev/null +++ b/Fluid_Mechanics-Fundamentals_&_Applications/Chapter05_2.ipynb @@ -0,0 +1,618 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:71d12a9b1f227d609ccf14dd1bc86e16c513b6cc011b9917a83e499f2febffe6" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "heading", + "level": 1, + "metadata": {}, + "source": [ + "Chapter 05:Mass, Bernoulli and Energy Equations" + ] + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-1, Page No:190" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Vairable Decleration\n", + "V=10#Volume of the bucket in Gal\n", + "r_in=1 #Radius of the hose in cm\n", + "r_e=0.4 #Radius of the hose at the nozzle exit in cm\n", + "t=50 #Time taken to fill the bucket in s\n", + "C_gl=3.7854 #Conversion factor gal to Lit\n", + "rho=1 #Denisty of water in kg/Lit\n", + "C_v=10**-3 #Conersion factor in m^3/lit\n", + "\n", + "#Calculations\n", + "\n", + "#Part (a)\n", + "V_dot=(V*C_gl)/t #Volume flow rate in Lit/s\n", + "m_dot=rho*V_dot #Mass flow rate in kg/s\n", + "\n", + "#Part(b)\n", + "A_e=pi*r_e**2*10**-4 #Cross-Sectional Area of the nozzle at exit in m^2\n", + "V_e=(V_dot*C_v)/A_e #Average Velocity of water at nozzle exit in m/s\n", + "\n", + "#Result\n", + "print \"The Volume Flow rate is\",round(V_dot,3),\"L/s and the mass flow rate is\",round(m_dot,3),\"kg/s\"\n", + "print \"The area of cross section at nozzle exit is\",round(A_e,5),\"m^2\"\n", + "print \"The Average Velocity of water is\",round(V_e,1),\"m/s\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The Volume Flow rate is 0.757 L/s and the mass flow rate is 0.757 kg/s\n", + "The area of cross section at nozzle exit is 5e-05 m^2\n", + "The Average Velocity of water is 15.1 m/s\n" + ] + } + ], + "prompt_number": 24 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-2, Page No:191" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "h_o=1.2 #Original Height in m\n", + "h_2=0.6 #Water level drop in m\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "D_tank=0.9 #Diameter of the tank in m\n", + "D_jet=0.013 #Diameter at the jet in m\n", + "\n", + "#Calculations\n", + "#After carrying out the theroetical calculations and integration we arrive to obtain\n", + "t_min=((h_o**0.5-h_2**0.5)/((g/2)**0.5))*((D_tank/D_jet)**2) #Time required to reach a level 0.6m in s\n", + "t=t_min/60 #Converting time from sec to min\n", + "\n", + "#Result\n", + "print \"The time it takes to half empty the tank is\",round(t,1),\"min\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The time it takes to half empty the tank is 11.6 min\n" + ] + } + ], + "prompt_number": 34 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-3, Page No:195" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "##Note:The symbols in the textbook are cumbersome to code hence a different one has been used in this coding\n", + "\n", + "#Variable Decleration\n", + "h=50 #Elevation difference in m\n", + "m_dot=5000 #Mass flow rate at which the water is to be supplied in kg/s\n", + "W_dot_out=1862 #Electric Power generated in kWh\n", + "n_generator=0.95 #Efficiency of the generator in fraction\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "C=10**-3 #Conversion factor in kJ/kg/m^2/s^2\n", + "\n", + "#Calculations\n", + "\n", + "#Part(a)\n", + "\n", + "#Calling e_mech_in-e_mech_out as del_e for convienence \n", + "del_e=g*h*C #Change in water's mechanical energy per unit mass in kJ/kg\n", + "delta_E_fluid=m_dot*del_e # Change in energy of the fluid in kW\n", + "\n", + "n_overall=W_dot_out/delta_E_fluid #Overall Efficiency in fraction\n", + "\n", + "#Part(b)\n", + "n_turbine_gen=n_overall/n_generator #Mechanical efficiency os the turbine in fraction\n", + "\n", + "#Part(c)\n", + "W_dot_shaft_out=n_turbine_gen*delta_E_fluid #Shaft power output in kW\n", + "\n", + "#Result\n", + "print \"The overall efficiency is\",round(n_overall,2)\n", + "print \"The mechanical efficiency of the turbine is\",round(n_turbine_gen,1)\n", + "print \"The shaft power output is\",round(W_dot_shaft_out,1),\"kW\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The overall efficiency is 0.76\n", + "The mechanical efficiency of the turbine is 0.8\n", + "The shaft power output is 1960.0 kW\n" + ] + } + ], + "prompt_number": 36 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-5, Page No:205" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "P1=400 #Pressure at upstream of the jet in kPa\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "rho=1000 #Density of water in kg/m^3\n", + "C1=1000 #Conversion factor in N/m^2.kPa\n", + "C2=1 #Conversion factor in kg.m/s^2.N\n", + "\n", + "#Calculations\n", + "#Applying the Bernoulli Equation\n", + "z2=(P1*C1*C2)/(rho*g) #maximun height the water jet reaches in m\n", + "\n", + "#Result\n", + "print \"The water jet rises up to\",round(z2,1),\"m\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The water jet rises up to 40.8 m\n" + ] + } + ], + "prompt_number": 39 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-6, Page No:206" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "h=5 #Height at which the water tank is filled in m\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "\n", + "#Calculations\n", + "z1=h #Decleration in terms of datum in m\n", + "#Applying the Bernoulli Equation\n", + "V2=(2*g*z1)**0.5 #Maximum velocity that the water jet can attain in m/s\n", + "\n", + "#Result\n", + "print \"The maximum velocity that the water jet can attain is\",round(V2,1),\"m/s\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The maximum velocity that the water jet can attain is 9.9 m/s\n" + ] + } + ], + "prompt_number": 40 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-7, Page No:207" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "P_atm=101.3 #Atmospheric pressure in kPa\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "rho=750 #Denisty of gasoline in kg/m^3\n", + "z1=0.75 #Location of point 2 w.r.t point 1\n", + "D=5*10**-3 #Diameter of the siphon pipe in m\n", + "V=4 #Volume of gasoline to be siphoned in Lit\n", + "z3=2.75 #Height of point 3 w.r.t to point 2 in m\n", + "C1=1 #conversion factor in N.s^2/kg.m\n", + "C2=10**-3 #Conversion factor in kPa.m^2/N\n", + "#Calculations\n", + "\n", + "#Part (a)\n", + "#Applying the Bernoulli Equation\n", + "V2=(2*g*z1)**0.5 #Velocity in m/s\n", + "A=(pi*D**2)/4 #Cross-Sectional Area in m^2\n", + "V_dot=V2*A*1000#Volume flow rate in L/s\n", + "delta_t=V/V_dot #Time required to siphon gasoline in s\n", + "\n", + "#Part(b)\n", + "#Applying Bernoulli Equations\n", + "P3=P_atm-(rho*g*z3*C1*C2) #Pressure at point 3 in kPa\n", + "\n", + "#result\n", + "print \"The time requires to siphon 4L gasoline is\",round(delta_t,1),\"s\"\n", + "print \"The pressure at point 3 is\",round(P3,1),\"kPa\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The time requires to siphon 4L gasoline is 53.1 s\n", + "The pressure at point 3 is 81.1 kPa\n" + ] + } + ], + "prompt_number": 45 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-8, Page No:208" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Declerations\n", + "g=9.81 #Acceleration due to Gravity in m/s^2\n", + "h3=0.12 #Difference in level in m\n", + "\n", + "#Calculations\n", + "#Applying Bernoulli Equations\n", + "V1=(2*g*h3)**0.5 #Velocity of Fluid in m/s\n", + "\n", + "#Result\n", + "print \"The velocity of fkuid is\",round(V1,2),\"m/s\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The velocity of fkuid is 1.53 m/s\n" + ] + } + ], + "prompt_number": 46 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-9, Page No:209" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "rho_hg=13600 #density of mercury in kg/m^3\n", + "rho_sw=1025 #density of sea-water in kg/m^3\n", + "rho_atm_air=1.2 #Density of air in kg/m^3\n", + "P_atm_air=762 #Atmospheric pressure 320km away from the eye in mm oh Hg\n", + "P_air=560 #Atmospheric pressure at the eye of the strom in mm og Hg\n", + "C=10**-3 #Conversion factor in m/mm\n", + "V_A=250 #Hurricane Wind Velocity in km/hr\n", + "C_k=1/3.6 #Conversion Factor from km/hr to m/s \n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "\n", + "#Calculations\n", + "\n", + "#part(a)\n", + "h3=(rho_hg*(P_atm_air-P_air)*C)/rho_sw #Pressure difference in m\n", + "\n", + "#Part(b)\n", + "#Applying Bernoulli Equations\n", + "h_air=(V_A**2*C_k**2)/(2*g) #Height of air column in m\n", + "rho_air=(P_air*rho_atm_air)/P_atm_air #Density of Air in the hurricane in kg/m^3\n", + "h_dynamic=(rho_air*h_air)/rho_sw #Sea-Water column equivalent to air-column in m\n", + "h2=h3+h_dynamic #Total storm surge at point 2 in m\n", + "\n", + "#Result\n", + "print \"The pressure difference between point's 1 and 3 in terms of sea-water column is\",round(h3,2),\"m\"\n", + "print \"The total Storm Surge at point2 is\",round(h2,2),\"m\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The pressure difference between point's 1 and 3 in terms of sea-water column is 2.68 m\n", + "The total Storm Surge at point2 is 2.89 m\n" + ] + } + ], + "prompt_number": 61 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-12, Page No:221" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "V_dot=50 #Volumetric Flow rate in L/s\n", + "rho=1 #Density of water \n", + "n_motor=0.9 #efficiency of the electric motor in fraction\n", + "W_dot_electric=15 #Power of the electric motor in kW\n", + "P2=300 #Absolute pressure at the outlet in kPa\n", + "P1=100 #Absolute pressure at the inlet in kPa\n", + "c=4.18 #Specific heat of water in kJ/kg C\n", + "#Calculations\n", + "\n", + "#Part(a)\n", + "m_dot=rho*V_dot #Mass flow rate in kg/s\n", + "W_dot_pump=n_motor*W_dot_electric #Mechanical shaft power delivered in kW\n", + "delta_E_dot_mech_fluid=(m_dot*((P2-P1)/rho))/1000 #Increase in mechanical energy in kW\n", + "n_pump=delta_E_dot_mech_fluid/W_dot_pump #Efficiency in fraction\n", + "\n", + "#part (b)\n", + "E_dot_loss=W_dot_pump-delta_E_dot_mech_fluid #Lost mechanical energy in kW\n", + "delta_T=(E_dot_loss)/(m_dot*c) #Temperature rise of water due to mechanical inefficiency in degree C\n", + "\n", + "#Result\n", + "print \"The Mechanical efficiency of the pump is\",round(n_pump,3)\n", + "print \"The temperature rise of water due to mechanical inefficiency is\",round(delta_T,3),\"Degree Centigrade\"\n", + "\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The Mechanical efficiency of the pump is 0.741\n", + "The temperature rise of water due to mechanical inefficiency is 0.017 Degree Centigrade\n" + ] + } + ], + "prompt_number": 14 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-13, Page No:222" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "V_dot=100 #Discharge through the power plant in m^3/s\n", + "rho=1000 #Density of water in kg/m^3\n", + "z1=120 #Elevation from which the water flows in m\n", + "h_l=35 #Elevation of point 2 in m\n", + "n_turbine_gen=0.8 #Overall efficiency of the generator in fraction\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "C=10**-3 #Conversion Factor\n", + "\n", + "#Calculations\n", + "m_dot=rho*V_dot #mass flow rate through the turbine in kg/s\n", + "\n", + "#Applying Bernoullis principle and taking point 2 as reference point z2=0\n", + "h_turbine=z1-h_l #extracted turbine head in m\n", + "W_dot_turbine=m_dot*g*h_turbine*C #Turbine Power in kW\n", + "W_dot_electric=C*n_turbine_gen*W_dot_turbine #Electrical Power Generated by the actual Unit in MW\n", + "\n", + "#Result\n", + "print \"The electrical Power generated is\",round(W_dot_electric,1),\"MW\"\n" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "The electrical Power generated is 66.7 MW\n" + ] + } + ], + "prompt_number": 15 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-14, Page No:223" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "void_fraction=0.5 #Void Fraction\n", + "l=12 #Dimension of the fan in cm\n", + "w=40 #Dimension of the fan in cm\n", + "h=40 #Dimension of the fan in cm\n", + "delta_t=1 #time in s\n", + "rho=1.2 #Ddensity of air in kg/m^3\n", + "D=0.05 #Diameter of opening in the case in m\n", + "alpha2=1.1 #kinetic correction factor\n", + "n_fan=0.3 #Efficiency of the fan-motor\n", + "#Calculations\n", + "#Part(a)\n", + "V=void_fraction*l*w*h #Volume in cm^3\n", + "V_dot=(V/delta_t)*10**-6 #Volumetric flow rate in m^3/s\n", + "m_dot=rho*V_dot #mass flow rate in kg/s\n", + "A=(pi*D**2)/4 #Area of the opening is the case in m^2\n", + "\n", + "#Notation has been changed to avoid conflict\n", + "Vel=V_dot/A #Velocity of the air thorught the opening in m/s\n", + "\n", + "#Applying Bernoullis principle\n", + "W_dot_fan=m_dot*alpha2*Vel**2*0.5 #Work done in W\n", + "W_dot_electric=W_dot_fan/n_fan #Electric Work done in W\n", + "\n", + "#Part(b)\n", + "#Applying Brnoullis principle\n", + "#Notation has been changed here\n", + "delta_P=(rho*W_dot_fan)/m_dot #Pressure rise across fan in Pa\n", + "\n", + "#Result\n", + "print \"Wattage of the fan to be purchased is\",round(W_dot_electric,4),\"W\"\n", + "print \"The pressure difference across the fan is\",round(delta_P,1),\"Pa\"" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Wattage of the fan to be purchased is 0.5049 W\n", + "The pressure difference across the fan is 15.8 Pa\n" + ] + } + ], + "prompt_number": 18 + }, + { + "cell_type": "heading", + "level": 2, + "metadata": {}, + "source": [ + "Example 5.5-15, Page No:225" + ] + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "import math\n", + "\n", + "#Variable Decleration\n", + "W_shaft=5 #Shaft Power in kW\n", + "n_pump=0.72 #Efficiency of the pump in fraction\n", + "g=9.81 #Acceleration due to gravity in m/s^2\n", + "h_l=4 #Head loss in m\n", + "z2=25 #Datum in m\n", + "rho=1000 #Density of water in kg/m^3\n", + "\n", + "#Calculations\n", + "W_dot_pump=n_pump*W_shaft #Useful mechanical power returned in kW\n", + "\n", + "#Applying Bernoullis Principle\n", + "m_dot=(W_dot_pump/(g*(z2+h_l)))*1000 #mass floe rate in kg/s\n", + "V_dot=(m_dot/rho) #Volumetric flow rate in m^3/s\n", + "delta_P=W_dot_pump/V_dot #Pressure difference in kPa\n", + "\n", + "#Result\n", + "print \"Discharge of water is\",round(V_dot,4),\"m^3/s\"\n", + "print \"The pressure difference across the pump is\",round(delta_P),\"kPa\"\n", + "#Answer in the coding is off by 1 kPa due to decimal point accuracy" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Discharge of water is 0.0127 m^3/s\n", + "The pressure difference across the pump is 284.0 kPa\n" + ] + } + ], + "prompt_number": 21 + } + ], + "metadata": {} + } + ] +} \ No newline at end of file -- cgit