Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

VVER-440 Thermal Hydraulics as a Computer Code Validation Challenge

This thesis concentrates on the validation of a generic thermal hydraulic computer code TRACE under the challenges of the VVER-440 reactor type. The code capability to model the VVER-440 geometry and thermal hydraulic phenomena specific to this reactor design has been examined and demonstrated acceptable. The main challenge in VVER-440 thermal hydraulics appeared in the modelling of the horizontal steam generator. The major challenge here is not in the code physics or numerics but in the formulation of a representative nodalization structure. Another VVER-440 specialty, the hot leg loop seals, challenges the system codes functionally in general, but proved readily representable. Computer code models have to be validated against experiments to achieve confidence in code models. When new computer code is to be used for nuclear power plant safety analysis, it must first be validated against a large variety of different experiments. The validation process has to cover both the code itself and the code input. Uncertainties of different nature are identified in the different phases of the validation procedure and can even be quantified. This thesis presents a novel approach to the input model validation and uncertainty evaluation in the different stages of the computer code validation procedure. This thesis also demonstrates that in the safety analysis, there are inevitably significant uncertainties that are not statistically quantifiable; they need to be and can be addressed by other, less simplistic means, ultimately relying on the competence of the analysts and the capability of the community to support the experimental verification of analytical assumptions. This method completes essentially the commonly used uncertainty assessment methods, which are usually conducted using only statistical methods.

Typenpoiston tehostaminen kunnallisessa jätevedenpuhdistuksessa membraanibioreaktorin avulla

Tämä diplomityö tehtiin Vihdin Vesihuoltolaitoksen Nummelan jäteveden puhdistamolle. Työssä tutkittiin typenpoistoa kunnallisista jätevesistä membraanibioreaktorin (MBR) avulla. MBR:ssä yhdistyvät perinteinen aktiivilieteprosessi ja kalvosuodatus. Työn tavoite oli päästä yli 95 % typenpoistoon. Aluksi typenpoisto oli yli 80 %, kun pilot-mittakaavan MBR-laitosta operoitiin perinteisen prosessin parametrein. Typenpoistoa onnistuttiin tehostamaan nostamalla nitraattipitoisen palautuslietteen kierrätystä prosessin alkupäähän (1600 L/h) ja lisäämällä aktiivista biomassaa reaktorissa. Yli 90 % typenpoisto edellytti myös pidempää viipymäaikaa (noin kaksinkertainen perinteiseen prosessiin verrattuna). Tutkimuksessa päästiin parhaimmillaan jopa 95 % typenpoistumaan operoimalla laitteistoa pienellä typpikuormalla (0,1 kg/vrk) ja alhaisemmalla lietepitoisuudella (10 g/L). Typpikuorman noustessa (0,3 kg/vrk) typenpoistoteho laski. Tätä onnistuttiin parantamaan (yli 90 %) nostamalla biomassan määrää reaktorissa (15 g/L). Hyvän typenpoiston saavuttaminen edellytti myös suurempaa metanolin ja hapen syöttöä.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.

Spin Polarization and Normal Reflection at Nanoscale Point Contacts

This thesis presents point-contact measurements between superconductors (Nb, Ta, Sn,Al, Zn) and ferromagnets (Co, Fe, Ni) as well as non-magnetic metals (Ag, Au, Cu, Pt).The point contacts were fabricated using the shear method. The differential resistanceof the contacts was measured either in liquid He at 4.2 K or in vacuum in a dilutionrefrigerator at varying temperature down to 0.1 K. The contact properties were investigatedas function of size and temperature. The measured Andreev-reflection spectrawere analysed in the framework of the BTK model – a three parameter model that describescurrent transport across a superconductor - normal conductor interface. Theoriginal BTK model was modified to include the effects of spin polarization or finitelifetime of the Cooper pairs. Our polarization values for the ferromagnets at 4.2 K agree with the literature data, but the analysis was ambiguous because the experimental spectra both with ferromagnets and non-magnets could be described equally well either with spin polarization or finite lifetime effects in the BTK model. With the polarization model the Z parametervaries from almost 0 to 0.8 while the lifetime model produces Z values close to 0.5. Measurements at lower temperatures partly lift this ambiguity because the magnitude of thermal broadening is small enough to separate lifetime broadening from the polarization. The reduced magnitude of the superconducting anomalies for Zn-Fe contacts required an additional modification of the BTK model which was implemented as a scaling factor. Adding this parameter led to reduced polarization values. However, reliable data is difficult to obtain because different parameter sets produce almost identical spectra.

