Sulasuolareaktorien historia

Kandidaatintyössä esitellään sulasuolareaktorien historian lisäksi lyhyesti reaktorityypin toimintaperiaate ja sulasuolareaktoritutkimuksen nykytilanne. Kandidaatintyöstä käy ilmi, että Oak Ridgen kansallinen laboratorio ORNL oli merkittävä tekijä sulasuolareaktoritutkimuksessa 1940-luvun lopulta aina 1980-luvun alkuun. Ainoat rakennetut sulasuolatekniikkaa käyttävät reaktorit olivat ORNL:n kehittämiä. Oak Ridgen kansallisen laboratorion lopetettua sulasuolatekniikan tutkimusohjelmansa, laantui sulasuolareaktoritutkimus pitkäksi aikaa. 2000-luvun alussa kiinnostus tekniikkaa kohtaan virkosi ja kirjoitushetkellä sulasuolareaktoreita kehitetään useissa tutkimusohjelmissa.

Novel technology for preparation of optically active chemicals

Asymmetric synthesis using modified heterogeneous catalysts has gained lots of interest in the production of optically pure chemicals, such as pharmaceuticals, nutraceuticals, fragrances and agrochemicals. Heterogeneous modified catalysts capable of inducing high enantioselectivities are preferred in industrial scale due to their superior separation and handling properties. The topic has been intensively investigated both in industry and academia. The enantioselective hydrogenation of ethyl benzoylformate (EBF) to (R)-ethyl mandelate over (-)-cinchonidine (CD)-modified Pt/Al2O3 catalyst in a laboratory-scale semi-batch reactor was studied as a function of modifier concentration, reaction temperature, stirring rate and catalyst particle size. The main product was always (R)-ethyl mandelate while small amounts of (S)-ethyl mandelate were obtained as by product. The kinetic results showed higher enantioselectivity and lower initial rates approaching asymptotically to a constant value as the amount of modifier was increased. Additionally, catalyst deactivation due to presence of impurities in the feed was prominent in some cases; therefore activated carbon was used as a cleaning agent of the raw material to remove impurities prior to catalyst addition. Detailed characterizations methods (SEM, EDX, TPR, BET, chemisorption, particle size distribution) of the catalysts were carried out. Solvent effects were also studied in the semi-batch reactor. Solvents with dielectric constant (e) between 2 and 25 were applied. The enantiomeric excess (ee) increased with an increase of the dielectric coefficient up to a maximum followed by a nonlinear decrease. A kinetic model was proposed for the enantioselectivity dependence on the dielectric constant based on the Kirkwood treatment. The non-linear dependence of ee on (e) successfully described the variation of ee in different solvents. Systematic kinetic experiments were carried out in the semi-batch reactor. Toluene was used as a solvent. Based on these results, a kinetic model based on the assumption of different number of sites was developed. Density functional theory calculations were applied to study the energetics of the EBF adsorption on pure Pt(1 1 1). The hydrogenation rate constants were determined along with the adsorption parameters by non-linear regression analysis. A comparison between the model and the experimental data revealed a very good correspondence. Transient experiments in a fixed-bed reactor were also carried out in this work. The results demonstrated that continuous enantioselective hydrogenation of EBF in hexane/2-propanol 90/10 (v/v) is possible and that continuous feeding of (-)-cinchonidine is needed to maintain a high steady-state enantioselectivity. The catalyst showed a good stability and high enantioselectivity was achieved in the fixed-bed reactor. Chromatographic separation of (R)- and (S)-ethyl mandelate originating from the continuous reactor was investigated. A commercial column filled with a chiral resin was chosen as a perspective preparative-scale adsorbent. Since the adsorption equilibrium isotherms were linear within the entire investigated range of concentrations, they were determined by pulse experiments for the isomers present in a post-reaction mixture. Breakthrough curves were measured and described successfully by the dispersive plug flow model with a linear driving force approximation. The focus of this research project was the development of a new integrated production concept of optically active chemicals by combining heterogeneous catalysis and chromatographic separation technology. The proposed work is fundamental research in advanced process technology aiming to improve efficiency and enable clean and environmentally benign production of enantiomeric pure chemicals.

