Laskentamalli PKL-koelaitteiston pienen vuodon kokeen simuloimiseksi Apros-ohjelmistolla

Työn tavoitteena oli tehdä Apros-laskentamalli PKL-koelaitteistosta ja testata kuinka hyvin Apros pystyy laskemaan PKL-koelaitteistolla suoritetun E2.2 pienen vuodon kokeen. Tavoitteena oli myös tarkastella boorittoman veden tulpan etenemistä pienen vuodon kokeen aikana. PKL-koelaitteisto vastaa saksalaista sähköteholtaan 1300 MW olevaa Philippsburg 2 painevesilaitosta. Koelaitteiston tilavuudet ja teho on skaalattu kertoimella 145. Työssä tehdyllä laskentamallilla tarkasteltiin boorittoman veden tulpan liikkumista pienen vuodon kokeen aikana. Kun malli oli valmis, laskenta suoritettiin Apros 5.05 versiolla. Boorittoman veden tulpan etenemisen laskennassa käytettiin toisen kertaluvun diskretointia, jolla booripitoisuuden muutokset säilyvät teräväreunaisina. Laskentamalli pystyi kuvaamaan koelaitteistolla suoritetussa pienen vuodon kokeessa tapahtuneet ilmiöt varsin hyvin. Eroa koetuloksiintuli pääkiertopiirien luonnonkiertojen alkamishetkistä ja primääripaineen käyttäytymisessä. Kokeen alkutilanne ei ollut stationääritila, joten alkutilanteen asettamisessa oli hankaluuksia. Varsinkin pääkiertopiirien veden pinnankorkeuksienasettamisessa oli vaikeuksia, koska veden pinnankorkeuksien erot pyrkivät tasoittumaan nopeasti kokeen aikana. Apros pystyi laskemaan PKL-koelaitteistolla suoritetun pienen vuodon kokeen hyvin. Mallilla tulisi kuitenkin laskea vielä toisentyyppisiäkin kokeita, ennen kuin voidaan varmuudella tietää mallin toimivuus. PKL-koelaitteisto vastaa pääpiirteiltään Suomeen rakennettavaa Olkiluoto 3 ydinvoimalaitosta. Tehty työ antaa lisävarmuutta, kun Olkiluoto 3 laitoksen turvallisuustarkasteluita tehdään.

The maximum voltage drop in an on-chip power distribution network: analysis of square, triangular and hexagonal power pad arrangements

A mathematical model of the voltage drop which arises in on-chip power distribution networks is used to compare the maximum voltage drop in the case of different geometric arrangements of the pads supplying power to the chip. These include the square or Manhattan power pad arrangement, which currently predominates, as well as equilateral triangular and hexagonal arrangements. In agreement with the findings in the literature and with physical and SPICE models, the equilateral triangular power pad arrangement is found to minimize the maximum voltage drop. This headline finding is a consequence of relatively simple formulas for the voltage drop, with explicit error bounds, which are established using complex analysis techniques, and elliptic functions in particular.

Direct-on-line axial flux permanent magnet synchronous generator static and dynamic performance

In distributed energy production, permanent magnet synchronous generators (PMSG) are often connected to the grid via frequency converters, such as voltage source line converters. The price of the converter may constitute a large part of the costs of a generating set. Some of the permanent magnet synchronous generators with converters and traditional separately excited synchronous generators couldbe replaced by direct-on-line (DOL) non-controlled PMSGs. Small directly networkconnected generators are likely to have large markets in the area of distributed electric energy generation. Typical prime movers could be windmills, watermills and internal combustion engines. DOL PMSGs could also be applied in island networks, such as ships and oil platforms. Also various back-up power generating systems could be carried out with DOL PMSGs. The benefits would be a lower priceof the generating set and the robustness and easy use of the system. The performance of DOL PMSGs is analyzed. The electricity distribution companies have regulations that constrain the design of the generators being connected to the grid. The general guidelines and recommendations are applied in the analysis. By analyzing the results produced by the simulation model for the permanent magnet machine, the guidelines for efficient damper winding parameters for DOL PMSGs are presented. The simulation model is used to simulate grid connections and load transients. The damper winding parameters are calculated by the finite element method (FEM) and determined from experimental measurements. Three-dimensional finite element analysis (3D FEA) is carried out. The results from the simulation model and 3D FEA are compared with practical measurements from two prototype axial flux permanent magnet generators provided with damper windings. The dimensioning of the damper winding parameters is case specific. The damper winding should be dimensioned based on the moment of inertia of the generating set. It is shown that the damper winding has optimal values to reach synchronous operation in the shortest period of time after transient operation. With optimal dimensioning, interferenceon the grid is minimized.

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.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Generating highly balanced sudoku problems as hard problems

