1 / 66

First year Application on 233,232 U(n,f) and 239 Pu(n,f) Preliminary results

IAEA Research Contract No. 15805 Prompt fission neutron spectrum calculations in the frame of extended Los Alamos and Point by Point models. First year Application on 233,232 U(n,f) and 239 Pu(n,f) Preliminary results. Professor Dr. Anabella TUDORA Bucharest University, Faculty of Physics

betty
Télécharger la présentation

First year Application on 233,232 U(n,f) and 239 Pu(n,f) Preliminary results

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. IAEA Research Contract No. 15805Prompt fission neutron spectrum calculations in the frame of extended Los Alamos and Point by Point models First year Application on 233,232U(n,f) and 239Pu(n,f) Preliminary results Professor Dr. Anabella TUDORA Bucharest University, Faculty of Physics Bucharest-Magurele, POB MG-11, R-76900, Romania IAEA-NDS, April 2010

  2. Content of the presentation I. Basic features of prompt neutron emission models used II. 239Pu(n,f), 233U(n,f), 232U(n,f) preliminary results 1. Point by Point model calculation of prompt neutron quantities: total average prompt neutron multiplicity and spectrum, ν(A,TKE), ν(A), P(ν), < ε >(A), <ν>(TKE) and so on - discussion about FF experimental distributions Y(A,TKE) - 2 methods of TXE partition between the 2 fragments forming a pair - different optical model parameterizations to calculate the compound nucleus cross-section of the inverse process σc(ε) of each fragment - average model parameter values obtained from PbP treatment 2. Most probable fragmentation approach: prompt neutron multiplicity and spectrum calculations 3. Discussion of different fission cross-section evaluations that can be used (as fission c.s. ratios) to calculate PFNM and PFNS at incident energies where multiple fission chances are involved

  3. Basic features of models based on neutron evaporation from fully-accelerated fission fragments (Los Alamos type models) SCM, one fragment SCM, one fragmentation, “most probable fragmentation” approach. :

  4. LS, one fragment Total spectrum Obs.:Most prob.fragmentation around AH=140 where νH = νL, fact observed for all fissioning systems with experimental sawtooth data. Madland & Nix assumption (Nucl.Sci.Eng.81(1982)213) is correct. CN cross-section of the inverse process: DI mechanism, SCAT2 code with optical model parameterizations appropriate for FF nuclei region: Becchetti-Greenless, Wilmore-Hodgson, Koning-Delaroche, as well as the simplified σcform proposed by Iwamoto.

  5. New form of FF residual nuclear temperature distribution (Vladuca and Tudora, Ann.Nucl.Energy 32 (2005) 1032-1046) • With conditions: • continuity in T=α • P(T=β)=0 • normalization on 1 • <T>=2Tm/3 parameterization in the computer code:

  6. Obviously for s=1 P(T) of LA Madland & Nix is re-obtained Anisotropy effect: the most important emission of prompt n. is from fully-accelerated FF but n. evaporation during fragment acceleration is also possible, leading to a non-isotropic spectrum in CMS. Another source of non-isotropic neutrons can be the emission at the scission moment (scission-neutrons). According to Terrel the anisotropy of neutron emission if present is symmetrical about 90o and the SCM spectrum could be described by: b = anisotropy parameter

  7. With the new P(T) and anisotropy taken into account, the prompt neutron spectrum in LS is :

  8. In the case where only one fissioning nucleus is involved: SF and • neutron induced fission in En range of the first chance • “Most probable fragmentation approach” with average values of • model parameters (<Er>, <TKE>, <Sn>, <Eγ>, <C>=Ac/<a>) • b) “Multi-modal” fission concept: total PFNM and PFNS calculated • as superposition of the multiplicity and spectrum of each mode • weighted with the modal branching ratios. Average model parameters • are determined for each mode. • c) Point by Point (PbP) model: the entire FF range covered by the • Y(A,TKE) distribution is taken into account. Total PFNM and PFNS • are calculated as superposition of the multiplicity and spectrum of • each pair weighted with the charge and mass distributions of FF. The PbP treatment is the most accurate because it takes into account the full range of possible fragmentations while the other two approaches consider only one ore few fragmentations (subsets) with average model parameters

  9. When more fission chances are involved: only “most prob. fragm”. approach is used because it is impossible to distinguish Y(A,TKE) of each chance and more, the secondary CN are formed at many excitation energies  too large amount of calculations. Fissioning nuclei of the main chain (formed by neutron evaporation from the precursor) – LA classical:

  10. Evaporation spectrum of neutrons emitted prior to the scission: • Weisskopf-Ewing spectrum as in the classical LA model • evaporation spectrum obtained from (n,xn) spectra provided by • GNASH-FKK (or other codes like Talys, Empire), from which the • contribution of neutrons leading to excitation energies of the • residual nucleus less than the fission-barrier height were substracted. • (Tudora et al., Nucl.Phys.A 756 (2005) 176) Most probab. fragm. at high En – the fission of secondary CN chains formed by charged particle emission is taken into account as following: (Tudora et al., Nucl.Phys.A 740 (2004) 33-58) 2) Protons way (fissioning nuclei formed by p. emission from CN of the main ch. 3) Neutrons via protons way (fiss.nuclei formed by n evap. from the nuclei formed by p. emission 4) Deuterons way (fiss. nuclei formed by d. emission from the nuclei of the main chain 5) Alpha way (fiss. nuclei formed by alpha emission form the nuclei of the main chain 6) Neutrons via alpha way (fiss. nuclei formed by n. evaporation from the nuclei formed by alpha emission

  11. Secondary nucleus chains and paths (ways) – Excitation energies, recursive formulae: Paths: k=2 (p) and k=5 (α) Sk p or α separation energy from “i” precursor of main ch. c=chain II (formed by p, nvp, d) and III (α and nv α) Paths: k=4 (d) Paths: k=3 (nvp) and k=6 (nvα)

  12. Total and partial fission cross-section ratios production cross sections of the j-th secondary nucleus by proton (σ2 j), neutron (σ3 j) and deuteron (σ4 j)emission and respectively by alpha and neutron via alpha emission Prompt neutron multiplicity: k=1(n), 3(nvp) and 6(nvα) i from n=1 for n, n=2 for nvp and nvα k=2(p), 4(d) and 5(α) i from n=1 for (p) and n=2 for (d)

  13. Total prompt neutron spectrum for k=1(n), 3(nvp) and 6(nvα) Total prompt neutron spectrum for k=2(p), 4(d) and 5(α) the individual PFNS of the i-th fissioning nucleus of the k way the evaporation spectrum of neutron emitted prior to the scission for secondary nucleus ways only Weisskopf-Ewing Details of the model in: A.Tudora, G.Vladuca, B.Morillon, Nucl.Phys.A 740 (2004) 33-58

  14. Computer codes: 1. SPECTRUM first version (2000) (Vladuca and Tudora, Computer Phys. Communic.125 (2000) 221-238) Program Library of CPC id. ADLH (2000) - Most prob. fragm, σc=variab, multiple fission chances, main CN chain - Average model parameters dependence on En (or E*) made for SF and neutron or proton induced fission (one incident energy) • 2. SPECTRUM second version (2002) • new P(T) form, anisotropy, possibility to use evaporation spectra • of neutrons emitted prior to the scission provided by GNASH-FKK • or other nuclear reaction codes (like Talys, Empire) • For SF and neutron or proton induced fission (multiple incident energies) 3. SPECTRUM extended version (2004) Including the extended model with fission of secondary CN chains and ways (paths) formed by charged particle emission (high En) All versions allow also the options: σc=const and σc=the simplified form of Iwamoto

  15. 4. SPECTRUM with multiple fragmentations (2002) Including 2 options: A. Multi-modal fission approach B. Point by Point approach Dimensioned for 300 pairs (meaning 600 fragments) • Auxiliary input files for all versions: • σc provided by SCAT-2 code with different optical model • parameterizations appropriated for FF nuclei • Fission c.s. ratios RF taken from evaluations (endf) or from • cross-section calculations (ECIS + STATIS + GNASH-FKK) • For the version with multiple fragmentations: • Model parameters of fragmentations provided by the codes: • MMPAR: for multi-modal fission approach • PAIRPAR: for PbP approach

  16. POINT BY POINT model FF pair range: all mass pairs {AL, AH} covered by Y(A, TKE) (with the step of 1 mass unit) are taken into account. For each mass pair 2Z/A, 4Z/A, 6Z/A… are taken with the charge number Z as the nearest integer values above and below the most probable charge Zp (UCD with charge polarization ΔZ). Usually 2Z/A are taken because P(Z) is a narrow distribution. For each FF pair “i” (meaning LFi, HFi) quantities as: number of prompt neutrons, prompt neutron spectrum, average prompt neutron energy in CMS and so on are calculated. PbP provides all multi-parametric data (quantities referring to each FF) these quantities do not depend on FF distributions Y(A, Z, TKE) And can be compared with experimental data. Such quantities are: the multi-parametric matrix ν(A,TKE), νpair(A), sawtooth ν(A), ε(A), Eγ(A), TXE(A) and so on.

  17. PbP provides all total average quantities (depending on FF distrib.) that can be compared with existing experimental data too. Such as < νp> , N(E), < νp>(TKE), <ε>(TKE), P(ν) etc. • Special mentions: • P(ν) very sensitive quantity - the PbP model is (to our knowledge) • the only able to provide P(ν) results in very good agreement with • experimental data for all fissioning systems having experimental P(ν) data (SF and neutron induced fission) : • 252Cf(SF), 248,244Cm(SF), 240,242Pu(SF), 235U(nth,f) and 239Pu(nth,f) • Tudora, Ann.Nucl.Energy 37 (4) (2010) 492-497 ( for 252Cf, 248,244Cm(SF)) • Tudora, Hambsch Ann.Nucl.Energy (2010) in press (7 fissioning systems) • PbP describes very well experimental < νp>(TKE) data too: • 252Cf(SF), 248,244Cm(SF), 233,235U(nth,f) • Tudora, Ann.Nucl.Energy 35 (2008) 1-10 PbP is a powerful model. The TXE partition between LF and HF forming a pair is almost avoided because the model is working under the concept of FF pair. The only one limitation is the need of Y(A,TKE)

  18. Because at this meeting many results concerning 235U(nth,f) and 252Cf(SF) were presented, I added few of my previous results concerning P(ν) and Eγ(A) Fig.3 from A.Tudora, F.-J.Hambsch, Ann.Nucl.Energy (2010) in press

  19. Fig.1 from A.Tudora, F.-J.Hambsch, Ann.Nucl.Energy (2010) in press

  20. Tudora A., Ann.Nucl.Energy 35 (2008) 1-10

  21. Tudora, Morillon, Vladuca, Hambsch…Nucl.Phys.A 756 (2005) 176-191

  22. PbP model parameters and average parameter values parami: Eri, TKEi, Sni, Ci Systematics of average param. Tudora, Ann.Nucl.Energy 36 (2009) 72 if equal T of LF and HF Level density parameter of each FF Super-fluid model (Ignatiuk) Δ(Z,A) pairing δW shell-corrections (Moller&Nix RIPL)

  23. TXE partition between the two fragments of the pair Lemaire et al. Phys.Rev.C 72 (2005) 024601 • Lemaire,Talou, Madland (Los Alamos) aL, aH super-fluid model, iterative procedure The two FF have the same Tm (thermodynamic equilibrium) • Tudora (Bucharest) Tudora Ann.Nucl.Energy 33 (2006) 1030-1038 aL, aH super-fluid model. νH/(νL+ νH) parameterization obtained on the basis of the systematic behaviour of experimental sawtooth data Not equal Tm of the two fragments

  24. Simulations to obtain Y(A) of FF (pre-neutron) in neutron induced fission Experimental P(Z) data of FF (pre or post neutron) can be used in the frame of the PbP treatment to obtain Y(A) of FF: A.Tudora, Ann.Nucl.Energy 37 (2010) 43-51 • Experimental P(Z) for 234U electromagnetic induced fission • (Schmidt et al., Nucl.Phys.A 665 (2000) 221-267) •  PbP treatment to obtain Y(A) of FF for 233U(n,f) • Experimental P(Z) of post-n FF for 239Pu(nth,f) measured by Gonnenwein (Cosi fan Tutte, Kaufmann et al., Proc.Int.Conf.NDT (1991), 133 ) •  PbP tratment to obtain Y(A) of FF for 239Pu(nth,f) • Exp.P(Z) of 233U electrmg.ind.fiss. and of 232U(nth,f) used to obtain Y(A) of FF PbP model to calculate prompt neutron emission data for232U(n,f)

  25. A.Tudora, Ann.Nucl.Energy 37 (2010) 43-51

  26. A.Tudora, Ann.Nucl.Energy 37 (2010) 43-51

  27. A.Tudora, Ann.Nucl.Energy 37 (2010) 43-51

  28. PROMPT NEUTRON EMISSION DATA FOR 239Pu(n,f), 233U(n,f), 232U(n,f) PRELIMINARY / PARTIAL RESULTS

  29. 239Pu(n,f) POINT BY POINT MODEL CALCULATION • Experimental FF distributions Y(A,TKE), from EXFOR: • more data sets at En=th: IRMM (Hambsch), Wagemans 1984/2010, • Surin 1971 CCPFEI, Akimov 1971 CCPFEI, Tsuchiya 2000 JPBKTO • Y(A) Akimov et al.: En=0.72 MeV, 1.72 MeV, 2.72 MeV, 4.48 MeV • TKE(A) obtained from experimental <TKE> as a function of En and • renormalization to the shape of TKE(A) at thermal En • CN cross-section of the inverse process: SCAT2 code with more • optical model potential parameterizations: Becchetti-Greenless, • Wilmore-Hodgson, Koning-Delaroche • TXE partition, two methods: νH/(νL+ νH) parameterization (Tudora, 2006) • Iterative procedure of Lemaire et al., 2005

  30. 239Pu(nth,f) multi-parametric representation of the PbP calculated matrix ν(A,TKE) A.Tudora, F.-J.Hambsch, Ann.Nucl.Energy (2010) in press

  31. A.Tudora, F.-J.Hambsch, Ann.Nucl.Energy (2010) in press

  32. PbP calculations at En=th: all experimental FF distributions were used, best agreement of PbP results with multiplicity and spectrum experimental data were obtained for FF distrib of IRMM and Surin. (Tudora, Ann.Nucl.Energy 37 (2010) 43-51) • PbP calculation at other En: • σc(ε) Becchetti-Greenless parameterization • TXE partition – the iterative procedure Lemaire et al., 2005

  33. 239Pu(n,f) average model parameter values obtained from PbP treatment

  34. 239Pu(n,f) MOST PROBABLE FRAGMENTATION APPROACH Using average model parameter values obtained by PbP treatment

  35. Iwamoto spectrum shape close to the case σc=const

More Related