Fluid flow and dispersion predictions in porous media by utilizing advanced numerical methods

Diplomityön tavoitteena oli tarkastella numeerisen virtauslaskennan avulla virtaukseen liittyviä ilmiöitä ja kaasun dispersiota. Diplomityön sisältö on jaettu viiteen osaan; johdantoon, teoriaan, katsaukseen virtauksen mallinnukseen huokoisessa materiaalissa liittyviin tutkimusselvityksiin, numeeriseen mallinnukseen sekä tulosten esittämiseen ja johtopäätöksiin. Diplomityön alussa kiinnitettiin huomiota erilaisiin kokeellisiin, numeerisiin ja teoreettisiin mallinnusmenetelmiin, joilla voidaan mallintaa virtausta huokoisessa materiaalissa. Kirjallisuusosassa tehtiin katsaus aikaisemmin julkaistuihin puoliempiirisiin ja empiirisiin tutkimusselvityksiin, jotka liittyvät huokoisen materiaalin aiheuttamaan painehäviöön. Numeerisessa virtauslaskenta osassa rakennettiin ja esitettiin huokoista materiaalia kuvaavat numeeriset mallit käyttäen kaupallista FLUENT -ohjelmistoa. Työn lopussa arvioitiin teorian, numeerisen virtauslaskennan ja kokeellisten tutkimusselvitysten tuloksia. Kolmiulotteisen huokoisen materiaalinnumeerisessa mallinnuksesta saadut tulokset vaikuttivat lupaavilta. Näiden tulosten perusteella tehtiin suosituksia ajatellen tulevaa virtauksen mallinnusta huokoisessa materiaalissa. Osa tässä diplomityössä esitetyistä tuloksista tullaan esittämään 55. Kanadan Kemiantekniikan konferenssissa Torontossa 1619 Lokakuussa 2005. ASME :n kansainvälisessä tekniikan alan julkaisussa. Työ on hyväksytty esitettäväksi esitettäväksi laskennallisen virtausmekaniikan (CFD) aihealueessa 'Peruskäsitteet'. Lisäksi työn yksityiskohtaiset tulokset tullaan lähettämään myös CES:n julkaisuun.

ADVANCED NUMERICAL METHODS FOR SIMULATING MICROBUBBLE IN AN AERATED MIXING TANK

Fluid mixing in mechanically agitated tanks is one of the major unit operations in many industries. Bubbly flows have been of interest among researchers in physics, medicine, chemistry and technology over the centuries. The aim of this thesis is to use advanced numerical methods for simulating microbubble in an aerated mixing tank. Main components of the mixing tank are a cylindrical vessel, a rotating Rushton turbine and the air nozzle. The objective of Computational Fluid Dynamics (CFD) is to predict fluid flow, heat transfer, mass transfer and chemical reactions. The CFD simulations of a turbulent bubbly flow are carried out in a cylindrical mixing tank using large eddy simulation (LES) and volume of fluid (VOF) method. The Rushton turbine induced flow is modeled by using a sliding mesh method. Numerical results are used to describe the bubbly flows in highly complex liquid flow. Some of the experimental works related to turbulent bubbly flow in a mixing tank are briefly reported. Numerical simulations are needed to complete and interpret the results of the experimental work. Information given by numerical simulations has a major role in designing and scaling-up mixing tanks. The results of this work have been reported in the following scientific articles: ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Large eddy simulations and PIV experiments of a two-phase air-water mixer, in Proceedings of ASME Fluids Engineering Summer Conference (2005). ·Honkanen M., Koohestany A., Hatunen T., Saarenrinne P., Zamankhan P., Dynamical States of Bubbling in an Aerated Stirring Tank, submitted to J. Computational Physics.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

Methods and Models for Accelerating Dynamic Simulation of Fluid Power Circuits

