Science and Technology of Nuclear Installations

Science and Technology of Nuclear Installations / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 584256 | 8 pages | https://doi.org/10.1155/2011/584256

Calculation of the Effective Delayed Neutron Fraction by Deterministic and Monte Carlo Methods

Academic Editor: Antonio Alvim
Received07 Mar 2011
Accepted05 May 2011
Published22 Aug 2011

Abstract

The studies on Accelerator-Driven Systems (ADSs) have renewed the interest in the theoretical and computational evaluation of the main integral parameters characterizing subcritical systems (e.g., reactivity, effective delayed neutron fraction š›½eff, and mean prompt neutron generation time). In particular, some kinetic parameters, as the effective delayed neutron fraction, are evaluated in Monte Carlo codes by formulations which do not require the calculation of the adjoint flux. This paper is focused on a theoretical and computational analysis about how the different š›½eff definitions are connected and which are the approximations inherent to the Monte Carlo definition with respect to the standard definition involving weighted integrals. By means of a refined transport computational analysis carried out in a coherent and consistent way, that is, using the same deterministic code and neutron data library for the š›½eff evaluation in different ways, the theoretical analysis is numerically confirmed. Both theoretical and numerical results confirm the effectiveness of the Monte Carlo š›½eff evaluation, at least in cases where spectral differences between total and prompt fluxes are negligible with respect to the value of the functionals entering the classical š›½eff formulation.

1. Introduction

The studies on Accelerator-Driven Systems (ADSs) have renewed the interest in the theoretical and computational evaluation of the main integral parameters characterizing subcritical systems (e.g., reactivity, effective delayed neutron fraction, and mean neutron generation time [1]). In particular, the extensive use of Monte Carlo codes for the analysis of the ADS neutronic behaviour is challenging the deterministic codes, used in the past for the analysis of critical systems, for what concerns the capability to reproduce Monte Carlo results, usually taken as reference [2ā€“4]. Nevertheless, some particular parameters, as the effective delayed neutron fraction, are evaluated in Monte Carlo codes by formulations which do not require the calculation of the adjoint flux. The assessment of the various formulations of the effective delayed neutron fraction is crucial for the system evaluation, since it plays an important role in determining its dynamic characteristics [5]. This paper is focused on a theoretical and computational analysis on the connections among the different š›½eļ¬€ definitions, with special attention to the approximations inherent to the Monte Carlo definition with respect to the standard definition involving weighted integrals. Theoretical results show how the Monte Carlo formulation of š›½eļ¬€ may be related to the classical definition, interpreting the classical one through a reactivity evaluation based on an ā€œimprovedā€ first-order approach of perturbation theory. The computational analysis is carried out in a coherent and consistent way, using the same deterministic code and neutron data library for the š›½eļ¬€ evaluations. A simplified system allows to investigate the features of the various procedures and to use the results to obtain a physical insight. The GUINEVERE system is then selected as a relevant test case for ADS technology. The GUINEVERE experience [6, 7], mainly devoted to the issues concerning online reactivity monitoring in ADS, is analysed by using a modified layout of the VENUS critical facility located at the SCKā€¢CEN site in Mol, Belgium, coupling the subcritical core facility to a deuteron accelerator delivering 14ā€‰MeV neutrons by deuterium-tritium reactions, by a continuous or pulsed beam. For the GUINEVERE experience the ERANOS system [8] and neutron data library JEFF 3.1 [9] have been used to perform transport calculations with 49 energy groups for a cylindrical schematization of the GUINEVERE start-up (at critical) configuration. The analysis of the results allows to draw some conclusions on the merits and the limits of the various formulations.

2. Theoretical Analysis

2.1. The General Case

