Improving Transportation and Warehousing Efficiency with Simulation-Based Decision Support Systems

Transportation and warehousing are large and growing sectors in the society, and their efficiency is of high importance. Transportation also has a large share of global carbondioxide emissions, which are one the leading causes of anthropogenic climate warming. Various countries have agreed to decrease their carbon emissions according to the Kyoto protocol. Transportation is the only sector where emissions have steadily increased since the 1990s, which highlights the importance of transportation efficiency. The efficiency of transportation and warehousing can be improved with the help of simulations, but models alone are not sufficient. This research concentrates on the use of simulations in decision support systems. Three main simulation approaches are used in logistics: discrete-event simulation, systems dynamics, and agent-based modeling. However, individual simulation approaches have weaknesses of their own. Hybridization (combining two or more approaches) can improve the quality of the models, as it allows using a different method to overcome the weakness of one method. It is important to choose the correct approach (or a combination of approaches) when modeling transportation and warehousing issues. If an inappropriate method is chosen (this can occur if the modeler is proficient in only one approach or the model specification is not conducted thoroughly), the simulation model will have an inaccurate structure, which in turn will lead to misleading results. This issue can further escalate, as the decision-maker may assume that the presented simulation model gives the most useful results available, even though the whole model can be based on a poorly chosen structure. In this research it is argued that simulation- based decision support systems need to take various issues into account to make a functioning decision support system. The actual simulation model can be constructed using any (or multiple) approach, it can be combined with different optimization modules, and there needs to be a proper interface between the model and the user. These issues are presented in a framework, which simulation modelers can use when creating decision support systems. In order for decision-makers to fully benefit from the simulations, the user interface needs to clearly separate the model and the user, but at the same time, the user needs to be able to run the appropriate runs in order to analyze the problems correctly. This study recommends that simulation modelers should start to transfer their tacit knowledge to explicit knowledge. This would greatly benefit the whole simulation community and improve the quality of simulation-based decision support systems as well. More studies should also be conducted by using hybrid models and integrating simulations with Graphical Information Systems.

Painevesityyppisen ydinvoimalaitoksen sekundääripiirin mallinnus

Tutkimuksen tavoitteena oli kehittää entistä tarkempi ja luotettavampi työkalu Loviisan ydinvoimalaitoksen ensimmäisen reaktoriyksikön sekundääripiirin 50-linjan stationääritilan toiminnan arviointiin. Toisena tavoitteena oli kartoittaa ja raportoida mallintamiseen käytetyn SOLVO-ohjelman kehitysmahdollisuuksia. Tutkimusten perusteella kehitettiin testatusti toimiva prosessimalli, joka dokumentoitiin ja validoitiin yksityiskohtaisesti. Tulevaisuudessa tapahtuvaa kehitystyötä, ja mallin käyttöä varten kehitettiin lisäksi Excel-pohjainen työkalu, jonka avulla SOLVO:n ja Excelin välinen integraatio voidaan myöhemmin julkaistavassa SOLVO:n versiossa viedä nykyistä pidemmälle. Mallintamistyön ohella tutkimuksen aikana selvitettiin olennaisimmat SOLVO:n erityispiirteet ja vahvuudet sekä kartoitettiin sen käyttöön liittyvät kehitystarpeet. Tärkeimpänä kehitysehdotuksena nousi esiin yksittäisissä komponenteissa suoritetun laskennan läpinäkyvyyden parantaminen. Seuraavassa kehitysvaiheessa myös komponenttikohtaiset laskentayhtälöt olisi suositeltavaa asettaa avoimiksi käyttäjäkohtaisille muutoksille. Työn aikana saavutettiin myös muita merkittäviä tuloksia, jotka liittyivät pääosin rinnak-kaisten 10- ja 50-linjojen välisiin yhteyksiin. Linjojen välisiä vaikutuksia analysoitaessa huomattiin niiden olevan olennaisessa asemassa erityisesti sarja-ajon aikana. Mikäli mallilla halutaan kuvata sekä sarja- että rinnanajoa, sen on käsitettävä molemmat linjat ja kaikki niihin liittyvät komponentit. Edellä mainitun lisäksi mallipohjaisen tarkastelun tuloksena tehtiin havaintoja, joiden perusteella nykyistä prosessia voidaan edelleen kehittää. Näistä havainnoista tärkeimpänä merivesipumppujen optimaalisen säätölämpötilan todettiin asettuvan 4,5 – 4,6 °C välille. Toinen huomio liittyi matalapaineturbiinien ulosvirtaushäviöihin, joihin hukataan juoksusta riippuen keskimäärin noin 10 kJ/kg enemmän entalpiaa kuin parhaassa mahdollisessa tapauksessa. Validoinnin yhteydessä havaituista pienistä poikkeamista huolimatta kehitetty malli vastaa hyvin laitokselta saatuja mittaustuloksia sekä muita samassa yhteydessä käytettyjä luotet-tavuuden arviointikriteerejä.