The objective of this dissertation is to improve the dynamic simulation of fluid power circuits. A fluid power circuit is a typical way to implement power transmission in mobile working machines, e.g. cranes, excavators etc. Dynamic simulation is an essential tool in developing controllability and energy-efficient solutions for mobile machines. Efficient dynamic simulation is the basic requirement for the real-time simulation. In the real-time simulation of fluid power circuits there exist numerical problems due to the software and methods used for modelling and integration. A simulation model of a fluid power circuit is typically created using differential and algebraic equations. Efficient numerical methods are required since differential equations must be solved in real time. Unfortunately, simulation software packages offer only a limited selection of numerical solvers. Numerical problems cause noise to the results, which in many cases leads the simulation run to fail. Mathematically the fluid power circuit models are stiff systems of ordinary differential equations. Numerical solution of the stiff systems can be improved by two alternative approaches. The first is to develop numerical solvers suitable for solving stiff systems. The second is to decrease the model stiffness itself by introducing models and algorithms that either decrease the highest eigenvalues or neglect them by introducing steady-state solutions of the stiff parts of the models. The thesis proposes novel methods using the latter approach. The study aims to develop practical methods usable in dynamic simulation of fluid power circuits using explicit fixed-step integration algorithms. In this thesis, twomechanisms whichmake the systemstiff are studied. These are the pressure drop approaching zero in the turbulent orifice model and the volume approaching zero in the equation of pressure build-up. These are the critical areas to which alternative methods for modelling and numerical simulation are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. Furthermore, the second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by small fluid volumes, a pseudo-dynamic solver is proposed. Instead of integration of the pressure in a small volume, the pressure is solved as a steady-state pressure created in a separate cascade loop by numerical integration. The hydraulic capacitance V/Be of the parts of the circuit whose pressures are solved by the pseudo-dynamic method should be orders of magnitude smaller than that of those partswhose pressures are integrated. The key advantage of this novel method is that the numerical problems caused by the small volumes are completely avoided. Also, the method is freely applicable regardless of the integration routine applied. The superiority of both above-mentioned methods is that they are suited for use together with the semi-empirical modelling method which necessarily does not require any geometrical data of the valves and actuators to be modelled. In this modelling method, most of the needed component information can be taken from the manufacturer’s nominal graphs. This thesis introduces the methods and shows several numerical examples to demonstrate how the proposed methods improve the dynamic simulation of various hydraulic circuits.

Simulation methods in the modelling of bioaffinity assays

Computational model-based simulation methods were developed for the modelling of bioaffinity assays. Bioaffinity-based methods are widely used to quantify a biological substance in biological research, development and in routine clinical in vitro diagnostics. Bioaffinity assays are based on the high affinity and structural specificity between the binding biomolecules. The simulation methods developed are based on the mechanistic assay model, which relies on the chemical reaction kinetics and describes the forming of a bound component as a function of time from the initial binding interaction. The simulation methods were focused on studying the behaviour and the reliability of bioaffinity assay and the possibilities the modelling methods of binding reaction kinetics provide, such as predicting assay results even before the binding reaction has reached equilibrium. For example, a rapid quantitative result from a clinical bioaffinity assay sample can be very significant, e.g. even the smallest elevation of a heart muscle marker reveals a cardiac injury. The simulation methods were used to identify critical error factors in rapid bioaffinity assays. A new kinetic calibration method was developed to calibrate a measurement system by kinetic measurement data utilizing only one standard concentration. A nodebased method was developed to model multi-component binding reactions, which have been a challenge to traditional numerical methods. The node-method was also used to model protein adsorption as an example of nonspecific binding of biomolecules. These methods have been compared with the experimental data from practice and can be utilized in in vitro diagnostics, drug discovery and in medical imaging.

Numerical Simulation of Stochastic Di erential Equations

Stochastic differential equation (SDE) is a differential equation in which some of the terms and its solution are stochastic processes. SDEs play a central role in modeling physical systems like finance, Biology, Engineering, to mention some. In modeling process, the computation of the trajectories (sample paths) of solutions to SDEs is very important. However, the exact solution to a SDE is generally difficult to obtain due to non-differentiability character of realizations of the Brownian motion. There exist approximation methods of solutions of SDE. The solutions will be continuous stochastic processes that represent diffusive dynamics, a common modeling assumption for financial, Biology, physical, environmental systems. This Masters' thesis is an introduction and survey of numerical solution methods for stochastic differential equations. Standard numerical methods, local linearization methods and filtering methods are well described. We compute the root mean square errors for each method from which we propose a better numerical scheme. Stochastic differential equations can be formulated from a given ordinary differential equations. In this thesis, we describe two kind of formulations: parametric and non-parametric techniques. The formulation is based on epidemiological SEIR model. This methods have a tendency of increasing parameters in the constructed SDEs, hence, it requires more data. We compare the two techniques numerically.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Fluid flow in T-junction of pipes