Let us take as reference system the following eigenvalue problem:š‹šœ‘=šœ”š…šœ‘,(1) where š‹ is the loss operator, š… the fission operator, and šœ”=1/š‘˜, inverse of the classical multiplication eigenvalue. Then, a perturbation š›æš… is introduced in the system. As a consequence, the eigenvalue problem for the perturbed system isš‹šœ‘āˆ—=[]šœ‘šœ”+š›æšœ”][š…+š›æš…āˆ—,(2) with šœ‘āˆ— as the perturbed flux. Following the exact approach of the Perturbation Theory (PT) [10], that is, considering the adjoint problem corresponding to the reference systemš‹+šœ‘+=šœ”š…+šœ‘+,(3) one obtains (the symbol āŸØā‹ÆāŸ© denotes integration over the full phase space, while, in the following, āŸØā‹ÆāŸ©š‘„ is denoting integration only over the specified variable š‘„):š›æšœ”ī«šœ‘šœ”+š›æšœ”=āˆ’+š›æš…šœ‘āˆ—ī¬āŸØšœ‘+š…šœ‘āˆ—āŸ©(4) or equivalentlyš›æšœ”šœ”ī«šœ‘=āˆ’+š›æš…šœ‘āˆ—ī¬āŸØšœ‘+(š…+š›æš…)šœ‘āˆ—āŸ©.(5) If šœ”+š›æšœ”=1/š‘˜āˆ—, we can write (4) in terms of the š‘˜ eigenvalue:š‘˜1āˆ’āˆ—š‘˜ī«šœ‘=āˆ’+š›æš…šœ‘āˆ—ī¬āŸØšœ‘+š…šœ‘āˆ—āŸ©.(6) If we write (2) asš‹[]=[]šœ‘+š›æšœ‘šœ”+š›æšœ”][š…+š›æš…][šœ‘+š›æšœ‘(7) and follow the first-order PT approach, while retaining the second-order term š›æšœ”š›æš…šœ‘, we obtain the ā€œimprovedā€ first-order formulationš›æšœ”ī«šœ‘šœ”+š›æšœ”ā‰…āˆ’+ī¬š›æš…šœ‘āŸØšœ‘+š…šœ‘āŸ©(8) or, in terms of the š‘˜ eigenvalue,š‘˜1āˆ’āˆ—š‘˜ī«šœ‘ā‰…āˆ’+ī¬š›æš…šœ‘āŸØšœ‘+š…šœ‘āŸ©.(9)

2.2. The ā€œPromptā€ Perturbation

If the perturbation is assumed to be equal to the delayed neutrons fission operator š…š‘‘, hence š›æš…=āˆ’š…š‘‘, the perturbed system given in (2) can be written asš‹šœ‘š‘=šœ”š‘ī€ŗš…āˆ’š…šī€»šœ‘š‘,(10) with the perturbed flux given by the prompt flux šœ‘š‘ and šœ”š‘=1/š‘˜š‘. Then (6) becomesš‘˜1āˆ’š‘š‘˜=ī«šœ‘+š…ššœ‘š‘ī¬ī«šœ‘+š…šœ‘š‘ī¬(11) and (9) takes the following form:š‘˜1āˆ’š‘š‘˜ā‰…ī«šœ‘+š…ššœ‘ī¬āŸØšœ‘+š…šœ‘āŸ©=š›½eļ¬€.(12) Thus, the effective delayed neutron fraction š›½eļ¬€ is an ā€œimprovedā€ PT first-order formulation of the relationship 1āˆ’(š‘˜š‘/š‘˜), widely used in Monte Carlo codes as š›½eļ¬€ estimator. The pure PT first-order formulation provides:š›æšœ”šœ”=š‘˜š‘˜š‘ī«šœ‘āˆ’1ā‰…āˆ’+ī¬š›æš…šœ‘āŸØšœ‘+=ī«šœ‘š…šœ‘āŸ©+š…ššœ‘ī¬āŸØšœ‘+š…šœ‘āŸ©=š›½eļ¬€.(13) To pass from (11) to (12), it is sufficient to replace šœ‘š‘ by šœ‘.

3. Standard and Monte Carlo š›½ Formulations

The effectiveness of delayed neutrons is normally obtained by one of the following two relationships:š›½eļ¬€=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–Ī¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)š‘‘,š‘–ī‚¬Ī£(šø)š‘“(š‘š)ī‚­(š«,šøā€²)Ī¦(š«,šøā€²)šøī…žī‚­šøī‚­š«āˆ‘š‘šĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)(šøā€²)Ī£š‘“(š‘š)ī‚­(š«,šøā€²)Ī¦(š«,šøā€²)šøī…žī‚­šøī‚­š«,(14)Ģƒš›½eļ¬€=ī“š‘šī“š‘–š›½š‘–(š‘š)Ī¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)š‘‘,š‘–ī‚¬šœˆ(šø)(š‘š)(šøā€²)Ī£š‘“(š‘š)ī‚­(š«,šøā€²)Ī¦(š«,šøā€²)šøī…žī‚­šøī‚­š«āˆ‘š‘šĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)(šøā€²)Ī£š‘“(š‘š)ī‚­(š«,šøā€²)Ī¦(š«,šøā€²)šøī…žī‚­šøī‚­š«,(15)where š‘š is the fissile isotope index, š‘– the delayed neutron family index,šœ’(š‘š)š‘‘,š‘– the delayed neutron spectrum for fissile isotope š‘š and delayed neutron family š‘– (basic data),šœ’(š‘š) the total neutron spectrum for fissile isotope š‘š,šœˆ(š‘š)š‘‘,š‘– the average value of delayed neutrons emitted from fissile isotope š‘š and delayed neutron family š‘–, for a given incident neutron spectrum (from basic data),šœˆ(š‘š) the total neutrons emitted from fissile isotope š‘š, and Ī£š‘“(š‘š) the macroscopic fission cross section for fissile isotope š‘š.

