{
 "metadata": {
  "name": "",
  "signature": "sha256:2259086c8cee958c316a471676c823626713e37889548ae1aff2070252fca121"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter12:Power System Transients"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E1 - Pg 30"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "rc=.5e-2;\n",
      "rs=1.5e-2;\n",
      "u=4\n",
      "\n",
      "L=2e-7 * math.log(rs/rc);\n",
      "print '%s %.2f %s' %(\"\\nL=\",L*1e7,\"e-7H/m\")\n",
      "C=u*1e-9/(18 * math.log(rs/rc))\n",
      "print '%s %.2f %s' %(\"\\nC=\", C*1e9,\"e-9F/m\")\n",
      "v=1/math.sqrt(L*C);\n",
      "print '%s %.2f %s' %(\"\\nv=\", v*1e-8,\"e8m/s\")\n",
      "Zc=math.sqrt(L/C)\n",
      "print '%s %.2f %s' %(\"\\nZc=\",Zc,\"ohm\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "L= 2.20 e-7H/m\n",
        "\n",
        "C= 0.20 e-9F/m\n",
        "\n",
        "v= 1.50 e8m/s\n",
        "\n",
        "Zc= 32.96 ohm\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E2 - Pg 32"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "ef=100.;\n",
      "Zc=400.;\n",
      "Z=50.;\n",
      "et=2.*ef*Z/(Z+Zc)\n",
      "print '%s %.2f %s' %(\"Surge transmitted=\",et,\"kV\\n\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Surge transmitted= 22.22 kV\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E3 - Pg 32"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "ef=200.;\n",
      "Zc=400.;\n",
      "Z1=500.;\n",
      "Z2=300.;\n",
      "et=2.*ef*(Z1*Z2/(Z1+Z2))/((Z1*Z2/(Z1+Z2))+Zc)\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage transmitted=\",et,\"kV\\n\")\n",
      "it1=et/Z1;\n",
      "print '%s %.2f %s' %(\"\\nSurge Current transmitted=\",it1,\"kA\\n\")\n",
      "it2=et/Z2;\n",
      "print '%s %.2f %s' %(\"\\nSurge Current transmitted=\",(it2*100.)/100.,\"kA\\n\")\n",
      "er=et-ef;\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage Reflected=\",er,\"kV\\n\")\n",
      "ir=it1+it2-(ef/Zc)\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Reflected=\",ir,\"kA\\n\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "Surge Voltage transmitted= 127.66 kV\n",
        "\n",
        "\n",
        "Surge Current transmitted= 0.26 kA\n",
        "\n",
        "\n",
        "Surge Current transmitted= 0.43 kA\n",
        "\n",
        "\n",
        "Surge Voltage Reflected= -72.34 kV\n",
        "\n",
        "\n",
        "Surge Current Reflected= 0.18 kA\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E4 - Pg 35"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "E=100.\n",
      "Zc=400.\n",
      "L=4000.\n",
      "\n",
      "print '%s %.2f %s %.2f' %(\"et=\",2*E,\"math.exp( - %.1f t) KV\\n\",Zc/L)\n",
      "print '%s %.2f %s %.2f' %(\"er=\",E,\"(2*math.exp( - %.1f t) -1) KV\\n\",Zc/L)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "et= 200.00 math.exp( - %.1f t) KV\n",
        " 0.10\n",
        "er= 100.00 (2*math.exp( - %.1f t) -1) KV\n",
        " 0.10\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E7 - Pg 36"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "V=300e3\n",
      "R=400.\n",
      "k=1.5e-27\n",
      "\n",
      "E=10.\n",
      "x=1.\n",
      "e=1e-5\n",
      "while (E>e):\n",
      "    f=(k*R*x**6.) +x -(2.*V)\n",
      "    df=(6.* k*R*x**5.) +1.\n",
      "    x1=x-(f/df)\n",
      "    E=abs(x1-x)\n",
      "    x=x1\n",
      "eA=round(x)\n",
      "IA=k*eA**6.\n",
      "\n",
      "print '%s %.2f %s %.2f' %(\"eA=\",eA,\"Ia=\",IA)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "eA= 97101.00 Ia= 1257.28\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E8 - Pg 36"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "V=300e3\n",
      "R1=400.\n",
      "R2=50.\n",
      "R=1.+(400./50.)\n",
      "k=1.5e-27\n",
      "\n",
      "E=10.\n",
      "x=1.\n",
      "e=1e-5\n",
      "while (E>e) :\n",
      "    f=(k*R1*x**6.) +(R*x) -(2.*V)\n",
      "    df=(6.* k*R1*x**5.) +R\n",
      "    x1=x-(f/df)\n",
      "    E=abs(x1-x)\n",
      "    x=x1\n",
      "eA=round(x)\n",
      "IA=k*eA**6.\n",
      "\n",
      "print '%s %.2f %s %.2f' %(\"eA=\",eA,\"Ia=\",IA)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "eA= 62640.00 Ia= 90.62\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E9 - Pg 45"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "ef=3000.;\n",
      "Zc=300.;\n",
      "ea=1700.;\n",
      "iF=ef/Zc\n",
      "print '%s %.2f %s' %(\"\\nCurrent in line=\",iF,\"kA\\n\")\n",
      "Ia=((2.*ef)-ea)/Zc\n",
      "print '%s %.2f %s' %(\"\\nCurrent through Arrester=\",Ia,\"kA\\n\")\n",
      "Ia=round(Ia *1000.)/1000.\n",
      "R=ea/Ia\n",
      "print '%s %.2f %s' %(\"\\nresistance of arrester=\",R,\"ohm\\n\")\n",
      "er=ea-ef;\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage Reflected=\",er,\"kV\\n\")\n",
      "Cr=er/ef;\n",
      "CR=ea/ef;\n",
      "print '%s %.2f %s %.2f' %(\"\\nCoeff of Reflection =\",Cr,\"Coeff of Refraction=\",CR)\n",
      "Cr=(R-Zc)/(R+Zc);\n",
      "CR=(R*2)/(R+Zc);\n",
      "print '%s %.2f %s %.2f' %(\"\\nVerification: Coeff of Reflection =\",Cr,\"Coeff of Refraction=\",CR)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "Current in line= 10.00 kA\n",
        "\n",
        "\n",
        "Current through Arrester= 14.33 kA\n",
        "\n",
        "\n",
        "resistance of arrester= 118.61 ohm\n",
        "\n",
        "\n",
        "Surge Voltage Reflected= -1300.00 kV\n",
        "\n",
        "\n",
        "Coeff of Reflection = -0.43 Coeff of Refraction= 0.57\n",
        "\n",
        "Verification: Coeff of Reflection = -0.43 Coeff of Refraction= 0.57\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E10 - Pg 45"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "ef=10000.;\n",
      "Zc=400.;\n",
      "iF=ef/Zc\n",
      "print '%s %.2f' %(\"\\n(a)\\nIncident Wave magnitude= A\",iF)\n",
      "\n",
      "R=1000\n",
      "et=ef*(R*2)/(R+Zc);\n",
      "it=et/R;\n",
      "er=et-ef;\n",
      "print '%s %.2f' %(\"\\n(b)\\nSurge Voltage Reflected= KV\",er/1000)\n",
      "ir=-1*er/Zc\n",
      "print '%s %.2f' %(\"\\nSurge Current Reflected= A\",ir)\n",
      "edr=et*it;\n",
      "print '%s %.2f' %(\"\\nRate of dissipation of energy= KW\",edr/1000)\n",
      "err=er*-ir;\n",
      "print '%s %.2f' %(\"\\nRate of reflection of energy= KW\",err/1000)\n",
      "\n",
      "print '%s %.2f' %(\"\\n(c)\\nfor complete dissipation, R=Zc= ohm\",Zc);\n",
      "\n",
      "R=50\n",
      "et=ef*(R*2)/(R+Zc);\n",
      "print '%s %.2f' %(\"\\n(d)\\nSurge Voltage Transmitted= KV\",et/1000)\n",
      "it=et/R;\n",
      "print '%s %.2f' %(\"\\nSurge Current Transmitted= A\",it)\n",
      "er=et-ef;\n",
      "print '%s %.2f' %(\"\\nSurge Voltage Reflected= kV\",er/1000)\n",
      "ir=-1*er/Zc\n",
      "print '%s %.2f' %(\"\\nSurge Current Reflected= A\",ir)\n",
      "edr=et*it;\n",
      "print '%s %.2f' %(\"\\nRate of dissipation of energy= KW\",edr/1000)\n",
      "err=er*-ir;\n",
      "print '%s %.2f' %(\"\\nRate of reflection of energy= KW\",err/1000)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "(a)\n",
        "Incident Wave magnitude= A 25.00\n",
        "\n",
        "(b)\n",
        "Surge Voltage Reflected= KV 4.29\n",
        "\n",
        "Surge Current Reflected= A -10.71\n",
        "\n",
        "Rate of dissipation of energy= KW 204.08\n",
        "\n",
        "Rate of reflection of energy= KW 45.92\n",
        "\n",
        "(c)\n",
        "for complete dissipation, R=Zc= ohm 400.00\n",
        "\n",
        "(d)\n",
        "Surge Voltage Transmitted= KV 2.22\n",
        "\n",
        "Surge Current Transmitted= A 44.44\n",
        "\n",
        "Surge Voltage Reflected= kV -7.78\n",
        "\n",
        "Surge Current Reflected= A 19.44\n",
        "\n",
        "Rate of dissipation of energy= KW 98.77\n",
        "\n",
        "Rate of reflection of energy= KW 151.23\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E11 - Pg 46"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "Zc=400.\n",
      "ef=20.\n",
      "z1=150.;\n",
      "z2=200.\n",
      "z=round((z1*z2/(z1+z2))*100.)/100.\n",
      "\n",
      "et=2*ef*z/(Zc+z)\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage Transmitted=\",et,\"kV\\n\")\n",
      "\n",
      "it1=et*1000./z1;\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Transmitted in line 1=\",it1,\"A\\n\")\n",
      "\n",
      "it2=et*1000./z2;\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Transmitted in line 2=\",it2,\"A\\n\")\n",
      "\n",
      "er=et-ef\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage Reflected=\",er,\"kV\\n\")\n",
      "ir=-1*er*1000./Zc\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Reflected=\",ir,\"A\\n\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "Surge Voltage Transmitted= 7.06 kV\n",
        "\n",
        "\n",
        "Surge Current Transmitted in line 1= 47.06 A\n",
        "\n",
        "\n",
        "Surge Current Transmitted in line 2= 35.29 A\n",
        "\n",
        "\n",
        "Surge Voltage Reflected= -12.94 kV\n",
        "\n",
        "\n",
        "Surge Current Reflected= 32.35 A\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E12 - Pg 49"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "ef=100.\n",
      "Zc=400.\n",
      "z=50.\n",
      "R=z+Zc;\n",
      "E=(2.*ef/(Zc+z+R))**2. *R\n",
      "E=round(E*100.)/100.\n",
      "print '%s %.2f %s %s %.2f %s' %(\"\\n(a)Energy transfer max when R=\",R,\"ohm\\n\",\"energy=\",E,\"KW\\n\");\n",
      "etB=2.*ef*z/(z+Zc+R);\n",
      "etB=round(etB*100.)/100.\n",
      "print '%s %.2f %s' %(\"\\n(b)Surge Voltage Transmitted=\",etB,\"kV\\n\")\n",
      "it=etB*1000./z;\n",
      "it=round(it*100.)/100.\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Transmitted =\",it,\"A\\n\")\n",
      "etA=2.*ef*(z+R)/(z+Zc+R);\n",
      "etA=round(etA*100.)/100.\n",
      "erA=etA-ef\n",
      "print '%s %.2f %s' %(\"\\n(c)Surge Voltage Reflected=\",erA,\"kV\\n\")\n",
      "irA=-1.*erA*1000./Zc\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Reflected=\",irA,\"A\\n\")\n",
      "iF=ef*1000./Zc\n",
      "Pi=ef*iF\n",
      "print '%s %.2f %s' %(\"\\n(d)Power Incident=\",Pi,\"kW\\n\")\n",
      "Pr=erA*-irA\n",
      "print '%s %.2f %s' %(\"\\nPower Reflected=\",Pr,\"kW\\n\")\n",
      "Pt=erA*it\n",
      "print '%s %.2f %s' %(\"\\nPower Transmitted=\",Pt,\"kW\\n\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "(a)Energy transfer max when R= 450.00 ohm\n",
        " energy= 22.22 KW\n",
        "\n",
        "\n",
        "(b)Surge Voltage Transmitted= 11.11 kV\n",
        "\n",
        "\n",
        "Surge Current Transmitted = 222.20 A\n",
        "\n",
        "\n",
        "(c)Surge Voltage Reflected= 11.11 kV\n",
        "\n",
        "\n",
        "Surge Current Reflected= -27.77 A\n",
        "\n",
        "\n",
        "(d)Power Incident= 25000.00 kW\n",
        "\n",
        "\n",
        "Power Reflected= 308.58 kW\n",
        "\n",
        "\n",
        "Power Transmitted= 2468.64 kW\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E14 - Pg 50"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "I=5.;\n",
      "z1=400.\n",
      "z2=50.\n",
      "V=I * z1* z2/(z2+z1)\n",
      "print '%s %.2f %s' %(\"\\nSurge Voltage Transmitted=\",V,\"kV\\n\")\n",
      "ic=V/z2\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Transmitted in cable=\",ic,\"kA\\n\")\n",
      "\n",
      "io=-V/z1;\n",
      "print '%s %.2f %s' %(\"\\nSurge Current Transmitted in OH line=\",io,\"kA\\n\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "Surge Voltage Transmitted= 222.22 kV\n",
        "\n",
        "\n",
        "Surge Current Transmitted in cable= 4.44 kA\n",
        "\n",
        "\n",
        "Surge Current Transmitted in OH line= -0.56 kA\n",
        "\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example E16 - Pg 51"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "i=100.\n",
      "L=4e-3\n",
      "C=300e-12\n",
      "E=i* math.sqrt(L/C)\n",
      "T=1./ math.sqrt(L*C)\n",
      "print '%s %.2f %s %.2f' %(\"e=\",E/1e3,\"*1e3 sin( %.3f *1e6 t) kV\",T/1e6)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "e= 365.15 *1e3 sin( %.3f *1e6 t) kV 0.91\n"
       ]
      }
     ],
     "prompt_number": 12
    }
   ],
   "metadata": {}
  }
 ]
}