The aim of this work is to study flow properties at T-junction of pipe, pressure loss suffered by the flow after passing through T-junction and to study reliability of the classical engineering formulas used to find head loss for T-junction of pipes. In this we have compared our results with CFD software packages with classical formula and made an attempt to determine accuracy of the classical formulas. In this work we have studies head loss in T-junction of pipes with various inlet velocities, head loss in T-junction of pipes when the angle of the junction is slightly different from 90 degrees and T-junction with different area of cross-section of the main pipe and branch pipe. In this work we have simulated the flow at T-junction of pipe with FLUENT and Comsol Multiphysics and observed flow properties inside the T-junction and studied the head loss suffered by fluid flow after passing through the junction. We have also compared pressure (head) losses obtained by classical formulas by A. Vazsonyi and Andrew Gardel and formulas obtained by assuming T-junction as combination of other pipe components and observations obtained from software experiments. One of the purposes of this study is also to study change in pressure loss with change in angle of T-junction. Using software we can have better view of flow inside the junction and study turbulence, kinetic energy, pressure loss etc. Such simulations save a lot of time and can be performed without actually doing the experiment. There were no real life experiments made, the results obtained completely rely on accuracy of software and numerical methods used.

On Topological Solitons in The Faddeev-Skyrme Model and its Extensions

The topological solitons of two classical field theories, the Faddeev-Skyrme model and the Ginzburg-Landau model are studied numerically and analytically in this work. The aim is to gain information on the existence and properties of these topological solitons, their structure and behaviour under relaxation. First, the conditions and mechanisms leading to the possibility of topological solitons are explored from the field theoretical point of view. This leads one to consider continuous deformations of the solutions of the equations of motion. The results of algebraic topology necessary for the systematic treatment of such deformations are reviewed and methods of determining the homotopy classes of topological solitons are presented. The Faddeev-Skyrme and Ginzburg-Landau models are presented, some earlier results reviewed and the numerical methods used in this work are described. The topological solitons of the Faddeev-Skyrme model, Hopfions, are found to follow the same mechanisms of relaxation in three different domains with three different topological classifications. For two of the domains, the necessary but unusual topological classification is presented. Finite size topological solitons are not found in the Ginzburg-Landau model and a scaling argument is used to suggest that there are indeed none unless a certain modification to the model, due to R. S. Ward, is made. In that case, the Hopfions of the Faddeev-Skyrme model are seen to be present for some parameter values. A boundary in the parameter space separating the region where the Hopfions exist and the area where they do not exist is found and the behaviour of the Hopfion energy on this boundary is studied.

Testate amoebae (thecamoebians) as indicators of aquatic mine impact

