{
 "metadata": {
  "name": "",
  "signature": "sha256:5804a8cf03ec2dc9fe48e7282fc5bb2ad9c57bb262c9f4fea9720f853cbc6882"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Ch-9, Nuclear Power Stations"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.1 Page 169"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "m=1*10**-3#mass of 1 grm in kgs\n",
      "c=3*10**8\n",
      "e=m*c**2 \n",
      "E=e/(1000*3600)\n",
      "print \"energy equivalent of 1 gram is %dkWh\"%E"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "energy equivalent of 1 gram is 25000000kWh\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.2 Page 169"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "amu=1.66*10**-27#mass equvalent in kgs\n",
      "c=3*10**8\n",
      "j=6.242*10**12\n",
      "e=amu*c**2\n",
      "E=e*j \n",
      "print \" energy evalent in joules is %ejoules \\n energy equvalent in Mev is %dMeV \\n hense shown\"%(e,E)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " energy evalent in joules is 1.494000e-10joules \n",
        " energy equvalent in Mev is 932MeV \n",
        " hense shown\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.3 Page 169"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "hm=2.0141\n",
      "hp=1.007825\n",
      "hn=1.008665\n",
      "nm=58.9342\n",
      "np=28\n",
      "nn=59\n",
      "um=235.0439\n",
      "up=92\n",
      "un=235\n",
      "hmd=hp+hn-hm ;nmd=np*hp+(nn-np)*hn-nm ;umd=up*hp+(un-up)*hn-um \n",
      "hbe=931*hmd ;nbe=931*nmd; ube=931*umd \n",
      "ahbe=hbe/2 ;anbe=nbe/nn ;aube=ube/un \n",
      "print \"\\t(a)\\n mass defect is for hydrogen %famu \\n total binding energy for hydrogens %fMev \\n average binding energy for hydrogen is %fMeV\"%(hmd,hbe,ahbe)\n",
      "print \"\\n\\t(b)\\n mass defect is for nickel %famu \\n total binding energy for nickel is %fMev \\n average binding energy for nickelis %fMeV\"%(nmd,nbe,anbe)\n",
      "print \"\\n\\t(c)\\n mass defect of uranium is %famu \\n total binding energy uranium is %fMev \\n average binding energy uranium is %fMeV\"%(umd,ube,aube)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\t(a)\n",
        " mass defect is for hydrogen 0.002390amu \n",
        " total binding energy for hydrogens 2.225090Mev \n",
        " average binding energy for hydrogen is 1.112545MeV\n",
        "\n",
        "\t(b)\n",
        " mass defect is for nickel 0.553515amu \n",
        " total binding energy for nickel is 515.322465Mev \n",
        " average binding energy for nickelis 8.734279MeV\n",
        "\n",
        "\t(c)\n",
        " mass defect of uranium is 1.915095amu \n",
        " total binding energy uranium is 1782.953445Mev \n",
        " average binding energy uranium is 7.587036MeV\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.4 Page 170"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from math import exp\n",
      "no=1.7*10**24\n",
      "hl=7.1*10**8\n",
      "t=10*10**8\n",
      "lm=0.693/(hl)\n",
      "lmda=lm/(8760*3600)\n",
      "ia=lmda*no\n",
      "n=no*(exp(-lm*t))\n",
      "print \"(lamda) disintegrations per sec is %ebq \\n initial activity is lamda*na is %ebq \\n final number of atoms is %eatoms\"%(lmda,ia,n)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "(lamda) disintegrations per sec is 3.095054e-17bq \n",
        " initial activity is lamda*na is 5.261592e+07bq \n",
        " final number of atoms is 6.405500e+23atoms\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.5 Page 170"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "um=5\n",
      "owp=2.6784*10**15\n",
      "an=6.023*10**23\n",
      "na1g=an/235\n",
      "na5g=an*5/235\n",
      "p=na5g/owp\n",
      "print \" 1 watt power requvires %efussions per day \\n number of atoms in 5 gram is %eatoms \\n power is %eMW \"%(owp,na5g,p)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " 1 watt power requvires 2.678400e+15fussions per day \n",
        " number of atoms in 5 gram is 1.281489e+22atoms \n",
        " power is 4.784533e+06MW \n"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.6 Page 171"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "pp=235\n",
      "pe=0.33\n",
      "lf=1\n",
      "teo=pp*8760*3600*10**6\n",
      "ei=teo/pe\n",
      "nfr=3.1*10**10#fessions required\n",
      "tnfr=nfr*ei\n",
      "t1gu=2.563*10**21 #total uranium atoms in 1 grm\n",
      "fure=tnfr/t1gu\n",
      "print \"total energy input %eWatt sec \\n energy input is %eWatt-sec\\n total number of fissions required is %efissions \\n fuel required is %e grams %dkg\"%(teo,ei,tnfr,fure,fure/1000)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "total energy input 7.410960e+15Watt sec \n",
        " energy input is 2.245745e+16Watt-sec\n",
        " total number of fissions required is 6.961811e+26fissions \n",
        " fuel required is 2.716274e+05 grams 271kg\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "example 9.7 Page 171"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "from math import log\n",
      "en=3*10**6\n",
      "a=12\n",
      "fen=0.1\n",
      "Es=2/(12+2/3)\n",
      "re=exp(Es)\n",
      "print \"(a)\\nratio of energies per collision is %f\"%(re)\n",
      "rietf=en/fen\n",
      "ldie=log(rietf)\n",
      "nc=ldie/Es\n",
      "print \"(b)\\npatio of iniial to final energies is %e \\nlogarithemic decrement in energy is %f \\nnumber of collisions is %d\"%(rietf,ldie,nc)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "(a)\n",
        "ratio of energies per collision is 1.171043\n",
        "(b)\n",
        "patio of iniial to final energies is 3.000000e+07 \n",
        "logarithemic decrement in energy is 17.216708 \n",
        "number of collisions is 109\n"
       ]
      }
     ],
     "prompt_number": 16
    }
   ],
   "metadata": {}
  }
 ]
}