Age-related differences on event-related potentials and brain rhythm oscillations during working memory activation.

Previous functional imaging studies have pointed to the compensatory recruitment of cortical circuits in old age in order to counterbalance the loss of neural efficiency and preserve cognitive performance. Recent electroencephalographic (EEG) analyses reported age-related deficits in the amplitude of an early positive-negative working memory (PN(wm)) component as well as changes in working memory (WM)-load related brain oscillations during the successful performance of the n-back task. To explore the age-related differences of EEG activation in the face of increasing WM demands, we assessed the PN(wm) component area, parietal alpha event-related synchronization (ERS) as well as frontal theta ERS in 32 young and 32 elderly healthy individuals who successfully performed a highly WM demanding 3-back task. PN(wm) area increased with higher memory loads (3- and 2-back > 0-back tasks) in younger subjects. Older subjects reached the maximal values for this EEG parameter during the less WM demanding 0-back task. They showed a rapid development of an alpha ERS that reached its maximal amplitude at around 800 ms after stimulus onset. In younger subjects, the late alpha ERS occurred between 1,200 and 2,000 ms and its amplitude was significantly higher compared with elders. Frontal theta ERS culmination peak decreased in a task-independent manner in older compared with younger cases. Only in younger individuals, there was a significant decrease in the phasic frontal theta ERS amplitude in the 2- and 3-back tasks compared with the detection and 0-back tasks. These observations suggest that older adults display a rapid mobilization of their neural generators within the parietal cortex to manage very low demanding WM tasks. Moreover, they are less able to activate frontal theta generators during attentional tasks compared with younger persons.

The timing of exploratory decision-making revealed by single-trial topographic EEGanalyses.

Decision-making in an uncertain environment is driven by two major needs: exploring the environment to gather information or exploiting acquired knowledge to maximize reward. The neural processes underlying exploratory decision-making have been mainly studied by means of functional magnetic resonance imaging, overlooking any information about the time when decisions are made. Here, we carried out an electroencephalography (EEG) experiment, in order to detect the time when the brain generators responsible for these decisions have been sufficiently activated to lead to the next decision. Our analyses, based on a classification scheme, extract time-unlocked voltage topographies during reward presentation and use them to predict the type of decisions made on the subsequent trial. Classification accuracy, measured as the area under the Receiver Operator's Characteristic curve was on average 0.65 across 7 subjects. Classification accuracy was above chance levels already after 516 ms on average, across subjects. We speculate that decisions were already made before this critical period, as confirmed by a positive correlation with reaction times across subjects. On an individual subject basis, distributed source estimations were performed on the extracted topographies to statistically evaluate the neural correlates of decision-making. For trials leading to exploration, there was significantly higher activity in dorsolateral prefrontal cortex and the right supramarginal gyrus; areas responsible for modulating behavior under risk and deduction. No area was more active during exploitation. We show for the first time the temporal evolution of differential patterns of brain activation in an exploratory decision-making task on a single-trial basis.

Early auditory processing

