File:Puconcentrationsplot.jpg
Relative concentrations $C$ of various isotopes of Plutonium, measured at the unload of fuel from various reactors, versus burn-out $B$, measured in teraJoules per kilogram of the fuel loaded.
Data are by Plutonium Topic Strategy, 2008.[1].
In the original table, the burn-out parameter $B$ is measured in MW d/t that seems to be MegaWatt days per tonne. This parameter is converted to the International system, taking into account, that day is 24*3600 seconds, and tonne is 1000kg.
References
- ↑ {https://web.archive.org/web/20110525233310/http://www.nda.gov.uk/documents/upload/Plutonium-Options-for-Comment-August-2008.pdf} NDA Plutonium Options. SMS Plutonium Topic Strategy, August 2008 - October 2008. Table 1: Examples of the types of variation in plutonium composition produced from different sources
Table of the input
RType 86400 238 239 240 241 242 100
Magnox 3000 0.1 80 16.9 2.7 0.3 82.7
Magnox 5000 -1 68.5 25.0 5.3 1.2 73.8
CANDU 7500 -1 66.6 26.6 5.3 1.5 71.9
AGR 18000 0.6 53.7 30.8 9.9 5.0 63.6
BWR 27500 2.6 59.8 23.7 10.6 3.3 70.4
BWR 30400 -1 56.8 23.8 14.3 5.1 71.1
PWR 33000 1.3 56.6 23.2 13.9 4.7 70.5
PWR 43000 2.0 52.5 24.1 14.7 6.2 67.2
PWR 53000 2.7 50.4 24.1 15.2 7.1 65.6
Reactor Type & Mean fuel burn-up (MW d/t) & Percentage of & u Isotopes at Discharge & Fissile Content
C++ generator of curves
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#define DB double
void ado(FILE *O, int X, int Y)
{ fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%');
fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y);
fprintf(O,"/M {moveto} bind def\n");
fprintf(O,"/L {lineto} bind def\n");
fprintf(O,"/S {stroke} bind def\n");
fprintf(O,"/s {show newpath} bind def\n");
fprintf(O,"/C {closepath} bind def\n");
fprintf(O,"/F {fill} bind def\n");
fprintf(O,"/o {.01 0 360 arc C S} bind def\n");
fprintf(O,"/O {.03 0 360 arc C S} bind def\n");
fprintf(O,"/times-Roman findfont 20 scalefont setfont\n");
fprintf(O,"/W {setlinewidth} bind def\n");
fprintf(O,"/RGB {setrgbcolor} bind def\n");}
int main(){int n,m;
char name[10][10]; for(n=0;n<10;n++) sprintf(name[n],"00000\n");
DB E[10],e[10];
DB p238[10];
DB p239[10];
DB p240[10];
DB p241[10];
DB p242[10];
DB p243[10]; // (this is hust sum)
DB MWt=1.e6; //Watt
//DB Year=365.24*24*60*60; // second
DB Day=24*60*60; //second
printf("Day=%8.4lf %8.4le \n",Day, 53*Day);
DB tonne=1000.; //kg
FILE *o; o=fopen("table.txt","r");
for(n=0;n<10;n++){ m=fscanf(o,"%s%lf%lf%lf%lf%lf%lf%lf",name[n], &E[n], &p238[n],&p239[n],&p240[n],&p241[n],&p242[n],&p243[n]);
e[n]=E[n]*(MWt*Day/tonne);
}
fclose(o);
n=0;
printf("%2d %6s %6.0lf %8.2le %4.0lf %4.0lf %4.0lf %4.0lf %4.0lf %4.0lf\n",n,name[n], E[n], e[n], p238[n],p239[n],p240[n],p241[n],p242[n],p243[n]);
for(n=1;n<10;n++){
printf("%2d %6s %6.0lf %8.2le %4.1lf %4.1lf %4.1lf %4.1lf %4.1lf %4.1lf\n",n,name[n], E[n], e[n],p238[n],p239[n],p240[n],p241[n],p242[n],p243[n]);
}
o=fopen("21.eps","w"); ado(o,504,204);
#define M(x,y) fprintf(o,"%6.4lf %6.4lf M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%6.4lf %6.4lf L\n",0.+x,0.+y);
#define O(x,y) fprintf(o,"%6.4lf %6.4lf O\n",0.+x,0.+y);
fprintf(o,"2 2 translate 100 100 scale 2 setlinecap 1 setlinejoin\n");
DB sy=2;
M(5,0) L(0,0) L(0,2) fprintf(o,".007 W S\n");
for(n=1;n<6;n++){ M(n,0)L(n,2) }
for(n=1;n<11;n++){ M(0,.2*n)L(5,.2*n)}
fprintf(o,".005 W S\n");
DB s=1.e12; // teraJoule
DB p=.02; // 2 percent
M(0,2) for(n=1;n<10;n++) L(e[n]/s,p*p239[n]) fprintf(o,".01 W S\n");
M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p240[n]) fprintf(o,"1 0 0 RGB .01 W S\n");
M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p241[n]) fprintf(o,"0 .7 0 RGB .01 W S\n");
M(0,0) for(n=1;n<10;n++) L(e[n]/s,p*p242[n]) fprintf(o,"0 0 .9 RGB .01 W S\n");
M(0,0) for(n=1;n<10;n++) {if(p238[n]>0) L(e[n]/s,p*p238[n])} fprintf(o,".8 0 .9 RGB .01 W S\n");
fprintf(o,"0 0 0 RGB .01 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p239[n])
fprintf(o,"1 0 0 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p240[n])
fprintf(o,"0 .6 0 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p241[n])
fprintf(o,"0 0 .9 RGB .02 W\n"); for(n=1;n<10;n++) O(e[n]/s,p*p242[n])
fprintf(o,".8 0 .8 RGB .02 W\n");for(n=1;n<10;n++){if(p238[n]>0){O(e[n]/s,p*p238[n])}}
fprintf(o,"showpage\n");
fclose(o);
system("epstopdf 21.eps");
system( "open 21.pdf");
}
Latex generator of labels
\documentclass[12pt]{article}
\usepackage{geometry}
\topmargin -96pt
\oddsidemargin -70pt
\paperwidth 518pt
\paperheight 216pt
\textwidth 1800pt
\textheight 1720pt
\usepackage{amsmath}
\usepackage{graphicx}
\usepackage{rotating}
\parindent 0pt
\newcommand{\sx}{\scalebox}
\newcommand \ing {\includegraphics}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\begin{document}
\begin{picture}(510,201)
\put(10,10){\ing{21}}
%\put(0,924){\sx{6}{Bu, \%}}
\put(-2,204){\sx{1.1}{$C,$}}
\put(1,188){\sx{1.1}{$\%$}}
\put(-1,168){\sx{1.}{$80$}}
\put(-1,128){\sx{1}{$60$}}
\put(-1,88){\sx{1}{$40$}}
\put(-1,48){\sx{1}{$20$}}
\put(2, 8){\sx{1}{$0$}}
\put(9, 1){\sx{1.1}{$0$}}
\put(109, 1){\sx{1.1}{$1$}}
\put(209, 1){\sx{1.1}{$2$}}
\put(309, 1){\sx{1.1}{$3$}}
\put(409, 1){\sx{1.1}{$4$}}
%\put(406,-6){\sx{1.3}{$4$}}
\put(454,1){\sx{1.1}{$B$, TJ/kg}}
%
\put(40,76){\sx{1.2}{\rot{90}\bf Magnox \ero}}
\put(60,76){\sx{1.2}{\rot{90}\bf Magnox \ero}}
\put(81,76){\sx{1.2}{\rot{90}\bf CANDU\ero}}
\put(172,78){\sx{1.2}{\rot{90}\bf AGR\ero}}
\put(254,72){\sx{1.2}{\rot{90}\bf BWR\ero}}
\put(280,72){\sx{1.2}{\rot{90}\bf BWR\ero}}
\put(303,72){\sx{1.2}{\rot{90}\bf PWR\ero}}
\put(391,68){\sx{1.2}{\rot{90}\bf PWR\ero}}
\put(478,68){\sx{1.2}{\rot{90}\bf PWR\ero}}
%
\put(403,111){\sx{1.3}{Pu-239}}
\put(403,56){\sx{1.3}{Pu-240}}
\put(403,39){\sx{1.3}{Pu-241}}
\put(403,22){\sx{1.3}{Pu-242}}
\put(403,10){\sx{1.3}{Pu-238}}
%
\end{picture}
\end{document}
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 06:13, 1 December 2018 | | 1,075 × 448 (125 KB) | Maintenance script (talk | contribs) | Importing image file |
- You cannot overwrite this file.
File usage
There are no pages that link to this file.
