Toms-ilmiön vaikutuksen mallinnus turbulenssivirtauksessa

Diplomityön tavoitteena oli laatia putkivirtaukselle laskentamalli, joka huomioi turbulenssin vaimentamiseen käytetyn lisäaineen vaikutuksen. Kirjallisuusosassa käsitellään turbulenssiin ja putkiston painehäviön laskentaan liittyviä asioita. Lisäksi käydään läpi virtausvastusta alentavien aineiden (DRA-aineiden) ominaisuuksia, vaikutusmekanismia ja niiden vaikutusta virtausominaisuuksiin sekä laskentamalleihin. Kokeellisessa osassa laadittiin kaksi Excel-pohjaista laskentamallia, joilla huomioitiin lisäaineen vaikutus putkiston virtausmäärien parantumiseen ja painehäviöihin. Malli I laskee ensin virtausvastuksen alenemisen (DR-efektin) annetuista lähtötiedoista. DR-efektien perusteella saadaan laskettua putkiston painehäviötiedot. Malli II laskee ensin virtausvastuskertoimet ilman DRA-ainetta ja DRA-aineen kanssa. Tämän jälkeen malli II laskee DR-efektin. Mallien vaatimat lähtötiedot koostuvat virtausmääristä, fluidin ominaisuuksista, polymeerin ominaisuuksista, lisäaineen annostelutiedoista ja putkiston tiedoista.

Toimittajavarastointi teollisuusyrityksen kilpailuetuna sekä sen aiheuttamien kustannusten tarkastelu

Työn tavoitteena on tutkia toimittajavarastointia sekä siitä aiheutuvia kus-tannuksia. Työssä perehdytään tarkemmin toimitusketjun eri osapuolien väliseen yhteistyömalliin sekä tarkastellaan varastoinnin aiheuttamia kus-tannuksia yleisellä tasolla, sekä toimittajavarastointimallin mukanaan tuo-mien erillisten kustannuserien käsittelyä. Lisäksi tutkitaan, miten toimitta-javarastointimalli tuo lisähyötyä toimittajan toimintaan kilpailijoihin nähden. Tutkimus sisältää teoriaosuuden ja empiirisen osuuden, joissa esitellään aiheeseen liittyviä lainalaisuuksia ja erään teollisuusyrityksen toimittajavarastointia kahden esimerkkiasiakkaan kesken, sekä näiden asiakkaiden toimittajavarastoinnista esimerkkiyritykselle aiheutuvia kustannuksia. Teoriaosuudessa selvitetään toimittajavarastointiin liittyviä teoreettisia lähtökohtia ja sen vaikutuksia yritysten toimintaan. Kyseessä on kvalitatiivinen tapaustutkimus, jossa empiirinen aineisto ke-rättiin havainnoinnin, haastattelujen ja valmiiden dokumenttien perusteella sekä esimerkkiyrityksen tietokannoista. Teoriaosuuden lähteenä käytettiin aikaisempaa alan kirjallisuutta sekä artikkeleita. Tuloksena esitetään kustannusten laskentamalli, jonka avulla voidaan selvittää toimittajavarastoinnista aiheutuvat kustannukset ja niiden ositusperusteet.

Ballistisen simuloinnin liikeidean kehittämismahdollisuudet

Asevaikutusta laukaisusta kohteeseen simuloiva integroitu laskentaketju koostuu sisä-, ulko- ja maaliballistiikan malleista. Ulkoballistiikka kattaa laskentamallit radanlaskennan, sääkorjauksen ja ammusaerodynamiikan alueilla. Graafisella käyttöliittymällä toteutetulla, fysikaalisesti tarkkaan mallinnukseen perustuvalla ja kokonaisuuden kattavalla laskentajärjestelmällä on kasvavaa tarvetta teknisiä ja koulutuksellisia tarkoituksia varten. Erikoisesti, jos laskentaketjuun lisätään räjähdysvaikutuksen mallintaminen, voidaan simuloida asejärjestelmien vaikutusta kohteessa käyttäjien arvostamalla tavalla. Tietointensiiviset ballistiikan laskentamallit ovat välttämättömiä työkaluja teknisen suunnitteluosaamisen kattamiseksi ja kilpailuedun luomiseksi verkostoituneessa yritysympäristössä. Yliopistotutkimuksen tuottamien laskennallisten menetelmien hyötykäyttö yritysten suunnittelujärjestelmissä syventää teknistä osaamista, jolla on myös henkilöstöä motivoiva vaikutus teknisesti vaikeutuvilla markkinoilla. Työssä arvioidaan toimialaa analysoimalla eri käyttötarpeita samoille tietokantoihin tukeutuville laskentamalleille. Tarkastellaan teknisiä perusteita, käyttöympäristöjä ja markkinoita liiketoimintamahdollisuuksien tunnistamiseksi. Työn tuloksena syvennetään näkemystä ydinosaamisista ja visioidaan liikeidean erottumista kilpailijoista, markkinoita ja sen kehittämistä.

