tesis_uam/Codigos para el manejo de un reactor ... - cedip

UNIVERSIDAD AUTONOMA METROPOLITANAUNIDAD IZTAPALAPAc. 4. I.SEMINARIO DE PROYECTOS (I Y TI)JRPLIC*CTON DE UN PAQUETE DECODIGOS CON LOS FUNDAMENTOSDE INGENIERIA NUCLEAR!iSALVADOR &SASERRANOMEXICO D + F+ MARZO 1984,, I

I N D I C ECONTENIDOPAGINA tIt1 OBJETIVOS DEL SEMINARIO t t t t t t t t t t t t t t t t t 2CAPITULO 3t CODIGO THERMAt t t t t t t t t t t t t t t t t t t t t t t t t t t313.1 PROPOSITO DEL PROGRAMA t t t t t t t t t t t t t t t t t t323t2 TEORIA EN QUE SE BASAt t t t t t t t t t t t t t t t t t t323t3 DIAGRAMA DE FLUJOt t t t t t t t * t t t t t t t t t t t t t t t383t4 USO DEL PROGRAMA (ENTRADAS Y SALIDAS)t t t t393t5 SIMULACION t t t t ~ ~ t t t t t t t t t ~ t t t ~ t t t t t t t t40t t3t5t1 CAHACTERISTICAS DE LOS NUCLEOS A SIMULAR413 , s RESULATDOS'~t t t t t t t t t * t t t t t t t t t t t t t t t t t t41

6.5 SIMULACION t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t 114CAPITULO 7t CODICO DECAY t t t t t t t t * t t * * t * t 4 4 , . * 4 * * * * * * * 1287tl PHOYOSITO DEL PROGRAMA 4 4 * 4 4 t t 4 4 t t t 4 t e t t t 1297t2 TEORIA EN QUE SE BASA t t t t t t t t t t t t e t t t t t e 1297t3 DIAGRAMA DE FLUJO t t t t t t t t t t * t t t t t t t t t t t t 1337t4 USO DEL PROGRAMA (ENTRADAS Y SALIDAS) t t t 1377t5 SIMULACION t t t t t t t t t t t t t t t t t t * t t t t t t t t t t t 1407tSt1 CADENAS UTILIZADAS PARA LA SIHULACION t t 1407t5t2 RESULTADOS t t t t t t t t t t t t t t t t t t t t * t t t t t t t 1417t5t3 GRAFICAS t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t 1417t6 CONCLUSIONES t t t t t t t t t t t t t t t t t t t t t t t * t t t t 157CAPITULO 8t CODIGO PLOT t t t t t t t t t t t t t t t t t t * t t t t t t t t t t t 1388tl PROPOSITO DEL PROGRAHA t t t t t t t t t t t t t t t t t t 1598.2 CAPACIDADES DE PLOT t t t t t t t t t t t t t t t t t t t t t 1598t3 USO DEL PROGRAMA (ENTRADAS Y SALIDAS) t t t 1598t4 DIAGRAMA DE FLUJO t t t t t t t t t t t t t t t t t t t t t t t 1628t5 PLOT COMO SUBRUTINA t t t t t t t t t t t t t t t t t . t t t 1628t6 SIMULACION t t t t t t t t t t t t t t t t t t t t t t t t t t t t t t 1648t7 CONCLUSIONES t t * t t t t t t t t t t t t t t t t t t t 4 t t t t t 174A P E N D I C E S

APENOICE A t * t t t t t t * t t t t t t t t t t t * t t * t t t t t * * * 175APENOICE B * * t t t t t * t t t t t * t t t t t t t t t t t t t * t t t * 180APENDICE C t t t t t t t t ~ t t t t t , t * t t t t t * t t t t t t t t t 189AFENDICE D t t * t t * t t t * t t t t t t t t t t * t t * t t t t t t t * 193APENOICE E t t * * t t t t t t t t t t t t t t t t t t t t t t t t t t * * 296AFENDICE F t t t t t t * t t t t t t t ~ * t t * t t * t t t t t * t t t t 204APENDICE G t t t t t t t * * t * t t t t t t t t t t t * t t t t t t t t t 211

CAF:'ITlll..O I +INTRODUCCTQNEl ... OE:JETIUO DE ESTE SEMXNARIO ES HACER UN MANllAl..., PARA EL.USO DE 1...08 CODTGOS ARCJI PRESENTRDOC 9 LOS CUf4LE.C SERVIRAN PARAREFCIRZAF;: LAS MATEi:RIAS IMPARTIDA!? POR EL AREh DE RECIJRCO$$ENERGETICOS BE: LA UNIVERSInAD AUTCINOMA METROPOLITANA ...''-,EN LA RAMA NUCLEAR,RI

EN'T'E:RCI:!; : 2 ,:i CI 2 + ,:3 +

2.5

I ’FS’SUHA 2T6;‘Y[PBBA“EUBATROH/2. .’ 1iFX GCiA

I $ +!... h E: i:: t.l R i:: :I: i:f N D lif: L. R I:' R i:i1:) U i::f; :I I:i f4 1:)E: i: Ai., O i:i Y I:) E: F' I. IJ ,.I O 13 1 : i:A i. I:) I?to

2.8

J3/:1C1,I!B'YPBBATuBA1CENTRODZCANAL\ REFRIGERANTE 1IDISTANCIA DENTRO32I.1 CANAL.FIGURA 2e3DESTRIBUCION DE TEJIPERATUiM D2 UNY,iiTALICO Y ¡?~-Ni3!T!fiICO BAJO ZUJOREPRI GZR,Q!TBT imBULZiT0.I3

Nu= h be2. i3?r=QPK2.\4

IUw Y V.P CON LOS UOl,UMENi~:C DEL AGUA Y IJEL COMBUS'Y%CliSLETIVANENTE: EN EL. ENREJADU ESTAS CON I...ACi CONSI'AN'I'EC' QUEL.JZAN EN ESTE COD'JI.E[J+

1I INICIA I1 .'\'

E N '1' F( \3 A s.."........"_................................ ........................................

p.cl (4 1, :f 1. 5....................................................................................

dlnrrlLnrlNrlu)Orlui03OOLnrlO0(oILlrlOOrl rl-tmc-u)InrlLFrlimi3wHRou-

__I_ ___I_____ I ^.---_____....'.O... .,....O....,.... O....,....O..".,....O... '.""b . "~' v . . ' . 9.. . '... . t!LO**CANAL REFRIGERANTEA =.RADIO DE LA PASTILLA DEL COMBUSTIBLE (PULG)R = ESPESOR DEL ENCAMISADO (PULGI?V = VELOCIDAO DEL REFRIGERANTE (FT/SEG)KC = CONDUCTIVIDAD TERMICA VAINA [BTU/LB-F)KF = CONOUCTIVIDAD TERMICA 8-F 1ANOTA RHOpMU,CP,K,KC,KF :0 .TUS DATOS SON :5 = ,02250c = .49600164,39999TBO S6OoOOORHO = 44,60000MU Z ,22200KF = 1~50000*** SALIDAS ***,B50,206 ITCtMAX) = 2898,682 .129Q (MAX) = ,532E+08

2 CFf) TMAX (COMB10 56080 516080 56080.I 560.4 565-8 820-2-5.8 561.8 572.4 1077.9-5.3 563-9 1329.6-488 €57281-4-3 802.2-3,8 01689-3.3 2213-4-2.8 586 3 623.3 2389-0592-7 632.7 2541.412668882769.312841 872885-0.2 629-5 676-1 2898-73.2 671.6 706.6 22375.16e2 693-5 701.0 1056.4798.4S~HtJL~CiOIu

-6.8500 ! C-3.8500 !GRA018 0.9 1.----+-----+-----+II V C !t v C i! V zt V !1 V C IHV C 1CtCCiiC 1C !-.8500 ! RV C!CCcC!1,6500 ! RVC !I2.1500 I V C4.6500 ! RV5.1500 ! RV5,6500 ! RVCCCIIII0.0 0.1 0.13 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.FACTOR DE NQRMALIZACION =.290E+04C=TEMP DE COMB sIhuuc3oiu 1,FACTOR DE NORMALIZACION = ,290E+OYS = ESPACIAMIENTO ENTRE BARRAS PITCH (PULG)H = ALTURA DEL NUCLEO ACTIVO (PULG)RHO = DENSIDAD DEL REFRIGERANTE (LB/FT**3)1

0.8 0.9 la-5.8375 ! C IC 1I11-2.8375 ! V C !HV C 1R c 1R c 1RR----+-----+----o+,1625 ! RV CC!c 1!II3.1625 ! V C 1ICICILI1IiI+ ~ ~ ~ ~ ~ ~ o ~ ~ ~ ~ + ~ ~ ~ ~ ~ + ~ - o ~ ~ + ~ .1.V= TEMP DE VAINA!IQUIERES ANOTAR OTROS DATOS (S/N)?A = RADIO DE LA PASTILLA DEL COMBUSTIBLE (PULG)!3 = ESPESOR DEL ENCAMISADO [PULG)V = VELOCIDAD DEL REFRIGERANTE (FT/SEG)KC = CONDUCTIVIOAD TERMICA VAINA (BTU/LB-F)

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0FACTOR DE NORMALIZACION =.297E+04V= TEMP DE VAINA ~ILWIKACibn) 3FACTOR DE NORMALIZACION =.297E+04- _ - -*** ENTRADAS *t*A = RADIO DE LA PASTILLA DEL COMBUSTIBLE (PULG)B = ESPESOR DEL ENCAMISADO (PULGIV = VELOCIDAD DEL REFRIGERANTE (FT/SEG)(F IKC = CONDllCTIVIDAD TERMICA VAINA (BTU/LB-F)125

~ ~_ ...... ^. .. ...... :........ .... .................. .. .... l...I.. _.._.,. ............I ......... 1- ..... . . .. . . . . .

-6,7257 ! CA0.6 O,? 0.8 0.9 I o-.*-*+*----+o.-.-+-3.7257 ! R VCIR V CR C 1R 1-07257 ! R vCCCII11III8IIIc 1c1CC2.2743 ! R V C4,7743 ! R V C I5,2743 ! R V C IRV C 111f1R= TEMP DF REFRFACTOR DE NORMALIZACION =.715E+04RRR ENTRADAS ***01S = ESPACIAMIENTO ENTRE BARRAS PITCH (PULGIH = ALTURA DEL NUCLEO ACTIVO (PULGIN = NUMERO DE ELEMEhTOS COMBUSTIBLEMIJ = VISCOSIDAD DEL REFRIGERANTE (LBIHR-FT)n A I

' V " "___O....r....O,...r....O....r....0....r....O....r....0.,..r....O..,.r....~- -*** ENTRADAS ***A = RADIO DE LA PASTILLA DEL COMBUSTIBLE (PULG)P = POTENCIA DEL REACTOR (MWT)K z CONDUCTIVIDAO TERMICA REF, (BTUILB-FIANOTA PIVpTBO :ANOTA RHO,MU,CP,K,KC,KF :A = ,24600

Af 0.8 009 I m-3,8292 ! CIV CRVLR V 1R 1R C tR V C !R C!CCIc 1C !2.1708 ! R V C !!1IOm0 001 0.2 Oe3 0.4 0.5 Om6 0.7 Om8 O09 1.I*+I?+TEMP O€ RfFRFICPDR DE NORHICIZACION fVI T€MP DE VAINA.1?6E+04FACTOR DE NORMALIZACION = . 176€+04C=TEMP DE COMBFACTOR DE NORMALIZACION = .176€+04SIWLAClOEIS = ESPACIAMIENTO ENTRE BARRAS PITCH (PULGIRHO = DENSIDAD DEL REFRIGERANTE [tB/FT**3)KF = CONDUCTIVIDAD TERMICA COMB. (BTU/LB-F)

31070105

IiII32

’rI1\IIIL-/--BarraVombusti bleuM G W 3.1 DISTRIBUCIGN DEL FLTJJO EN u1J NUCLEOCILINDRIC O.33

Formaciondel nucleode' ebullicionPeliculaparcial deebullicionConveccion IFIGURA 3.2Ií/IT1Tfl. 'FIGURA 3.335

III .;1II'SFIGU RA 304 SUBEXFBXAXI'?IErJTO O BAJA CALIDAD DE-DMB.

u:' N .y r< i'j 6, f3.........."".. ........"..If: N 'r' R A I:) 61 !j 13 1:: I. A c U N F' 1; I:; u 13 A [: :r:w N ci I:.::................................... "...._........^.I.I- N u (= I_. E 13............ .................................. .38

,39

8 A I... :I I:> A s....................................................- .I.........

7-I

I* .42

ANOTA LA VELOCIDAD viASICA DETEMP. DF ENTRADA DEL NUCLEU; PRESIOlJ DE EIJTHAijA:A FUENTE (BTü/HH-FT**3COMB Cl.0.INCH ) INCH 1It,

. v . . . *0 ......... 0 ....,..-.. 0 . . . . , . . . . o ....,.." 0 " " v " " o . * ".'... o . .I . , . . . . o.A LA EiJTf?ADA UEL NiICLEU..............................................................1 -6.507 575.6 - 57'5.2 572.4 568.3 562.7561 .O 560.1 652.5 6655.1 68.42 .o04 -4.453 669.6 666.8 647.0 618.8 579.276 6 QQ e9ejRilas 00_ I rY8bo&*$'586.9 576.3 652.2 6733.7 69.76 -9u13 1.712 776.7 771.4 734.6 682.0 b08.4 1588 . 3 578.2 652-2 4742-4 69-93 - O 016 3.767 719.9 716.2 690.9 bS4.6590.1 563.0 652.1 6763-7 70-376U4.2I19 5.822 633.3 632.1 623.3 610.8 593.3588.2 585.6 652.0 6775-1 70-62 - o o

.r,...O....r....O....,....0.........0.........o...... . ..o.. .‘”.O. .......O-6.5075 ! F2 1-5.8225 ! F32l I-5,1375 ! I-2,397s ! F 3 21c !-1.7125 ! F 3 2 l !-1,0275 ! F 3 2 l!FACTOR DE NORt~ALIZACION = .789E+03F= TEMP UN RADIO2234.9 586.0

,....O....,....O...,,..~..O....,....O..........O....~..,.o....*....o....,....o.A LA ENTRABA OEL kUCLE015.57 307922.44 .O147..............................................................?da >Z4&0 79 ''I!W)a. * Tf*'Tk P R€!Lk€ION CffFt VAI; P.1 -5.704 586.6 586.2 583.6 579.8 s74.5572.9 572.1 652.4 6357.6 53.79 .o04 -3.903 673.9 671.3 653.0 626.710 -.300 781.7 776.4 739.6 656.9593.1 582.4 6402.3 54.5013 1 .5Ol 771 -7 766.9 732.8 684.3s97.4 587.5 652.1 6422.7 54.87 00016 3.302 718.8 715.4 692.0 658.5 611.6 1598.5 591.6 652.0 6438.3 5S.18 . o019 5,104 638.2 637 . 1 628.9 617.3 601.1596.3 593.6 652.0 6446.7 55.35 .u0

-5,1035 ! ' F 31 I-4.5031 ! F 321 I-1.5010 ! F 3 2 l !w.9006 ! F 3 2 l!2.1015 ! F 3 2 l !FACTOR DE NOHr4ALXZACION = .784€+031= TEM A 1/4 DEL f?AI)IOFACTOR DE hORMALfZACIUtd =F= TEMP iJN RADIu.f64€+03I2235.4 594.2TEPPERATURA DE SALIíjA* < .64

' V " "O....r....O....,....O.,..,....~....~....0....,....0....,....0....,....0ALA ENTRADA bEL IUIJCLEO1S.84 309829.13 .O1471 -5.657572.9z-5,064 6JPI,?3 14 'w4576.7586.1 585.7 583.2 579.5 574.4572.1 652.4 b415.4 55.39 .o0614.6 61.3.5 594.8 579.4i.%O572.6 652.3 6416.9 55.44 . o o596.7 ' I590.2 652.0 6491.8 56.71

' V " "o. .. .,.... O....r....O....r....O....r....~....~....~....,....~,...,....~F3t-5,6569 !-5.0614 ! F32C IFACTOR DE NOHMALIZACIOrJ = .775€+031= TEM A 114 DEL RADIOFACTOR DE hOHMALIZACIDN =F= TEMP UN RADIOFACTORPRES I orJi.7fCE+03I2234-6 592.6

’ V “ ” O....r....O....r....0....r....~....r....O...r.......~....,....~....,....o15.58 296732.94 .O148..............................................................1 -5.938 579.1 578.7 576.1 572.5 5b7 .456’5.7 564.9 652.4 b333.9 58.27 .o04 -4.063 664.6 662.1 644.2 618.7 5653.07 -2. 1 BB 734.7 730.4 700.3 657 . 3 597.25as.5 575.4 652.2 6382.0 59.02 .o0 I13 1.563 760.8 756.0 723.0 675.7 609.6 ~589.9 580.6 652.1 6404.3 S9.41 .o0I6 3.438 709.2 706.0 683.2 6’30.7 605.1591 -3 584-8 652-0 6421 -6 59-74 - ~- _o 0602,2\ . 594b.6*QB ’590.4 586.6 651.9 6420.8 59.88 .o019 5.313 630.6 429.4 621 .S 610.2 594.4589-5 587 - 1 651.9 6430-8 59-93 o o

-5.3125 !-4.6875 !-4.062s I-3.4375 I-2.1875 !F2IF 32IF 3?1tF 3 216 , 1' f 3 i?&C t*-2.BA2S L ,f 3 2 3c 'iF 3 2 1 !t 1- -FACTOR DE NORMALIZACION =1= TEM A 114 DEL HA010.772E+03, * -I3= TEM A 314 DEL RADIO

Í I -_I_' V " "~....r....~....V....0....r....O....~....O....V....~....V....~,...V....~A L A €PrTRADn DEL NUCLEOq Z.4AC................................................................. *70% k f 3&4 7Ff'VAf. *Y RW. :f@&,A e REkw8m e##1 -6 . 389 486.4 486 . 3 485.5 484.4 482.91'496.2 486.4 581 -4 552.6 -43 .o04 -4.372 530 . 4 529.6 524 . 5 517.1 506.8FORMACIUN 1'JETA DE VACIOI

'T""~....T....O....T....~....T....O....T....O....~....~....T....~....T....~ALA EVTHADA UEL NUCLEO3.59 136451 -63 ,0173..............................................................1 -1 .6?1 bl8.S 616.2 599.9 575.6 544.1537.2 520.1 605.1 1939.23 -1,280 1190-6 1086-3 986-5 843-9 644.2EBULL IC IOiJ SUBEivFR I ADA4.IE B ULL I C I ON S iJFJ EMF R I A U A-*EBULL IC ION SUBENFRIADA9 0.256 2411-5 2357.7 1982.4 1445-9 694-910 0.085 2474.1 2410.4 2029.9 1474.5 697.0613.4 540.5 605.1 7982.5 .59 044_.I -- ---e-.-L- n ... A . . . .n

1I- I613.3 573.4 605.0 10307.9 . 6 0EBULLICIOb S1JHENFHIADA695.0 lri15R.2 . 34EBULL C I ON SUBENFR I ADAEBULLICIOh SUBENFRIAOA? ., . .,. I '.n. . .n

‘ V ” ” O....r....O....r....0....r....0....V....~.........~....r....~....V....~-1.4503 ! F 3 21-1 ,2797 ! F 3 ? lf- - ---‘1,1090 ! f 1-09384 1 P t-.76?8 L F i !-05972 ! F 3 2 1c I-. Y266 ! F 3 1c . I-.I2559 ! F,5972 ! F 3 ? ir I,7678 !FACTOR DE NOR+lALIZACIOPd =.247E+041= TEM A 1/4 DEL RADIO -ON q .Este I3= TEM A 3/4 DEL RADIOFACTOR DE NORI1ALItACTDi~ = .247€+04F= TEMP IJhJ RADIO1582.0 603.1”.A I

A LA ENTRADA OEL kiJCLEü189.11 3514035.00 .O0901 -3.638 533.0 532.1 525.1 515.3 501.5496.4 496.0 627.5 49981.9 3604.47 .o04 -2,489 764.5 757.6 706.7 634.3 533.116 2.106 837.0 770.8 676.9 545.4502.7 499 . 7 599.0 50121.0 3614.3819 3.255 613.8 610.9 589.5 559.1 51b.6iHUl&m7r56- n .....

.'T""O....r....O....r....O....r....O.........O...,.....o....r....o....T....o-3.2548 ! F 3 21 I-2.8719 ! F 3 2 1 1*i?04890 tF . t1-2.t060 tíF? 1=lD'T231 1 v 1'-1,3402 ! F 3 2 1 1-09573 I F 3 a 1c !00573 ! F 3 2 1c !1.3402 ! F 3 I? 1 I1.7231 ! F 3 2 1c 1a.lO6O t2.4890 t2r87t9 t3,2548 ! F 3 213.6377 ! 5¿1 t1tua9 $.C= TEM E:J EL CENTRUFACTOR DE NORMALIZACION = . 108E+041= TEM A 1/4 DEL HAiIIOFACTOR DE IIORMALIZACION = . 108E+04F= TEMP IJni RADIO1048.0 500.3. . I. .. . . . .

' V " "0....r..,.o....,....0....r....O....r....O....r....O....,....O....r....O..........................1 -5.700676.5 673.3 651.9 621.3 578.44 -3.9001752.5 1719.8 1497.2 11711.0 731.1 ',

:I , I., c: I:] N C: I, 1.J C 2: I5 N I: t;\5%

CAPITULO 4r CODIGO ROOTS4+1 PROPOSITO DEL PROGRAMAROOTS ES UN CODIGO ESCRITO EN FORTRAN, PARA EL USOINTERACTIVO DE COMPUTADORA, PARA EXAMINAR LAS SOLUCIONES DE LAECUACION DE REACTIVIDAD+4+2 TEORIA EN QUE SE BASASABEMOS QUE LA OPERACTON DE LOS REACTORES NUCLEARES NOSIEMPRE ESTA EN ESTADO CRITICO, POR EJEMPLO, PARA UN REACTORt CIUE SE ARRANCA O AL ELEVAR SU NIVEL DE POTENCIA, ESSUPERCRITTCO, DE OTRO MODO, AL PARAR EL REACTOR (SHUT DOHN) OAL REDUCIR SU POTENCIA, ESTE ES SUBCRITICO+ EL ESTUDIO DELCOMPORTAMINETO DE LA POBLACION DE NEUTRONES EN UN REACTOR NOCRITICO SE LE DENOMINA 'CINETICA DEL REACTOR'+FEL GRADO DE CRITICIDAD DE UN REACTOR GENERALMENTE SE REGULACON EL USO DE BARRAS DE CONTRO O DE ELEMENTOS QUIMICOSI LAS 1BARRAS DE CONTROL SON PIEZAS O ENSAMBLES HECHAS DE MATERIALESABSORBENTES DE NEUTRONES, CUYO MOVIMIENTO TIENE EFECTO EN EL IFACTOR DE MULTIPLICACION DEL SISTEMA+ ASI, SI LAS BARRAS DE ICONTROL SE RETIRAN DE UN REACTOR CRITICO, EL REACTOR TIENDE A ISER UN REACTOR SUPERCRITICO, PERO SI LAS BARRAS SE INCERTAN, IEL SISTEMA CAE A SUBCRITICOI CON ELEMENTOS QUIMICOS, EL ICONTROL SE REALIZA VARIANDO LA CONCENTRACION DE UN ABSOBEDORQUIMICO DE NEUTRONES, GENERALMENTE ACID0 BORIC0 (H3B03) ENREACTORES MODERADOS O ENFRIADOS CON AGUA+LOS FACTORES MAS IMPORTANTES QUE AFECTAN LA CRITICIDAD ESLA TEMPERATURA DEL REACTOR+ UNO DE LOS PARAMETROS QUE SERELACIONAN CON EL VALOR K ES LA DEPENDENCIA DE LA TEMPERATURA,Y CAMBIOS DE T PROVOCAN CAMBIOS EN K,SABEMOS QUE MUCHOS DE LOS NEUTRONES EMITIDOS EN UNA FISIONAPARECEN VIRTUALMENTE EN EL INSTANTE EN QUE OCURRE LA FISION;ESTOS SON LOS NEUTRONES INMEDIATOS+ UNA FRACCION MUY PEQUENADE LOS NEUTRONES DE FISION APARECEN UN TIEMPO DESPUES DEOCURRIDO EL EVENTO DE LA FISION; ESTOS SON LOS NEUTRONESRETARDADOS. EL COMPORTAMIENTO DE UN REACTOR DEPENDE DE LASPROPIEDADES DE EBQS DOS TIPOS DE NEUTRONES+EN ESTE CASO CONSIDERAREMOS UN REACTOR CON NEUTRONESRETARDADOS+EN CALCULOS MAS PRECISOS DE LA CINETICA DEL REACTOR, SENECESITA CONSIDERAR EN DETALLE LA PRODUCCION Y LADESINTEGRACION DE CADA UNO DE LOS SEIS GRUPOS DE NEUTRONESRETARDADOS PRECURSORES+ PARA SXMPLIFICAR LA DISCUSION, POR ELMOMENTO SUPONDREMOS SOLO UN GRUPO DE NEUTRONES RETARDADOS QUEIII,I

070105FIGUiiA 40i GUPICA DE LA ECUACION DE Í3EACTIVIDAD PARAUN GRUPO DE NEUTRONES RETARDADOS.61

.APARECEN EN LA DESINTEGRACION DE UN PRECURSOR HIPOTETICOI PARAESTE CALCULO SE CONSIDERA UN REACTOR TERMICO HOMOGENEO QUEPUEDE SER O NO CRITICO, OBTENIENDO:4. IDONDE (\CP -l)/ \Le SE CONOCE COMO REACTIVIDAD, EN ESTE CASOPARA UN REACTOR INFINITO, SE DENOTA POR EL SIMBOLO p, AS1PERO PARA UN REACTOR FINITO SE DEFINE COMOENTONCESQUE SE CONOCE COMO ECUACION DE REACTIVIDAD, PARA UN GRUPO DENEUTRONES RETARDADOS+CONSIDERANDO EL SIGNIFICADO DE LA ECUACION DE REACTIVIDAD,SE DEBE NOTAR QUE SI UN REACTOR ES SUPERCRITICO, AS1 K 3 1, Y?ES POSITIVO+ EN ESTE CASO, EL REACTOR TIENE "REACTI'JIDADPOSITIVA"+ DE OTRO MODO, CUANDO EL REACTOR ES SUBCRITICO, y\

DONDE Al. Y A2 SON CONSTANTES+SE OBSERUA EN LA FIGURA(4+1), QUE CUANDO .f ES POSITIVA,wl, ES POSITIVA Y ~ 2 ES , NEGATIVA+ AS1 COMO EL TIEMPO EN ELSEGUNDO TERMINO DE LA EC+(4+,5) DECAE RAPIDAMENTE Y EL FLUJOCRECE COMO EXP(w1xt) + DE OTRA MANERA, CUANDO ES NEGATIVO,PERO wl Y ~2 SON NEGATIVOS, PERO WZ ES MAS NEGATIVO QUE wI.4AS1 EL SEGUNDO TERMINO DE LA EC+(4+5) DE NUEVO DECAE MASRAPIDAMENTE QUE EL PRIMERO Y EL FLUJO HORA DISMINUYE COMOEXP(wlxt)+ EN CADA CASO, p POSITIVA O ! NEGATIVA, EL FLUJOSE APROXIMA A EXP(wiXt), ES DECIR,EL INVERSO DE Wl, SE CONOCE COMO PERIODO DEL REACTOR OALGUNAS VECES PERIODO ESTABLE, Y SE DENOTA POR "T", AS1 QUET=l/wiEL COMPORTAMIENTO DEL FLUJO PUEDE SER ESCRITO COMO,AHORA AL INCLUIR LOS SEIS GRUPOS DE NEUTRONES RETARDADOS,EN LUGAR DE DOS, LA ECUACION RESULTANTE DE REACTIVIDAD, DEBETENER LA MISMA FORMA QUE LA EC+(4+5) PERO NECESARIAMENTE ESMAS COMPLICADA+ ESTA ECUACION GENERAL DE REACTIVIDAD ES,6DONDE Y %; SE REFIEREN AL iESJMO GRUPO DE NEUTRONESRETARDADOS, OBSERVANDO EL TERMINO DEL LADO DERECHO DE LAEC+(4+6), FACILMENTE SE MUESTRA QUE HAY SIETE RAICES DE LAECUACION PARA CADA VALOR POSITIVO O NEGATIVO DEL FLUJO, DENUEVO ESTA DADO POR UNA SUMA DE EXPONENCIALES, PERO HAY AHORASIETE EN LUGAR DE DOC+ SIN EMBARGO, SIMILAR AL CASO DE UNGRUFO DE NEUTRONES RETARDADOS, AL INCREMENTARSE EL TIEMPO #SE APROXIMA MAS A EXP(Wlxt), DONDE Wl, ES LA PRIMERA RAIZ DELA ECUACION, SI ES POSITIVA O NEGATIVA+ EL PERIODO DELREACTOR ES DE NUEVO EL INVERSO DE wi+AHORA PARA LA CINETICA DE REACTORES RAPIDOS, SABEMOS QUE ELTIEMPO DE VIDA DE LOS NEUTRONES INMEDIATOS EN ESTE TIPO DEREACTOR ES MUY CORTO, DEL ORDEN DE 10xr-7 SEGUNDOS+ POR LOTANTO EXCEPTO PARA INSERCIONES DE REACTIVIDAD MUY GRANDE QUERESULTA EN GRANDES VALORES DE w l (CORTOS PERIODOS DELREACTOR), LA CANTIDAD w++ EN LA EC+(4+6) PUEDE DESPRECIARSE,AS1 LA ECUACION DE REATJUIDAD SE APROXIMA A63

443 DIAGRAMA DE FLUJOPARA ESTE CODIGO EL DIAGRAMA DE FLUJO UNICAMENTRE MUESTRALAS OPERACIONES PRINCIPALES, PERO PARA MAYORES DETALLESCONSULTAR EL LISTADO, QUE ESTA EN EL APENDICE C+4+4 USO DEL PROGRAMA (ENTRADAS Y SALIDAS)EL CODIGO SE PUEDE USAR DE DOS MANERAS:1) PARA CALCULAR LAS RAICES DE LA ECUACION DE REACTIVIDAD+ii) PARA CALCULAR SOLAMENTE LA REACTIVIDAD+ADEMAS SE PUEDE USAR PARA UNO O SEIS GRUPOS DE NEUTRONESRETARDADOS+E N T R A D A S+====E===========EL USUARIO DEBE ESPECIFICAR LOS MODOS DE CALCULO+- EL NUMERO DE GRUPOS DE NEUTRONES RETARDADOS+- LAS CONSTANTES DE DESINTEGRACION ( 1 DE CADA GRUPO EN(l/SEG+), SE DEBEN ANOTAR EMPEZANDO POR LA MAS PEQUENA+- LA FRACCION DE NEUTRONES RETARDADOS ( 1 PARA CADAGRUPO +-. EL TIEMPO DE GENERACION DE LOS NEUTRONES (SEG+)EL USUARIO DEBE ANOTAR LOS SIGUIENTES DATOS PARA EL MODOCORRESPONDIENTE+PARA CALCULAR LAS RAICES (MODO i), DEBERA ANOTAR CERO, PARAEL CALCULO DE REACTIVIDAD (MODO ii) SE DEBERA ANOTAR i+LA SIGUIENTE PREGUNTA QUE HACE EL CODIGO ES "TUS DATOSESTAN CORRECTOS (S/N)?", EN CASO DE QUE ESTEN CORRECTOS LOSDATOS, EL USUARIO ANOTARA "N", Y EL CODIGO PREGUNTARA DE NUEVOTODAS LOS DATOS, PERO LAS LINEAS QUE ESTEN CORRECTAS BASTARACON PRESIONAR "RETURN", Y LOS DATOS ANTES ESCRITOS NO SERANALTERADOS, ENTONCES SOLO SE MODIFSCARAN LOS DATOS QUE NO ESTENERRONEOS, UNA VEZ QUE LOS DATOS ESTAN CORRECTOS "Si', EL CODIGOPROCEDE SEGUN EL MODO DE OPERACION INDICADO+MODO 1)

- EL USUARIO DEBERA ANOTAR LA REACTIVIDAD, DE LA QUE SEQUIERAN SABER SUS RAICESiMODO ii)- EL USUARIO ANOTARA EL PERIODO DEL REACTOR EN SEGUNDOS,QUE DEBE SER MAYOR QUE -1/ DEL PRIMER GRUPO,EN CADA MODO, EL CODIGO PREGUNTA PARA NUEVAS ENTRADAS (YASEA REACTIVIDAD O PERIODO.) DESPUES DE CADA CALCULO, LOSCALCULOS PUEDEN TERMINARSE ANOTANDO LOS SIGUIENTES DATOSMODO i)UNA REACTIVIDAD MAYOR O IGUAL A 1,MODO ii)-. PARA UN PERIODO DE CERO+MODO i)LAS SALIDAS PARA CADA MODO SONSE IMPRIMEN LAS RAICES DE LA ECUACION DE RACTIVIDAD PARA LAREACTSVIDAD QUE SE ANOTO (EL NUMERO VARIA SEGUN LOS GRUPOS DENEUTRONES RETARDADOS), Y EL PERIODO DEL REACTOR PARA ESAREACTIVIDAD,MODO ii)SE IMPRIME LA REACTIVIDAD PARA EL PERIODO DADO+4+5 SIMULACION4+5+1 DATOS PARA LA SIMULACIONA CONTINUACION SE PRESENTAN LOS DATOS QUE SE USARON PARALA SIMULACION DE ESTE CODIGO, ESTOS DATOS SE OBTUVIERON DE LAFIGURA 7i2 DEL LIBRO DE LAMARSH, INTRODUCCION A LA INGENIERIANUCLEAR +PARA EL MODO (i), SE ANOTARON REACTIVIDADES NEGATIVAS QUESE MUESTRAN EN LA TABLA 4ilj PARA SEIS GRUPOS DE NEUTRONES YCON DOS DIFERENTES TIEMPOS DE GENERACION DE NEUTRONESi(EST0SDATOS NO SE ENCUENTRAN EN LA FIGURA ANTES MENCIONADA, DEBIDO AQUE LOS VALORES SON NEGATIVOS, Y EN LA FIGURA SE TOMARONVALORES ABSOLUTOS)PARA EL MODO (ii), SE ANOTARON LOS PERIODOS QUE ESTAN EN LA66

