diff options
Diffstat (limited to 'Turbomachines_by_A_V_Arasu/7-DIMENSIONAL_AND_MODEL_ANALYSIS.ipynb')
| -rw-r--r-- | Turbomachines_by_A_V_Arasu/7-DIMENSIONAL_AND_MODEL_ANALYSIS.ipynb | 357 |
1 files changed, 357 insertions, 0 deletions
diff --git a/Turbomachines_by_A_V_Arasu/7-DIMENSIONAL_AND_MODEL_ANALYSIS.ipynb b/Turbomachines_by_A_V_Arasu/7-DIMENSIONAL_AND_MODEL_ANALYSIS.ipynb new file mode 100644 index 0000000..5a9869a --- /dev/null +++ b/Turbomachines_by_A_V_Arasu/7-DIMENSIONAL_AND_MODEL_ANALYSIS.ipynb @@ -0,0 +1,357 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 7: DIMENSIONAL AND MODEL ANALYSIS" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.10: SPEED_DISCHARGE_AND_POWER.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"DmDp=1/10//The model ratio to prototype\n", +"Hm=5//The head developed by the model in m\n", +"Hp=8.5//The head developed by the prototype in m\n", +"Pp=8000*10^3//The power developed by the prototype in W\n", +"Np=120//The speed of running of the prototype in rpm\n", +"d=1000//density of the water in kg/m^3\n", +"g=9.81//Acceleration due to gravity in m/s^2\n", +"n0=0.85//Overall efficiency of the prototype\n", +"\n", +"//calculations\n", +"Nm=((Hm/Hp)^(1/2))*(1/DmDp)*(Np)//Speed of the mpdel in rpm\n", +"Qp=Pp/(d*g*n0*Hp)//Discharge from the prototype in m^3/s\n", +"Qm=((DmDp)^(3))*(Nm/Np)*(Qp)//Discharge from the model in m^3/s\n", +"Pm=((DmDp)^(5))*((Nm/Np)^(3))*(Pp)*10^-3//Power of the model in kW\n", +"\n", +"//output\n", +"printf('(a)Speed of the model is %3.1f rpm\n(b)Discharge from the model is %3.3f m^3/s\n(c)Power of the model is %3.1f kW',Nm,Qm,Pm)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.11: SPEED_AND_POWER_DEVELOPED.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"P1=6600//Initial power developed by the turbine in kW\n", +"N1=100//Initial speed of the turbine in rpm\n", +"H1=30//Initial head of the turbine in m\n", +"H2=18//Final head of the turbine in m\n", +"\n", +"//calculations\n", +"N2=N1*((H2/H1)^(1/2))//The final speed of the turbine in rpm\n", +"P2=P1*((H2/H1)^(3/2))//The final power developed by the turbine in kW\n", +"\n", +"//output\n", +"printf('(1)The final speed of the turbine is %3.2f rpm\n(2)The final power developed by the turbine is %3i kW',N2,P2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.12: PERFORMANCE_OF_TURBINE.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"H1=25//The initial head on the turbine in m\n", +"N1=200//The initial speed of the turbine in rpm \n", +"Q1=9//The initial discharge of the turbine in m^3/s\n", +"n0=0.9//Overall efficiency of the turbine \n", +"H2=20//The final head on the turbine in m\n", +"d=1000//density of the water in kg/m^3\n", +"g=9.81//Acceleration due to gravity in m/s^2\n", +"\n", +"//calculations\n", +"N2=N1*((H2/H1)^(1/2))//The final speed of the turbine in rpm\n", +"Q2=Q1*((H2/H1)^(1/2))//The final discharge of the turbine in m^3/s\n", +"P1=n0*d*g*Q1*H1*10^-3//Power produced by the turbine initially in kW\n", +"P2=P1*((H2/H1)^(3/2))//Power produced by the turbine finally in kW\n", +"\n", +"//output\n", +"printf('(a)The final speed of the turbine is %3.2f rpm\n(b)The final discharge of the turbine is %3.2f m^3/s\n(c)Power produced by the turbine initially is %3.3f kW\n(d)Power produced by the turbine finally is %3.2f kW',N2,Q2,P1,P2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.13: SPECIFIC_SPEED.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"P1=5000*10^3//The initial power produced in W\n", +"H1=250//The initial head produced in m\n", +"N1=210//The initial speed of turbine in rpm\n", +"n0=0.85//Overall efficiency of the turbine \n", +"H2=160//The final head produced in m\n", +"d=1000//density of the water in kg/m^3\n", +"g=9.81//Acceleration due to gravity in m/s^2\n", +"\n", +"\n", +"//calculations\n", +"Nu=N1/((H1)^(1/2))//The unit speed of the turbine \n", +"Pu=P1/((H1)^(3/2))*10^-3//The unit power of the turbine \n", +"Q1=P1/(d*g*n0*H1)//The initial discharge of the turbine in m^3/s\n", +"Qu=Q1/((H1)^(1/2))//The unit discharge of the turbine \n", +"Q2=Qu*((H2)^(1/2))//The final discharge of the turbine in m^3/s\n", +"N2=Nu*((H2)^(1/2))//The final speed of the turbine in rpm\n", +"P2=Pu*((H2)^(3/2))//The final power of the turbine in kW\n", +"Ns=(N2*((P2)^(1/2)))/((H2)^(5/4))//The specific speed of the turbine\n", +"\n", +"//output\n", +"printf('(a)The unit speed of the turbine is %3.2f\n(b)The unit power of the turbine is %3.3f\n(c)The unit discharge of the turbine is %3.3f\n(d)The final discharge of the turbine is %3.2f m^3/s\n(e)The final speed of the turbine is %3.2f rpm\n(f)The final power of the turbine is %3.1f kW\n(g)The specific speed of the turbine is %3.2f',Nu,Pu,Qu,Q2,N2,P2,Ns)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.5: SPEED_OF_PROTOTYPE.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"Nm=1000//Speed of the model in rpm\n", +"Hm=8//Head of the model in m\n", +"Pm=30//Power of the model in kW\n", +"Hp=25//Head of the prototype in m\n", +"DmDp=1/5//The scale of the model to original\n", +"\n", +"//calculations\n", +"Np=((Hp/Hm)^(1/2))*(DmDp)*(Nm)//Speed of the prototype in rpm\n", +"Pp=(Pm)*((1/DmDp)^(5))*(Np/Nm)^(3)//Power developed by the prototype in kW\n", +"QpQm=((1/DmDp)^(3))*(Np/Nm)//Ratio of the flow rates of two pump(model and prototype)\n", +"\n", +"//output\n", +"printf('(1)Speed of prototype pump is %3.1f rpm\n(2)Power developed by the prototype pump is %3i kW\n(3)Ratio of the flow rates of two pumps is %3.4f',Np,Pp,QpQm)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.6: HEAD_SPEED_AND_SCALE_RATIO.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"Hp=85//Head of the prototype in m\n", +"Qp=(20000/3600)//Flow rate of the prototype in m^3/s\n", +"Np=1490//Speed of the prototype in rpm\n", +"Dp=1.2//Diameter of the prototype in m\n", +"dp=714//Density of the prototype fluid in kg/m^3\n", +"Pp=4//Power of the prototype in MW\n", +"Pm=500*10^-3//Power of the model in MW\n", +"Qm=0.5//Flow rate of the prototype in m^3/s\n", +"dm=1000//Density of the model fluid (water) in kg/m^3\n", +"\n", +"//calculations\n", +"NpNm=(Qp/Qm)//Ratio of the speeds of the prototype and the model in terms of (Dm/Dp)^(3)\n", +"DmDp=1/(((NpNm)^(3))*(dp/dm)*(Pm/Pp))^(1/4)//The ratio of the diameters of model and the prototype or the scale ratio \n", +"NmNp=1/(NpNm*((DmDp)^(3)))//The speed ratio or the ratio of speeds of the model and the prototype\n", +"HmHp=((1/NmNp)^(2))*((1/DmDp)^(2))//The head ratio or the ratio of heads of the model and the prototype \n", +"\n", +"//output\n", +"printf('(1)The head ratio of the model is %3.1f\n(2)The speed ratio of the model is %3.1f\n(3)The scale ratio of the model is %3.1f',HmHp,NmNp,DmDp)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.7: SPEED_AND_DISCHARGE_OF_THE_MODEL.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"Np=400//The speed of the prototype in rpm\n", +"Qp=1.7//The discharge of the prototype in m^3/s\n", +"Hp=36.5//The head of the prototype in m\n", +"Pp=720//The power input of the prototype in kW\n", +"Hm=9//The head of the model in m\n", +"DmDp=1/6//The scale of model to prototype \n", +"\n", +"//calculations\n", +"Nm=((Hm/Hp)^(1/2))*(1/DmDp)*Np//Speed of the model in rpm\n", +"Qm=((DmDp)^(3))*(Nm/Np)*(Qp)//Discharge of the model in m^3/s\n", +"Pm=((DmDp)^(5))*((Nm/Np)^(3))*Pp//Power required by the model in kW\n", +"\n", +"//output\n", +"printf('(a)Speed of the model is %3.2f rpm\n(b)Discharge of the model is %3.4f m^3/s\n(c)Power required by the model is %3.2f kW',Nm,Qm,Pm)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.8: IMPELLER_DIAMETER_OF_PUMP2.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"N1=1000//The running speed of the pump-1 in rpm\n", +"D1=0.3//The impeller diameter of pump-1 in m\n", +"Q1=0.02//The discharge of pump-1 in m^3/s\n", +"H1=15//The head developed by the pump-1 in m\n", +"N2=1000//The running speed of the pump-2 in rpm\n", +"Q2=0.01//The discharge of pump-2 in m^3/s\n", +"\n", +"//calculations\n", +"D2=(((Q2/Q1)*(N1/N2))^(1/3))*(D1)//Impeller diameter of the pump-2 in m\n", +"H2=(((D2/D1)*(N2/N1))^(2))*(H1)//Head developed by the pump-2 in m\n", +"\n", +"//output\n", +"printf('(a)Impeller diameter of the pump-2 is %3.3f m\n(b)Head developed by the pump-2 is %3.2f m',D2,H2)" + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 7.9: SPECIFIC_SPEEDS.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"clc\n", +"clear\n", +"//input data\n", +"DmDp=1/10//The model ratio to prototype \n", +"Pm=1.84//Power developed by the model in kW\n", +"Hm=5//Head developed by the model in m\n", +"Nm=480//Speed of the model in rpm\n", +"Hp=40//Head developed by the prototype in m\n", +"\n", +"//calculations\n", +"Np=((Hp/Hm)^(1/2))*(DmDp)*(Nm)//Speed of the prototype in rpm\n", +"Pp=((1/DmDp)^(5))*((Np/Nm)^(3))*Pm//Power developed by the prototype in kW\n", +"Nsp=((Np*((Pp)^(1/2)))/((Hp)^(5/4)))//Specific speed of the prototype\n", +"Nsm=((Nm*((Pm)^(1/2)))/((Hm)^(5/4)))//Specific speed of the prototype\n", +"\n", +"//output\n", +"printf('(a)Power developed by the prototype is %3i kW\n(b)Speed of the prototype is %3.2f rpm\n(c)Specific speed of the prototype is %3.1f\n(d)Specific speed of the model is %3.1f\n Thus the specific speed of the model is equal to the prototype and thus it is verified',Pp,Np,Nsp,Nsm)" + ] + } +], +"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 +} |