Bio-ethanol valorization towards C4s including 1-butanol over metal modified alumina and zeolite catalysts

Bio-ethanol has been used as a fuel additive in modern society aimed at reducing CO2-emissions and dependence on oil. However, ethanol is unsuitable as fuel supplement in higher proportions due to its physico-chemical properties. One option to counteract the negative effects is to upgrade ethanol in a continuous fixed bed reactor to more valuable C4 products such as 1-butanol providing chemical similarity with traditional gasoline components. Bio-ethanol based valorization products also have other end-uses than just fuel additives. E.g. 1-butanol and ethyl acetate are well characterised industrial solvents and platform chemicals providing greener alternatives. The modern approach is to apply heterogeneous catalysts in the investigated reactions. The research was concentrated on aluminium oxide (Al2O3) and zeolites that were used as catalysts and catalyst supports. The metals supported (Cu, Ni, Co) gave very different product profiles and, thus, a profound view of different catalyst preparation methods and characterisation techniques was necessary. Additionally, acidity and basicity of the catalyst surface have an important role in determining the product profile. It was observed that ordinary determination of acid strength was not enough to explain all the phenomena e.g. the reaction mechanism. One of the main findings of the thesis is based on the catalytically active site which originates from crystallite structure. As a consequence, the overall evaluation of different by-products and intermediates was carried out by combining the information. Further kinetic analysis was carried out on metal (Cu, Ni, Co) supported self-prepared alumina catalysts. The thesis gives information for further catalyst developments aimed to scale-up towards industrially feasible operations.

Flexible multibody dynamics and intelligent control of a hydraulically driven hybrid redundant robot machine

The assembly and maintenance of the International Thermonuclear Experimental Reactor (ITER) vacuum vessel (VV) is highly challenging since the tasks performed by the robot involve welding, material handling, and machine cutting from inside the VV. The VV is made of stainless steel, which has poor machinability and tends to work harden very rapidly, and all the machining operations need to be carried out from inside of the ITER VV. A general industrial robot cannot be used due to its poor stiffness in the heavy duty machining process, and this will cause many problems, such as poor surface quality, tool damage, low accuracy. Therefore, one of the most suitable options should be a light weight mobile robot which is able to move around inside of the VV and perform different machining tasks by replacing different cutting tools. Reducing the mass of the robot manipulators offers many advantages: reduced material costs, reduced power consumption, the possibility of using smaller actuators, and a higher payload-to-robot weight ratio. Offsetting these advantages, the lighter weight robot is more flexible, which makes it more difficult to control. To achieve good machining surface quality, the tracking of the end effector must be accurate, and an accurate model for a more flexible robot must be constructed. This thesis studies the dynamics and control of a 10 degree-of-freedom (DOF) redundant hybrid robot (4-DOF serial mechanism and 6-DOF 6-UPS hexapod parallel mechanisms) hydraulically driven with flexible rods under the influence of machining forces. Firstly, the flexibility of the bodies is described using the floating frame of reference method (FFRF). A finite element model (FEM) provided the Craig-Bampton (CB) modes needed for the FFRF. A dynamic model of the system of six closed loop mechanisms was assembled using the constrained Lagrange equations and the Lagrange multiplier method. Subsequently, the reaction forces between the parallel and serial parts were used to study the dynamics of the serial robot. A PID control based on position predictions was implemented independently to control the hydraulic cylinders of the robot. Secondly, in machining, to achieve greater end effector trajectory tracking accuracy for surface quality, a robust control of the actuators for the flexible link has to be deduced. This thesis investigates the intelligent control of a hydraulically driven parallel robot part based on the dynamic model and two schemes of intelligent control for a hydraulically driven parallel mechanism based on the dynamic model: (1) a fuzzy-PID self-tuning controller composed of the conventional PID control and with fuzzy logic, and (2) adaptive neuro-fuzzy inference system-PID (ANFIS-PID) self-tuning of the gains of the PID controller, which are implemented independently to control each hydraulic cylinder of the parallel mechanism based on rod length predictions. The serial component of the hybrid robot can be analyzed using the equilibrium of reaction forces at the universal joint connections of the hexa-element. To achieve precise positional control of the end effector for maximum precision machining, the hydraulic cylinder should be controlled to hold the hexa-element. Thirdly, a finite element approach of multibody systems using the Special Euclidean group SE(3) framework is presented for a parallel mechanism with flexible piston rods under the influence of machining forces. The flexibility of the bodies is described using the nonlinear interpolation method with an exponential map. The equations of motion take the form of a differential algebraic equation on a Lie group, which is solved using a Lie group time integration scheme. The method relies on the local description of motions, so that it provides a singularity-free formulation, and no parameterization of the nodal variables needs to be introduced. The flexible slider constraint is formulated using a Lie group and used for modeling a flexible rod sliding inside a cylinder. The dynamic model of the system of six closed loop mechanisms was assembled using Hamilton’s principle and the Lagrange multiplier method. A linearized hydraulic control system based on rod length predictions was implemented independently to control the hydraulic cylinders. Consequently, the results of the simulations demonstrating the behavior of the robot machine are presented for each case study. In conclusion, this thesis studies the dynamic analysis of a special hybrid (serialparallel) robot for the above-mentioned special task involving the ITER and investigates different control algorithms that can significantly improve machining performance. These analyses and results provide valuable insight into the design and control of the parallel robot with flexible rods.