FIGURA 7+2 ANTES MENCIQU@M PURA SEIS GRUPOS DE NEUTRONES YTRES DIFERENTES TIEMPOS DE GENERWION DEL NEUTRON, CON LO QUESE OBSERVO QUE EL COt>T’&-* F CfdNABA CORRECTAMENTE Y SEPROCEDIO A ANOTAR LOS DATOS PARA SIMULAR EL MODO ti), QUE NOHABIA FORMA DE COMPARAR SU RESULTADOS+4+5+2 GKAFICASLOS DATOS MOSTRADOS ANTERIOR MENTE FUERON GRAFICADOS CONAYUDA DEL CODOGO PLOT+ DONDE SE PUEDE OBSERVAR QUE LASGRAFICAS DE PERIODO CONTRA REACTXUIDAD SIGUEN ELCOMPORTAMIENTE DE LA FIGURA 7+2 DEL LIBRO (DEBE OBSERVAREiE QUELAS GRAFICAS DEL LIBRO ESTAN EL LOG-LOG), AS1 TAMBLEN LAGRAFECA DE REACTZVIDADES NEGATIVAS (QUE DA PERIODOSNEGATIVOS), ES SEMEJANTE A LA MOSTRADA, Y ES PRACTICAMENTE LAMISMcl CURVA PARA LOS TIEMPOS DE GENERACION UTILIZADOS PARA LASIMULACION+

PARA U-i GRUPOI0.1I 0.0065 IILLAIJDA (h) BETA ( p )0.0124 O. 000215O. 0305 0.0014240.1110 0,0012740.3010 O. 002568'\1,14CO 0.0007483 010 0.000273

~~1PEZLEC30 'PERI OD00.001 o. 00010.001 o. O001-8e-3 -83.11-4283 -0 . o1 -82.62 -82.62-4e-5 I -2180 -2159 -0.02 -81.57 -81.57-6e-5 I -1465 -1452 . -004 -81. og -81.09-8e-5 1' -1108-1097 -0.06-80.94-le-4 I - 894.0 - 865.6 -0.08 -80, &7 -80.87-2e-4 I - 466.6 - 462.0 -0.1 -80.82 -80.82-4e-4 I - 253.3 - 251.3 -0.2 -80.73 -m. 73-6.e-4 I - 183.7 - 182.7 -0.4 -8C.69 -eC.69- 8e-4 I - 149.9 -149.2 -0.6 -80.67-80.67-le-3 1 - 130.6 - 130.0 -0.8 -80.67 -V0.5?-2-3 I - 97.33 - 97.21 -0.1 el -80.66 -50.66- e6.68 - 2 -W,65-2- 84.24 -2 -80.65 -80.65

~~RIZ/iCTIVIDAf>rEACTI VI DADO, O01 o * O001 e e O01 0.0001.L . 82e-4 8,23e-5 8.20e-5 6.13e-3 5.23e-32 . 26e-4 1 o We-4 1 .'Ole-46.56e-3 5.44e-33 , 00e-4 1 . 35-4 1 . 34e-4I .6 I -1476 5 . 69e-3L45e-4 1.98e-4 1.95e-4 8,256-3 5 . 03e-35.63e-3 3.6Se-3 3.63-4 I 02 1.3374 0110 5.57e-3L . 64e-3 6-53e-4 6-44e-4 7.26~-3?.Ole-32.60e-37.75e-49.56e-47.643-4 1.08 I ,5584 . 0186 7.55e-30226 8,SCZ-3w3.73-3 1 . 26e-3 1 , 24e-3 7161 0306 8. e2e-36 . 82e-33. 90e-31.86e-3.1.02 A344 -0538o114O1 25. O1522.7Ce-32.98e-32.6Ge-32.87e-35032 C%S05584,116901%3 o 36e-33.2Le-36274.14e4o 02903 93- 34-962-33.70-34b5le-3. 71612Q520 3377

