diff options
Diffstat (limited to 'mechanics_of_fluid/Chapter13_1.ipynb')
| -rwxr-xr-x | mechanics_of_fluid/Chapter13_1.ipynb | 166 |
1 files changed, 166 insertions, 0 deletions
diff --git a/mechanics_of_fluid/Chapter13_1.ipynb b/mechanics_of_fluid/Chapter13_1.ipynb new file mode 100755 index 00000000..4bd1ca7e --- /dev/null +++ b/mechanics_of_fluid/Chapter13_1.ipynb @@ -0,0 +1,166 @@ +{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:6ab6ca7ec50b12b4ad5664b1412855d8b1bafd30fee8faf6c1773c772ff215d4"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Chapter13-Fluid Machines"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex1-pg618"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate Hydraulic Efficiency and Overall Efficiency and Outlet angles of the guide vanes and Rotor blade angle at inlet and oulet\n",
+ "\n",
+ "## Maximum hydraulic efficiency occurs for minimum pressure loss, that is, when\n",
+ "\n",
+ "## dp1/dQ=2.38Q-1.43=0\n",
+ "\n",
+ "Q_opt=1.43/2.38;\n",
+ "\n",
+ "p1_min=1.19*Q_opt**2-1.43*Q_opt+0.47; ## MPa\n",
+ "\n",
+ "rho=1000.; ## kg/m**3\n",
+ "g=9.81; ## m/s**2\n",
+ "w=69.1; ## rad/s\n",
+ "P=200.*10.**3.; ## W\n",
+ "Ohm_P=0.565; ## rad\n",
+ "d=0.5; ## m\n",
+ "h=0.06; ## m\n",
+ "\n",
+ "p1=p1_min*10.**6./(rho*g); ## mH2O, coversion of units\n",
+ "\n",
+ "H=(w*P**(1/2.)/(rho**(1/2.)*Ohm_P))**(4/5.)/g;\n",
+ "\n",
+ "Hydraulic_efficiency=(H-p1)/H;\n",
+ "print'%s %.3f %s'%(\"Hydraulic Efficiency =\",Hydraulic_efficiency,\"\")\n",
+ "\n",
+ "\n",
+ "Overall_efficiency=P/(Q_opt*rho*g*H);\n",
+ "print'%s %.3f %s'%(\"Overall Efficiency =\",Overall_efficiency,\"\")\n",
+ "\n",
+ "\n",
+ "H_Euler=H-p1;\n",
+ "\n",
+ "u1=w*0.25;\n",
+ "v_w1=g*H_Euler/u1;\n",
+ "A=math.pi*d*h*0.95; \n",
+ "v_r=Q_opt/A;\n",
+ "\n",
+ "alpha1=math.atan(v_r/v_w1);\n",
+ "print'%s %.3f %s'%(\"Outlet angles of the guide vanes =\",alpha1,\"degrees\")\n",
+ "\n",
+ "beta1=math.atan(v_r/(v_w1-u1));\n",
+ "print'%s %.2f %s'%(\"Rotor blade angle at inlet =\",beta1,\"degrees\")\n",
+ "\n",
+ "u2=w*0.325/2;\n",
+ "beta2=math.atan(v_r/u2);\n",
+ "print'%s %.3f %s'%(\"Rotor blade angle at outlet =\",beta2,\"degrees\")\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Hydraulic Efficiency = 0.896 \n",
+ "Overall Efficiency = 0.855 \n",
+ "Outlet angles of the guide vanes = 0.321 degrees\n",
+ "Rotor blade angle at inlet = 1.16 degrees\n",
+ "Rotor blade angle at outlet = 0.539 degrees\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex2-pg622"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "import math\n",
+ "#calculate Overall efficiency and Limiting value for the height of the draft tube above\n",
+ "w=6.25;\n",
+ "D=0.75; ## m\n",
+ "gv_angle=15; ## guide vane angle in degrees\n",
+ "g=9.81; ## m/s^2\n",
+ "H=27.5; ## m\n",
+ "A1=0.2; ## m^2\n",
+ "rho=1000.; ## kg/m^3\n",
+ "p_atm=101.3*10**3.;\n",
+ "p_min=35.*10.**3.;\n",
+ "\n",
+ "u1=math.pi*w*D;\n",
+ "v1=u1*math.sin(105.)/math.sin(60.);\n",
+ "v_r1=v1*math.sin(gv_angle);\n",
+ "v_w1=v1*math.cos(gv_angle);\n",
+ "v_w2=0.;\n",
+ "\n",
+ "n_hydraulic=u1*v_w1/g/H;\n",
+ "\n",
+ "n_overall=0.97*n_hydraulic;\n",
+ "print'%s %.1f %s'%(\"Overall efficiency =\",n_overall,\"\")\n",
+ "\n",
+ "\n",
+ "Q=A1*v_r1;\n",
+ "\n",
+ "P=n_overall*Q*rho*g*H;\n",
+ "Ohm_P=w*2.*math.pi/(g*H)**(5/4)*(P/rho)**(1/2);\n",
+ "\n",
+ "## sigma > 0.119*(0.5)^(1.84) = 0.0331\n",
+ "\n",
+ "sigma=0.0331;\n",
+ "\n",
+ "##((p_atm-p_min)/(rho*g)-z0)/H > 0.0331\n",
+ "\n",
+ "z0=((p_atm-p_min)/(rho*g))-sigma*H;\n",
+ "print'%s %.2f %s'%(\"Limiting value for the height of the draft tube above =\",z0,\"m\")\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Overall efficiency = -1.9 \n",
+ "Limiting value for the height of the draft tube above = 5.85 m\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
\ No newline at end of file |