Öljynjalostamon tyhjöyksikön esilämmitysuunin optimointi

Työn tarkoituksena oli vertailla kahta kaupallista prosessiuunisimulointiohjelmaa sekä selvittää onko toinen ohjelma soveltuva Neste Jacobs Oy:n käyttöön. Tavoitteena oli selvittää ohjelmien ominaisuuksien eroja ja rajoitteita sekä tuloksien yhdenmukaisuutta. Kirjallisuusosassa perehdyttiin prosessiuunien toimintaperiaatteisiin, yleisiin mitoitusperiaatteisiin ja säteilylämmönsiirron perusteisiin. Kirjallisuusosassa käsiteltiin myös prosessiuunien polttimien toimintaa ja NOx-päästöjä vähentäviä polttotekniikoita. Molemmilla laskentaohjelmilla laadittiin kuumaöljy- ja tyhjöuuni, joiden tuloksia vertailtiin tärkeimpien prosessiarvojen perusteella ja toisen ohjelman osalta myös laitevalmistajan ilmoittamiin arvoihin. Tällä ohjelmalla tutkittiin tyhjöuunin höyrynruiskutuksen ja koksaantumisen vaikutuksia prosessiin sekä säteilyosan leveyden muutoksen vaikutusta säteilylämmönsiirtoon. Ohjelmien ominaisuuksien keskeisimmät erot liittyivät kaksifaasivirtauksen laskentaan, prosessiaineen massavirran rajoituksiin, ilman esilämmittimen käyttöön ja laskentamallien paineen kontrollointiin. Tulokset erosivat savukaasun ja tuubien lämpötilojen ja painehäviön laskennassa. Tuloksien ja ominaisuuksien perusteella toisen ohjelman ei todettu olevan täysin soveltuva Neste Jacobs Oy:n käyttöön.

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Laskentamalli PKL-koelaitteiston pienen vuodon kokeen simuloimiseksi Apros-ohjelmistolla

Työn tavoitteena oli tehdä Apros-laskentamalli PKL-koelaitteistosta ja testata kuinka hyvin Apros pystyy laskemaan PKL-koelaitteistolla suoritetun E2.2 pienen vuodon kokeen. Tavoitteena oli myös tarkastella boorittoman veden tulpan etenemistä pienen vuodon kokeen aikana. PKL-koelaitteisto vastaa saksalaista sähköteholtaan 1300 MW olevaa Philippsburg 2 painevesilaitosta. Koelaitteiston tilavuudet ja teho on skaalattu kertoimella 145. Työssä tehdyllä laskentamallilla tarkasteltiin boorittoman veden tulpan liikkumista pienen vuodon kokeen aikana. Kun malli oli valmis, laskenta suoritettiin Apros 5.05 versiolla. Boorittoman veden tulpan etenemisen laskennassa käytettiin toisen kertaluvun diskretointia, jolla booripitoisuuden muutokset säilyvät teräväreunaisina. Laskentamalli pystyi kuvaamaan koelaitteistolla suoritetussa pienen vuodon kokeessa tapahtuneet ilmiöt varsin hyvin. Eroa koetuloksiintuli pääkiertopiirien luonnonkiertojen alkamishetkistä ja primääripaineen käyttäytymisessä. Kokeen alkutilanne ei ollut stationääritila, joten alkutilanteen asettamisessa oli hankaluuksia. Varsinkin pääkiertopiirien veden pinnankorkeuksienasettamisessa oli vaikeuksia, koska veden pinnankorkeuksien erot pyrkivät tasoittumaan nopeasti kokeen aikana. Apros pystyi laskemaan PKL-koelaitteistolla suoritetun pienen vuodon kokeen hyvin. Mallilla tulisi kuitenkin laskea vielä toisentyyppisiäkin kokeita, ennen kuin voidaan varmuudella tietää mallin toimivuus. PKL-koelaitteisto vastaa pääpiirteiltään Suomeen rakennettavaa Olkiluoto 3 ydinvoimalaitosta. Tehty työ antaa lisävarmuutta, kun Olkiluoto 3 laitoksen turvallisuustarkasteluita tehdään.

