Kaivoslastauskoneen hybridisointianalyysi virtuaali- ja in-loop -simuloinnin avulla

Työssä tutkitaan raskaiden työkoneiden hybridisointimitoitusta simuloimalla. Työssä esitetään simulation-in-the-loop-simulointiin perustuva järjestelmä, jolla esimerkkitapauksena oleva kaivoslastauskone työympäristöineen voidaan mallintaa mekaaniselta osaltaan monikappaledynamiikkaan perustuvalla ohjelmistolla ja hybridijärjestelmän osalta Simulinkissa. Yhdistetty simulointi mahdollistaa hybridityökoneen virtuaalimallin ohjaamisen käyttäjän toimesta reaaliajassa. Simuloinnista saadaan tuloksena mm. työsykli, jota voidaan käyttää hybridisointimitoitukseen. Hybridisointi toteutetaan kahdella erilaisella kokoonpanolla, joista analysoidaan suorituskykyä sekä polttoaineen kulutusta. Tuloksia verrataan pelkästään dieselmoottoria voimanlähteenä käyttävään lastauskoneeseen. Työssä tehty tutkimus osoittaa, että (sarja-) hybridisoinnilla voidaan saavuttaa merkittäviä etuja raskaiden työkoneiden polttoainetehokkuudessa. Dieselmoottoria voidaan ajaa sellaisessa staattisessa toimintapisteessä, jonka hyötysuhde on korkea riippumatta työkoneen kuormituksesta. Saavutettu hyöty on toteutetussa tutkimuksessa parhaimmillaan jopa 56 % vähennys polttoaineenkulutuksessa. Lisäksi tarvittava dieselin nimellisteho pienenee huomattavasti. Tutkimuksen osana esitellään myös Hardware-in-the-Loop -laitteisto, jonka avulla voidaan liittää oikea sähkömoottori ja taajuudenmuuttaja osaksi virtuaalisesti simuloitua työkonetta.

Sähköajoneuvon sähkökäytön ja mekaniikan reaaliaikaisen yhdistelmäsimulaation toteuttaminen virtuaaliympäristössä

Esitetyn tutkimuksen tavoite on selvittää yksinkertaisen sähköajoneuvosimulaation avulla sähkökäyttösimulaation ja mekaniikkasimulaation välisen ohjelmallisen rajapinnan toimivuutta ja samalla toteuttaa LUT Energian hybridilaboratorion ensimmäinen rajapintaa hyödyntävä sähkökäytön ja mekaniikan yhdistelmäsimulaatio. Tutkimus on osa Lappeenrannan teknillisen yliopiston hybridikäyttötutkimusta, jossa muun muassa virtuaalisimuloinnin ja in-loop –simulaatioiden avulla tutkitaan raskaiden työkoneiden hybridisoinnin vaikutuksia esimerkiksi niiden suorituskykyyn, energiatehokkuuteen ja käytettävyyteen. Tulokset osoittavat mainitun rajapinnan toimivuuden, mikä mahdollistaa tutkimuksen etenemisen. Lisäksi raportti selvittää virtuaalisimuloinnin keskeisiä käsitteitä (virtual prototyping, hardware-in-the-loop sekä human-in-the-loop –simulointi).

Separation of Gas Bubbles from Circulation Waters and Pulp Suspensions in Open Channel Flow

The objective of this thesis was to study the removal of gases from paper mill circulation waters experimentally and to provide data for CFD modeling. Flow and bubble size measurements were carried out in a laboratory scale open gas separation channel. Particle Image Velocimetry (PIV) technique was used to measure the gas and liquid flow fields, while bubble size measurements were conducted using digital imaging technique with back light illumination. Samples of paper machine waters as well as a model solution were used for the experiments. The PIV results show that the gas bubbles near the feed position have the tendency to escape from the circulation channel at a faster rate than those bubbles which are further away from the feed position. This was due to an increased rate of bubble coalescence as a result of the relatively larger bubbles near the feed position. Moreover, a close similarity between the measured slip velocities of the paper mill waters and that of literature values was obtained. It was found that due to dilution of paper mill waters, the observed average bubble size was considerably large as compared to the average bubble sizes in real industrial pulp suspension and circulation waters. Among the studied solutions, the model solution has the highest average drag coefficient value due to its relatively high viscosity. The results were compared to a 2D steady sate CFD simulation model. A standard Euler-Euler k-ε turbulence model was used in the simulations. The channel free surface was modeled as a degassing boundary. From the drag models used in the simulations, the Grace drag model gave velocity fields closest to the experimental values. In general, the results obtained from experiments and CFD simulations are in good qualitative agreement.

