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.

Real-Time Simulation of Mobile and Industrial Machines Using the Multibody Simulation Approach

This thesis introduces a real-time simulation environment based on the multibody simulation approach. The environment consists of components that are used in conventional product development, including computer aided drawing, visualization, dynamic simulation and finite element software architecture, data transfer and haptics. These components are combined to perform as a coupled system on one platform. The environment is used to simulate mobile and industrial machines at different stages of a product life time. Consequently, the demands of the simulated scenarios vary. In this thesis, a real-time simulation environment based on the multibody approach is used to study a reel mechanism of a paper machine and a gantry crane. These case systems are used to demonstrate the usability of the real-time simulation environment for fault detection purposes and in the context of a training simulator. In order to describe the dynamical performance of a mobile or industrial machine, the nonlinear equations of motion must be defined. In this thesis, the dynamical behaviour of machines is modelled using the multibody simulation approach. A multibody system may consist of rigid and flexible bodies which are joined using kinematic joint constraints while force components are used to describe the actuators. The strength of multibody dynamics relies upon its ability to describe nonlinearities arising from wearing of the components, friction, large rotations or contact forces in a systematic manner. For this reason, the interfaces between subsystems such as mechanics, hydraulics and control systems of the mechatronic machine can be defined and analyzed in a straightforward manner.

Acoustic environments in change. Five village soundscapes

Five village soundscapes toim./ed. R. M. Schafer. kirjan liitteinä olevista CD-levyistä kaksi liittyy teokseen Acoustic environments in change ja kaksi teokseen Five village soundscapes. Five villages soundscapes. 1st ed. Vancouver : A.R.C. Publications, cop. 1977.

A comparison of electricity spot prices simulation using ARMA-GARCH and mean-reverting models

The aim of this work is to compare two families of mathematical models for their respective capability to capture the statistical properties of real electricity spot market time series. The first model family is ARMA-GARCH models and the second model family is mean-reverting Ornstein-Uhlenbeck models. These two models have been applied to two price series of Nordic Nord Pool spot market for electricity namely to the System prices and to the DenmarkW prices. The parameters of both models were calibrated from the real time series. After carrying out simulation with optimal models from both families we conclude that neither ARMA-GARCH models, nor conventional mean-reverting Ornstein-Uhlenbeck models, even when calibrated optimally with real electricity spot market price or return series, capture the statistical characteristics of the real series. But in the case of less spiky behavior (System prices), the mean-reverting Ornstein-Uhlenbeck model could be seen to partially succeeded in this task.

3D Stripixel Detectors Simulation

Solid-state silicon detectors have replaced conventional ones in almost all recent high-energy physics experiments. Pixel silicon sensors don't have any alternative in the area near the interaction point because of their high resolution and fast operation speed. However, present detectors hardly withstand high radiation doses. Forthcoming upgrade of the LHC in 2014 requires development of a new generation of pixel detectors which will be able to operate under ten times increased luminosity. A planar fabrication technique has some physical limitations; an improvement of the radiation hardness will reduce sensitivity of a detector. In that case a 3D pixel detector seems to be the most promising device which can overcome these difficulties. The objective of this work was to model a structure of the 3D stripixel detector and to simulate electrical characteristics of the device. Silvaco Atlas software has been used for these purposes. The structures of single and double sided dual column detectors with active edges were described using special command language. Simulations of these detectors have shown that electric field inside an active area has more uniform distribution in comparison to the planar structure. A smaller interelectrode space leads to a stronger field and also decreases the collection time. This makes the new type of detectors more radiation resistant. Other discovered advantages are the lower full depletion voltage and increased charge collection efficiency. So the 3D stripixel detectors have demonstrated improved characteristics and will be a suitable replacement for the planar ones.

Past, present and future - assessing changes in aquatic environments using subfossil diatoms