Modeling and Simulation ofReal-Time System as a Tool to Improve Capacity Management - Case Fenix

Tämän tutkimustyön kohteena on TietoEnator Oy:n kehittämän Fenix-tietojärjestelmän kapasiteettitarpeen ennustaminen. Työn tavoitteena on tutustua Fenix-järjestelmän eri osa-alueisiin, löytää tapa eritellä ja mallintaa eri osa-alueiden vaikutus järjestelmän kuormitukseen ja selvittää alustavasti mitkä parametrit vaikuttavat kyseisten osa-alueiden luomaan kuormitukseen. Osa tätä työtä on tutkia eri vaihtoehtoja simuloinnille ja selvittää eri vaihtoehtojen soveltuvuus monimutkaisten järjestelmien mallintamiseen. Kerätyn tiedon pohjaltaluodaan järjestelmäntietovaraston kuormitusta kuvaava simulaatiomalli. Hyödyntämällä mallista saatua tietoa ja tuotantojärjestelmästä mitattua tietoa mallia kehitetään vastaamaan yhä lähemmin todellisen järjestelmän toimintaa. Mallista tarkastellaan esimerkiksi simuloitua järjestelmäkuormaa ja jonojen käyttäytymistä. Tuotantojärjestelmästä mitataan eri kuormalähteiden käytösmuutoksia esimerkiksi käyttäjämäärän ja kellonajan suhteessa. Tämän työn tulosten on tarkoitus toimia pohjana myöhemmin tehtävälle jatkotutkimukselle, jossa osa-alueiden parametrisointia tarkennetaan lisää, mallin kykyä kuvata todellista järjestelmää tehostetaanja mallin laajuutta kasvatetaan.

Primääri-sekundäärivuodon tunnistaminen aktiivisuusmittauksin Loviisan voimalaitoksella

Diplomityössä tutkittiin Loviisan voimalaitoksen primääri- ja sekundääripiirin aktiivisuusmittausten kykyä tunnistaa pienet primääri-sekundäärivuodot. Tarkasteltavat primääri-sekundäärivuotojen suuruudet valittiin laitoksen hätätilanne- ja häiriönselvitysohjeiden mukaisesti. Vuodon vaikutuksia arvioitiin erilaisilla primäärijäähdytteen ominaisaktiivisuuksilla. Ominaisaktiivisuudet primääripiirissä määritettiin nuklidikohtaisesti erilaisille polttoainevuototapauksille. Työssä huomioitiin myös transienteissa mahdollisesti esiintyvä primääripiirin aktiivisuustasoa nostava spiking-ilmiö. Vuodon tarkempaa tunnistamista varten työssä laskettiin tarkasteltaville mittareille kalibrointikertoimet. Primääri-sekundäärivuoto mallinnettiin APROS-simulointiohjelmalla laitoksen eri käyttötiloissa ja kahdella eri vuotokoolla. Varsinainen aktiivisuuslaskenta suoritettiin SEKUN-ohjelmalla. Työssä tätä aktiivisuus- ja päästölaskentaohjelmaa muokattiin ohjelmoimalla siihen tarkasteltavat aktiivisuusmittaukset sekä primääripiirin puhdistus ja ulospuhallus. Laskelmien tuloksena saatiin arviot kunkin tarkasteltavana olleen aktiivisuusmittauksen soveltuvuudesta primääri-sekundäärivuodon tunnistamiseen erilaisissa polttoainevuototapauksissa ja reaktorin eri tehotasoilla. Häiriönselvitysohje I3:n käyttöönottoa edellyttävät vuotokoot määritettiin aktiivisuusmittausten havaitseman perusteella. Erityisesti kuumavalmiustilassa tapauksissa, joissa reaktorisydämessä oletetaan olevan tiiveytensä menettäneitä polttoainesauvoja, spikingin vaikutus jäähdytteiden aktiivisuuspitoisuuksiin ja mittaustuloksiin oli merkittävä. Niiltä osin, kuin tulokset käsittelevät ohjeissa vuodon tunnistamiseen käytettyjä aktiivisuusrajoja, tulokset osoittivat aktiivisuusrajat oikeiksi. Kuumavalmiudessa aktiivisuusmittausten mittausalueet saattavat joissakin tapauksissa rajoittaa primääri-sekundäärivuodon tunnistamista.