ANOTA LA HEACTIVIDAD:REACTIVIOAD = -.3000E+Ol PERIOOO 5 -.8065€+0;!W ( 1 ) = -.1240€-01W(2) -,3049E-01W(3) -.111OE+OOw(4) = -,3007E+00W(5) -,1140E+01W(6) -.301OE+OlW(7) = -,7517E+03ANOTALA REACTIVIDAD:REACTIVIDAD = -.4OOOE+Ol PERIODO = -.8065E+O2)W(1) = -,1240E-01WC2) = -,3049€-01W(31 = -,1110E+00ANOTA LA REACTIVIDAD:vJ(4) -,3008Et00W(5) -,1140E+01W(6) -,301OE+OlW(71 = -,8014€+03QUIERES ANOTAR QTROS DATOS(S/N)?ANOTA EL NUMFRO DE GRUPO DE NEUTRONES RETARDADOS:ANOTA LAS CONSTANTES DE DESINTFGRACION EMPEZANDO POR LA MAS PEQUENA:ANOTA LA FRACCION DE NEUTRONES RETARDADOS PARA CAOA GRUPOANOTA EL TIEMPO DE GENERACION DE NEUTRON:ANOTA O PARA CALCULAR LAS RAICES: 1 PARA CALCULAR LA REACTIVIDADDATOS ENTRADA Y DE SALIDAEL NUMERO DE GRUPOS DE FJEUTROlJES RETARDADOS ES:6L A b1B D A,1240E-01BETA,2150E-03,3050E-01 14114E-02111 OE*OO .I 1274E-02-3010E+00 .ZS68E-02-11 40E+01 ,7480E-03.30 1 OEtO 1 .¿'730€-03LA FRACCION TOTAL DE NEUTRONES RETARDADOS =EL TIEMPO OE GENERACION DE NEUTRON ES:SE EVALUARAPJ LAS HAICESTUS DATOS SON CORRECTOS[S/NI?ANOTA LA REACTIVIDAD:REACTIVIDAD = -.1000€-04 PERIODO = -.8524E+O4W(1) -,1173E-O3N(2) -,1435E-01Wt3) -,6825E-01W(41 -,1947E+00W(5) = -.1023E+Olh(6) : -,2894€+01,6502E-O2,1000E-03A .

I____' V " " 0....r....0....,....0...-r....O....r....O....r....0....,....0... . v . ' . 0W(7) =: -.6553€+02Wll) -,2335E-O3WlZ) = -.14366-014W(7) 3: -.6563E+O2W(2) = -,lU38E-01W(3) f -.6849E-0104ANOTA LA REACTIVIDAD:W(3) =W(4)-,6865E-01-.1952E+OOANOTA LA REACTIVIDAD:W(4) W(5) =-.1954E+OO-.1024E+01REACTIVIDAD = -o1000&-03 PERIODO = -.8856E+03W(53 = -.10Z4E+OlWC61 = -.2895E+OlREACTIVIDAD = -.¿?000E-03 PERIODO = -.4620E+03.... 0 ......... 0 . _ . . A . . . . 3-

REACTIVIDAD = -,YOOOE+01 PERIODO = -,8065E+02ANOTA LA REACTIVIDAD:W(1) -.124OE-O1W(2) -.30496-01W(3) = -.lllOE+OOWí4) = -,3008€+00W(5) = -.1140€+01W(6) = -,3010€+01W(71 = -.8013E+04QUIERES ANOTAR OTHOS DATOSIS/N)?ANOTA EL NUMERO DE GRUPO DE NEllTRONES RETARDADOS:ANOTA LAS CONSTANTES DE DESINTEGRACION EMPEZANDO POR LA WAS PEOUENA:ANOTA LA FRACCION DE NEUTRONES RETARDADOS PARA CADA GRUPOANOTA EL TIEMPO DE GENERACIOFJ OE NEUTRON:ANOTA O PARA CALCiILAR LAS RAICES; 1 PARA CALCULAR LA REACTIVIDADDATOS OE ENTRADA YDE SALIDAEL NUMERO DE GRUPOS DELAE46DA,1240E-01.3050E-01,lllOE+00.3010E+OQ, 114OE+O1. 3010E+O 1NEUTRONES RETARDADOS ES:bBETA-2150E-03, 14246-02. 1274E-O2,2568E-02,74806-03, 27 30E-03LA FRACCION TOTAL DE NEUTRONES RETARDADOS = ,6502E-02EL TIEMPO DE GENERACION DE NEUTRON ES: 1000E+00SE CALCULARA LA REACTIVIDAD:TUS DATOS SON CORRECTOS(S/NI?AEL PERIODO DELPERIODO DEL REACTOR = ,1000€+04 REACTIVIDAD= ,18196-03ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(l) 1:PERIODO DEL REACTOR = .8000€+03 REACTIVIDAD= ,2265E-03ANOTA EL PERIODO DEL REACTOR ET > -i/LAMBDA(il 1:PERIOOO DEL REACTOR = .h000E+03 REACTIVIDAD= , 3001E-03ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(lI 1:PERIODO DEL REACTOR = .4000E+03 REACTIVIDAD= ,4449E-03ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBüA(l) 1:PERIODO DEL REACTOR = .2000E+03 REACTIVIDAD= ,86256-03ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(ll 1:$3

PERIODO DEL REACTOR = .1OOOE+O3 REACTIVIDAD= .1641E-02ANOTA EL PERIODO DEL REACTOR (1 > -1/LAMBDAtll 1:PERIODO DEL REACTOR = .8OOOE+Oi! REACTIVIDAD= , 2 O 1 OE-02ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBOA(l) 1:PERIODO DEL REACTOR = .6000E+02 REACTIVIDAD= ,2602E-O2ANOTA EL PERIODO DEL REACTOR (T >PERIODO DEL REACTOR = .4000E+02ANOTA EL PERIODO DEL REACTOR tT >PERIODO DEL REACTOR = ,20OOE+O2ANOTA EL PERIODO DEL REACTOR (TPERIODO DEL REACTOR = .1000E+02ANOTA EL PERIODO DEL REACTOR' (T >PERIODO DEL REACTOR = 8000E+O 1ANOTA EL PERIODO DEL REACTOR CT >PERIODO DEL REACTOR = .6000€+01ANOTA EL PERIODO DEL REACTOR (T >-l/LAMBDA(l) 1:REACTIVIDAD=-l/LAMBDA(iI 1:REACTIVIDAD=-1/LAMBOA(11 1:REACTIVIDAD=-1/LAMBDA(11 1 :REACTIVIDAD=-1/LAMBDA(11 1:REACTIVIDAD=-l/LAMBDA(l) 1:,3725E-02.6820E-02,1247E=01 .. 1 S16E-0 1.1953€*0 1PER1000 DEL REACTOR.4000E+01REACTIVIDAD=,2798E-01ANOTA EL PERIODO DEL REACTOR (T >PERIODO DEL REACTOR =,2000E+01ANOTA EL PERIODO DEL REACTOR (T *PERIODO DEL REACTOR = ,1000E+01ANOTA EL PERIODO DEL REACTOR (T >PERIODO DEL REACTOR =.0000E+00ANOTA EL PERIODO'OEL REACTOR (T >PERIODO DEL REACTOR =.6000E+00ANOTA EL PERIODO DEL REACTOR IT >PERIODO DEL REACTOR = ,4000E+00ANOTA EL PERIODO DEL REACTOR (T >PERIODO DEL REACTOR =,2000E+00ANOTA EL PERIODO DEL REACTOR (T >-l/LAMBDA(l) 1:REACTIVIDAD=-l/LAMBDA(lI 1:REACTIVIDAD=-l/LAMBDA(l) 1:REACTIVIDAD=-l/LAMBDA(il 1:REACTIVIDAD=-1/LAMBDA(11 1:REACTIVIDAD=-1/LAMBDA(11 1 :REACTIVIDAD==l/LAMBDA(ll 1:,5187E-0 1,9557E-Oí. 1158€+00, 1476E+00.2046E+00.3374E+00

-_- _IC_-- -.- - --v c ‘I I * t * rPERIODO UEL REACTOR = ,1000E+00 REACTIVIDAD= .5031E+00.-. ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(l) 1:PERIODO DEL REACTOR = .8000€-01 REACTIVIDAD= .S584E+OOANOTA EL PERIODO DEL HEACTOH (T > -l/LAMBDA(l) 1:PERIODO DEL REACTOR = .6000E-01 REACTIVIDAD+ .6274E+00ANOTA EL PERIODO DEL ‘REACTOR (T > -l/LAMBDA(l) I:PERIODO DEL REACTOR = ,400OE-01 REACTIVIDAD= ,716lE+00ANOTA EL PERIODO DEL REACTOR (T > -i/LAMBDA(l) I:PERIODO DEL REACTOR = .2000E-01 REACTIVIDAD= , 8344E+00ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(lI I:PERIODO DEL REACTOR = .1000E-02 REACTIVIDAD= , 9902E+00ANOTA EL PERIODO DEL REACTOR (T > -1/LAMBDA(lI 1:PERIODO DEL REACTOR = .8000E-03 REACTIVIDAD= .9921E+00ANOTA EL PERIODO DEL REACTOR [T > -l/LAMBDA(l) 1:PERIODO DEL REACTOR = ,6000E-03 REACTIVIDAD= -994 1E+00ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(l) 1:PERIODO DEL REACTOR = .4000E-03 REACTIVIDAD= ,9960€+00ANOTA EL PERIODO DEL REACTOR (T > -1/LAMBDA(lI ):PERIODO DEL REACTOR = .rIOOOE-03 REACTIVIDAD: ,9980€+00ANOTA EL PERIODO DEL REACTOR (1 > -1/LAMBDA(l) I:QUIERES ANOTAR OTROS DATOS(S/N)?ANOTA EL NUMERO DE GRUPO DE NEUTRONES RETARDADOS:ANOTA LAS CONSTANTES DE DESINTEGRACION EMPEZANDO FOR LA MAS PEQUENA:ANOTA LA FRACCION DE NEUTRONES RETARDADOS PARA CADA GRUPOANOTA EL TIEMPO DE GENERACION DE NEUTRON:ANOTA O PARA CALCULAR LAS RAICES: 1 PARA CALCULAR LA REACTIVIDADDATOS DE ENTRADA YDE SALIDAEL NUMERO DE GRUPOS DELAMBDA,124OE-01,3050E-01, 11 1 OE+OO,301 OE+OO.lluOE+Ol.3010E+01NEUTRONES RETARDADOS ES:6BETA, 2250E-03, 1424E-O2.1274E-02,2568E-02,748OE-O3,2730E-03LA FRACCION TOTAL OE fJEUTRONES RETARDADOS.6502E-02

PERIDLA GRAFICA NORMALIZADA0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.+ - - - ~ ~ + - - o ~ ~ ~ - - - - ~ + ~ . ~ - ~ ~ ~ ~ ~ ~ ~ + ~ ~ o -0.0000 3 io.ooo(1 3iO.OGO0 3 10,0000 3 10*0000 3 19.0000 3 13*0000 a 10.0000 3 10.0000 3 13,0000 3 13,0000. 31 13,0000 31 13.0000 31 11.0000 3 1 1?*O000 3 1 11.0000 3 1 1,8000 3 1 1.6000 3 1 1,4000 3 1 1.2ooo 32 1 1.loo0 32 1 1.O800 32 1 1,0600 32 1 1.O400 !32 1 1.0200 !3 2 1 1.o010 ! 32 l!,00011 ! 32 1.O006 ! 3¿? 1.O004 !2 1,0002 !2 1+ - - - - - + ~ - ~ ' - + - - - - ~ + ~ - ~ - - + o ~ ~ o ~ + ~ - ~ ~ ~ + ~ ~ ~0.0 0.1 0.2 0.3 0.4 0.S 0.6 0.7 0.9 0.9 1.I1 2: REAC olOOE+Ol2 = .10UE+O13 = . 1 r30E+01Tx1xi?x3 x4 XS., 100E+04.800E+03.600E+03.400E+03.200E+03..600€+02. 4 o O'€ + o 2.2OOE+O2. 100€+02,800E+Ol.600€+01. 400E+01.2OOF+01l82E-03oZ27E-03.309E-03.445E-G3,863E-03160 20e26,373E-02. b82E-02125E-0 10 t 52E-0 1.I o102 o1. 51 9E-0 1. 19SE-03,368E-03.957E-03. 126F-02,191E-02,270E-02,298E-02,336E-,393E-Oí.4Q7E-0;!,820E-O4, 102E-03. 134E-03, 19SE-U3,364E-03,644E-03-7646-03,942E-03124E-O2. 18hE-i)2-261 €-O2,287E-o2,322E-02.371€-02.4S3E-O2* . . x.. .

~- _.l___l-.... . a , ..., . .(J...',.. ..o. ..',...'~.....".. 0 "".....O... .,." . u ......... o .REACTLA GRAFICA NORMALIZADA-1 0.9 -08 0.7 0.6 0.5 0.4 7.3 9.2 -01 O.+ o ~ ~ o . + . . . o . + o ~ " o o + ~ ~ o o ~ ~ ~ o o * ~ + o = o o o + . w o ~ . + * ~ ~ o ~ + ~ o o o w + ~ " o ~ . +-.oooo 12 I-.o000 ! 2 t*.ouoo ! 2 I*.o001 ! 2 1*.o001 ! 2 I=.o001 ! 2 i-.o002 ! ¿ !-.O004 ! i ? !-.O006 ! 2!-,O008 ! 2!..o010 ! 2!-.o020 t 2!=,O040 !c'!-.O060 ! 2!*.O080 ! 2!..o100 ! c7!'.O200 ! 2!..O400 ! I-!*,Oh00 ! 2!*,9800 ! 2!0 . 1000 ! 2!*.2000 1 2!*.4QQO ! 2!..6000 í 2!..8000 ! ;>!.o000 ! i!'.O000 ! ¿!;,o000 ! 2!i.0000 ! 2!+o~~o~+.~o.*+o~o"*o+*~oo"+o~oo.+..~o~+ooooo+o~~w~+-1 9.9 -08 -.7 0.6 -0s -.4 0.3 -02 -01 0.1 = PERT AC E N = .86iE+04Z = PER1 AC E PJ = 61E+04T x1 X2 x3 x4 Y5-,100E-04-,200E-04-.400E-04-,6OOE-.800E-.100E-03-.861€+04 -.852E+O4-0433 04-0218 04-0147 04 .0.111 04-.894E+03 -,88hE+03-.U67€+03 0.-.-,1OOE-O2 -.131€+03 -,130E+03200~902 -*97?€+02 -.97i?Et02O. 4OOE-O2 -.8b7E*üZ -,867E+O200600E-02 -.842€+02 -.B42E+Oil-,800E*Oc7 *.832E+02 -.8-.100E-01-.2OOE-O1-.826E+Ot-.BlhE+O2=.8-.81óE+O213')

4*6 CONCLUSIONESSE COMPROBO QUE EL CODIGO ESTUVIERA FUNCIONANDOCORRECTAMENTE j AL DARLE DIFERENTES VALORES DE REACTIVIDADjPERIODO, Y TIEMPO DE GENERACION DEL NEUTRON, CONOCIENDO LOSRESULTADOS Y COMPARANDO CON LOS OHTENIDOS POR EL CODIGO,TENIENDO EN CUENTA SU FUNCIONAMIENTO CQRRECTOj SE REALIZARONGRAFICAS (USANDO EL CODIGO PLOT) DE RACTIUIDADES NEGATIVASCONTRA PERIODOS NEGATIVOS PARA DOS TIENPOS DE GENARACLON DELNEUTRON +

CAPITULO 5+ CODIGO ENRICHSil PROPOSITOS DEL PROGRAMAENRICH ES UN CODIGO ESCRITO EN FORTRAN PARA CREAR TABLAS DELOS REQUERIMIENTOS DE LOS COSTOS DE ALIMENTACION DE MATrRIALPARA EL ENRIQUECIMIENTO DEL URANIO POR EL MFTciDO DE DIFUSIONGASEOSA +5.2 TEORIA EN QUE SE BASALOS REACTORES DE ACbA LIGERA REQUIEREN QUE EL URANIO-235CONTENIDO EN EL URANIO NATURAL (Ot711X) SEA INCREMENTADO AL2-4 X, TP'.- ENRIQUECIMIENTO SE HACE FRECUENTEMENTE POR DIFUSIONGttSEOSA DEL HEXAFLUORURO DE URANIO A TRAVES DE MILES DE FASESQUE CONTIENEN BARRERAS POROSAS+ UNA PORCION MAYOR DEL COSTO DELOS ENSAMBLES DEL COMBUSTIBLE NUCLEAR PARA LA PRODUCCION DEPOTENCXA NUCLEAR COMERCIAL ES ATRIBUIBLE A LA ALIMENTACION DELMATERIAL Y LOS REQUERIMIENTOS DE TRABAJO DE SEPARACIONf LACUBIERTA DEL COMBUSTIBLE NUCLEAR ENRIQUECIDO DENTRO DE TUBOSDE CIRCONIO, SU COSTO ES SORPRENDENTEMENTE BAJO,REQUERIMIENTOS DEL COMBUSTIBLE PARA REACTORES DEAGUA LIGERAPOR SUPUESTO HAY MUCHOS DISENOS POSIBLES DE REACTORES QUEREPRESENTA DIFERENTES MEDIOS DE CUMPLIR EL MISMO FIN SENALADO,USANDO LA ENERGIA LIGERADA POR LA FISION NUCLEAR, PARAPRODUCIR VAPOR Y CONDUCIRLO A LAS TURBINAS-GENERADORAC PARALA PRODUCCION DE POTENCIA ELECTRICA, LAS VENTAJAS YDESVE:NTAJAS DE UN DISENO PARTICULAR DE REACTOR, NO DEPENDEN DECOMO LOS COSTOS DEL CICLO DE COMGUSTIBLE SON ESTIMADOS,PARA ILUSTRAR LA INTER-RELACION DE LOS COSTOS DEL CICLO DECOMBUSTIBLE, CONSIDERAREMOS EL BALANCE DE MATERIAL TIPICO DELDIAGRAMA DE FLUXJO MOSTRADO EN LA FIGURA(S+l)+ LA RELACION DEFLUJO ESTA DADO EN TONELADAS WETRICAS POR ANO DE CONTENIDO DEURANIO+ L.0S REACTORES DE AGUA A PRESTON SUPONEN QUE PARAOPERAR El... REACTOR AL 85 X DE FACTOR DE POTENCIA, ES DECIR A850 MNe POR TODO EL ANOtDE 1000 Ml(e AL ANO* LA EFICIENCIATERMICA ES CONSIDERADA COMO 1000/3077 O 32,s %+ EL REACTOR DADE 40+782 TON O 46,782 KG DE URANIO AL 2+548 Y. DEENRIQUECIMIENTOfLA DESCARGA DE 45,448 KG DE URANIO AL 1,016 XDE ENRIQUECIMIENTO +NOSOTROS DESEAMOS CALCULAR EL MATERIAL ALIMENTADO Y TRABAJOREQUERIDO QUE MUESTRA EL DIAGRAMA DE FLUJO (DE LA FIFURA S*I),DESDE QUE ESTOS REQUERIMIENTOS SON ANOTADOS PARA LA ESTIMACJONDEL COSTO DEL COMBUSTIBLE NUCLEAR, LA RESPUESTA QUE NOSOTROSBUSCAMOS PARA CORROBORARLO SON:145j659 KG DE URANIO NATURAL ALIMENTADO A LA DIFUSIONGASEOSA EN CASCADA, Y 139+38 TON POR ANO DE TRAEAJO SEPARATIVO

fr’ *-. ---

EFECTUADO POR EL PROCESO DE ENRIQUECIMIENTO+ENRIQUECIMIENTO POR EL METODO DE DIFUSION GASEOSA+EL METODO DE DIFUSION GASEOSA PARA LA SEPRACION DE LOSISOTOPOS HACE USO DE LAS RELACIONES DIFERENTES A QUE LOS GASES(O VAPORES) DE DIFERENTES PESOS MOLECULARES SE DIFUNDEN ATRAVES DE UNA BARRERA POROSA+ LA BARRERAS USADAS EN EL PROCESODE SrmwAcIoN DEL URANIO CONTIENEN MILLONES DE POROS PORMILIMETRO CUADRADO, EL DIAMETRO PROMEDIO DE LOS POROS SONAPROXIMADAMENTE 0+05 + COMO SE SABE LOS ISOTOPOS MOLECULARESLIGEROS SE DIFUNDEN A TRAVES DE TALES BARRERAS MAS RAPIDAMENTEQUE L.0S ISOTOPOS PESADOS, AS1 QUE UNA SEPARACION PARCIAL DEL.0S ISOTOPOS OCURRE+ EL GAS QUE PASA PRIMERO A TRAVES DE LABARRERA SERAN RELATIVAMENTE MAS LIGEROS, CONSIDERANDO EL GASSOBRANTE to QUE NO PASA) CONTIENE UNA PORCION MAYOR DEISOTOPOS PESADOS,ELUNICO COMPUESTO UTIL DE URANIO PARA SU USO EN EL PROCESODE DIFUSION GASEOSA, ES EL HEXAFLUORUROI AUQUE ESTE COMPUESTOES SOLIDO A TEMPERATURA AMBIENTE, SE VAPORIZA (SUBLIMA)FACILMENTE APROXIMADAMENTE A UNA TEMPERATURA DE 56 GRADOSCENTIGRADOS; LAS PLANTAS DE DIFUSION OPERAN A TEMPERATURAS YPRESIONES TALES RUE EL HEXÁFLUORURO DE URANIO ESTE EN LA FASEGASEOSA (O VAPOR)+ PARA LA DIFUSION GASEOSA, EL HEXAFLUORUROTIENE AI ... GlINAS VENTAJAS t PRIMERIO , EL PESO ATOMIC0 DEL FLUOR ESRELATIVAMENTE PEQUENOj ASS QUE EL PESO MOLECULAR DE LOS GASESA DIFUNDIRSE, RUE ES INEVXTAELEMENTE ALTO PORQUE EL PESOATOMICO DEL URANIO ES GRANDE, NO SERA DEMASIADC) GRANDE+ ASI,CUANDO El ... PESO ME1 . ECUL..AR NC) ES DEMASIADO GRANDE , HAY UNA MAYOREFICIENCIA DE SEPARACION EN ESTAS CONDICIONES+ SEGUNDO, ELF.I...IJC)R CONSTA DE UN SOLO ISCITOF'O, POR LO TANTO 9 EN EL. PROCESC)DE: DIFUS.:ION SOLO SON SI-PARAD(lS 1-0s ISOTOPOS DE URANIO + ADEMAS 9El ... t.IEXAI::'l.-lJORURO DE lJl?ANI(l ES UN COMPUESTO ESTAE:LE + !31NEMBARGO ESTA MEZCL.A REACCIONA MUY FACIL.MENTE , Y ES MUYCíIWR(J8IVA PARA LA MAYCIRIA DE LOS METALES Y ATACA A MATERIALES~RGAN:X:COS TAI... ES COMO ncE:rrE:s LUBRICANTES + ESTA ACT:IVIDADR u :I M I c: n 0 c ASIQNA EL. USC) DE QUE LOS PíATERIAL..ES ESTRUCTURALES PARA LOSC~í3MI::'ONIEN'TES DE L.AS F'L.AN'TAS DE DIFUSION SEIAN DE LOS METAI...ESNX(1iJEi:i ... Y AI..lIMINJO t

I...A REL.ACSON DE TRANSPORTE DE UNA ESPECIE MOLECULAR ATRAVES DE UNA RARREHA PQI3QSA POR DIE-USION ES PROPORCIONAL A LAFRECI.lENCIA CON QIJE I...AS MOLECUL.AS DE ESAS ESPECIES GOLPEAN LARARRERA, Y ESTA FHE:CUENCIA ES PROPORCIONAL AL INVERSO DE LAVELOCIDAD PRQMEDIO DE LAS MOI.,ECULAS+ ENTONCES EN UNA MEZCLA DEDOS ESPECIES MOLECULARES TALES COMO U(235)Fb Y U(2381F6, EN ELHEXAFI.,IJOR\JRQ DE URANIO NATURAL, EL ISOTOPO MOLECULAR LIGEROI.I(235)Fili, CON LA VELOCIDAD PROMEDIO MAYOR, TIENDE A DIFUNDIRSEA TRAVES DE UNA BARRERA POROSA MAS RAPIDAMENTE QUE LOSISOTOPOS MOLECULARES PESADOS U(238)F6+ DE ESTA MANERA SE PUEDEREALIZAR UNA SEPARACION PARCIAL DE LOS ISOTOPOS DE URANIO+UNA INDICACION DEL GRADO DE SEPARACION ( O ENRIQUECIMIENTOISOTROPIC0 1 QUE PUEDE SER REALIZADO POR DIFUSION A TRAVES DEUNA BARRERA PUEDE SER OBTENIDA DE LA SIGUIENTE MANERA, SI MHES EL PESO MOLECUL.AR DE LA ESPECIE PESADA U(238)Fót ES DECIR352 UMAS, Y VH ES SU VELOCIDAD MOLECULAR PROMEDIO, Y SI ML ESEL PESO MOLECULAR DE LA ESPECIE LIGERA U(235)Fbt ES DECIR 349UMAS, Y VL ES SU VELOCIDAD MOLECULAR PROMEDIO, ENTONCES LASENERGIAS CINETICAS EN LA MEZCLA A UNA TEMPERATURA UNIFORME SONIGUALES, POR LO TANTOLA RELACION DE LAS VELOCIDADES MOLECULARES DE DOS ESPECIESISOTROPICAS, QUE ES IGUAL A LA RELACION DE LAS VELOCIDADESRESPECTIVAS A TRAVES DE LA DIFUSION DE LA BARRERA POROSA, ESLLAMADA EL FACTOR DE SEPARACION TEORICA d* ; AS1DEESTE FACTOR ES CERCANO A LA UNIDAD, QUE INDICA QUE EL GRADOENRIQUECIMIENTO DE URANIO-233 COMO RESULTADO DE LA DIFUSIONA TRAVES DE UNA BARRERA ES MUY PEQUERA, POR LO TANTO PARATENER UN ENRIQUECIMIENTO SUSTANCIAL ES NECESARIO EMPLEAR UNA"CASCADA" QUE CONSISTE DE MUCHAS ETAPAS CON UN INCREMENTO DEENRIQUECIMIENTO DE UNA ETAPA A LA SIGUIENTE+ UN GRAN NUMERO DECOMPRESORES SON NECESARIOS PARA HACER CIRCULAR EL GASHEXAFLUORURO+ UN GRAN NUMERO DE INTERCAMBIADORES DE CALOR SONNECESARIOS PARA MOVER EL CALOR DE LA COMPRESION, AS1 LASPLANTAS DE DIFUSION GASEOSA SON COMPLEJAS Y MUY GRANDES,CARACTERIZADOLAS LA ENORME DEMANDA DE POTENCIA Y AGUA FRIA,LA OPERACION EFICIENTE DE UNA PLANTA DE DIFUSION GASEOSAOBLIGA A UNA ELECCION PARA UNA COMFOSICION DE COLAS (OMATERIAL DESECHADO)+ LA FIGURA 1 MUESTRA 143,234 KG DEMATERIAL AGOTADO DESECHADO POR LA DIFUSION GASEOSA EN CASCADA83

i.. . ..I ..., . .-. ..1 ., . .. ...,. , .. .. . . . . o U-235 0 .U-238F1Gm.A 5.2 UNA ETAPA 32L PROCZSO 33 SZFAMCICIjPO3 3IF'JSION GASEOSA..-- . . .- - , . -a-----. .i . '. . .FIG'J2A 5.3 PRGCZSC 92 SZPAZ43IOM ?ORL)IFU3IOX GASEOSA.8W.U

Ecuacion del Balance de Materiales.F = P + WXf F = Xp? + XwiVFIGUHA 5.4 RBCICLO lQ? CASCADA.85

A UNA COMPOSICION DE COLAS DE 0+2 X DE URANIO-239+ LA PLANTADEBE SER OPERADA AL 011 X DE COLAS, RECUPERANDO UNA GRANFRACCION DEL ISOTOPO VALIOSO, PERO EL COSTO DE TRABAJO DESEPARACION ES ADICIONAL+SE INVOLUCRAN DOS BALANCES ESENCIALES EN LA ESTIMACION DELOS COSTOS DE ENRIQUECIMIENTO+ UN BALANCE DE MATERIAL PARA ELURANIO-239, PARA ESTABLECER LA RELACION DE ALIMENTACION (F)-- :. PRODUCTO (Pi, DE LAS COMPOSICIONES DE LA ALIMENTACION(XF), PRODUCTOS (XP) Y DESPERDICIOS (XI.0, CONVENCIONALMENTEEXPRESADO EN PORCENTAJE DE PESO DE URANIO-235+ASI, SI SOLO URANIO NATURAL SE ALIMENTA EN LA CASCADA DE LAFIGURA(l)+PARA PRODUCIR 47,284 KG DE PRODUCTO SE REQUIERE ALIMENTAR4+595(47,284)= 217 9 266 KG DE URANIO NATURAL+ LOSREQUERIMIENTOS ACTUALES DE LA ALIMENTACIQN SON MENORES QUEESTO PORQUE DE 44,899 KG DE COMBUSTIBLE RECICLADO POR PROCESOQUIMICO +iILOS EFECTOS DE TALES FLUJOS RECICLADOS PUEDEN SER IMPUESTOS ICALCULANDO EL PRODUCTO ENRIQUECIDO DE URANIO QUE DEBA SERPRODUCIDO DE LA ALIMENTACION AL RECICLO+POR ESTO, P=44,859/2+877-19,592 KG+ LA CANTIDAD NETA DEURANIO NATURAL ALIMENTADO REQUIERE UNA REDUCCION 9CONFIRMANDO LOS 145,659 KG QUE SE MUESTRAN EN EL DIAGRAMA+UN SEGUNDO BALANCE -UN BALANCE DEL VALOR DE L.A FUNCION- ESNECESARIO PARA EVALUAR LOS REQUERIMXENTOS DE TRABAJOSEPARATIVO PARA PRODUCIR URANIO ENRIQUECIDO DE UNO O MASFLUJOS ALIMENTADOS+ EL VALOR DE LA FUNCION SE ESCOGE ENFUNCION DE : V(X)=(2X-I) LN(X/l-X')+ ASI, PARA CUALQUIERCOMPOSICION X, EN FRACCION DE PESO O URANIO-235, HAY UNNUMERO, V(X), QUE REPRESENTA NUMERICAMENTE EL VALOR DE LAFUNCION ASOCIADA CON SU COMPOSICION+ El. VALOR DE LA FUNCION NOREPRESENTA EL VALOR POR UNIDAD DE MATERIAL, DESDE QUE ESTAFALTA INCLUYE LA CANTIDAD DE MATERIAL ALIMENTADO Y ENERGIA

CONSUMIDA, YEL COSTO DEL EQUIPO UTILIZADO*LAS UNIDADES DE TRABAJO SEPARATIVO (SWU) DELIBERADAS PORUNA PLANTA DE DIFUSION, SON DEFINIDAS COMO EL INCREMENTO EN ELVALOR DE LA FUNCION, DEL FLUJO EFLUENTE DE LA CASCADARELATIVA, PARA EL FLUJO ALIMENTADO A LA CASCADA+ PARA UNAPRODUCCION EN CASCADA DE 2*548 X DE PRODUCTO Y +2% DEDESPERDICIO DEL +711 ALIMENTADO:DESDE QUE EL VALOR DE LA FUNCION ES DIMENSIONADA, SWU TIENELAS MISMAS UNIDADES COMO P+W Y F (KG DE URANIO)*REESCRIBIENDO, PARA UNA UNIDAD DE PRODUCTO (KG), tN EL CUALW=F-P:SWU/F' = U(XP)+(F-P/F)mU(XW)-F/PxV(XF)DONDE: XP=0+02548 ; V(XP)13*43850XF10+00711 ; U(XF)=4+86888Xbl~0+00200 ; V(XW)=6+18776RECORDANDO QUE F/P TIENE VALOR DE 41595 PARA ESTASCOMPOSICIONES,SWU/PiB + 33099 I

USANDO ESTOS DOS BALANCES ESCENCIALES -EL BALANCE DEMATERIAL Y EL VALOR DE LA FUNCION DEL BALANCE- NOSOTROSPODEMOS ESTIMAR EL MATERIAL ALIMENTADO Y REQUERIMIENTOS DELTRAE:AJO SEPARATIVO PARA ENRIQUECIMIENTO DE URANIO POR ELMETaDo DE IIIFUSION GfiSEasA, EL COSTO ENTONCES SE CALCULASIMPLEMENTE MULTIPLICANDO LA ALIMENTACION REQUERIDA POR LOSCOSTOS POR UNIDAD DE ALIMENTACION, Y EL TRABAJO SEPARATIVOREQUERIDO POR LA UNIDAD DE COSTOS DE TRABAJO SEPARATIVO:COSTOS= F($ ALIMENTACION) + (SWU($ SWU)LOS COSTOS REPRESENTATIVOS ACTUALES SON DE s5*00 roR KGDE URANIO NATURAL ALIMENTADO Y $ 36,OO POR KG SWU+ASI, EL ENRIQUECIMIENTO DE HEXAFLUORURO DE URANIO UTILIZADOPARA LA FABRICACION DE COMBUSTIBLE ES CONSIDERADO $ 8,659,000/47,284 O APROXIMADAMENTE $ 180*OO POR KG DE URANIO CONTENIDOTENIENDO CREDIT0 INMEDLATAMENTE PARA RECICLAR UN FLUJO QUEPUEDE SER ACTUALMENTE APLAZADO EN UN TIEMPO DE 5 ANOS, DESPUESDE QUE EL REACTOR SE CONSTRUYA* EN UNA PLANTA DE DIFUSIONGASEOSA, LA DISPONIBILIDAD DE URANIO BARATO FAVORECE ALTASCOMPOSICIONES DE COLAS, LA DIPONIBILIDAD DE POTENCIA BARATA YBAJAS CARGAS DE CAPITAL FAVORECE BAJAR LAS COMPOSICIONES DECOLASI HA RESULTADO UN PROGRAMA POLITICO COMUN DE IMPONER PARAEL TRABAJO SEPARATIVO CON LA BASE DE UNA COMPOSICION DE COLASDEL 0,s '31EL CALCULO DE FASES, SE BASA EN EL FACTOR DE SEPARACION DE0,004 roR CADA FASE*513 DIAGKAYA DE FLUJOEN EL DIAGRAMA DE FLUJO SOLO SE PRESENTAN LAS OPERACIONESPRINCIPALES, PERO PARA MAYORES DETALLES CON RESPECTO A ESTECODIGO VER EL APENDICE D, DONDE SE ENCUENTRA EL LISTADO+~". ' ..-.. . , ,., .

I'I

5*4 USO DEL PROGRAMA (ENTRADAS Y SALIDAS)*E N T R A D A S===============ESTE CODIGO CALCULA HASTA 10 TABLAS PARA 10 CASOSDIFERENTES, QUE EL USUARIO DEBE ANOTAR*EL USUARIO DEBE ANOTAR (TODO EN FORMATO LIBRE)*- EL COSTO DE MATERIAL ALIMENTADO Y EL COSTO DE UNA UNIDADDE TRABAJO SEPARATIVO (SWU) EN $/KGtEL ENRIQUECIMIENTO QUE TIENE EL URANIO ALIMENTADO (1%=1),Y EL NUMERO DE CASOS QUE SERAN CORRIDOSt (ESTOS DOS DATOS ENLA MISMA LINEA)*- EL ENRIQUECIMIENTO INICIAL Y FINAL, EL INCREMENTO DEENRIQUECIMIENTO, Y EL ENRIQUECIMIENTO DE LAS COLAS, PARA CADACASO (ESTOS TRES DATOS EN LA MISMA LINEA)*- UNA VEZ QUE SE HICIERON LOS' CALCULOS, EL CODIGO PREGUNTA,"SI QUIERES ANOTAR OTROS DATOS (S/N)?", SI EL USUARIO ANOTA"S" (S=SI), EL CODIGO PREGUNTA DE NUEVO TODOS LOS DATOS, ENCASO CONTRARIO "N" (N-NO), EL CODIGO TERMINAtCADA TABLA CONTIENE LA INFQRMACIQN DEL MATERIAL ALIMENTADO,LA UNIDAD DE TRABAJO SEPARATIVO Y LOS COSTOS DE ALIMENTACION,UNIDAD DE TRABAJO SEPARATIVO Y TOTALES EN FUNCION DE UN KG DEURANIO Y EN FUNCION DE URANIO-235, ADEMAS EL ENRIQUECIMI€NTODE LAS ETAPAS Y LAS ETAPAS TOTALES, PARA ALCANZAR ELENRIQUECIMIENTO DESEADO, CON EL INCREMENTO DE ENRIQUECIMIENTOEL CUclL ESPECIFICA EL USUARIOt5t5 SIMULACION5t5.1 DATOS PARA LA SIMULACION Y RESULTADOS,A CONTINUACION SE PRESENTAN LOS DATOS QUE SE ANOTARON ALCODIGO, AS1 TAMBIEN COMO LOS RESULTADOS QUE SE OBTUVIERONtS*St2 GRAFICASEN ESTA SECCION SE PRESENTAN LAS GRAFICAS REALIZADASTOMANDO COMO DATOS LOS RESULTADOS OBTENIDOS EN LAS TABLASANTERIQRESt LAS GRAFICAS MUESTRAN UNA VARIACION DE LOS COSTOS

DE LA UNIDAD DE TRABAJO SEPARATIVO, LOS COSTOS EN KG DE URANIOY LOS COSTOS TOTALES CONTRA EL ENRIQUECIMIENTO DE LAS COLAS.OBSERVANDOSE QUE LOS COSTOS DE ALIMENTACTON AUMENTAN ALAUMENTAR EL ENRIQUECIMIENTO DE LAS COLAS, LOS COSTOS DE LAUNIDAD DE TRABAJO SEPARATIVO TIENEN EL COMPORTAMIENTO INVERSOQUE LOS COSTOS DE ALIMENTACION, ENTONCES LOS COSTOS TOTALES(QUE ES LA SUMA DE LOS COSTOS DE ALIMENTACION Y DE LA UNIDADDE TRAEAJO SEPARATIVO), TIENEN UN COMPORTAMIENTO QUE TIENE UNPUNTO MINIMO, QUE INDICA EL OPTIMO DE LOS COSTOS ANTESMENCIONADOS.

I . “ ‘ O....r....O....T....O....r....O,...v....O... ,..... O .... .... o...‘.MATEf?ftCiiCTCIi

dO444u3tT W'Uz aa c:+Ym u3i.iuu3-7oc'/?3r,99c)aJrI-.IUc:73!-I4J

14OI-hv>@4ot- 2dUUUW3Ciy7w.-Uaaf-lzr,3u>LzwmWCm3!-I 4+I313xO4cnww2Qcnc)f-cn3I-o3 H$W4xUziOYVw>YWf-Q3CKbVLf3 3t3'BaJ3'J-IW3a-,4I--.-Ck33zc3wh,E?33 z:n .3J'Oac.4I.-cb-x3L)cc-2. -r(3.Jxu-Ov>occir3oaJhv;MPinMryNt3a .13aY :.

ICU.c:-.b+,b43Ob4atbciaunuaZtww4aCb3.............MNMl-bMm5M00aOmO=f+OMmbmdNMLncJG.U'U'000004~d+""MMMMMMMMM4Cb3L n nmI-MV)N0-VIL4cb4c33(3YoOtbd4w3bDb1toO Wt- tu3Q3 a>l-ltL)w az wwa2:Ua3HuH>Ywmw ar% Qtzc4C4CbVIMeI4Db'n.4J9. uVCIa wz ac*LLObO............. -. f T..-,+n,vc,330OOiO~33000h < L- ..-UQmpriNI3.*..................... .--.0'3 r=) a o 3 o o 3 o 3 .a o 0 o o o o o o. ..._ .. -- - Ir> ,.A e,*-u C'k;r)9axouUA

O44o4c43oY3t-zwOI 4c .-cl Huv: w w r)c3v) HO WI- !-- o^,Lu w 3c3ao:eaXaCL3UuU>iyuwC-A-H@A-d n9-ULLC O n3I-aCl) Yo-ULLc42 .'3 cY .o .a .- 4m .w .3 ..n9Jr c - 3u.o w7 D2 ..-3 JLL __ -cac2 L3....................3 I> 3 3 O *2 G 3 3 0 O O O 3 3 3 O 3 0 Oi-.I n, M rt r, h n ,P -> - n. m t in .I\ A CT

c->b43cb0InM3b3bab3c43c3b1bE3c33'3y:cdcnw# U 3o W?- bo z3 Uz4a>ulMNI-3I-3CYl>t(I-UWav)wnLPMNI3>3uYw3c3'3 z%, 3a 4u uc'J:c)ozoHO30300000 ~000~0033.................. .. - -- . .--- ..

.. - .. .-. . .~ .. , . . . o . . . .. . . . () . .' Y ' ' ~! - v . " J ' vt.loo0 ! O I:) D.2000 1 o O o !,3000 ! o O D I.YO00 ! 0 o D IO = COSTO DE ALIM. (KGU) FACTOR DE NORMALIZACION = , 155E+050 = COSTü DE SVJU (KGU) FACTOR DE NORMALIZACION = .lSSE+OS0 = COSTO TOTAL (KGU) FACTOR DE NORMALIZACION = . iSSE+OST X 1 X2 X3 x4 x5 X6 !.100E+00 .162€+04 .138E+05 .155F+O5.200E+00 .192E+O4 .104€+05 .123E+OS,300€+00 ,236€+04 .864€+04 .11OE+O5.400€+00 .3Oc)E+04 .747E+04 .106E+05.SOOE+00 .Y50E+OY .6h2F+04 .111E+OSQUIERES ANOTAR OTROS DATOS (SIN)?ANOTA EL #DE CtJRVAS Y EL NUMERO DE PUNTOS PARA CADA CURVA:ANOTA LA CONSTANTE DE NbRMALTZACION PARA CADA CURVA:ANTITA LA OPCIürJ PARA LOCALIZAR EL ORIGEN:ANOTA EL CETRERO PARA EL EJE T MAX CINCO LETRAS.:ANOTA EL NIVEL PARA LA LONGITUD DEL EJE(1"60;2~1001:ANOTA EL LETRERO PARA CADA CIJRVA (MAX 20 LETRAS):ANOTA XI,1 --->.r..i XE.I, 12 o-->3 --->4 -o->!j 99-w--o>7 --->a --->9 -o->$0 .-o>:RE3 MODIFICAR ALGUN PUNTO(S/NI?

1'0 o o ! I1 D.1500 ! O 0 O,2000 ! 0 O 0.25riu ! í7 o O,3000 ! O O 003500 ! 0 O D.4000 ! 0 o o,4500 ! O O O f.5oou ! !I n t0 = COSTO ALI'.l. (KGU) FACTOR DE NIJRMALIZACION = .2056+04O = COSTO S24U (KGIJI FACTOR DE NORMALIZACION = .205€+04O COSTO TOTAL {KGUI FACTOR DE NORMALIZACION = .20SE+04T X1' XS x3 x4 xs X6.100E+OO .'i7OE+O3 .598E+03 .117E+04.150E+00 .610E+03 .498E+O3 .lllE+04m200EtO0 .6SRE+03 .431E+03 ,109€+04.ZSOE+OO .7166+03 ,381€+03 .llOE+O4.300E+00 .788€+03 .342E+03 .113€+04.3SOE+0ü .8816+03 .31iF+O3 .119€+04.400E+00 . 100E+u4 .285E+O3 .129€+04.450E+00 .117€+04 .262Et03 .143E+OY.SOOE+00 .412E+04 .243€+03 .166E+O4.55OE+00 .183€+04 .226E+03 .205€+04QUIERES ANOTAR OTROS DATOS (C/M]?ANOTA EL #DE CIIRVAS Y EL NUMERO DE PUNTOS PARA CAUA CURVA:ANOTA LA CONSTANTE DE NURMACIZACION PARA CADA CURVA:ANOTA LA OPCIOFJ PARA LOCALXZAR EL ORIGEN:ANOTA EL CETRERO PARA EL EJE T MAX CINCO LETRAS.:ANOTA EL bIVEL PARA LA LONGITUD DEL EJE(1"60;2"100):ANOTA EL CETRERO PARA CADA CURVA (MAX 20 LETRAS):ANOTA X1, ..... Xf.1, T1 -*o>2 --o>3 ---a4 --->r, --->

5+6 CONCLUSIONESOBSERVANDO LAS TABLAS DADAS POR ESTE CODIGO SE PUEDECONCLUIR LO SIGUIENTE,LOS COSTOS TOTALES, LOS COSTOS DE UNA UNIDAD DE TRABAJOSEPARATIVO (SWU) Y LOS COSTOS DE ALIMENTACION EN FUNCION DEU-235, VARIAN GRANDEMENTE PARA BAJOS ENRIQUECIMIENTOS(APROXIMADAMENTE AL 7%) Y AL SEGUIR INCREMENTANDO ELENRIQUECIMIENTO DEL URANIO, LOS COSTOS VARIAN MUY POCO, ENCAMBIO LOS COSTOS TOTALES, DE ALIMENTACION Y DE SWU EN FUNCIONDE KG-U, TIENEN UNA UARIACION QUE SE PODRIA CONSIDERAR LINEAL,CON RESPECTO AL ENRIQUECIMIENTO DEL URANIO, AS1 TAMEIEN COMOLOS KG DE URANIO ALIMENTADO Y LAS UNIDADES DE TRABAJOSEPARATIVO, TODO ESTO CUANDO EL ENRIQUECIMIENTO DE LAS COLASSE MANTIENE CONSTANTE, A UN PORCENTAJE DE U-235 EN EL MATERIALALIMENTADO+AHORA CUANDO CUANDO SE VARIA EL ENRIQUECIMIENTO DE LASCOLAS, PARA UN MISMO PORCENTAJE DE ENRIQUECIMIENTO DEL URANIODE SALIDA, AS1 COMO EL DEL URANIO ALIMENTADO, LOS KG DE URANIOALIMENTADO AUMENTAN Y LA UNIDAD DE TRABAJO SEPARATIVODISMINUYE CON FORME SE AUMENTA EL ENRIQUECIMIENTO DE LASCOLAS, POR LO TANTO LOS COSTOS DE ALIMENTACION Y DE LA UNIDADDE TRABAJO SEPARATIVO AUMENTAN Y DISMINUYEN RESPECTIVAMENTE,QUE ES EL COMPORTAMIENTO ESPERADO, YA QUE A MAYORENRIQUECIMIENTO DE LAS COLAS, SE REQUIERE ALIMENTAR UNACANTIDAD MAYOR DE URANIO PORQUE SE PIERDE UNA BUENA CANTIDADEN EL MATERIAL DE DESECHO, Y LA CANTIDAD DE UNIDADES DETRARAJO SEPARATIVO DISMINUYEN PORQUE EL ENRIQUECIMIENTO DE LASCOLAS AUMENTA (POR LA EC+ SWU = b(U(XP) - FU(XF)) PERO CONRESPECTO A LOS COSTOS TOTALES, ESTOS DISMINUYEN HASTA UNPUNTO, A ESTE PUNTO DONDE LOS COSTOS TOTALES SON MINIMOS LEPODRIAMOS LLAMAR PUNTO OPTIMO CECONOMICAMENTE HABLANDO) DELENRIQUECIMIENTO DE LAS COLAS, DESPUES AUMENTAN AL AUMENTAR ELENRIQUECIMIENTO DE LAS COLAS (VER GRAFICAS)+

c $1 I::':x: 'I' i.1 I... i:1 6 *

CAPITULO 6 + CODOGO STAT6+i PROPOSITO DEL PROGRAMASTAT ES UN CODIGO PARA COMPUTADORA ESCRITO EN FORTRAN PARAHACER UN ANALICIS ESTADICTICO DE LAS MEDIDAS RADIACTIUASI PARAESTE ANALISIS SE SUPONE QUE LAS MEDIDAS ESTAN DISTRIBUIDAS DEACUERDO A LA DISTRIBUCION NORMAL+6*2+ TEORIA EN QUE SE BASA EL CODIGO STAT+EL DECAIMENTO RADIACTIUO ES UN PROCESO ALEATORIO*CONSECUENTEMENTE, CUALQUIER MEDIDA QUE SE BASA EN LAOBSERUACION DE LA RADIACION EMITIDA EN DECAIMIENTO RADIACTIUOESTA SUJETA A ALGUNA FLUCTUACION ESTADISTICA+ ESTASFLUCTUACIONES INHERENTES REPRESENTAN UNA FUENTE INEVITABLE DEINCERTIDUMBRE EN TODAS LAS MEDIDAS NUCLEARES, Y FRECUENTEMENTEPUEDEN SER LA FUENTE PREDOMINANTE DE ERROR O IMPRESICION* ELTERMINO "ESTADISTICA DE CONTEO" INCLUYE EL ANALISJSESTADISTICO REQUERIDO PARA LOS RESULTADOS DEL PROCESO DEEXPERIMENTOS DE CONTEO NUCLEAR Y PARA HACER LA PREDICCIONACERCA DE LA PRECISION ESPERADA DE LAS CANTIDADES DERIVADAS DEESTOS EXPERIMENTOS+EL UALOR DE LA ESTADISTICA DE CONTEO ESTA DENTRO DE DOSCATEGORIAS GENERALES* LA PRIMERA ES QUE SIRVE PARA CHECAR LAFLUCTUACION NORMAL DE ALGUNA PARTE DEL EQUIPO DE CONTEONUCLEAR+ AS1 UN CONJUNTO DE MEDIDAS SE REGISTRAN EAJOCONDICIONES EN QUE TODOS LOS ASPECTOS DEL EXPERIMENTO SEMANTIENEN CONSTANTES+ PORQUE LA INFLUENCIA DE LASFLUCTUACIONES ESTADISTICAS, DE ESTAS MEDIDAS NO SERAN TODASLAS MISMAS PERO MOSTRARAN ALGUN GRADO DE UARIACION INTERNA* LACANTIDAD DE ESTAS FLUCTUACIONES PUEDEN SER CUANTIFICADAS YCOMPARADAS CON LAS PREDICCIONES DE MODELOS ESTADISTICOS+ SI LACANTIDAD DE LAS FLUCTUACXONES OBSERVADAS NO ES CONSISTENTE CONLAS PREDICCIONES UNO PUEDE CONCLUIR QUE EXISTE ALGUNAANORMALIDAD EN EL SISTEMA DE CONTEOI LA SEGUNDA APLICACION ESGENERALMENTE MAS COSTOSA Y TRATA DE LA SITUACION EN QUENOSTROS SOLO TENEMOS UNA MEDIDA+ PODEMOS USAR AS1 LAESTADISTICA DE CONTEO PARA PREDECIR UNA INCERTIDUMBREESTADISTICA INHERENTE Y AS1 ESTIMA UNA PRECISION QUE DEBERIASER ASOCIADA SIMPLEMENTE CON LA MEDICION+CARACTERIZACION DE DATOSPODEMOS EMPEZAR SUPONIENDO QUE TENEMOS UNA COLECCION DE MMEDIDAS INDEPENDIENTES DE LAS MISMAS CANTIDADES FISICAS:

ADEMAS SUPONEMOS QUE UN VALOR SIMPLE DE X i TIPICO DE ESTECONJUNTO PUEDE TOMAR VALORES ENTEROS SOLAMENTE, ASI QUE LOSDATOS PUEDEN REPRESENTAR, POR EJEMPLO, UN NUMERO DE LECTURASSUCESIVAS DE UN CONTADOR DE RADIACION PARA INTERVALOS DETIEMPO (REPETIDOS) DE IGUAL LONGITUD+ LAS PROPIEDADESELEMENTALES DE ESTE CONJUNTO DE DATOS SON:nSUMAMEDIA EXPERIMENTAL -&= *h6b Ik.2ES FRECUENTEMETE CONVENIENTE REPRESENTAR EL CONJUNTO DEDATOS POR UNA FUNCION DE DISTRIBUCION DE FRECUENCIA F(X)CORRESPONDIENTE+ EL VALOR DE F(X) ES LA FRECUENCIA RELATIVA ENQUE EL NUMERO APARECE EN EL CONJUTO DE DATOS. POR DEFINICIONLA DISTRIBUCION AUTOMATICA NORMALIZADA, ES3z FCK)=t 6.4*=OLA FORMA RELATIVA DE LA FUNCION DE DISTRIBUCION INDICACUANTITATIVAMENTE LA CANTIDAD DE FLUCTUACION INTERNA EN ELCONJUNTO DE DATOS, Y EL ANCHO ES UNA MEDIDA RELATIVA DE LACANTIDAD DE FLUCTUACION O DISPERSION DE LA MEDIA EN UNCONJUNTO DE DATOS DADO+ES POSIBLE CALCULAR LA MEDIA EXPERIMENTAL USANDO LA FUNCIONDE DISTRIBUCION A LOS DATOS, PORQUE LA MEDIA DE CUALQUIERDISTRIBUCION ES SIMPLEMENTE:006.5TAMBIEN ES POSIBLE DERIVAR OTRO PARAMETRO, CONOCIDO COMO LAUARIANZA, QUE SIRVE PARA CUANTIFICAR LA CANTIDAD DEFLUCTUACION INTERNA EN EL CONJUNTO DE DATOS+ EL PRIMER PASOPARA DEFINIR LA DESUIACION DE CUALQUIER DATO COMO LA CANTIDADQUE DIFIERE DEL VALOR MEDIO ES6.6ESTA DEBE TENER UNA CONTRIBUCION IGUAL A LAS DESVIACIONESPOSITIVAS Y NEGATIVAS, AS1 QUE- h2 €

~SI DAMOS EL CUADRADO DE CADA DESVIACION, DEBE RESULTAR UNNUMERO POSITIVOt AHORA PODEMOS INTRODUCIR LA UARIANZA SxK.2COMOnQUE AHORA SIRVE COMO UNA INDICACION DE LA FLUCTUACIONINHERENTE DEL PUNTO ORIGINALt SI EL NUMERO DE PUNTOS, M ESRAZONABLEMENTE GRANDE, LA VARIANZA ES ESCENCIALMENTE EL VALORDEL PROMEDIO DE LA DESVIACION CUADRADA DE CADA PUNTO. PARA SERPRECISO, LA VARIANZA ESTA FUNDAMENTALMENTE DEFINIDA COMO ELVALOR PROMEDIO DE LA DESVIACION DE CADA PUNTO DEL VALOR MEDIOVERDADERO t QUE DEBERIA SER DERIVADO SI UN NUMERO XNFINITODE PUNTOS ES ACUMULADO, rlCOMO CONOCEMOS PARA-UN VALOR INFINITO DE HEDIDCIS,NOSOTROS USAMOS EL VALOR k DERIVADO DEL MISMO CONJUNTO DEDATOS EN QUE SE VAN A CALCULAR LOS VALORES DE LASDESVIACIONES* EL USO DE LA MEDIA EXPERIMENTAL, EN VEZ DE LAMEDIA TEORICA, TIENDE A REDUCIR LA DESVIACION PROMEDIO, POR LO 1TANTO, RESULTA UNA VARIANZA MENOR QUE LA NORMAL, ENIESTADISTICA EL NUMERO DE GRADOS DE LIBERTAD DEL SISTEMA TIENEQUE SER REDUCIDO UNA UNIDAD (M-l.)tSE PUEDE CALCULAR LA UARIANZA DIRECTAMENTE DE LA FUNCION DEDISTRIBUCION F(X)t PORQUE LA ECt(6.9) INDICA QUE SXX2 ESSIMPLEMENTE EL VALOR PROMEDIO DE (X-)1 )XX2 QUE SE PUEDEESCRIBIR EL MISMO PROMEDIO COMO'F=¿SI LA ECUACION (6t10) SE EXPANSIONA, OBTENEMOS6. IIIDE LO DESCRXTO ANTERIORMENTE SE PUEDEN OETENER DOSCONCLUSIONES IMPORTANTEStit CUALQUIER CONJUNTO DE DATOS PUEDE DESCRIBIRSECOMPLETAMENTE POR LA FUNCION DE DISTRIBUCION DE FRECUENCIAF(X) tZt DOS PROPIEDADES DE ESTA FUNCION DE DISTRIEUCION DEFRECUENCIA SON DE PARTICULAR INTERES: LA MEDIA EXPERIMENTAL YLA VARIANZAI.LA MEDIA EXPERIMENTAL ESTA DADA POR LA ECt(6t5), Y ES ELVALOR APROXIMADO DEL CENTRO DE LA DISTRIEUCIONt LA VARIANZAESTA DADA POR LA ECt(6tl.O) Y E8 LA MEDIAD DEL ANCHO DE LA

DISTRIBUCION, O LA CANTIDAD DE FLUCTUACION INTERNA DELCONJUNTO DE DATOS,M O D E L O S E S T A D I S T I C O S tBAJO CIERTAS CONDICIONES, SE PUEDE PREDECIR LA FUNCION DEDISTRIBUCION QUE DESCRIBIRA EL RESULTADO DE MUCHASREPETICIONES DE UNA MEDIDA DADA. NOSOTROS DEFINIREMOS UNAMEDIDA COMO EL NUMERO DE CUENTAS DE UN SUCESO RESULTANTE DE UNNUMERO DE EVENTOS DADOS+ EN CADA EVENTO SUPONEMOS QUE ES UNPROCESO BINARLO, EN QUE SOLO DOS RESULTADOS SON POSIBLES: ELEVENTO OCURRE O NO OCURREt PARA TODO LO QUE SIGUE SUPONDREMOSQUE LA PROBABILIDAD DE QUE OCURRA ES CONSTANTE EN TODOS LOSEVENTOS +EN EL CASO DE EVENTOS DE CONTEO DE RADIACION NUCLEAR, UNEXPERIMENTO CONSISTE EN OBSERVAR UN NUCLEO RADIACTIVO DADO PORUN PERIODO DE TIEMPO t, EL NUMERO DE EVENTOS ES EQUIVALENTE AEL NUMERO DE NUCLEOS EN OBSERVACION, Y LA MEDIDA CONSISTE ENCONTAR AQUELLO NUCLEOS QUE SE DESINTEGRANt SI IDENTIFICAMOS LAPROBABILIDAD DE QUE CUALQUIER EVENTO OCURRA ES Pt EN EL CASODE DESINTEGRACION RADIACTIVA, LA PROBABILIDAD ES IGUAL A(l-EXP(- ht), DONDE '& ES LA CONSTANTE DE DESINTECRACION DELNUCLEO RADIACTIVO+EXISTEN ALGUNOS MODELOS ESTADISTICOS COMO, LA DISTRIBUCIONBINOMIAL, LA DISTRIEXJCION DE POISSON Y LA DISTRIBUCION NORMALtEN ESTE CODIGO SE SUPONE QUE LOS EVENTOS ESTAN DISTRIBUIDOS DEACUERDO A LA DISTRIBUCION NORMALtD I S T K I B U C I O N N O R M A LtMUCHOS PROCESOS BINARIOS PUEDEN SER REPRECENTADOS POR UNASPROBABILIDADES PEQUENAS DE QUE OCURRANt AQUI SE INCLUYENMUCHOS EXPERIMENTOS DE CONTEO NUCLEAR EN QUE UN NUMERO GRANDEDE NUCLEOS HACE QUE EL NUMERO DE EVENTOS, SEAN CONSIDERADOSCOMO UNA FEQUENA FRRCCION RELATIVAMENTE PEQUENA DE LAS CUENTASREGISTRADAS, SIMILARMENTE EN EXPERIMENTOS NUCLEARES, MUCHASPARTICULAS CON UNA GRAN ACELEKACION GOLPEA A UN NUCLEO BLANCOPARA QUE CADA PRODUCTO DE LA REACCION SEA REGISTRADO, BAJO 1ESTAS CONDICIONES LA APROXIHACION DE QUE P .

Z+ LA DISTRIBUCION ESTA REPRESENTADA POR UN SIMPLEPARAMETRO Cjt j QUE ESTA DADO POR EL PRODUCTO MP+3+ EL VALOR PREDICHO POR LA VARIANZA d XXZ SE DEFINE PORSE PUEDEN HACER DOS OBSERVACIONES IMPORTANTES CON RESPECTOA LA DISTRIBUCION NORMAL:I+ LA DISTRIBUCION ES SIMETRICA CERCA DEL VALOR MEDIO 3*POR LO TANTO P(X) DEPENDE SOLO DEL VALOR ABSOLUTO DE LADESVIACION DE CUALQUIER VALOR X DEL VALOR MEDIO, DEFINIDO COMO6% \x-W11, DEBIDO A QUE EL VALOR MEDIO DE ES GRANDE, LOSVALORES DE F'(X) PARA VALORES ADYACENTES DE X NO SON MUYDIFERENTES UNO DEL OTRO. EN OTRAS PALABRAS, LA DISTRIBUCIONVARIA SUAVEMENTE+PARA CHECAR QUE LAS FLUCTUACIONES OBSERVADAS CON EL SITEMADE CONTEO SON CONSISTENTES CON LAS FLUCTUACIONES ESPERADAS DELMODELO ESTADISTICOj SE INTRODUCE LA PRUEBA DE " CHI-CUADRADA"+CHI-CUADRADA ES SIMPLEMENTE OTRO PARAMETRO DE LADXSTRIBUCION DE LOS DATOS EXPERIMENTALES Y SE DEFINE COMODONDE LA SUMATORIA ESTA DADA SOBRE CADA PUNTO XNDIVIDUAL X i+CHI-CUADRADA ESTA RELACIONADA CON LA VARIACION MOSTRADA PORx-Xc6.15AHORA SI LA FLUCTUACION PRESENTE EN LOS DATOS ES ANGOSTA,MODELADA POR LA DISTRIBUCION NORMAL, AS1 S2ZG1+ PERO NOSOTROSSABEMOS QUE PARA LA DISTRIBUCION NORMAL @?= ADEMAS, EES IGUAL A POR LO TANTO, EL GRADO DE QUE LA RELACION S*xZ/)('SE DESVIE DE LA UNDAD ES MEDIDO DIRECTAMENTE LA MAGNITUD PARAQUE LA VARIANZA MOSTRADA OBSERVADA, SEA DIFERENTE DE LAVARIANZA PREDICHA+I6+3+ DIAGRAMA DE FLUJO+EL DIAGRAMA DE FLUJO QUE SE PRESENTA A CONTINUACION ESTALIMITADO A LAS OPERACIONES PRINCIPALES, PERO SI SE REQUIEREMAYORES DETALLES EL LISTADO SE ENCUENTRA EN EL APENDICE E+

I og

6t4t USO DEL PROGRAMA ( ENTRADAS Y SALIDAS I tCUANDO EL CODIGO CORRE IMPRIME UNA COLUMNA EN LA CUAL ESTANLOS PRIMEROS CINCO CALCULOS QUE SE DESCRIBEN EN LAS "SALIDAS",EL USUARIO DEBE ANOTAR LOS DATOS EN FORMATO LIBRE (SOLO UNOPOR LINEA), CUANDO TODOS LOS PUNTOS YA SE ANOTARON, EL USUARIODEBERA ANOTAR 99t99 PARA SENALAR QUE YA SE ANOTARON TODOS LOSDATOS tEL PROGRAMA ACEPTA HASTA 200 PUNTOS, YERO ESTE NUMERO SEPUEDE INCREMENTAR CAMBIANOOLO EN LAS DECLARACIONES DEDIMENSIONtSI SUPONEMOS QUE SE TIENEN M MEDIDAS DE UNA FUENTERADIACTIVA, TENIENDO CADA MEDIDA Xi CUENTAS CON i=1j2j3,tttMtEL USUARIO AL ANOTAR CADA MEDIOA, EL CODIGO CALCULA E IMPRIMEINMEDIATAMENTE LOS SIGUIENTES CINCO RESULTADOSt1 t M: NUMERO DE MEDIDAS ANOTADAS*Zt X(M):LA ULTIMA MEDIDA ANOTADAt13t XAVR(X): EL PROMEDIO DE LOS PUNTOS (CON LA EC 6t2)t4* DEL2(X-Xi)XX2: PARA LA ULTIMA MEDIDA ANOTADAt3t SIGMA: LA DESVIACION ESTANDAR (SXXZ) DE LAS M MEDIDAS 1ANOTADAS (CON LA EC+ 6t8)tPARA M=lj SIGMA SE CALCULA CON SQRT(Xi)tYA QUE EL USUARIO ANOTO TODAS LAS MEDIDAS A 4NALISARj SE 1IMPRIMEN LA MEDIA Y LA DESVIACION ESTANDARtLOS CALCULOS SIGUIENTES SE HACEN CON EL 75Xj 90Xj 99% DELINTERVALO DE CONFIANZAt UN INTERVALO DE CONFIANZA DEL 75%REPRECENTA EL INTERVALO EN QUE LA MEDIA VERDADERA SEENCUENTRA, ESTOS LIMITES SE CALCULAN POR MEDIO DE LA ECUACIONIX(+/-) =: x +/-tsDONDE X+ Y X- SON LOS LIMITES SUPERIOR E INFERIOR DELINTERVALO, t SE OBTIENE DE LA DISTHIBUCION-t PARA M(21 Y SESUPONE DISTRIBUCION NORMAL PARA CASOS CON MAS DE 120 MEDIDAS*LA TABLA DE PORCENTAJES DE LAS MEDIDAS DE LA DISTRUBUCION-tjSE MUESTRAN EN LA TABLA 6tIj DONDE SON LOS GRADOS DELIBERTAD QUE SE DEFINEN COMO M-1 PARA UN CONJUNTO DE MMEOIDASt NOTESE QUE SI CORRESPONDE A LA DISTRIBUCIONNORMAL t

EL TERCER CONJUNTO DE RESULTADOS SON PARA LA PRUEBA DE"CHI-CUADRADA" (XXX21, ESTA PRUEBA SE REALIZA PARA VER SI TODO .EL CONJUNTO DE MEDIDAS REPRECENTA UN "CONJUNTO TIPICO",LLAMADO AS1 SI MUCHOS CONJUNTOS DE MEDIDAS SE DAN BAJO LASMISMAS CONDICIONES, AS1 EL PRIMER CONJUNTO CONCUERDA CON LOSOTROS+ PARA ESTE PROPOCITO, SE CALCULA Xmm.! CON LA EC+(á+15)+CABE NOTAR QUE XXX2 DA LA DISPERSION EN EL CONJUNTO DEMEDIDAS+ ENTONCES SE DEBERA GUARDAR EL CONJUNTO QUE .TENGA UN"PROMEDIO" DE DISPERSION; Y RECHAZAR AQUELLOS QUE TIENENDISPERSIONES MUY PEQUENAS O MUY GRANDES+EN LA TABLA 6+II SE DA EL PORCENTAJE DE PUNTOS DE XXm2+ LAINTERPRETACION ES LA SIGUIENTE: DADO XXX2 PARA UN CONJUNTO DEMEDIDAS 9 QUE ES LA PROBABILIDAD DE QUE OTRO CONJUNTO DEMEDIDAS SE OBTENGAN BAJO LAS MISMAS CONDICIONES, DEEERA TENERXXX2 (DE DISPERSION) MAYOR QUE EL PRIMER CONJUNTO, CON LAPROBABILIDAD DE P/lOO+COMO UN CRITERIO DE ACEPTACION DEL CONJUNTO DE MEDIDAS, ENSTAT SE SUPONE QUE XxrZ DEE€ ESTAR DENTRO DEL INTERVALO DE1-99X; CUALQUIER CONJUNTO FUERA DE ESTE INTERVALO SERARECHAZADO+ EN LA TABLA 6+11, SE REFIERE A LOS GRADOS DELIBERTAD M-l+SI EL CONJUNTO DE DATOS CONTIENE MAS DE 101 MEDIDAS, SESUPONE QUE Y= SQRT(ZXXZ2) SE APROXIMA A UNA DISTRIBUCIONNORMAL CON MEDIA Y= SQRT(3, -1) Y S=l, AS1Y USANDO EL PORCENTAJE PARA UNA DISTRIEUCION NORMAL'(t=2+33PARA EL 99x1 COMO EL CRITERIO PARA RECHAZAR MEDIDAS+PARA FINALIZAR, STAT EXAMINA CADA MEDIDA DEL CONJUNTO YDETERMINA CUALES MEDIDAS ESTAN "MUY DESVIADAS" DE LA MEDIDA+SI SE DETECTAN ESTOS PUNTOS SE IMPRIMEN+ EL USUARIO DEBERAELIMINAR ESTOS PUNTOS Y VOLVER A CORRER EL CODIGO,EL CRITERIO PARA EL RECHAZO DE LOS PUNTOS ES EL "CRITERIOMODIFICADO DE CHAUUENET"+ PARA CADA PUNTO EN EL CONJUNTOSE CALCULA, Y SE COMPARA CON LA t CORRESPONDIENTE ALPORCENTAJE DE MEDIDAS CALCULADO POR (1/2M)+ (PARA M100)+

Pwecnta~e poinfs cf tho i-dístr;J&nP U 10 5 2 1 GZ 0.1TABLA - 6.1,. IIcLr- t 2.41 6-31 12.71 31.82 63-66 318.3 636-62 1-60 292 4.30 696 9.92 22-33 31.603 1.42 235 3.18 454 J-84 10-3 12.924 PW 213 2.78 3.75 4.60 7.17 8.615 1-30 24?3-36 4-03 5-69 6.87- ci6 1.27 1-94 1.45 3.14 3-71 5-21 5967 1.25 1.89 236 3-00 3.50 4.1 5.41- 8 . -1.24 1-86 . . 1-31 2.90 3.36 4-50 5-049 1.23 1.83 226 2-62 3.2S 4-30 478- 13 1.12 - 141 2-23 . -2.76 . '3.17 4-14 .+59 .12 1-21 1-78 2-18 2-68 3.05 3-93 4-3215 1-20 1.75 2.13 1-60 2-95 3-73 4-0721) 1.18 1-72 2-09 2-53 2,85 335 3-8524 1.18 1-71 2-06 2-49 280 34 3-7530 1-17 1-70 2-04 2-46 2.75 3.39 3-6540 1-17 I -68 2-02 2.42 270 3-31 3.5560 1-16 1-67 200 1-39 2.66 3.23 3-46120 1-16 1-66 1.98 2-36 2-62 3-16 3.37m 1.1s I -64 1.96 233 258 3-09 3-29P 99.5 99 97.5 9s10 . 5 25 1 0.5 0.1r-1 00'3Y3 0-0'157 0.0'932 O.Cí3932 001üG WEOl 00506 Olüj3 00717 0-115 0216 12-324 0.207 0-237 0-484 0.7112.7 14-616-257.783-845.997.8;s-495.02 6-637.38 9-2f9.35 11.3411-14 i3.37-88io-úo12.84I4-8Ó1g-8313-8116-2118.47TABLA 6 .SI5 0-412 0.554 0-831 1-136 0676 OS72 1-24 i-6:7 0.989 1-24 1.69 >!7S 1.34 1-65 1-18 2.739 1-73 2-09 270 3.3310 2-16 2-56 3-25 3-9411 2-60 305 3.82 4-S712 3-07 3-57 4-40 5-2313 3.s7 4-11 5.01 5.5914 497 4-66 5.63 6-5715 4-60 513 6-26 7-26I6 5-14 5-81 6-91 7-9617 5-70 6-41 7-56 8.6718 6-26 7-01 5-23 93919 6.84 ?--63 8.91 lei220 7.43 8-26 9-59 10-8521 8-03 8-90 10.28 11-5922 8-64 9-54 10.98 12.3423 9.26 1020 11-69 13.09U 9-89 10.S6 1240 13.8525 1052 ri-s- 13-12 1.1-6126 11.16 I220 13.04 15.327 11-81 E-SS 14-57 i6.JS23 12-46 13-56 15-31 16.9329 13.12 14-26 16.05 17.7133 13.79 14-35 16-79 18.49do 20.71 E.:6 24.13 26.51S!) 27-99 29-31 32.36 34-1669 35-53 37-48 @4S 43-1970 43-25 4á-43 48.76 51.7480 51.17 53-54 57-15 tii.3990 59.20 GI 15 65-65 69.13100 67.33 10-96 74.22 77.939-2410.6412.0213-361465114712.914-07IS-SI1692i2-s3 15.391445 16-8116-01 18413-53 20.0919-02 21-671ó-75:s.SS70.2821.9523-5920.5222.4624.3226.1227.8815.99 18.31 20.48 2321 25.19 3-5917-28 19-65 21.92 24.73 26-76 21.2618.55 21-03 23.31 25.22 28.30 32.9119.81 2-36 24.74 27.69 29-32 34.5321-06 23.1% 26-12 29.14 31.32 36.122-31 3.00 27-49 3658 32-80 37.7023-54 26.30 3.65 32-OO 34-27 39-2521-77 27-39 30.19 33-41 35-72 40.7925-99 73.87 31.53 3.81 37-16 42-3127~030-1-132.85 x.19 38.58 43.8228-41 31.41 34-17 37.57 40-00 45.3129.62 32.67 3S.48 3-93 41-40 468030.S1 33-92 36-7s 4-29 42.30 48.273i.01 35-17 3s.o~ 41-61 u-ia 49.7333.20 36.42 39-36 42-98 45-56 51.1834-35 31-65 40.65 4.31 46-93 52-6235456 38.59 41.92 45-64 49-29 54-0536-74 40-11 43.19 d5.96 29-63 55-4837-92 41-33 4.46 4.28 50-99 SU.E?39-09 42-5ó 45-72 49.59 52-34 58.3040.26 93-77 46.94 53-89 53.67 59-7051.81 55.76 59.34 63-63 65.77 73463-17 67.50 71.42 -6.15 79.49 S1.567440 79.0s 83-3C IIa 33 91-93 39-61

TABA 6.111LINEAR INTERPOLATION OF t VALUES FOR DATA POINT REJECTIONDegrees of Freedomti v g 0 2.231 ’Pc10 < -v < 25 ____ - _-2.23 + (v-10) * 0.017c__--.--25-

TABLA 6.IVCumulativc Normal Distribution'..1 A.O1 .O2 .O3-.o4-5160s557SS48.633 I.67W-.os-SI995596S987.5368-673632395636.6026.6jo6.6772-.O7S279S675.w.w3.6808-.O8S31957 14.6 I03.&SO.6844-5359s7s3.6141.65 17.6879.6915 .6?50 .69SS .701Y.7257 -7291 .7324 I .7357.7SSO .76l1 .7612 .7673-7881 .7910 -7939 I .7Y67.?OS4.7389.77&.7995.I3264.TO88.742Z.:734.SO23.a289.7 I23.I454.7764.SO5 I-8315.7 I57.748h.7794.a078.!?;a-7 193.75 I7.7S?3.SI06.S365.7224.7549.7352.813333891.49192 9107 .9222 I .9236.8538.373.w25.!?¡Y39.92s 1.8531.87493944.9115.Y265.Y5543770.iiObl.9131S2?935773790.XOXO.9 147.9292.8599.8810.s9w.9 I429 3%.862 I.8830.yo15.9:77.9319' I1.5', 1.61.1131.9-9332 9345 9357 .9370.942 9463 .Y474 .Y4849SW -9% .9573 .4582.9MI .9hM .VhM9713 .Y719 .Y726 .9732.9383.9495.9591.967 I.9738.9394.9505.9s99,9078s744.94%.95 I5.9MX!JiNt.9750.MI8.!IS25.9ó I ó.969 3-9756.9a9.9535.9625.96W.9761.se? I9545.96330706.97672.0212.22.32.4,9772 .9778 .9783 .9783.9821 .9926 .990 .9334-9861 .YS64 .9865 .9871.9593 .9896 .9998 1 .99@1N I Y .y93 .Y922 .!XI25.9793.9S38,9875.9m.9927.9798.9842.9878.9906.!I929.9803.9E*.98S I9909.993 I.98G?.9550.963.99l I.99?2.98 12.933.9987.99 I3.99x.98173857.989!1.99 1 t.-933625262.72.82.99938 .YYSO .YY~ s.13.9953 .9955 .9956 99579965 .9Y66 .Y967 .9Y6Y974 ,9975 .9976 99779981 .9992 .9982 9983.99:5.995c.9969.99779984-9946.%o03970.9978.9384.9C48.Y96 I.Q97 I.9979.9985.99:9.9961.997?99799985.99519.263,9973.9930.9966-9052.99&9974.998 I9863.03.13.2333.49997 .99Y7 9987 .W88.9990 .9991 .99YI '.9991.9993 .!%Y3 .FV% .Y993.9995 .9995 .9995 .9996.W7 .9997 I .W97 .*Y7II9988.9991.9994.99963997-.9989.9991.9991.9996.9997-.9999.9992.Y934.99%.93079989.999?.Y995.9996.9Y9779990.99939995.99969997--.99909993.99?5.9937.9w87,Ir 3

6+5+ SIMULACION6+5+1 DATOS PARA LA SIMULACIONPARA OBTENER LOS DATOS PARA SIMILAR ESTE CODIGO FUENECESARIO HACER ALGUNOS EXPERIMENTOS DE CONTEO NUCLEAR, QUE SEREALIZARON EN EL LABORATORIO DE FISICA NUCLEAR DE LA UAM-I+ YEL SISTEMA UTILIZADO SE MUESTRA EN LA FIGURA 6+1, Y LOSRESULTADOS SE PRESENTAN EN LA SIGUIENTE SECCION+615.3, RESULTADOSLOS RESULTADOS QUE SE OBTUVIERON TANTO DE LOS EXPERIMENTOSREALIZADOS COMO LOS OBSERVADOS AL CORRER EL CODIGO SEPRESENTAN A CONTINUACIONII I y

FIGURA 6.1DIAGEIAX-9 A BLOQUSS DEL SISTEMA UTILIZADOPARA CONTZO NUCLSAB..*....I ’*: 2 ‘4 5T .i 4b 1.7 8i4c3 691) Fuentes ( Cs-137 de 8 microcuries y eo-60 de 34microcuries).2) Detector (HARSKAÍ’J Na(T1) ).3) Fuente de alto voltaje t ORTEC 456 ).4) Preamplificador ( HARSHX.7 NB-11).5) Amplificador ( HARSHAJ NA-17 ) O6) Osciloscopio7) Discriminador ( ORTEC 420A ).8) Escalador ( HARSHAfl NS-12 ).9) Reloj ( HARSHA‘J NT-29)...... &. - .< ”.

12T.>il5b74'21 II111;"131 '1113s171819¿C212223242§110. 110.00r351 O¿ . si)11;.Ql1fib.33102 . 50116. 105.2092. 103.00115. 104.71107 105.00120 . 1 0 ~ ~ 6 7111. 1 O 7 , 1 :I99 . 106.368s. 104.6712L).97.106.15lOS,l41 i'i, 1 c! 5 di?1U7. 195.01)102. lUS.t35122. 106.56101. 106.26105. 1 ob ¿d97. 105.5i?111. 105-77117. 10b.26Y7. 5.86103. 5.96.o0'36.2558.78132.25116.64121.00105.804.0u177.7615.2154.22348.44318.49172.7364.h41.2713.30238.5327.701.44162.R927.32115.3378.774.1610048810,6071 O . 0 1711.21011.43211.55911.48510.6b411.1S8100b0910,35611,49212.24012,35312.1731 1 -76411,42911.73811.47911.17611.32511.11411.10811.0¿710,803* ORARA P LOS DATOSESTOS DATOS ES 26.434

TúDUS LUC PüPiTUS SATISFACEN LA MODIFICACIUN DE CHAVEVET1i??a5f>7891 '1I tl?131415161715192821¿?252 4252627í? 3293031JZ333 435Jh37383940fJ I4 2434445'4 ó1 1 Y6 . 0 O1176.501191.671206.001215.20l¿lb,b31221,uo12228 . 501230.001234.001235 .'9 11239.921240.921240.791244 . 071246.Úu1247.62lc'4b.5012463 . o51246.351246.671í47.4512¿46.171246.381 245 . o 01240.óQ1 ZJb .1245.451¿46.06124'5.151 ¿ 4 b . 2 61247.391249.36124d 11124d . 111247.2812 4 7 . II 01,248 . 341249 . 621249.231248.1249.12441 . 46. O0930.259¿!0,11I649.0013S4.2406.70625.002756.25144.001296.003bii.461943.33145.853.192109.88841 .OG651.27506.2364.8531 09240.11273.75793.77(31 -391918.441.57525.60625.00.2238.27362.197304.56358.63849.821349.493037.5d2hbb.582034.75.o1978.56.?OC). u4LioQ.042743.76243.50142.91384.166-4733.85343.13440.25343.61243.00038.67536.98740.28137.94737.94736 . 55337.51536.10035.40255.09354.50533,59032.72131.92531.37431 -26230.59131.30631 -21 131.01731.78q31.44231.06330.77230.73630 . 95s31 -96932.72833.15932.70732,68132.425-32.20532.86932.57432.24b32.01731 -662

120Q.99- 1295.971191.i8- 1314.781154.0ó- 1351.407 S% 1.16690% 1.67699% 2.682DbTERES AhJTAG UTRiIS üATGS (C/IJ)?UN IvrAXIMr) DE ZIIEDEFJ ANOTARSEPARA TEHI4INAHENTRADA, ANOTAbrx I r-1 1 XAVR DELL) SIGMA~L.oIoIwIa.mooII ~ ~ ~ ~ ~ o I I . o I I ~ I I o * I o ~ ~ . ~ ~ ~ . ~ . ~ a * ~ . o .1 8 1 I 81 ,u0 , o 0 9.000? 71. 76.00 25.00 7.0713 Th. 75.00 . o o 5,0004 77 14.06 4.7875 77 .o4 4.1476 76.83 66.69 5.4567 52. 80.71 127.37 7.041J Y .-I. 82 3íf 135.14 8.035J dj. 82.b7 53-44 7.5h610 93. 83-70 86.49 7.84b11 84.09 15.28 7.55612 8Y.SB 29.34 7.40313 17, G3.63 114.33 7.78211 8

7';.!4(j- Q3.tir) 75% 1.164f 1 .oh - 77.92 90% 1.7¿66?.99- 106.61 99% 2.870* PRUtBA DE Ct4I-CUADRADQ PARA LOS DATOS!x PF(Ut5A DE LA MODIFTCALIOM DE GHAiIVEfJET PAHA LOS DATOSTODOS LljS PUUTOS SATISFACEN LA gvtODIFXCACION DE CHAVEVET914. 914.00004. 909.0335: . 9¿'s . 3333 1. 927.259ü7.913.937 . 929.43Si2'2. 926.50'277 . 4Z3.890 0 *3 . 93J.30914. 929.73970. 933.0898s. 937 .o7n. 979..o025.ni)1024.004'j05b2149.64230.0357.3342.251858.57542.89247 -351362.842296 . b29 3'1 . h 130.2227.07128.16023.429%2.12921.14219,58618.323~23.56123.57423,05724 8b927 . 82520.144

I_____.,J.."' .". 0 ..'.v.... o... .r....O....l....O....r....O....r....O... '.. ..o.v15 '201 I 340.d7 405.35 27.tJY71 (7 -2 I! 1 YJ(1.30 1394.89 2G.54417 o42 . 438.59 11.64 27 . 652189480 939.11 79.01 c'h.9181s 972. 940.84 970.81 Zf . 22620 q54. 'i 4 1 . 50 1Sb.25 ¿6 -663dl 907 933,Yb 1079.59 ¿7 .OS6¿2 "37 . 539.73 7.4G 26 4 1 1I23 97'1. 941 ,43 1411.15 27 . 07224 963. 942.33 427.11 i 26.84025 907 . c)40.92 1150.57 27.2092 (> 013. 939.85 720.7¿ 27.¿1ó27 0 ct 7 9 'I(1. 1 i 47.46 2ti.723CY 9b(i . 9';d.Hh 1146.31 27,05029 QZ(i. 938 . Y0 109.89 25.63930 379. 939.83 1534.03 27.20131 404, 938.68 1202.52 27.50752 843, Y37.25 ' 1958,Ob 28.239'33 3?¡1 936 , 8') 1bS.08 27.8i3934 ad7 930.52 2'369.41 2t3.7793s 973. 939.31 1134.72 28 95236 912. 938 . 56 705.20 ZN . 89637 933. 938.43 19.65 28.50238 917. 937 . b8 766.42 28.4903íJ '>4h. 537.90 bT.bb ¿#. 1444c ql?. 937.25 637.56 l28.08141 9270 337.00 100.00 27.77442 suo. 936.12 1304.58 28,02143 896. 935.1'1 1535.5-i 28.35344 946, 935.30 22.13 28.0310 5 be!'). 93t2.93 10502.1s 31 .92145 Q4f7 . 933 . 0'2 47.79 51 . 53247 931. 933.04 4.17 31 -25848 375. 931.83 3230.03 32,020~ ~ ~ ~ o i ~ o ~ w ~ c ~ ~ ~ ~ ~ ~ ~ ~ ~ rcPARA 47 GRATJOS DE LIBERTADILOS LIMITES DE CHI-CUADRARACORR NU E S AL INTERVALO ESTAPJ ENTRE 27.44 Y 72-41* PkI'F'IA I,E LA i~~i'!nTFICkCTdr' C'iArlVEiiiET PARA LüS DAT9S Id0

__I__._.__I ~.. . . . o . . ......O....,.... O.........O....,....O.... ..... 0 ' . . . v . . . . 0 . . ' , v . . . ' o .1 lllb9. 11169.002 11161, 11165.007 11 tlA-5. 111271.00J 11 li8J7 . 11305.(105 ilQ77. 114 39 . 406 1174s. 11491.007 11928. 11553.43R 117(17. 11572.63c3 i 1 ¿5,204;) 110lP. 11554.Y543 11981 ilSb4.Hh44 11Qb7. 11574.00.o016.00105.6835 . 55744q44.01 183.6411 OU04.02 164.6412899 13.56 332.h3666564.13 323,25414U304.50 338.17018056.66 317.7S71 O 3 1 ¿I 4 . 13 30415.4632 0 7 16308 . 521147b1.3b 295.447200.68 281.73349011.21 277,82049094-28 274.4331 ?d¿;3? 42 2n2.4735 1075 . 0 1 279,4751714.47 270.dll5207.94 263.34212542R.GS 269.3101159b.52 263.942137147 -47 270.891121h10.05 275.5981753115.75 284.3117825.00 278,6991767'33.31 6,70049403.75 264,542156434.?6 2119 -3971ISST2.61 242.48720'527 1 . 13 3 í) U . 1 441149.24 294.9935567.15 290.3653115.04 285.8259832Li.34 2b6.928808,%2 2d?.!59i39453.41 280 54411778.21 277 . 134173843 -66; 282.19347122.92 280,693174580.7a 2H5.49714 3 fj 7 Q. 59 283.1535H O ó 4 , 23 Y 265.2b8¿qY79d.31 294.548173171.75 298.105154449 . 36 300.869

75%90%99%1.1691.679r'. 694CHI CUAUgADA PAíIA ESTOS DATOS ES 336.309PARA 45 Gi7ADiJC DE LIBERTADtL HI-CUADRADACOHRESPUIJD T t i4TE 3 AL INTERVALO -42 Y 67.431 61.952234 5.06 48,1975 3341:. 722s.ou 614.246h 677'3. 5017.39 6h.3647 386%. 269.90 61.0133 3892. 1444.00 58.5379 38990 1600.00 Sb.774In 3:j7?. 136.89 53.58411 3"1';. 536s 55 2751 -49 53.A1917 3QU6. 3YbCi . 9% 1375.19 52.62613 3867 3868.77 3-13 50.38914 3845. 3867 O7 487.14 48.82715 3920 3670,óO 2440.35 43.19516 3796) . 38b5,9:, 4841.2b 50.97617 3273. 39b6.13 14.6¿ 49.2701,3 5h27, 3bi4.00 1369 . O0 48.b8319 5449. 3868 . 47 646405i) 51.17220 3813, 3865.9 2299.22 51.07021 5894. 36b7-3 765.44 50.179¿ 2 Jlii37. 3 d 7 3 . 7f 1 HZdb . 33 s7.535

.i....O....r....O....,......r....,....~....~....~....r....~.........~....i....~.23 3951. 31377 . 13 5456.73 58.47323 3892. 3877.75 203.06 57.26829 3863. 3883.79 -63 55.86730 3851. 1004.9035 3712. 3875.34 26680.87 59.93540 3801, 3876.02 5628.74 58.267a1 3800. 3674.17 5501.33 58.74642 3923. 3875.48 2864.02 58.63946 3867, 3875.6% 128.77 56.31547 3957 . 3877.38 6338.91 56 . 9504% 3850. 3876.81 718.91 56.479INTERVALO CONFIDENC PROB.ClUE MEDIA VERDAD TT3784.2b- 3970.59 90% 1.674 I99% 2.678ESTUS DATOSESTAN DENTRO DEL RANGO DE 1. - 9Q. x* PRUEBA i)E LA i-~OEjI~ICACIOfJ DE CHAUVENET PARA LOSDATOS

.-... ..O....r....O....r....O....r....O....v....O... .r....O....r....~....r....O.__1"; x ( '41 XAVR DEL¿ SIGMA9 9 1 1 0 ~ ~ ~ ~ ~ ~ 9 " ~ ~ ~ 0 1 ~ L ~ ~ ~1 7707, 7303.00 .o0 85.4582 72jg. 72hb.50 1190.25 48,790I, 7 '43 , 7332.09 1b129.02 115,2694 7717. 7328.25 126.56 94.4155 7413. 7345.20 4596.83 90.1S3ii 7400. 7354.33 85. 83,6567 7 54s. 7 35 3 'J O 64, 76.4485 7213* 7335.50 15066.29 bó.36d3 7 :I 9 (1 7353-67 ¿1121.70 97 . 45410 7270, 7345,30 5670.06 95,61411 7352, 37.10 90.73012 7?7?, 3761 -74 88.63813 724 ?. 7332.85 8 7 2 24 69,054I?) 7 4 113 7337,db 4243.55 87.59115 7307. 7341.13 ¿103.71 85.?5316 7?37. 7334.63 9530.65 86.47117 72b6. 2094.41 s4.5511 (9 7270, 2578.34 83.00014 7 ?¿,) .) 73?3,05 10619.8d d4.4342 17 7-70, 752ij.70 1993.11 t12.853c' 1 7,345 73 1 7 . 1 d S157.80 82 . 42722 7333, 151.74 80.48723 7462; 64.19824 7441. 85.743¿'i 7753- d3.947;rh 7L:Ofi. 7332.12 5758.51 83,694¿72829303152 7 I' 33 7 3 33 . 5 3 101 ab. 54 42.14933 7 32h 7353.12 172,ib 90.7¿934 7209. 7337.65 ¿3b6 58 69,756-3 s 7743. 7337,biO 27 . 04 80.43136 7378. 7338.92 1527,48 ~37.41557 7x50. 7s.zi5,r>;J 75-27 84.2053 I! 7iiti1 7342-42 19204,lb 88.11139 7 j:7

IiiiTtRVALO COidF IClLivC PPOB.flUE E;EDIA VkRDAD TTts Eri ESTE INTERVALO* PHUERA DE CMI-CUADQ&IlA PARA LOS DATOS7 5% l.lbb99% 1.67699% 2,682EST~JS DATUS ESTAN DEiJTRU DEL RANGIl DE 1 , = 93,ESTc) ES ACEPTABLEXXDEL.I -o-"o--.--.o-1 i 3 8 . 836.03 .no 2f3.948P a7 4. lr58,SG 420.2s 28.9913 891, 862 . b7 69-45 21,7334 690, 869, 420,25 22,3985 804, 856, 2745.76 35,1335 dd'i. bbl . e> 9 3 . 4 s 33.963-7I Ro7 . db2

75% 1.17290% 1,7029QX 2.756I* PRUEBA uE CHI-CUAURADA PARA LUS DATOSESTOS 3ATílf. C5Tk;J i)E'JTY(! DEL hA~JfJ0 DE 1. - 93. XESTO ES ACFYTAULEI)€VEVET

6.6 CONCLUSIONESEN ESTE CODIGO DE ESTADISTICA DE CONTEO, ES IMPORTANTECUANDO SE REALIZAN MEDIDAS NUCLEARES, DEBIDO A QUE SE PUEDEDAR LA PRECISION ESPERADA DE LOS DATOS OBTENIDOS, Y SIRVETAMBIEN PARA CHECAR SI EL SISTEMA DE CONTEO ESTA OPERANDOCORRECTAMENTE,EN EL CASO DE LAS MEDIDAS QUE SE ANALIZARON CON ESTE CODIGOY QUE SE OE3TUUIERON EXPERIMENTALMENTE, MUESTRAN QUE ESTANDENTRO DEL INTERVALO ESPERADO DE ERROR, PERO EN UN CONJUNTO DEMEDIDAS, SE OBSERVA QUE LAS NO SON ACEPTABLES, ESTO SE DEBE AQUE AL TOMAR ESTAS MEDIDAD SE DEJO QUE SE SATURARA ELDETECTOR, POR LO CUAL LAS DISPERSIONES ENCONTRADAS FUERON MUYGRANDES +

i?~:cli'i' E:S IJN CCII:ICIC;O E?~CRITCI EN FORTRAN PARA EL uso.Lbl.rFr* .,. ..I. IVIL 1 XU0 DIE: CClMF'IJ'r'ADClRAC QUE CALCULA EL. NUMERO DE: ATOMOS YI...A i^JC'T:X:VXT.)AFS DE DIFERENTES ELEMENTOS DE UNA CADENA QUE SEl~~I~~:SXN'T'I~~Ri? RADIACTI'JAMENTE 9 COMO UNA FUNCION DEL TIEMPO +TE:ORIA EN QUE SE BASAC:I.~ANOO I...OS Amims TIENEN MAS NEUTRONES QUE PROTONES ENNIJCI..EO.j CC)N Z MhYUR QUE 20, ESTOS NEUTRONES EXTRA SONNEt:CE:CARIO!j PARA LA ESTABILIDAD DE LOS NIJCLEOS PESADOS + ELE:XCEC13 DE NEUTRONES AClUA COMO UNA "GOMA" NUCLEAR, MANTENIENDOUNIDO Al, NUCLEO COMPENSANDO LAS FUENZAS ELECTRICAS REPULSIVASENTRE LOS PROTONES CARGADOS POSITIVAMENTE+ GENERALMENTE HAYCACOS EN QUE ALGUNOS NUCLEOS CON EL MISMO NUMERO ATOMIC0 PERODIFERENTE NUMERO DE NEUTRONES ( A ESTOS ATOMOS SE LES CONOCECOMO ISOTOPOS DEL ELEMENTO ), AUN SI ESTOS SON MUCHOS OALGUNOS NEUTRONES PARA UN NUMERO DE PROTONES DADO, ELRESULTADO SON NUCLEOS NO ESTABLES Y SUFREN UNA DESINTEGRACIONRADIACTIVAI LOS ISOTOPOS DEL OXIGENO ( 2-8 1 CON N = 819 Y 10SON ESTABLES, PERO LOS ISOTOPOS CON N I: 5,697911 Y 12 SONRADIACTIVOS+ EL CASO DE LOS ISOTOPOS CON N I: 596 Y 7 NO TIENENLOS NEUTRONES SUFICIENTES PARA LA ESTABILIDAD, Y LOS ISOTOPOSCON N 11 Y 13, TIENEN DEMASIADOS NEUTRONES+NUCLEOS COMO EL 0-15, QUE ESTA FALTO DE NEUTRONES, SUFREDESINTEGRACION B+ EN ESTE PROCESO, UNO DE LOS PROTONES EN ELNUCLEO, ES TRANSFORMADO EN UN NEUTRON, Y SE EMITE UN POSITRONY UN NEUTRON+ AS1 EL NUMERO DE PROTONES SE REDUCE DE 8 A 7,AS1 RUE EL NUCLEO RESULTANTE ES UN ISOTOPO DE NITROGEN0 N-15,QUE ES ESTABLE+ ESTA TRANSFORMACION SE ESCRIBLE COMODONDE (st SIGNIFICA EL POSITRON EMITIDO, QUE ES COMUNMENTELLAMADO RAYO-BETA, Y 4 ES EL NEUTRINO+ COMPARANDOLO, CON UNNUCLEO DE 0-19, QUE TIENE NEUTRONES EXCESIVOS, SE DESINTEGRAPOR DESINTEGRACION /3' 9 EMITIENDO UN ELECTRON NEGATIVO Y UNANTINEUTRINODONDE 3 ES EL ANTINEUTRINO+ EN ESTE CASO, UN CAMBIO DELNEUTRON EN UN PROTON, Y EL NUMERO ATOMIC0 CRECE EN UNA UNIDAD,SE DEEE NOTAR QUE EN AMBOS CASOS DE DESINTEGRACION fiC O @', EL NUMERO DE MASA ATOMICA RUE QUEDA ES EL MISMO+EN AMBAS FORMAS DE DESINTEGRACION (3 p LA EMISION DE

ELECTRONES APARECE EN UN ESPECTRO CONTINUO DE ENERGIA QUE SEMUESTRA EN LA FIGURA 7tl.tLA ORDENADA EN LA FIGURA, N(E), ES IGUAL A EL NUMERO DEELECTRONES EMITIDOS POR UNIDAD DE ENERGIA CON LA ENERGIACINETICA Et AS1 EL NUMERO ACTUAL DE ELECTRONES EMITIDOS CONENERGIA CINETICA ENTRE E Y E+dE ES N(E)dE* SE DEBE NOTAR EN LAFIGURA QUE HAY DEFINIDA UNA ENERGIA MAXIMA, EH~X,APROXIMADAMENTE CUANDO LOS ELECTRONES NO SON OBSERVADOS* ESTOMUESTRA QUE LA ENERGIA PROMEDIO DE LOS ELECTRONES c ESAPROXIMADflMENTE IGUAL AeO+3 EH~XEN EL CASO DE DESINTEGRACIONEN DESINTEGRACION , E 5 O .4 EH~XFRECUENTEMENTE, EL NUCLEO DESCENDIENTE, QUE ES EL NUCLEOFORMADO EN LA DESINTEGRACION fi 9 TAMBIEN ES INESTABLE Y SUFREDECAIMIENTO b ESTE CONTINUA LA CADENA DE DESINTEGRACION COMOYUN NUCLEO AL QUE LE FALTAN NEUTRONES TAMBIEN PUEDEINCREMENTAR SU NUMERO DE NEUTRONES POR “CAPTURA ELECTRONICA”+EN ESTE PROCESO, UN ELECTRON INTERACTUA CON UNO DE LOSYROTQNES EN EL NUCLEO, Y UN NEUTRO SE FORMA DE LA UNIONt ESTOPERMITE UNA VACANTE EN EL ELECTRON QUE DEBERA SER OCUPADA POROTRO ELECTRONI ESTO CONDUCE A LA EMISION DE RAYOS-X, QUE SON,NECESARIAMENTE CARACTERISTICAS DEL ELEMENTO DESCENDIENTE, OPOR LA EMISION DE UN ELECTRON . HUECO,GENERALMENTE EL ELECTRON QUE ES CAPTURADO POR EL NUCLEO ESEL ELECTRON MAS INTIMO O ELECTRON -DE LA CAPA Y\, Y AS1 ESTEMODO DE DESINTEGRACION ES LLAMADO CAPTURA- KI ENTONCES ELNUCLEO DESCENDIENTE PRODUCIDO EN LA CAPTURA ES EL MISMO QUE ELNUCLEO FORMADO EN LA DESINTEGRACION (3OTRA MANERA POR LA CUAL ALGUNOS NUCLEOS INESTABLES SUFRENDESINTEGRACION RADIACTIUA E.S POR LA EMISION DE UNA YARTICULAALFA t ESTA PARTICULA ES UN NUCLEO ALTAMENTE ESTABLE DELISOTOPO He-4, QUE CONSISTE DE DOS PROTONES Y DOS NEUTRONESI LAEMISION DE UNA PARTICULA ALFA REDUCE EL NUMERO ATOMIC0 POR DOSY EL NUMERO DE PIASA EN CUATRO* ASI, POR EJEMPLO, LADESINTEGRACION ALFA DEL U-238 CONDUCE AL The-234 DE ACUERDO CONLA ECUACION,a39u ?a -> “3qThQo =t :tieLA DESINTEGRACION POR LA EMISION DE PARTICULAS ALFA ESRELATIVAMENTE RARA EN NUCLEOS LIGEROS, PERO ES COMUN EN LOS

NUCL OS PESADOlSt EN CONTRASTE PON LA DESINTECRACION (3 9 LASCUCAS ALFA SON EHITIDAS a UN ESPECTRO DE EWERGIAPART 1DISC'kETO ( UNA LINEA 1.LOS NUCLEOS QUE SE FORMAN POR EL RESULTADO DEDESINTEGRACION 6 ( + O - 1 9 CAPTURA ELECTRONICA, ODES1 TEGRACION ALFA, SI TAMBSEN QUEDAN EN ESTADO EXCITAQO, LAFB TRAN FORMACION CONTINUARA, EL NUCLEO EXCITADO (DESCENDXENTE)GENERALMENTE SE ,DESINTEGRA CON LA EMISION DE UNO O MAS RAYOSGAMMht UN EJEHPLO DE ESTA SITUACION SE MUESTRA EN LA FIGURA7.2, PARA LA DESINTEGRACION DEL CO-60, EL DIAGRAMA QUE SEMUESWRA EN LA FIGURA 7.2, SE CONOCE COMO ESQUEMA DEDESIhTEGRACION* SE DEBE NOTAR QUE LA MAYORIA DE LOS RAYOS-GhMMA SON EMITIDOS POR LOS NUCLEOS DESCENDIENTES, EN ESTE CASOEL Ni-60,C A L C U L O S D E R A D I A C T I V I D A DLOS CALCULOS DE LA DESINTEGRACION RADIACTIUA DE NUCLEOSESTA GOBERNADA POR UNA LEY QUE MANIFIESTA, QUE LA PROBABfLIDADPOR UNXDAD DE TIEMPO QUE EL NUCLEO SE DESINTEGRE ES CONS'f'ANTEtESTA CONSTANTE SE LLAMA " CONSTANTE DE DESINTEGRACION " Y SEDENOTA POR ")r +CQNSIDERANDO DESINTEGRACION DE UN MATERIAL RADIACTIVOt SI AUN TIEMPO t HAY N(t) ATOMOS QUE NO SE HAN DESINTEQRADO,REVI$ANDü LA DEFINICION DE LA CONSTANTE DE DESINTEGRACION,N(t)dt ESTE SE DESINTEGRARA, EN EL INTERVALO DE TIEMPO DG dt Yt+dtr LA RELkCION EN QUE LOS ATOMOS SE DESINTEGRAN ES POR LOTANTO N(t) DESINTEGRACIONES POR UNIDAD DE TIEMPO, ESTARELACION DE DESINTEGRACION ES LLAMADA ACTIVIDAD DE LA MUESTRAY SE DENOTA POR A, AS1 LA ACTIVIDAD EN EL TIEMPO ESTA DADA PORLA UNIDAD DE MEDIAD DE LA ACTIVIDAD ES EL CURIE, DOMDE UNCURIE SE DENOTA POR Ci, Y SE DEFINE COMO 3JE+l0DESINTEGRACIONES POR SEGUNDO* PARA DESCRIBIR PEWENASACTIUIDADES HCIY EL MICROCURIE*^^, LE-^ CURIES,DENOTADO PORfCiDEESENTONCES N(t)dt, Y SE SIGUE QUE LA DISMINUCION DEL NUMERONUCLEOS QUE NO SE DESINTEGRAN EN LA MUESTRA EN EL TIEMPO dtINTEGRANDO LA ECUACION DADONDE NO ES EL NUMERO DE ATOMS A trot DONDE LA ECUACION PARALA ACTIVIDAO DE LA MUESTRA EN EL TIEWPO t, ES3.3

.IFIGURA 7.1ENERGIAIESPECTRO TIPIZO DE LOS XLECTRGNES ETiTTIDGSEN LA DI!10.05860.0co2.5062.1580.24 $1.332INi-60O0FIGURA 7.2ESQUEMA DE DESINTEGRACION DEL (20-60I( SE ILUSTRAN LAS RADIACIONES EMITIDAS)

DONDE A0 ES LA ACTIVIDAD A tsO+ AS1 LA ACTIVIDAD DISMINUYEEXPONENCIALMENTE CON EL TIEMPO+EL TIEMPO EN QUE LA ACTIVIDAD CAE POR UN FACTOR DE DOS SECONOCE COMO "VIDA MEDIA" Y ESTA DADA POR EL SIMBOLO T1/2+ CONESTA DEFINICIONASI, APLICANDO LOGARITMO Y RESOLVIENDO LA ECUACION PARA T1/2,f. YAHORA LA VIDA PROMEDIO is ,DE UN NUCLEO RADIACTIVO ESTAKELfiCIONADA CON LA CONSTANTE DE DESINTEGRACION POR LA FORMULAt = '1%DE LA EC+(7+1) SE PUEDE VER RUE LA VIDA PROMEDIO, LA ACTIVIDADCAE A l/e DE LA ACTIVIDAD INICIAL+ LA RELACION RUE EXISTEENTRE LA VIDA MEDIA Y LA VIDA PROMEDIO ES,EL COMPORTAMIENTO EXPONENCIAL DE LA DESINTEGRACIONRADIACTIVA SE MUESTRA EN LA FIGURA (7+3), DONDE SE INDICA LAVIDA MEDIA Y LA VIDA PROMEDIO+FRECUENTEMENTE SE CONSIDERAN PROBLEMAS EN QUE UN NUCLEORADI&TIVO ES PRODUCIDO A UN REACTOR NUCLEAR o POR UNACELERADOR DE PARTICULAS, SI SUPONEMOS POR SIMPLICIDAD QUE ELNUCLEO SE PRODUCE A UNA RELACION CONSTANTE R ATOMOS/SEG+ ELCAMBIO EN EL NUMERO DE ATOMOS DE UN NUCLEO EN EL TIEMPO dtESTA DADO PORESTA ECUACION AL SER INTEGRADA QUEDA,DONDE NO ES EL NUMERO DE ATOMOS KADIACTIUOS EN EL TIEMPO t z O +LA ECUACION CORRESPONDIENTE A LA ACTIVIDAD ES,

ACTIVIDAD AO/L&O/€"\T1/2 tFIGURA 703 COT.PORTN.'IENTO DE LA DESINTEGRACION DEU?A* IJUESTRA RRDIAC!JZVA* .

SI AOzOj SE OBSERVA QUE, A SE INCREMENTA DESDE CERO, Y PARAt - I -'*..**aj A SE APROXIMA AL VALOR MAXIM0 AMGBX= Rt SI AOf0, AS1LA ACTIVIDAD DEBIDO A LA DESINTEGRACION DE LOS ATOMOSPRESENTES ORIGINALMENTE SE SUMA A LA ACTIVIDAD DEL NUCLEOPRODUCIDO NUEVAMENTE+ EN AMBOS CASOS, LA ACTIVIDAD SE APROXIMAQAL VALOR Ama:

MUESTRA EL LISTADO+7.4 USO DEL. PROGRAMA (ENTRADAS Y SALIDAS)EL PROGRAMA PUEDE MANEJAR UNA CADENA DE DESINTEGRACION CONUN MAXIM0 DE CINCO ELEMENTOS RADIACTIUOS, SIN INCLUIR ALDESCENDIENTE ESTABLE,NOTA: TODOS LOS FORMATOS DE LECTURA SON LIBRES+LA PRIMERA PREGUNTA DEL CODIGQ ES;- rl: DE ELEMENTOS RADIACTIUOS EN LA CADENA (MAX 5)?+ AQUI ELUSUARIO ESPECIFICA LA LONGITUD DE LA CADENA (MI, DANDO COMORESULTADO LOS SIGUIENTES CASOS+PARA M-1, SE DEBERA ANOTARNits), NUMERO INICIAL DE CITOMOS PADRES+ IR, RELACION CONSTANTE DE PRODUCCION DE PADRE (ATOM/SEG), SIR ES CERO EL PRORLEMA SE REDUCE A UNA CADENA SIMPLE DEDESINTEGRACION+Tl, VIDA MEDIA DEL ELEMENTO PADRE+IPARA M=2 (CASO ESPECIAL)+SI SON DOS LOS ELEMENTOS RADIACTIUOS EN LA CADENA, ELUSUARIO PUEDE ESTUDIAR PROBLEMAS MAS COMPLICADOS QUE TIENENDOS CARACTERISTICAS ADICIONALES, DE LAS INDICADAS PARA LOSDEMAS CASOS:a) INICIALMENTE EL DESCENDIENTE TIENE N2(0) ATOMOSPRESENTES+PUEDE TOMAR EL VALOR DE CERO,b) EL DESCENDIENTE TAMBIEN SE PRODUCE A UNA RELACIONCONSTANTE DE P ( CITOM/SEG ), P PUEDE SER CERO,EN ESTE CASO SE DEBERA ANOTARN1(0), NUMERO INICIAL DE ATOMOS PADRES,N2(0), NUMERO INICIAL DE ATOMOS DESCENDIENTES,Rj RELACION CONSTANTE DE PRODUCCION DEL PADRE,P, RELACION CONSTANTE DE PRODUCCION DEL DESCENDIENTE+

Tí, VIDA MEDIA DEL ELEMENTO PAORE*T2.9 VIDA MEDIA DEL ELEMENTO DESCENDIENTEtPARA LOS CASOS M=3,4,5 UNO DEBE SUPONER QUE A t=O, SOLOEXISTEN LOS ATOMOS PADRES, Y QUE EL UNICO MEDIO DE PRODUCCIONEN LA CAOENA ES POR DESINTEGRACION RADIACTIVA, EXCEPTO PARALOS PADRES QUE SE PRODUCEN A UNA RELACION CONSTANTE, R (ATOMOS/SEGPARA M-3N(O), NUMERO INICIAL DE ATOMOS PADRES*R, RELACION CONSTANTE DE PRODUCCION DEL PADRE*T1, VIDA MEDIA DEL ELEMENTO PADRE.T2, VIDA MEDIA DEL PRIMER DESCENDIENTE,73, VIDA MEDIA DEL SEGUNDO DESCENDIENTE+PARA M=4LO MISMO QUE PARA Ms3,PERO ADEMAS,T4, VIDA MEDIA DEL CUARTO DESCENDIENTEtPARA MISLO MISMO QUE EN Mz4, PERO ADEMAS,ITJ, VIDA MEDIA DEL QUINTO DESCENDIENTEt133 CODIGO DE LA UNIOAD DE TIEMPO

MODIFICARAN LOS DATOS MAL ANOTADOS, PARA QUE EL CODIGOCONTINUE SE DEBERA ANOTAR “S“+PARA ADQUIRTR LOS CALCULOS HECHOS, EL USUARIO DEEE ANOTARCUATRO PARAMETROS:1+ EL CODIGO DE LA UNIDAD DE TIEMPO PARA LOS CALCULOS,EL USUARIO DEBERA ANOTAR EL CODIGO DE TIEMPO ( ANTESMENCIONADODESEADAS) PARA TENER LOS CALCULOS EN LAS UNIDADES DE TIEMPOZ+ TIEMPO INICIAL DE CALCULOScESTE ES EL TIEMPO ( EN LAS UNIDADES ESPECIFICADASANTERIOMENTE 1 PARA QUE EL PRIMER CALCULO SE REALICE+ ESTEPUEDE SER CERO O DIFERENTE DE CERO+3+ INCREMENTO DE TIEMPO PARA LOS CALCULOS,EL INCREMENTO DE TIEMPO ( EN LAS UNIDADES YA ESPECIFICADAS), ES EL TIEMPO QUE SE INCREMENTA EN CADA CALCULO,4+ NUMERO DE CALCULOS+ESTE ES EL NUMERO DE CALCULOS QUE EL CODIGO HARA,COMENZANDO CON EL TIEMPO INICIAL, E UNCREMENTANDO CON ELINCREMENTO DADO, HASTA COMPLETAR EL NUMERO DE CALCULOS,SOLO 20 CALCULOS SE HAHAN+YA QUE TODOS LOS CALCULOS FUERON HECHOS, EL CODIGO MUESTRALA SIGUIENTE PREGUNTA “QUIERES ANOTAR OTROS DATOS (S/N)?”, YDEPENDE DE LO QUE EL USUARIO DESEE, EL PROGRAMA WOLVERA ACORRER O NO+S A L I D A S.............................................LAS SALIDAS DESPUES DE ANOTAR LOS DATOS SON, LOS CALCULOSDE N ( NUMERO DE ATOMOS 1 Y A ( ACTIVIDAD 1 PARA CADA ELEMENTODE: LA CADENA Y PARA CADA INCREMENTO DE TIEMPO, ENTONCES SEMUESTRA LA SIGUIENTE TABLA+TIEMPO(SEG) N A0+10 El0 O + 763E91.0N1Al39

DESPUES DE QUE EL NUMERO DESEADO DE CALCULOS SE COMPLETAjEL PROGRAMA GRAFICA LA ACTIVIDAD CONTRA TIEMPO Y REGRESA ALPUNTO DONDE PREGUNTA UN NUEVO CODIGO DE UNIDAD DE TIEMPO PARANUEVOS CALCULOS, DE ESTA HANERA EL USUARIO PUEDE EXAMINARDIFERENTES SECCIONES DE TIEMPO, INTRODUCIENDO NUEVOS CODIGOSDE TIEMPO, O PUEDE CONTINUAR CON EL ANTERIOR Y VARIAR ELTIEMPO INICIAL Y EL INCREMENTO DE TIEMPO. SI EL USUARIONECESITA TERMINAR EL PROGRAMA, DEBERA ANOTAR CERO PARA ELCODIGO DE LA UNIDAD DE TIEMPO*7,s SIMULACION7+5tl CADENAS UTILIZADAS EN LA SIMULACIONILAS CADENAS QUE SE MUESTRAN A CONTINUACION FUERONUTILIZADAS 'PARA LA SIMULACION DE DECAY, TAMBIEN SE OESER'JANLOS DATOS PARA CADA UNA,N10 = 1E6R 2t775 EllT = 64+8, 3 (HORAS)Ai-1-198 -----'.# Ai~-198 (ESTABLE)PARA ESTA REACCION SE UTILIZAN DIFERENTES DATOS DEPENDIENDODEL COMBUSTIBLE UTILIZADO+N1QN20RPT1T2U-233 U-235 PU-2391E6 lE6 1E6O O O7*8Ell 1 14ElZ 1 +4E121+7Ell 4t22E10 2+4E12ár7j 3 6+7, 3 6,7939*2, 3 9*2, 3 9,293PARA EL U-235 SE HICIERON DOS CORRIDAS CON OIFERENTE CODIGODE TIEMPO,SM-157 ---- :> E~J :* Gd {ESTABLE)N10N20RPU-235IE3O1 + 25E9O

T1T2o t5,2XCit2,3~ ~ - 1 4---_-9>a Sfi-149(ESTABLE)U-233 u-235N10 1E6 1E6R1 t3Ell 1+9El1TL.c44,3 55,37t5+2 RESULTADOSEN LOS RESULTADOS DADOS A CONTXNUACION, SE OBSERVA QUEPARA LA PRIMERA SIMULACION SE DA LA CORRIDA COMPLETA, PARA QUEEL USUARIO LA VEA, PERO PARA LAS DEMAS SIMULACION QUE SEHICIERON SE MUESTRAN SOLAMENTE LAS GRAFICASt7+5t3 GRAFICASA ESTE CODIGO SE LE ADAPTO PLOT COMO SUBRUTINA, PARA QUE SEOBSERVARAN DE UNA FORMA GRAFICA LOS RESULTADOS, COMO SE PODRAOEtSERVAH, LAS GRAFICAS SON DE ACTIVIDAD DEL PADRE AS1 TAMBIENLA DE SUS DESCENDIENTES, CONTRA EL TIEMPO* NO SE GRAFICO ELNUMERO DE ATOMOS CONTRA EL TIEMPO, OEBIDO A QUE EL NUMERO DEATOMOS ES LA ACTIVIDAD DIVIDIDA POR UNA COSNTANTE LA DEDESINTEGRACION I+

___..,....O....r....O....,....O....,....O....p....O... ', .." 0 . " . v . " . 0 ....,...'O

.LI

I_-..,....O....r....0....,....0....,....O....T....O...,.....~....T....O....r....Of

;i-- - 1 c..T....p .. .,."'O ..",...'~ ....,.... 0 ....,.... 057 - i l l ) ?05,,? + l 7 o 1157r +I 7~-".,.' .O....,....O... IT""

on, I-& .A- .,,L~ .-.- .A. - . . n . . . . A . . . . 0 . . . . A

-_, . '7"" o ....r....~....,....~....T....~....,....~....,......~....T....~...' V " ' .C

......O.........O....r....O....r....O....r....O...'r.....O....r.."O....r....0.PIJANTT

P 1ci 1Pfr!1' !PIFPP?

?....o . .. .,....o.. '.,.... o . . . ......0....,..'.0.......,.O....,....O....,....[A - . n . . A . .

_-___-_. .. .o. .. ..'..'O . . " ~ . . . ......... ' ~0 ......... 0 " ' . v . " . . O ..'...... o . . . .I ... i-_AN

AI!iI

7t6 CONCLUSIONESLAS CONCLUSIONES QUE SE PUEDEN OBTENER DE ESTE CODIGO ESQUE SIGUE EL COMPORTAMIENTO DE LA LEY DE DESINTEGRACIONES,PERO ADEMAS SE PUEDE UTILIZAR PARA REACCIONES QUE RESULTAN DELA FISION NUCLEARt

ICAPITULO 8+ CODIGO PLOT8+1+ PROPOSITO DEL PROGRAMAPLOT ES UN CODIGO ESCRITO EN FORTRAN PARA TRAZAR EN UNPLANO VARIAS CURVAS EN EL MISMO SISTEMA DE EJES+ UNARESTRICCION IMPORTANTE DE ESTE CODIGO ES QUE LOS VALORES DE LAVARIABLE INDEPENDIENTE DEBE SER IGUAL PARA TODAS LAS CURVASTRAZADAS POR ESTE CODIGO+SUPONIENDO QUE X1 DEPENDE DE T EN ALGUNA FUNCION Y EN UNDIAGRAMA DE T ~5 X1 QUE SEA TRAZADO, LOS VALORES DE T DEBERANESTAN IGUALMENTE ESPACIADOS+PLOT NORMALIZA LOS VALORES DE X1, AS1 QUE EL VALOR ABSOLUTOMAXIM0 DE DE X i SON TRAZADOS CONTRA T+8+2+ CAPACIDADES DE PLOT1) PLOT PUEDE COLOCAR HASTA 10 CURVAS EN EL MISMO SISTEMADE EJES,2) PUEDEN GRAFICARSE hASTA 200 PUNTOS POR CURVA+3) EL USUARIO PUEDE ELIGIR EL LUGAR DEL ORIGEN TAL QUE SEMUESTRE:3) 301.,0 EL PRIMERO Y SEGUNDO CUADRANTE+h) SOLO EL TERCER Y CUARTO CUADRANTE+C) TODOS LOS CUADRANTES,COMO SE MUESTRA EN LA FIGURA 8+1+4) LA LONGITUD DEL EJE X PUEDE SER TAL QUE:EL) SE AJUSTA A UN PAPEL ESTANDAR DE 8-112 PULG+ (OPCIONI)+ EN ESTE CASO SON 60 SUBDIVISIONES A LO LARGO DEL EJE X+b) PUEDEN TENERSE 100 SUBDIVISIONES (OPCION 2) A LO LARGODEL EJE XI EN ESTE CASO LA LONGITUD ES DE 11-1/2 PULG+5) EL USUARIO PUEDE ANOTAR SUS CONSTANTES DE NORMALIZACIONPARTICULARES PARA CADA CURVA O PUEDE DEJAR LA QUE EL CODIGCICALCULA; EN CADA CASO LA CURVA SERA NORMALIZADA POR su VALORABSOLUTO +8+3+ USO DEL PROGRAMA (ENTRADAS Y SALIDAS),ENTRADAS