The deterioration of surface waters is one of the most important issues in the environmental management of the European Union. Thus, the EU Water Framework Directive 2000/60/EC (WFD) requires “good ecological and chemical status” of surface waters by 2015 allowing only a slight departure from ecological reference conditions characterized by the biological communities typical for the conditions of minimal anthropogenic impact. The WFD requires the determination of ecological reference conditions and the present ecological status of surface waters. To meet this legislative demand, sedimentary diatom assemblages were used in these studies with various methods 1) to assess natural and human activity induced environmental changes, 2) to characterize background conditions 3) to evaluate the present ecological status and 4) to predict the future of the water bodies in the light of palaeolimnological data. As the WFD refers to all surface waters, both coastal and inland sites were included. Two long and two short sediment cores from the Archipelago Sea in the northern Baltic Sea were examined for their siliceous microfossils in order to assess (1) the Holocene palaeoenvironmental history and (2) the recent eutrophication of the area. The diatom record was divided into local diatom assemblage zones (LDAZ, long cores) and diatom assemblage zones (DAZ, short cores). Locally weighted weighted averaging regression and calibration (LWWA) was applied for the quantitative reconstruction of past TN concentrations (short cores). An age model for the long cores was constructed by using independent palaeomagnetic and AMS-14C methods. The short cores were dated using radiometric (210Pb, 226Ra and 137Cs) methods. The long cores date back to the early history of the Archipelago Sea, which was freshwater – no salinity increase referable to the brackish phase of the Yoldia Sea is recognized. The nutrient status of the lacustrine phase was slightly higher in the Archipelago Sea than in the Baltic Proper. Initial brackish-water influence is observed at 8 150 ±80 cal. BP (LDAZ4), but fully brackish conditions were established at 7 700 ±80 cal. BP (LDAZ5). The diatom assemblages indicate increasing salinity, warming climate and possible eutrophic conditions during the lacustrine to brackish-water transition. The decreasing abundance of Pseudosolenia calcar-avis (Schultze) Sundström and the increasing abundance of the ice-cover indicator species Pauliella taeniata (Grunow) Round and Basson indicate decreasing salinity and climatic cooling after ca. 5 000 cal. BP. Signs of eutrophication are visible in the most recent diatom assemblage zones of both short cores. Diatom-inferred total nitrogen (DI-TN) reconstructions partially fail to trace the actual measured total nitrogen concentrations especially from the late 1980s to the mid 1990s. This is most likely due to the dominating diatom species Pauliella taeniata, Thalassiosira levanderi Van Goor and Fragilariopsis cylindrus (Grunow) W. Krieger being more influenced by factors such as the length of the ice-season rather than nutrient concentrations. It is concluded that the diatom assemblages of the study sites are principally governed by climate fluctuations, with a slight influence of eutrophication visible in the most recent sediments. There are indications that global warming, with reduced ice cover, could impact the spring blooming diatom species composition in the Archipelago Sea. In addition, increased sediment accumulation in the early 90s coincides with the short ice-seasons suggesting that warming climate with decreasing ice-cover may increase sedimentation in the study area. The diverse diatom assemblages dominated by benthic species (54 %) in DAZ1 in the Käldö Fjärd core can be taken as background diatom assemblages for the Archipelago Sea. Since then turbidity has increased and the diatom assemblages have been dominated by planktonic diatoms from around the mid 1800s onwards. The reconstructed reference conditions for the total nitrogen concentrations fluctuate around 400 μg l-1. Altogether two short sediment cores and eight short cores for top-bottom analysis were retrieved from Lake Orijärvi and Lake Määrjärvi to assess the impact of the acid mine drainage (AMD) derived metals from the Orijärvi mine tailings on the diatom communities of the lakes. The Cu (Pb, Zn) mine of Orijärvi (1757 – 1956) was the first one in Finland where flotation techniques (1911 – 1955) were used to enrich ore and large quantities of tailings were produced. The AMD derived metal impact to the lakes was found to be among the heaviest thus far recorded in Finland. Concentrations of Cu, Pb and Zn in Lake Orijärvi sediments are two to three orders of magnitude higher than background values. The metal inputs have affected Lake Orijärvi and Lake Määrjärvi diatom communities at the community levels through shifts in dominant taxa (both lakes) and at the individual level through alteration in frustule morphology (Lake Orijärvi). At present, lake water still has elevated heavy metal levels, indicating that the impact from the tailings area continues to affect both lakes. Lake Orijärvi diatom assemblages are completely dominated by benthic species and are lacking planktonic diatoms. In Lake Määrjärvi the proportion of benthic and tychoplanktonic diatoms has increased and the planktonic taxa have decreased in abundance. Achnanthidium minutissimum Kützing and Brachysira vitrea (Grun.) R. Ross in Hartley were the most tolerant species to increased metal concentrations. Planktonic diatoms are more sensitive to metal contamination than benthic taxa, especially species in the genus Cyclotella (Kützing) Brébisson. The ecological reference conditions assessed in this study for Lake Orijärvi and Lake Määrjärvi comprise diverse planktonic and benthic communitites typical of circumneutral oligotrophic lakes, where the planktonic diatoms belonging to genera Cyclotella , Aulacoseira Thwaites, Tabellaria Ehrenberg and Asterionella Hassall dominate in relative abundances up to ca. 70%. The benthic communities are more diverse than the planktonic consisting of diatoms belonging to the genera Achnanthes Bory, Fragilaria Lyngbye and Navicula St. Vincent. This study clearly demonstrates that palaeolimnological methods, especially diatom analysis, provide a powerful tool for the EU Water Frame Work Directive for defining reference conditions, natural variability and current status of surface waters. The top/bottom approach is a very useful tool in larger-scale studies needed for management purposes. This “before and after” type of sediment sampling method can provide a very time and cost effective assessment of ecological reference conditions of surface waters.