Development of a simulation tool for torsional vibration of branched power transmission system

In this master's thesis a mechanical model that is driven with variable speed synchronous machine was developed. The developed mechanical model simulates the mechanics of power transmission and its torsional vibrations. The mechanical model was developed for the need of the branched mechanics of a rolling mill and the propulsion system of a tanker. First, the scope of the thesis was to clarify the concepts connected to the mechanical model. The clarified concepts are the variable speed drive, the mechanics of power transmission and the vibrationsin the power transmission. Next, the mechanical model with straight shaft line and twelve moments of inertia that existed in the beginning was developed to be branched considering the case of parallel machines and the case of parallel rolls. Additionally, the model was expanded for the need of moreaccurate simulation to up to thirty moments of inertia. The model was also enhanced to enable three phase short circuit situation of the simulated machine. After that the mechanical model was validated by comparing the results of the developed simulation tool to results of other simulation tools. The compared results are the natural frequencies and mode shapes of torsional vibration, the response of the load torque step and the stress in the mechanical system occurred by the permutation of the magnetic field that is arisen from the three phase short circuit situation. The comparisons were accomplished well and the mechanical model was validated for the compared cases. Further development to be made is to develop the load torque to be time-dependent and to install two frequency converters and two FEM modeled machines to be simulated parallel.

Numerical Modelling of Combustion and Emission Formation in a Diesel Engine with Different Fuel Injection Characteristics

Päästöjen vähentäminen on ollut viime vuosina tärkeässä osassa polttomoottoreita kehitettäessä.Monet viralliset tahot asettavat uusia tiukempia päästörajoituksia. Päästörajatovat tyypillisesti olleet tiukimmat autoteollisuuden valmistamille pienille nopeakäyntisille diesel-moottoreille, mutta viime aikoina paineita on kohdistunut myös suurempiin keskinopeisiin ja hidaskäyntisiin diesel-moottoreihin. Päästörajat ovat erilaisia riippuen moottorin tyypistä, käytetystä polttoaineesta ja paikasta missä moottoria käytetään johtuen erilaisista paikallisista laeista ja asetuksista. Eniten huomiota diesel-moottorin päästöissä täytyy kohdistaa typen oksideihin, savun muodostukseen sekä partikkeleihin. Laskennallisen virtausmekaniikan (CFD) avulla on hyvät mahdollisuudet tutkia diesel-moottorin sylinterissä tapahtuvia ilmiöitä palamisen aikana. CFD on hyödyllinen työkalu arvioitaessa moottorin suorituskykyä ja päästöjen muodostumista. CFD:llä on mahdollista testata erilaisten parametrien ja geometrioiden vaikutusta ilman kalliita moottorinkoeajoja. CFD:tä voidaan käyttää myös opetustarkoituksessa lisäämään paloprosessin tuntemusta. Tulevaisuudessa palamissimuloinnit CFD:llä tulevat epäilemättä olemaan tärkeä osa moottorin kehityksessä. Tässä diplomityössä on tehty palamissimuloinnit kahteen erilaisilla poittoaineenruiskutuslaitteistoilla varustettuun Wärtsilän keskinopeaan diesel-moottoriin. W46 moottorin ruiskutuslaitteisto on perinteinen mekaanisesti ohjattu pumppusuutin ja W46-CR moottorissa on elektronisesti ohjattu 'common rail' ruiskutuslaitteisto. Näiden moottorien ja käytössä olevien ruiskutusprofiilien lisäksi on simuloinneilla testattu erilaisia uusia ruiskutusprofiileja, jotta erityyppisten profiilien hyvät ja huonot ominaisuudet tulisivat selville. Matalalla kuormalla kiinnostuksen kohteena on nokipäästöjen muodostus ja täydellä kuormalla NOx-päästöjen muodostus ja polttoaineen kulutus. Simulointien tulokset osoittivat, että noen muodostusta matalalla kuormalla voidaan selvästi vähentää monivaiheisella ruiskutuksella, jossa yksi ruiskutusjakso jaetaan kahteen tai useampaan jaksoon. Erityisen tehokas noen vähentämisessä vaikuttaa olevan ns. jälkiruiskutus (post injection). Matalat NOx-päästöt ja hyvä polttoaineen kulutus täydellä kuormalla on mahdollista saavuttaaasteittain nostettavalla ruiskutusnopeudella.