Particle model for simulating limestone reactions in novel fluidised bed energy applications

The development of carbon capture and storage (CCS) has raised interest towards novel fluidised bed (FB) energy applications. In these applications, limestone can be utilized for S02 and/or CO2 capture. The conditions in the new applications differ from the traditional atmospheric and pressurised circulating fluidised bed (CFB) combustion conditions in which the limestone is successfully used for SO2 capture. In this work, a detailed physical single particle model with a description of the mass and energy transfer inside the particle for limestone was developed. The novelty of this model was to take into account the simultaneous reactions, changing conditions, and the effect of advection. Especially, the capability to study the cyclic behaviour of limestone on both sides of the calcination-carbonation equilibrium curve is important in the novel conditions. The significances of including advection or assuming diffusion control were studied in calcination. Especially, the effect of advection in calcination reaction in the novel combustion atmosphere was shown. The model was tested against experimental data; sulphur capture was studied in a laboratory reactor in different fluidised bed conditions. Different Conversion levels and sulphation patterns were examined in different atmospheres for one limestone type. The Conversion curves were well predicted with the model, and the mechanisms leading to the Conversion patterns were explained with the model simulations. In this work, it was also evaluated whether the transient environment has an effect on the limestone behaviour compared to the averaged conditions and in which conditions the effect is the largest. The difference between the averaged and transient conditions was notable only in the conditions which were close to the calcination-carbonation equilibrium curve. The results of this study suggest that the development of a simplified particle model requires a proper understanding of physical and chemical processes taking place in the particle during the reactions. The results of the study will be required when analysing complex limestone reaction phenomena or when developing the description of limestone behaviour in comprehensive 3D process models. In order to transfer the experimental observations to furnace conditions, the relevant mechanisms that take place need to be understood before the important ones can be selected for 3D process model. This study revealed the sulphur capture behaviour under transient oxy-fuel conditions, which is important when the oxy-fuel CFB process and process model are developed.

Ydinvoimalaitosten vesikemia ja korroosionesto