Flexible multibody approach in bone strain estimation during physical activity:

Bone strain plays a major role as the activation signal for the bone (re)modeling process, which is vital for keeping bones healthy. Maintaining high bone mineral density reduces the chances of fracture in the event of an accident. Numerous studies have shown that bones can be strengthened with physical exercise. Several hypotheses have asserted that a stronger osteogenic (bone producing) effect results from dynamic exercise than from static exercise. These previous studies are based on short-term empirical research, which provide the motivation for justifying the experimental results with a solid mathematical background. The computer simulation techniques utilized in this work allow for non-invasive bone strain estimation during physical activity at any bone site within the human skeleton. All models presented in the study are threedimensional and actuated by muscle models to replicate the real conditions accurately. The objective of this work is to determine and present loading-induced bone strain values resulting from physical activity. It includes a comparison of strain resulting from four different gym exercises (knee flexion, knee extension, leg press, and squat) and walking, with the results reported for walking and jogging obtained from in-vivo measurements described in the literature. The objective is realized primarily by carrying out flexible multibody dynamics computer simulations. The dissertation combines the knowledge of finite element analysis and multibody simulations with experimental data and information available from medical field literature. Measured subject-specific motion data was coupled with forward dynamics simulation to provide natural skeletal movement. Bone geometries were defined using a reverse engineering approach based on medical imaging techniques. Both computed tomography and magnetic resonance imaging were utilized to explore modeling differences. The predicted tibia bone strains during walking show good agreement with invivo studies found in the literature. Strain measurements were not available for gym exercises; therefore, the strain results could not be validated. However, the values seem reasonable when compared to available walking and running invivo strain measurements. The results can be used for exercise equipment design aimed at strengthening the bones as well as the muscles during workout. Clinical applications in post fracture recovery exercising programs could also be the target. In addition, the methodology introduced in this study, can be applied to investigate the effect of weightlessness on astronauts, who often suffer bone loss after long time spent in the outer space.

Testing and simulation of fractionary electromechanical rotative drives

This paper concerns the development of drives that use electromechanical rotative motor systems. It is proposed an experimental drive test structure integrated to simulation softwares. The objective of this work is to show that an affordable model validation procedure can be obtained by combining a precision data acquisition with well tuned state-of-the-art simulation packages. This is required for fitting, in the best way, a drive to its load or, inversely, to adapt loads to given drive characteristics.

Acceleration of Solar Energetic Particles in coronal shocks through self-generated turbulence