Improving distributed simulation in a workstation environment

Technological development brings more and more complex systems to the consumer markets. The time required for bringing a new product to market is crucial for the competitive edge of a company. Simulation is used as a tool to model these products and their operation before actual live systems are built. The complexity of these systems can easily require large amounts of memory and computing power. Distributed simulation can be used to meet these demands. Distributed simulation has its problems. Diworse, a distributed simulation environment, was used in this study to analyze the different factors that affect the time required for the simulation of a system. Examples of these factors are the simulation algorithm, communication protocols, partitioning of the problem, distributionof the problem, capabilities of the computing and communications equipment and the external load. Offices offer vast amounts of unused capabilities in the formof idle workstations. The use of this computing power for distributed simulation requires the simulation to adapt to a changing load situation. This requires all or part of the simulation work to be removed from a workstation when the owner wishes to use the workstation again. If load balancing is not performed, the simulation suffers from the workstation's reduced performance, which also hampers the owner's work. Operation of load balancing in Diworse is studied and it is shown to perform better than no load balancing, as well as which different approaches for load balancing are discussed.

Soodakattilan vesi-höyrykierron simulointi

Työssä tutkitaan soodakattilan vesi-höyrykierron laskentaa APROS-simulointiohjelman avulla sekä mahdollisuutta rakentaa soodakattilalle oma mallikirjasto. Lisäksi tarkasteltiin mahdollisuuksia räätälöidä APROS:ista yrityksen omaa tarvetta vastaava simulointiohjelma. Kirjallisuusosuudessa esitellään soodakattilan rakenne pääpiirteittäin sekä APROS-simulointiohjelman ja Nowa-kiertolaskuohjelman laskentaperusteet. Tässä osuudessa kerrotaan myös yleisesti ohjelmien sisällöstä, rakenteesta sekä niiden käytöstä. Työn kokeellisessa osassa on kerrottu soodakattilan vesi-höyrykierron mallinnuksesta APROS:in avulla sekä mallikirjaston luomisen perusteet. Lisäksi tässä osuudessa kerrotaan alimallien tekemisestä ja niiden tarpeellisuudesta sekä simuloinnin ongelmakohdista. Työstä saatujen tulosten perusteella dynaamisen simuloinnin ja mallikirjaston avulla saadaan kiertolaskuja laskettua helposti ja luotettavasti. Alimallikirjastoa ylläpitämällä nopeutetaan kiertolaskuvariaatioden tekoa huomattavasti.

Sahatavaran kuivauksen kehittäminen