Ydinvoimalaitosten vesikemian optimointi ja korroosionesto on välttämätöntä laitosten taloudellisen ja turvallisen käytön kannalta. Eri laitoksiin liittyvää vesikemiaa ja järjestelmissä havaittavia korroosion muotoja on tutkittu laajasti ja tutkitaan yhä edelleen. Monien prosessien ymmärtäminen vaatii usean eri tieteenalan osaamista, kuten kemiantekniikan, energiatekniikan sekä materiaalitekniikan. Tässä työssä kerrotaan yksinkertaistaen vesikemiaan ja korroosioon liittyviä prosesseja ja reaktioita. Työssä käsitellään kevytvettä jäähdytteenä sekä moderaattorina käyttävien ydinvoimalaitosten eri korroosiomuotoja sekä säteilyn vaikutusta näihin suoraan tai vesikemian kautta. Työssä kerrotaan korroosio- ja aktivoitumistuotteiden muodostumisesta ja kulkeutumisesta sekä näiden tuotteiden vaikutuksista laitosten toimintaan. Korroosion ja materiaalien aktivoitumisen pohjalta tarkastellaan kattavasti ydinvoimalaitosten tyypillisimpiä vesikemian muokkauskeinoja sekä korroosionhallintaa. Tärkeimpiin asioihin syvennytään hieman lähemmin. Tarkastelun kohteena ovat eniten käytetyt ydinvoimalaitokset, eli länsimaiset paine- ja kiehutusvesilaitokset sekä venäläisvalmisteiset VVER-laitokset. Tarkoituksena on ollut luoda tiivis tietopaketti opiskelijoiden käyttöön muun opintomateriaalin tueksi.

Decentralized management of organic household wastes in the Kathmandu Valley using small-scale composting reactors

The sustainable management of municipal solid waste in the Kathmandu Valley has always been a challenging task. Solid waste generation has gone rapidly high in the Kathmandu Valley over the last decade due to booming population and rapid urbaniza-tion. Finding appropriate landfill sites for the disposal of solid wastes generated from the households of the Kathmandu Valley has always been a major problem for Nepalese government. 65 % of total generated wastes from the households of Nepal consist of organic materials. As large fractions of generated household wastes are organic in na-ture, composting can be considered as one of the best sustainable ways to recycle organ-ic wastes generated from the households of Nepal. Model Community Society Development (MCDS), a non-governmental organization of Nepal carried out its small-scale project in five households of the Kathmandu Valley by installing composting reactors. This thesis is based on this small-scale project and has used secondary data provided by MCDS Nepal for carrying out the study. Proper man-agement of organic wastes can be done at household levels through the use of compost-ing reactors. The end product compost can be used as soil conditioners for agricultural purposes such as organic farming, roof-top farming and gardening. The overall average organic waste generation in the Kathmandu Valley is found to be 0,23 kg/person/day and the total amount of organic household wastes generated in the Kathmandu Valley is around 210 Gg/yr. Produced composts from five composting reac-tors contain high amount of moistures but have sufficient amount of nutrients required for the fertility of land and plant growth. Installation of five composting reactors in five households have prevented 2,74 Mg of organic wastes going into the landfills, thus re-ducing 107 kg of methane emissions which is equivalent to 2,7 Mg of carbondioxide.

Methyl and ethyl chloride synthesis in microreactors

