A Soft Contact Collision Method For Real-Time Simulation of Triangularized Geometries in Multibody Dynamics

When modeling machines in their natural working environment collisions become a very important feature in terms of simulation accuracy. By expanding the simulation to include the operation environment, the need for a general collision model that is able to handle a wide variety of cases has become central in the development of simulation environments. With the addition of the operating environment the challenges for the collision modeling method also change. More simultaneous contacts with more objects occur in more complicated situations. This means that the real-time requirement becomes more difficult to meet. Common problems in current collision modeling methods include for example dependency on the geometry shape or mesh density, calculation need increasing exponentially in respect to the number of contacts, the lack of a proper friction model and failures due to certain configurations like closed kinematic loops. All these problems mean that the current modeling methods will fail in certain situations. A method that would not fail in any situation is not very realistic but improvements can be made over the current methods.

Adaptive Routing Approaches for Networked Many-Core Systems

Through advances in technology, System-on-Chip design is moving towards integrating tens to hundreds of intellectual property blocks into a single chip. In such a many-core system, on-chip communication becomes a performance bottleneck for high performance designs. Network-on-Chip (NoC) has emerged as a viable solution for the communication challenges in highly complex chips. The NoC architecture paradigm, based on a modular packet-switched mechanism, can address many of the on-chip communication challenges such as wiring complexity, communication latency, and bandwidth. Furthermore, the combined benefits of 3D IC and NoC schemes provide the possibility of designing a high performance system in a limited chip area. The major advantages of 3D NoCs are the considerable reductions in average latency and power consumption. There are several factors degrading the performance of NoCs. In this thesis, we investigate three main performance-limiting factors: network congestion, faults, and the lack of efficient multicast support. We address these issues by the means of routing algorithms. Congestion of data packets may lead to increased network latency and power consumption. Thus, we propose three different approaches for alleviating such congestion in the network. The first approach is based on measuring the congestion information in different regions of the network, distributing the information over the network, and utilizing this information when making a routing decision. The second approach employs a learning method to dynamically find the less congested routes according to the underlying traffic. The third approach is based on a fuzzy-logic technique to perform better routing decisions when traffic information of different routes is available. Faults affect performance significantly, as then packets should take longer paths in order to be routed around the faults, which in turn increases congestion around the faulty regions. We propose four methods to tolerate faults at the link and switch level by using only the shortest paths as long as such path exists. The unique characteristic among these methods is the toleration of faults while also maintaining the performance of NoCs. To the best of our knowledge, these algorithms are the first approaches to bypassing faults prior to reaching them while avoiding unnecessary misrouting of packets. Current implementations of multicast communication result in a significant performance loss for unicast traffic. This is due to the fact that the routing rules of multicast packets limit the adaptivity of unicast packets. We present an approach in which both unicast and multicast packets can be efficiently routed within the network. While suggesting a more efficient multicast support, the proposed approach does not affect the performance of unicast routing at all. In addition, in order to reduce the overall path length of multicast packets, we present several partitioning methods along with their analytical models for latency measurement. This approach is discussed in the context of 3D mesh networks.

On Energy Efficient Computing Platforms

In accordance with the Moore's law, the increasing number of on-chip integrated transistors has enabled modern computing platforms with not only higher processing power but also more affordable prices. As a result, these platforms, including portable devices, work stations and data centres, are becoming an inevitable part of the human society. However, with the demand for portability and raising cost of power, energy efficiency has emerged to be a major concern for modern computing platforms. As the complexity of on-chip systems increases, Network-on-Chip (NoC) has been proved as an efficient communication architecture which can further improve system performances and scalability while reducing the design cost. Therefore, in this thesis, we study and propose energy optimization approaches based on NoC architecture, with special focuses on the following aspects. As the architectural trend of future computing platforms, 3D systems have many bene ts including higher integration density, smaller footprint, heterogeneous integration, etc. Moreover, 3D technology can signi cantly improve the network communication and effectively avoid long wirings, and therefore, provide higher system performance and energy efficiency. With the dynamic nature of on-chip communication in large scale NoC based systems, run-time system optimization is of crucial importance in order to achieve higher system reliability and essentially energy efficiency. In this thesis, we propose an agent based system design approach where agents are on-chip components which monitor and control system parameters such as supply voltage, operating frequency, etc. With this approach, we have analysed the implementation alternatives for dynamic voltage and frequency scaling and power gating techniques at different granularity, which reduce both dynamic and leakage energy consumption. Topologies, being one of the key factors for NoCs, are also explored for energy saving purpose. A Honeycomb NoC architecture is proposed in this thesis with turn-model based deadlock-free routing algorithms. Our analysis and simulation based evaluation show that Honeycomb NoCs outperform their Mesh based counterparts in terms of network cost, system performance as well as energy efficiency.

Development of porous glass-fiber reinforced composite for bone implants. Evaluation of antimicrobial effect and implant fixation