Furthermore,š›½š‘–(š‘š)=šœˆ(š‘š)š‘‘,š‘–šœˆšœ“(š‘š),(16) with šœˆšœ“(š‘š) being the average value of the number of total neutrons emitted from the fissile isotope š‘š for a given incident neutron spectrum šœ“ (šœ“ has to be calculated by a cell or system neutronic calculation). In (14) and (15) the direct flux and the adjoint flux are solutions for the following eigenvalue equations, respectively:š›€ā‹…āˆ‡šœ‘(š«,š›€,šø)+Ī£š‘”=1(š«,šø)šœ‘(š«,š›€,šø)š‘˜ī“š‘ššœ’(š‘š)(šø)<šœˆ(š‘š)(šøā€²)Ī£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦ī€·š«,šøī…žī€ø>šøā€²+<Ī£š‘ ī€·š«,šøī…žī€øĪ¦ī€·āŸ¶šøš«,šøī…žī€ø>šøī…ž(17)āˆ’š›€ā‹…āˆ‡šœ‘+(š«,š›€,šø)+Ī£š‘”(š«,šø)šœ‘+=1(š«,š›€,šø)š‘˜ī“š‘ššœˆ(š‘š)(šø)Ī£š‘“(š‘š)ī«šœ’(š«,šø)(š‘š)ī€·šøī…žī€øĪ¦+ī€·š«,šøī…žī€øī¬šøā€²+ī«Ī£š‘ ī€·š«,šøāŸ¶šøī…žī€øĪ¦+ī€·š«,šøī…žī€øī¬šøā€².(18) Current Monte Carlo (MC) calculations approximate š›½eļ¬€ by the following formula:š›½MCeļ¬€š‘˜=1āˆ’š‘š‘˜,(19) where š‘˜š‘ is the multiplication eigenvalue of the equationš›€ā‹…āˆ‡šœ‘š‘(š«,š›€,šø)+Ī£š‘”(š«,šø)šœ‘š‘=1(š«,š›€,šø)š‘˜š‘ī“š‘šīƒ¬šœ’(š‘š)ī“(šø)āˆ’š‘–šœ’(š‘š)š‘‘,š‘–(šø)š›½š‘–(š‘š)īƒ­Ć—ī‚¬šœˆ(š«)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²+ī«Ī£š‘ ī€·š«,šøī…žī€øĪ¦āŸ¶šøš‘ī€·š«,šøī…žī€øī¬šøā€².(20) Equation (20) is written in a general form, taking into account the spatial dependence of the delayed emission fraction š›½ on position š«. In this case š›½š‘–(š‘š)(š«) is given byš›½š‘–(š‘š)(š«)=šœˆ(š‘š)š‘‘,š‘–šœˆī‚ƒī‚¬(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²Ī£ī‚‚ī‚¬š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚„.(21)If we consider the system governed by (20) as a perturbation of the system governed by (17), we obtain as exact perturbation (11), that is,š‘˜1āˆ’š‘š‘˜=ī“š‘šī“š‘–Ī¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)š‘‘,š‘–(šø)š›½š‘–(š‘š)ī‚¬šœˆ(š«)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚­šøī‚­š«āˆ‘š‘šĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚­šøī‚­š«,(22)where šœ‘+ is the eigenfunction of (18) and šœ‘š‘ the eigenfunction of (20). When comparing (22) with (15), it can be noticed that, besides the replacement of šœ‘ by šœ‘š‘, the differences lie in the presence of the terms š›½š‘–(š‘š)(š«) (under integration) replacing the terms š›½š‘–(š‘š). If we assume in both (20) and (22) the following approximation:š›½š‘–(š‘š)(š«)ā‰ˆš›½š‘–(š‘š)=šœˆ(š‘š)š‘‘,š‘–šœˆšœ“(š‘š),(23) equation (20) takes the following form:š›€ā‹…āˆ‡šœ‘š‘(š«,š›€,šø)+Ī£š‘”(š«,šø)šœ‘š‘=1(š«,š›€,šø)Ģƒš‘˜š‘ī“š‘šīƒ¬šœ’(š‘š)ī“(šø)āˆ’š‘–šœ’(š‘š)š‘‘,š‘–(šø)š›½š‘–(š‘š)īƒ­Ć—ī‚¬šœˆ(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²+ī«Ī£š‘ ī€·š«,šøī…žī€øĪ¦āŸ¶šøš‘ī€·š«,šøī…žī€øī¬šøā€²(24) and (22) becomesĢƒš›½MCeļ¬€Ģƒš‘˜=1āˆ’š‘š‘˜=ī“š‘šī“š‘–š›½š‘–(š‘š)ƗĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)š‘‘,š‘–ī‚¬šœˆ(šø)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøī…žī‚­šøī‚­š«āˆ‘š‘šĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚­šøī‚­š«.(25)