The acceleration of solar energetic particles (SEPs) by flares and coronal mass ejections (CMEs) has been a major topic of research for the solar-terrestrial physics and geophysics communities for decades. This thesis discusses theories describing first-order Fermi acceleration of SEPs through repeated crossings at a CME-driven shock. We propose that particle trapping occurs through self-generated Alfvén waves, leading to a turbulent trapping region in front of the shock. Decelerating coronal shocks are shown to be capable of efficient SEP acceleration, provided seed particle injection is sufficient. Quasi-parallel shocks are found to inject thermal particles with good efficiency. The roles of minimum injection velocities, cross-field diffusion, downstream scattering efficiency and cross-shock potential are investigated in detail, with downstream isotropisation timescales having a major effect on injection efficiency. Accelerated spectra of heavier elements up to iron are found to exhibit significantly harder spectra than protons. Accelerated spectra cut-off energies are found to scale proportional to (Q/A)^1.5, which is explained through analysis of the spectral shape of amplified Alfvénic turbulence. Acceleration times to different threshold energies are found to be non-linear, indicating that self-consistent time-dependent simulations are required in order to expose the full extent of acceleration dynamics. The well-established quasilinear theory (QLT) of particle scattering is investigated by comparing QLT scattering coefficients with those found via full-orbit simulations. QLT is found to overemphasise resonance conditions. This finding supports the simplifications implemented in the presented coronal shock acceleration (CSA) simulation software. The CSA software package is used to simulate a range of acceleration scenarios. The results are found to be in agreement with well-established particle acceleration theory. At the same time, new spatial and temporal dynamics of particle population trapping and wave evolution are revealed.

Numerical Simulation of Stochastic Di erential Equations

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

Häiriöanalyysejä tehonkorotetun kiehutusvesivoimalaitoksen käytettävyyden varmistamiseksi

TVO suunnittelee reaktoritehon 10 %:n korotusta Olkiluoto 1 ja 2 -voimalaitoksille. Reaktoriteho nostetaan 2500 MW:sta 2750 MW:iin polttoaineen rikastusastetta nostamalla ja pääkiertovirtausta kasvattamalla. Samalla syöttövesivirtaus reaktoriin ja tuorehöyryvirtaus turpiineille kasvaa. Lauhteenpuhdistusjärjestelmän kapasiteettia ei voida kuitenkaan kasvattaa, joten massavirran lisäys toteutetaan ottamalla käyttöön korkeapainesivulauhteen eteenpäinpumppaus. Lauhteen esilämmityslinjojen, lauhteenpuhdistuksen ja syöttövesipumppujen massavirta säilyy siten nykyisellään. Muita merkittäviä tehonkorotukseen liittyviä laitosmuutoksia ovat pääkiertopumppujen uusinta ja korkeapaineturpiinin muutokset. Tehonkorotetun prosessin käytettävyyden varmistamiseksi tehdään häiriöanalyysejä Apros-prosessisimulointiohjelmistoa käyttäen. OL1 ja OL2 -laitoksista on olemassa validoitu 2500 MW:n laitosmalli, josta muokatulla 2750 MW:n laitosmallilla simuloinnit tehdään. Häiriöanalyysien avulla selvitetään säätöjärjestelmien kyky pitää prosessin tila hallinnassa ilman suojausautomaation laukeamista. Simuloituihin tapauksiin kuuluu pumppujen ja venttiilien vikaantumistapauksia sekä turpiini- ja reaktoripuolen pikasulku- ja osittaispikasulkutapauksia. Myös meriveden lämpötilan vaikutusta häiriötilanteisiin tarkastellaan. Analyysien perusteella voimalaitosten ohjaus- ja suojausautomaatio toimivat hyvin myös korotetulla teholla. Tehonkorotuksen jälkeiset suuremmat massavirrat aiheuttavat kuitenkin voimakkaampia reaktoripaineen ja -tehon vaihteluita varsinkin venttiilien sulkeutumistapauksissa. Simuloinnit osoittivat, että tehonkorotus 2500 MW:sta 2750 MW:iin on mahdollinen, mutta aiheuttaa pieniä muutoksia laitoksen suojausjärjestelmien laukaisurajoihin.

Integration of Torrefaction with Steam Power Plant

Torrefaction is one of the pretreatment technologies to enhance the fuel characteristics of biomass. The efficient and continuous operation of a torrefaction reactor, in the commercial scale, demands a secure biomass supply, in addition to adequate source of heat. Biorefinery plants or biomass-fuelled steam power plants have the potential to integrate with the torrefaction reactor to exchange heat and mass, using available infrastructure and energy sources. The technical feasibility of this integration is examined in this study. A new model for the torrefaction process is introduced and verified by the available experimental data. The torrefaction model is then integrated in different steam power plants to simulate possible mass and energy exchange between the reactor and the plants. The performance of the integrated plant is investigated for different configurations and the results are compared.