The environmental impacts of a single mine often remain local, but acidic and metal-rich acid mine drainage (AMD) from the waste materials may pose a serious threat to adjacent surface waters and their ecosystems. Testate amoebae (thecamoebian) analysis was used together with lake sediment geochemistry to study and evaluate the ecological effects of sulphidic metal mines on aquatic environments. Three different mines were included in the study: Luikonlahti Cu-mine in Kaavi, eastern Finland, Haveri Cu-Au mine in Ylöjärvi, southern Finland and Pyhäsalmi Zn-Cu-S mine in Pyhäjärvi, central Finland. Luikonlahti and Haveri are closed mines, but Pyhäsalmi is still operating. The sampling strategy was case specific, and planned to provide a representative sediment sample series to define natural background conditions, to detect spatial and temporal variations in mine impacts, to evaluate the possible recovery after the peak contamination, and to distinguish the effects of other environmental factors from the mining impacts. In the Haveri case, diatom analyses were performed alongside thecamoebian analysis to evaluate the similarities and differences between the two proxies. The results of the analyses were investigated with multivariate methods (direct and indirect ordinations, diversity and distance measure indices). Finally, the results of each case study were harmonized, pooled, and jointly analyzed to summarize the results for this dissertation. Geochemical results showed broadly similar temporal patterns in each case. Concentrations of ions in the pre-disturbance samples defined the natural baseline against which other results were compared. The beginning of the mining activities had only minor impacts on sediment geochemistry, mainly appearing as an increased clastic input into the lakes at Haveri and Pyhäsalmi. The active mining phase was followed by the metallic contamination and, subsequently, by the most recent change towards decreased but still elevated metal concentrations in the sediments. Because of the delay in the oxidation of waste material and formation of AMD, the most intense, but transient metal contamination phase occurred in the post-mining period at Luikonlahti and Haveri. At Pyhäsalmi, the highest metal contamination preceded effluent mitigation actions. Spatial gradients were observed besides the temporal evolution in both the pre-disturbance and mine-impacted samples from Luikonlahti and Pyhäsalmi. The geochemical gradients varied with distance from the main source of contaminants (dispersion and dilution) and with water depth (redox and pH). The spatial extent of the highest metal contamination associated with these mines remained rather limited. At Haveri, the metallic impact was widespread, with the upstream site in another lake basin found to be contaminated. Changes in thecamoebian assemblages corresponded well with the geochemical results. Despite some differences, the general features and ecological responses of the faunal assemblages were rather similar in each lake. Constantly abundant strains of Difflugia oblonga, Difflugia protaeiformis and centropyxids formed the core of these assemblages. Increasing proportions of Cucurbitella tricuspis towards the surface samples were found in all of the cases. The results affirmed the indicator value of some already known indicator forms, but such as C. tricuspis and higher nutrient levels, but also elicited possible new ones such as D. oblonga ‘spinosa’ and clayey substrate, high conductivity and/or alkalinity, D. protaeiformis ‘multicornis’ and pH, water hardness and the amount of clastic material and Centropyxis constricta ‘aerophila’ and high metal and S concentrations. In each case, eutrophication appeared to be the most important environmental factor, masking the effects of other variables. Faunal responses to high metal inputs in sediments remained minor, but were nevertheless detectable. Besides the trophic state of the lake, numerical methods suggested overall geochemical conditions (pH, redox) to be the most important factor at Luikonlahti, whereas the Haveri results showed the clearest connection between metals and amoebae. At Pyhäsalmi, the strongest relationships were found between Ca- and S-rich present loading, redox conditions and substrate composition. Sediment geochemistry and testate amoeba analysis proved to be a suitable combination of methods to detect and describe the aquatic mine impacts in each specific case, to evaluate recovery and to differentiate between the effects of different anthropogenic and natural environmental factors. It was also suggested that aquatic mine impacts can be significantly mitigated by careful design and after-care of the waste facilities, especially by reducing and preventing AMD. The case-specific approach is nevertheless necessary because of the unique characteristics of each mine and variations in the environmental background conditions.

Design and characterization of large area position sensitive radiation detectors

Planar, large area, position sensitive silicon detectors are widely utilized in high energy physics research and in medical, computed tomography (CT). This thesis describes author's research work relating to development of such detector components. The key motivation and objective for the research work has been the development of novel, position sensitive detectors improving the performance of the instruments they are intended for. Silicon strip detectors are the key components of barrel-shaped tracking instruments which are typically the innermost structures of high energy physics experimental stations. Particle colliders such as the former LEP collider or present LHC produce particle collisions and the silicon strip detector based trackers locate the trajectories of particles emanating from such collisions. Medical CT has become a regular part of everyday medical care in all developed countries. CT scanning enables x-ray imaging of all parts of the human body with an outstanding structural resolution and contrast. Brain, chest and abdomen slice images with a resolution of 0.5 mm are possible and latest CT machines are able to image whole human heart between heart beats. The two application areas are presented shortly and the radiation detection properties of planar silicon detectors are discussed. Fabrication methods and preamplifier electronics of the planar detectors are presented. Designs of the developed, large area silicon detectors are presented and measurement results of the key operating parameters are discussed. Static and dynamic performance of the developed silicon strip detectors are shown to be very satisfactory for experimental physics applications. Results relating to the developed, novel CT detector chips are found to be very promising for further development and all key performance goals are met.

Modelling of potentiometric ion sensors