FIGURA 8.19 .OIX0.0I* OPCIONT0.0X-1.0ZeOPCIONT1 .oX' -1.039 OPCIONT

EL CODIGO PREGUNTA AL USUARIO PARA SUMINISTRAR LAS ENTRADASSIGUIENTES ( TODOS LOS DATOS DE ENTRADA ESTAN EN FORMATO LIEREY DEBEN SER SEPARADOS POR UNA COMA):I1, ANOTE EL NUMERO DE CURVAS Y EL NUMERO DE PUNTOS PARACADA CURVA, EL USUARIO DEBE SUMINISTRAR ENTEROS, EL PRIMERO(MI DEBE ESTAR EN EL INTERVALO DE 1-5: EL SEGUNDO ENTERO (N)DEBE ESTAR EN EL INTERVALO 1-200+ CADA CURVA SERA TRAZADASIMULTANEAMENTE DEBIENDO TENER EL MISMO NUMERO DE PUNTOS+2+ ANOTE LA CONSTANTE DE NORMALIZACION PARA CADA CURVA, ELUSUARIO DEBE ANOTAR LOS VALORES EN COLUMNA, EL PROGRAMAUTILIZA EL PRIMER VALOR PARA NORMALIZAR LA PRIMERA CURVADIVIDIENDO TODOS LOS VALORES DE X l POR LA CONSTANTE DENORMALIZACION DADA, Y AS1 SUCESIVAMENTE+ EL USUARIO PUEDEUTILIZAR LA OPCION DEL CODIGO, ANOTANDO CEROS PARA CADACONSTANTE DE NORMALIZACION, EN ESTE CASO EL VALOR MAXXMOABSOLUTO DE CADA CURVA ES ENCONTRADO Y ES USADO PARANORMALIZAR+3, ANOTE LA OPCION DE TRAZO, PARA LOCALIZACION DEL ORIGEN+EL USUARIO DEBE ANOTAR 192 O 31) EL CODIGO IMPRIME X DE O+O --> 1+02) EL CODIGO IMPRIME X DE -1tO e-2 . 0,o3) EL CODIGO SMPRIME X DE -1,O -->a 1tO4, ANOTE EL LETRERO PARA EL EJE T EN CINCO LETRAS+ ELUSUARIO DEBE ANOTAR CINCO SSMBOLOS CUALQUIERA PARA EL EJE T+3+ ANOTA EL CODIGO PARA LA LONGITUD DEL EJE (1~60: 2=100)+EL USUARIO DEBE ANOTAR 1 O 2+LA PRIMERA OPCION TRAZA EL EL EJE X DE LONGITID DE 60COLUMNAS TAL QUE LA GRAFICA PUEDE AJUSTARSE A UN PAPELESTANDAR DE 8-112 PULG+ POR 11 PULG+ PARA SER IMPRESO A LASALIDA +LA SEGUNDA OPCION ES DE 11-112 PULG+ ; EN ESTE CASO LAIMPRESION DE LA SALIDA DEBE SER EN UN PAPEL GRANDE DECOMPUTADORA,6 t ANOTA EL LETRERO PARA CADA CURVA (MAXIM0 20 LETRAS PARACADA CURVA)+ EL USUARIO PUEDE ANOTAR EL LETRERO PARA CADACURVA GRAFICADA, ANOTANDO HASTA 20 LETRAS+7, ANOTA Xl,X2,++rtXM,T+ EL USUARIO DEBE ANOTAR N NUMEKOSDE LINEAS QUE DEBE CONTENER MI1 COLUMNAS+ EL PRIMER VALORAPUNTADO EN CADA LINEA SE REFIERE AL VALOR DE X1 PARA LAPRIMERA CURVA, Y AS1 SUCESIVAMENTE, EL ULTIMO VALOR APUNTADOEN CADA LINEA ES EL VALOR DE LA VARIABLE INDEPENDIENTE T PARAEL PUNTO EN ESTUDIO,8, QUIERES MODIFICAR ALGUN PUNTO (S/N)t EN CASO DE QUE ELUSUARIO NO COMETIERA ERROR EN ALGUNA LINEA DEBERA ANOTAR N,CON LO QUE EL PROGRAMA CONTINUARA PARA OBTENER LA GRAFICA, EN

