Difference between revisions of "File:G17t1694x1281.png"

Power of relaxation heat power per kg of mixture of Am-241 and Pu-241 versus time for varios cases.

The mixture refers to initial concentation 85% of Am-241 and 15% of Pu-241

Description

C++ generator of curves

// file ado.cin should be loaded in order to compile the code below.

#include <iostream>
#include <cmath>
//#include <iomanip>
//using namespace std;
#define DB double
#include "ado.cin"

DB ln2 = log(2.0);
DB T1 = 14.33; // Pu-241 half-life (years) https://www.chemlin.org/isotope/plutonium-241
DB T2 = 432.6; // Am-241 half-life (years) https://www.chemlin.org/isotope/americium-241
DB L1 = ln2 / T1; // Year^-1
DB L2 = ln2 / T2; // Year^-1
DB MeV=1.60218e-13; // Joule
DB Year=365.242374*24*3600; //Second
DB amu=1.660539040e-27 ; //kg
DB DecayEnergy=5.63782*MeV; // https://www.chemlin.org/isotope/americium-241
DB EA=DecayEnergy*L2/Year / (241.05683*amu); // Watt/kg for Am-241
DB DecayEnergyPu=0.02078*MeV; // https://www.chemlin.org/isotope/Plutonium-241
DB EB=DecayEnergyPu*L1/Year / (241.05685*amu); // Watt/kg for Pu-241, https://www.chemlin.org/isotope/plutonium-241
DB A(DB t) {return EA*exp(-L2*t); }
DB B(DB t) {return EA* (L1/(L2-L1)) * ( exp(-L1*t) - exp(-L2*t) ); }
DB C(DB t) {return EB*exp(-L1*t); }

int main() { int m,n; DB d,t;
printf("EA=%9.5lf EB=%9.5lf (W/kg) , Initial energy generation by Am and by Pu\n",EA,EB);
//for(n=0;n<201; n+=10) {t=n; printf("%03.0f %8.4lf %8.4lf %8.4lf\n",t,A(t),B(t),A(t)-B(t));}
for(n=0;n<201;n+=10) {t=n; 
  printf("%03.0f %8.4lf %8.4lf %8.4lf %8.4lf\n", t,  A(t),  B(t),  B(t)+.3*C(t), .85 * A(t) + .15 * (B(t)+.4*C(t)) );
}

printf("A(0)=%4.2f\n",A(0.));
printf("A(100)=%4.2f\n",A(100.));
printf("(A(0)-A(100))/A(100)=%4.2f\n",(A(0.)-A(100.))/A(100.));

DB M32=.85*A(32)+0.15*(B(32)+.3*C(32));
DB M100=.85*A(100)+0.15*(B(100)+.3*C(100));
DB MM=(M32-M100)/M100;
printf("M32=%4.2f\n",M32);
printf("M100=%4.2f\n",M100);
printf("(M32-M100)/M100=%4.2f\n",MM);

FILE *o; o=fopen("g17.eps","w"); 
ado(o,1620,1240);
fprintf(o,"10 10 translate 10 10 scale 2 setlinecap 1 setlinejoin\n");
#define M(x,y) fprintf(o,"%3.0lf %5.3lf M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%3.0lf %5.3lf L\n",0.+x,0.+y);
M(0,124) L(0,0) L(162,0)
M(0,100) L(160,100)
fprintf(o,".5 W S\n");

for(n=10; n<121; n+=10){ M(0,n) L(160,n)}
for(m=10; m<161; m+=10){ M(m,0) L(m,120)}
fprintf(o,".3 W S\n");
M(0,125) L(100,125)
M(0,115) L(100,115)
M(0,105) L(100,105)
M(0, 95) L(100, 95)
fprintf(o,".12 W S\n");

for(n=0;n<164;n+=4){t=n; d=.3*C(t); if(n==0) M(t,d) else L(t,d)} fprintf(o,"0 0 0 RGB .2 W S\n"); // beta by Pu
for(n=0;n<164;n+=4){t=n; d=A(t); if(n==0) M(t,d) else L(t,d)} fprintf(o,"0 0 1 RGB .3 W S\n"); // alpha by Am
for(n=0;n<164;n+=2){t=n; d=B(t)+.3*C(t); if(n==0) M(t,d) else L(t,d)} fprintf(o,"1 0 0 RGB .4 W S\n"); // alpha by Am by Pu	
for(n=0;n<164;n+=4){t=n; d=.85*A(t)+0.15*(B(t)+.3*C(t)); if(n==0) M(t,d) else L(t,d)} fprintf(o,"0 .6 0 RGB .4 W S\n"); // mixture, no beta	
for(n=0;n<164;n+=1){t=n; d=B(t); if(n%2==0) M(t,d) else L(t,d)} fprintf(o,".5 0 0 RGB .2 W S\n"); // alpha by Am by Pu	

fprintf(o,"showpage\n");
fclose(o);

system("epstopdf g17.eps");
system( "open g17.pdf");
return 0;
}

