Statistical Inversion Theory in X-ray Tomography

In this thesis the X-ray tomography is discussed from the Bayesian statistical viewpoint. The unknown parameters are assumed random variables and as opposite to traditional methods the solution is obtained as a large sample of the distribution of all possible solutions. As an introduction to tomography an inversion formula for Radon transform is presented on a plane. The vastly used filtered backprojection algorithm is derived. The traditional regularization methods are presented sufficiently to ground the Bayesian approach. The measurements are foton counts at the detector pixels. Thus the assumption of a Poisson distributed measurement error is justified. Often the error is assumed Gaussian, altough the electronic noise caused by the measurement device can change the error structure. The assumption of Gaussian measurement error is discussed. In the thesis the use of different prior distributions in X-ray tomography is discussed. Especially in severely ill-posed problems the use of a suitable prior is the main part of the whole solution process. In the empirical part the presented prior distributions are tested using simulated measurements. The effect of different prior distributions produce are shown in the empirical part of the thesis. The use of prior is shown obligatory in case of severely ill-posed problem.

Inversion and assimilation methods with applications in geophysical remote sensing

Stratospheric ozone can be measured accurately using a limb scatter remote sensing technique at the UV-visible spectral region of solar light. The advantages of this technique includes a good vertical resolution and a good daytime coverage of the measurements. In addition to ozone, UV-visible limb scatter measurements contain information about NO2, NO3, OClO, BrO and aerosols. There are currently several satellite instruments continuously scanning the atmosphere and measuring the UVvisible region of the spectrum, e.g., the Optical Spectrograph and Infrared Imager System (OSIRIS) launched on the Odin satellite in February 2001, and the Scanning Imaging Absorption SpectroMeter for Atmospheric CartograpHY (SCIAMACHY) launched on Envisat in March 2002. Envisat also carries the Global Ozone Monitoring by Occultation of Stars (GOMOS) instrument, which also measures limb-scattered sunlight under bright limb occultation conditions. These conditions occur during daytime occultation measurements. The global coverage of the satellite measurements is far better than any other ozone measurement technique, but still the measurements are sparse in the spatial domain. Measurements are also repeated relatively rarely over a certain area, and the composition of the Earth’s atmosphere changes dynamically. Assimilation methods are therefore needed in order to combine the information of the measurements with the atmospheric model. In recent years, the focus of assimilation algorithm research has turned towards filtering methods. The traditional Extended Kalman filter (EKF) method takes into account not only the uncertainty of the measurements, but also the uncertainty of the evolution model of the system. However, the computational cost of full blown EKF increases rapidly as the number of the model parameters increases. Therefore the EKF method cannot be applied directly to the stratospheric ozone assimilation problem. The work in this thesis is devoted to the development of inversion methods for satellite instruments and the development of assimilation methods used with atmospheric models.

Fractional slot permanent magnet synchronous motors for low speed applications

This study compares different rotor structures of permanent magnet motors with fractional slot windings. The surface mounted magnet and the embedded magnet rotor structures are studied. This thesis analyses the characteristics of a concentrated two-layer winding, each coil of which is wound around one tooth and which has a number of slots per pole and per phase less than one (q < 1). Compared to the integer slot winding, the fractional winding (q < 1) has shorter end windings and this, thereby, makes space as well as manufacturing cost saving possible. Several possible ways of winding a fractional slot machine with slots per pole and per phase lessthan one are examined. The winding factor and the winding harmonic components are calculated. The benefits attainable from a machine with concentrated windingsare considered. Rotor structures with surface magnets, radially embedded magnets and embedded magnets in V-position are discussed. The finite element method isused to solve the main values of the motors. The waveform of the induced electro motive force, the no-load and rated load torque ripple as well as the dynamic behavior of the current driven and voltage driven motor are solved. The results obtained from different finite element analyses are given. A simple analytic method to calculate fractional slot machines is introduced and the values are compared to the values obtained with the finite element analysis. Several different fractional slot machines are first designed by using the simple analytical methodand then computed by using the finite element method. All the motors are of thesame 225-frame size, and have an approximately same amount of magnet material, a same rated torque demand and a 400 - 420 rpm speed. An analysis of the computation results gives new information on the character of fractional slot machines.A fractional slot prototype machine with number 0.4 for the slots per pole and per phase, 45 kW output power and 420 rpm speed is constructed to verify the calculations. The measurement and the finite element method results are found to beequal.

Development of advanced silicon radiation detectors for harsh radiation environment