Potentiometric ion sensors are a very important subgroup of electrochemical sensors, very attractive for practical applications due to their small size, portability, low-energy consumption, relatively low cost and not changing the sample composition. They are investigated by the researchers from many fields of science. The continuous development of this field creates the necessity for a detailed description of sensor response and the electrochemical processes important in the practical applications of ion sensors. The aim of this thesis is to present the existing models available for the description of potentiometric ion sensors as well as their applicability and limitations. This includes the description of the diffusion potential occurring at the reference electrodes. The wide range of existing models, from most idealised phase boundary models to most general models, including migration, is discussed. This work concentrates on the advanced modelling of ion sensors, namely the Nernst-Planck-Poisson (NPP) model, which is the most general of the presented models, therefore the most widely applicable. It allows the modelling of the transport processes occurring in ion sensors and generating the potentiometric response. Details of the solution of the NPP model (including the numerical methods used) are shown. The comparisons between NPP and the more idealized models are presented. The applicability of the model to describe the formation of diffusion potential in reference electrode, the lower detection limit of both ion-exchanger and neutral carrier electrodes and the effect of the complexation in the membrane are discussed. The model was applied for the description of both types of electrodes, i.e. with the inner filling solution and solidcontact electrodes. The NPP model allows the electrochemical methods other than potentiometry to be described. Application of this model in Electrochemical Impedance Spectroscopy is discussed and a possible use in chrono-potentiometry is indicated. By combining the NPP model with evolutionary algorithms, namely Hierarchical Genetic Strategy (HGS), a novel method allowing the facilitation of the design of ion sensors was created. It is described in detail in this thesis and its possible applications in the field of ion sensors are indicated. Finally, some interesting effects occurring in the ion sensors (i.e. overshot response and influence of anionic sites) as well as the possible applications of NPP in biochemistry are described.

FE-analyysin soveltaminen S960 QC teräksisen I-profiilin äärikestävyyden määrittämisessä

Tässä työssä tutkittiin FE-analyysin soveltamista S960 QC teräksisen I-profiilin kestävyyden määrittämisessä. Työn tavoitteena oli tarkastella nykyisten suunnitteluohjeiden soveltuvuutta ultralujille teräksille ja koota ohjemateriaali I-profiilin optimoimisesta sekä FE-analyysin hyö-dyntämisestä I-profiilin staattisen ja dynaamisen kestävyyden määrittämisessä. I-profiili mitoitettiin ja optimoitiin Eurokoodi 3:ssa esitettyjen PL3 mukaisten mitoitusohjeiden avulla. Rakenteelle suoritettiin Eurokoodi 3:n ja IIW:n mukaiset lommahdus-, kiepahdus- ja vä-symiskestävyystarkastelut. Väsymistarkastelussa sovellettiin nimellisen jännityksen, rakenteelli-sen jännityksen ja tehollisen lovijännityksen menetelmiä sekä murtumismekaniikkaa. Rakenteel-lisen jännityksen menetelmässä sovellettiin lisäksi lineaarista ja parabolista pintaa pitkin ekstra-polointia, paksuuden yli linearisointia sekä Dong:in menetelmää. Lommahdus-, kiepahdus- ja väsymistarkasteluissa hyödynnettiin analyyttistä laskentaa, FE-analyysiä sekä Frank2d sovellusta. Tarkastelujen perusteella voidaan todeta, että analyyttisillä menetelmillä saadaan numeerisia me-netelmiä varmemmalla puolella olevia tuloksia. Lommahdustarkastelussa ero tulosten välillä on suurimmillaan 8 % ja kiepahdustarkastelussa suurimmillaan 20 % mutta väsymistarkastelussa saadut tulokset eroavat keskenään huomattavasti. Väsymistarkastelussa tehollisen lovijännityksen menetelmällä sekä rakenteellisen jännityksen menetelmän Dong:in menetelmällä saadaan huo-mattavasti muita menetelmiä pidempiä kestoikiä, kun taas yksinkertaisemmilla menetelmillä saa-dut kestoiät ovat lyhyempiä. Rakenteen kestävyyden määrittäminen analyyttisillä menetelmillä on melko helppoa, mutta tu-lokset ovat monesti liian konservatiivisia. FE-analyysillä saadaan puolestaan hyvin tarkkoja tu-loksia mallin ollessa yksityiskohtainen. Mallintaminen on kuitenkin aikaa ja resursseja vievää ja vaatii käyttökokemusta. FE-analyysin mahdolliset hyödyt on aina arvioitava tapauskohtaisesti tarkasteltavan geometrian, kuormitusten ja reunaehtojen perusteella.