As shown in the previous section, the classical š›½eļ¬€ definition given in (15) is an ā€œimprovedā€ PT first-order approximation of the relationship 1āˆ’(š‘˜š‘/š‘˜) given in (25). Along the point of view of the perturbation approach, leading to a relationship among system integral properties like the eigenvalues and the delayed neutron effectiveness, (24) and (25) are intimately coupled, in the sense that the coherent prompt flux šœ‘š‘ to be used in the š›½eļ¬€ definition given by (25) is the eigenfunction of (24). We recall that (24) is an approximation of (20), owing to the assumption given in (23). If we want to remove this assumption, we have to make explicit the š›½(š«) definition given in (21) and insert it into (20), thus obtainingš›€ā‹…āˆ‡šœ‘š‘(š«,š›€,šø)+Ī£š‘”(š«,šø)šœ‘š‘=1(š«,š›€,šø)š‘˜š‘ī“š‘šīƒ¬šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²āˆ’ī“š‘–šœ’(š‘š)š‘‘,š‘–(šø)šœˆ(š‘š)š‘‘,š‘–ī‚¬Ī£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²īƒ­+ī«Ī£š‘ ī€·š«,šøī…žī€øĪ¦āŸ¶šøš‘ī€·š«,šøī…žī€øī¬šøā€².(26) Following the same perturbation approach, that is, perturbed system given by (26) and unperturbed system given by (17), we findš›½MCeļ¬€š‘˜=1āˆ’š‘š‘˜=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–Ī¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)š‘‘,š‘–ī‚¬Ī£(šø)š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚­šøī‚­š«āˆ‘š‘šĪ¦ī‚¬ī‚¬+(š«,šø)šœ’(š‘š)ī‚¬šœˆ(šø)(š‘š)ī€·šøī…žī€øĪ£š‘“(š‘š)ī€·š«,šøī…žī€øĪ¦š‘ī€·š«,šøī…žī€øī‚­šøā€²ī‚­šøī‚­š«.(27)

Analogously, the classical š›½eļ¬€definition given in (14) is an ā€œimprovedā€ PT first-order approximation of the relationship 1āˆ’(š‘˜š‘/š‘˜) given in (27). Also in this case (26) and (27) are intimately coupled, in the sense that the coherent prompt flux šœ‘š‘ to be used in the š›½eļ¬€ definition given by (27) is the eigenfunction of (26). The previous elaboration can be summarized in the following way:(a)as introduced in the previous section, the classical š›½eļ¬€ definitions, like those given in (14) and (15), which do not require any prompt flux calculation, are ā€œimprovedā€ PT first-order approximations of the relationship 1āˆ’(š‘˜š‘/š‘˜) used by Monte Carlo calculations;(b)the relationship 1āˆ’(š‘˜š‘/š‘˜) provides a š›½eļ¬€ estimate coherent with the assumptions made when performing the prompt flux calculation, that is, (24) provides the estimate given by (25), and (26) provides the estimate given by (27).

4. ERANOS Formalism and Calculation Set up