Methyl chloride is an important chemical intermediate with a variety of applications. It is produced today in large units and shipped to the endusers. Most of the derived products are harmless, as silicones, butyl rubber and methyl cellulose. However, methyl chloride is highly toxic and flammable. On-site production in the required quantities is desirable to reduce the risks involved in transportation and storage. Ethyl chloride is a smaller-scale chemical intermediate that is mainly used in the production of cellulose derivatives. Thus, the combination of onsite production of methyl and ethyl chloride is attractive for the cellulose processing industry, e.g. current and future biorefineries. Both alkyl chlorides can be produced by hydrochlorination of the corresponding alcohol, ethanol or methanol. Microreactors are attractive for the on-site production as the reactions are very fast and involve toxic chemicals. In microreactors, the diffusion limitations can be suppressed and the process safety can be improved. The modular setup of microreactors is flexible to adjust the production capacity as needed. Although methyl and ethyl chloride are important chemical intermediates, the literature available on potential catalysts and reaction kinetics is limited. Thus the thesis includes an extensive catalyst screening and characterization, along with kinetic studies and engineering the hydrochlorination process in microreactors. A range of zeolite and alumina based catalysts, neat and impregnated with ZnCl2, were screened for the methanol hydrochlorination. The influence of zinc loading, support, zinc precursor and pH was investigated. The catalysts were characterized with FTIR, TEM, XPS, nitrogen physisorption, XRD and EDX to identify the relationship between the catalyst characteristics and the activity and selectivity in the methyl chloride synthesis. The acidic properties of the catalyst were strongly influenced upon the ZnCl2 modification. In both cases, alumina and zeolite supports, zinc reacted to a certain amount with specific surface sites, which resulted in a decrease of strong and medium Brønsted and Lewis acid sites and the formation of zinc-based weak Lewis acid sites. The latter are highly active and selective in methanol hydrochlorination. Along with the molecular zinc sites, bulk zinc species are present on the support material. Zinc modified zeolite catalysts exhibited the highest activity also at low temperatures (ca 200 °C), however, showing deactivation with time-onstream. Zn/H-ZSM-5 zeolite catalysts had a higher stability than ZnCl2 modified H-Beta and they could be regenerated by burning the coke in air at 400 °C. Neat alumina and zinc modified alumina catalysts were active and selective at 300 °C and higher temperatures. However, zeolite catalysts can be suitable for methyl chloride synthesis at lower temperatures, i.e. 200 °C. Neat γ-alumina was found to be the most stable catalyst when coated in a microreactor channel and it was thus used as the catalyst for systematic kinetic studies in the microreactor. A binder-free and reproducible catalyst coating technique was developed. The uniformity, thickness and stability of the coatings were extensively characterized by SEM, confocal microscopy and EDX analysis. A stable coating could be obtained by thermally pretreating the microreactor platelets and ball milling the alumina to obtain a small particle size. Slurry aging and slow drying improved the coating uniformity. Methyl chloride synthesis from methanol and hydrochloric acid was performed in an alumina-coated microreactor. Conversions from 4% to 83% were achieved in the investigated temperature range of 280-340 °C. This demonstrated that the reaction is fast enough to be successfully performed in a microreactor system. The performance of the microreactor was compared with a tubular fixed bed reactor. The results obtained with both reactors were comparable, but the microreactor allows a rapid catalytic screening with low consumption of chemicals. As a complete conversion of methanol could not be reached in a single microreactor, a second microreactor was coupled in series. A maximum conversion of 97.6 % and a selectivity of 98.8 % were reached at 340°C, which is close to the calculated values at a thermodynamic equilibrium. A kinetic model based on kinetic experiments and thermodynamic calculations was developed. The model was based on a Langmuir Hinshelwood-type mechanism and a plug flow model for the microreactor. The influence of the reactant adsorption on the catalyst surface was investigated by performing transient experiments and comparing different kinetic models. The obtained activation energy for methyl chloride was ca. two fold higher than the previously published, indicating diffusion limitations in the previous studies. A detailed modeling of the diffusion in the porous catalyst layer revealed that severe diffusion limitations occur starting from catalyst coating thicknesses of 50 μm. At a catalyst coating thickness of ca 15 μm as in the microreactor, the conditions of intrinsic kinetics prevail. Ethanol hydrochlorination was performed successfully in the microreactor system. The reaction temperature was 240-340°C. An almost complete conversion of ethanol was achieved at 340°C. The product distribution was broader than for methanol hydrochlorination. Ethylene, diethyl ether and acetaldehyde were detected as by-products, ethylene being the most dominant by-product. A kinetic model including a thorough thermodynamic analysis was developed and the influence of adsorbed HCl on the reaction rate of ethanol dehydration reactions was demonstrated. The separation of methyl chloride using condensers was investigated. The proposed microreactor-condenser concept enables the production of methyl chloride with a high purity of 99%.