Towards sustainable Iron- and steelmaking with economic optimization

The iron and steelmaking industry is among the major contributors to the anthropogenic emissions of carbon dioxide in the world. The rising levels of CO2 in the atmosphere and the global concern about the greenhouse effect and climate change have brought about considerable investigations on how to reduce the energy intensity and CO2 emissions of this industrial sector. In this thesis the problem is tackled by mathematical modeling and optimization using three different approaches. The possibility to use biomass in the integrated steel plant, particularly as an auxiliary reductant in the blast furnace, is investigated. By pre-processing the biomass its heating value and carbon content can be increased at the same time as the oxygen content is decreased. As the compression strength of the preprocessed biomass is lower than that of coke, it is not suitable for replacing a major part of the coke in the blast furnace burden. Therefore the biomass is assumed to be injected at the tuyere level of the blast furnace. Carbon capture and storage is, nowadays, mostly associated with power plants but it can also be used to reduce the CO2 emissions of an integrated steel plant. In the case of a blast furnace, the effect of CCS can be further increased by recycling the carbon dioxide stripped top gas back into the process. However, this affects the economy of the integrated steel plant, as the amount of top gases available, e.g., for power and heat production is decreased. High quality raw materials are a prerequisite for smooth blast furnace operation. High quality coal is especially needed to produce coke with sufficient properties to ensure proper gas permeability and smooth burden descent. Lower quality coals as well as natural gas, which some countries have in great volumes, can be utilized with various direct and smelting reduction processes. The DRI produced with a direct reduction process can be utilized as a feed material for blast furnace, basic oxygen furnace or electric arc furnace. The liquid hot metal from a smelting reduction process can in turn be used in basic oxygen furnace or electric arc furnace. The unit sizes and investment costs of an alternative ironmaking process are also lower than those of a blast furnace. In this study, the economy of an integrated steel plant is investigated by simulation and optimization. The studied system consists of linearly described unit processes from coke plant to steel making units, with a more detailed thermodynamical model of the blast furnace. The results from the blast furnace operation with biomass injection revealed the importance of proper pre-processing of the raw biomass as the composition of the biomass as well as the heating value and the yield are all affected by the pyrolysis temperature. As for recycling of CO2 stripped blast furnace top gas, substantial reductions in the emission rates are achieved if the stripped CO2 can be stored. However, the optimal recycling degree together with other operation conditions is heavily dependent on the cost structure of CO2 emissions and stripping/storage. The economical feasibility related to the use of DRI in the blast furnace depends on the price ratio between the DRI pellets and the BF pellets. The high amount of energy needed in the rotary hearth furnace to reduce the iron ore leads to increased CO2 emissions.

Höyryjärjestelmän pelisimulaattorikonsepti

Tässä diplomityössä on suunniteltu konsepti myöhemmin toteutettavaa pelisimulaat-toria varten. Pelisimulaattorista tulee osa öljynjalostamon tuotanto-operaattoreiden höyryjärjestelmän koulutusmoduulia. Työssä höyryjärjestelmän simulaattoriin on laadittu selkeä malli ja rakenne, joka arvioi käyttäjän suoritusta sekä antaa palautteen suorituksen onnistumisesta. Höyryjärjestelmän toiminnan ymmärtäminen on tärkeää, sillä sen avulla höyryntuotantoa voidaan optimoida. Kirjallisuusosassa on käsitelty koulutussimulaattoreita sekä niiden käyttökohteita. Lisäksi on esitetty höyryntuottajat sekä höyrynkuluttajat, jotka huomioitiin pelisimu-laattorikonseptissa. Jokaiselle yksikköprosessille on laadittu toimintakuvaus ja ope-rointiohjeet, joiden perusteella pelaaja saa käsityksen simulaattorin toimintaperiaat-teesta. Soveltavassa osassa on käsitelty kaikkien konseptissa mallinnettujen laitteiden aine- ja energiataseet sekä muut laskennassa tarvittavat teknilliset korrelaatiot. Soveltavas-sa osassa esitettiin myös pelisimulaattorin arviointiperusteet, arviointiin tarvittavien kustannusten korrelaatiot, häiriötilanne sekä esimerkkisimulointi ja suoritusesimerkki häiriön korjaamiseksi. Pelisuorituksen arvioinnin suunnitteluun kiinnitettiin huomiota, jolloin pelisimulaat-tori vastaa mahdollisimman hyvin todellista prosessia sekä säilyttää motivaation si-mulaattoriharjoitteluun.