Utilization of AMI-Technology in a Multilayered Demand Response Service Environment

This thesis studied the performance of Advanced metering infrastructure systems in a challenging Demand Response environment. The aim was to find out what kind of challenges and bottlenecks could be met when utilizing AMI-systems in challenging Demand Response tasks. To find out the challenges and bottlenecks, a multilayered demand response service concept was formed. The service consists of seven different market layers which consist of Nordic electricity market and the reserve markets of Fingrid. In the simulations the AMI-systems were benchmarked against these seven market layers. It was found out, that the current generation AMI-systems were capable of delivering Demand Response on the most challenging market layers, when observed from time critical viewpoint. Additionally, it was found out, that to enable wide scale Demand Response there are three major challenges to be acknowledged. The challenges hindering the utilization of wide scale Demand Response were related to poor standardization of the systems in use, possible problems in data connectivity solutions and the current electricity market regulation model.

Novel virtual environment and real-time simulation based methods for improving life-cycle efficiency of non-road mobile machinery

Virtual environments and real-time simulators (VERS) are becoming more and more important tools in research and development (R&D) process of non-road mobile machinery (NRMM). The virtual prototyping techniques enable faster and more cost-efficient development of machines compared to use of real life prototypes. High energy efficiency has become an important topic in the world of NRMM because of environmental and economic demands. The objective of this thesis is to develop VERS based methods for research and development of NRMM. A process using VERS for assessing effects of human operators on the life-cycle efficiency of NRMM was developed. Human in the loop simulations are ran using an underground mining loader to study the developed process. The simulations were ran in the virtual environment of the Laboratory of Intelligent Machines of Lappeenranta University of Technology. A physically adequate real-time simulation model of NRMM was shown to be reliable and cost effective in testing of hardware components by the means of hardware-in-the-loop (HIL) simulations. A control interface connecting integrated electro-hydraulic energy converter (IEHEC) with virtual simulation model of log crane was developed. IEHEC consists of a hydraulic pump-motor and an integrated electrical permanent magnet synchronous motorgenerator. The results show that state of the art real-time NRMM simulators are capable to solve factors related to energy consumption and productivity of the NRMM. A significant variation between the test drivers is found. The results show that VERS can be used for assessing human effects on the life-cycle efficiency of NRMM. HIL simulation responses compared to that achieved with conventional simulation method demonstrate the advances and drawbacks of various possible interfaces between the simulator and hardware part of the system under study. Novel ideas for arranging the interface are successfully tested and compared with the more traditional one. The proposed process for assessing the effects of operators on the life-cycle efficiency will be applied for wider group of operators in the future. Driving styles of the operators can be analysed statistically from sufficient large result data. The statistical analysis can find the most life-cycle efficient driving style for the specific environment and machinery. The proposed control interface for HIL simulation need to be further studied. The robustness and the adaptation of the interface in different situations must be verified. The future work will also include studying the suitability of the IEHEC for different working machines using the proposed HIL simulation method.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

Preparative ion exchange for salt solutions

Työssä tutkittiin moniarvoisten metalliformiaattien valmistusta ioninvaihto-menetelmällä. Kirjallisuustutkimus käsitteleetunnettuja alumiiniformiaatin ja rautaformiaatin valmistusmenetelmiä, kationinvaihtohartsien ominaisuuksia, ioninvaihtohartsien selektiivisyyttä ja alumiinin, raudan, magnesiumin ja sinkin vesikemiaa. Laboratoriokokeiden avulla tutkittiin sinkki-, magnesium-, rauta(II)- ja alumiiniformiaattien valmistusta ioninvaihdolla. Kokeet suoritettiin kolonnissa, joka oli pakattu makrohuokoisella tai geelimäisellä vahvalla kationin-vaihtohartsilla. Hartsi vaihdettiin natriummuodosta metallimuotoon metallikloridi- tai metallisulfaattiliuoksella.Metalli eluoitiin hartsista natriumformiaatilla. Formiaattien valmistus onnistui makrohuokoista vahvaa kationinvaihtohartsia käyttämällä. Rauta(II)formiaatin valmistus oli vaikeampaa kuin muiden formiaattien, koska rauta(II) hapettui osittain rauta(III):ksi valmistuksen aikana. Alumiiniformiaattia valmistettiin käyttäen sekä makrohuokoista että geelimäistä hartsia. Makrohuokoisen hartsin havaittiin soveltuvan geelimäistä hartsia paremmin alumiiniformiaatin valmistukseen. Kungeelimäistä hartsia käytettiin, noin 30 % alumiinista jäi kiinni hartsiin eikä siten eluoitunut. Ioninvaihdon selektiivisyyskertoimien saamiseksi suoritettiin tasapainokokeita. Selektiivisyyskertoimia käytettiin ioninvaihtokolonnin dynaamisessa simuloinnissa. Ioninvaihdon simuloiminen dynaamisella kolonnimallilla onnistui hyvin.