The following š›½eļ¬€ formulation is adopted in ERANOS [11] (given in the energy multigroup scheme, with g as energy group index):š›½eļ¬€=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘‘,š‘–,š‘”āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²ī‚­(š«)š«ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘”āˆ‘š‘šāˆ‘š‘”ā€²šœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²ī‚­(š«)š«,(28) where šœ’š‘‘,š‘–,š‘” is a unique delayed neutron spectrum (derived from basic data) for each fissile isotope š‘š, delayed neutron family i, andšœ’š‘”=āˆ‘š‘ššœ’š‘”(š‘š)āˆ‘š‘”ā€²šœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€²āˆ‘š‘šāˆ‘š‘”ā€²šœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€².(29) The quantity šœ’š‘” is a unique fission spectrum, used by ERANOS for spatial calculations, for all the fissile isotopes. The weighting flux šœ“š‘” is obtained by the fuel cell cross section calculation. In (28) the direct flux and the adjoint flux are solutions for the following eigenvalue equations, respectively:š›€ā‹…āˆ‡šœ‘š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘š‘”=1(š«,š›€)š‘˜šœ’š‘”ī“š‘šī“š‘”ī…žšœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²ī“(š«)+š‘”ī…žĪ£š‘ ,š‘”ā€²ā†’š‘”(š«)Ī¦š‘”ā€²(š«);(30)āˆ’š›€ā‹…āˆ‡šœ‘+š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘+š‘”=1(š«,š›€)š‘˜ī“š‘ššœˆš‘”(š‘š)Ī£(š‘š)š‘“,š‘”ī“(š«)š‘”ā€²šœ’š‘”ā€²Ī¦+š‘”ā€²ī“(š«)+š‘”ā€²Ī£š‘ ,š‘”ā†’š‘”ā€²(š«)Ī¦+š‘”ī…ž(š«).(31)

It can be seen that, besides the approximation on the delayed neutron spectrum and on the total fission spectrum, (28) is the energy-multigroup version of (14). Considering as perturbed state with respect to the reference state in (30) the solution for the following equation:š›€ā‹…āˆ‡šœ‘š‘,š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘š‘,š‘”=1(š«,š›€)š‘˜š‘ī“š‘šī“š‘”ā€²īƒ¬šœ’š‘”šœˆš‘”(š‘š)ā€²āˆ’ī“š‘–šœ’š‘‘,š‘–,š‘”šœˆ(š‘š)š‘‘,š‘–īƒ­Ī£(š‘š)š‘“,š‘”ī…ž(š«)Ī¦š‘,š‘”ā€²(+ī“š«)š‘”ī…žĪ£š‘ ,š‘”ā€²ā†’š‘”(š«)Ī¦š‘,š‘”ā€²(š«),(32) we obtain as exact perturbationš›½MCeļ¬€š‘˜=1āˆ’š‘š‘˜=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–Ć—ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘‘,š‘–,š‘”āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²ī‚­(š«)š«ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘”āˆ‘š‘šāˆ‘š‘”ā€²šœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²ī‚­(š«)š«.(33)

On the basis of the theory outlined in the previous section, (28) is the ā€œimprovedā€ PT first-order approximation of (33). Unfortunately, ERANOS cannot calculate the solution for (32) and consequently (33) cannot be applied because, as mentioned before, it uses a unique fission spectrum for all the fissile isotopes (and in (32) both šœ’š‘” and the various šœ’š‘‘,š‘–,š‘” appear). In addition, to use (for our scopes) the š›½eļ¬€ modules built-in in ERANOS we need a further assumption, that is, we have to evaluatešœˆ(š‘š)=āˆ‘š‘”ā€²šœˆš‘”(š‘š)ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€²āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ī…žšœ“š‘”ā€²,(34) with šœ“š‘” as weighting flux obtained by the cell cross-section calculation. Thus, our new reference state is given by the following equation:š›€ā‹…āˆ‡šœ‘š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘š‘”=1(š«,š›€)š‘˜šœ’š‘”ī“š‘ššœˆ(š‘š)ī“š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²+ī“(š«)š‘”ā€²Ī£š‘ ,š‘”ā€²ā†’š‘”(š«)Ī¦š‘”ā€²(š«),(35) with the associated adjointāˆ’š›€ā‹…āˆ‡šœ‘+š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘+š‘”=1(š«,š›€)š‘˜ī“š‘ššœˆ(š‘š)Ī£(š‘š)š‘“,š‘”ī“(š«)š‘”ā€²šœ’š‘”ā€²Ī¦+š‘”ā€²ī“(š«)+š‘”ā€²Ī£š‘ ,š‘”ā†’š‘”ā€²(r)Ī¦+š‘”ā€²(š«),(36) where (of course, because of (34), (37) provides the same result as (29)):šœ’š‘”=āˆ‘š‘ššœ’š‘”(š‘š)šœˆ(š‘š)āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€²āˆ‘š‘ššœˆ(š‘š)āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€².(37) Then, settingšœ’š‘,š‘”=āˆ‘š‘šī‚ƒšœ’š‘”šœˆ(š‘š)āˆ’āˆ‘š‘–šœ’š‘‘,š‘–,š‘”šœˆ(š‘š)š‘‘,š‘–ī‚„āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€²āˆ‘š‘šī‚ƒšœˆ(š‘š)āˆ’āˆ‘š‘–šœˆ(š‘š)š‘‘,š‘–ī‚„āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²šœ“š‘”ā€²,(38) Equation (32) becomesš›€ā‹…āˆ‡šœ‘š‘,š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘š‘,š‘”=1(š«,š›€)š‘˜š‘šœ’š‘,š‘”ī“š‘šīƒ¬šœˆ(š‘š)āˆ’ī“š‘–šœˆ(š‘š)š‘‘,š‘–īƒ­ī“š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²(+ī“š«)š‘”ā€²Ī£š‘ ,š‘”ā€²ā†’š‘”(š«)Ī¦š‘,š‘”ī…ž(š«).(39)