Role of angiotensin II and vasopressin receptors within the supraoptic nucleus in water and sodium intake induced by the injection of angiotensin II into the medial septal area

In this study we investigated the effects of the injection into the supraoptic nucleus (SON) of non-peptide AT1- and AT2-angiotensin II (ANG II) receptor antagonists, DuP753 and PD123319, as well as of the arginine-vasopressin (AVP) receptor antagonist d(CH2)5-Tyr(Me)-AVP, on water and 3% NaCl intake induced by the injection of ANG II into the medial septal area (MSA). The effects on water or 3% NaCl intake were assessed in 30-h water-deprived or in 20-h water-deprived furosemide-treated adult male rats, respectively. The drugs were injected in 0.5 µl over 30-60 s. Controls were injected with a similar volume of 0.15 M NaCl. Antagonists were injected at doses of 20, 80 and 180 nmol. Water and sodium intake was measured over a 2-h period. Previous administration of the AT1 receptor antagonist DuP753 into the SON decreased water (65%, N = 10, P<0.01) and sodium intake (81%, N = 8, P<0.01) induced by the injection of ANG II (10 nmol) into the MSA. Neither of these responses was significantly changed by injection of the AT2-receptor antagonist PD123319 into the SON. On the other hand, while there was a decrease in water intake (45%, N = 9, P<0.01), ANG II-induced sodium intake was significantly increased (70%, N = 8, P<0.01) following injection of the V1-type vasopressin antagonist d(CH2)5-Tyr(Me)-AVP into the SON. These results suggest that both AT1 and V1 receptors within the SON may be involved in water and sodium intake induced by the activation of ANG II receptors within the MSA. Furthermore, they do not support the involvement of MSA AT2 receptors in the mediation of these responses.

Analysis and modelling of chemical looping combustion process with and without oxygen uncoupling

Effective control and limiting of carbon dioxide (CO₂) emissions in energy production are major challenges of science today. Current research activities include the development of new low-cost carbon capture technologies, and among the proposed concepts, chemical combustion (CLC) and chemical looping with oxygen uncoupling (CLOU) have attracted significant attention allowing intrinsic separation of pure CO₂ from a hydrocarbon fuel combustion process with a comparatively small energy penalty. Both CLC and CLOU utilize the well-established fluidized bed technology, but several technical challenges need to be overcome in order to commercialize the processes. Therefore, development of proper modelling and simulation tools is essential for the design, optimization, and scale-up of chemical looping-based combustion systems. The main objective of this work was to analyze the technological feasibility of CLC and CLOU processes at different scales using a computational modelling approach. A onedimensional fluidized bed model frame was constructed and applied for simulations of CLC and CLOU systems consisting of interconnected fluidized bed reactors. The model is based on the conservation of mass and energy, and semi-empirical correlations are used to describe the hydrodynamics, chemical reactions, and transfer of heat in the reactors. Another objective was to evaluate the viability of chemical looping-based energy production, and a flow sheet model representing a CLC-integrated steam power plant was developed. The 1D model frame was succesfully validated based on the operation of a 150 kWth laboratory-sized CLC unit fed by methane. By following certain scale-up criteria, a conceptual design for a CLC reactor system at a pre-commercial scale of 100 MWth was created, after which the validated model was used to predict the performance of the system. As a result, further understanding of the parameters affecting the operation of a large-scale CLC process was acquired, which will be useful for the practical design work in the future. The integration of the reactor system and steam turbine cycle for power production was studied resulting in a suggested plant layout including a CLC boiler system, a simple heat recovery setup, and an integrated steam cycle with a three pressure level steam turbine. Possible operational regions of a CLOU reactor system fed by bituminous coal were determined via mass, energy, and exergy balance analysis. Finally, the 1D fluidized bed model was modified suitable for CLOU, and the performance of a hypothetical 500 MWth CLOU fuel reactor was evaluated by extensive case simulations.