Deoxygenation of fatty acids for production of fuels and chemicals

The decreasing fossil fuel resources combined with an increasing world energy demand has raised an interest in renewable energy sources. The alternatives can be solar, wind and geothermal energies, but only biomass can be a substitute for the carbon–based feedstock, which is suitable for the production of transportation fuels and chemicals. However, a high oxygen content of the biomass creates challenges for the future chemical industry, forcing the development of new processes which allow a complete or selective oxygen removal without any significant carbon loss. Therefore, understanding and optimization of biomass deoxygenation processes are crucial for the future bio–based chemical industry. In this work, deoxygenation of fatty acids and their derivatives was studied over Pd/C and TiO2 supported noble metal catalysts (Pt, Pt–Re, Re and Ru) to obtain future fuel components. The 5 % Pd/C catalyst was investigated in semibatch and fixed bed reactors at 300 °C and 1.7–2 MPa of inert and hydrogen–containing atmospheres. Based on extensive kinetic studies, plausible reaction mechanisms and pathways were proposed. The influence of the unsaturation in the deoxygenation of model compounds and industrial feedstock – tall oil fatty acids – over a Pd/C catalyst was demonstrated. The optimization of the reaction conditions suppressed the formation of by–products, hence high yields and selectivities towards linear hydrocarbons and catalyst stability were achieved. Experiments in a fixed bed reactor filled with a 2 % Pd/C catalyst were performed with stearic acid as a model compound at different hydrogen–containing gas atmospheres to understand the catalyst stability under various conditions. Moreover, prolonged experiments were carried out with concentrated model compounds to reveal the catalyst deactivation. New materials were proposed for the selective deoxygenation process at lower temperatures (~200 °C) with a tunable selectivity to hydrodeoxygenation by using 4 % Pt/TiO2 or decarboxylation/decarbonylation over 4 % Ru/TiO2 catalysts. A new method for selective hydrogenation of fatty acids to fatty alcohols was demonstrated with a 4 % Re/TiO2 catalyst. A reaction pathway and mechanism for TiO2 supported metal catalysts was proposed and an optimization of the process conditions led to an increase in the formation of the desired products.

Loistehon säätö ja kompensointi

Kantaverkkoyhtiö Fingrid Oy on tarkastellut kantaverkosta otettua tai tuotettua loistehoa alueittain Fingridin määrittelemällä maantieteellisellä tavalla. Alue koostuu liittymispisteistä, jotka ovat yksittäisten verkkoyhtiöiden liittymispisteitä. Yhden alueen sisällä voi siis olla monia eri verkkoyhtiöiden liittymispisteitä. Vuodesta 2016 alkaen loistehon tarkastelu vaihtuu aluetarkastelusta liittymispistetarkasteluun. Liittymispistetarkasteluun siirtymisen myötä tulee verkkoyhtiöiden investoida erinäisiin kompensointiratkaisuihin, jotta yksittäisten liittymispisteiden loistehot pysyvät Fingridin asettaman loistehon tarkasteluikkunan sisällä. Työssä perehdytään erinäisiin kompensointiratkaisuihin ja niiden kytkentöihin sekä kytkentöjen vaikutuksiin olemassa olevaan sähköverkkoon. Lisäksi työssä käsitellään tuulivoimaloiden hyödyntämistä loistehon säädössä. Esimerkkitapauksessa tarkastellaan erästä Caruna Oy:n verkkoa ja pohditaan teknistä ja mahdollisimman taloudellista kompensointiratkaisua kyseiseen verkkoon. Työn lopputulos on, että maakaapeloinnin määrän kasvaessa loistehon kompensoinnin tarve kasvaa. Näin ollen Fingrid Oy:ltä tarvitaan selkeitä linjauksia moneen ongelmakohtaan, kuten tuotanto-kulutuspisteen loistehon tarkasteluun. Lisäksi Energiavirastolta tarvitaan kannusteita loistehon kompensointiin, sillä verkkoyhtiöiden liiketoiminta perustuu regulaatiomalliin ja kompensointilaitteiston yksikköhinnasto on hyvin suppea.

