2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti...

39
Fare clic per modificare lo stile del sottotitolo Progetto Progetto neutronico neutronico di di massima dell massima dell impianto impianto sperimentale sperimentale . . Carlo Maria Viberti [email protected]

Transcript of 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti...

Page 1: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Fare clic per modificare lo stile del sottotitolodello schema

““Progetto Progetto neutronico neutronico didimassima dellmassima dell’’impiantoimpiantosperimentalesperimentale””..

Carlo Maria [email protected]

Page 2: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Punto di partenza Punto di partenza I:I:il il fuel assemblyfuel assembly

Fuel Assembly

L FA ≅

9.7

cm

∼30 kgUO2-enr./FA

∼45 kgPb/FA

x

y z=0

(hattiva = 90 cm)hFA ≅ 130 cm

Page 3: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

x

y z=0Restrizioni:– keff ~ 0.95

Punto di partenza Punto di partenza II:II:ll’’assemblaggio sotto-criticoassemblaggio sotto-critico

quanti FAs?

quanto riflettore?

(r tot?) vessel

riflettore

FAs

Beamline & target

rtot

Page 4: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Configurazione Configurazione ““pienapiena””

MCNPX 2.6.0, kcode, ENDF/B.-VII

Page 5: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Configurazione Configurazione ““pienapiena””

MCNPX 2.6.0, kcode, ENDF/B.-VII

Page 6: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

rtot=120 cm

Errore sistematicoErrore sistematico

keff = 0.946

+0.007-0.003

!"#

$%&

∼1.8 tUO2-enr.

∼2.7 tPb in FAs

∼84 tPb in rifl.

x

y z=0

Page 7: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Studio del Studio del targettarget

Ottimizzata per p @ 70 MeV, 0.75 mA

p_i

Sbarra di Berillio,

concezione TRADE

∼200 gBe

Page 8: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

(Z(X)=20·X-9.5)

Ottimizzazione dimensioniOttimizzazione dimensionibersaglio bersaglio II

Dimensioneproposta

Page 9: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

(Z(X)=20·X-9.5)

Ottimizzazione dimensioniOttimizzazione dimensionibersaglio bersaglio II

Dimensioneproposta

Ginocchio: oltre il ginocchiosfuggono solo i protonisecondari. Infatti…

Page 10: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Ottimizzazione dimensioniOttimizzazione dimensionibersaglio bersaglio IIII

Energia (MeV)

Partenza

Start

Spettro energeticodei protoni in fugadal target perd i f f e r e n t idimensioni delbersaglio (A, B, …).

Page 11: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Ottimizzazione Ottimizzazione spotspot fasciofascio

3 cm

~ 0.6 cm

r

R = 1.5 cmDimensioneproposta

(Z(X)=20·X-9.5)

Page 12: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Caratterizzazione sorgenteCaratterizzazione sorgente

dE d!

4"!( )E! " 0.091#p _ i $ 4.3 #1014 neutroni sec .

dV

V! " 49.3 kW.

Fondamentale per progettare il sistema direfrigerazione dedicato al target

( 1 MeV/cm3/p_i = 0.75 kWatt/cm3 )

(Corrente fascio: 0.75 mA)

1 20-1-2-10

-50

510

1520

(Corrente fascio: 0.75 mA)

Yield totale dei neutroni di sorgente:

Page 13: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Collocazione assialeCollocazione assialesorgente sorgente II

δδ

Centro del target

Centro dell’active core

Target

δδ Distanza tra i due centriDistanza tra i due centri

qualequale δ δ è è il migliore? il migliore?

BeamlineActive core

hattiva =

90 cm

Page 14: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

δ ≅ 5 cmDimensioneproposta

Collocazione assialeCollocazione assialesorgente sorgente IIII

Page 15: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Deutoni Deutoni ??

Page 16: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Deutoni Deutoni ??

Yield totale: p @ 70 MeV ≈ d @ 50 MeV

Page 17: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

In un A.D.S. la “sorgente” deve essereben conosciuta: vi è la necessità di

studiare sperimentalmente p+9Be perverificare/(eventualmente)correggere i

risultati forniti da MCNPX.

Conoscere Conoscere la la sorgente sorgente !!

“Resa neutronica di protoni da 62 MeV su bersagliodi Berillio: esperimento preliminare”.

M. Ripani, INFN-Ge

Page 18: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Andamento radiale potenzaAndamento radiale potenza

Page 19: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

ENDF/B.-VII:

Ptot ~190 kWth

Andamento radiale potenzaAndamento radiale potenza

JENDL-3.2:

Ptot ~230 kWth Ptot = 199

+29-10

!"#

$%&

kWthTotale:

Page 20: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Distribuzione assiale densitàDistribuzione assiale densitàdi potenzadi potenza

Page 21: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Flusso neutronico barraFlusso neutronico barrainternainterna

Flusso integrale barra interna: φtot ~ 6·1012 n/cm2/sec( p @ 70 MeV, 0.75 mA )

ENDF/B.-VII

Page 22: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Flusso neutronico barraFlusso neutronico barrainternainterna

ENDF/B.-VII

Flusso “lento”:φ<0.5 MeV ≅ 4.4·1012 n/cm2/sec ⇒ φ<0.5 MeV / φtot = (69.8±0.5)%

