Validation Studies of Thermal-hydraulic Code for Safety Analysis of Nuclear Power Plants

This thesis gives an overview of the validation process for thermal hydraulic system codes and it presents in more detail the assessment and validation of the French code CATHARE for VVER calculations. Three assessment cases are presented: loop seal clearing, core reflooding and flow in a horizontal steam generator. The experience gained during these assessment and validation calculations has been used to analyze the behavior of the horizontal steam generator and the natural circulation in the geometry of the Loviisa nuclear power plant. The cases presented are not exhaustive, but they give a good overview of the work performed by the personnel of Lappeenranta University of Technology (LUT). Large part of the work has been performed in co-operation with the CATHARE-team in Grenoble, France. The design of a Russian type pressurized water reactor, VVER, differs from that of a Western-type PWR. Most of thermal-hydraulic system codes are validated only for the Western-type PWRs. Thus, the codes should be assessed and validated also for VVER design in order to establish any weaknesses in the models. This information is needed before codes can be used for the safety analysis. Theresults of the assessment and validation calculations presented here show that the CATHARE code can be used also for the thermal-hydraulic safety studies for VVER type plants. However, some areas have been indicated which need to be reassessed after further experimental data become available. These areas are mostly connected to the horizontal stem generators, like condensation and phase separation in primary side tubes. The work presented in this thesis covers a large numberof the phenomena included in the CSNI code validation matrices for small and intermediate leaks and for transients. Also some of the phenomena included in the matrix for large break LOCAs are covered. The matrices for code validation for VVER applications should be used when future experimental programs are planned for code validation.

Experimental insights into deformation dynamics and intermittency in rapid granular shear flows

This work concerns the experimental study of rapid granular shear flows in annular Couette geometry. The flow is induced by continuous driving of the horizontal plate at the top of the granular bed in an annulus. The compressive pressure, driving torque, instantaneous bed height and rotational speed of the shearing plate are measured. Moreover, local stress fluctuations are measured in a medium made of steel spheres 2 and 3 mm in diameter. Both monodisperse packing and bidisperse packing are investigated to reveal the influence of size diversity in intermittent features of granular materials. Experiments are conducted in an annulus that can contain up to 15 kg of spherical steel balls. The shearing granular medium takes place via the rotation of the upper plate which compresses the material loaded inside the annulus. Fluctuations of compressive force are locally measured at the bottom of the annulus using a piezoelectric sensor. Rapid shear flow experiments are pursued at different compressive forces and shear rates and the sensitivity of fluctuations are then investigated by different means through monodisperse and bidisperse packings. Another important feature of rapid granular shear flows is the formation of ordered structures upon shearing. It requires a certain range for the amount of granular material (uniform size distribution) loaded in the system in order to obtain stable flows. This is studied more deeply in this thesis. The results of the current work bring some new insights into deformation dynamics and intermittency in rapid granular shear flows. The experimental apparatus is modified in comparison to earlier investigations. The measurements produce data for various quantities continuously sampled from the start of shearing to the end. Static failure and dynamic shearing ofa granular medium is investigated. The results of this work revealed some important features of failure dynamics and structure formation in the system. Furthermore, some computer simulations are performed in a 2D annulus to examine the nature of kinetic energy dissipation. It is found that turbulent flow models can statistically represent rapid granular flows with high accuracy. In addition to academic outcomes and scientific publications our results have a number of technological applications associated with grinding, mining and massive grain storages.

Thermal energy storage with phase change materials in building envelopes

Els materials de canvi de fase (PCM) han estat considerats per a l’emmagatzematge tèrmic en edificis des de 1980. Amb la inclusió dels PCM en plaques de guix, guix, formigó o altres materials que s’utilitzen per a cobrir les parets, l’emmagatzematge tèrmic pot ser part de les estructures fins i tot en edificis lleugers. Les noves tècniques de microencapsulació han obert moltes possibilitats en aplicacions per a edificis. El treball que es presenta és el desenvolupament d’un formigó innovador mesclat amb PCM microencapsulat, amb un punt de fusió de 26 oC i una entalpia de canvi de fase de 110 kJ/kg. El primer experiment va ser la inclusió del PCM microencapsulat dins del formigó i la construcció d’una caseta amb aquest nou formigó-PCM. Es va construir una segona caseta al costat de la primera amb les mateixes característiques i orientació però amb formigó convencional que serveix com a referència. Durant els anys 2005 i 2006 es va analitzar el comportament d’ambdues casetes i més tard es va edificar un mur Trombe a la paret sud de totes dues per investigar la seva influència durant la tardor i l’hivern.