Lämpötilakerrostumien syntyminen ydinvoimalaitoksen ensiöpiirissä luonnonkierron olosuhteissa

Tässä diplomityössä on tutkittu lämpötilakerrostumien syntymistä RENATA-koelaitteistolla, joka muistutti geometrialtaan painevesireaktorin paineastian ylätilaa. Kokeet tehtiin siten, että aluksi RENATA täytettiin lämpimällä vedellä, jonka jälkeen koelaitteistoon juoksutettiin pohjasta käsin kylmää vettä. Kokeiden tuloksia verrattiin kirjallisuudessa esitettyyn korrelaatioon. Koetilanne mallinnettiin myös Fluent-virtauslaskentaohjelmalla, jolloin saatiin tietoa ohjelman kyvystä käsitellä lämpötilakerrostumia. Kokeiden tuloksissa havaittiin olevan selvää yhteyttä korrelaatioon. Korrelaation kriittistä rajaa suuremmilla arvoilla kylmä vesi kerrostui lämpimän veden alapuolelle. Lämpimän ja kylmän veden väliin muodostui muutaman senttimetrin paksuinen rajakerros, lämpötilakerrostuma, jossa lämpötilan muutos oli suurimmillaan parinkymmenen asteen luokkaa. Tämä lämpötilakerrostuma nousi hitaasti ylöspäin kokeen edetessä. Vastaavasti korrelaation kriittistä rajaa pienemmillä arvoilla lämmin ja kylmä vesi sekoittuivat keskenään. Myös Fluentilla lasketuissa simuloinneissa kylmä vesi kerrostui lämpimän veden alapuolelle. Lämpötilakerrostuma ei kuitenkaan noussut ylöspäin niin kuin kokeessa tapahtui, vaan se seisahtui koelaitteiston yläosaan.

Influence of mixed fuels combustion on boiler operation

The aim of this work is to study the effect of different fuel mixtures on the operation of circulating fluidized bed (CFB) boiler. The applicability of heat balance modeling software IPSEpro to simulate CFB boiler operation is also investigated. The work discusses various types of boilers and methods of boiler operation. The fuel properties and the possible fuel influence on the boiler efficiency are described. Various biofuel types that are possible to use in combination with other fuels are presented. Some examples of the fuel mixtures use are given. A CFB boiler model has been constructed using IPSEpro and applied to analyze boiler operation outside design conditions. In the simulations, the effect of different load levels and moisture contents for the fuel mixture has been studied.

Virtual technology and haptic interface solutions for design and control of mobile working machines

Commercially available haptic interfaces are usable for many purposes. However, as generic devices they are not the most suitable for the control of heavy duty mobile working machines like mining machines, container handling equipment and excavators. Alternative mechanical constructions for a haptic controller are presented and analysed. A virtual reality environment (VRE) was built to test the proposed haptic controller mechanisms. Verification of an electric motor emulating a hydraulic pump in the electro-hydraulic system of a mobile working machine is carried out. A real-time simulator using multi-body-dynamics based software with hardware-in-loop (HIL) setup was used for the tests. Recommendations for further development of a haptic controller and emulator electric motor are given.

Studies of Solar Activity with Emphasis on Quasi-Periodic Oscillations

In this thesis mainly long quasi-periodic solar oscillations in various solar atmospheric structures are discussed, based on data obtained at several wavelengths, focussing, however, mainly on radio frequencies. Sunspot (Articles II and III) and quiet Sun area (QSA) (Article I) oscillations are investigated along with quasi-periodic pulsations (QPP) in a flaring event with wide-range radio spectra (Article IV). Various oscillation periods are detected; 3–15, 35–70 and 90 minutes (QSA), 10-60 and 80-130 minutes (in sunspots at various radio frequencies), 3-5, 10-23, 220-240, 340 and 470 minutes (in sunspots at photosphere) and 8-12 and 15-17 seconds (in a solar flare at radio frequencies). Some of the oscillation periods are detected for the first time, while some of them have been confirmed earlier by other research groups. Solar oscillations can provide more information on the nature of various solar structures. This thesis presents the physical mechanisms of some solar structure oscillations. Two different theoretical approaches are chosen; magnetohydrodynamics (MHD) and the shallow sunspot model. These two theories can explain a wide range of solar oscillations from a few seconds up to some hours. Various wave modes in loop structures cause solar oscillations (<45 minutes) both in sunspots and quiet Sun areas. Periods lasting more than 45 minutes in the sunspots (and a fraction of the shorter periods) are related to sunspot oscillations as a whole. Sometimes similar oscillation periods are detected both in sunspot area variations and respectively in magnetic field strength changes. This result supports a concept that these oscillations are related to sunspot oscillations as a whole. In addition, a theory behind QPPs at radio frequencies in solar flares is presented. The thesis also covers solar instrumentation and data sources. Additionally, the data processing methods are presented. As the majority of the investigations in this thesis focus on radio frequencies, also the most typical radio emission mechanisms are presented. The main structures of the Sun, which are related to solar oscillations, are also presented. Two separate projects are included in this thesis. Solar cyclicity is studied using the extensively large solar radio map archieve from Metsähovi Radio Observatory (MRO) at 37 GHz, between 1978 and 2011 (Article V) covering two full solar cycles. Also, some new solar instrumentation (Article VI) was developed during this thesis.