Sudoku problems are some of the most known and enjoyed pastimes, with a never diminishing popularity, but, for the last few years those problems have gone from an entertainment to an interesting research area, a twofold interesting area, in fact. On the one side Sudoku problems, being a variant of Gerechte Designs and Latin Squares, are being actively used for experimental design, as in [8, 44, 39, 9]. On the other hand, Sudoku problems, as simple as they seem, are really hard structured combinatorial search problems, and thanks to their characteristics and behavior, they can be used as benchmark problems for refining and testing solving algorithms and approaches. Also, thanks to their high inner structure, their study can contribute more than studies of random problems to our goal of solving real-world problems and applications and understanding problem characteristics that make them hard to solve. In this work we use two techniques for solving and modeling Sudoku problems, namely, Constraint Satisfaction Problem (CSP) and Satisfiability Problem (SAT) approaches. To this effect we define the Generalized Sudoku Problem (GSP), where regions can be of rectangular shape, problems can be of any order, and solution existence is not guaranteed. With respect to the worst-case complexity, we prove that GSP with block regions of m rows and n columns with m = n is NP-complete. For studying the empirical hardness of GSP, we define a series of instance generators, that differ in the balancing level they guarantee between the constraints of the problem, by finely controlling how the holes are distributed in the cells of the GSP. Experimentally, we show that the more balanced are the constraints, the higher the complexity of solving the GSP instances, and that GSP is harder than the Quasigroup Completion Problem (QCP), a problem generalized by GSP. Finally, we provide a study of the correlation between backbone variables – variables with the same value in all the solutions of an instance– and hardness of GSP.

Generating hard SAT/CSP instances using expander graphs

In this paper we provide a new method to generate hard k-SAT instances. We incrementally construct a high girth bipartite incidence graph of the k-SAT instance. Having high girth assures high expansion for the graph, and high expansion implies high resolution width. We have extended this approach to generate hard n-ary CSP instances and we have also adapted this idea to increase the expansion of the system of linear equations used to generate XORSAT instances, being able to produce harder satisfiable instances than former generators.

Otolith shape analysis as a tool for population discrimination of the white grunt (Haemulon plumieri) stock in the northern coast of the Yucatan Peninsula, Mexico

In this study, the population structure of the white grunt (Haemulon plumieri) from the northern coast of the Yucatan Peninsula was determined through an otolith shape analysis based on the samples collected in three locations: Celestún (N 20°49",W 90°25"), Dzilam (N 21°23", W 88°54") and Cancún (N 21°21",W 86°52"). The otolith outline was based on the elliptic Fourier descriptors, which indicated that the H. plumieri population in the northern coast of the Yucatan Peninsula is composed of three geographically delimited units (Celestún, Dzilam, and Cancún). Significant differences were observed in mean otolith shapes among all samples (PERMANOVA; F2, 99 = 11.20, P = 0.0002), and the subsequent pairwise comparisons showed that all samples were significantly differently from each other. Samples do not belong to a unique white grunt population, and results suggest that they might represent a structured population along the northern coast of the Yucatan Peninsula

A parameterization method for the computation of invariant tori and their whiskers in quasi-periodic maps: explorations and mechanisms for the breakdown of hyperbolicity