On the basis of (16), (39) is equivalent to:š›€ā‹…āˆ‡šœ‘š‘,š‘”(š«,š›€)+Ī£š‘”,š‘”(š«)šœ‘š‘,š‘”=1(š«,š›€)š‘˜š‘šœ’š‘,š‘”ī“š‘šīƒ¬ī“1āˆ’š‘–š›½š‘–(š‘š)īƒ­šœˆ(š‘š)ī“š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²(+ī“š«)š‘”ā€²Ī£š‘ ,š‘”ā€²ā†’š‘”(š«)Ī¦š‘,š‘”ā€²(š«)(40) Considering now (39) as the perturbed state with respect to the reference state in (35), taking into account (38), we obtainš›½MCeļ¬€š‘˜=1āˆ’š‘š‘˜=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–Ć—ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘‘,š‘–,š‘”āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²ī‚­(š«)š«ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘”āˆ‘š‘ššœˆ(š‘š)āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘,š‘”ā€²ī‚­(š«)š«.(41) The corresponding ā€œimprovedā€ PT first-order approximation of (41) isš›½eļ¬€=ī“š‘šī“š‘–šœˆ(š‘š)š‘‘,š‘–ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘‘,š‘–,š‘”āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²ī‚­(š«)š«ī‚¬āˆ‘š‘”Ī¦+š‘”(š«)šœ’š‘”āˆ‘š‘ššœˆ(š‘š)āˆ‘š‘”ā€²Ī£(š‘š)š‘“,š‘”ā€²(š«)Ī¦š‘”ā€²ī‚­(š«)š«.(42) and Equation (39), both sides of (41) and (42) can be calculated by ERANOS.

5. ERANOS Results

The GUINEVERE system, shown in Figure 1, is selected to perform the numerical analysis. The fuel fissile zone is made of Uranium 30% (weight) enriched in U235. The fuel matrix has the following volumetric fractions: Uranium 17%, Stainless Steel 16%, Lead 60% and Air 7%. The zone Air + SS Sheet has the following volumetric fractions: Stainless Steel 93% and Air 7%. The RZ system is a simplified schematization of the GUINEVERE start-up (at critical) configuration (after start-up the central zone is replaced by the beam tube of the deuteron accelerator delivering, in continuous or pulsed mode, 14ā€‰MeV neutrons by deuterium-tritium reactions).

ERANOS transport calculations are performed with 49 energy groups (see the appendix), P0 transport approximation, and angular quadrature S4 (module BISTRO [8]). The neutron data library JEFF 3.1 [9] is used for the calculations. In Table 1 the delayed neutron data šœˆ(š‘š)š‘‘,š‘– [12] used for the analysis are shown. Delayed neutron spectra šœ’š‘‘,š‘–,š‘” are the ones built-in in ERANOS and taken from [12].


Family12345678Sum

U2350.000550.002450.001620.003260.005090.001520.001420.000390.01630
U2380.000390.004840.001740.006370.013670.009210.005950.004330.04650

The following šœˆ(š‘š) values were obtained by (28): šœˆ(235)=2.49 and šœˆ(238)=2.75. This high šœˆ(238) value is due to the fact that about 61% of U(238) total fissions occur in the energy range 2 Ć· 8ā€‰MeV, where š‘£(238) values span from 2.67 to 3.48 (only about 5% of U(238) total fissions occur above 2ā€‰MeV). In Figure 2 a comparison between the total fission spectrum šœ’š‘” obtained by (37) and the prompt fission spectrum šœ’š‘,š‘” obtained by (38) is shown. In Table 2 a synthesis of the results is reported.


