Geostatistical prediction of clay percentage based on soil survey data

Selostus: Viljelymaiden savespitoisuuden alueellistaminen geostatistiikan ja pistemäisen tiedon avulla

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.

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.

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).

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.

Mikrohuokoisen pinnoitteen kehittäminen suodatinkankaan päälle

Mikrohuokoisia membraaneja käytetään monissa suodatussovelluksissa kuten vedenpuhdistuksessa. Tämän työn tavoite oli pinnoittaa suodatinkankaita mikrohuokoisella membraanilla. Membraanimateriaalin ja suodatuskankaan yhdistelmällä tavoiteltiin komposiittimateriaalia, jolla on membraanin selektiivinen erotuskyky ja kankaan mekaaniset ominaisuudet. Tämän työn kirjallisuusosa käsittelee membraanin valmistusta faasi-inversiomenetelmällä, pinnoitteeksi soveltuvia membraanimateriaaleja sekä pinnoitteen lujittamista. Kokeellisessa osassa suodatuskankaalle kehitettiin sopivaa pinnoitusmenetelmää. Pinnoitemateriaaleina käytettiin akryylidispersioainetta, polyamidia ja polysulfonia. Tuloksista huomattiin, ettei akryylidispersioaine soveltunut pinnoitemateriaaliksi. Kangasta onnistuttiin pinnoittamaan polyamidilla ja pinnoitusolosuhteita kontrolloimalla saatiin aikaiseksi pinnoite, joka pysyi kankaassa kiinni ja jonka rakenne koostui noin 12 μm huokosista. Polymeeriliuoksen pitoisuuden huomattiin vaikuttavan pinnoitteen rakenteeseen, mutta siihen jäi edelleen suuria avoimia huokosia. Suuret yksittäiset huokoset laskivat kuplapistepainetta ja tekivät huokoskokojakaumasta leveämmän. Polyamidipitoisuudella 14 ja 16 m–% pinnoite oli tunkeutunut kankaan sisään ja sen rakenne oli huokoinen, mutta pintakerrokseen muodostunut tiivis ihokerros tukki materiaalia. Tulosten perusteella kankaiden pinnoittaminen membraanikerroksella on mahdollista mutta pinnoitusprosessia täytyy edelleen kehittää, jotta päästään pienempään huokoskokoon ja kapeampaan huokoskokojakaumaan pinnoitteessa.

Conventional and nonconventional Kramers-Kronig analysis in optical spectroscopy

This Master’s Thesis is dedicated to the investigation and testing conventional and nonconventional Kramers-Kronig relations on simulated and experimentally measured spectra. It is done for both linear and nonlinear optical spectral data. Big part of attention is paid to the new method of obtaining complex refractive index from a transmittance spectrum without direct information of the sample thickness. The latter method is coupled with terahertz tome-domain spectroscopy and Kramers-Kronig analysis applied for testing the validity of complex refractive index. In this research precision of data inversion is evaluated by root-mean square error. Testing of methods is made over different spectral range and implementation of this methods in future is considered.

Membrane development for oil contaminated water treatments

The aim of this Master’s thesis study was to develop a membrane for oil contaminated water treatments. Oily wastewaters are a big problem to environment and therefore it is important to find an efficient method for their treatment. There are several treatment methods, but one of the most promising methods is membrane filtration. In the theoretical part of this study the membrane technology and polymeric membrane preparation with phase inversion and membrane modification methods was discussed. It was also told about the most important properties of the membranes. Oily waters, their treatment methods and oily wastewater sources were discussed more specifically. In the experimental part membranes from cellulose acetate were prepared and membranes were modified with two different methods. Modification methods were surface modification and polymer mixing. The modification purpose was to make membranes more hydrophilic and increase surface charge, which can reduce fouling. Membranes were characterized by determining zeta potential, contact angle, oil retention, pure water permeability, pressure-normalized flux and fouling. It were used both synthetic and real spent oil-water emulsion in membrane filtration. Surface modification resulted membranes, which had better properties than unmodified membrane. The amount of substance used in surface modification affected a lot to membrane properties, so it would be necessary to try different amounts of substance to develop the best membrane for oil-water emulsion treatment.

Spectral retinal image processing and analysis for ophthalmology

Diabetic retinopathy, age-related macular degeneration and glaucoma are the leading causes of blindness worldwide. Automatic methods for diagnosis exist, but their performance is limited by the quality of the data. Spectral retinal images provide a significantly better representation of the colour information than common grayscale or red-green-blue retinal imaging, having the potential to improve the performance of automatic diagnosis methods. This work studies the image processing techniques required for composing spectral retinal images with accurate reflection spectra, including wavelength channel image registration, spectral and spatial calibration, illumination correction, and the estimation of depth information from image disparities. The composition of a spectral retinal image database of patients with diabetic retinopathy is described. The database includes gold standards for a number of pathologies and retinal structures, marked by two expert ophthalmologists. The diagnostic applications of the reflectance spectra are studied using supervised classifiers for lesion detection. In addition, inversion of a model of light transport is used to estimate histological parameters from the reflectance spectra. Experimental results suggest that the methods for composing, calibrating and postprocessing spectral images presented in this work can be used to improve the quality of the spectral data. The experiments on the direct and indirect use of the data show the diagnostic potential of spectral retinal data over standard retinal images. The use of spectral data could improve automatic and semi-automated diagnostics for the screening of retinal diseases, for the quantitative detection of retinal changes for follow-up, clinically relevant end-points for clinical studies and development of new therapeutic modalities.