Calculation of the eigenfunctions of the two-group neutron diffusion equation and application to modal decomposition of BWR instabilities

In this thesis, numerical methods aiming at determining the eigenfunctions, their adjoint and the corresponding eigenvalues of the two-group neutron diffusion equations representing any heterogeneous system are investigated. First, the classical power iteration method is modified so that the calculation of modes higher than the fundamental mode is possible. Thereafter, the Explicitly-Restarted Arnoldi method, belonging to the class of Krylov subspace methods, is touched upon. Although the modified power iteration method is a computationally-expensive algorithm, its main advantage is its robustness, i.e. the method always converges to the desired eigenfunctions without any need from the user to set up any parameter in the algorithm. On the other hand, the Arnoldi method, which requires some parameters to be defined by the user, is a very efficient method for calculating eigenfunctions of large sparse system of equations with a minimum computational effort. These methods are thereafter used for off-line analysis of the stability of Boiling Water Reactors. Since several oscillation modes are usually excited (global and regional oscillations) when unstable conditions are encountered, the characterization of the stability of the reactor using for instance the Decay Ratio as a stability indicator might be difficult if the contribution from each of the modes are not separated from each other. Such a modal decomposition is applied to a stability test performed at the Swedish Ringhals-1 unit in September 2002, after the use of the Arnoldi method for pre-calculating the different eigenmodes of the neutron flux throughout the reactor. The modal decomposition clearly demonstrates the excitation of both the global and regional oscillations. Furthermore, such oscillations are found to be intermittent with a time-varying phase shift between the first and second azimuthal modes.

La metrica di Agmon ed il decadimento esponenziale delle soluzioni di equazioni alle derivate parziali

The aim of this master thesis is to study the exponential decay of solutions of elliptic partial equations. This work is based on the results obtained by Agmon. To this purpose, first, we define the Agmon metric, that plays an important role in the study of exponential decay, because it is related to the rate of decay. Under some assumptions on the growth of the function and on the positivity of the quadratic form associated to the operator, a first result of exponential decay is presented. This result is then applied to show the exponential decay of eigenfunctions with eigenvalues whose real part lies below the bottom of the essential spectrum. Finally, three examples are given: the harmonic oscillator, the hydrogen atom and a Schrödinger operator with purely discrete spectrum.