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.

Surface studies of limestones and dolostones : characterisation using various techniques and batch dissolution experiments with hydrochloric acid solutions

In this thesis, stepwise titration with hydrochloric acid was used to obtain chemical reactivities and dissolution rates of ground limestones and dolostones of varying geological backgrounds (sedimentary, metamorphic or magmatic). Two different ways of conducting the calculations were used: 1) a first order mathematical model was used to calculate extrapolated initial reactivities (and dissolution rates) at pH 4, and 2) a second order mathematical model was used to acquire integrated mean specific chemical reaction constants (and dissolution rates) at pH 5. The calculations of the reactivities and dissolution rates were based on rate of change of pH and particle size distributions of the sample powders obtained by laser diffraction. The initial dissolution rates at pH 4 were repeatedly higher than previously reported literature values, whereas the dissolution rates at pH 5 were consistent with former observations. Reactivities and dissolution rates varied substantially for dolostones, whereas for limestones and calcareous rocks, the variation can be primarily explained by relatively large sample standard deviations. A list of the dolostone samples in a decreasing order of initial reactivity at pH 4 is: 1) metamorphic dolostones with calcite/dolomite ratio higher than about 6% 2) sedimentary dolostones without calcite 3) metamorphic dolostones with calcite/dolomite ratio lower than about 6% The reactivities and dissolution rates were accompanied by a wide range of experimental techniques to characterise the samples, to reveal how different rocks changed during the dissolution process, and to find out which factors had an influence on their chemical reactivities. An emphasis was put on chemical and morphological changes taking place at the surfaces of the particles via X-ray Photoelectron Spectroscopy (XPS) and Scanning Electron Microscopy (SEM). Supporting chemical information was obtained with X-Ray Fluorescence (XRF) measurements of the samples, and Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) and Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) measurements of the solutions used in the reactivity experiments. Information on mineral (modal) compositions and their occurrence was provided by X-Ray Diffraction (XRD), Energy Dispersive X-ray analysis (EDX) and studying thin sections with a petrographic microscope. BET (Brunauer, Emmet, Teller) surface areas were determined from nitrogen physisorption data. Factors increasing chemical reactivity of dolostones and calcareous rocks were found to be sedimentary origin, higher calcite concentration and smaller quartz concentration. Also, it is assumed that finer grain size and larger BET surface areas increase the reactivity although no certain correlation was found in this thesis. Atomic concentrations did not correlate with the reactivities. Sedimentary dolostones, unlike metamorphic ones, were found to have porous surface structures after dissolution. In addition, conventional (XPS) and synchrotron based (HRXPS) X-ray Photoelectron Spectroscopy were used to study bonding environments on calcite and dolomite surfaces. Both samples are insulators, which is why neutralisation measures such as electron flood gun and a conductive mask were used. Surface core level shifts of 0.7 ± 0.1 eV for Ca 2p spectrum of calcite and 0.75 ± 0.05 eV for Mg 2p and Ca 3s spectra of dolomite were obtained. Some satellite features of Ca 2p, C 1s and O 1s spectra have been suggested to be bulk plasmons. The origin of carbide bonds was suggested to be beam assisted interaction with hydrocarbons found on the surface. The results presented in this thesis are of particular importance for choosing raw materials for wet Flue Gas Desulphurisation (FGD) and construction industry. Wet FGD benefits from high reactivity, whereas construction industry can take advantage of slow reactivity of carbonate rocks often used in the facades of fine buildings. Information on chemical bonding environments may help to create more accurate models for water-rock interactions of carbonates.

Modelling the super-equilibria in thermal biomass conversion : applications and limitations of the constrained free energy method