CASO CONTRARIO DEBERA ANOTAR S Y EL CODIGO RESPONDERA CON LASSIGUIENTES PREGUNTAS*9* ANOTA LA LINEA QUE MODIFICARAS+ EN ESTE CASO EL USUARIODEBERA ANOTAR EL NUMERO QUE APARECE A LA IZQUIERDA DE LASLINEAS DE DATOS, DE LA LINEA QUE MODIFICARA*10+ ANOTA TODA LA LINEA MODIFICADA* EL USUARIO DEBERA ANOTARLA LINEA COMPLETA DE DATOS CORRECTOS DE LA MISMA FORMA COMO SEANOTARON LAS LINEAS ANTERIORES+ POR EJEMPLOEN ESTA LINEA SE DESEA CAMBIAR EL NUMERO 27 POR EL 28,ENTONCES SE PROCEDE DE LA SIGUIENTE FORMA, EN LA PREGUNTA 8 SEANOTA Si EN LA PREGUNTA 9 SE ANOTA 11 EN LA PREGUNTA 10 SEANOTA 34,28912 ;Y SI ES LA UNICA MODIFICACION A LOS DATOS,CUANDO EL CODIGO HACE DE NUEVO LA PREGUNTA 8, EL USUARIO ANOTAN, EN CASO CONTRARIO CONTINUA CON LA PREGUNTA 9 ,SALIDAS-....-.".---LA SALIDA APARECE EN LA PANTALLA, EN CASO DE QUE LA SEGUNDAOPCION SEA USADA, POR EL TAMANO DE PLOT, LA IMPRESION DE LASALIDA DEBE SER OBTENIDA A TRAVES DEL IMPRESOR EN UN PAPEL DECOMPUTADORA*SI EL USUARIO NECESITA QUE EL FACTOR DE NORMALIZACION PARATODAS LAS CURVAS EN LA MISMA GRAFICA SEA EL MISMO, SE DEBERAANOTAR EL MISMO FACTOR DE NQRMALIZACION PARA CADA CARVAI DEOTRA MANERA CADA CURVA SERA NORMALIZADA DE ACUERDO CON ELVALOR MAXIMO*I8*4+ DIAGRAMA DE FLUJO,EL DIAGRAMA DE FLUJO SE MUESTRA EN LA FIGURA SIGUIENTE, YESTA LIMITADO SOLAMENTE A LAS OPERACIONES IMPORTANTES, PERO SISE QUIERE MAYORES DETALLES VER APENDICE GI8*5+ PLOT COMO SUBRUTINACOMO PLOT ES UN PROGRAMA QUE GRAFICA VARIAS CURVAS EN UNMISMO SISTEMA DE EJES, SE PODRIA USAR SATISFACTORIAMENTE COMOSUBRUTINA DE ALGUN PROGRAMA PRINCIPAL PARA OBSERVAR Y ANALIZARSUS RESULTADOS MAS FACILMENTE, PARA ADAPTAR PLOT COMOSUBRUTINA SE DEBE TOMAR EN CUENTA EL FUNCIONAMIENTO DE ESTE,LOS PUNTOS DE LAS CURVAS DEBEN ESTAR EN FORMA MATRICIAL "EJEY" ( X(10,200) DONDE 10 ES EL NUMERO MAXIMO QUE PUEDE GRAFICARPLOT EN EL MISMO SISTEMA DE EJES Y 200 ES EL NUMERO MAXIMO DE