Prosessiautomaatio pinnankorkeusmittauksen koelaitteistoissa

Pinnankorkeuden tunteminen kiehutusvesireaktorin painesäiliössä on erittäin tärkeää sen turvallisuusvaikutusten takia. Pinnankorkeutta mitataan vesipatsaiden korkeutta havaitsevien paine-eromittausten avulla. Säteilyturvakeskuksen YVL-ohjeiden mukaan turvallisuuteen vaikuttavien mittausten täytyy noudattaa moninkertaistus- ja erilaisuusperiaatteita. Yleensä erilaisuusperiaatetta on toteutettu käyttämällä erityyppisiä paine-eromittareita, mutta erilaisella fysikaalisella toimintaperiaatteella oleva mittaus olisi parempi ja toteuttaisi paremmin erilaisuusperiaatetta. Uimurikytkin olisi tällainen fysikaalisesti eri periaatteeseen perustuva pinnankorkeuden mittauslaite. Ydinvoimalaan tarkoitettu teknologia tulee kelpoistaa riippumattoman tahon toimesta ennen käyttöönottoa. Kelpoistamiskokeita varten Lappeenrannan teknillisen yliopiston Ydinturvallisuuden tutkimusyksikköön rakennettiin vuosina 2011–2013 kaksi koelaitteistoa. Näillä koelaitteistoilla tutkittiin uimurikytkimien toimintaa ja ominaisuuksia erilaisissa kiehutusvesireaktorin käyttötilanteissa. Koelaitteistot tarvitsivat toimiakseen automaatiojärjestelmät, jotka suunniteltiin pääosin noudattamalla suunnittelun elinkaarimallia sekä automaatiosuunnittelun sisältökokonaisuuksia. Automaatiojärjestelmien suunnittelu aloitettiin määrittelemällä koejärjestelyjen asettamat vaatimukset, jonka jälkeen tehtiin teknologiavalinnat. Seuraavaksi suunniteltiin automaatiojärjestelmien logiikkaohjelmistot, joiden kuvaukseen tämä työ pääasiassa keskittyy. Logiikkaohjelmistot toteutettiin graafisella National Instruments LabView -ohjelmointikielellä. Logiikkaohjelmistojen tuli hoitaa tiedonkeruuta, käyttöautomaatiota, turvallisuustehtäviä sekä kokeisiin liittyviä erikoistehtäviä. Ohjelmistot saatiin esikokeiden aikana toimimaan halutusti, ja varsinaiset kokeet voitiin suorittaa ilman merkittäviä ongelmia.

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.

Heterogeneous mass transfer in fluidized beds by computational fluid dynamics