Possibility of press forming wood films into primary packaging

In recent times the packaging industry is finding means to maximize profit. Wood used to be the most advantageous and everyday material for packaging, worktables, counters, constructions, interiors, tools and as materials and utensils in the food companies in the world. The use of wood has declined vigorously, and other materials like plastic, ceramic, stainless steel, concrete, and aluminum have taken its place. One way that the industry could reduce its cost is by finding possibilities of using wood for primary packaging after which it can be safely recycled or burned as a carbon source for energy. Therefore, the main objective of this thesis is to investigate the possibility of press-forming a wood film into primary packaging. In order to achieve the stated objectives, discussion on major characteristics of wood in terms of structure, types and application were studied. Two different wood species, pine and birch were used for the experimental work. These were provided by a local carpentry workshop in Lappeenranta and a workshop in Ruokolahti supervised by Professor Timo Kärki. Laboratory tests were carried out at Lappeenranta University of Technology FMS workshop on Stenhøj EPS40 M hydraulic C-frame press coupled with National Instruments VI Logger and on the Adjustable packaging line machine at LUT Packaging laboratory. The tests succeeded better on the LUT packaging line than on the Stenhoj equipment due to the integrated heating system in the machine. However, there is much work to be done before the quality of a tray produced from the wood film is comparable to that of the wood plastic composite tray.

Novel virtual environment and real-time simulation based methods for improving life-cycle efficiency of non-road mobile machinery

Virtual environments and real-time simulators (VERS) are becoming more and more important tools in research and development (R&D) process of non-road mobile machinery (NRMM). The virtual prototyping techniques enable faster and more cost-efficient development of machines compared to use of real life prototypes. High energy efficiency has become an important topic in the world of NRMM because of environmental and economic demands. The objective of this thesis is to develop VERS based methods for research and development of NRMM. A process using VERS for assessing effects of human operators on the life-cycle efficiency of NRMM was developed. Human in the loop simulations are ran using an underground mining loader to study the developed process. The simulations were ran in the virtual environment of the Laboratory of Intelligent Machines of Lappeenranta University of Technology. A physically adequate real-time simulation model of NRMM was shown to be reliable and cost effective in testing of hardware components by the means of hardware-in-the-loop (HIL) simulations. A control interface connecting integrated electro-hydraulic energy converter (IEHEC) with virtual simulation model of log crane was developed. IEHEC consists of a hydraulic pump-motor and an integrated electrical permanent magnet synchronous motorgenerator. The results show that state of the art real-time NRMM simulators are capable to solve factors related to energy consumption and productivity of the NRMM. A significant variation between the test drivers is found. The results show that VERS can be used for assessing human effects on the life-cycle efficiency of NRMM. HIL simulation responses compared to that achieved with conventional simulation method demonstrate the advances and drawbacks of various possible interfaces between the simulator and hardware part of the system under study. Novel ideas for arranging the interface are successfully tested and compared with the more traditional one. The proposed process for assessing the effects of operators on the life-cycle efficiency will be applied for wider group of operators in the future. Driving styles of the operators can be analysed statistically from sufficient large result data. The statistical analysis can find the most life-cycle efficient driving style for the specific environment and machinery. The proposed control interface for HIL simulation need to be further studied. The robustness and the adaptation of the interface in different situations must be verified. The future work will also include studying the suitability of the IEHEC for different working machines using the proposed HIL simulation method.

Electrical parameters and water permeability properties of monolayers formed by T84 cells cultured on permeable supports