Bioenergi ses som en viktig del av det nu- och framtida sortimentet av inhemsk energi. Svartlut, bark och skogsavfall täcker mer än en femtedel av den inhemska energianvändningen. Produktionsanläggningar kan fungera ofullständigt och en mängd gas-, partikelutsläpp och tjära produceras samtidigt och kan leda till beläggningsbildning och korrosion. Orsaken till dessa problem är ofta obalans i processen: vissa föreningar anrikas i processen och superjämviktstillstånd är bildas. I denna doktorsavhandling presenteras en ny beräkningsmetod, med vilken man kan beskriva superjämviktstillståndet, de viktigaste kemiska reaktionerna, processens värmeproduktion och tillståndsstorheter samtidigt. Beräkningsmetoden grundar sig på en unik frienergimetod med bivillkor som har utvecklats vid VTT. Den här så kallade CFE-metoden har tidigare utnyttjats i pappers-, metall- och kemiindustrin. Applikationer för bioenergi, vilka är demonstrerade i doktorsavhandlingen, är ett nytt användingsområde för metoden. Studien visade att beräkningsmetoden är väl lämpad för högtemperaturenergiprocesser. Superjämviktstillstånden kan uppstå i dessa processer och det kemiska systemet kan definieras med några bivillkor. Typiska tillämpningar är förbränning av biomassa och svartlut, förgasning av biomassa och uppkomsten av kväveoxider. Också olika sätt att definiera superjämviktstillstånd presenterades i doktorsavhandlingen: empiriska konstanter, empiriska hastighetsuttryck eller reaktionsmekanismer kan användas. Resultaten av doktorsavhandlingen kan utnyttjas i framtiden i processplaneringen och i undersökning av nya tekniska lösningar för förgasning, förbränningsteknik och biobränslen. Den presenterade metoden är ett bra alternativ till de traditionella mekanistiska och fenomenmodeller och kombinerar de bästa delarna av både. --------------------------------------------------------------- Bioenergia on tärkeä osa nykyistä ja tulevaa kotimaista energiapalettia. Mustalipeä, kuori ja metsätähteet kattavat yli viidenneksen kotimaisesta energian kulutuksesta. Tuotantolaitokset eivät kuitenkaan aina toimi täydellisesti ja niiden prosesseissa syntyy erilaisia kaasu- ja hiukkaspäästöjä, tervoja sekä prosessilaitteita kuluttavia saostumia ja ruostumista. Usein syy näihin ongelmiin on prosessissa esiintyvä epätasapainotila: tietyt yhdisteet rikastuvat prosessissa ja muodostavat supertasapainotiloja. Väitöstyössä kehitettiin uusi laskentamenetelmä, jolla voidaan kuvata nämä supertasapainotilat, tärkeimmät niihin liittyvät kemialliset reaktiot, prosessin lämmöntuotanto ja tilansuureet yhtä aikaa. Laskentamenetelmä perustuu VTT:llä kehitettyyn ainutlaatuiseen rajoitettuun vapaaenergiamenetelmään. Tätä niin kutsuttua CFE-menetelmää on aiemmin sovelluttu onnistuneesti muun muassa paperi-, metalli- ja kemianteollisuudessa. Väitöstyössä esitetyt bioenergiasovellukset ovat uusi sovellusalue menetelmälle. Työ osoitti laskentatavan soveltuvan hyvin korkealämpöisiin energiatekniikan prosesseihin, joissa kemiallista systeemiä rajoittavia tekijöitä oli rajallinen määrä ja siten super-tasapainotila saattoi muodostua prosessin aikana. Tyypillisiä sovelluskohteita ovat biomassan ja mustalipeän poltto, biomassan kaasutus ja typpioksidipäästöt. Työn aikana arvioitiin myös erilaisia tapoja määritellä super-tasapainojen muodostumista rajoittavat tekijät. Rajoitukset voitiin tehdä teollisiin mittauksiin pohjautuen, kokeellisia malleja hyödyntäen tai mekanistiseen reaktiokinetiikkaan perustuen. Tulevaisuudessa väitöstyön tuloksia voidaan hyödyntää prosessisuunnittelussa ja tutkittaessa uusia teknisiä ratkaisuja kaasutus- ja polttotekniikoissa sekä biopolttoaineiden tutkimuksessa. Kehitetty menetelmä tarjoaa hyvän vaihtoehdon perinteisille mekanistisille ja ilmiömalleille yhdistäen näiden parhaita puolia.

Electromagnetic and thermal design of a multilevel converter with high power density and reliability