Cranial bone reconstructions are necessary for correcting large skull bone defects due to trauma, tumors, infections and craniotomies. Traditional synthetic implant materials include solid or mesh titanium, various plastics and ceramics. Recently, biostable glass-fiber reinforced composites (FRC), which are based on bifunctional methacrylate resin, were introduced as novel implant solution. FRCs were originally developed and clinically used in dental applications. As a result of further in vitro and in vivo testing, these composites were also approved for clinical use in cranial surgery. To date, reconstructions of large bone defects were performed in 35 patients. This thesis is dedicated to the development of a novel FRC-based implant for cranial reconstructions. The proposed multi-component implant consists of three main parts: (i) porous FRC structure; (ii) bioactive glass granules embedded between FRC layers and (iii) a silver-polysaccharide nanocomposite coating. The porosity of the FRC structure should allow bone ingrowth. Bioactive glass as an osteopromotive material is expected to stimulate the formation of new bone. The polysaccharide coating is expected to prevent bacterial colonization of the implant. The FRC implants developed in this study are based on the porous network of randomly-oriented E-glass fibers bound together by non-resorbable photopolymerizable methacrylate resin. These structures had a total porosity of 10–70 volume %, of which > 70% were open pores. The pore sizes > 100 μm were in the biologically-relevant range (50-400 μm), which is essential for vascularization and bone ingrowth. Bone ingrowth into these structures was simulated by imbedding of porous FRC specimens in gypsum. Results of push-out tests indicated the increase in the shear strength and fracture toughness of the interface with the increase in the total porosity of FRC specimens. The osteopromotive effect of bioactive glass is based on its dissolution in the physiological environment. Here, calcium and phosphate ions, released from the glass, precipitated on the glass surface and its proximity (the FRC) and formed bone-like apatite. The biomineralization of the FRC structure, due to the bioactive glass reactions, was studied in Simulated Body Fluid (SBF) in static and dynamic conditions. An antimicrobial, non-cytotoxic polysaccharide coating, containing silver nanoparticles, was obtained through strong electrostatic interactions with the surface of FRC. In in vitro conditions the lactose-modified chitosan (chitlac) coating showed no signs of degradation within seven days of exposure to lysozyme or one day to hydrogen peroxide (H2O2). The antimicrobial efficacy of the coating was tested against Staphylococcus aureus and Pseudomonas aeruginosa. The contact-active coating had an excellent short time antimicrobial effect. The coating neither affected the initial adhesion of microorganisms to the implant surface nor the biofilm formation after 24 h and 72 h of incubation. Silver ions released to the aqueous environment led to a reduction of bacterial growth in the culture medium.

Höyryturbiinin kanavistojen numeerinen virtauslaskenta

Diplomityössä tutkittiin Fortumin Loviisan ydinvoimalaitoksen ulosvirtauskanaviston ja suurnopeuskosteudenerottimen toimintaa, sekä selvitettiin taustalla olevaa teoriaa ja aiemmin tehtyjä tutkimuksia. Tavoitteena oli ymmärtää ja esittää laitteiden toimintaa, sekä tutkia voiko ulosvirtauskanaviston suorituskykyä parantaa geometrian muutoksilla. Työssä luotiin tutkittaville kohteille geometriat ja laskentahilat, joiden avulla simuloitiin niiden toimintaa eri käyttötilanteissa numeerisen virtauslaskennan avulla. Laskennan reunaehdot saatiin olemassa olevasta prosessimallista ja aiemmista turbiiniselvityksistä. Ulosvirtauskanaviston suorituskyky laskettiin kolmella eri lauhdutinpaineella neljällä eri geometrialla. Geometrian muutokset vaikuttivat selkeästi ulosvirtauskanaviston suorituskykyyn ja sitä saatiin parannettua. Kaksi kolmesta muutoksesta, lisäkanavat ja oikaistu vesilippa, pa-ransivat suorituskykyä. Lokinsiipien poistaminen heikensi ulosvirtauskanaviston toi-mintaa. Suurnopeuskosteudenerottimen mallintaminen jäi lähtötietojen ja ajan puutteen takia hieman tavoitteesta. Sekä ulosvirtauskanaviston että suurnopeuskosteudenerotti-men jatkotutkimusta ja mahdollisia toimenpiteitä varten saatiin arvokasta uutta tietoa.

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.

Two-phase flow measurements with spacer grid in a rod bundle geometry