T84 is an established cell line expressing an enterocyte phenotype whose permeability properties have been widely explored. Osmotic permeability (P OSM), hydraulic permeability (P HYDR) and transport-associated net water fluxes (J W-transp), as well as short-circuit current (I SC), transepithelial resistance (R T), and potential difference (deltaV T) were measured in T84 monolayers with the following results: P OSM 1.3 ± 0.1 cm.s-1 x 10-3; P HYDR 0.27 ± 0.02 cm.s-1; R T 2426 ± 109 omega.cm², and deltaV T 1.31 ± 0.38 mV. The effect of 50 µM 5,6-dichloro-1-ethyl-1,3-dihydro-2H-benzimidazol-2-one (DCEBIO), a "net Cl- secretory agent", on T84 cells was also studied. We confirm the reported important increase in I SC induced by DCEBIO which was associated here with a modest secretory deltaJ W-transp. The present results were compared with those reported using the same experimental approach applied to established cell lines originating from intestinal and renal epithelial cells (Caco-2, LLC-PK1 and RCCD-1). No clear association between P HYDR and R T could be demonstrated and high P HYDR values were observed in an electrically tight epithelium, supporting the view that a "water leaky" barrier is not necessarily an "electrically leaky" one. Furthermore, the modest secretory deltaJ W-transp was not consistent with previous results obtained with RCCD-1 cells stimulated with vasopressin (absorptive fluxes) or with T84 cells secreting water under the action of Escherichia coli heat stable enterotoxin. We conclude that, while the presence of aquaporins is necessary to dissipate an external osmotic gradient, coupling between water and ion transport cannot be explained by a simple and common underlying mechanism.

Metal additive manufacturing of internal flow channels with powder bed fusion processes

This thesis studies the advantages, disadvantages and possibilities of additive manufacturing in making components with internal flow channels. These include hydraulic components, components with cooling channels and heat exchangers. Processes studied in this work are selective laser sintering and selective laser melting of metallic materials. The basic principles of processes and parameters involved in the process are presented and different possibilities of internal channel manufacturing and flow improvement are introduced

Remote sensing for three-dimensional modelling of hydromorphology

Successful management of rivers requires an understanding of the fluvial processes that govern them. This, in turn cannot be achieved without a means of quantifying their geomorphology and hydrology and the spatio-temporal interactions between them, that is, their hydromorphology. For a long time, it has been laborious and time-consuming to measure river topography, especially in the submerged part of the channel. The measurement of the flow field has been challenging as well, and hence, such measurements have long been sparse in natural environments. Technological advancements in the field of remote sensing in the recent years have opened up new possibilities for capturing synoptic information on river environments. This thesis presents new developments in fluvial remote sensing of both topography and water flow. A set of close-range remote sensing methods is employed to eventually construct a high-resolution unified empirical hydromorphological model, that is, river channel and floodplain topography and three-dimensional areal flow field. Empirical as well as hydraulic theory-based optical remote sensing methods are tested and evaluated using normal colour aerial photographs and sonar calibration and reference measurements on a rocky-bed sub-Arctic river. The empirical optical bathymetry model is developed further by the introduction of a deep-water radiance parameter estimation algorithm that extends the field of application of the model to shallow streams. The effect of this parameter on the model is also assessed in a study of a sandy-bed sub-Arctic river using close-range high-resolution aerial photography, presenting one of the first examples of fluvial bathymetry modelling from unmanned aerial vehicles (UAV). Further close-range remote sensing methods are added to complete the topography integrating the river bed with the floodplain to create a seamless high-resolution topography. Boat- cart- and backpack-based mobile laser scanning (MLS) are used to measure the topography of the dry part of the channel at a high resolution and accuracy. Multitemporal MLS is evaluated along with UAV-based photogrammetry against terrestrial laser scanning reference data and merged with UAV-based bathymetry to create a two-year series of seamless digital terrain models. These allow the evaluation of the methodology for conducting high-resolution change analysis of the entire channel. The remote sensing based model of hydromorphology is completed by a new methodology for mapping the flow field in 3D. An acoustic Doppler current profiler (ADCP) is deployed on a remote-controlled boat with a survey-grade global navigation satellite system (GNSS) receiver, allowing the positioning of the areally sampled 3D flow vectors in 3D space as a point cloud and its interpolation into a 3D matrix allows a quantitative volumetric flow analysis. Multitemporal areal 3D flow field data show the evolution of the flow field during a snow-melt flood event. The combination of the underwater and dry topography with the flow field yields a compete model of river hydromorphology at the reach scale.