Diplomityö tehtiin UPM-Kymmene puuteollisuuteen kuuluvalle Heinolan sahalle, joka valmistaa korkealaatuisia erikois- ja perussahatavaroita kotimaan ja ulkomaiden markkinoille. Sahatavarat kuivataan sahan lämminilmakuivaamoissa sahatavaroiden käyttötarkoitusta vaatimaan tasoon. Asiakkaiden vaatimusten kasvaessa on sahatavaran onnistunut kuivaus noussut erittäin tärkeäksi osaksi hyvää sahatavaratuotantoa. Diplomityön tavoitteena oli kehittää sahan kuivaustoimintaa. Kuivaustoimintaa kehitettiin ohjaamalla tavaravirtoja sopiviin kuivaamoihin kuivattavan sahatavaran ehdoilla, sekä kehittämällä uusia optimoituja kuivauskaavoja kamarikuivaamoille. Tämän lisäksi työssä selvitettiin sahan nykyinen kuivauskapasiteetti ja arvioitiin sen riittävyys nykyisiin ja tulevaisuuden tuotantomääriin. Lopuksi pohdittiin, miten ja millä keinoilla kuivauskapasiteettia voitaisiin jatkossa kehittää nykyistä suuremmaksi ilman mittavia investointeja. Kuivauskaavoja kehittämällä, sekä materiaalivirtoja ohjaamalla oikeisiin kuivaamoihin, pystyttiin kuivauskapasiteettia lisäämään selvästi aikaisemmasta tasosta. Kuivauskapasiteetin suuruus todettiin riittäväksi nykyisiin sahausmääriin verrattuna, mutta suunniteltuihin suurempiin sahausmääriin verrattuna, kapasiteetti ei ollut riittävän suuri. Kuitenkin, kehittämällä kuivaustoimintaa edelleen, voidaan kuivauskapasiteettia kasvattaa selvästi aikaisemmin lasketusta tasosta. Kuivauskaavojen optimoinnin apuvälineenä käytettiin LAATUKAMARI-simulointiohjelmaa, jonka todettiin soveltuvan hyvin kuivauskokeissa käytettyihin kamarikuivaamoihin. Simulointi on helpoin ja suositeltavin tapa kehittää kuivauskaavoja, varsinkin uusilla ja hyväkuntoisilla kamarikuivaamoilla. Samalla vältetään työläiden ja hitaiden kuivauskokeiden tekeminen. Hyvin tärkeää on käyttää simuloinnissa todellisia lähtöarvoja, oikean lopputuloksen aikaansaamiseksi.

Prosessihyötysuhteen parantamiskohteiden kartoitus painevesireaktorityyppisessä ydinvoimalaitoksessa

Työn tavoitteena on kartoittaa painevesireaktorityyppisen ydinvoimalaitoksen prosessihyötysuhteen parantamiskohteita. Aluksi kirjallisuudesta etsitään hyötysuhteen parantamiskeinoja ideaalisessa höyryvoimalaitosprosessissa. Näistä valitaan sopivimmat tarkastelun kohteeksi todellisessa voimalaitoksessa: syöttöveden esilämmityksen tehostaminen väliottohöyryvirtausta kasvattamalla ja syöttöveden esilämmittimen lämmönsiirtopintaa lisäämällä. Tarkastelussa pyritään löytämään paras mahdollinen hyötysuhde väliottohöyrylinjojen putkikokoa sekä esilämmittimien putkien lukumäärää muuttamalla. Diskreetin optimoinnin iteraatioaskel määritetään hyötysuhteen osittaisderivaattojen avulla. Tehtäviä muutoksia simuloidaan APROS-simulointiohjelmalla, jossa käytetään Loviisan voimalaitoksesta tehtyä mallia VVER-440. Työssä havaittiin, että pelkkiä väliottohöyrylinjojen putkikokoja – ja massavirtaa – kasvattamalla Loviisan voimalaitoksen hyötysuhdetta voidaan parantaa parhaimmillaan 32,75%:sta 32,85%:iin. Syöttöveden esilämmittimien lämmönsiirtopintaa lisäämällä saadaan suurempi parannus hyötysuhteeseen: 32,75%:sta 32,99%:iin. Näissä tapauksissa muutettiin kaikkia väliottohöyrylinjoja tai syöttöveden esilämmittimien lämpöpintoja. Työssä tarkasteltiin myös joitakin pienempiä muutoskohteita, joista paras hyötysuhteen kasvu saatiin korkeapaine-esilämmittimien lämmönsiirtopintaa kasvattamalla sekä toisen väliottohöyrylinjan (RD12) ja sitä vastaavan syöttöveden esilämmittimen muutosten yhteisvaikutuksena.

Sähkökäyttöisen konejärjestelmän vuorovaikutteinen simulointi