This thesis describes the development of advanced silicon radiation detectors and their characterization by simulations, used in the work for searching elementary particles in the European Organization for Nuclear Research, CERN. Silicon particle detectors will face extremely harsh radiation in the proposed upgrade of the Large Hadron Collider, the future high-energy physics experiment Super-LHC. The increase in the maximal fluence and the beam luminosity up to 1016 neq / cm2 and 1035 cm-2s-1 will require detectors with a dramatic improvement in radiation hardness, when such a fluence will be far beyond the operational limits of the present silicon detectors. The main goals of detector development concentrate on minimizing the radiation degradation. This study contributes mainly to the device engineering technology for developing more radiation hard particle detectors with better characteristics. Also the defect engineering technology is discussed. In the nearest region of the beam in Super-LHC, the only detector choice is 3D detectors, or alternatively replacing other types of detectors every two years. The interest in the 3D silicon detectors is continuously growing because of their many advantages as compared to conventional planar detectors: the devices can be fully depleted at low bias voltages, the speed of the charge collection is high, and the collection distances are about one order of magnitude less than those of planar technology strip and pixel detectors with electrodes limited to the detector surface. Also the 3D detectors exhibit high radiation tolerance, and thus the ability of the silicon detectors to operate after irradiation is increased. Two parameters, full depletion voltage and electric field distribution, is discussed in more detail in this study. The full depletion of the detector is important because the only depleted area in the detector is active for the particle tracking. Similarly, the high electric field in the detector makes the detector volume sensitive, while low-field areas are non-sensitive to particles. This study shows the simulation results of full depletion voltage and the electric field distribution for the various types of 3D detectors. First, the 3D detector with the n-type substrate and partial-penetrating p-type electrodes are researched. A detector of this type has a low electric field on the pixel side and it suffers from type inversion. Next, the substrate is changed to p-type and the detectors having electrodes with one doping type and the dual doping type are examined. The electric field profile in a dual-column 3D Si detector is more uniform than that in the single-type column 3D detector. The dual-column detectors are the best in radiation hardness because of their low depletion voltages and short drift distances.

Kuormituksen vaikutus kromatografiseen erotukseen alikriittisellä vedellä

Teollisessa kromatografiassa kolonnia pyritään kuormittamaan mahdollisimman paljon, jotta saataisiin maksimoitua erotetun komponentin määrä aikayksikköä kohden. Tässä työssä kuormitusta tutkittiin nostamalla syöttöliuoksen, synteettisen melassin, näyteväkevyyttä 80-125 ºC:ssa. Eluenttina oli paineistettu kuumaa vesi ja hartsina vahva Na-muotoinen PS-DVB pohjainen vahva kationinvaihtohartsi. Lämpötilaa nostamalla piikit kapenivat ja tulivat symmetrisemmiksi, erotus nopeutui sekä suola erottui usein paremmin sokereista. Syöttöliuoksen kuiva-ainetta lisättiin asteittain 55 p-% saakka, jolloin ei vielä havaittu ongelmia erotuksessa. Lämpötilassa 125 ºC havaittiin erotuksen aikana kuormituksesta riippumatonta sakkaroosin invertoitumista. Vertailtaessa eri stationäärifaaseja havaittiin Na-muotoisen PS-DVB pohjaisen kationinvaihtohartsin erottavan yleensä sokereita, sokerialkoholeja, oligosakkarideja ja betaiinia lähes poikkeuksetta paremmin alhaisilla pitoisuuksilla kuin neutraalihartsi ja Na-muotoinen zeoliitti. Erottuminen ei yleensä parantunut lämpötilaa nostamalla, mutta piikit kapenivat ja erotus nopeutui. Monosakkaridien erotus huononi 125 ºC:ssa kationinvaihtohartsilla. Tutkittaessa terveysvaikutteisten ksylo-oligosakkaridien soveltuvuutta alikriittiseen erotukseen, niiden havaittiin huomattavasti hydrolysoituvan happamissa olosuhteissa koeputkessa 100 ºC:ssa kahdessa tunnissa. Näytteessä olevien epäpuhtauksien havaittiin katalysoineen hydrolyysiä. Hydrolysoituminen oli hitaampaa neutraaleissa olosuhteissa korotetussa lämpötilassa. Tästä voitiin tehdä johtopäätös, että alikriittiset olosuhteet eivät sovi ksylo-oligosakkaridien erotukseen.

DTC-induktiomoottorikäytön soveltaminen keskitaajuusalueelle