Novel control, haptic and calibration methods for teleoperated electrohydraulic servo systems

The aim of this thesis is to propose a novel control method for teleoperated electrohydraulic servo systems that implements a reliable haptic sense between the human and manipulator interaction, and an ideal position control between the manipulator and the task environment interaction. The proposed method has the characteristics of a universal technique independent of the actual control algorithm and it can be applied with other suitable control methods as a real-time control strategy. The motivation to develop this control method is the necessity for a reliable real-time controller for teleoperated electrohydraulic servo systems that provides highly accurate position control based on joystick inputs with haptic capabilities. The contribution of the research is that the proposed control method combines a directed random search method and a real-time simulation to develop an intelligent controller in which each generation of parameters is tested on-line by the real-time simulator before being applied to the real process. The controller was evaluated on a hydraulic position servo system. The simulator of the hydraulic system was built based on Markov chain Monte Carlo (MCMC) method. A Particle Swarm Optimization algorithm combined with the foraging behavior of E. coli bacteria was utilized as the directed random search engine. The control strategy allows the operator to be plugged into the work environment dynamically and kinetically. This helps to ensure the system has haptic sense with high stability, without abstracting away the dynamics of the hydraulic system. The new control algorithm provides asymptotically exact tracking of both, the position and the contact force. In addition, this research proposes a novel method for re-calibration of multi-axis force/torque sensors. The method makes several improvements to traditional methods. It can be used without dismantling the sensor from its application and it requires smaller number of standard loads for calibration. It is also more cost efficient and faster in comparison to traditional calibration methods. The proposed method was developed in response to re-calibration issues with the force sensors utilized in teleoperated systems. The new approach aimed to avoid dismantling of the sensors from their applications for applying calibration. A major complication with many manipulators is the difficulty accessing them when they operate inside a non-accessible environment; especially if those environments are harsh; such as in radioactive areas. The proposed technique is based on design of experiment methodology. It has been successfully applied to different force/torque sensors and this research presents experimental validation of use of the calibration method with one of the force sensors which method has been applied to.

Simulation model of wind turbine gearbox and lubrication system

As increasing efficiency of a wind turbine gearbox, more power can be transferred from rotor blades to generator and less power is used to cause wear and heating in the gearbox. By using a simulation model, behavior of the gearbox can be studied before creating expensive prototypes. The objective of the thesis is to model a wind turbine gearbox and its lubrication system to study power losses and heat transfer inside the gearbox and to study the simulation methods of the used software. Software used to create the simulation model is Siemens LMS Imagine.Lab AMESim, which can be used to create one-dimensional mechatronic system simulation models from different fields of engineering. When combining components from different libraries it is possible to create a simulation model, which includes mechanical, thermal and hydraulic models of the gearbox. Results for mechanical, thermal, and hydraulic simulations are presented in the thesis. Due to the large scale of the wind turbine gearbox and the amount of power transmitted, power loss calculations from AMESim software are inaccurate and power losses are modelled as constant efficiency for each gear mesh. Starting values for simulation in thermal and hydraulic simulations were chosen from test measurements and from empirical study as compact and complex design of gearbox prevents accurate test measurements. In further studies to increase the accuracy of the simulation model, components used for power loss calculations needs to be modified and values for unknown variables are needed to be solved through accurate test measurements.

Risk Assessment of Modern Landfill Structures in Finland