g17.eps

the file will loaded upon request, 6264 bytes.

Output

The code above generates the picture, but also the text output that may help to interpret the result:

EA=114.57715 EB= 12.74890 (W/kg) , Initial energy generation by Am and by Pu
000 114.5772  -0.0000   3.8247  98.1555
010 112.7559  43.5626  45.9205 102.8485
020 110.9637  69.7263  71.1800 105.0688
030 109.1999  85.1748  86.0709 105.7753
040 107.4641  94.0280  94.5805 105.5592
050 105.7560  98.8261  99.1667 104.7846
060 104.0750 101.1347 101.3446 103.6759
070 102.4207 101.9188 102.0482 102.3713
080 100.7927 101.7732 101.8530 100.9557
090  99.1906 101.0645 101.1137  99.4815
100  97.6139 100.0184 100.0487  97.9807
110  96.0623  98.7741  98.7928  96.4728
120  94.5354  97.4170  97.4286  94.9700
130  93.0328  95.9999  96.0070  93.4793
140  91.5540  94.5549  94.5593  92.0050
150  90.0987  93.1018  93.1045  90.5497
160  88.6666  91.6527  91.6544  89.1149
170  87.2572  90.2149  90.2159  87.7011
180  85.8703  88.7926  88.7932  86.3087
190  84.5054  87.3884  87.3888  84.9379
200  83.1621  86.0038  86.0041  83.5884
A(0)=114.58
A(100)=97.61
(A(0)-A(100))/A(100)=0.17
M32=105.75
M100=97.98
(M32-M100)/M100=0.08

g17t.tex

\documentclass[12pt]{article}
\usepackage{geometry}
\topmargin -96pt
%\oddsidemargin -64pt
\oddsidemargin -44pt
\paperwidth 1700pt
\paperheight 1286pt
\textwidth 1600pt
\textheight 1300pt
\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}(1700,1254)
\put(50,30){\ing{g17}}
\put(-9,1226){\sx{5}{$P,$}}
\put(-5,1158){\sx{5}{$\rm \frac{W}{kg}$}}
\put(-32,1024){\sx{5}{$100$}}
%\put(-8,924){\sx{5}{$90$}}
\put(-8,824){\sx{5}{$80$}}
%\put(-8,724){\sx{5}{$70$}}
\put(-8,624){\sx{5}{$60$}}
%\put(-8,524){\sx{5}{$50$}}
\put(-8,424){\sx{5}{$40$}}
%\put(-8,324){\sx{5}{$30$}}
\put(-8,224){\sx{5}{$20$}}
%\put(-8,124){\sx{5}{$10$}}
\put(4, 24){\sx{5}{$0$}}
\put(46,-9){\sx{5}{$0$}}
%\put(128,-9){\sx{5}{$10$}}
\put(228,-9){\sx{5}{$20$}}
\put(428,-9){\sx{5}{$40$}}
\put(628,-9){\sx{5}{$60$}}
\put(828,-9){\sx{5}{$80$}}
\put(1020,-9){\sx{5}{$100$}}
\put(1220,-9){\sx{5}{$120$}}
\put(1420,-9){\sx{5}{$140$}}
\put(1584,-6){\sx{5}{$t,$ y}}
\end{picture}
\end{document}

Conversion

In order to improve the potability, compatibility with various versions of wiki, the picture is converted to png8 with the command below:

convert g17t.pdf png8:G17t1694x1281.png

Similar pictures

The idea is borrowed from the result by Ruvim Pechor, personage of the sci-fi Tartaria, see https://mizugadro.mydns.jp/t/index.php/File:Pu241heatPlo800x517.jpg

The draft of the description and similar pictures (perhaps with some mistakes) are loaded at https://mizugadro.mydns.jp/w/ChatGPT/TRY05/

References

Keywords

Am-241 Pu-241