Induction Motor Drive Energy Efficiency - Simulation and Analysis

A coupled system simulator, based on analytical circuit equations and a finite element method (FEM) model of the motor has been developed and it is used to analyse a frequency-converterfed industrial squirrel-cage induction motor. Two control systems that emulate the behaviour of commercial direct-torque-controlled (DTC) and vector-controlled industrial frequency converters have been studied, implemented in the simulation software and verified by extensive laboratory tests. Numerous factors that affect the operation of a variable speed drive (VSD) and its energy efficiency have been investigated, and their significance in the simulation of the VSD results has been studied. The dependency of the frequency converter, induction motor and system losses on the switching frequency is investigated by simulations and measurements at different speeds for both the vector control and the DTC. Intensive laboratory measurements have been carried out to verify the simulation results.

CFD Simulation of Two-phase and Three-phase Flows in Internal-loop Airlift

Airlift reactors are pneumatically agitated reactors that have been widely used in chemical, petrochemical, and bioprocess industries, such as fermentation and wastewater treatment. Computational Fluid Dynamics (CFD) has become more popular approach for design, scale-up and performance evaluation of such reactors. In the present work numerical simulations for internal-loop airlift reactors were performed using the transient Eulerian model with CFD package, ANSYS Fluent 12.1. The turbulence in the liquid phase is described using κ- ε the model. Global hydrodynamic parameters like gas holdup, gas velocity and liquid velocity have been investigated for a range of superficial gas velocities, both with 2D and 3D simulations. Moreover, the study of geometry and scale influence on the reactor have been considered. The results suggest that both, geometry and scale have significant effects on the hydrodynamic parameters, which may have substantial effects on the reactor performance. Grid refinement and time-step size effect have been discussed. Numerical calculations with gas-liquid-solid three-phase flow system have been carried out to investigate the effect of solid loading, solid particle size and solid density on the hydrodynamic characteristics of internal loop airlift reactor with different superficial gas velocities. It was observed that averaged gas holdup is significantly decreased with increasing slurry concentration. Simulations show that the riser gas holdup decreases with increase in solid particle diameter. In addition, it was found that the averaged solid holdup increases in the riser section with the increase of solid density. These produced results reveal that CFD have excellent potential to simulate two-phase and three-phase flow system.

Simulation of hydrogen peroxide direct synthesis in a microreactor