Electric energy demand has been growing constantly as the global population increases. To avoid electric energy shortage, renewable energy sources and energy conservation are emphasized all over the world. The role of power electronics in energy saving and development of renewable energy systems is significant. Power electronics is applied in wind, solar, fuel cell, and micro turbine energy systems for the energy conversion and control. The use of power electronics introduces an energy saving potential in such applications as motors, lighting, home appliances, and consumer electronics. Despite the advantages of power converters, their penetration into the market requires that they have a set of characteristics such as high reliability and power density, cost effectiveness, and low weight, which are dictated by the emerging applications. In association with the increasing requirements, the design of the power converter is becoming more complicated, and thus, a multidisciplinary approach to the modelling of the converter is required. In this doctoral dissertation, methods and models are developed for the design of a multilevel power converter and the analysis of the related electromagnetic, thermal, and reliability issues. The focus is on the design of the main circuit. The electromagnetic model of the laminated busbar system and the IGBT modules is established with the aim of minimizing the stray inductance of the commutation loops that degrade the converter power capability. The circular busbar system is proposed to achieve equal current sharing among parallel-connected devices and implemented in the non-destructive test set-up. In addition to the electromagnetic model, a thermal model of the laminated busbar system is developed based on a lumped parameter thermal model. The temperature and temperature-dependent power losses of the busbars are estimated by the proposed algorithm. The Joule losses produced by non-sinusoidal currents flowing through the busbars in the converter are estimated taking into account the skin and proximity effects, which have a strong influence on the AC resistance of the busbars. The lifetime estimation algorithm was implemented to investigate the influence of the cooling solution on the reliability of the IGBT modules. As efficient cooling solutions have a low thermal inertia, they cause excessive temperature cycling of the IGBTs. Thus, a reliability analysis is required when selecting the cooling solutions for a particular application. The control of the cooling solution based on the use of a heat flux sensor is proposed to reduce the amplitude of the temperature cycles. The developed methods and models are verified experimentally by a laboratory prototype.

Permanent magnet synchronous generator design solution for large directdrive wind turbines: Thermal behavior of the LC DD-PMSG

Wind is one of the most compelling forms of indirect solar energy. Available now, the conversion of wind power into electricity is and will continue to be an important element of energy self-sufficiency planning. This paper is one in a series intended to report on the development of a new type of generator for wind energy; a compact, high-power, direct-drive permanent magnet synchronous generator (DD-PMSG) that uses direct liquid cooling (LC) of the stator windings to manage Joule heating losses. The main param-eters of the subject LC DD-PMSG are 8 MW, 3.3 kV, and 11 Hz. The stator winding is cooled directly by deionized water, which flows through the continuous hollow conductor of each stator tooth-coil winding. The design of the machine is to a large degree subordinate to the use of these solid-copper tooth-coils. Both steady-state and timedependent temperature distributions for LC DD-PMSG were examined with calculations based on a lumpedparameter thermal model, which makes it possible to account for uneven heat loss distribution in the stator conductors and the conductor cooling system. Transient calculations reveal the copper winding temperature distribution for an example duty cycle during variable-speed wind turbine operation. The cooling performance of the liquid cooled tooth-coil design was predicted via finite element analysis. An instrumented cooling loop featuring a pair of LC tooth-coils embedded in a lamination stack was built and laboratory tested to verify the analytical model. Predicted and measured results were in agreement, confirming the predicted satisfactory operation of the LC DD-PMSG cooling technology approach as a whole.

Characterisation of Cladonia rangiferina transcriptome and genome, and the effects of dehydration and rehydration on its gene expression

Lichens are symbiotic organisms, which consist of the fungal partner and the photosynthetic partner, which can be either an alga or a cyanobacterium. In some lichen species the symbiosis is tripartite, where the relationship includes both an alga and a cyanobacterium alongside the primary symbiont, fungus. The lichen symbiosis is an evolutionarily old adaptation to life on land and many extant fungal species have evolved from lichenised ancestors. Lichens inhabit a wide range of habitats and are capable of living in harsh environments and on nutrient poor substrates, such as bare rocks, often enduring frequent cycles of drying and wetting. Most lichen species are desiccation tolerant, and they can survive long periods of dehydration, but can rapidly resume photosynthesis upon rehydration. The molecular mechanisms behind lichen desiccation tolerance are still largely uncharacterised and little information is available for any lichen species at the genomic or transcriptomic level. The emergence of the high-throughput next generation sequencing (NGS) technologies and the subsequent decrease in the cost of sequencing new genomes and transcriptomes has enabled non-model organism research on the whole genome level. In this doctoral work the transcriptome and genome of the grey reindeer lichen, Cladonia rangiferina, were sequenced, de novo assembled and characterised using NGS and traditional expressed sequence tag (EST) technologies. RNA extraction methods were optimised to improve the yield and quality of RNA extracted from lichen tissue. The effects of rehydration and desiccation on C. rangiferina gene expression on whole transcriptome level were studied and the most differentially expressed genes were identified. The secondary metabolites present in C. rangiferina decreased the quality – integrity, optical characteristics and utility for sensitive molecular biological applications – of the extracted RNA requiring an optimised RNA extraction method for isolating sufficient quantities of high-quality RNA from lichen tissue in a time- and cost-efficient manner. The de novo assembly of the transcriptome of C. rangiferina was used to produce a set of contiguous unigene sequences that were used to investigate the biological functions and pathways active in a hydrated lichen thallus. The de novo assembly of the genome yielded an assembly containing mostly genes derived from the fungal partner. The assembly was of sufficient quality, in size similar to other lichen-forming fungal genomes and included most of the core eukaryotic genes. Differences in gene expression were detected in all studied stages of desiccation and rehydration, but the largest changes occurred during the early stages of rehydration. The most differentially expressed genes did not have any annotations, making them potentially lichen-specific genes, but several genes known to participate in environmental stress tolerance in other organisms were also identified as differentially expressed.

