summaryrefslogtreecommitdiff
path: root/Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb
diff options
context:
space:
mode:
authorTrupti Kini2016-04-12 23:30:26 +0600
committerTrupti Kini2016-04-12 23:30:26 +0600
commitbde93a104da9d4510f751a003cdf51ddaa45da1c (patch)
tree033e6b60c456cc12fdda7e3317851624e28fabb5 /Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb
parent999a3252fe20da210975277eb235d4b132592ab6 (diff)
downloadPython-Textbook-Companions-bde93a104da9d4510f751a003cdf51ddaa45da1c.tar.gz
Python-Textbook-Companions-bde93a104da9d4510f751a003cdf51ddaa45da1c.tar.bz2
Python-Textbook-Companions-bde93a104da9d4510f751a003cdf51ddaa45da1c.zip
Added(A)/Deleted(D) following books
A Power_Plant_Engineering_by_P._K._Nag/Ch1.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch10.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch11.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch12.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch2.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch3.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch4.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch5.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch6.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch7.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch8.ipynb A Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb A Power_Plant_Engineering_by_P._K._Nag/screenshots/LoadCurve.png A Power_Plant_Engineering_by_P._K._Nag/screenshots/LoadCurve2.png A Power_Plant_Engineering_by_P._K._Nag/screenshots/LoadDurationCurve.png A sample_notebooks/FameethaAlamBasha/Chapter7_1.ipynb A sample_notebooks/Hrituraj/Ch-6_1.ipynb A sample_notebooks/NishthaRani/Chapter12.ipynb A sample_notebooks/SaleemAhmed/Chapter10_1.ipynb A sample_notebooks/VivekMaindola/chap3.ipynb
Diffstat (limited to 'Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb')
-rw-r--r--Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb474
1 files changed, 474 insertions, 0 deletions
diff --git a/Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb b/Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb
new file mode 100644
index 00000000..b7e3e772
--- /dev/null
+++ b/Power_Plant_Engineering_by_P._K._Nag/Ch9.ipynb
@@ -0,0 +1,474 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter 9 : Nuclear Power Plant"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.1 Pg: 648"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 2,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ " The mass defect = 0.13703 amu \n",
+ " The binding energy per nucleon = 7.97 MeV \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "mp=1.007277##Atomic Mass of proton in amu\n",
+ "mn=1.008665##Atomic Mass of neutron in amu\n",
+ "me=0.00055##Atomic Mass of electron in amu\n",
+ "mo=15.99491##Atomic Mass of oxygen in amu\n",
+ "np=8##Number of protons in oxygen\n",
+ "ne=8##Number of electrons in oxygen\n",
+ "nn=8##Number of neutrons in oxygen\n",
+ "a=931##One amu in MeV\n",
+ "No=16##Number of nucleons in oxygen\n",
+ "\n",
+ "#Calculations\n",
+ "m=(np*mp)+(ne*me)+(nn*mn)-mo##The mass defect in amu\n",
+ "B=m*a##Binding energy in MeV\n",
+ "Bn=B/No##Binding energy per nucleon\n",
+ "\n",
+ "#Output\n",
+ "print \" The mass defect = %3.5f amu \\n The binding energy per nucleon = %3.2f MeV \"%(m,Bn)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.2 Pg: 649"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 24,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The decay constant = 1.36e-11 s**-1\n",
+ "The initial activity of 1 g of radium 226 = 3.61e+10 dis/s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "amr=226.095##Atomic mass of radium in amu\n",
+ "AC=6.023*10**23##Avogadro constant in molecules/g.mol\n",
+ "h=1620##Half life of radium in years\n",
+ "\n",
+ "#Calculations\n",
+ "D=(0.6931/(h*365*24*3600))##The decay constant in 1/s\n",
+ "Na=AC/amr##Number of atoms per gram of radium \n",
+ "Ao=D*Na##Initial activity in dis/s\n",
+ "\n",
+ "#Output\n",
+ "print \"The decay constant = %0.2e s**-1\"%D\n",
+ "print \"The initial activity of 1 g of radium 226 = %0.2e dis/s\"%Ao,"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.3 Pg: 649"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 25,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The fuel consumed of U-235 per day = 3.9 g/day \n"
+ ]
+ }
+ ],
+ "source": [
+ "from __future__ import division\n",
+ "#Input data\n",
+ "F=190##Each fission of U-235 yeilds in MeV\n",
+ "a=85##Assuming the Neutrons absorbed by U-235 cause fission in percentage\n",
+ "b=15##Non fission capture to produce an isotope U-236 in percentage\n",
+ "Q=3000##The amount of thermal power produced in MW\n",
+ "\n",
+ "#Calculations\n",
+ "E=F*1.60*10**-13##Each fission yields a useful energy in J\n",
+ "N=1/E##Number of fissions required \n",
+ "B=((10**6)*(N*86400))/(a/100)##One day operation of a reactor the number of U-235 nuclei burned is in absorptions per day\n",
+ "M=(B*235)/(6.023*10**23)##Mass of U-235 consumed to produce one MW power in g/day\n",
+ "M1=M*3##Mass of U-235 consumed to produce 3000 MW power in g/day\n",
+ "\n",
+ "#Output\n",
+ "print \"The fuel consumed of U-235 per day = %3.1f g/day \"%(M1)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.4 Pg: 650"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 9,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The cross section of neutrons = 5.02 barns \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "sa1=10##Cross section of nucleus in barns\n",
+ "N=2200##Neutrons in m/s\n",
+ "En1=0.1##Kinetic energy of neutrons increases in eV\n",
+ "En2=0.02525##Kinetic energy of neutron in eV\n",
+ "\n",
+ "#Calculations\n",
+ "sa2=sa1/((En1/En2)**0.5)##The cross section of neutrons in barns\n",
+ "\n",
+ "#Output\n",
+ "print \"The cross section of neutrons = %3.2f barns \"%(sa2)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.5 Pg: 650"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The microscopic absorption cross section of natural uranium = 7.6 barns \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "U1=99.285##Uranium consists of U-238 in percentage \n",
+ "U2=0.715##Uranium consists of U-235 in Percentage\n",
+ "E=0.025##The energy of neutrons in eV\n",
+ "sc=2.72##Capture cross section for U-238 in barns\n",
+ "sf=0##fission cross section for U-238 in barns\n",
+ "sc1=101##Capture cross section for U-235 in barns\n",
+ "sf1=579##fission cross section for U-235 in barns\n",
+ "\n",
+ "#Calculations\n",
+ "sa=(U1/100)*(sc+sf)+(U2/100)*(sc1+sf1)##The microscopic absorption cross section of natural uranium in barns\n",
+ "\n",
+ "#Output\n",
+ "print \"The microscopic absorption cross section of natural uranium = %3.1f barns \"%(sa)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.6 Pg: 650"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 12,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The microscopic capture cross section of water = 0.0222 cm**-1 \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "p=1##The density of water in g/cm**3\n",
+ "sch=0.332##The microscopic capture cross section of hydrogen in barn\n",
+ "sco=0.0002##The microscopic capture cross section of oxygen in barn\n",
+ "\n",
+ "#Calculations\n",
+ "N=(6.023*10**23)*p/18##Number of molecules of water per cm**3\n",
+ "scw=(2*N*sch*10**-24)+(N*sco*10**-24)##The microscopic capture cross section of water in cm**-1\n",
+ "\n",
+ "#output\n",
+ "print \"The microscopic capture cross section of water = %3.4f cm**-1 \"%(scw)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.7 Pg: 650"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 15,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The thermal neutron flux = 9.10e+11 cm**-2s**-1\n",
+ "The average distance that a neutron travels before it is absorbed = 0.010 cm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "m=230##The amount of boron piece in g\n",
+ "mw=10##The molecular weight of boron \n",
+ "R=9.57*10**13##Reaction rate in cm**-3s**-1\n",
+ "d=2.3##Density of boron in g/cm**3\n",
+ "sa=755##Absorbption cross section in barns\n",
+ "ss=4##Elastic scattering cross section in barns\n",
+ "\n",
+ "#Calculations\n",
+ "st=sa+ss##The total cross section in barns\n",
+ "N=(d*6.023*10**23)/mw##The number density of neutrons in cm**-3\n",
+ "S=N*st*10**-24##Number density of neutrons for total in cm**-1\n",
+ "F=R/S##Neutron flux in cm**-2s**-1\n",
+ "L=1/S##Average distance a neutron travels before it is absorbed in cm\n",
+ "\n",
+ "#Output\n",
+ "print \"The thermal neutron flux = %0.2e cm**-2s**-1\"%F\n",
+ "print \"The average distance that a neutron travels before it is absorbed = %0.3f cm\"%L"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.8 Pg: 651"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 18,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ " The logarithmic energy decrement representing the neutron energy loss per elastic collision = 0.158 \n",
+ " The number of collisions necessary = 121 \n"
+ ]
+ }
+ ],
+ "source": [
+ "from math import log\n",
+ "#Input data\n",
+ "Eni=4.8##The energy of the newly born electron in MeV\n",
+ "Enf=0.025##The energy of the electron after slow down in eV\n",
+ "A=12##The mass number of the graphite (carbon)\n",
+ "\n",
+ "#Calculations\n",
+ "L=1-(((A-1)**2/(2*A))*log((A+1)/(A-1)))##The logarithmic energy decrement\n",
+ "n=(log(Eni*10**6/Enf))/L##The number of collisions required to slowdown the neutron \n",
+ "\n",
+ "#Output\n",
+ "print \" The logarithmic energy decrement representing the neutron energy loss per elastic collision = %3.3f \\n The number of collisions necessary = %3.0f \"%(L,n)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.9 Pg: 651"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 26,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ " (a) The rating of the reactor = 3.36 MW/tonne \n",
+ " (b)The rate of consumption of U-235 per day = 0.353 kg/day (or) 353 g/day \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "f=100##The reactor is fuelled of natural uranium in tonnes\n",
+ "A=238.05##The atomic mass of natural uranium \n",
+ "F=10**13##The average thermal neutron flux in neutrons/cm**2s\n",
+ "A1=235.04##The atomic mass of U-235\n",
+ "sf=579##The fission cross section of U-235 in barns\n",
+ "sc=101##The capture cross section of U-235 in barns\n",
+ "E=200##The energy released per fission in MeV\n",
+ "P=0.715##U-235 in natural uranium in percentage\n",
+ "N=2200##The average thermal neutron in m/s\n",
+ "\n",
+ "#Calculations\n",
+ "n=((f*1000)*6.023*10**26*(P/100))/A##The number of U-235 atoms in the reactor in atoms\n",
+ "R=(sf*10**-24)*F*n##The rate of fission in the reactor in fissions/s\n",
+ "T=R*E*1.602*10**-19##Thermal power of the reactor in MW\n",
+ "Rr=T/f##Rating the reactor MW/tonne\n",
+ "Rc=(((R*A1*60*60*24))/(6.023*10**26))##The rate of consumption of U-235 by fission in kg/day\n",
+ "Rcc=Rc*1000##The rate of consumption of U-235 by fission in g/day\n",
+ "\n",
+ "#Output\n",
+ "print \" (a) The rating of the reactor = %3.2f MW/tonne \\n (b)The rate of consumption of U-235 per day = %3.3f kg/day (or) %3.0f g/day \"%(Rr,Rc,Rcc)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.10 Pg: 652"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 27,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The specific energy release rate for a light water uranium reactor = 138.13 W/cm**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Input data\n",
+ "f=3.5##Mass fraction of U-235 in the fuel in percentage\n",
+ "G=180##Energy per fission in Mev\n",
+ "F=10**13##The neutron flux in neutrons/cm**2s\n",
+ "sf=577##Fission cross section of U-235 in barns\n",
+ "M=1.602*10**-13##One MeV in J\n",
+ "\n",
+ "#Calculations\n",
+ "N=2.372*(f/100)*10**22##The fuel density for a uranium oxide fuel in nuclei/cm**3\n",
+ "q=G*N*sf*10**-24*F##The rate of energy release in MeV/cm**3s\n",
+ "qg=q*M##The rate of energy release in W/cm**3\n",
+ "\n",
+ "#Output\n",
+ "print \"The specific energy release rate for a light water uranium reactor = %3.2f W/cm**3\"%(qg)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex: 9.11 Pg: 653"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 28,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "The reactor power level at the end of 1s is 3.196 MW\n"
+ ]
+ }
+ ],
+ "source": [
+ "from math import exp\n",
+ "#Input data\n",
+ "P=1##The operating power of a reactor in W\n",
+ "K=1.0015##The effective multiplication factor of Reactor becomes suppercritical \n",
+ "t=0.0001##The average neutron life in s\n",
+ "t1=1.0001##Neutron life time in s\n",
+ "\n",
+ "#Calculations \n",
+ "d=(K-1)/K##The reactivity \n",
+ "Z=(d*P)/t##The number of neutrons\n",
+ "n=exp(Z)/10**6##Neutron density * 10**6\n",
+ "\n",
+ "#Output\n",
+ "print \"The reactor power level at the end of 1s is %3.3f MW\"%(n)"
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.9"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}