Using modern methods of reactor physics we have performed full-scale calculations of the natural reactor Oklo. For reliability we have used recent version of two Monte Carlo codes: Russian code MCU REA and world wide known code MCNP (USA). Both codes produce similar results. We have constructed a computer model of the reactor Oklo zone RZ2 which takes into account all details of design and composition. The calculations were performed for three fresh cores with different uranium contents. Multiplication factors, reactivities and neutron fluxes were calculated. We have estimated also the temperature and void effects for the fresh core. As would be expected, we have found for the fresh core a significant difference between reactor and Maxwell spectra, which was used before for averaging cross sections in the Oklo reactor. The averaged cross section of 149 62 Sm and its dependence on the shift of resonance position (due to variation of fundamental constants) are significantly different from previous results.Contrary to results of some previous papers we find no evidence for the change of the fine structure constant in the past and obtain new, most accurate limits on its variation with time: −4 · 10 −17 year −1 ≤α/α ≤ 3 · 10 −17 year −1 . A further improvement in the accuracy of the limits can be achieved by taking account of the core burnup. These calculations are in progress.

No abstract

No abstract

The time-dependent bounce-averaged CQL3D flux-conservative finite-difference Fokker-Planck equation (FPE) solver has been upgraded to include finite-orbit-width (FOW) capabilities which are necessary for an accurate description of neoclassical transport, losses to the walls, and transfer of particles, momentum, and heat to the scrape-off layer. The FOW modifications are implemented in the formulation of the neutral beam source, collision operator, RF quasilinear diffusion operator, and in synthetic particle diagnostics. The collisional neoclassical radial transport appears naturally in the FOW version due to the orbitaveraging of local collision coefficients coupled with transformation coefficients from local (R, Z) coordinates along each guiding-center orbit to the corresponding midplane computational coordinates, where the FPE is solved. In a similar way, the local quasilinear RF diffusion terms give rise to additional radial transport of orbits. We note that the neoclassical results are obtained for 'full' orbits, not dependent on a common small orbit-width approximation. Results of validation tests for the FOW version are also presented.

The behaviour of the fast particle population during 18 keV hydrogen and 26 keV deuterium neutral beam injection in deuterium plasmas is investigated. Experiments reveal large fast ion losses. The experimental results are confirmed using different types of modelling: simulation using the NUBEAM module, solution of the Boltzmann kinetic equation and solution of the 3D fast ion tracking algorithm. The dynamics of the energetic particle redistribution and losses during sawtooth oscillation and toroidal Alfvén eigenmodes are studied. A method to decrease fast ion losses under the current conditions (0.4 T, 0.2 MA) is shown. The influence of the plasma parameters on the energetic ion confinement rate is investigated. Modelling for the Globus-M2 conditions (1 T, 0.5 MA) is performed.