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.

Origin of cold gas and dust in massive elliptical galaxies

The recent availability of multi-wavelength data revealed the presence of large reservoirs of warm and cold gas and dust in the innermost regions of the majority of massive elliptical galaxies. To prove an internal origin of cold and warm gas, the investigation of the spatially distributed cooling process which occurs because of non-linear density perturbations and subsequent thermal instabilities is of crucial importance. The first goal of this work of thesis is to investigate the internal origin of warm and cold phases. Numerical simulations are the powerful tool of analysis. The way in which a spatially distributed cooling process originates has been examined and the off-centre amount of gas mass which cools when different and differently characterized AGN feedback mechanisms operate has been quantified. This thesis demonstrates that the aforementioned non-linear density perturbations originate and develop from AGN feedback mechanisms in a natural fashion. An internal origin of the warm phase from the once hot gas is shown to be possible. Computed velocity dispersions of ionized and hot gas are similar. The cold gas as well can originate from the cooling process: indeed, it has been estimated that the surrounding stellar radiation, which is one of the most feasible sources of ionization of the warm gas, does not manage to keep ionized all the gas at 10^4 K. Therefore, cooled gas does undergo a further cooling which can lead the warm phase to lower temperatures. However, the gas which has cooled from the hot phase is expected to be dustless; nonetheless, a large fraction of early type galaxies has detectable dust in their cores, both concentrated in filamentary and disky structures and spread over larger regions. Therefore a regularly rotating disk of cold and dusty gas has been included in the simulations. A new quantitative investigation of the spatially distributed cooling process has therefore been essential: the contribution of the included amount of dust which is embedded in the cold gas does have a role in promoting and enhancing the cooling. The fate of dust which was at first embedded in cold gas has been investigated. The role of AGN feedback mechanisms in dragging (if able) cold and dusty gas from the core of massive ellipticals up to large radii has been studied.