Työn tavoitteena oli toteuttaa simulointimalli, jolla pystytään tutkimaan kestomagnetoidun tahtikoneen aiheuttaman vääntömomenttivärähtelyn vaikutuksia sähkömoottoriin liitetyssä mekaniikassa. Tarkoitus oli lisäksi selvittää kuinka kyseinen simulointimalli voidaan toteuttaa nykyaikaisia simulointiohjelmia käyttäen. Saatujen simulointitulosten oikeellisuus varmistettiin tätä työtä varten rakennetulla verifiointilaitteistolla. Tutkittava rakenne koostui akselista, johon kiinnitettiin epäkeskotanko. Epäkeskotankoon kiinnitettiin massa, jonka sijaintia voitiin muunnella. Massan asemaa muuttamalla saatiin rakenteelle erilaisia ominaistaajuuksia. Epäkeskotanko mallinnettiin joustavana elementtimenetelmää apuna käyttäen. Mekaniikka mallinnettiin dynamiikan simulointiin tarkoitetussa ADAMS –ohjelmistossa, johon joustavana mallinnettu epäkeskotanko tuotiin ANSYS –elementtimenetelmäohjelmasta. Mekaniikan malli siirrettiin SIMULINK –ohjelmistoon, jossa mallinnettiin myös sähkökäyttö. SIMULINK –ohjelmassa mallinnettiin sähkökäyttö, joka kuvaa kestomagnetoitua tahtikonetta. Kestomagnetoidun tahtikoneen yhtälöt perustuvat lineaarisiin differentiaaliyhtälöihin, joihin hammasvääntömomentin vaikutus on lisätty häiriösignaalina. Sähkökäytön malli tuottaa vääntömomenttia, joka syötetään ADAMS –ohjelmistolla mallinnettuun mekaniikkaan. Mekaniikan mallista otetaan roottorin kulmakiihtyvyyden arvo takaisinkytkentänä sähkömoottorin malliin. Näin saadaan aikaiseksi yhdistetty simulointi, joka koostuu sähkötoimilaitekäytöstä ja mekaniikasta. Tulosten perusteella voidaan todeta, että sähkökäyttöjen ja mekaniikan yhdistetty simulointi on mahdollista toteuttaa valituilla menetelmillä. Simuloimalla saadut tulokset vastaavat hyvin mitattuja tuloksia.

CFD simulation of complex phenomena containing suspensions and flow throughporous media

There is an increasing reliance on computers to solve complex engineering problems. This is because computers, in addition to supporting the development and implementation of adequate and clear models, can especially minimize the financial support required. The ability of computers to perform complex calculations at high speed has enabled the creation of highly complex systems to model real-world phenomena. The complexity of the fluid dynamics problem makes it difficult or impossible to solve equations of an object in a flow exactly. Approximate solutions can be obtained by construction and measurement of prototypes placed in a flow, or by use of a numerical simulation. Since usage of prototypes can be prohibitively time-consuming and expensive, many have turned to simulations to provide insight during the engineering process. In this case the simulation setup and parameters can be altered much more easily than one could with a real-world experiment. The objective of this research work is to develop numerical models for different suspensions (fiber suspensions, blood flow through microvessels and branching geometries, and magnetic fluids), and also fluid flow through porous media. The models will have merit as a scientific tool and will also have practical application in industries. Most of the numerical simulations were done by the commercial software, Fluent, and user defined functions were added to apply a multiscale method and magnetic field. The results from simulation of fiber suspension can elucidate the physics behind the break up of a fiber floc, opening the possibility for developing a meaningful numerical model of the fiber flow. The simulation of blood movement from an arteriole through a venule via a capillary showed that the model based on VOF can successfully predict the deformation and flow of RBCs in an arteriole. Furthermore, the result corresponds to the experimental observation illustrates that the RBC is deformed during the movement. The concluding remarks presented, provide a correct methodology and a mathematical and numerical framework for the simulation of blood flows in branching. Analysis of ferrofluids simulations indicate that the magnetic Soret effect can be even higher than the conventional one and its strength depends on the strength of magnetic field, confirmed experimentally by Völker and Odenbach. It was also shown that when a magnetic field is perpendicular to the temperature gradient, there will be additional increase in the heat transfer compared to the cases where the magnetic field is parallel to the temperature gradient. In addition, the statistical evaluation (Taguchi technique) on magnetic fluids showed that the temperature and initial concentration of the magnetic phase exert the maximum and minimum contribution to the thermodiffusion, respectively. In the simulation of flow through porous media, dimensionless pressure drop was studied at different Reynolds numbers, based on pore permeability and interstitial fluid velocity. The obtained results agreed well with the correlation of Macdonald et al. (1979) for the range of actual flow Reynolds studied. Furthermore, calculated results for the dispersion coefficients in the cylinder geometry were found to be in agreement with those of Seymour and Callaghan.