Heat transfer effectiveness in nuclear rod bundles is of great importance to nuclear reactor safety and economics. An important design parameter is the Critical Heat Flux (CHF), which limits the transferred heat from the fuel to the coolant. The CHF is determined by flow behaviour, especially the turbulence created inside the fuel rod bundle. Adiabatic experiments can be used to characterize the flow behaviour separately from the heat transfer phenomena in diabatic flow. To enhance the turbulence, mixing vanes are attached to spacer grids, which hold the rods in place. The vanes either make the flow swirl around a single sub-channel or induce cross-mixing between adjacent sub-channels. In adiabatic two-phase conditions an important phenomenon that can be investigated is the effect of the spacer on canceling the lift force, which collects the small bubbles to the rod surfaces leading to decreased CHF in diabatic conditions and thus limits the reactor power. Computational Fluid Dynamics (CFD) can be used to simulate the flow numerically and to test how different spacer configurations affect the flow. Experimental data is needed to validate and verify the used CFD models. Especially the modeling of turbulence is challenging even for single-phase flow inside the complex sub-channel geometry. In two-phase flow other factors such as bubble dynamics further complicate the modeling. To investigate the spacer grid effect on two-phase flow, and to provide further experimental data for CFD validation, a series of experiments was run on an adiabatic sub-channel flow loop using a duct-type spacer grid with different configurations. Utilizing the wire-mesh sensor technology, the facility gives high resolution experimental data in both time and space. The experimental results indicate that the duct-type spacer grid is less effective in canceling the lift force effect than the egg-crate type spacer tested earlier.

Simulation model of wind turbine gearbox and lubrication system

As increasing efficiency of a wind turbine gearbox, more power can be transferred from rotor blades to generator and less power is used to cause wear and heating in the gearbox. By using a simulation model, behavior of the gearbox can be studied before creating expensive prototypes. The objective of the thesis is to model a wind turbine gearbox and its lubrication system to study power losses and heat transfer inside the gearbox and to study the simulation methods of the used software. Software used to create the simulation model is Siemens LMS Imagine.Lab AMESim, which can be used to create one-dimensional mechatronic system simulation models from different fields of engineering. When combining components from different libraries it is possible to create a simulation model, which includes mechanical, thermal and hydraulic models of the gearbox. Results for mechanical, thermal, and hydraulic simulations are presented in the thesis. Due to the large scale of the wind turbine gearbox and the amount of power transmitted, power loss calculations from AMESim software are inaccurate and power losses are modelled as constant efficiency for each gear mesh. Starting values for simulation in thermal and hydraulic simulations were chosen from test measurements and from empirical study as compact and complex design of gearbox prevents accurate test measurements. In further studies to increase the accuracy of the simulation model, components used for power loss calculations needs to be modified and values for unknown variables are needed to be solved through accurate test measurements.

Elementtikoon vaikutus väsymismenetelmien tarkkuuteen

Kandidaatintyössä tutkittiin kolmen hitsatun palstalevyliitoksen väsymiskestävyyttä äärellisten elementtien menetelmällä. Tutkimuksen pääpaino oli 2D-mallin antamien tulosten suppenevuudessa kolmella eri elementtikoolla ja kolmessa eri mittasuhdetapauksessa. Tutkimustuloksissa havaittiin noin 5 % muutos harvimman ja tiheimmän elementtiverkon tulosten välillä. Tutkimuksessa havaittiin myös että muutos harvimman ja tiheimmän elementtiverkon tuloksissa pieneni hiukan mitä suurempia kappaleet ja niiden hitsausliitokset olivat.

Simulation and material calibration of ultra high strength steel (UHSS) S960 welded joint under static tensile test utilizing finite element method (FEM)

The aim of this work was to calibrate the material properties including strength and strain values for different material zones of ultra-high strength steel (UHSS) welded joints under monotonic static loading. The UHSS is heat sensitive and softens by heat due to welding, the affected zone is heat affected zone (HAZ). In this regard, cylindrical specimens were cut out from welded joints of Strenx® 960 MC and Strenx® Tube 960 MH, were examined by tensile test. The hardness values of specimens’ cross section were measured. Using correlations between hardness and strength, initial material properties were obtained. The same size specimen with different zones of material same as real specimen were created and defined in finite element method (FEM) software with commercial brand Abaqus 6.14-1. The loading and boundary conditions were defined considering tensile test values. Using initial material properties made of hardness-strength correlations (true stress-strain values) as Abaqus main input, FEM is utilized to simulate the tensile test process. By comparing FEM Abaqus results with measured results of tensile test, initial material properties will be revised and reused as software input to be fully calibrated in such a way that FEM results and tensile test results deviate minimum. Two type of different S960 were used including 960 MC plates, and structural hollow section 960 MH X-joint. The joint is welded by BöhlerTM X96 filler material. In welded joints, typically the following zones appear: Weld (WEL), Heat affected zone (HAZ) coarse grained (HCG) and fine grained (HFG), annealed zone, and base material (BaM). Results showed that: The HAZ zone is softened due to heat input while welding. For all the specimens, the softened zone’s strength is decreased and makes it a weakest zone where fracture happens while loading. Stress concentration of a notched specimen can represent the properties of notched zone. The load-displacement diagram from FEM modeling matches with the experiments by the calibrated material properties by compromising two correlations of hardness and strength.