In two previous papers [J. Differential Equations, 228 (2006), pp. 530 579; Discrete Contin. Dyn. Syst. Ser. B, 6 (2006), pp. 1261 1300] we have developed fast algorithms for the computations of invariant tori in quasi‐periodic systems and developed theorems that assess their accuracy. In this paper, we study the results of implementing these algorithms and study their performance in actual implementations. More importantly, we note that, due to the speed of the algorithms and the theoretical developments about their reliability, we can compute with confidence invariant objects close to the breakdown of their hyperbolicity properties. This allows us to identify a mechanism of loss of hyperbolicity and measure some of its quantitative regularities. We find that some systems lose hyperbolicity because the stable and unstable bundles approach each other but the Lyapunov multipliers remain away from 1. We find empirically that, close to the breakdown, the distances between the invariant bundles and the Lyapunov multipliers which are natural measures of hyperbolicity depend on the parameters, with power laws with universal exponents. We also observe that, even if the rigorous justifications in [J. Differential Equations, 228 (2006), pp. 530-579] are developed only for hyperbolic tori, the algorithms work also for elliptic tori in Hamiltonian systems. We can continue these tori and also compute some bifurcations at resonance which may lead to the existence of hyperbolic tori with nonorientable bundles. We compute manifolds tangent to nonorientable bundles.

Right ventricular dysfunction in children and adolescents conceived by assisted reproductive technologies.

Assisted reproductive technologies (ART) predispose the offspring to vascular dysfunction, arterial hypertension, and hypoxic pulmonary hypertension. Recently, cardiac remodeling and dysfunction during fetal and early postnatal life have been reported in offspring of ART, but it is not known whether these cardiac alterations persist later in life and whether confounding factors contribute to this problem. We, therefore, assessed cardiac function and pulmonary artery pressure by echocardiography in 54 healthy children conceived by ART (mean age 11.5 ± 2.4 yr) and 54 age-matched (12.2 ± 2.3 yr) and sex-matched control children. Because ART is often associated with low birth weight and prematurity, two potential confounders associated with cardiac dysfunction, only singletons born with normal birth weight at term were studied. Moreover, because cardiac remodeling in infants conceived by ART was observed in utero, a situation associated with increased right heart load, we also assessed cardiac function during high-altitude exposure, a condition associated with hypoxic pulmonary hypertension-induced right ventricular overload. We found that, while at low altitude cardiac morphometry and function was not different between children conceived by ART and control children, under the stressful conditions of high-altitude-induced pressure overload and hypoxia, larger right ventricular end-diastolic area and diastolic dysfunction (evidenced by lower E-wave tissue Doppler velocity and A-wave tissue Doppler velocity of the lateral tricuspid annulus) were detectable in children and adolescents conceived by ART. In conclusion, right ventricular dysfunction persists in children and adolescents conceived by ART. These cardiac alterations appear to be related to ART per se rather than to low birth weight or prematurity.

Kuumalanka-anemometri virtausmittauksissa

Tässä työssä selostetaan kuumalanka-anemometrin käyttö virtausmittauksissa. Kuumalanka-anemometrilla saadaan mitattua virtausnopeuden ja -suunnan lisäksi nopeusheilahteluja. Mittaustaajuus on tyypillisesti useita kymmeniä tuhansia mittauksia sekunnissa ja signaali on jatkuva. Nykytekniikalla pystytään helposti tallentamaan mittauslaitteistolta saatu viesti tietokoneelle ja muuntamaan se nopeudeksi. Hetkellisten nopeuksien avulla voidaan laskea turbulenttisen virtauksen ominaisuuksia, kuten turbulenssin intensiteetti ja spektri. Kuumalanka-anemometrissa lämmitetään sähköisesti ohutta lankaa, joka on mitattavassa virtauksessa. Langan sähköteho on suunnilleen yhtäsuuri kuin langasta konvektiolla siirtyvä lämpöteho. Tällöin on teoreettisesti mahdollista laskea virtausnopeus lämpötehosta lämmönsiirtokorrelaatioilla. Käytännössä laitteisto joudutaan kuitenkin erikseen kalibroimaan, mutta sähkötehon teoreettista riippuvuutta konvektiosta käytetään hyväksi. Kuumalangan lämmitettävä osuus on tyypillisesti halkaisijaltaan 5 µm ja pituudeltaan noin 1 mm. Sitä käytetään pääasiassa kaasuvirtausten mittaamiseen ja valtaosassa mittauksissa virtausaineena on ilma. Kuumalanka voi olla toteutettu kuumakalvotekniikalla, jossa halkaisijaltaan noin 50 - 70 µm paksuinen kuitu on päällystetty ohuella sähköä johtavalla kalvolla. Kuumakalvoanturin ei tarvitse olla muodoltaan sylinterimäinen, se voi olla mm. kartiomainen tai kiilamainen. Erikoispäällystetyllä kuumakalvoanturilla on mahdollista mitata myös nestevirtauksia. Mitattaessa kaasuvirtauksia kuumakalvon etuna on selvästi parempi kestävyys verrattuna kuumalankaan. Nimitystä kuumalanka-anemometri käytetään yleisesti molemmista anturityypeistä Tämän työn alussa käsitellään sylinterin yli tapahtuvaan virtaukseen liittyvää virtausmekaniikkaa ja lämmönsiirtoa. Anemometrin sähköinen osa, laitteisto ja sen kalibrointi käydään läpi. Langan suuntariippuvuuden laskentaan esitetään tarvittavat yhtälöt. Työssä esitellään kolme laitteistolla tehtyä perusmittausta: anturin kohtauskulman muuttaminen, pyörähdyssymmerisen suihkun nopeuskenttä ja tuulitunnelin rajakerros. Lisäksi esitellään yksi käytännöllinen ja vaativampi mittaus, jossa on mitattu nopeusprofiili radiaalikompressorin diffuusorin loppuosassa.