Comparison of Different Concentrated Solar Power Collector Designs and Development of a Linear Fresnel Solar Collector Model

Concentrated solar power (CSP) is a renewable energy technology, which could contribute to overcoming global problems related to pollution emissions and increasing energy demand. CSP utilizes solar irradiation, which is a variable source of energy. In order to utilize CSP technology in energy production and reliably operate a solar field including thermal energy storage system, dynamic simulation tools are needed in order to study the dynamics of the solar field, to optimize production and develop control systems. The object of this Master’s Thesis is to compare different concentrated solar power technologies and configure a dynamic solar field model of one selected CSP field design in the dynamic simulation program Apros, owned by VTT and Fortum. The configured model is based on German Novatec Solar’s linear Fresnel reflector design. Solar collector components including dimensions and performance calculation were developed, as well as a simple solar field control system. The preliminary simulation results of two simulation cases under clear sky conditions were good; the desired and stable superheated steam conditions were maintained in both cases, while, as expected, the amount of steam produced was reduced in the case having lower irradiation conditions. As a result of the model development process, it can be concluded, that the configured model is working successfully and that Apros is a very capable and flexible tool for configuring new solar field models and control systems and simulating solar field dynamic behaviour.

Delayed Schwann cell and oligodendrocyte remyelination after ethidium bromide injection in the brainstem of Wistar rats submitted to streptozotocin diabetogenic treatment

Schwann cell disturbance followed by segmental demyelination in the peripheral nervous system occurs in diabetic patients. Since Schwann cell and oligodendrocyte remyelination in the central nervous system is a well-known event in the ethidium bromide (EB) demyelinating model, the aim of this investigation was to determine the behavior of both cell types after local EB injection into the brainstem of streptozotocin diabetic rats. Adult male Wistar rats received a single intravenous injection of streptozotocin (50 mg/kg) and were submitted 10 days later to a single injection of 10 µL 0.1% (w/v) EB or 0.9% saline solution into the cisterna pontis. Ten microliters of 0.1% EB was also injected into non-diabetic rats. The animals were anesthetized and perfused through the heart 7 to 31 days after EB or saline injection and brainstem sections were collected and processed for light and transmission electron microscopy. The final balance of myelin repair in diabetic and non-diabetic rats at 31 days was compared using a semi-quantitative method. Diabetic rats presented delayed macrophage activity and lesser remyelination compared to non-diabetic rats. Although oligodendrocytes were the major remyelinating cells in the brainstem, Schwann cells invaded EB-induced lesions, first appearing at 11 days in non-diabetic rats and by 15 days in diabetic rats. Results indicate that short-term streptozotocin-induced diabetes hindered both oligodendrocyte and Schwann cell remyelination (mean remyelination scores of 2.57 ± 0.77 for oligodendrocytes and 0.67 ± 0.5 for Schwann cells) compared to non-diabetic rats (3.27 ± 0.85 and 1.38 ± 0.81, respectively).

Flow in porous media of oil-water systems

Fluid flow behaviour in porous media is a conundrum. Therefore, this research is focused on filtration-volumetric characterisation of fractured-carbonate sediments, coupled with their proper simulation. For this reason, at laboratory rock properties such as pore volume, permeability and porosity are measured, later phase permeabilities and oil recovery in function of flow rate are assessed. Furthermore, the rheological properties of three oils are measured and analysed. Finally based on rock and fluid properties, a model using COMSOL Multiphysics is built in order to compare the experimental and simulated results. The rock analyses show linear relation between flow rate and differential pressure, from which phase permeabilities and pressure gradient are determined, eventually the oil recovery under low and high flow rate is established. In addition, the oils reveal thixotropic properties as well as non-Newtonian behaviour described by Bingham model, consequently Carreau viscosity model for the used oil is given. Given these points, the model for oil and water is built in COMSOL Multiphysics, whereupon successfully the reciprocity between experimental and simulated results is analysed and compared. Finally, a two-phase displacement model is elaborated.