Työssä tarkastellaan olemassa olevan suoraa vääntömomentin säätöä käyttävän taajuudenmuuttajan soveltamista keskitaajuusalueella toimiviin induktiomoottorikäyttöihin. Keskinopeusalueen sovellusten pyörimisnopeudet ovat tyypillisesti 6000…30000rpm. Tällöin invertterin lähtötaajuuden on nelinapaista moottoria ohjattaessa ulotuttava 1000Hz:iin. ABB:n ACS600 taajuudenmuuttajan nykyinen syöttötaajuus ulottuu noin 400Hz:iin ja sen keskimääräinen kytkentätaajuus on luokkaa 3kHz. Taajuudenmuuttajan keskimääräistä kytkentätaajuutta ei haluta tästä merkittävästi nostaa, koska tällöin pääteasteen mitoitusta ja rakennetta on muutettava. Tarkastelussa keskitytään täten jännitemodulointiin, joka määrittelee invertterin lähtöjännitteen käyrämuodon ja pääteasteen tehokytkimien kytkentätaajuuden. Työssä esitetään suoran käämivuon säädön periaatetta soveltava 30-kulmainen modulointimenetelmä (30-modulointi), jolla moottorin syöttövirrasta voidaan eliminoida 5. ja 7. yliharmoninen komponentti. Yliharmonisten komponenttien eliminointi mahdollistaa passiivisilla komponenteilla toteutetun alipäästösuodattimen asentamisen invertterin lähtöön, jolloin moottorissa tapahtuvat tehohäviöt saadaan pieniksi. 30-moduloinnin ohjaukseen ja säätöön esitetään menetelmät, jotka mahdollistavat sen toteuttamisen nykyiseen taajuudenmuuttajaan ohjelmallisesti ilman suuria laitteistomuutoksia. 30-moduloinnin ominaisuuksia tarkastellaan analyyttisin menetelmin ja sen toimintaa testataan esitettyjen teorioiden perusteella simuloimalla. 30-moduloinnin ohjelmallinen implementointi nykyiseen ACS600:een mahdollistaa periaatteessa noin 600Hz:n syöttötaajuuden saavuttamisen. Tällöin invertterin keskimääräinen kytkentätaajuus voidaan säätää koko pyörimisnopeusalueella alle 4kHz:n tasolle.

Adaptive MCMC methods with applications in environmental and geophysical models

This work presents new, efficient Markov chain Monte Carlo (MCMC) simulation methods for statistical analysis in various modelling applications. When using MCMC methods, the model is simulated repeatedly to explore the probability distribution describing the uncertainties in model parameters and predictions. In adaptive MCMC methods based on the Metropolis-Hastings algorithm, the proposal distribution needed by the algorithm learns from the target distribution as the simulation proceeds. Adaptive MCMC methods have been subject of intensive research lately, as they open a way for essentially easier use of the methodology. The lack of user-friendly computer programs has been a main obstacle for wider acceptance of the methods. This work provides two new adaptive MCMC methods: DRAM and AARJ. The DRAM method has been built especially to work in high dimensional and non-linear problems. The AARJ method is an extension to DRAM for model selection problems, where the mathematical formulation of the model is uncertain and we want simultaneously to fit several different models to the same observations. The methods were developed while keeping in mind the needs of modelling applications typical in environmental sciences. The development work has been pursued while working with several application projects. The applications presented in this work are: a winter time oxygen concentration model for Lake Tuusulanjärvi and adaptive control of the aerator; a nutrition model for Lake Pyhäjärvi and lake management planning; validation of the algorithms of the GOMOS ozone remote sensing instrument on board the Envisat satellite of European Space Agency and the study of the effects of aerosol model selection on the GOMOS algorithm.

A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Preparation of Polysulfone Membranes

As it is known, a huge part of all commercially available membranes are prepared by immersion precipitation. This way is the primary way to get flat membranes. The advantages of immersion precipitation are: wide field of the polymers, which can be used (polymer must be soluble in a solvent or a solvent mixture) and ease of performing. The literature part of this work deals with phase inversion membrane preparation methods and casting parameters affecting membrane performance. Also some membrane types and materials are discussed. In the experimental part of this work 73 membrane samples were made with different casting parameters (polymer concentration in the casting solution and precipitation time) and tested for the retention and permeability. The results of these experiments are collected and combined into the figures and tables which are presented in this thesis. This work showed and confirmed connection between membrane performance and casting parameters (concentration of polymer in the casting solution and precipitation time).