Development a distributed simulation environment on a cluster of workstations

Simulation has traditionally been used for analyzing the behavior of complex real world problems. Even though only some features of the problems are considered, simulation time tends to become quite high even for common simulation problems. Parallel and distributed simulation is a viable technique for accelerating the simulations. The success of parallel simulation depends heavily on the combination of the simulation application, algorithm and message population in the simulation is sufficient, no additional delay is caused by this environment. In this thesis a conservative, parallel simulation algorithm is applied to the simulation of a cellular network application in a distributed workstation environment. This thesis presents a distributed simulation environment, Diworse, which is based on the use of networked workstations. The distributed environment is considered especially hard for conservative simulation algorithms due to the high cost of communication. In this thesis, however, the distributed environment is shown to be a viable alternative if the amount of communication is kept reasonable. Novel ideas of multiple message simulation and channel reduction enable efficient use of this environment for the simulation of a cellular network application. The distribution of the simulation is based on a modification of the well known Chandy-Misra deadlock avoidance algorithm with null messages. The basic Chandy Misra algorithm is modified by using the null message cancellation and multiple message simulation techniques. The modifications reduce the amount of null messages and the time required for their execution, thus reducing the simulation time required. The null message cancellation technique reduces the processing time of null messages as the arriving null message cancels other non processed null messages. The multiple message simulation forms groups of messages as it simulates several messages before it releases the new created messages. If the message population in the simulation is suffiecient, no additional delay is caused by this operation A new technique for considering the simulation application is also presented. The performance is improved by establishing a neighborhood for the simulation elements. The neighborhood concept is based on a channel reduction technique, where the properties of the application exclusively determine which connections are necessary when a certain accuracy for simulation results is required. Distributed simulation is also analyzed in order to find out the effect of the different elements in the implemented simulation environment. This analysis is performed by using critical path analysis. Critical path analysis allows determination of a lower bound for the simulation time. In this thesis critical times are computed for sequential and parallel traces. The analysis based on sequential traces reveals the parallel properties of the application whereas the analysis based on parallel traces reveals the properties of the environment and the distribution.

Enhancing the efficiency of freight transport by using simulation

The role of transport in the economy is twofold. As a sector of economic activity it contributes to a share of national income. On the other hand, improvements in transport infrastructure create room for accelerated economic growth. As a means to support railways as a safe and environmentally friendly transportation mode, the EU legislation has required the opening of domestic railway freight for competition from beginning of year 2007. The importance of railways as a mode of transport has been great in Finland, as a larger share of freight has been carried on rails than in Europe on average. In this thesis it is claimed that the efficiency of goods transport can be enhanced by service specific investments. Furthermore, it is stressed that simulation can and should be used to evaluate the cost-efficiency of transport systems on operational level, as well as to assess transportation infrastructure investments. In all the studied cases notable efficiency improvements were found. For example in distribution, home delivery of groceries can be almost twice as cost efficient as the current practice of visiting the store. The majority of the cases concentrated on railway freight. In timber transportation, the item with the largest annual transport volume in domestic railway freight in Finland, the transportation cost could be reduced most substantially. Also in international timber procurement, the utilization of railway wagons could be improved by combining complementary flows. The efficiency improvements also have positive environmental effects; a large part of road transit could be moved to rails annually. If impacts of freight transport are included in cost-benefit analysis of railway investments, up to 50 % increase in the net benefits of the evaluated alternatives can be experienced, avoiding a possible inbuilt bias in the assessment framework, and thus increasing the efficiency of national investments in transportation infrastructure. Transportation systems are a typical example of complex real world systems that cannot be analysed realistically by analytical methods, whereas simulation allows inclusion of dynamics and the level of detail required. Regarding simulation as a viable tool for assessing the efficiency of transportation systems finds support also in the international survey conducted for railway freight operators; operators use operations research methods widely for planning purposes, while simulation is applied only by the larger operators.