CaseEquationkkp š›½ e ļ¬€ (pcm)1āˆ’kp/k (pcm)

(a)(28) for š›½ e ļ¬€ 723
(b)(31) for Ļ†1.01955
(c)(31) for Ļ† +1.01955
(d)(35) for Ļ†1.02149
(e)(36) for Ļ† +1.02149
(f)(39) for Ļ† p1.01402
(g)(41) for 1āˆ’kp/k 732
(h)(41) for š›½ e ļ¬€ (using Ļ† p)732
(i)(42) for š›½ e ļ¬€ (using Ļ†)731


Upper energy (eV)Group #Upper energy (eV)Group #Upper energy (eV)Group #Upper energy (eV)Group #

1.96E + 0712.470E + 06141.11E + 05271.49E + 0240
1.73E + 0722.23E + 06156.74E + 04289.17E + 0141
1.49E + 0732.02E + 06164.09E + 04296.79E + 0142
1.38E + 0741.65E + 06172.48E + 04304.02E + 0143
1.16E + 0751.35E + 06181.50E + 04312.26E + 0144
1.00E + 0761.22E + 06199.12E + 03321.37E + 0145
8.19E + 0671.11E + 06205.53E + 03338.32E + 0046
6.70E + 0681.00 + 06213.35E + 03344.000E + 0047
6.07E + 0699.07E + 05222.03E + 03355.40E āˆ’ 0148
5.49E + 06108.21E + 05231.23E + 03361.00E āˆ’ 0149
4.49E + 06114.98E + 05247.49E + 0237
3.68E + 06123.02E + 05254.54E + 0238
3.04E + 06131.83E + 05263.04E + 0239

Case (a) is the classical š›½eļ¬€ formulation adopted in ERANOS, with šœ‘ and šœ‘+ from cases (b) and (c), respectively. Cases (d) and (e) differ from (b) and (c) for the presence of šœˆ(š‘š) in the fission production terms, and this approximation provides a difference of +186ā€‰pcm with respect to the standard treatment of the fission production terms. The prompt case (f), again with the presence of šœˆ(š‘š) in the fission production terms, provides a difference of āˆ’721ā€‰pcm (~š›½eļ¬€) with respect to cases (d) and (e).

Case (g), representing the approach used by Monte Carlo calculations, provides exactly the same value as in case (h), that is, the last term in (41). This numerical result confirms the theoretical analysis described above, linking the Monte Carlo formulation to the PT approach. Finally, case (i) shows that the replacement of the prompt flux in case (h) by the total flux does not affect appreciably the š›½eļ¬€ value (at least for the fast system considered).

Actually, the differences between šœ‘ and šœ‘š‘, from cases (d) and (f) in Table 2, respectively, are not remarkable, as it can be seen from Figure 3 (where šœ‘ and šœ‘š‘ are undistinguishable). šœ‘ and šœ‘š‘ in Figure 3 are spatial average values over the core volume, same area under šœ‘ and šœ‘š‘ curves.

6. Conclusions

When evaluating š›½eļ¬€ by the relationship 1āˆ’(š‘˜š‘/š‘˜), the quality of the obtained results depends on the quality of the description of the delayed neutron emissions assumed in the š‘˜š‘ prompt calculation. It could seem obvious, but if we look at the standard route used in ERANOS (and in general in deterministic codes) to calculate š›½eļ¬€, by a formulation which does not require any prompt flux calculation, we realize that great detail is given to the characteristics of the delayed neutron emission through the š›½eļ¬€ formulation, and without the need to transpose such a detail into a full system prompt neutron calculation. For example, in the ERANOS case, it is not possible to set up a prompt calculation having the same quality of the delayed neutron emission description available in the š›½eļ¬€ formulation. And, in fact, in the present case we are obliged to lower the level of the information about the delayed neutron emission in the š›½eļ¬€ definition in order to run the corresponding prompt calculation.

Following the perturbation theory approach, a rigorous relationship may be established between the š›½eļ¬€ evaluation by the formula 1āˆ’(š‘˜š‘/š‘˜), used in Monte Carlo codes, and the corresponding š›½eļ¬€ calculation by the classical formulation involving direct and adjoint fluxes, and it has been shown how the classical formulation can be considered an ā€œimprovedā€ PT first-order formulation of the formula 1āˆ’(š‘˜š‘/š‘˜). By means of a transport computational analysis using a 49 energy group structure and carried out in a coherent and consistent way, that is, using the same deterministic code ERANOS and neutron data library JEFF 3.1 for the š›½eļ¬€ evaluation in different ways, the theoretical analysis is numerically confirmed.