Determination of thermal diffusivity in papaya pulp as a function of maturation stage

In order to determine the penetration of the thermal wave in the papaya fruit pulp (Carica papaya L.), cv. Golden, thermal diffusivity of the pulp was obtained measuring temperature at four different depths. Measurements were carried out initially with the fruit on the first stage of maturity. The changes of the thermal diffusivity were expressed as a function of ripening. A temporal decrease of the thermal diffusivity was observed. Chemical (pH, soluble solids and total titratable acidity) and physical (pulp firmness) properties were measured as well and the results were compared to the thermal diffusivity change.

Comblike poly(a-alkyl g glutamate)s: computer simulation studies of an intermediate thermal phase.

Monte Carlo (MC) simulations have been used to study the structure of an intermediate thermal phase of poly(R-octadecyl ç,D-glutamate). This is a comblike poly(ç-peptide) able to adopt a biphasic structure that has been described as a layered arrangement of backbone helical rods immersed in a paraffinic pool of polymethylene side chains. Simulations were performed at two different temperatures (348 and 363 K), both of them above the melting point of the paraffinic phase, using the configurational bias MC algorithm. Results indicate that layers are constituted by a side-by-side packing of 17/5 helices. The organization of the interlayer paraffinic region is described in atomistic terms by examining the torsional angles and the end-to-end distances for the octadecyl side chains. Comparison with previously reported comblike poly(â-peptide)s revealed significant differences in the organization of the alkyl side chains.

Relaxation and derelaxation of pure and hydrogenated amorphous silicon during thermal annealing experiments

The structural relaxation of pure amorphous silicon a-Si and hydrogenated amorphous silicon a-Si:H materials, that occurs during thermal annealing experiments, has been analyzed by Raman spectroscopy and differential scanning calorimetry. Unlike a-Si, the heat evolved from a-Si:H cannot be explained by relaxation of the Si-Si network strain but it reveals a derelaxation of the bond angle strain. Since the state of relaxation after annealing is very similar for pure and hydrogenated materials, our results give strong experimental support to the predicted configurational gap between a-Si and crystalline silicon.

Dynamic and quasistatic trajectories in quasifission reactions and particle emission

We show that the quasifission paths predicted by the one-body dissipation dynamics, in the slowest phase of a binary reaction, follow a quasistatic path, which represents a sequence of states of thermal equilibrium at a fixed value of the deformation coordinate. This establishes the use of the statistical particle-evaporation model in the case of dynamical time-evolving systems. Pre- and post-scission multiplicities of neutrons and total multiplicities of protons and α particles in fission reactions of 63Cu+92Mo, 60Ni+100Mo, 63Cu+100Mo at 10 MeV/u and 20Ne+144,148,154Sm at 20 MeV/u are reproduced reasonably well with statistical model calculations performed along dynamic trajectories whose slow stage (from the most compact configuration up to the point where the neck starts to develop) lasts some 35×10−21 s.

Development of rotating drum under mechanical and thermal load

A variable temperature field sets exacting demands to the structure under mechanical load. Most of all the lifetime of the rotating drum structure depends on temperature differences between parts inside the drum. The temperature difference was known because of the measurements made before. The list of demands was created based on customers’ needs. The limits of this paper were set to the inner structure of the drum. Creation of ideas for the inner structure was started open minded. The main principle in the creation process was to create new ideas for the function of the product with the help of sub-functions. The sub-functions were created as independent as possible. The best sub-functions were combined together and the new working principles were created based on them. Every working principle was calculated separately and criticized at the end of the calculation process. The main objective was to create the new kind of structure, which is not based too much to the old, inoperative structure. The affect of own weight of the inner structure to the stress values was quite small but it was also taken into consideration when calculating the maximum stress value of the structure. Because of very complex structures all of the calculations were made with the help of the ProE – Mechanica software. The fatigue analyze was made also for the best structure solution.