Multilevel converter modulation: implementation and analysis

Multilevel converters provide an attractive solution to bring the benefits of speed-controlled rotational movement to high-power applications. Therefore, multilevel inverters have attracted wide interest in both the academic community and in the industry for the past two decades. In this doctoral thesis, modulation methods suitable especially for series connected H-bridge multilevel inverters are discussed. A concept of duty cycle modulation is presented and its modification is proposed. These methods are compared with other well-known modulation schemes, such as space-vector pulse width modulation and carrier-based modulation schemes. The advantage of the modified duty-cycle modulation is its algorithmic simplicity. A similar mathematical formulation for the original duty cycle modulation is proposed. The modified duty cycle modulation is shown to produce well-formed phase-to-neutral voltages that have lower total harmonic distortion than the space-vector pulse width modulation and the duty cycle modulation. The space-vector-based solution and the duty cycle modulation, on the other hand, result in a better-quality line-to-line voltage and current waveform. The voltage of the DC links in the modules of the series-connected H-bridge inverter are shown to fluctuate while they are under load. The fluctuation causes inaccuracies in the voltage production, which may result in a failure of the flux estimator in the controller. An extension for upper-level modulation schemes, which changes the switching instants of the inverter so that the output voltage meets the reference voltage accurately regardless of the DC link voltages, is proposed. The method is shown to reduce the error to a very low level when a sufficient switching frequency is used. An appropriate way to organize the switching instants of the multilevel inverter is to make only one-level steps at a time. This causes restrictions on the dynamical features of the modulation schemes. The produced voltage vector cannot be rotated several tens of degrees in a single switching period without violating the above-mentioned one-level-step rule. The dynamical capabilities of multilevel inverters are analyzed in this doctoral thesis, and it is shown that the multilevel inverters are capable of operating even in dynamically demanding metal industry applications. In addition to the discussion on modulation schemes, an overvoltage in multilevel converter drives caused by cable reflection is addressed. The voltage reflection phenomenon in drives with long feeder cables causes premature insulation deterioration and also affects the commonmode voltage, which is one of the main reasons for bearing currents. Bearing currents, on the other hand, cause fluting in the bearings, which results in premature bearing failure. The reflection phenomenon is traditionally prevented by filtering, but in this thesis, a modulationbased filterless method to mitigate the overvoltage in multilevel drives is proposed. Moreover, the mitigation method can be implemented as an extension for upper-level modulation schemes. The method exploits the oscillations caused by two consecutive voltage edges so that the sum of the oscillations results in a mitigated peak of the overvoltage. The applicability of the method is verified by simulations together with experiments with a full-scale prototype.

Impedanssisovituksen tarkkailu jännitevälipiirillisessä taajuusmuuttajakäytössä

Invertterikäyttöjen pitkissä kuormakaapeleissa jännitteen aaltomuodossa on havaittavissa ylijännitteitä. Ylijännitteet johtuvat invertterikäytön käyttämien kytkinten ohjaamisesta pulssinleveysmodulaatiolla, niiden nopeista kytkeytymisistä ja itse järjestelmästä. Ylijännitteet ovat haitallisia järjestelmässä, joten niitä pyritään estämään esim. suodattamalla. Yksi käytännön sovellus sähkökoneen suotimena on RC-suodin. RC-suotimen toimintakyvyn aleneminen on kuitenkin oletettavaa, joten sähkökone altistuu tällöin ylijännitteille. Mittaamalla RC-suotimen kuntoa voidaan valvoa generaattorin kokemia ylijännitteitä. Työssä tarkastellaan järjestelmän häiriöitä ja ylijännitteitä, RC-suodinta sekä suotimen vajaatoimintatiloja, ja kehitetään mittaus tarkkailemaan vajaatoimintaa.

Vuotasapainon hallinta kokosiltasyötetyssä muuntajassa

Keskitaajuudella toimivia muuntajia käytetään laajalti tehoelektroniikkasovelluksissa kuten DC/DC-konverttereissa ja muissa hakkuriteholähteissä. Muuntaja on induktiivinen komponentti, jonka magneettisen tasapainon säilyttäminen hakkuriteholähteissä on laitteen virheettömän toiminnan kannalta tärkeää. Muuntajaa syöttävän virtapiirin on muodostettava symmetrinen syöttöjännite, jotta muuntajan vuo ei ajaudu positiiviseen tai negatiiviseen kyllästykseen. Tässä diplomityössä esitetään muuntajan sähkömagneettinen toimintaperiaate, kyllästymisen syyt hakkuriteholähteissä sekä kehitetään aktiivinen ohjaus vuotasapainon säilyttämiseksi. Hakkuriteholähteissä käytettävissä muuntajissa on monesti useampi kuin kaksi käämiä. Tässä työssä tutkittavassa muuntajassa on useita ensiöitä ja useita toisioita ja muuntajaa syötetään keskitaajuudella. Tämä tuo uusia ongelmia verrattuna perinteiseen yksivaiheiseen DC/DC-konvertteriin. Näihin ongelmiin esitetään ratkaisut diplomityön tutkimuksessa.