Yliopiston sähköverkon mallintaminen

Lappeenrannan teknillinen yliopisto tutkii älykkäiden sähköverkkojen kehittämistä. Yliopisto on hankkinut sähköverkkoonsa tuuliturbiinin ja aurinkopaneeleita, joilla pystytään tuottamaan sähköenergiaa sähköverkkoon. Näitä tuotantoja voidaan käyttää myös tutkimuksessa. Tässä työssä luodaan simulaatiomalli yliopiston sähköverkosta Matlab® Simulink® -ohjelmalla. Simulaatiomalliin mallinnetaan yliopiston sisäinen keskijänniteverkko ja osa pienjänniteverkosta. Simulaatiomalli toteutetaan ohjelman valmiilla komponenteilla, joihin lasketaan tarvittavat parametrit. Tuuliturbiinin ja aurinkopaneelien sähköntuotantotehot määritetään säätiladatojen avulla. Verkon komponenteille lasketaan arvot komponenttien tyyppitietojen perusteella ja asetetaan simulaatiomallin parametreiksi. Simulaatiomalli luodaan yliopiston sisäisen verkon tehonjaon tarkastelemiseksi. Työssä selvitetään myös mahdollisuuksia luodun simulaatiomallin käyttämiseen vikatilanteiden tarkastelussa.

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.

Visualization of thermal-hydraulic nodes

The purpose of this thesis was to develop a program that can illustrate thermal-hydraulic node dimensions used in SMABRE simulations. These created node illustrations are used to verify the correctness of the designed simulation model and in addition they can be included in scientific reports. This thesis will include theory about SMABRE and relevant programs that were used to achieve the ending results. This thesis will give explanations for different modules that were created and used in the finished program, and it will present the different problems encountered and provide the solutions. The most important objective in this thesis is to display the results of generic VVER-1000 node dimensions and verify the correctness in the displayed part. The finished program was created using code language Python.

Probing MST1 Kinase Sequence and Structure for Rational Drug Discovery (Targeting diabetes by blocking protein-protein interactions)

Apoptotic beta cell death is an underlying cause majorly for type I and to a lesser extent for type II diabetes. Recently, MST1 kinase was identified as a key apoptotic agent in diabetic condition. In this study, I have examined MST1 and closely related kinases namely, MST2, MST3 and MST4, aiming to tackle diabetes by exploring ways to selectively block MST1 kinase activity. The first investigation was directed towards evaluating possibilities of selectively blocking the ATP binding site of MST1 kinase that is essential for the activity of the enzymes. Structure and sequence analyses of this site however revealed a near absolute conservation between the MSTs and very few changes with other kinases. The observed residue variations also displayed similar physicochemical properties making it hard for selective inhibition of the enzyme. Second, possibilities for allosteric inhibition of the enzyme were evaluated. Analysis of the recognized allosteric site also posed the same problem as the MSTs shared almost all of the same residues. The third analysis was made on the SARAH domain, which is required for the dimerization and activation of MST1 and MST2 kinases. MST3 and MST4 lack this domain, hence selectivity against these two kinases can be achieved. Other proteins with SARAH domains such as the RASSF proteins were also examined. Their interaction with the MST1 SARAH domain were evaluated to mimic their binding pattern and design a peptide inhibitor that interferes with MST1 SARAH dimerization. In molecular simulations the RASSF5 SARAH domain was shown to strongly interact with the MST1 SARAH domain and possibly preventing MST1 SARAH dimerization. Based on this, the peptidic inhibitor was suggested to be based on the sequence of RASSF5 SARAH domain. Since the MST2 kinase also interacts with RASSF5 SARAH domain, absolute selectivity might not be achieved.