summaryrefslogtreecommitdiff
path: root/Basic_Fluid_Mechanics_by_Peerless
diff options
context:
space:
mode:
Diffstat (limited to 'Basic_Fluid_Mechanics_by_Peerless')
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/10-External_flows.ipynb256
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/11-elementary_analysis.ipynb338
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/12-TURBOMACHINES.ipynb222
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/13-Hydraulic_power_transmission.ipynb152
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/14-Further_Devolopments.ipynb111
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/2-Similarity.ipynb244
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/3-Dimensional_Analaysis.ipynb132
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/5-Control_Volume_Analysis.ipynb170
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/6-Steady_one_dimensional_reversible_flow.ipynb311
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/7-Steady_one_dimensional_Irreversible_flow.ipynb402
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/8-analysis_of_dimensional_constant_density_laminar_flow.ipynb125
-rw-r--r--Basic_Fluid_Mechanics_by_Peerless/9-analysis_of_dimensional_constant_density_turbulent_flow.ipynb109
12 files changed, 2572 insertions, 0 deletions
diff --git a/Basic_Fluid_Mechanics_by_Peerless/10-External_flows.ipynb b/Basic_Fluid_Mechanics_by_Peerless/10-External_flows.ipynb
new file mode 100644
index 0000000..84141e4
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/10-External_flows.ipynb
@@ -0,0 +1,256 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 10: External flows"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"g= 32.2 //ft/sec^2\n",
+"u= 3.6*10^-5 //lbf sec/ft^2\n",
+"d= 64 //lbm/ft^2\n",
+"l= 20 //ft\n",
+"a= 0.5\n",
+"//CALCULATIONS\n",
+"sw= u*g/(a*d)\n",
+"sw1= u^2*g*l/(2*a*d)\n",
+"Re=[1 2 3 4 5 6 7 8 9 10]*10^5\n",
+"Vinf=Re*u*g/(d*a)\n",
+"Cd=[1.2 1.15 0.94 0.68 0.305 0.31 0.32 0.33 0.34 0.35]\n",
+"cdre=Cd.*Re^2\n",
+"D=sw1*cdre\n",
+"//RESULTS\n",
+"printf ('velocity = %.2e ft/sec',sw)\n",
+"printf ('\n Force = %.2e lbf',sw1)\n",
+"disp(Vinf)\n",
+"disp(D)\n",
+""
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.2: ex2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"g= 32.2 //ft/sec^2\n",
+"u= 3.6*10^-5 //lbf sec/ft^2\n",
+"d= 64 //lbm/ft^2\n",
+"l= 20 //ft\n",
+"a= 0.5\n",
+"//CALCULATIONS\n",
+"sw= u*g/(a*d)\n",
+"sw1= u^2*g*l/(2*a*d)\n",
+"Re=[1 2 3 4 5 6 7 8 9 10]*10^5\n",
+"Vinf=Re*u*g/(d*a)\n",
+"Cd=[1.2 1.15 0.94 0.68 0.305 0.31 0.32 0.33 0.34 0.35]\n",
+"cdre=Cd.*Re^2\n",
+"D=sw1*cdre\n",
+"//RESULTS\n",
+"plot(Vinf,D)\n",
+"xtitle('','Vinf, ft/sec', 'D, lbf') \n",
+"\n",
+"//data for curves b,c,d is not given"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"v1= 10 //ft/sec\n",
+"v2= 9 //ft/sec\n",
+"a= 1.02\n",
+"r= 5.95\n",
+"//CALCULATIONS\n",
+"ca= (v1/v2)^2\n",
+"Cd= r*(ca-1+2-2*ca)+2*a*ca\n",
+"//RESULTS\n",
+"printf ('Drage coeffcieicnt = %.2f',Cd)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"A= 320 //ft/^2\n",
+"w= 18000 //lbf\n",
+"v= 230 //ft/sec\n",
+"ad= 0.0765 //lbm/ft^3\n",
+"p= 5 //per cent\n",
+"c= 0.055\n",
+"n= 1.75\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"CL= 2*w*(1-(p/100))*g/(ad*v^2*A)\n",
+"D= w*(1-(p/100))*c*n/CL\n",
+"//RESULTS\n",
+"printf ('lift coeefieicnt = %.2f',CL)\n",
+"printf ('\n Drage force = %.f',D)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"bi= 70 //degrees\n",
+"i= 8 //degrees\n",
+"bo= 130 //degrees\n",
+"s= 5 //degrees\n",
+"vi= 1200 //ft/sec\n",
+"g= 32.2 //ft/sec^2\n",
+"a= 0.48\n",
+"s1= 1.4 //in\n",
+"b= 5 //in\n",
+"Cx= 0.06\n",
+"//CALCULATIONS\n",
+"O= bo-s-bi+i\n",
+"Vo= vi*sind(bi-i)/sind(bo-s)\n",
+"Fy= -a*vi*sind(bi-i)*(s1/12)*(b/12)*(Vo*cosd(bo-s)-vi*cosd(bi-i))/g\n",
+"dp= a*(Vo^2*(1+Cx)-vi^2)/(2*g)\n",
+"//RESULTS\n",
+"printf ('Fluid deflection angle = %.f degrees',O)\n",
+"printf ('\n Vo = %.f ft/sec',Vo)\n",
+"printf ('\n Force on each blade = %.f lbf',Fy)\n",
+"printf ('\n Pressure difference = %.f lbf/ft^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 10.6: ex_6.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"ari= 62 //degrees\n",
+"aro= 125 //degrees\n",
+"vri= 1200 //ft/sec\n",
+"vro= 1294 //ft/sec\n",
+"vrr= 550 //ft/sec\n",
+"//CALCULATIONS\n",
+"v1= vri*sind(ari)\n",
+"v2= vrr+vri*cosd(ari)\n",
+"vi= sqrt(v1^2+v2^2)\n",
+"ai= atand(v1/v2)\n",
+"vo= vro*sind(aro)\n",
+"vo1= vro*cosd(aro)+vrr\n",
+"vo2= sqrt(vo^2+vo1^2)\n",
+"ao= atand(vo/vo1)+180\n",
+"//RESULTS\n",
+"printf ('absolute velocity = %.f ft/sec',vi)\n",
+"printf ('\n direction = %.1f degrees',ai)\n",
+"printf ('\n absolute velocity = %.f ft/sec',vo2)\n",
+"printf ('\n direction = %.1f degrees',ao)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/11-elementary_analysis.ipynb b/Basic_Fluid_Mechanics_by_Peerless/11-elementary_analysis.ipynb
new file mode 100644
index 0000000..79ca59c
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/11-elementary_analysis.ipynb
@@ -0,0 +1,338 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 11: elementary analysis"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"rt= 1.3 //ft\n",
+"rr= 0.6 //ft\n",
+"Q= 75 //ft^3\n",
+"rm= 0.95\n",
+"w1= 40 //rev/sec\n",
+"bim= 153 //degrees\n",
+"bom= 147 //degrees\n",
+"w= 62.4 //lb/ft^3\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"A= %pi*(rt^2-rr^2)\n",
+"Va= Q/A\n",
+"Vbm= rm*w1\n",
+"a= acotd(-Vbm/Va)\n",
+"im= a-bim\n",
+"vwm= Vbm+Va*cotd(bom)\n",
+"dvwm= rm*vwm\n",
+"C= w*Q*dvwm/g\n",
+"Cw= C*w1\n",
+"dp= Cw/Q\n",
+"//RESULTS\n",
+"printf ('Incidence = %.1f degrees',im) \n",
+"printf ('\n Oulet velocity = %.2f ft/sec',vwm)\n",
+"printf ('\n Change of whirl at the mean radius = %.2f ft^2/sec',dvwm)\n",
+"\n",
+"printf ('\n Torque = %.f lbf/ft',C)\n",
+"printf ('\n Rate of working = %.f ft lbf/sec',Cw)\n",
+"printf ('\n Workdone by the rotor = %.f lbf/ft^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"vbm= 38 //ft/sec\n",
+"va= 17.94 //ft/sec\n",
+"a= 147.5 //degrees\n",
+"vwm= 10.37 //ft/sec\n",
+"C= 1430 //lbf/ft\n",
+"P= 763 //lbf/ft^2\n",
+"//CALCULATIONS\n",
+"vwm1= vbm+va*cotd(a)\n",
+"p= (vwm-vwm1)/vwm\n",
+"C1= C*(1-p)\n",
+"P1= P*(1-p)\n",
+"//RESULTS\n",
+"printf ('Oulet Velocity = %.2f ft/sec',vwm1) \n",
+"printf ('\n Torque = %.f lbf/ft',C1)\n",
+"printf ('\n Workdone by the rotor = %.f lbf/ft^2',P1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"a= 154 //degrees\n",
+"vbm= 38 //ft/sec\n",
+"bom= 147 //degrees\n",
+"vwm= -7.78 //ft/sec\n",
+"w= 62.4 //lbf/ft^3\n",
+"g= 32.2 //ft/sec^2\n",
+"vb= 38 //ft/sec\n",
+"A= 4.18 //ft^2\n",
+"e= 0.95\n",
+"//CALCULATIONS\n",
+"vat= (vwm-vb)*tand(bom)\n",
+"Q= vat*A\n",
+"a1= acotd(-vbm/vat)\n",
+"imt= a1-a\n",
+"C= w*Q*vwm*e/g\n",
+"//RESULTS\n",
+"printf ('Flow rate = %.1f ft^3',Q)\n",
+"printf ('\n Incidence angle= %.f degrees',imt)\n",
+"printf ('\n Torque= %.f lbf ft',C)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"rt= 0.5 //ft\n",
+"rr= 0.16 //ft\n",
+"dv1= 88.3 //ft/sec\n",
+"b= 150 //degrees\n",
+"r= [0.16 0.3 0.5]\n",
+"vw= [2.5 5 7.5]\n",
+"vb= [46.6 88.3 132.5]\n",
+"vrb= [44.16 88.3 132.5]\n",
+"v1= [-1.154 -0.385]\n",
+"//CALCULATIONS\n",
+"A= %pi*(rt^2-rr^2)\n",
+"Va= -dv1*tand(b)\n",
+"Q= Va*A\n",
+"a= atand(v1)+180\n",
+"i= b-a\n",
+"//RESULTS\n",
+"printf ('Velocity = %.2f ft/sec',Va)\n",
+"printf ('\n Flow rate = %.1f ft^3',Q)\n",
+"disp(v1)\n",
+"disp(a)\n",
+"disp(i)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"rt= 0.5 //ft\n",
+"rr= 0.16 //ft\n",
+"dv1= 88.3 //ft/sec\n",
+"b= 150 //degrees\n",
+"a= 5 //degrees\n",
+"v1= [-0.933 -0.311]\n",
+"i= [1.0 5.0 6.7]\n",
+"//CALCULATIONS\n",
+"b1= b+a\n",
+"A= %pi*(rt^2-rr^2)\n",
+"Va= -dv1*tand(b1)\n",
+"Q= Va*A\n",
+"a1= atand(v1)+180\n",
+"\n",
+"//RESULTS\n",
+"printf ('Velocity = %.2f ft/sec',Va)\n",
+"printf ('\n Flow rate = %.1f ft^3/sec',Q)\n",
+"\n",
+"disp(v1)\n",
+"disp(a1)\n",
+"disp(i)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.6: ex_6.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1 //in\n",
+"b= 0.75 //in\n",
+"w= 180 //rev/sec\n",
+"B= 120 //degrees\n",
+"Bo= 150 //degrees\n",
+"ro= 3 //ft\n",
+"bo= 0.5 //ft\n",
+"Vbo= 180 //ft/sec\n",
+"w1= 62.4 //lbf/ft^3\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"Q= -2*%pi*(r/12)^2*(b/12)*w*tand(B)\n",
+"Vfo= Q/(2*%pi*(ro/12)*(bo/12))\n",
+"Vwo= Vbo*(ro/12)+Vfo*cotd(Bo)\n",
+"C= w1*Q*Vwo*(ro/12)/g\n",
+"dp= w1*Vwo*w*(ro/12)/g\n",
+"ari= atand(-Q*0.8/(2*%pi*(r/12)^2*(b/12)*w))+180\n",
+"i1= ari-B\n",
+"//RESULTS\n",
+"printf ('Flow rate = %.2f ft^3/sec',Q)\n",
+"printf ('\n radial velocity= %.2f ft/sec',Vfo)\n",
+"printf ('\n outlet whirl velocity= %.2f ft/sec',Vwo)\n",
+"printf ('\n Torque= %.2f lbf ft',C)\n",
+"printf ('\n Stagnant pressure = %.f lbf/ft^2',dp)\n",
+"printf ('\n Incidence angle = %.1f degrees',i1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 11.7: ex_7.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.4\n",
+"Mai= 0.5 //ft/sec\n",
+"T= 582 //R\n",
+"psi= 3040 //lbf/in^2\n",
+"R= 53.3 //ft lbf/lbm\n",
+"g= 32.2 //ft/sec^2\n",
+"Vwi= 300 //ft/sec\n",
+"m= 35 //lb/sec\n",
+"rm= 0.7 //ft\n",
+"rp= 4.25\n",
+"w= 1200 //rev/sec\n",
+"cp= 0.24\n",
+"J= 778 //lb\n",
+"//CALCULATIONS\n",
+"tr= 1+0.5*(r-1)*Mai^2\n",
+"Ti= T/tr\n",
+"pr= tr^(r/(r-1))\n",
+"pi= psi/pr\n",
+"ai= pi/(R*Ti)\n",
+"Vi= Mai*(r*R*g*Ti)^0.5\n",
+"Vai= sqrt(Vi^2-Vwi^2)\n",
+"h= m/(2*%pi*ai*rm*Vai)\n",
+"pr1= rp^(1/12)\n",
+"Vwo= Vwi+(pr1^((r-1)/r)-1)*(cp*J*g*T/(rm*w))\n",
+"BO= acotd((Vwo-w*rm)/Vai)\n",
+"//RESULTS\n",
+"printf ('Absolute air velocity = %.f ft/sec',Vi)\n",
+"printf ('\n air velocity = %.f ft/sec',Vai)\n",
+"printf ('\n Blade height = %.3f ft',h)\n",
+"printf ('\n velocity = %.f ft/sec',Vwo)\n",
+"printf ('\n outlet balde angle = %.1f degrees',BO)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/12-TURBOMACHINES.ipynb b/Basic_Fluid_Mechanics_by_Peerless/12-TURBOMACHINES.ipynb
new file mode 100644
index 0000000..76886a1
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/12-TURBOMACHINES.ipynb
@@ -0,0 +1,222 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 12: TURBOMACHINES"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"d= 0.0764 //lbm/ft^3\n",
+"u= 3.74*10^-7 //lbf sec/ft^2\n",
+"D= 15 //in\n",
+"g= 32.2 //ft/sec^2\n",
+"p= 14.7 //lb/in^2\n",
+"r1= [0.02 0.04 0.06 0.08 0.1 1.15]\n",
+"r2= [0.0338 0.0267 0.0199 0.0159 0.0132 0.0100]\n",
+"r3= [0.46 0.92 1.38 1.84 2.3 2.64]\n",
+"r4= [2.97 2.35 1.75 1.4 1.16 0.88]\n",
+"r5= [0.0206 0.0163 0.0121 0.0097 0.0081 0.0061]\n",
+"//CALCULATIONS\n",
+"re= (d/u)*(p*100*2*%pi/60)*(D/12)^2/g\n",
+"//RESULTS\n",
+"printf ('Reynolds Number = %.2e ',re) \n",
+"xtitle('','m lbm/sec', 'dPs lbf/ft^2') \n",
+"\n",
+"disp(r1)\n",
+"disp(r2)\n",
+"disp(r3)\n",
+"disp(r4)\n",
+"disp(r5)\n",
+"plot(r3,r5)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"psif= 10.2 //lbf/in^2\n",
+"usit= 3.8*10^-7 //lbf sec/ft^2\n",
+"usif= 3.52*10^-7 //lbf sec/ft^2\n",
+"Tsit= 530 //R\n",
+"Tsif= 480 //R\n",
+"wf= 15000 //rev/min\n",
+"//CALCULATIONS\n",
+"Psit= psif*usit*sqrt(Tsit/Tsif)/usif\n",
+"wt= wf*sqrt(Tsit/Tsif)\n",
+"//RESULTS\n",
+"printf ('Pressure in the test cell = %.1f lbf/in^2',Psit) \n",
+"printf ('\n Compressor speed = %.f rev.min',wt) "
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 62.3 //lbf/ft^3\n",
+"d= 0.375 //in\n",
+"ro= 0.75 //ft\n",
+"l= 1.25 //ft\n",
+"b= 120 //degrees\n",
+"do= 0.25 //in\n",
+"p= 750 //lbf/in^2\n",
+"g= 32.1 //ft/sec^2\n",
+"f= 0.03\n",
+"f1= 0.9\n",
+"f2= 0.3\n",
+"w1= 60 //rad/sec\n",
+"//CALCULATIONS\n",
+"Q= sqrt(((p/w)+((60*ro)^2/(2*g))+do)*%pi^2*g*(d/12)^4/((d/do)^4-1+(l*f/(d/12))+f1+f2))*0.353\n",
+"Vwo= w1*ro+(4*Q/(%pi*(do/12)^2))*cosd(b)\n",
+"C= w*Q*Vwo*ro/g\n",
+"//RESULTS\n",
+"printf ('Flow Rate = %.4f ft^3/sec',Q) \n",
+"printf ('\n Vwo = %.2f ft/sec',Vwo) \n",
+"printf ('\n Driving Torque = %.3f lbf ft',C) "
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"W= 38 //rev/sec\n",
+"w= 62.4 //lbf/ft^3\n",
+"m= 2000 //lbm/sec\n",
+"g= 32.2 //ft/sec^2\n",
+"ps= 5000 //lbf/ft^2\n",
+"S3= 4.6\n",
+"e= 0.91\n",
+"//CALCULATIONS\n",
+"S1= W*(w*m^2/(g*ps)^3)^0.25\n",
+"D= S3*(m^2/(w*g*ps))^0.25\n",
+"//RESULTS\n",
+"printf ('S1 = %.3f',S1) \n",
+"printf ('\n Diameter = %.2f ft',D) \n",
+"printf ('\n efficiency = %.2f ',e)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 12.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"d= 6 //in\n",
+"f= 0.25\n",
+"l= 1200 //ft\n",
+"p= 55 //lbm/ft^3\n",
+"w= 740 //rev/min\n",
+"g= 32.2 //ft/sec^2\n",
+"n= 0.87\n",
+"d1= 1.78 //ft\n",
+"//CALCULATIONS\n",
+"D= (0.13*%pi^2*(d/12)^5/(8*f*l*0.012^2))^0.25*d1\n",
+"m= 0.012*p*(w*2*%pi/60)*D^3\n",
+"dps= 0.13*p*(w*2*%pi*D/60)^2/g\n",
+"P= m*10*dps/(p*n)\n",
+"//RESULTS\n",
+"printf ('Diameter = %.2f ft',D) \n",
+"printf ('\n Mass flow rate = %.1f lbm/sec',m) \n",
+"printf ('\n pressure rise = %.1f lbf/ft^2',dps)\n",
+"printf ('\n shaft power = %.2e ft lbf/sec',P)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/13-Hydraulic_power_transmission.ipynb b/Basic_Fluid_Mechanics_by_Peerless/13-Hydraulic_power_transmission.ipynb
new file mode 100644
index 0000000..8d88e4f
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/13-Hydraulic_power_transmission.ipynb
@@ -0,0 +1,152 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 13: Hydraulic power transmission"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 13.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"nop= 0.88\n",
+"nom= 0.88\n",
+"Pm= 75 //hp\n",
+"p= 3000 //lb/in^2\n",
+"d= 54.5 //lbm/ft^3\n",
+"u= 1.05*10^-4\n",
+"d1= 0.5 //in\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"nt= (7/11)*nop*nom\n",
+"pp= Pm/nt\n",
+"Q= nop*pp*550/(p*144)\n",
+"Re= 4*d*Q/(%pi*u*(d1/12)*g)\n",
+"//RESULTS\n",
+"printf ('n trans = %.3f ',nt)\n",
+"printf ('\n Input power = %.f hp',pp)\n",
+"printf ('\n Flow rate = %.3f ft^3/sec',Q)\n",
+"printf ('\n Reynolds Number = %.1e ',Re)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 13.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"lc= 0.25\n",
+"a= 90 //degrees\n",
+"p= 3000 //lb/in^2\n",
+"g= 32.2 //ft/sec^2\n",
+"d1= 0.5 //in\n",
+"Q= 0.171 //ft^3/sec\n",
+"d= 54.5 //lbm/ft^3\n",
+"n1= 2 \n",
+"n2= 6\n",
+"lc1= 0.9\n",
+"nop= 0.88\n",
+"nom= 0.88\n",
+"//CALCULATIONS\n",
+"P1= 4*p*144/11\n",
+"P2= 8*d*Q^2*(n1*lc+n2*lc1)/(%pi^2*(d1/12)^4*g)\n",
+"pt= P1+P2\n",
+"dpm= (p*144-pt)\n",
+"ntrans= nop*nom*dpm/(p*144)\n",
+"//RESULTS\n",
+"printf ('Frictional pressure drop = %.2e lbf/ft^2',P1) \n",
+"printf ('\n Extra Frictional pressure drop = %.2e lbf/ft^2',P2) \n",
+"printf ('\n Total pressure drop = %.2e lbf/ft^2',pt)\n",
+"printf ('\n Motor pressure drop = %.2e lbf/ft^2',dpm)\n",
+"printf ('\n Overall transmission coefficiency = %.3f',ntrans)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 13.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"bip= 135 //degrees\n",
+"bop= 150 //degrees\n",
+"bot= 140 //degrees\n",
+"bos= 137 //degrees\n",
+"r= 1.8\n",
+"r1= 1.8\n",
+"r2= 0.7\n",
+"r3= 0.95\n",
+"//CALCULATIONS\n",
+"R= (1+(cotd(bip)/cotd(bos)))*r^2-r1*(cotd(bop)/cotd(bos))\n",
+"R1= r2*r3^2*(1+(cotd(bip)/cotd(bos)))-(cotd(bot)/cotd(bos))\n",
+"R2= (R1-R)/(R-1)\n",
+"//RESULTS\n",
+"printf ('R1 = %.2f',R) \n",
+"printf ('\n R2 = %.2f',R1) \n",
+"printf ('\n Torque ratio = %.2f',R2) "
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/14-Further_Devolopments.ipynb b/Basic_Fluid_Mechanics_by_Peerless/14-Further_Devolopments.ipynb
new file mode 100644
index 0000000..e7d39ba
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/14-Further_Devolopments.ipynb
@@ -0,0 +1,111 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 14: Further Devolopments"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 14.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"a= 60.5\n",
+"Q= 0.2 //ft^3/sec\n",
+"d= 3 //in\n",
+"u= 0.0325\n",
+"g= 32.2 //ft/sec^2\n",
+"T= [50.0 60.0 70.0 80.0 90.0 100.0]\n",
+"Ep= [294.5 188.6 113.2 60.4 37.7 24.5]\n",
+"Eh= [0 69.9 139.8 209.7 279.5 349.4]\n",
+"Et= [295 258 253 270 317 374]\n",
+"//CALCULATIONS\n",
+"re= a*4*Q/(%pi*(d/12)*u*g)\n",
+"//RESULTS\n",
+"printf ('Reynolds Number = %.1f ',re)\n",
+"disp(T)\n",
+"disp(Ep)\n",
+"disp(Eh)\n",
+"disp(Et)\n",
+"plot(T,Ep)\n",
+"plot(T,Eh)\n",
+"plot(T,Et)\n",
+"\n",
+"xtitle('','T (F)', 'Eh,Ep,Eh&Ep (kW)')\n",
+"\n",
+""
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 14.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"wcb= 2 //ton\n",
+"wc= 100 //ton\n",
+"wa= 6.5 //ton\n",
+"wca= 20 \n",
+"r= 0.8\n",
+"r1= 1.2\n",
+"//CALCULATIONS\n",
+"wca1= wc/wa\n",
+"wca2= wcb*(wca1/wca)^1.5\n",
+"Wca= wcb*r^(9/4)*(1/r1)^(9/4)*(wca1/wca)^1.5\n",
+"//RESULTS\n",
+"printf ('(Wc/W)a = %.2f ',wca1)\n",
+"printf ('\n Wc,a = %.2f ton',wca2)\n",
+"printf ('\n Wc,a = %.2f ton',Wca)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/2-Similarity.ipynb b/Basic_Fluid_Mechanics_by_Peerless/2-Similarity.ipynb
new file mode 100644
index 0000000..1e153fd
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/2-Similarity.ipynb
@@ -0,0 +1,244 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 2: Similarity"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.1: Force_applied.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 4\n",
+"l1= 4 //units\n",
+"l2= 10 //units\n",
+"//CALCULATIONS\n",
+"sxy= (4/r)\n",
+"sxy1= l1^2\n",
+"sxy2= l2^2\n",
+"//RESULTS\n",
+"printf ('x^2+4*y^2 = %.f ',sxy)\n",
+"printf ('\n x^2+4*y^2 = %.f ',sxy1)\n",
+"printf ('\n x^2+4*y^2 = %.f ',sxy2)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.3: force_required.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"vo= 10 //ft/sec\n",
+"a= 0.5 //ft^-1\n",
+"b= 1 //ft\n",
+"x= -2 //ft\n",
+"y= 2 //ft\n",
+"b1= 2\n",
+"a1= 3/5 //ft\n",
+"//CALCULATIONS\n",
+"Vx= vo/(a*x^2+b)\n",
+"Vy= -2*a*b*vo*x*y/(a*x^2+b)^2\n",
+"V= sqrt(Vx^2+Vy^2)\n",
+"fx= -2*a*b^2*vo^2*x/(a*x^2+b)^3\n",
+"fy= 2*a*b^2*vo^2*y*(b-a*x^2)/(a*x^2+b)^4\n",
+"f= sqrt(fx^2+fy^2)\n",
+"r= b1^2/a1\n",
+"f1= f*r\n",
+"//RESULTS\n",
+"printf ('Vx = %.2f ft/sec',Vx)\n",
+"printf ('\n Vx = %.2f ft/sec',Vy)\n",
+"printf ('\n V = %.2f ft/sec',V)\n",
+"printf ('\n fx = %.2f ft/sec^2',fx)\n",
+"printf ('\n fy = %.2f ft/sec^2',fy)\n",
+"printf ('\n f = %.2f ft/sec^2',f)\n",
+"printf ('\n r = %.2f in the present case',r)\n",
+"printf ('\n f1 = %.2f ft/sec^2',f1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.4: force_at_the_end.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1/5\n",
+"b1= 2 //ft\n",
+"a1= 3/5 //ft\n",
+"//CALCULATIONS\n",
+"r= (a1*b1)^2*r\n",
+"//RESULTS\n",
+"printf ('ratio of resultant forces acting on coorresponding fluid elements = %.2f ',r)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.5: air_speed.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"vos= 70 //ft/sec\n",
+"as= 78 //ft\n",
+"am= 72 //ft\n",
+"ls1= 6 //ft\n",
+"lm= 2 //ft\n",
+"um= 386 //ft/sec\n",
+"us= 372 //ft/sec\n",
+"dm= 0.4\n",
+"//CALCULATIONS\n",
+"vom= vos*as*ls1*um/(am*lm*us)\n",
+"Ds= dm*(am/as)*(us/um)^2\n",
+"//RESULTS\n",
+"printf ('Air speed = %.f ft/sec',vom)\n",
+"printf ('\n Ds = %.3f lbf',Ds)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.6: ratio_of_resultant_forces.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"vom= 236 //ft/sec\n",
+"as= 0.072 //ft\n",
+"am= 62.4 //ft\n",
+"ls1= 2 //ft\n",
+"lm= 8 //ft\n",
+"um= 248 //ft/sec\n",
+"us= 3.86 //ft/sec\n",
+"r= 0.4/3.3\n",
+"//CALCULATIONS\n",
+"voh= vom*as*ls1*um/(am*lm*us)\n",
+"Ds= r*(as/am)*(um/us)^2*(ls1/lm)*(lm-ls1)\n",
+"//RESULTS\n",
+"printf ('Air speed = %.2f ft/sec',voh)\n",
+"printf ('\n Drag force = %.3f lbf',Ds)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 2.7: Temperature_a_exit.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"To1= 540 //R\n",
+"po3= 12.6 //lbf/in^2\n",
+"l3= 3 //ft\n",
+"po1= 14.7 //lbf/in^2\n",
+"l1= 1 //ft\n",
+"vo1= 500 //ft/sec\n",
+"r= 0.83\n",
+"P1= 1 //lbf/in^2\n",
+"//CALCULATIONS\n",
+"To3= To1*(po3*l3/(po1*l1))^r\n",
+"Vo3= vo1*sqrt(To3/To1)\n",
+"P3= P1*po3*l3/(po1*l1)\n",
+"//RESULTS\n",
+"printf ('To3 = %.f R',To3)\n",
+"printf ('\n Vo3 = %.f ft/sec',Vo3)\n",
+"printf ('\n P3 = %.2f lbf/ft',P3)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/3-Dimensional_Analaysis.ipynb b/Basic_Fluid_Mechanics_by_Peerless/3-Dimensional_Analaysis.ipynb
new file mode 100644
index 0000000..ce10c7c
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/3-Dimensional_Analaysis.ipynb
@@ -0,0 +1,132 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 3: Dimensional Analaysis"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 3.1: dimensions_of_g.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"g= 32.2 //ft/sec^2\n",
+"t= 1 //hr\n",
+"g1= 32.2 //ft/sec^2\n",
+"g2= 32.2 //lbm ft/lbf\n",
+"u= 2.4*10^-5 //lbf sec/ft^2\n",
+"//CALCULATIONS\n",
+"q2= g*(t*60*60)^2\n",
+"go= g*(t*60*60)^2\n",
+"q3= g/g2\n",
+"u1= u/(t*60*60)\n",
+"//RESULTS\n",
+"printf ('q2 = %.2e lbm ft/lbf hr^2',q2)\n",
+"printf ('\n go = %.2e lbm ft/lbf hr^2',go)\n",
+"printf ('\n go = %.f slug ft/lbf sec^2',q3)\n",
+"printf ('\n viscosity = %.2e lbf hr/ft^2',u1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 3.2: dimensions_of_vectors.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"g= 32.2 //ft/sec^2\n",
+"m= 1 //lb\n",
+"//CALCULATIONS\n",
+"m1= g/m\n",
+"//RESULTS\n",
+"printf ('1 lbf/sec ft^2 = %.1f lbm/ft sec',m1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 3.5: dimensions_of_velocity.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"n1=1\n",
+"n2= 3\n",
+"n3=2\n",
+"//CALCULATIONS\n",
+"a1= -n1\n",
+"a2= -n3\n",
+"a3= -n1-a2+3*a1\n",
+"b1= -n1\n",
+"b2= -n1\n",
+"b3= n1+3*b1-b2\n",
+"//RESULTS\n",
+"printf ('a1 = %.f ',a1)\n",
+"printf ('\n a2 = %.f ',a2)\n",
+"printf ('\n a3 = %.f ',a3)\n",
+"printf ('\n b1 = %.f ',b1)\n",
+"printf ('\n b2 = %.f ',b2)\n",
+"printf ('\n b3 = %.f ',b3)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/5-Control_Volume_Analysis.ipynb b/Basic_Fluid_Mechanics_by_Peerless/5-Control_Volume_Analysis.ipynb
new file mode 100644
index 0000000..b82ea2e
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/5-Control_Volume_Analysis.ipynb
@@ -0,0 +1,170 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 5: Control Volume Analysis"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 4.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 20 //lbm/sec\n",
+"sh= 0.004\n",
+"m1= 0.12 //lbm/sec\n",
+"m2= 12.2 //lbm/sec\n",
+"m3= 0.130 //lbm/sec\n",
+"//CALCULATIONS\n",
+"mw1= w/((1/sh)+1)\n",
+"ma1= w-mw1\n",
+"ma4= ma1-m2\n",
+"mw4= mw1+m1-m3\n",
+"mr= ma4+mw4\n",
+"sh1= mw4/ma4\n",
+"//RESULTS\n",
+"printf ('mw1 = %.4f lbm/sec',mw1)\n",
+"printf ('\n ma1 = %.2f lbm/sec',ma1)\n",
+"printf ('\n ma4 = %.2f lbm/sec',ma4)\n",
+"printf ('\n mw4 = %.2f lbm/sec',mw4)\n",
+"printf ('\n mr = %.2f lbm/sec',mr)\n",
+"printf ('\n specific humidity = %.5f lbm/sec',sh1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 5.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 62.4 //lbf/ft^3\n",
+"g= 32.2 //ft/sec^2\n",
+"v= 86.5 //ft/sec\n",
+"d2= 3 //in\n",
+"d1= 6 //in\n",
+"dp= 50 //lbf/in^2\n",
+"//CALCULATIONS\n",
+"Fb= -((%pi*(w/g)*v^2*(1/d1)^2*(1-(d2/d1)^2)*0.25)-dp*144*(%pi/4)*(1/d2)^2)\n",
+"//RESULTS\n",
+"printf ('Load on the bolts = %.f lbf',Fb)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 5.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"F1= 237 //lb\n",
+"dp= 50 //lbf/in^2\n",
+"D= 6 //in\n",
+"//CALCULATIONS\n",
+"F2= dp*144*(%pi/4)*(D/12)^2\n",
+"Fb= F1-F2\n",
+"//RESULTS\n",
+"printf ('Load on the bolts = %.f lbf',Fb)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 5.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w1= 0.0286 //lbm/ft^3\n",
+"v= 2500 //ft/sec\n",
+"A= 2.5 //ft^3\n",
+"k= 0.015\n",
+"p2= 700 //lbf/ft^2\n",
+"p1= 628 //lbf/ft^2\n",
+"v2= 3500 //ft/sec\n",
+"g= 32.17 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"ma= w1*v*A\n",
+"mf= k*ma\n",
+"mt= ma+mf\n",
+"F= (p2-p1)*A+(mt*v2/g)-(ma*v/g)\n",
+"//RESULTS\n",
+"printf ('air mass flow rate = %.2f lbm/sec',ma)\n",
+"printf ('\n Fuel flow rate = %.2f lbm/sec',mf)\n",
+"printf ('\n Fuel flow rate at station 2 = %.2f lbm/sec',mt)\n",
+"printf ('\n Thrust force = %.f lbf',F)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/6-Steady_one_dimensional_reversible_flow.ipynb b/Basic_Fluid_Mechanics_by_Peerless/6-Steady_one_dimensional_reversible_flow.ipynb
new file mode 100644
index 0000000..38d269e
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/6-Steady_one_dimensional_reversible_flow.ipynb
@@ -0,0 +1,311 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 6: Steady one dimensional reversible flow"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.10: ex_10.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"Q= 400 //ft^3/sec\n",
+"b1= 25 //ft\n",
+"b2= 20 //ft\n",
+"h1= 6 //ft\n",
+"z1= 2.5 //ft\n",
+"z2= 3.3 //ft\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"hc1= (Q^2/(g*b1^2))^(1/3)\n",
+"hc2= (Q^2/(g*b2^2))^(1/3)\n",
+"r= (hc1/hc2)*((h1/hc1)+0.5*(hc1/h1)^2)+((z1-z2)/hc2)\n",
+"//RESULTS\n",
+"printf ('hc1 = %.3f ft',hc1)\n",
+"printf ('\n hc2 = %.3f ft',hc2)\n",
+"printf ('\n Ratio = %.3f ',r)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"g= 32.2 //ft/sec^2\n",
+"h= 4 //ft\n",
+"d2= 0.16 //ft\n",
+"d1= 0.3 //ft\n",
+"dp= 12.6 //lbf/in^2\n",
+"//CALCULATIONS\n",
+"Q= (%pi/4)*sqrt(2*g*dp*h/((1/d2^4)-(1/d1^4)))\n",
+"//RESULTS\n",
+"printf ('Volumetric flow rate = %.2f ft^3/sec',Q)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 0.0765 //lbm/ft^3\n",
+"v1= 120 //ft/sec\n",
+"g= 32.2 //ft/sec^2\n",
+"//CALCULATIONS\n",
+"dp= w*v1^2/(2*2*g)\n",
+"//RESULTS\n",
+"printf ('Difference in pressure= %.2f lbf/ft^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r=1.4\n",
+"g= 32.2 //ft/sec^2\n",
+"R= 53.3 //lbf ft/lbm\n",
+"T1= 760 //R\n",
+"p2= 2 //lbf/in^2\n",
+"p1= 3 //lbf/in^2\n",
+"//CALCULATIONS\n",
+"V2= sqrt(2*r*R*g*T1*(1-(p2/p1)^((r-1)/r))/(r-1))\n",
+"//RESULTS\n",
+"printf ('Velocity in working section = %.f ft/sec',V2)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.4\n",
+"g= 32.2 //ft/sec^2\n",
+"R= 53.3 //lbf ft/lbm\n",
+"T= 32 //C\n",
+"T1= 2000 //R\n",
+"r1= 1.32\n",
+"p= 1440 //lbf/in^2\n",
+"v1= 1.2306 //ft^3/lbm\n",
+"v2= 1.2546 //ft^3/lbm\n",
+"bm= 3.13*10^5 //lbf/in^2\n",
+"w= 62.4 //lbf/ft^3\n",
+"//CALCULATIONS\n",
+"a1= sqrt(r*R*(460+T)*g)\n",
+"a2= sqrt(r1*R*T1*g)\n",
+"r2= p/(v1-v2)\n",
+"a3= sqrt(-g*(v1+v2)^2*0.5^2*r2)\n",
+"a4= sqrt(bm*144*g/w)\n",
+"//RESULTS\n",
+"printf ('Acoustic veloctiy = %.f ft/sec',a1)\n",
+"printf ('\n Acoustic veloctiy = %.f ft/sec',a2)\n",
+"printf ('\n Acoustic veloctiy = %.f ft/sec',a3)\n",
+"printf ('\n Acoustic veloctiy = %.f ft/sec',a4)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.4\n",
+"ma2= 2.5 //ft/sec\n",
+"g= 32.17 //ft/sec^2\n",
+"p2= 1 //lbf/in^2\n",
+"ps= 17.08 //lbf/in^2\n",
+"ps2= 75 //lbf/in^2\n",
+"Ts= 720 //R\n",
+"R= 53.3 //lbf ft/lbm\n",
+"A= 4 //ft^2\n",
+"ps3= 0.4 //lbf/in^2\n",
+"A2= 0.685 //ft^2\n",
+"P= 5 //per cent\n",
+"//CALCULATIONS\n",
+"R1= (1+0.5*(r-1)*ma2^2)^(r/(r-1))\n",
+"R2= (2*(r/(r-1))*(p2/ps)^(2/(r))*(1-(p2/ps)^((r-1)/r)))^0.5\n",
+"m2= R2*ps2*144*(g/(R*Ts))^0.5*0.1\n",
+"m= m2*A\n",
+"At= A*R2/A2\n",
+"m1= m*(1-(P/100))\n",
+"mrp= (1-(P/100))*R2\n",
+"//RESULTS\n",
+"printf ('Mass flow rate= %.1f lbm/sec',m)\n",
+"printf ('\n Area of throat= %.3f ft^2',At)\n",
+"printf ('\n Mass flow rate= %.1f lbm/sec',m1)\n",
+"printf ('\n Mass flow rate parameter = %.4f',mrp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.7: ex_7.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r1= 10 //ft\n",
+"r2= 0.2 //miles\n",
+"w= 0.0765 //lbm/ft^2\n",
+"g= 32.2 //ft/sec^2\n",
+"V1= 1 //ft/sec\n",
+"//CALCULATIONS\n",
+"k= r2*5280*V1 \n",
+"dp= w*k^2*10*((1/r1)^2-(1/(5280*r2))^2)/(2*g)\n",
+"//RESULTS\n",
+"printf ('k = %.f ft^2/sec',k)\n",
+"printf ('\n pressure difference = %.1f lbf/ft^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 6.9: ex_9.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 12 //ft\n",
+"q= 300 //ft^3/sec\n",
+"h= 10 //ft\n",
+"g= 32.2 //ft/sec^2\n",
+"R= 2.6\n",
+"//CALCULATIONS\n",
+"hc= ((q/12)^2/g)^(1/3)\n",
+"r= h/hc\n",
+"h1= hc*(((h/hc)+0.5*(hc/h)^2)-0.5*R^2)\n",
+"//RESULTS\n",
+"printf ('hc = %.2f ft',hc)\n",
+"printf ('\n stream depth = %.2f ft',h1)\n",
+""
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/7-Steady_one_dimensional_Irreversible_flow.ipynb b/Basic_Fluid_Mechanics_by_Peerless/7-Steady_one_dimensional_Irreversible_flow.ipynb
new file mode 100644
index 0000000..49cabdd
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/7-Steady_one_dimensional_Irreversible_flow.ipynb
@@ -0,0 +1,402 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 7: Steady one dimensional Irreversible flow"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.10: ex_10.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"p2= 67.2 //lbf/in^2\n",
+"p1= 63 //lbf/in62\n",
+"r= 1.4\n",
+"n= 0.6\n",
+"T1= 870 //R\n",
+"ma1= 0.8 //ft/sec\n",
+"//CALCULATIONS\n",
+"dt= (p2/p1)^((r-1)/r)-1\n",
+"dt1= dt/n\n",
+"T2= T1*(1+dt1)\n",
+"Ts1= T1*(1+0.5*(r-1)*ma1^2)\n",
+"ps1= p1*(1+0.5*(r-1)*ma1^2)^(r/(r-1))\n",
+"ps2= p2*(Ts1/T2)^(r/(r-1))\n",
+"dp= ps1-ps2\n",
+"//RESULTS\n",
+" printf ('dT = %.5f ',dt)\n",
+" printf ('\n dT1 = %.5f ',dt1)\n",
+" printf ('\n Temperature = %.f R',T2)\n",
+" printf ('\n Temperature = %.1f R',Ts1)\n",
+" printf ('\n Pressure = %.1f lbf/in^2',ps1)\n",
+" printf ('\n Pressure = %.1f lbf/in^2',ps2)\n",
+" printf ('\n pressure difference = %.1f lbf/in^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.11: ex_11.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.4\n",
+"ma3= 3 //ft/sec\n",
+"ps= 80 //lbf/ft^2\n",
+"Ts= 840 //R\n",
+"r1= 53.3 //ft lbm/ft^3\n",
+"A3= 2 //in^2\n",
+"g= 32.2 //ft/sec^2\n",
+"ma1= 1.6\n",
+"//CALCULATIONS\n",
+"R= (1+(r-1)*0.5*ma3^2)^(r/(r-1))\n",
+"p3= ps/R\n",
+"R1= 1+0.5*(r-1)*ma3^2\n",
+"T3= Ts/R1\n",
+"w3= p3*144/(r1*T3)\n",
+"V3= ma3*sqrt(r*r1*g*T3)\n",
+"m= w3*V3*A3/144\n",
+"ra= ((r+1)/(2*r*ma1^2-(r-1)))^(1/(r-1))*(0.5*(r+1)*ma1^2/(1+0.5*(r-1)*ma1^2))^(r/(r-1))\n",
+"ps2= ps*ra\n",
+"dp= ps-ps2\n",
+"//RESULTS\n",
+" printf ('outlet pressure = %.2f lbf/in^2',p3)\n",
+" printf ('\n outlet temperature = %.f R',T3)\n",
+" printf ('\n mass flow rate = %.3f lbm/sec',m)\n",
+" printf ('\n mass flow rate = %.3f lbm/sec',m)\n",
+" printf ('\n ps2 = %.1f lbf/in^2',ps2)\n",
+" printf ('\n pressure difference = %.1f lbf/in^2',dp)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.5\n",
+"f= 0.025\n",
+"//CALCULATIONS\n",
+"r1= (2/f)*(r^2-1)\n",
+"//RESULTS\n",
+"printf ('ratio L/D2 = %.f',r1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"a= 6 //degrees\n",
+"r= 1.5\n",
+"l= 100 //ft\n",
+"f= 0.025\n",
+"K= 0.15\n",
+"//CALCULATIONS\n",
+"R= r^4-1\n",
+"R1= cotd(a/2)*(1-(1/r))\n",
+"p1= f*l\n",
+"p2= 2.5*(l-p1)/l\n",
+"p3= (1-r^2)^2\n",
+"p4= K*p3\n",
+"pt= p4+p2\n",
+"//RESULTS\n",
+"printf ('lowest ratio = %.2f',R)\n",
+"printf ('\n contribtuion of friction in pipe = %.1f lbf/ft^2',p1)\n",
+"printf ('\n contribtuion of diffuser in pipe = %.3f lbf/ft^2',p2)\n",
+"printf ('\n stagnant pressure drop = %.3f lbf/ft^2',p3)\n",
+"printf ('\n contribtuion of friction in pipe after reduction = %.3f lbf/ft^2',pt)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"d= 4 //in\n",
+"q= 0.5 //ft^3/sec\n",
+"w= 62.4 //lb/ft^3\n",
+"u= 2.7*10^-5 //lbf sec/ft^2\n",
+"e= 0.0005 //ft\n",
+"g= 32.1 //ft/sec^2\n",
+"f= 0.0235\n",
+"lt= 400 //ft\n",
+"//CALCULATIONS\n",
+"V= 4*q/(%pi*(d/12)^2)\n",
+"Re= w*V*(d/12)/(u*g)\n",
+"r= e/(d/12)\n",
+"dz= (V^2/(2*g))*(1.7+f*lt/(d/12))\n",
+"//RESULTS\n",
+"printf ('mean flow velocity = %.2f ft/sec',V)\n",
+"printf ('\n Reynolds number = %.2e',Re)\n",
+"printf ('\n Relative roughness = %.4f',r)\n",
+"printf ('\n difference in the levels of water = %.1f ft',dz)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"d= 4 //in\n",
+"v= 6.64 //ft/sec\n",
+"//CALCULATIONS\n",
+"Q= %pi*0.25*(d/12)^2*v\n",
+"//RESULTS\n",
+"printf ('Flow rate= %.3f ft^3/sec',Q)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.5: ex_5.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"d= 0.366 //ft\n",
+"i= 12\n",
+"//CALCULATIONS\n",
+"pd= d*i\n",
+"//RESULTS\n",
+"printf ('Required pipe diameter = %.2f in',pd)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.6: ex_6.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"Ps1= 1050 //lbf/ft^2\n",
+"fr= 10.7\n",
+"p= 36.6 //lbf/ft^2\n",
+"p1= 195 //lbf/ft^2\n",
+"fr1= 16\n",
+"fr2= 1.8\n",
+"//CALCULATIONS\n",
+"p2= fr*p\n",
+"dp= Ps1-p2\n",
+"lc= dp/p\n",
+"sp= Ps1+p1-p*(fr1+fr2)\n",
+"lc1= sp/p\n",
+"//RESULTS\n",
+"printf ('Pressure = %.f lbf/ft^2',p1)\n",
+"printf ('\n pressure difference = %.f lbf/ft^2',dp)\n",
+"printf ('\n Loss coefficient = %.f ',lc)\n",
+"printf ('\n Loss coefficient = %.1f ',lc1)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.7: ex_7.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"p1= 50 //lbf/in^2\n",
+"R= 96.3 //ft lbf/lbm R\n",
+"T= 80 //F\n",
+"p2= 20 //lbf/in^2\n",
+"r= 1.31\n",
+"u= 2.34*10^-7 //lbf sec/ft^2\n",
+"e= 0.00005 //ft\n",
+"m= 5*10^4 //lbm/sec\n",
+"d= 1.5 //ft\n",
+"g= 32.2 //ft/sec^2\n",
+"f= 0.113\n",
+"//CALCULATIONS\n",
+"w1= p1*144/(R*(460+T))\n",
+"V1= 4*(m/3600)/(%pi*w1*d^2)\n",
+"Ma1= V1/(r*R*g*(460+T))^0.5\n",
+"Re= w1*V1*d/(u*g)\n",
+"dx= (((1/(r*Ma1^2))*10*(1-(p2/p1)^2))+log(p2/p1))*d/f\n",
+"//RESULTS\n",
+"printf ('density = %.3f lbm/ft^3',w1)\n",
+"printf ('\n mean flow velocity = %.1f ft/sec',V1)\n",
+"printf ('\n Match number = %.4f ',Ma1)\n",
+"printf ('\n Reynolds number = %.2e ',Re)\n",
+"printf ('\n Length of pipe = %.2e ft',dx)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 7.9: ex_9.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r= 1.4\n",
+"R= 53.3 //ft lbf/lbm R\n",
+"g= 32.2 //ft/sec^2\n",
+"T1= 410 //R\n",
+"v= 2500 //ft/sec\n",
+"P1= 628 //lbf/in^2\n",
+"//CALCULATIONS\n",
+"v1= sqrt(r*g*R*T1)\n",
+"Ma1= v/v1\n",
+"Ts1= T1*(1+0.5*(r-1)*Ma1^2)\n",
+"Ps1= P1*(1+0.5*(r-1)*Ma1^2)^(r/(r-1))\n",
+"Ps2= Ps1*((r+1)/(2*r*Ma1^2-r+1))^(1/(r-1))*(0.5*(r+1)*Ma1^2/(1+0.5*(r-1)*Ma1^2))^(r/(r-1))\n",
+"//RESULTS\n",
+"printf ('acoustic velocity = %.f ft/sec',v1)\n",
+"printf ('\n Match number = %.2f ',Ma1)\n",
+"printf ('\n Stagnition temperature = %.f R',Ts1)\n",
+"printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps1)\n",
+"printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps2)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/8-analysis_of_dimensional_constant_density_laminar_flow.ipynb b/Basic_Fluid_Mechanics_by_Peerless/8-analysis_of_dimensional_constant_density_laminar_flow.ipynb
new file mode 100644
index 0000000..41f3ea4
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/8-analysis_of_dimensional_constant_density_laminar_flow.ipynb
@@ -0,0 +1,125 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 8: analysis of dimensional constant density laminar flow"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 8.2: ex_2.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"w= 78.9 //lbf.ft^3\n",
+"d= 0.01 //in\n",
+"u= 8.67*10^-9 //lbf/ hr ft^2\n",
+"h= 18 //ft\n",
+"l= 10 //ft\n",
+"//CALCULATIONS\n",
+"Q= %pi*w*(d/12)^4*(h+l)/(l*128*u)\n",
+"//RESULTS\n",
+"printf ('Flow rate = %.2e ft^3/hr',Q)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 8.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"x= 0.1 //ft\n",
+"w= 62.4 //lbf/ft^3\n",
+"v1= 10 //ft/sec\n",
+"u= 2.4*10^-5 //lbf/ft\n",
+"g= 32.2 //ft/sec^2\n",
+"k= 4.91\n",
+"//CALCULATIONS\n",
+"s= k*x*(w*v1*x/(u*g))^-0.5\n",
+"Tw= 0.332*w*v1^2*(u*g/(w*x*v1))^0.5/g\n",
+"R= 0.332*6*Tw\n",
+"//RESULTS\n",
+"printf ('Thickness = %.2e*ft',s)\n",
+"printf ('\n Shear stress = %.3f lbf/ft^2',Tw)\n",
+"printf ('\n Shear stress = %.3f lbf/ft',R)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 8.4: ex_4.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"r=1\n",
+"r1=1\n",
+"//CALCULATIONS\n",
+"e1= r+r1\n",
+"e2= r-r1\n",
+"//RESULTS\n",
+"printf ('vorticity = %.f*k',e1)\n",
+"printf ('\n vorticity = %.f',e2)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Basic_Fluid_Mechanics_by_Peerless/9-analysis_of_dimensional_constant_density_turbulent_flow.ipynb b/Basic_Fluid_Mechanics_by_Peerless/9-analysis_of_dimensional_constant_density_turbulent_flow.ipynb
new file mode 100644
index 0000000..7e6aecc
--- /dev/null
+++ b/Basic_Fluid_Mechanics_by_Peerless/9-analysis_of_dimensional_constant_density_turbulent_flow.ipynb
@@ -0,0 +1,109 @@
+{
+"cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 9: analysis of dimensional constant density turbulent flow"
+ ]
+ },
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 9.1: ex_1.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"n=7\n",
+"w= 62.4 //lbf/ft^3\n",
+"v= 6 //ft/sec\n",
+"d= 2 //in\n",
+"u= 2.34*10^-5 //lbf/ft^3\n",
+"f= 0.0178\n",
+"g= 32.2 //ft/sec^2\n",
+"R= 1.224\n",
+"R1= 8 //ft/sec\n",
+"//CALCULATIONS\n",
+"r= (n+1)*(2*n+1)/(2*n^2)\n",
+"Red= w*v*(d/12)/(u*g)\n",
+"C= (d/Red)^(1/7)*R*(R1/f)^(4/7)\n",
+"V = v*sqrt(f/8)\n",
+"//RESULTS\n",
+"printf ('Vmax/V = %.3f',r)\n",
+"printf ('\n Red = %.2e',Red)\n",
+"printf ('\n C = %.2f',C)\n",
+"printf ('\n Velocity = %.3f ft/sec',V)"
+ ]
+ }
+,
+{
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Example 9.3: ex_3.sce"
+ ]
+ },
+ {
+"cell_type": "code",
+ "execution_count": null,
+ "metadata": {
+ "collapsed": true
+ },
+ "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"clear\n",
+"Re= 5\n",
+"g= 32.2 //ft/sec^2\n",
+"u= 2.34*10^-5 //lbf/ft sec\n",
+"w= 62.4 //lbf/ft^3\n",
+"v= 0.283 //ft/sec\n",
+"Re1= 70\n",
+"v1= 0.0374 //ft/sec\n",
+"//CALCULATIONS\n",
+"y= Re*u*g/(w*v)\n",
+"y1= Re1*u*g/(w*v)\n",
+"y2= Re*u*g/(w*v1)\n",
+"y3= Re1*u*g/(w*v1)\n",
+"//RESULTS\n",
+"printf ('y = %.6f ft',y)\n",
+"printf ('\n y = %.5f ft',y1)\n",
+"printf ('\n y = %.5f ft',y2)\n",
+"printf ('\n y = %.4f ft',y3)"
+ ]
+ }
+],
+"metadata": {
+ "kernelspec": {
+ "display_name": "Scilab",
+ "language": "scilab",
+ "name": "scilab"
+ },
+ "language_info": {
+ "file_extension": ".sce",
+ "help_links": [
+ {
+ "text": "MetaKernel Magics",
+ "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+ }
+ ],
+ "mimetype": "text/x-octave",
+ "name": "scilab",
+ "version": "0.7.1"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
generated by cgit (git 2.34.1) at 2024-12-22 15:27:32 +0000