Klas Jareteg

Behaviour and Stability of the Two-Fluid Model for Fine-Scale Simulations of Bubbly Flow in Nuclear Reactors

Abstract In the present work, we formulate a simplistic two-fluid model for bubbly steam-water flow existing between fuel pins in nuclear fuel assemblies. Numerical simulations are performed in periodic 2D domains of varying sizes. The appearance of a non-uniform volume fraction field in the form of meso-scales is investigated and shown to be varying with the bubble loading and the domain size, as well as with the numerical algorithm employed. These findings highlight the difficulties involved in interpreting the occurrence of instabilities in two-fluid simulations of gas-liquid flows, where physical and unphysical instabilities are prone to be confounded. The results obtained in this work therefore contribute to a rigorous foundation in on-going efforts to derive a consistent meso-scale formulation of the traditional two-fluid model for multiphase flows in nuclear reactors.

Development of a point-kinetic verification scheme for nuclear reactor applications

In this paper, a new method that can be used for checking the proper implementation of time- or frequency-dependent neutron transport models and for verifying their ability to recover some basic reactor physics properties is proposed. This method makes use of the application of a stationary perturbation to the system at a given frequency and extraction of the point-kinetic component of the system response. Even for strongly heterogeneous systems for which an analytical solution does not exist, the point-kinetic component follows, as a function of frequency, a simple analytical form. The comparison between the extracted point-kinetic component and its expected analytical form provides an opportunity to verify and validate neutron transport solvers. The proposed method is tested on two diffusion-based codes, one working in the time domain and the other working in the frequency domain. As long as the applied perturbation has a non-zero reactivity effect, it is demonstrated that the method can be successfully applied to verify and validate time- or frequency-dependent neutron transport solvers. Although the method is demonstrated in the present paper in a diffusion theory framework, higher order neutron transport methods could be verified based on the same principles.

A numerical framework for bubble transport in a subcooled fluid flow

In this paper we present a framework for the simulation of dispersed bubbly two-phase flows, with the specific aim of describing vapor–liquid systems with condensation. We formulate and implement a framework that consists of a population balance equation (PBE) for the bubble size distribution and an Eulerian–Eulerian two-fluid solver. The PBE is discretized using the Direct Quadrature Method of Moments (DQMOM) in which we include the condensation of the bubbles as an internal phase space convection. We investigate the robustness of the DQMOM formulation and the numerical issues arising from the rapid shrinkage of the vapor bubbles. In contrast to a PBE method based on the multiple-size-group (MUSIG) method, the DQMOM formulation allows us to compute a distribution with dynamic bubble sizes. Such a property is advantageous to capture the wide range of bubble sizes associated with the condensation process. Furthermore, we compare the computational performance of the DQMOM-based framework with the MUSIG method. The results demonstrate that DQMOM is able to retrieve the bubble size distribution with a good numerical precision in only a small fraction of the computational time required by MUSIG. For the two-fluid solver, we examine the implementation of the mass, momentum and enthalpy conservation equations in relation to the coupling to the PBE. In particular, we propose a formulation of the pressure and liquid continuity equations, that was shown to correctly preserve mass when computing the vapor fraction with DQMOM. In addition, the conservation of enthalpy was also proven. Therefore a consistent overall framework that couples the PBE and two-fluid solvers is achieved.

Phoswich scintillator for proton and gamma radiation of high energy

We present here a Phoswich scintillator design to achieve both high resolution gamma ray detection, and good efficiency for high energy protons. There are recent developments of new high resolution scintillator materials. Especially the LaBr3(Ce) and LaCl3(Ce) crystals have very good energy resolution in the order of 3% for 662 keV gamma radiation. In addition, these materials exhibit a very good light output (63 and 32 photons/keV respectively).A demonstrator detector in the form of an Al cylinder of 24 mm diameter and a total length of 80 mm with 2 mm wall thickness, containing a LaBr3(Ce) crystal of 20 mm diameter and 30 mm length directly coupled to a LaCl3(Ce) crystal of 50 mm length, and closed with a glass window of 5 mm, was delivered by Saint Gobain. To the glass window a Hamamatsu R5380 Photomultiplier tube (PMT) was coupled using silicon optical grease.