.,.I,/I63... I ..-__I.- - .... --*

PUNTOS PARA CADA CURVA), Y EL "EJE Xi',POR UN VECTOR ( -ruoo) 200 ES EL NUMERO MAXIMODEBE ESTAR REPRECENTADODE PUNTOS)+MUCHAS VECES SE TRATA DE ANALIZAR LAS CURVAS DE ALGUNPROGRAMA PRINCIPAL, CON EL MISMO FACTOR DE NORMALIZACION PARATODAS LAS CURVAS, POR LO RUE SUGIERE HACER LA SIGUIENTEMODIFICACION +YA QllE SE CALCULO EL VECTOR DE LOS FACTORES DENORMALIZACION ( XMAX(1) ) PARA LAS CURVAS, ADICIONARXMAXX=XMAX(i)DO ## 1:=2,N## IF (XMAXX+LT+XMAX(I))XMAXX~XMAX(I)8+6+ BXMULACIONPARA LA SIMILACXON DE ESTE CODIGO, SE UTILIZARON DATOS DECODIGOS COMO CHANNEL, THERMA, ROOTS, DECAY, CUYAS GRAFICAS SEMUESTRAN, ADEMAS SE GRAFICO UNA ESTROFOIDE, CUYA ECUACSON ESTADADA POR,DONDE A ES UNA CONSTANTE Y SE TOMO EL VALOR DE i00+ UNAOBSERVACION IMPORTANTE ES QUE PARA PODER GRAFICAR ESTAECUACIONtSE' UTILSZARON DOS GRAFICAS, YA QUE SE TIENEN DOSVALORES DE "y" POR CADA VALOR DE "x"+