Thermal modelling of commercial lithium-ion batteries

The accelerating adoption of electrical technologies in vehicles over the recent years has led to an increase in the research on electrochemical energy storage systems, which are among the key elements in these technologies. The application of electrochemical energy storage systems for instance in hybrid electrical vehicles (HEVs) or hybrid mobile working machines allows tolerating high power peaks, leading to an opportunity to downsize the internal combustion engine and reduce fuel consumption, and therefore, CO2 and other emissions. Further, the application of electrochemical energy storage systems provides an option of kinetic and potential energy recuperation. Presently, the lithium-ion (Li-ion) battery is considered the most suitable electrochemical energy storage type in HEVs and hybrid mobile working machines. However, the intensive operating cycle produces high heat losses in the Li-ion battery, which increase its operating temperature. The Li-ion battery operation at high temperatures accelerates the ageing of the battery, and in the worst case, may lead to a thermal runaway and fire. Therefore, an appropriate Li-ion battery cooling system should be provided for the temperature control in applications such as HEVs and mobile working machines. In this doctoral dissertation, methods are presented to set up a thermal model of a single Li-ion cell and a more complex battery module, which can be used if full information about the battery chemistry is not available. In addition, a non-destructive method is developed for the cell thermal characterization, which allows to measure the thermal parameters at different states of charge and in different points of cell surface. The proposed models and the cell thermal characterization method have been verified by experimental measurements. The minimization of high thermal non-uniformity, which was detected in the pouch cell during its operation with a high C-rate current, was analysed by applying a simplified pouch cell 3D thermal model. In the analysis, heat pipes were incorporated into the pouch cell cooling system, and an optimization algorithm was generated for the estimation of the optimalplacement of heat pipes in the pouch cell cooling system. An analysis of the application of heat pipes to the pouch cell cooling system shows that heat pipes significantly decrease the temperature non-uniformity on the cell surface, and therefore, heat pipes were recommended for the enhancement of the pouch cell cooling system.

Development of laboratory exercises for basic nuclear thermal hydraulics measurements

The thesis focuses on light water reactors (pressurized water reactors, boiling water reactors) and measurement techniques for basic thermal hydraulics parameters that are used in a nuclear power plant. The goal of this work is a development of laboratory exercises for basic nuclear thermal hydraulics measurements.

Determination of the thermal parameters of high-power batteries by local heat

A new approach to the determination of the thermal parameters of high-power batteries is introduced here. Application of local heat flux measurement with a gradient heat flux sensor (GHFS) allows determination of the cell thermal parameters in di_erent surface points of the cell. The suggested methodology is not cell destructive as it does not require deep discharge of the cell or application of any charge/discharge cycles during measurements of the thermal parameters of the cell. The complete procedure is demonstrated on a high-power Li-ion pouch cell, and it is verified on a sample with well-known thermal parameters. A comparison of the experimental results with conventional thermal characterization methods shows an acceptably low error. The dependence of the cell thermal parameters on state of charge (SoC) and measurement points on the surface was studied by the proposed measurement approach.

Visualization of thermal-hydraulic nodes

The purpose of this thesis was to develop a program that can illustrate thermal-hydraulic node dimensions used in SMABRE simulations. These created node illustrations are used to verify the correctness of the designed simulation model and in addition they can be included in scientific reports. This thesis will include theory about SMABRE and relevant programs that were used to achieve the ending results. This thesis will give explanations for different modules that were created and used in the finished program, and it will present the different problems encountered and provide the solutions. The most important objective in this thesis is to display the results of generic VVER-1000 node dimensions and verify the correctness in the displayed part. The finished program was created using code language Python.