The main objective of this research is to estimate and characterize heterogeneous mass transfer coefficients in bench- and pilot-scale fluidized bed processes by the means of computational fluid dynamics (CFD). A further objective is to benchmark the heterogeneous mass transfer coefficients predicted by fine-grid Eulerian CFD simulations against empirical data presented in the scientific literature. First, a fine-grid two-dimensional Eulerian CFD model with a solid and gas phase has been designed. The model is applied for transient two-dimensional simulations of char combustion in small-scale bubbling and turbulent fluidized beds. The same approach is used to simulate a novel fluidized bed energy conversion process developed for the carbon capture, chemical looping combustion operated with a gaseous fuel. In order to analyze the results of the CFD simulations, two one-dimensional fluidized bed models have been formulated. The single-phase and bubble-emulsion models were applied to derive the average gas-bed and interphase mass transfer coefficients, respectively. In the analysis, the effects of various fluidized bed operation parameters, such as fluidization, velocity, particle and bubble diameter, reactor size, and chemical kinetics, on the heterogeneous mass transfer coefficients in the lower fluidized bed are evaluated extensively. The analysis shows that the fine-grid Eulerian CFD model can predict the heterogeneous mass transfer coefficients quantitatively with acceptable accuracy. Qualitatively, the CFD-based research of fluidized bed process revealed several new scientific results, such as parametrical relationships. The huge variance of seven orders of magnitude within the bed Sherwood numbers presented in the literature could be explained by the change of controlling mechanisms in the overall heterogeneous mass transfer process with the varied process conditions. The research opens new process-specific insights into the reactive fluidized bed processes, such as a strong mass transfer control over heterogeneous reaction rate, a dominance of interphase mass transfer in the fine-particle fluidized beds and a strong chemical kinetic dependence of the average gas-bed mass transfer. The obtained mass transfer coefficients can be applied in fluidized bed models used for various engineering design, reactor scale-up and process research tasks, and they consequently provide an enhanced prediction accuracy of the performance of fluidized bed processes.

Simulation of the Interaction Mean Free Path between Neutrons and TRISO Particles

The interaction mean free path between neutrons and TRISO particles is simulated using scripts written in MATLAB to solve the increasing error present with an increase in the packing factor in the reactor physics code Serpent. Their movement is tracked both in an unbounded and in a bounded space. Their track is calculated, depending on the program, linearly directly using the position vectors of the neutrons and the surface equations of all the fuel particles; by dividing the space in multiple subspaces, each of which contain a fraction of the total number of particles, and choosing the particles from those subspaces through which the neutron passes through; or by choosing the particles that lie within an infinite cylinder formed on the movement axis of the neutron. The estimate from the current analytical model, based on an exponential distribution, for the mean free path, utilized by Serpent, is used as a reference result. The results from the implicit model in Serpent imply a too long mean free path with high packing factors. The received results support this observation by producing, with a packing factor of 17 %, approximately 2.46 % shorter mean free path compared to the reference model. This is supported by the packing factor experienced by the neutron, the simulation of which resulted in a 17.29 % packing factor. It was also observed that the neutrons leaving from the surfaces of the fuel particles, in contrast to those starting inside the moderator, do not follow the exponential distribution. The current model, as it is, is thus not valid in the determination of the free path lengths of the neutrons.

Thermal reactions of the major hydrocarbon components of biomass gasification gas

Gasification of biomass is an efficient method process to produce liquid fuels, heat and electricity. It is interesting especially for the Nordic countries, where raw material for the processes is readily available. The thermal reactions of light hydrocarbons are a major challenge for industrial applications. At elevated temperatures, light hydrocarbons react spontaneously to form higher molecular weight compounds. In this thesis, this phenomenon was studied by literature survey, experimental work and modeling effort. The literature survey revealed that the change in tar composition is likely caused by the kinetic entropy. The role of the surface material is deemed to be an important factor in the reactivity of the system. The experimental results were in accordance with previous publications on the subject. The novelty of the experimental work lies in the used time interval for measurements combined with an industrially relevant temperature interval. The aspects which are covered in the modeling include screening of possible numerical approaches, testing of optimization methods and kinetic modelling. No significant numerical issues were observed, so the used calculation routines are adequate for the task. Evolutionary algorithms gave a better performance combined with better fit than the conventional iterative methods such as Simplex and Levenberg-Marquardt methods. Three models were fitted on experimental data. The LLNL model was used as a reference model to which two other models were compared. A compact model which included all the observed species was developed. The parameter estimation performed on that model gave slightly impaired fit to experimental data than LLNL model, but the difference was barely significant. The third tested model concentrated on the decomposition of hydrocarbons and included a theoretical description of the formation of carbon layer on the reactor walls. The fit to experimental data was extremely good. Based on the simulation results and literature findings, it is likely that the surface coverage of carbonaceous deposits is a major factor in thermal reactions.