Difference between revisions of "File:Puconcentrationsplot.jpg"

From TORI
Jump to navigation Jump to search
(Importing image file)
 
 
Line 1: Line 1:
  +
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.
−
Importing image file
  
  +
  +
Data are by Plutonium Topic Strategy, 2008.<ref>{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
  +
</ref>.
  +
  +
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==
  +
<references/>
  +
  +
==Table of the input==
  +
<poem>
  +
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
  +
</poem>
  +
  +
==[[C++]] generator of curves==
  +
<poem><nomathjax><nowiki>
  +
#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");
  +
}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] generator of labels==
  +
<poem><nomathjax><nowiki>
  +
\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}
  +
</nowiki></nomathjax></poem>
  +
  +
[[Category:C++]]
  +
[[Category:Explicit plot]]
  +
[[Category:Latex]]
  +
[[Category:Nuclear waste]]
  +
[[Category:Plutonium]]
  +
[[Category:Tartaria]]

Latest revision as of 08:47, 1 December 2018

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

  1. ↑ {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/TimeThumbnailDimensionsUserComment
current06:13, 1 December 2018Thumbnail for version as of 06:13, 1 December 20181,075 × 448 (125 KB)Maintenance script (talk | contribs)Importing image file

You cannot overwrite this file.

There are no pages that use this file.

Metadata

This file contains additional information, probably added from the digital camera or scanner used to create or digitize it.

If the file has been modified from its original state, some details may not fully reflect the modified file.

Retrieved from "https://mizugadro.mydns.jp/t/index.php?title=File:Puconcentrationsplot.jpg&oldid=20943"