NH90 helikopterin moottoreiden syklilaskennan tehostaminen ja elinkaarihallinta

Tutkimustyö liittyy NH90 Helikopterin RTM-322 moottoreiden syklivaihdettavi-en osien huollonsuunnitteluun ja sen kehittämiseen. Tavoitteena oli löytää nykyi-sen käyttösyklilaskennan tilalle säästävämpi ratkaisumalli ja saada määritettyä suihkumoottorin osakokonaisuuksille lisälentoaikaa sekä huollonsuunnitteluun järkevämmät käyttöjaksot. Tähän päästiin selvittämällä käytännön lentotoiminnan avulla uudet laskentakertoimet ja -mallit moottorin eri osien eliniän määrittämi-seksi. Konkreettisena tavoitteena oli saada koneen käyttöhistoriaan perustuen kohdeorganisaatiolle esitys, miten moottoreiden käyttösyklilaskentaa tulisi muut-taa ja tehostaa. Diplomityön aihe muodostuu kolmesta tutkimusongelmakokonaisuudesta; moot-torin käyttöhistorian tutkimisesta, moottorivalmistajan hyväksymien syklilasken-tamenetelmien tutkimisesta ja moottorin syklien säästövaikutustarkastelusta tut-kimuksen määrittämillä uusilla laskentamenetelmillä. Tutkimuksen tuloksena muodostui NH90 huollonsuunnitteluun arvio mahdollisis-ta moottorisyklisäästöistä sekä moottoreiden näkökulmasta priorisoitu lista heli-kopterien moottoreiden laitevaihtohuoltojen järjestyksestä seuraavan viiden vuo-den aikana. Syklien säästötarkastelun pohjalta suoritettiin yhden koneen osalta koko käyttöhistorian ajalta takaisinlaskentatarkastelu. Tutkimuksessa selvitettiin Ilmavoimien meneillään olevien projektien valossa moottorisyklien tarkanlasken-nan toteuttamismahdollisuudet tulevaisuudessa. Tutkimuksessa saatujen tulosten avulla vältyttäisiin tulevaisuudessa RTM-322 moottorin ennenaikaisilta laitevaih-doilta koko NH90 kaluston suhteen ja näin saataisiin aikaan merkittävän suuria säästöjä vuoteen 2020 mennessä.

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.

Welding time models for cost calculations in the early stages of the design process

It is generally accepted that between 70 and 80% of manufacturing costs can be attributed to design. Nevertheless, it is difficult for the designer to estimate manufacturing costs accurately, especially when alternative constructions are compared at the conceptual design phase, because of the lack of cost information and appropriate tools. In general, previous reports concerning optimisation of a welded structure have used the mass of the product as the basis for the cost comparison. However, it can easily be shown using a simple example that the use of product mass as the sole manufacturing cost estimator is unsatisfactory. This study describes a method of formulating welding time models for cost calculation, and presents the results of the models for particular sections, based on typical costs in Finland. This was achieved by collecting information concerning welded products from different companies. The data included 71 different welded assemblies taken from the mechanical engineering and construction industries. The welded assemblies contained in total 1 589 welded parts, 4 257 separate welds, and a total welded length of 3 188 metres. The data were modelled for statistical calculations, and models of welding time were derived by using linear regression analysis. Themodels were tested by using appropriate statistical methods, and were found to be accurate. General welding time models have been developed, valid for welding in Finland, as well as specific, more accurate models for particular companies. The models are presented in such a form that they can be used easily by a designer, enabling the cost calculation to be automated.

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.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Long-term fertilizer field trials: comparison of three mathematical response models

Selostus: Lannoituksen pitkäaikaiset kenttäkokeet: kolmen matemaattisen mallin vertailu

Generalizing sample tree information with semiparametric and parametric models

Summary

Predicting the growth response to thinning for Scots pine stands using individual-tree growth models

Summary

Calculation and comparison of different permanent sample plot types

Summary

The application of agricultural land rating and crop models to CO2 and climate change issues in Northern regions : the Mackenzie Basin case study

Selostus: Viljelyvyöhykkeiden ja kasvumallien soveltaminen ilmastonmuutoksen tutkimisessa: Mackenzien jokialue, Kanada

Individual-tree basal area growth models for scots pine, pubescent birch and Norway spruce on drained peatlands in Finland

Summary