The purpose of this thesis was to investigate environmental permits of landfills with respect to the appropriateness of risk assessments focusing on contaminant migration, structures capable to protect the environment, waste and leachate management and existing environmental impacts of landfills. According to the requirements, a risk assessment is always required to demonstrate compliance with environmental protection requirements if the environmental permit decision deviates from the set requirements. However, there is a reason to doubt that all relevant risk factors are identified in current risk assessment practices in order to protect people end environment. In this dissertation, risk factors were recognized in 12 randomly selected landfills. Based on this analysis, a structural risk assessment method was created. The method was verified with two case examples. Several development needs were found in the risk assessments of the environmental permit decisions. The risk analysis equations used in the decisions did not adequately take into account all the determining factors like waste prospects, total risk quantification or human delineated factors. Instead of focusing on crucial factors, the landfill environmental protection capability is simply expressed via technical factors like hydraulic conductivity. In this thesis, it could be shown, that using adequate risk assessment approaches the most essential environmental impacts can be taken into account by consideration of contaminant transport mechanisms, leachate effects, and artificial landfill structures. The developed structural risk analysing (SRA) method shows, that landfills structures could be designed in a more cost-efficient way taking advantage of recycled or by-products. Additionally, the research results demonstrate that the environmental protection requirements of landfills should be updated to correspond to the capability to protect the environment instead of the current simplified requirements related to advective transport only.

Analysis of PACTEL small break LOCA experiment using TRACE

A small break loss-of-coolant accident (SBLOCA) is one of problems investigated in an NPP operation. Such accident can be analyzed using an experiment facility and TRACE thermal-hydraulic system code. A series of SBLOCA experiments was carried out on Parallel Channel Test Loop (PACTEL) facility, exploited together with Technical Research Centre of Finland VTT Energy and Lappeenranta University of Technology (LUT), in order to investigate two-phase phenomena related to a VVER-type reactor. The experiments and a TRACE model of the PACTEL facility are described in the paper. In addition, there is the TRACE code description with main field equations. At the work, calculations of a SBLOCA series are implemented and after the calculations, the thesis discusses the validation of TRACE and concludes with an assessment of the usefulness and accuracy of the code in calculating small breaks.

Coolability of porous core debris beds : Effects of bed geometry and multi-dimensional flooding

This thesis addresses the coolability of porous debris beds in the context of severe accident management of nuclear power reactors. In a hypothetical severe accident at a Nordic-type boiling water reactor, the lower drywell of the containment is flooded, for the purpose of cooling the core melt discharged from the reactor pressure vessel in a water pool. The melt is fragmented and solidified in the pool, ultimately forming a porous debris bed that generates decay heat. The properties of the bed determine the limiting value for the heat flux that can be removed from the debris to the surrounding water without the risk of re-melting. The coolability of porous debris beds has been investigated experimentally by measuring the dryout power in electrically heated test beds that have different geometries. The geometries represent the debris bed shapes that may form in an accident scenario. The focus is especially on heap-like, realistic geometries which facilitate the multi-dimensional infiltration (flooding) of coolant into the bed. Spherical and irregular particles have been used to simulate the debris. The experiments have been modeled using 2D and 3D simulation codes applicable to fluid flow and heat transfer in porous media. Based on the experimental and simulation results, an interpretation of the dryout behavior in complex debris bed geometries is presented, and the validity of the codes and models for dryout predictions is evaluated. According to the experimental and simulation results, the coolability of the debris bed depends on both the flooding mode and the height of the bed. In the experiments, it was found that multi-dimensional flooding increases the dryout heat flux and coolability in a heap-shaped debris bed by 47–58% compared to the dryout heat flux of a classical, top-flooded bed of the same height. However, heap-like beds are higher than flat, top-flooded beds, which results in the formation of larger steam flux at the top of the bed. This counteracts the effect of the multi-dimensional flooding. Based on the measured dryout heat fluxes, the maximum height of a heap-like bed can only be about 1.5 times the height of a top-flooded, cylindrical bed in order to preserve the direct benefit from the multi-dimensional flooding. In addition, studies were conducted to evaluate the hydrodynamically representative effective particle diameter, which is applied in simulation models to describe debris beds that consist of irregular particles with considerable size variation. The results suggest that the effective diameter is small, closest to the mean diameter based on the number or length of particles.