Enhanced Fission Probability of Even-Z Fragments in the Decay of Hot and Rotating 210Rn* Compound System

Mass and charge distribution of the cross-section for the fission fragments obtained in the decay of hot and rotating compound system formed in the reaction 48Ca + 162Dy → 210Rn* at an incident energy 139.6 MeV has been calculated using the dynamical cluster-decay model. Isotopic composition for each element belonging to the symmetric mass region has been obtained. The shell closure at N=50 for light and at Z=50 for heavy mass binary fragments gives a deep minima in the fragmentation potential at touching configuration and governs the fission partition of the compound system. The fission fragments of the symmetric mass region have their dominating presence along with strong odd-even staggering i.e., even-Z fission fragments are more probable than the odd ones, similar to the observed trends of the yield.

Comparison of Direct and Adjoint k and α-Eigenfunctions

Modal expansions based on k-eigenvalues and α-eigenvalues are commonly used in order to investigate the reactor behaviour, each with a distinct point of view: the former is related to fission generations, whereas the latter is related to time. Well-known Monte Carlo methods exist to compute the direct k or α fundamental eigenmodes, based on variants of the power iteration. The possibility of computing adjoint eigenfunctions in continuous-energy transport has been recently implemented and tested in the development version of TRIPOLI-4®, using a modified version of the Iterated Fission Probability (IFP) method for the adjoint α calculation. In this work we present a preliminary comparison of direct and adjoint k and α eigenmodes by Monte Carlo methods, for small deviations from criticality. When the reactor is exactly critical, i.e., for k0 = 1 or equivalently α0 = 0, the fundamental modes of both eigenfunction bases coincide, as expected on physical grounds. However, for non-critical systems the fundamental k and α eigenmodes show significant discrepancies.

ANALYSIS OF TIME-EIGENVALUE AND EIGENFUNCTIONS IN THE CROCUS BENCHMARK

Time-dependent neutron transport in non-critical state can be expressed by the natural mode equation. In order to estimate the dominant eigenvalue and eigenfunction of the natural mode, CEA had extended the α-k method and developed the generalized iterated fission probability method (G-IFP) in the TRIPOLI-4® code. CRIEPI has chosen to compute those quantities by a time-dependent neutron transport calculation, and has thus developed a time-dependent neutron transport technique based on k-power iteration (TDPI) in MCNP-5. In this work, we compare the two approaches by computing the dominant eigenvalue and the direct and adjoint eigenfunctions for the CROCUS benchmark. The model has previously been qualified for keffs and kinetic parameters by TRIPOLI-4 and MCNP-5. The eigenvalues of the natural mode equations by α-k and TDPI are in good agreement with each other, and closely follow those predicted by the inhour equation. Neutron spectra and spatial distributions (flux and fission neutron emission) obtained by the two methods are also in good agreement. Similar results are also obtained for the adjoint fundamental eigenfunctions. These findings substantiate the coherence of both calculation strategies for natural mode.

SUPER-HISTORY METHODS FOR ADJOINT-WEIGHTED TALLIES IN MONTE CARLO TIME EIGENVALUE CALCULATIONS

Time eigenvalues emerge in several key applications related to neutron transport problems, including reactor start-up and reactivity measurements. In this context, experimental validation and uncertainty quantification would demand to assess the variation of the dominant time eigenvalue in response to a variation of nuclear data. Recently, we proposed the use of a Generalized Iterated Fission Probability method (G-IFP) to compute adjoint-weighted tallies, such as kinetic parameters, perturbations and sensitivity coefficients, for Monte Carlo time (or alpha) eigenvalue calculations. With the massive use of parallel Monte Carlo calculations, it would be therefore useful to trade the memory burden of the G-IFP method (which is comparable to that of the standard IFP method for k-eigenvalue problems) for computation time and to rely on history-based schemes for such adjoint-weighted tallies. For this purpose, we investigate the use of the super-history method as applied to estimating adjoint-weighted tallies within the α-k power iteration, based on previous work on k-eigenvalue problems. Verification of the algorithms is performed on some simple preliminary tests where analytic solutions exist. In addition, the performances of the proposed method are assessed by comparing the super-history and the G-IFP methods for the same sets of benchmark problems.

Fission fragment mass and total kinetic energy distributions of spontaneously fissioning plutonium isotopes

The fission-fragment mass and total kinetic energy (TKE) distributions are evaluated in a quantum mechanical framework using elongation, mass asymmetry, neck degree of freedom as the relevant collective parameters in the Fourier shape parametrization recently developed by us. The potential energy surfaces (PES) are calculated within the macroscopic-microscopic model based on the Lublin-Strasbourg Drop (LSD), the Yukawa-folded (YF) single-particle potential and a monopole pairing force. The PES are presented and analysed in detail for even-even Plutonium isotopes with A = 236–246. They reveal deep asymmetric valleys. The fission-fragment mass and TKE distributions are obtained from the ground state of a collective Hamiltonian computed within the Born-Oppenheimer approximation, in the WKB approach by introducing a neck-dependent fission probability. The calculated mass and total kinetic energy distributions are found in good agreement with the data.