Teollisuusverkon stabiilisuus saarekekäyttöön siirryttäessä

Työn aluksi on käsitelty sähkövoimajärjestelmien stabiilisuutta yleisesti. Myös teollisuuden sähköverkkojen erityispiirteitä on tarkasteltu. Teollisuusverkkojen komponentteja koskevassa selvityksessä eri laitteet on käsitelty siten, että niiden tärkeimmät ominaisuudet verkon erilaisissa muutostilanteissa tulevat esille. Tahtigeneraattorit ja epätahtikoneet on käsitelty yksityiskohtaisesti. Tutkimuksessa on esitelty teollisuusverkon saarekekäyttöön johtavia syitä, saarekkeeseen siirtymisen kriteereitä ja saarekekäytön toteutustapoja. Tehonvajaussuojaukseen liittyen on tutustuttu kuormanpudotusjärjestelyihin. Myös höyryverkon kuormien tasausten mahdollisuuksiin on luotu katsaus. Verkostolaskentaa käsittelevässä luvussa on keskitytty sopivan laskentamallin luomiseen liittyviin menetelmiin ja rajoituksiin. Lisäksi on kerrottu työn käytännön osan laskentamallin luomisen keskeiset periaatteet. Työn lopuksi on suoritettu simulointeja, joissa on tutkittu erään teollisuusverkon stabiilisuutta saarekekäyttöön siirryttäessä Calpos-verkostolaskentaohjelmistolla.

Characteristics of 2D convective structures in Catalonia (NE Spain): an analysis using radar data and GIS

Flood simulation studies use spatial-temporal rainfall data input into distributed hydrological models. A correct description of rainfall in space and in time contributes to improvements on hydrological modelling and design. This work is focused on the analysis of 2-D convective structures (rain cells), whose contribution is especially significant in most flood events. The objective of this paper is to provide statistical descriptors and distribution functions for convective structure characteristics of precipitation systems producing floods in Catalonia (NE Spain). To achieve this purpose heavy rainfall events recorded between 1996 and 2000 have been analysed. By means of weather radar, and applying 2-D radar algorithms a distinction between convective and stratiform precipitation is made. These data are introduced and analyzed with a GIS. In a first step different groups of connected pixels with convective precipitation are identified. Only convective structures with an area greater than 32 km2 are selected. Then, geometric characteristics (area, perimeter, orientation and dimensions of the ellipse), and rainfall statistics (maximum, mean, minimum, range, standard deviation, and sum) of these structures are obtained and stored in a database. Finally, descriptive statistics for selected characteristics are calculated and statistical distributions are fitted to the observed frequency distributions. Statistical analyses reveal that the Generalized Pareto distribution for the area and the Generalized Extreme Value distribution for the perimeter, dimensions, orientation and mean areal precipitation are the statistical distributions that best fit the observed ones of these parameters. The statistical descriptors and the probability distribution functions obtained are of direct use as an input in spatial rainfall generators.

Otsonin käyttö paperin ja kartongin valmistuksen kiertovesien mikrobiologisen tilan hallintaan