Galvaanisen erotuksen toteuttaminen tasavirtasähkönjakelun asiakaspäätelaitteessa

Suomessa sähkönjakeluverkon pääasialliset jännitetasot ovat 20 kV ja 400 V. 20 kV jännitetasolla sähkö viedään lähelle kuluttajia ja muunnetaan pienemmäksi lähellä asiakkaita. Haittapuolena on se, että 20 kV avojohtosähkönsiirtoverkko on hyvin vika-altis ja usein yhden haaran vikaantuessa monta muutakin jää ilman sähköä. Lisäksi hintavien ja suurien jakelumuuntajien määrä on suuri. Vaihtoehtona on toteuttaa osa sähkönjakelusta tasajännitteellä, jolloin tehollinen pienjännite olisi 400 V:a suurempi. Tällöin sähköä voitaisiin siirtää pidempiä matkoja ilman, että asiakaskohtaisia tai muutaman asiakkaan kattavia 20 kV siirtolinjoja tarvitsisi käyttää. Tämä edellyttää vaihtosuuntauksen toteuttamista kuluttajan päässä. Tässä työssä esiteltävällä 1 kVA:n tehoisella vaihtosuuntaajalla muodostetaan tasasähköjakeluverkosta saatavasta 750 V tasasähköstä yksivaiheista (230 VRMS, 50 Hz) verkkojännitettä. Laite on suunniteltu toteuttamaan galvaaninen erotus mahdollisimman hyvällä hyötysuhteella. Suurtaajuusmuuntajan mitoitus mahdollisimman hyvälle hyötysuhteelle on haastava tehtävä, koska vaatimuksia sille asettavat sekä syöttävä resonanssikonvertteri että syötettävä syklokonvertteri. Mitoituksessa on pyrittävä löytämään mahdollisimman hyvä hyötysuhde kustannusten ja toteutettavuuden suhteen.

Power cycling lifetime estimation of IGBT power modules based on chip temperature modeling

In this doctoral thesis, methods to estimate the expected power cycling life of power semiconductor modules based on chip temperature modeling are developed. Frequency converters operate under dynamic loads in most electric drives. The varying loads cause thermal expansion and contraction, which stresses the internal boundaries between the material layers in the power module. Eventually, the stress wears out the semiconductor modules. The wear-out cannot be detected by traditional temperature or current measurements inside the frequency converter. Therefore, it is important to develop a method to predict the end of the converter lifetime. The thesis concentrates on power-cycling-related failures of insulated gate bipolar transistors. Two types of power modules are discussed: a direct bonded copper (DBC) sandwich structure with and without a baseplate. Most common failure mechanisms are reviewed, and methods to improve the power cycling lifetime of the power modules are presented. Power cycling curves are determined for a module with a lead-free solder by accelerated power cycling tests. A lifetime model is selected and the parameters are updated based on the power cycling test results. According to the measurements, the factor of improvement in the power cycling lifetime of modern IGBT power modules is greater than 10 during the last decade. Also, it is noticed that a 10 C increase in the chip temperature cycle amplitude decreases the lifetime by 40%. A thermal model for the chip temperature estimation is developed. The model is based on power loss estimation of the chip from the output current of the frequency converter. The model is verified with a purpose-built test equipment, which allows simultaneous measurement and simulation of the chip temperature with an arbitrary load waveform. The measurement system is shown to be convenient for studying the thermal behavior of the chip. It is found that the thermal model has a 5 C accuracy in the temperature estimation. The temperature cycles that the power semiconductor chip has experienced are counted by the rainflow algorithm. The counted cycles are compared with the experimentally verified power cycling curves to estimate the life consumption based on the mission profile of the drive. The methods are validated by the lifetime estimation of a power module in a direct-driven wind turbine. The estimated lifetime of the IGBT power module in a direct-driven wind turbine is 15 000 years, if the turbine is located in south-eastern Finland.