A set of models in Aspen plus was built to simulate the direct synthesis process of hydrogen peroxide in a micro-reactor system. This process model can be used to carry out material balance calculation under various experimental conditions. Three thermodynamic property methods were compared by calculating gas solubility and Uniquac-RK method was finally selected for process model. Two different operation modes with corresponding operation conditions were proposed as the starting point of future experiments. Simulations for these two modes were carried out to get the information of material streams. Moreover, some hydrodynamic parameters such as gas/liquid superficial velocity, gas holdup were also calculated with improved process model. These parameters proved the proposed experimental conditions reasonable to some extent. The influence of operation conditions including temperature, pressure and circulation ratio was analyzed for the first operation mode, where pure oxygen was fed into dissolving tank and hydrogen-carbon dioxide mixture was fed into microreactor directly. The preferred operation conditions for the system are low temperature (2°C) and high pressure (30 bar) in dissolving tank. High circulation ratio might be good in the sense that more oxygen could be dissolved and fed into reactor for reactions, but meanwhile hydrodynamics of microreactor should be considered. Furthermore, more operation conditions of reactor gas/liquid feeds in both of two operation modes were proposed to provide guidance for future experiment design and corresponding hydrodynamic parameters were also calculated. Finally, safety issue was considered from thermodynamic point of view and there is no explosion danger at given experimental plan since the released reaction heat will not cause solvent vaporization inside the microchannels. The improvement of process model still needs further study based on the future experimental results.

Microscopic simulation of pigment coating consolidation

The control of coating layer properties is becoming increasingly important as a result of an emerging demand for novel coated paper-based products and an increasing popularity of new coating application methods. The governing mechanisms of microstructure formation dynamics during consolidation and drying are nevertheless, still poorly understood. Some of the difficulties encountered by experimental methods can be overcome by the utilisation of numerical modelling and simulation-based studies of the consolidation process. The objective of this study was to improve the fundamental understanding of pigment coating consolidation and structure formation mechanisms taking place on the microscopic level. Furthermore, it is aimed to relate the impact of process and suspension properties to the microstructure of the coating layer. A mathematical model based on a modified Stokesian dynamics particle simulation technique was developed and applied in several studies of consolidation-related phenomena. The model includes particle-particle and particle-boundary hydrodynamics, colloidal interactions, Born repulsion, and a steric repulsion model. The Brownian motion and a free surface model were incorporated to enable the specific investigation of consolidation and drying. Filter cake stability was simulated in various particle systems, and subjected to a range of base substrate absorption rates and system temperatures. The stability of the filter cake was primarily affected by the absorption rate and size of particles. Temperature was also shown to have an influence. The consolidation of polydisperse systems, with varying wet coating thicknesses, was studied using imposed pilot trial and model-based drying conditions. The results show that drying methods have a clear influence on the microstructure development, on small particle distributions in the coating layer and also on the mobility of particles during consolidation. It is concluded that colloidal properties can significantly impact coating layer shrinkage as well as the internal solids concentration profile. Visualisations of particle system development in time and comparison of systems at different conditions are useful in illustrating coating layer structure formation mechanisms. The results aid in understanding the underlying mechanisms of pigment coating layer consolidation. Guidance is given regarding the relationship between coating process conditions and internal coating slurry properties and their effects on the microstructure of the coating.

Tutorials for MeVEA Simulation Software

The objective of the thesis was to create three tutorials for MeVEA Simulation Software to instruct the new users to the modeling methodology used in the MeVEA Simulation Software. MeVEA Simulation Software is a real-time simulation software based on multibody dynamics. The simulation software is designed to create simulation models of complete mechatronical system. The thesis begins with a more detail description of the MeVEA Simulation Software and its components. The thesis presents the three simulation models and written theory of the steps of model creation. The first tutorial introduces the basic features which are used in most simulation models. The basic features include bodies, constrains, forces, basic hydraulics and motors. The second tutorial introduces the power transmission components, tyres and user input definitions for the different components in power transmission systems. The third tutorial introduces the definitions of two different types of collisions and collision graphics used in MeVEA Simulation Software.