___.r...,O....r....@....r.... O....,....O....,....0... ,... .o ...'v..,' o... ' V " " 926

EJE XLA GRAFICA NORMALIZADA-1 0 -08 -06 4 -02 0.0 0.2 o, 0.6 0.8 1.+ ~ o ~ ~ ~ + ~ ~ ~ ~ ~ + -+----,+---em-+ ~ ~ ~ o ~ + ~ ~0.0000 ! 2 1'4.0000 ! ? ! l Id.0000 ! ¿ ! ! 1 I;',o000 ! 2 ! 1 Ib.0000 ! 2 ! 1 I0,0000 ! 2 ! 10,0000 ! 7 . 1 1 Ir3.0000 ! % ! 1. I¿?.o000 ! ? ! 16,0000 i Z ! 1 I0,0000 ! 2 ! ! 1 I4,0000 ! 2 ! 1 I8.0000 ! 2 ! 1 I;I00U30 ! 2 ! 1 1c,neoo ! 2 ! 1 I0,0000 ! 2!1 14,0000 ! 2 I2,0000 L 2 It\,oooo ! i?! 1 14.0000 ! 2 ! 1II0.9000 ! 2 . I 16oOOOO ! 2 1 1 I2.0000 ! 2 I 1 Ia,oooo ! 2 I 1 I1 111.0000 ! 2 1I0.001)o ! 2 1 It-~.000O ! 1 1 I2,0000 I 2 I 1 It3.(7000 t- 1I14,0000 2 t 1+-~---+~~-~-+~--..+--~--+-~---+-----~---w-+--w--+.-~--+---~-+-10 -,Y -.6 -.Y -02 9.0 0.2 0.4 U.6 3.5 1.IT X l XF! x3 x4 x5N = .192€+03N = .192E+03

QUIERES ANOTAR OTROS DATOS (S/N)?

.... . . _. . . . ...- . - .. .. .. . ... .. . . (1. , . .., . . . ,j ~. .. .. ., , v . , . . i; . . . ' V ' 'L' . o fEJE XLA GHAFICA NORMALIZADA-1. -.El -.6 0.4 -.2 0.0 0.2 0.4 0.6 0.8 1.+-----+-----+-----+-----+-----+---+.----+-----+-----+-----+~----+* I00.0000 :94.0000 ! * * I88.0000 ! * * 1* ! * I82,0000 !76.0000 ! I I70,0000 ! 1 Ih4.0000 ! I I58,0000 ! I II52.0000 !I46.0000 1 * I * I40.0000 ! * I * I34.0000 ! I2a.0000 ! * * I2% . o o o o ! 116.0000 ! * ! * I10~0000 ! *!* I-4.0000 ! I%.O000 ! * 18.0000 ! *!* 114.0000 ! * ! * 1120.00QO * . *I26,0000 ! * I * I32.0000 ! * I * I38.0000 ! * 1 * I44.0000 ! * I * 15O.r)000 ! *I*I56.0000 ! * I * t92*0000 ! * I * I58.0000 ! * I * I74.0000 * 1 *+-----+-~-~-+-----+----.+-----+-----+----+----+-----+~.---+-----+-1. -.rj 0.6 -.4 -.% 0.0 0.2 0.4 0.6 0.8 1.* = CURVA DE %$@- FACTOR DE NDRMALIZACION = . 192E+03* = A= 100 FACTOR UE MURMALIZACION = o 19SE+03-. 100E+03-.940€+02-.8 8 O E + O 2-.820E+U2-.7hOE+02-. 70OE+02-.640€+02-.580E+O2-.520E+02-.460E+02-.4OOE+Ot?-.340E+02-.280E+02-,¿?2OE+O2-.160€+rJ2-.100E+Oi-T X l XI? x3 x4 xs X6..OOOE+001 bSE*O;!.?22€+02.258€+0228 1E+O;Z. 2

*** SALIDAS ***8,0004.3781415Si (MAX) = ,3íiE*08H = 2124,448W= 345 I b92R .I645'E-02RH = 32OE-03DIST ( TT 1GR AFICA NORMAL IZADA0,O 0.1 0.2 0.3 0.4 0.5 0.4 017 0,8 O,? 1,+. -I--- + ----. ..(-..----.+.----.-+-----+-- I._.t.-.+.------+----.-+--- --$.-*6,0000 ! C !-s,sooo ! R v C !-5,0000 ! V c !-*4,5000 ! R ir C !~"4,0000 ! R iJ c !-3, so00 ! R V C !-*3,0000 ! RV C!-2,5000 ! R V C!-2,0000 ! R V C !-1,5000 ! R V C !,-1,0000 ! R V C !-,5000 ! R V C!,0000 ! R V C*so00 ! R V C:,0000 ! R iJ C!í,SOOO ! R V C !2.0000 ! R 02,5000 ! R V C !3,0000 ! RV C !3.5000 ! R V I C !4,0000 ! R V C !4,5000 ! RV C !c,oooo ! V C !5.5000 ! V C !6,0000 ! C !+_....I...._ + -.-.-I- +.." ---.-. 4 ...-- I- + -----. .t --"..-- +". -..".".-*-----* -.".--- + -----.. .C !0.0 0,í 0.2 0.3 0.4 0.5 0.6 0.7 0.8 O.? í,fiN= TEHP DE REFRFACTOR DE NORMALIZACION =U= TEMP DE VAINAFACTOR DE NORMALIZACION =C=TEMP DE COMEFACTOR DE NORMALIZACION =.314E+04.314E*04.314E+04

OISTLA GRAFICA MOHMALIZAOA0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 I.+-----+-----+-----+-o---+--"o-t-~---+-~---+---~~+----~+-----+-1.h495 ! FZ I-105356 ! F3 2C 1-1.4216 ! F 3 21 I-1,3081 ! F 3 2 1 1-1.1943 ! F 3 2 1 1-1.0b06 ! F 3 2 1c I-.9669 ! F 3 2 1c I-.853I : F 3 2 1c I0,7394 F 3 2 1c L-.625ó ! F 3 i? 1c I-05119 ! F 3 2 1 c I-.3951 ! F 3 2 1c . I-02844 ! F 3 2 1c !-01706 ! F 3 2 lC!-.O569 ! F 3 2 1c.O569 ! F 3 2 1C.1706 ! F 3 i? lC!02844 ! F 3 2 1c !03981 ! F 3 2 1c !05119 ! F 3 2 1 c I.625b ! F 3 2 1c !07394 ! F 3 2 1c 1-8531 ! F 3 2 It I.9669 ! F 3 2 1c I1.0606 ! F 3 2 1c I1.1943 ! F 3 2 1 11.3081 ! F 3 2 1 C I1.4218 1 F 3 21I1.5356 1 -F 321 t1.6493 ! F2 f+-----+-~---+-----+----~+--"~-+--~--+-~~--+-----+---~-+--~--+0.0 0.1 0.2 0.3 0.4 0.5 0,h 0.7 0.8 0.9 1.C= TEM EM EL CENTROFACTOR DE MOt?!íALIZACION = . 248Et041= TEM A í/Y ¡)EL RADIOFACTUR ')E NOH~4ALIZACION = .24AEtO42= TEM A 1/2 DEL RAOIOFACTOR DE NORMALIZACION = ,248E+043= TEM A 3/Y UEL RADIOFACTOR DE miOhllALIIAC1í)IJ = .248E+O4F= TEMP IJN &'A010FACTOR DE iJOKI1ALIZACION = .248E+04PRESIOPJ DE SALIDA TEMPERATURA DE SALIUA1582.9 602.9POTENCIA POR IJIUIDAD DE LONGITIJD (KW/FTI 8.57

(I REACTiVIDAD (T=.l) FACTOR DE NORMALIZACXOW = 100E+01It15 (T=.OolI FACTQI? DE NORVALIZACION . 1 OOE+O 11) = 'I (T=.Oúul) FACTOR DE NClHMALIZACIOiJ = .100E+01. 100E*04,8OüE+í)3,600F+03.4OOE+O3,20OE+O3100€+03.¿300E+O2, 6 O O E + O 2.40i)E+Oa,206E+G?100E+Oi?.1)OOE+01.600E+Q1,4OOF+OI0 aOOE+OlT X l xi?082 4. 1 o 3o 135E-03.957E-03. 1ZhE-O2.191E-d227 OE-02,298E-O¿?336E-02, 393E-02.YY7E-O2x3 x4 x5 X6,820E-04. 102E-o3. 134E-03,195E-03.364€-03,644E-03.764€-03.94S€-O3. 124E-02. 1A6E-O2,261E-02,287E-02032iEE-02,371E-02,453E-02

.loo0 !,1500 !.io00 !.2500 !,3000 !.3500 !.4000 !o O O IO O o I0 U D Io o D 1O0 D 1I) o D !n O D I,4500 IOO O I.5ooo !oO D IO = COSTO ALTf.1. IKGU) FACTOR DE NURMALIZACION = .987E+030 = COSTO SiJCl (KGU) FACTOR DE NORMALIZACION = ,9a7~+03D COSTO TOTAL (KGUI FACTOR DE NORMALIZACION .987E+0311!T x1 X2 X3 x4 x5 X6 1.100EtQü .237€+03 ,59AE+03 .853E+03.lSOEtC)ü .254€+03 .4911Et03 ,752EtO3.200E+00 .;>74€+03 .43fEt03 .7OSE+O3O ,6796*03,300E+QO .328€+03 ,34ZE+03 .6716+03,350~tnu .367~+03 .3ii~to3 .678~+03,400Et00 .418E+03 .285€+03 .703E+O3,450Et00 .48qEtO? .262€+03 .751E+ü3.500E+Oü .592€+03 ,243Et03 .835€+03.5SOE+00 .761€+03 .;226€+03 ,987E+03I1 o-.>2 o-->3 o-->(JCIIiiRES ANOTAR UTRCIS DATOS (S/N)?ANOTA EL #[)E CIlRVAS Y EL NIJblERO DE PUNTOS PARA CADA CURVA:ANOTA LA CONSTANTE DE NORMALIZACION PARA CADA CURVA:ANOTA LA OPCION PARA LOCALIZAR EL ORIGEN:ANOTA EL LFTRERO PARA EL EJE T MAX CINCO LETRAS,:Aa1iOTA EL NIVEL PARA LA LONGITUD DEL EJE(1n60;2n100):ANOTA EL LETREPO PARA CADA CURVA (MAX 20 LETRAS):AFJOTA X l r ..*or XM, T

~.-_____. o . ..r.,..()....*....O....,....O... ,,.. . o . . . ...... 0 ..'.

8+7+ CONCLUSIONESPLOT ES UN CODIGO QUE GRAFICA VARIAS CURVAS EN EL MISMOSISTEMA DE EJES, SE DEBE OBSERVAR QUE EL PAPEL SE DEBE GIRAR90 GRADOS PARA VER CORRECTAMENTE LA GRAFICA, POR QUE t (EJE:"x" 1 j ESTA COLOCADO VERTICALMENTE +ESTE CODIGO ES IMPORTANTE DEBIDO A QUE SE PUEOE UTILIZARCOMO SUE3RU'TINA EN OTROS PROGRAMASj Y MUCHAS VECES LOSRESULTADOS GRAFICOS SE INTERPRETAN MEJOR,.

¡,':¡I.. x 1.;. . .)I"./::: /:.' I ' !' . 'i I'; :;:' (>¡

_____7-E l1?3U567d9101 112131415IC'1738192c21¿i?2324252 h27282930313233363536373 9394 ('414243444s4h4748495051S2535455r -1 bs 7CcCCCCCC&CC177329793560977rkiA'JIiFL Cl C3!,lPilTEP C:iIUE IPJ FORTRArJ Ft)k CALCULATIiJG TEMPEHATURfl)ISTHIí!iJTIl7iJ ALIIIUG A COOLkrlT CHANiJEL IN A PllR USING THFSI;.1PLfl '100EL PHESENTED IFU CtlAPTER 8 dF LAIlAKSH, INTEO.T O EJUC L E A H E NG I fJE E H I NG .THE PJOTATIOFJ IJSEü IS IDENTICAL TO THAT OF THE CHAPTER.TijE l1r;lITS USFD ARE FTp HR, BTUt DEG F, P4wJT.TÍiE I JPilT iifllTS FbH A,B,S,H ARE IiJ INCHES, JUT THEY AK(.I'ii~lE9IhTELY TnNVEHTED TO FT. THE UNITS FOR V INPUT APEFT/SEC. IT IS ALScj It~~GiEDIATELY LClrJVERTFD TO FT/iiR,(SELIM SkPlCAKTAR, APPIL 13, L981)'CP'K ='KC'kF

581; 9bQh Ib20364G 'iub57nsb970717273747576777879(308182&3Lj 4o 5136fi7oe89rs o'$192939495'4 6979899O0O1o 2u304o506C! 7I, rjri 91011lé?151 4

_.._. -. v- - ._'E 31slh17le19202122232425i?h2728293cI3132335455?h373839L( 04 10243444546u7484950s1S?c-3 3'; 4555657LjH59GO61bi?b3t> 4b íb h676869707115220777253304004546606364

~7.. ... .;I 47?737075767775736081Hi?8384tr5Yti57bB9990'3 19733949596974'899O0o10263o 4o 5U60708!I 91 o11121314751 OS12012112315016115slb0300222223124025591 1913727873927***

.: 11274567li91 i t11l?131415ik17IS1929¿ ld?232425Ch2 7262930315?3334353637;e3940414243$1 (445dh47as49L'"3 o'i11-12535455'? fY'j 7cCCCI:CCCCCCCCCcCCCCCCC20SATURATION TEMPERATURE AS FUNCTION OF PRESSURE,F 1F7 Y A 7 P'TERMOMIDRAULICA DE UN PWR',/)299 WRITE(br3Ol)301 FoHMAT(~X,'ANOTA COMB 0,D.: VAI. 0.D.: ESPACIADO : ',/,15Y, 'ALTIIQA DEL IJLlCLEO (EN PULGADAS) : ')R t CI I) ( 3, * 1 !IF U, 1) C L , S , ti

Y ~ -35 CI 6 C [:Ai-CIIl ATIiiiJ OF AXI4L C[iOHf)IrJATE i)F CEDTEti OF FIRST IiJCHEf4EtJT,17 C1619 C20 C7c = -1-4 /? o t CiEL L / 2 , oCALCIILATIOIJ OF AVERAGE BULK TEMPERATURE FOR INCREMENT,21 C¿2¿33.i I TF (h,42 14? FJRF*lAT(/,lOX, 'A LA EIJTRADA DEL ljUCLEO *,I)IZQ25 45 F3RMAT[1OX, 'VEL. PROM. DEL CAidAL *,SX,*REYNOLDS d'p3Xr'FACTOR DI26

CCCCCHECK FOR POSSIBILITY OF SURFACE BOILIIJG, IF DETECTED, "SU89COOLEO BíIILItdG'' PRIFJTEO OlJT AND BOILING BURNOUT CHECK INITIATEDCCCCALCtJL ATIOi-J OF CRITICAL M A T FLUX RATIO-CCCCCHECK F09 SOILING SURNOIJT.P R I Ir) T F O 13 11T A 1.1O C A L c UL A T I ON $ TERM I rJ A T E 0IF DETECTED, "SUBCOOLED BURNOUT"cCCCCCCCCCCI:CCCCRECALC'üLATIOiJ OF TEMPERATUHE AT OUTSIDE OF SURFACE CLADDING ANDHEAT 7iiA.lSF-EH CIIEFFICIENT FUI? SlJbCOOLFD HC)ILING CASES-11C1 TSAV=TS4TAV+I ,9*SHFAV**O,?S*EXPHTC=SIfF ABAV)(-PAV/900,0)CALCIILATEOIJ OF TEMPERATURF AT OUTSIOE SURFACE OF FUEL FOR IN-CRF-SizrErlT, i)ETFi?ii"I~dfZr) i3Y AD!IIrJG TEMPEKATIJí7E QISE ACROSS CLADDINGTU Tt'ivlPE ;?ATLiF!t AT OiITSIDE SlJUFACE OF CLAUDING, AN ITERATIONTECH iI9i:t- IS iICEb TO OC3TATrJ THE CLADDING THERMAL CONOUCTIVITYTO WITHIIJ 1X,115 KCLl=AKCL(TSAVI12 O T F = T SA V + iJ TTCLAV- (TSAV+TF)/?.OK C 1.2 = A Y C L ( T C L A V 1125 KCLl=KCL2G'j TCi 120A V * O F II * * 2 O / 8 O / K C L 1 *AL OG (DC L 1 D F U 1-KCLP)/KCLt?)10) 130,130,12SC~LC~)LAT lOiJ OF TEVPERAT!JRE AT CENTER OF FUEL FOR IIJCHEtIEPiT, ANITERATICIrJ TECHMIDUE IS USED TO OBTAIN FlJEL THERMAL CONDUCTIVITYMITHIN 1%.