ABSTRACT (English)An accurate processing of the order between sensory events at the millisecond time scale is crucial for both sensori-motor and cognitive functions. Temporal order judgment (TOJ) tasks, is the ability of discriminating the order of presentation of several stimuli presented in a rapid succession. The aim of the present thesis is to further investigate the spatio-temporal brain mechanisms supporting TOJ. In three studies we focus on the dependency of TOJ accuracy on the brain states preceding the presentation of TOJ stimuli, the neural correlates of accurate vs. inaccurate TOJ and whether and how TOJ performance can be improved with training.In "Pre-stimulus beta oscillations within left posterior sylvian regions impact auditory temporal order judgment accuracy" (Bernasconi et al., 2011), we investigated if the brain activity immediately preceding the presentation of the stimuli modulates TOJ performance. By contrasting the electrophysiological activity before the stimulus presentation as a function of TOJ accuracy we observed a stronger pre-stimulus beta (20Hz) oscillatory activity within the left posterior sylvian region (PSR) before accurate than inaccurate TOJ trials.In "Interhemispheric coupling between the posterior sylvian regions impacts successful auditory temporal order judgment" (Bernasconi et al., 2010a), and "Plastic brain mechanisms for attaining auditory temporal order judgment proficiency" (Bernasconi et al., 2010b), we investigated the spatio-temporal brain dynamics underlying auditory TOJ. In both studies we observed a topographic modulation as a function of TOJ performance at ~40ms after the onset of the first sound, indicating the engagement of distinct configurations of intracranial generators. Source estimations in the first study revealed a bilateral PSR activity for both accurate and inaccurate TOJ trials. Moreover, activity within left, but not right, PSR correlated with TOJ performance. Source estimations in the second study revealed a training-induced left lateralization of the initial bilateral (i.e. PSR) brain response. Moreover, the activity within the left PSR region correlated with TOJ performance.Based on these results, we suggest that a "temporal stamp" is established within left PSR on the first sound within the pair at early stages (i.e. ~40ms) of cortical processes, but is critically modulated by inputs from right PSR (Bernasconi et al., 2010a; b). The "temporal stamp" on the first sound may be established via a sensory gating or prior entry mechanism.Behavioral and brain responses to identical stimuli can vary due to attention modulation, vary with experimental and task parameters or "internal noise". In a fourth experiment (Bernasconi et al., 2011b) we investigated where and when "neural noise" manifest during the stimulus processing. Contrasting the AEPs of identical sound perceived as High vs. Low pitch, a topographic modulation occurred at ca. 100ms after the onset of the sound. Source estimation revealed activity within regions compatible with pitch discrimination. Thus, we provided neurophysiological evidence for the variation in perception induced by "neural noise".ABSTRACT (French)Un traitement précis de l'ordre des événements sensoriels sur une échelle de temps de milliseconde est crucial pour les fonctions sensori-motrices et cognitives. Les tâches de jugement d'ordre temporel (JOT), consistant à présenter plusieurs stimuli en succession rapide, sont traditionnellement employées pour étudier les mécanismes neuronaux soutenant le traitement d'informations sensorielles qui varient rapidement. Le but de cette thèse est d'étudier le mécanisme cérébral soutenant JOT. Dans les trois études présentées nous nous sommes concentrés sur les états du cerveau précédant la présentation des stimuli de JOT, les bases neurales pour un JOT correct vs. incorrect et sur la possibilité et les moyens d'améliorer l'exécution du JOT grâce à un entraînement.Dans "Pre-stimulus beta oscillations within left posterior sylvian regions impact auditory temporal order judgment accuracy" (Bernasconi et al., 2011),, nous nous sommes intéressé à savoir si l'activité oscillatoire du cerveau au pré-stimulus modulait la performance du JOT. Nous avons contrasté l'activité électrophysiologique en fonction de la performance TOJ, mesurant une activité oscillatoire beta au pré-stimulus plus fort dans la région sylvian postérieure gauche (PSR) liée à un JOT correct.Dans "Interhemispheric coupling between the posterior sylvian regions impacts successful auditory temporal order judgment" (Bernasconi et al., 2010a), et "Plastic brain mechanisms for attaining auditory temporal order judgment proficiency" (Bernasconi et al., 2010b), nous avons étudié la dynamique spatio-temporelle dans le cerveau impliqué dans le traitement du JOT auditif. Dans ses deux études, nous avons observé une modulation topographique à ~40ms après le début du premier son, en fonction de la performance JOT, indiquant l'engagement des configurations de générateurs intra- crâniens distincts. La localisation de source dans la première étude indique une activité bilatérale de PSR pour des JOT corrects vs. incorrects. Par ailleurs, l'activité dans PSR gauche, mais pas dans le droit, est corrélée avec la performance du JOT. La localisation de source dans la deuxième étude indiquait une latéralisation gauche induite par l'entraînement d'une réponse initialement bilatérale du cerveau. D'ailleurs, l'activité dans la région PSR gauche corrèlait avec la performance de TOJ.Basé sur ces résultats, nous proposons qu'un « timbre-temporel » soit établi très tôt (c.-à-d. à ~40ms) sur le premier son par le PSR gauche, mais module par l'activité du PSR droite (Bernasconi et al., 2010a ; b). « Le timbre- temporel » sur le premier son peut être établi par le mécanisme neuronal de type « sensory gating » ou « prior entry ».Les réponses comportementales et du cerveau aux stimuli identiques peut varier du à des modulations d'attention ou à des variations dans les paramètres des tâches ou au bruit interne du cerveau. Dans une quatrième expérience (Bernasconi et al. 2011B), nous avons étudié où et quand le »bruit neuronal« se manifeste pendant le traitement des stimuli. En contrastant les AEPs de sons identiques perçus comme aigus vs. grave, nous avons mesuré une modulation topographique à env. 100ms après l'apparition du son. L'estimation de source a révélé une activité dans les régions compatibles avec la discrimination de fréquences. Ainsi, nous avons fourni des preuves neurophysiologiques de la variation de la perception induite par le «bruit neuronal».

Primääri-sekundäärivuodon tunnistaminen aktiivisuusmittauksin Loviisan voimalaitoksella