Diplomityössä tutkittiin otsonin käytön soveltuvuutta ja etsittiin optimaalista soveltamistapaa kartonki- ja paperikoneen kiertoveden mikrobiologisen tilan hallintaan. Kirjallisuusosassa selvitettiin otsonin tuotannon nykytekniikkaa sekä otettiin selvää saatavilla olevista kapasiteetiltaan riittävistä otsonigeneraattoreista. Kirjallisuuden perusteella selvitettiin myös otsonin reaktioita paperi- ja kartonkikoneen kiertovesissä olevien orgaanisten ja epäorgaanisten aineiden kanssa. Lisäksi käsiteltiin otsonin käytön eri sovelluksia vedenpuhdistuksessa sekä perehdyttiin aikaisemmin tehtyihin tutkimuksiin otsonin käytöstä mikrobitorjuntaan paperi- ja kartonkikoneiden kiertovesistä. Kokeellisessa osassa määritettiin sopiva otsoniannos valkaisematonta sellua käyttävän koneen kiekkosuotimen kirkkaalle ja samealle suodokselle. Vesien otsonointi suoritettiin otsonipilotlaitteistolla. Työssä tutkittiin otsonoinnin vaikutusta kiertoveden mikrobiologisiin ja kemiallisiin ominaisuuksiin. Lisäksi tarkasteltiin kiertoveden otsonoinnin vaikutuksia koepaperikoneella valkaisemattomasta sel-lusta ja hylkymassasta valmistettujen lopputuotteiden laatuun. Työssä arvioitiin myös otsonin käytön mahdollisuutta osittain korvaamaan jo käytössä olevia hapettavia biosidejä. Samealle suodokselle sopivaksi otsoniannokseksi osoittautui 120 mg O3/l (0,3 mg O3/mg CODCr) ja kirkkaalle suodokselle riittävä annos olisi 60 mg O3/l (0,2 mg/mg CODCr). Näillä annoksilla kiertoveden mikrobireduktiot olivat 98–100 %. Kiertoveden hapetuskyky kasvoi 9-14 % ja pH laski yhdellä yksiköllä otsonoinnin vaikutuksesta. Kiertoveden kemiallisiin ominaisuuksiin otsonoinnilla oli edullinen vaikutus, sillä veden sameus, kiintoaine- ja CODCr-pitoisuudet laskivat. Koepaperikoneella havaittiin perälaatikon ilmapitoisuuden vähentyneen jopa 49 % käytettäessä ot-sonoitua kiertovettä. Mikrobiologiset analyysit suoritettiin koepaperikoneen massalle, viiravedelle ja lopputuotteelle. Tutkimusten perusteella mikrobientorjunnan kannalta parhaimmaksi otsonin annos-telupaikaksi osoittautui samea suodos. Massan, viiraveden ja lopputuotteen mikrobimäärät alenivat selvästi. Kemialliset analyysit osoittivat, että kirkasta suodosta käytettäessä oli koepaperikoneen vii-raveden kemiallinen tila parantunut LK-aineen muuttuessa hydrofiilisemmäksi ja sen sitoutumisky-vyn parantuessa. Samealla suodoksella ei muutosta havaittu. Valkaisemattomasta sellusta ja päällystetystä hylystä koepaperikoneella valmistettujen lopputuotteiden lujuus- ja aistinvaraisissa ominaisuuksissa ei tapahtunut merkittävää muutosta. Tulosten perusteella otsonia voisi käyttää korvaamaan hinnaltaan kalliimpia hapettavia biosidejä. Otsonin etuna on, ettei sille resistenttejä mikrobikantoja ole. Kirjallisuuden perusteella otsonin heikkoutena on sen tehon paikallisuus eli otsonin lisänä on käytettävä jotain muuta sekundääristä biosidia ylläpitämään matala mikrobitaso otsonoinnin jälkeen.

Draft genome sequence of Aeromonas molluscorum strain 848TT, isolated from bivalve molluscs

We report here the draft genome sequence of Aeromonas molluscorum 848T, the type strain of this Aeromonas species, which was isolated from wedge shells (Donax trunculus) obtained from a retail market in Barcelona, Spain, in 1997.