.r....U....r....O....,....O....,....u... .,. . . $ .' ' V ' -' ' . v!). . . '.5 IiEhrLETT-PACKARU 32201A.7.13 EDIT/3000 FHI, MAR 16, 1584, 2:03 PM'2Q? fjI32'3 3'34'3536373839' 4 o'4 1Q?'4:'44'45'40'47'48'U9s 051*5?53' 5 u:,'i, C'2657cjaSQ5 oG1o ?#6364'65'bh'n 7'58'69'7 o'7 1'72+7 3'7L!'75'76'7717 f3'7ci1 I t f''& 1' 8 2'b 1;; Li:, i

~ -.. ._. . .- .... .. _ .. , , < . .. v . , . . o . . . . v . . . . u . . . . v . . . . o . . . ... .,o....v....u.... V ' . y ' ? ' V ' I

12345678910111213141516171819202122232'4252627282s30313233343s3637383940414243444s464748495051523354555657$CONTROL NULISTC ROOTS: A CODE Tri CALCULATE THE ROOTS OF T:iE REACTIVITY EOüATIíJi\,C (S.5. FEH 24,L981)CDTMENSI~~~J ~d~(7)COMWPJ XLAM(6),BETA(6) ,XtP,FdGRUATA DATU/'tJ'/605 WRITE(6,lO)10FOHi4AT(/7Y, 'kilOT4 EL tJUVEH0 DE GHUP;J DE iJEUTROI"JES RETAHUAULC: '1RFIAr) (S, *I iJGH11FIiRMAT (I 5)WRITE(6,20)20FORMAT (7X, "ANOTTAfJTES DE DESINTEGRACION 'I,*"EMPEZANDO POR LANA:")READ (5,*) (XLAM( I H 121 F OR f.14 T ( b t 15 . 7 )fiRITE(6r¿?S)25FORMAT(ÍX," A LA FRACCIQN DE NEUTRONES RETARDADOS",*Ic PARA CADA G 'I 1REAE)(T,*? (BETAII)rf=lrNGRIWKI TE (6 30)30FORllAT(7X, 'ANUTA EL TIE:IPQ DE GEr!JEI4ACIOiI DE NE'JTtiUiJ: '1REAb(S,*) XLPWRITE (6,35135FORMAT(7X~"AMOTA O PARA CALCULAR LAS RAICES; 1 " I*"PARA CALCULAR LA REACTIVIDAOn)HEAL) (S, *) NXNRITE(6,SZ)52F~HtlAT(/,¿!OX,'DATOS DE ENTRADA YWRITE(6153) NDE S/iLIDA',/)53FQRIYAT(?X, "ELWRITE (6rS4)IMERQ DE GRUPOS DE MEUTWONES RETARDADOS ES:',Il54FORtiAT (?OX, 'LAMBDA', lox, 'BFTA')DO 56 I=l,rJSH5657),BETA(IIX,E11.41159WRITE (brhc?) THETA62FORtlAT (/, 7X, 'LA FR4CCIdF.l TOTAL DE PJEUTK$~:JE~ dETAHGAPnS =*,3E11.4)WRITE(b,43) XLP63FORMAT(7Xr'EL TIEMPO DE GENERACION DE NEUTR3Pd ES: ',7X,E11.4,/]IF(EJX.NE.1)GO TO 813600 FnRi.1AT (7X, "SE EVALUARAM LAS HAICES", /I'A R I TE ( 6 , b 0 1 )601 F OR PI A T I 7 X , I' SE CLILARA LA REACTIVIDAD:",/)823 WRITE (6,602)602 FOHi'iAT (7X, "TIIS DATOS StJN CDRRECTifS(S/N) 3")READ (5,603) DATIIJ603 FORMAT (AI)IF(DATIN.EQ,DATU) TU 605IF(NXoEQo1) GUT0Ct

_- __. v 'E a5859hO616263646566670069707172737475767778798081828304858687808990919293949596979899O0o1O2030405Ob07O80910111213145060C110115120130150C200210300C376379381C500510WRITE (br401FORMAT(/,?X,"ANOTA LA REACTIVIDAD: ''1READ (5,* 1 HOTF(9ü,GE,1.) GOTO 949IF(RO,GT,OI GOTO 110IF(H0,LT.O) GOTO 120GOTO SOYJF4AX=10Q.WP I I i = oROMA X =E O U I hF1 A Y 1IF(HOMAX.GT,HOI GOTO 130WMAX=WMAX+100,GOTfl 115WMAX=O.WMIU=-XLAil(llWGU= (dMAX+WFlIN) /2.R OG U = E Q U ( WG IJ 1DELTA=ROGU/RU*IF(ABS(0ELTA) . GOTO 200IF(HOtiU.LT.RO) GOTO 150W M A X = 'Ai G tlGOT0 130W M I t i = i4 G UGUT11 130PEHIDP=1 ./kv (11b't!RITE(61379) ROIPERIODDO 378 x=~,NGR+~WRITE(6r381) 1,FORMAT (/7X I 'HE ACTIVIDAD = 'I E 1 1 -4 I SX I 'PERIODO = ' I E 11.4)=F ClkblA T ( 2 OX, ' W ( ' v I 1 I ' 1 ' I E 1 1 . 4 )GOTO 501'ANOTA EL PERIODO DEL REACTOR (T > -l/LAMBDA(l) 1:')GUIF(I;JGU.EG,ü) GOTO 999WGU= 1 i./WI;LIIFCWGU,LT.*XLAMI13) GOT0 SO0WGU=l,/WGUW R I TE (6I5 SO ) kiG I I , f?d

1516171819¿?O21222324252627282930313233550 FORFIAT(/?Xp "PERIODO DEL REACTOR = 'rE11,4,3X, 'REACTIVIDAD= *,i3E11.4)GOTri '300999 WHITk(6,700)700 FORf4AT(/r7X, "IJLJIERES ANOTAR OTROS UAldS (S/PJ)'?")READ(T,IOl) OATUR701 FORMAT (Ai)IF(DATUR,NE.DATU) GO TO 605STUPFPJOFIJNCTIOIJ EOU(rJGIJ)COWI~lON XLAM(6) r5ETAC6) r XLPpNGRDY=QeY=WEU*XLP/(l.*WGU*XLDO 12 J=lrNGR12 DY=DY+dETACJ)/(L~JGU+XLAM(Jl)EQU=Y+üY*!4GU/ (1,+VJGU*XLP)RET 11 R IdEND

I1 $CII1?ITRUL fJOL I ST? C F#Ji?ICli: FtEij [IATEÍIIAL AiJO SEPARATIII’I g.IURY, CALCULATIONS FOR3 C Fdh’ICriIiJG IJRAtIItJF-1 B Y GASEüiJS UIFF IJSIOiJ. (S.S. FEd. 27,4 C 1-96115 C P(T) - AFlOUlJT OF PHODlJCT DESIRED IN KGU6 C WX(I1 - z u-23s IiJ COLAS7 C XPP(1) IVITIAL % OF U-235 IN PROOlJCTtl I: XPE(1) - FIrJAL % OF iJ-2-55 IiJ PRí)UUCT4, C kINC(1) - 1:1CiiE:i1ErJT OF % FROM IiJITIAL TO FINAL %l U C XF - i dF !I-235 Ii? URAIJIUM FEED11 C CF - COST üF FEEL) / KGlJ12 C CS - COST OF ONE SEPARATIVE NORK UNIT I KGU13 C A1d - IJUi4REH OF CASES TO BE STUDIED14 c IF ‘-1cIQE 1iiA.J TEiI CASES ARE TO RE STiJOIEú, IiJCREASE THE DIMENSION1 5 C ST4TEI’IE3JT T:) THE ijECIQEI1 IJU:*IBER OF CASES,Ib111’4EhJSIdiJ F ( 101 hX ( 1 O) r XPP ( 1 O) ,E I?JC ( 1 O) XPE ( 1 O)17\)AT A D AT I ad/ ’ fJ /18 V ( X ) = (2,*X-i. ) * ALOG (X/ I: 1 .OX 1 19 WHITE (6131¿U 913 dR I TF (6,5 9 1 )21 5Ül Fi7R .1AT(SX, ‘AiJíjTA EL C!ISTO DEL MATERIAL ALIMENTA00 Y ’,122 *‘C:lbT 1 O t 1 Silk! Eli S/KG:’) 23 HEAI)(Sr*) CFFCS 124 WR I TE (6 t 502)25 502 FURMAT(SX, ‘ANOTA EL ENRIQUECIMIENTO DE ALIMENTACION(1X = 1.0) ‘i26 */,5X ‘ Y EL liuflEYO DE CASUS QUE SEíiAiJ COKRIDUS:’)!27 HEAII (S,*I XF-, At\ir2hI\I=AiJ29 130 131 503 FORI?lAT(SX, ‘ANOTA EL ENRIQUECIMIENTO DE COLASILOS LIMITES DE ‘ , I32 *’ FNVIilkJEC I:~iIF_;iTt)’,33 *‘ Y tI, I;ICiiEiíit_iJTú’,/,5%, ’i,E EIJRIWUECI~IIE~TO PARA LL CASO’,543535 P(I)=l*(I) rEINC(1)37 10 C!JNT I fJl.IE! ,3c< tJqTTE (hrl-l) XF,CF,CS39 N4ITE (or5140 r)’i is I=~,;J41 15 di=xP,= ( I 144 5 20 C 3KTIiJiJEJb I)) GO TCi 4047 ,O.Ol €INC(II=100,48 .I GO TO 40495 o51,., !d b k 1 =vvx ( I ) / 1 O0r: t> 1 = x f’/ 1 do./, 5X, ’ IlJFERIOH Y SUPERIOR DE LA TABLA ’,brJ C53 CS - fGiJr4dEit UF E;JHICHIiJG STAGESSN - ;JCIMBER OF STRIPPING STAGES54 Cs5 C ON SEPARATION FACTOR OF 0.0045 b l))/(XFl/(l-XFl)))!I 7 S.j=590.*ALUti( (XF1/(1-XF 111 / (kXl/(l-bJXl)) 1- I. Z.A” --&-‘,-’*.“-*&---.-5- - .. *- i_ _”_-

I u95 8 T=SiJtS'j Y F IC'= ( X ? 1 $1 1 / ( XF 1 -'AX 1 r) 3S ,iP= ( I/ ( X fJ 11 -V(i~X 111 -FP* ( V ( X F 1 1 -V (YJX 11 )(J 1 FSFr)=FP*P (I Ib;? s JLI=S.;.JP*P (I 1h3 C CFF - CUST UF FEED, BASED ON KGU64 C CSw - COST QF SAU, BASED UN KGUo '3 C TE - Ti-lTAL CílST, BASEL) ON KGUu h c CFLt - CfiST dF FEE[), BASED OiJ 1 GM U0235(3 7 c CSCl - LUST IF SblJIJ, OASED ON 1 GI1 U-23Sb8 C TCU - T3TAL COST, BASE0 OiJ 1 GM U-23569 CFF=CF*FEED70c s w = c s * s w lJ71T L =CF F+C Ski72 CFIi=CFF/(lO*XP)73 Csll=ssbJ/ (lO*XP)74T C U=C F U t C SI175 IIIS=S76 IIS=T77 TF (1IS.LT.T) IIS=IIStl7u Ir (IíI5.LT.S) IIIS=IIIC+l79 NB!ITE (br 12) XP, FEEU, S?tlJ,CFF,CS\J, TC ,CFU,CSU, TCU, IS, 11s80 30 C8NTINUE81 XP=XP+fINC(I)62 GO TU 20!5 7, 40 C 1.) ($1 T I J IJ E$4 s (3 C:iiuT T'li.lE:'i 5 3 FiittIiAT (3X,'CALCULüS U€ LOS COSTOS PARA EL ENRIQUECIMIENTO',86 * ' DE URA POR BIFUSION GASEOSA',/)8 7 4 FdRI.1AT (/ 'MATERIAL ALIMENTADO ='rF6.3,' WT X U-235',//,IbS . lIOX,'COSTO DE MATERIAL ALIMENTADO = S',F6-2,//, !8 9 ¿?IdX,'Cl.)STO Dt UlJA IJfJIUAG [>E TRABAJO SEPARATIVO = $ 'rF6.2,/)í) 5 Fli&*fhT (gX,"CASj;" ,2X,'iruX (% U-235)',2X,'XPP cx U-235)',31 lax, 'EIrd192 I FORMnTI93 9 FORMAT Ii?r//,lOX,'LA COMPOSICION DE COLAS'rFSi34 1' 4T % U-2 'rF4.1,' KGU'r//,lOX,'COSTOS COMO *,'i 5 2'FIiíJCIUIi L>E KÚU Y U-235 HESPECTIVAMENTE8,//,20X,9 h *' TAtiLA íJE LdS COSTOS DE SERVICIO PARA ENHIQUECIMIENTO',//,'$7 3SX, 'Xí",BXI'ALIi1',7X,'S~~U',SX, 'COST ALIM COST SkdU COST 'TOTAL'98 42Xt'COST ALIM COST SHU TOTAL CDST'r4X,'ENRIU 'rYX,'TOTAL'99 5,/,1X,'(NT X U0235)',2X,'(KGU)',lSX,'F(KGU)(~) F(KGU)($) ' e.o0 ó*F(KGUI ($)'2X, ' F(2351 ($1 F(235l ($1 F(235) ($1'r3X,01 7' COCAS'r5X'COLAS'r //ILO¿? 12 F kJK',lA T (4 X F

. _.,. ... .

7. . . .i).. . . . [, ....... ,;u . . v . . . . o . . . ..."O... .,....O... ...'-'u. ' V;t 1123(15b7h4I1112131u15Ih171819¿O212?232425?h27?O2930315?335435363736394 o414s43J4454 í,474 a4950'i 1'1253545s5t.') 7

o . . . ... . . . o . . . 'T .... 0 ... ......o... . ~ . . [ l . ..". . . , . . . .,. . . :)ili-:~iiLETT-1'kLt:#i.'i) 322OlA.7.13 EOIT/3000 FHI, MAR 16, 1984, 2:08 PF1 (C)4204 5 !I460C4 b 7611SSO625630h 3 16506517 (J 07 o 17 1 o711720C730771r\ 0 0HOOO

15115%1334567691 0

727374770777879í) 0U18283A44 5Qb37888990$192'3 39495Yh37'3 H9 c4'O0'O 102(J 3Ii4d5'06I o 7'Obb 1) 9*I o'1 1'1213'1 4'15'Ih17161920212?2324I2526272 811121314151617181920?1222324262728C?q" A

'29I30'5 1'32'33'34'35' 5h'373 ti'39'40'4 1' k i ?'4?'4 4'45'Y6 .!474 H)J9'50'5 1'52'5 3'S4'55l%h'57'58'59' f> o'6 1'h2'63'6 4'65'h6'6 7'bM'69'7 O'7 1)7?'7 3'7 4'7 5'76'77'?ti'7 9'4 o'8 1'82'83Ii441 +' L13 o40s il61)7980901002 o o4467?O01s1

.,.,.. o___~-..'."'... o . . . ......o. .. .."..O... ."." . o - . . . . . . 2 ......... t).- ..-.. .. .. . . . .. . . . . .. .. . __._ .-.-. . - .' V . '.:,123145c,79910111213A

1415? (i24304 0601 2 o1 so200n A n

E l123456789101 112131415161718192021222324252627282930313233343536373839404142434445464748495051Si?5354555657$CONTROL NOLISTCC535507273747s7780

. V . . ' 'E 2585960616263646566b768697071727374757677787980818283H4858687888990919293949596979899O0O1020304os060708091011121314

v€ 315161718192021222324252627232930313233343s3637383940414243444546474849SO5152535455565758596061b26364666760697071b5130613513140142144146148150200210

'

_--- ~.____..._- ___.*.: O . . .-.. .O.........O....,....0....,....~....~....0....,....~... .... o .E 5293031323334353637383940414243444s464748495051525354555b575859606162635465666768697071727374757677787980818283848511315207772225304045466061641051201211231 so161155160200220

8607ut38990919293949596979899O0O1020304225513514515516518636517

~ '.,"_

1234567b9LOI1i21314:516!78.9?O?112?3!4?5?b?7!819so511233I4r5;65738391041+213141s1617+049505152335435563753 FDRI4AT(lOX, 'A;JOTA EL ¡'JIVEL PARA LA LOYJGITIJD DEL EJE'

'.E 2585960bl525354656607fJ 85970717273747576777e79BO61520350859h3788993031923394353b979899O0U1'32O33435363708091011121314303152022253040444660616263646671727374751 os120

1516171619?O?1?2-33435262728193031323334353637383940114243!I495167738$930515253541211231501611551602002102202252302.3 124024525025560060 1\

EIBLIOGRAFIA-. 1+ J+ 'CRUZ Y J+ AREAS, EXPERIMENTUS UE XLNSTRUMEN7'ACION YFISICA NUCLEAR, UAM-I, 1982,-. 3 + ,J ,.I UUIJERS'I'AD'l AND L + J + HAMIL.TON, NUCLEAR REACTURn E! !Ai I... :I: :; x : * E: D w 1: I." E 'Y f :I. '3 7 6 +. S j. + NUCLEAR TE:CI-INOl ... 1lGY 9 UCILIJME 6% 9 NlJMBER 2 AUGUST 1983,F'AG 139,.... I::+ SELXH sANcCit:wm, CHANNEL, A CODE TO CALCULATETI:MPERAT[JFW DISTRIBUTION ALONG A COOLANT CHANNEL IN A rm, ED,SCHOOL. OF AMNE, UNIVERSITY OF OKLAHOMA+13+ CELIM SANCAKTAR, DECA, A CODE IN FORTRAN FORRADIACTIUE DECAY CHAIN CALCULATIONS, ED+ SCHOOL OF AMNE,UNIVERSITY OF OKLAHOMA +14, SELIM SANCAKTAR, PLOT, A CODE IN FORTRAN FOR PLOTTING* MULTIPLE CURVES ON THE SAME SET OF AXES, ED+ SCHOOL OF AMNE,UNIVERSITY OF OKLAHOMA+-* 15, SELIM SANCAKTAR, ROOTS, A CODE TO EVALUATE THE ROOTSOF THE REACTIVITY ERUATIONj ED+ SCHOOL OF ANNE, UNIVERSITY OF2 15

OKLAHOMA +- 16+ SELIM SANCAKTAR, STAT, A CODE IN FORTRAN FOR STATICAL.ANALYSIS OF RADIACTIVITY MEASUREMENTS, ED+ SCHOOL OF AMNE,UNIVERSITY OF OKLAHOMA+I. -.-U_/