Diplomityössä tutkittiin Loviisan voimalaitoksen primääri- ja sekundääripiirin aktiivisuusmittausten kykyä tunnistaa pienet primääri-sekundäärivuodot. Tarkasteltavat primääri-sekundäärivuotojen suuruudet valittiin laitoksen hätätilanne- ja häiriönselvitysohjeiden mukaisesti. Vuodon vaikutuksia arvioitiin erilaisilla primäärijäähdytteen ominaisaktiivisuuksilla. Ominaisaktiivisuudet primääripiirissä määritettiin nuklidikohtaisesti erilaisille polttoainevuototapauksille. Työssä huomioitiin myös transienteissa mahdollisesti esiintyvä primääripiirin aktiivisuustasoa nostava spiking-ilmiö. Vuodon tarkempaa tunnistamista varten työssä laskettiin tarkasteltaville mittareille kalibrointikertoimet. Primääri-sekundäärivuoto mallinnettiin APROS-simulointiohjelmalla laitoksen eri käyttötiloissa ja kahdella eri vuotokoolla. Varsinainen aktiivisuuslaskenta suoritettiin SEKUN-ohjelmalla. Työssä tätä aktiivisuus- ja päästölaskentaohjelmaa muokattiin ohjelmoimalla siihen tarkasteltavat aktiivisuusmittaukset sekä primääripiirin puhdistus ja ulospuhallus. Laskelmien tuloksena saatiin arviot kunkin tarkasteltavana olleen aktiivisuusmittauksen soveltuvuudesta primääri-sekundäärivuodon tunnistamiseen erilaisissa polttoainevuototapauksissa ja reaktorin eri tehotasoilla. Häiriönselvitysohje I3:n käyttöönottoa edellyttävät vuotokoot määritettiin aktiivisuusmittausten havaitseman perusteella. Erityisesti kuumavalmiustilassa tapauksissa, joissa reaktorisydämessä oletetaan olevan tiiveytensä menettäneitä polttoainesauvoja, spikingin vaikutus jäähdytteiden aktiivisuuspitoisuuksiin ja mittaustuloksiin oli merkittävä. Niiltä osin, kuin tulokset käsittelevät ohjeissa vuodon tunnistamiseen käytettyjä aktiivisuusrajoja, tulokset osoittivat aktiivisuusrajat oikeiksi. Kuumavalmiudessa aktiivisuusmittausten mittausalueet saattavat joissakin tapauksissa rajoittaa primääri-sekundäärivuodon tunnistamista.

Symmetries of the free Schrödinger equation in the non-commutative plane

We study all the symmetries of the free Schr odinger equation in the non-commu- tative plane. These symmetry transformations form an infinite-dimensional Weyl algebra that appears naturally from a two-dimensional Heisenberg algebra generated by Galilean boosts and momenta. These infinite high symmetries could be useful for constructing non-relativistic interacting higher spin theories. A finite-dimensional subalgebra is given by the Schröodinger algebra which, besides the Galilei generators, contains also the dilatation and the expansion. We consider the quantization of the symmetry generators in both the reduced and extended phase spaces, and discuss the relation between both approaches.

Symmetries of the free Schrödinger equation in the non-commutative plane

We study all the symmetries of the free Schrödinger equation in the non-commu- tative plane. These symmetry transformations form an infinite-dimensional Weyl algebra that appears naturally from a two-dimensional Heisenberg algebra generated by Galilean boosts and momenta. These infinite high symmetries could be useful for constructing non-relativistic interacting higher spin theories. A finite-dimensional subalgebra is given by the Schröodinger algebra which, besides the Galilei generators, contains also the dilatation and the expansion. We consider the quantization of the symmetry generators in both the reduced and extended phase spaces, and discuss the relation between both approaches.

Pinnan ominaisuuksien karakterisointi kiiltomittalaitteella ja kiiltomittalaitteen soveltuvuus tuotannon laadunvarmistuksessa

Alahuulten, kärkilistojen, turboputkien ja telapinnoitteiden pinnankarheutta mitataan Rautpohjan tehtaalla nykyään pääosin koskettavalla pinnankarheusprofilometrilla. Profilometrin käyttö on hankalaa ja epäluotettavaa tai käytössä on muita ongelmia. Paperin kuituorientaation mittaamiseen käytetään nykyisin mm. TSO-menetelmää tai repäisymenetelmää. Tämän diplomityön tarkoituksena oli kartoittaa diffraktiiviseen optiseen elementtiin perustuvan kiiltomittalaiteprototyypin mahdollisuuksia karakteroida pinnan ominaisuuksia ja tutkia kiiltomittalaitteen soveltuvuutta tuotannon laadunvarmistuksessa. Alahuuli-, kärkilista-, turboputki-, telapinnoite- ja paperinäytteistä valmistettiinnäytesarjoja, jotka mitattiin nykymenetelmin ja kiiltomittalaitteella. Tulostenkorreloituvuutta tutkittiin ja taloudellista merkitystä arvioitiin investointilaskelmin. Paperin kuituorientaation, paperin pinnan kuituorientaation ja telapinnoitteiden mittaaminen kiiltomittalaitteella osoittautui taloudellisesti erittäin merkittäväksi. Ennen mahdollista kiiltomittalaitteen käyttöönottoa prototyyppiä on kehitettävä, kiillon korreloituvuutta tulee tutkia edelleen ja aiheesta on tehtävä lisätutkimuksia.

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.

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.

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.