Flusso “veloce”:φ>0.5 MeV ≅ 1.9·1012 n/cm2/sec ⇒ φ>0.5 MeV / φtot = (30.2±0.4)%

Flusso integrale barra interna: φtot ~ 6·1012 n/cm2/sec( p @ 70 MeV, 0.75 mA )

Page 23: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Simulazioni di bruciamento Simulazioni di bruciamento II

Barra di MA nellaposizione più interna.

MCB-1c

∼20 pcm

Page 24: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Flusso medioparte attiva

Flusso medioriflettore

σ(n, γ)del 99Tc

Flusso neutronico riflettoreFlusso neutronico riflettore

Flusso integrale medio riflettore: φtot ~ 1·1012 n/cm2/sec( p @ 70 MeV, 0.75 mA )

ENDF/B.-VII

Page 25: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Simulazioni di bruciamento Simulazioni di bruciamento IIII

Barra di LLFPnel riflettore.

MCB-1c

∼5 pcm

Page 26: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Ulteriori simulazioniUlteriori simulazioni

A seguito delle simulazione effettuate pervalutare gli O.d.G. delle caratteristichefondamentali reattore, sono stati eseguiti altristudi per comprendere la componente velocedei neutroni e l’effetto di “pozzetti” perl’inserimento di apparati sperimentali.

Page 27: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Dipendenza dallDipendenza dall’’acceleratoreacceleratore

Legnaro - Ciclotrone

p @ 70 MeV, 0.75 mA

Legnaro - IFMIF

p @ 20 MeV, 120 mA

(entrambi ENDF/B.-VII)

Pfascio|target ~ 50 kW

Ptot ~ 190 kWthφtot ~ 6·1012 n/cm2/sec

Percentuale flusso “lento”:φ<0.5 MeV / φtot = (69.8±0.5)%

Percentuale flusso “veloce”:φ>0.5 MeV / φtot = (30.2±0.4)%

Pfascio|target ~ 2.4 MW

Ptot ~ 3.4 MWthφtot ~ 1·1014 n/cm2/sec

Percentuale flusso “lento”:φ<0.5 MeV / φtot = (69.3±0.5)%

Percentuale flusso “veloce”:φ>0.5 MeV / φtot = (30.7±0.4)%

Bar

ra in

tern

a

Page 28: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Dipendenza dai materialiDipendenza dai materialistrutturalistrutturali: studio : studio sul flussosul flusso..

Med

io p

er b

arra

Page 29: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Studio Studio sul flussosul flusso: : tolto tolto He, Al,He, Al,cladding rodcladding rod e e cladding FAcladding FA

Med

io p

er b

arra

Rimane combustibile e matrice di piombo

Page 30: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Med

io p

er b

arra

Studio Studio sul flussosul flusso: : ripristinoripristinocladding FAcladding FA

Manca elio-gas, alluminio e cladding rod

Page 31: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Med

io p

er b

arra

Studio Studio sul flussosul flusso: : ripristino ilripristino ilcladding rodcladding rod

Manca elio-gas e alluminio

Page 32: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Med

io p

er b

arra

Studio Studio sul flussosul flusso: : ripristinoripristinoll’’alluminioalluminio

Manca elio-gas

Page 33: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Studio Studio sul flussosul flusso: : ripristinatoripristinatoanche anche HeHe

Med

io p

er b

arra

Tutto ripristinato

Page 34: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Studio Studio sul flussosul flusso: : ZircaloyZircaloyanziché anziché AISIAISI

Med

io p

er b

arra

AISI

Page 35: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Studio Studio sul flussosul flusso::combustibilecombustibile

Bar

ra in

tern

a

φ>0.5/ φtot ≅ 32.8%Frazione flusso veloce - noOX

φ>0.5/ φtot ≅ 30.1%Frazione flusso veloce - yesOX

Page 36: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

57 FA

keff ≅ 0.913

Ptot ≅ 117 kWth

Pozzi Pozzi per per esperimenti esperimenti I:I:

60 FA

keff = 0.946

+0.007-0.003

!"#

$%&

Ptot = 199

+29-10

!"#

$%&

kWth

keffpiena ! 0.95, keff

bucata " 0.95

Page 37: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

65 FA

keff ≅ 0.934

Ptot ≅ 161 kWth

68 FA

keff ≅ 0.964

Ptot = 298 kWth

Pozzi Pozzi per per esperimenti esperimenti II:II: keff

piena > 0.95, keffbucata < 0.95

Page 38: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

73 FA

keff ≅ 0.952

Ptot ≅ 217 kWth

76 FA

keff ≅ 0.978

Ptot ≅ 492 kWth

Pozzi Pozzi per per esperimenti esperimenti III:III: keff

piena ! 0.95, keffbucata " 0.95

Page 39: 2011 02 07-Conf-ANNricco/ansaldo/Progetto neutronico...dei protoni in fuga dal target per differenti dimensioni del bersaglio (A, B, …). Ottimizzazione spot fascio 3 cm ~ 0.6 cm

Pozzi Pozzi per per esperimenti esperimenti IVIV

Pozzo vuoto più interno

φ>0.5/ φtot~26%Frazione flusso veloce

φ<0.5/ φtot~74%Frazione flusso lento