Both theoretical and numerical results confirm the effectiveness of the š›½eļ¬€ evaluation by the relationship 1āˆ’(š‘˜š‘/š‘˜), at least in cases where spectral differences between total and prompt fluxes are negligible with respect to the value of the functionals entering the classical š›½eļ¬€ formulation. For other material configurations more investigations are needed to obtain a detailed quantification of the effects involved.

Appendix

The 49 energy group structure, presented in Table 3, has been adopted in order to capture the various spectral phenomena associated to the source-driven subcritical assembly. Different built-in energy group grids are available into ERANOS. In this work the base library at 1968 energy groups, together with two derived energy grids at 172 and 33 energy groups, has been used to produce the cross sections. The 49 energy group structure has been obtained using the standard 33 energy group below 0.82ā€‰MeV and the finest 172 energy group above this energy in order to better take into account the D-T neutron external source foreseen for the GUINEVERE subcritical configurations (neutron energy at about 14ā€‰MeV).

References

  1. S. Dulla, P. Picca, D. Tomatis, P. Ravetto, and M. Carta, ā€œIntegral parameters in source-driven systems,ā€ Progress in Nuclear Energy, vol. 53, no. 1, pp. 32ā€“40, 2011. View at: Publisher Site | Google Scholar
  2. R. K. Meulekamp and S. C. van der Marck, ā€œCalculating the effective delayed neutron fraction with Monte Carlo,ā€ Nuclear Science and Engineering, vol. 152, no. 2, pp. 142ā€“148, 2006. View at: Google Scholar
  3. B. Verboomen, W. Haeck, and P. Baeten, ā€œMonte Carlo calculation of the effective neutron generation time,ā€ Annals of Nuclear Energy, vol. 33, no. 10, pp. 911ā€“916, 2006. View at: Publisher Site | Google Scholar
  4. Y. Nagaya and T. Mori, ā€œCalculation of effective delayed neutron fraction with Monte Carlo perturbation techniques,ā€ Annals of Nuclear Energy, vol. 38, no. 2-3, pp. 254ā€“260, 2011. View at: Publisher Site | Google Scholar
  5. Y. Nagaya, G. Chiba, T. Moria, D. Irwanto, and K. Nakajima, ā€œComparison of Monte Carlo calculation methods for effective delayed neutron fraction,ā€ Annals of Nuclear Energy, vol. 37, no. 10, pp. 1308ā€“1315, 2010. View at: Publisher Site | Google Scholar
  6. P. Baeten, H. AĆÆt Abderrahim, T. Aoust et al., ā€œThe GUINEVERE project at the VENUS facility,ā€ in Proceedings of the International Conference on the Physics of Reactors (PHYSOR '08), vol. 4, pp. 2867ā€“2873, Interlaken, Switzerland, September 2008. View at: Google Scholar
  7. http://www.sckcen.be/en/Our-Research/Research-projects/EU-projects-FP6-FP7/GUINEVERE.
  8. G. Rimpault et al., ā€œThe ERANOS code and data system for fast reactor neutronic analyses,ā€ in Proceedings of the International Conference on the Physics of Reactors, (PHYSOR '02), Seoul, Korea, October 2002. View at: Google Scholar
  9. A. Santamarina, D. Bernard, P. Blaise et al., ā€œThe JEFF-3.1.1 Nuclear Data Libraryā€”Report 22,ā€ NEA 6807, Nuclear Energy Agency, 2009. View at: Google Scholar
  10. A. Gandini, ā€œA generalized perturbation method for bi-linear functionals of the real and adjoint neutron fluxes,ā€ Journal of Nuclear Energy, vol. 21, no. 10, pp. 755ā€“765, 1967. View at: Google Scholar
  11. V. Zammit and E. Fort, ā€œConstantes de Neutrons RetardĆ©s pour le calcul du Ī²eff avec ERANOS,ā€ CEA Technical Note NT-SPRC-LEPh-97-230, 1997. View at: Google Scholar
  12. A. D'Angelo and J. L. Rowlands, ā€œConclusions concerning the delayed neutron data for the major actinides,ā€ Progress in Nuclear Energy, vol. 41, no. 1ā€“4, pp. 391ā€“412, 2002. View at: Publisher Site | Google Scholar

Copyright © 2011 M. Carta et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

1944Ā Views | 690Ā Downloads | 6Ā Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.