Studies on integrating kinematic design method with mechanical systems simulation techniques

Theultimate goal of any research in the mechanism/kinematic/design area may be called predictive design, ie the optimisation of mechanism proportions in the design stage without requiring extensive life and wear testing. This is an ambitious goal and can be realised through development and refinement of numerical (computational) technology in order to facilitate the design analysis and optimisation of complex mechanisms, mechanical components and systems. As a part of the systematic design methodology this thesis concentrates on kinematic synthesis (kinematic design and analysis) methods in the mechanism synthesis process. The main task of kinematic design is to find all possible solutions in the form of structural parameters to accomplish the desired requirements of motion. Main formulations of kinematic design can be broadly divided to exact synthesis and approximate synthesis formulations. The exact synthesis formulation is based in solving n linear or nonlinear equations in n variables and the solutions for the problem areget by adopting closed form classical or modern algebraic solution methods or using numerical solution methods based on the polynomial continuation or homotopy. The approximate synthesis formulations is based on minimising the approximation error by direct optimisation The main drawbacks of exact synthesis formulationare: (ia) limitations of number of design specifications and (iia) failure in handling design constraints- especially inequality constraints. The main drawbacks of approximate synthesis formulations are: (ib) it is difficult to choose a proper initial linkage and (iib) it is hard to find more than one solution. Recentformulations in solving the approximate synthesis problem adopts polynomial continuation providing several solutions, but it can not handle inequality const-raints. Based on the practical design needs the mixed exact-approximate position synthesis with two exact and an unlimited number of approximate positions has also been developed. The solutions space is presented as a ground pivot map but thepole between the exact positions cannot be selected as a ground pivot. In this thesis the exact synthesis problem of planar mechanism is solved by generating all possible solutions for the optimisation process ¿ including solutions in positive dimensional solution sets - within inequality constraints of structural parameters. Through the literature research it is first shown that the algebraic and numerical solution methods ¿ used in the research area of computational kinematics ¿ are capable of solving non-parametric algebraic systems of n equations inn variables and cannot handle the singularities associated with positive-dimensional solution sets. In this thesis the problem of positive-dimensional solutionsets is solved adopting the main principles from mathematical research area of algebraic geometry in solving parametric ( in the mathematical sense that all parameter values are considered ¿ including the degenerate cases ¿ for which the system is solvable ) algebraic systems of n equations and at least n+1 variables.Adopting the developed solution method in solving the dyadic equations in direct polynomial form in two- to three-precision-points it has been algebraically proved and numerically demonstrated that the map of the ground pivots is ambiguousand that the singularities associated with positive-dimensional solution sets can be solved. The positive-dimensional solution sets associated with the poles might contain physically meaningful solutions in the form of optimal defectfree mechanisms. Traditionally the mechanism optimisation of hydraulically driven boommechanisms is done at early state of the design process. This will result in optimal component design rather than optimal system level design. Modern mechanismoptimisation at system level demands integration of kinematic design methods with mechanical system simulation techniques. In this thesis a new kinematic design method for hydraulically driven boom mechanism is developed and integrated in mechanical system simulation techniques. The developed kinematic design method is based on the combinations of two-precision-point formulation and on optimisation ( with mathematical programming techniques or adopting optimisation methods based on probability and statistics ) of substructures using calculated criteria from the system level response of multidegree-of-freedom mechanisms. Eg. by adopting the mixed exact-approximate position synthesis in direct optimisation (using mathematical programming techniques) with two exact positions and an unlimitednumber of approximate positions the drawbacks of (ia)-(iib) has been cancelled.The design principles of the developed method are based on the design-tree -approach of the mechanical systems and the design method ¿ in principle ¿ is capable of capturing the interrelationship between kinematic and dynamic synthesis simultaneously when the developed kinematic design method is integrated with the mechanical system simulation techniques.

Real-time simulation of multibody systems with applications for working mobile vehicles

This dissertation describes an approach for developing a real-time simulation for working mobile vehicles based on multibody modeling. The use of multibody modeling allows comprehensive description of the constrained motion of the mechanical systems involved and permits real-time solving of the equations of motion. By carefully selecting the multibody formulation method to be used, it is possible to increase the accuracy of the multibody model while at the same time solving equations of motion in real-time. In this study, a multibody procedure based on semi-recursive and augmented Lagrangian methods for real-time dynamic simulation application is studied in detail. In the semirecursive approach, a velocity transformation matrix is introduced to describe the dependent coordinates into relative (joint) coordinates, which reduces the size of the generalized coordinates. The augmented Lagrangian method is based on usage of global coordinates and, in that method, constraints are accounted using an iterative process. A multibody system can be modelled as either rigid or flexible bodies. When using flexible bodies, the system can be described using a floating frame of reference formulation. In this method, the deformation mode needed can be obtained from the finite element model. As the finite element model typically involves large number of degrees of freedom, reduced number of deformation modes can be obtained by employing model order reduction method such as Guyan reduction, Craig-Bampton method and Krylov subspace as shown in this study The constrained motion of the working mobile vehicles is actuated by the force from the hydraulic actuator. In this study, the hydraulic system is modeled using lumped fluid theory, in which the hydraulic circuit is divided into volumes. In this approach, the pressure wave propagation in the hoses and pipes is neglected. The contact modeling is divided into two stages: contact detection and contact response. Contact detection determines when and where the contact occurs, and contact response provides the force acting at the collision point. The friction between tire and ground is modelled using the LuGre friction model, which describes the frictional force between two surfaces. Typically, the equations of motion are solved in the full matrices format, where the sparsity of the matrices is not considered. Increasing the number of bodies and constraint equations leads to the system matrices becoming large and sparse in structure. To increase the computational efficiency, a technique for solution of sparse matrices is proposed in this dissertation and its implementation demonstrated. To assess the computing efficiency, augmented Lagrangian and semi-recursive methods are implemented employing a sparse matrix technique. From the numerical example, the results show that the proposed approach is applicable and produced appropriate results within the real-time period.

Studies on the Multifaceted Interaction of Atoms and an Electromagnetic Field

In this Thesis the interaction of an electromagnetic field and matter is studied from various aspects in the general framework of cold atoms. Our subjects cover a wide spectrum of phenomena ranging from semiclassical few-level models to fully quantum mechanical interaction with structured reservoirs leading to non-Markovian open quantum system dynamics. Within closed quantum systems, we propose a selective method to manipulate the motional state of atoms in a time-dependent double-well potential and interpret the method in terms of adiabatic processes. Also, we derive a simple wave-packet model, based on distributions of generalized eigenstates, explaining the finite visibility of interference in overlapping continuous-wave atom lasers. In the context of open quantum systems, we develop an unraveling of non-Markovian dynamics in terms of piecewise deterministic quantum jump processes confined in the Hilbert space of the reduced system - the non-Markovian quantum jump method. As examples, we apply it for simple 2- and 3-level systems interacting with a structured reservoir. Also, in the context of ion-cavity QED we study the entanglement generation based on collective Dicke modes in experimentally realistic conditions including photonic losses and an atomic spontaneous decay.

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.

Computational modeling of stented coronary arteries

The application of computational fluid dynamics (CFD) and finite element analysis (FEA) has been growing rapidly in the various fields of science and technology. One of the areas of interest is in biomedical engineering. The altered hemodynamics inside the blood vessels plays a key role in the development of the arterial disease called atherosclerosis, which is the major cause of human death worldwide. Atherosclerosis is often treated with the stenting procedure to restore the normal blood flow. A stent is a tubular, flexible structure, usually made of metals, which is driven and expanded in the blocked arteries. Despite the success rate of the stenting procedure, it is often associated with the restenosis (re-narrowing of the artery) process. The presence of non-biological device in the artery causes inflammation or re-growth of atherosclerotic lesions in the treated vessels. Several factors including the design of stents, type of stent expansion, expansion pressure, morphology and composition of vessel wall influence the restenosis process. Therefore, the role of computational studies is crucial in the investigation and optimisation of the factors that influence post-stenting complications. This thesis focuses on the stent-vessel wall interactions followed by the blood flow in the post-stenting stage of stenosed human coronary artery. Hemodynamic and mechanical stresses were analysed in three separate stent-plaque-artery models. Plaque was modeled as a multi-layer (fibrous cap (FC), necrotic core (NC), and fibrosis (F)) and the arterial wall as a single layer domain. CFD/FEA simulations were performed using commercial software packages in several models mimicking the various stages and morphologies of atherosclerosis. The tissue prolapse (TP) of stented vessel wall, the distribution of von Mises stress (VMS) inside various layers of vessel wall, and the wall shear stress (WSS) along the luminal surface of the deformed vessel wall were measured and evaluated. The results revealed the role of the stenosis size, thickness of each layer of atherosclerotic wall, thickness of stent strut, pressure applied for stenosis expansion, and the flow condition in the distribution of stresses. The thicknesses of FC, and NC and the total thickness of plaque are critical in controlling the stresses inside the tissue. A small change in morphology of artery wall can significantly affect the distribution of stresses. In particular, FC is the most sensitive layer to TP and stresses, which could determine plaque’s vulnerability to rupture. The WSS is highly influenced by the deflection of artery, which in turn is dependent on the structural composition of arterial wall layers. Together with the stenosis size, their roles could play a decisive role in controlling the low values of WSS (<0.5 Pa) prone to restenosis. Moreover, the time dependent flow altered the percentage of luminal area with WSS values less than 0.5 Pa at different time instants. The non- Newtonian viscosity model of the blood properties significantly affects the prediction of WSS magnitude. The outcomes of this investigation will help to better understand the roles of the individual layers of atherosclerotic vessels and their risk to provoke restenosis at the post-stenting stage. As a consequence, the implementation of such an approach to assess the post-stented stresses will assist the engineers and clinicians in optimizing the stenting techniques to minimize the occurrence of restenosis.

Matkatavaroiden käsittelyn asiakaslähtöinen kehittäminen kansainvälisellä lentoasemalla

Lentoasema toimii matkustajien ja matkatavaroiden kauttakulkureittinä vaihdettaessa kulkumuotoa maaliikennevälineen ja lentokoneen tai kahden lentokoneen välillä. Matkustajien kanssa samassa lentokoneessa pitää kulkea myös heidän lentokoneen ruumassa kuljetettavat matkatavaransa. Lentoasemalla suoritettavan matkatavaroiden käsittelyn tehtävänä on huolehtia tavaroiden siirtämisestä matkustajan ja lentokoneen ruuman sekä eri lentokoneiden välillä. Kasvavan kapasiteettitarpeen lisäksi matkatavaroiden käsittelyyn vaikuttavia tulevaisuuden haasteita ovat lentokonetyyppien kehitys, säännösten muutokset, muutokset vastuualueiden jaossa, teknologinen kehitys sekä palvelun laadun merkityksen kasvu. Työssä arvioidaan näiden muutosten merkitystä Helsinki-Vantaan lentoaseman matkatavaroiden käsittelyyn. Tämä on suoritettu esittelemällä käsittelyprosessien nykytilanne ja lentoasemilla tarjottavat palvelut sekä matkatavaroiden käsittelyn vaikutus kokonaispalvelun laatuun. Tämä on mahdollistanut mahdollisten kehityskohteiden löytämisen etenkin palvelun parantamisen kannalta. Vaikka matkatavaroiden käsittelystä vastuullinen osapuoli vaihtuu prosessien aikana monesti, työssä on pyritty keskittymään lentoaseman vastuualueeseen. Tämä tarkoittaa käsittelyssä tarvittavia tiloja sekä tieto- ja mekaanisia järjestelmiä, joiden avulla lentoaseman pitää tarjota kaikille lentoasemalla toimiville yhtiöille tasavertaiset toimintamahdollisuudet. Työssä havaitut merkittävimmät kehityskohteet koskevat tiloja ja saapuvien matkatavaroiden luovuttamista.

Treatment-related renal side-effects in pediatric cancer patients

Background: The long-term side-effects of cancer treatments are of growing importance, since the number of pediatric cancer survivors has considerably increased. Renal side-effects should be noted early to prevent further deterioration. Renal dysfunction may also develop long after cancer treatment. Easy and reliable methods for assessing renal function are needed. Aims: The aims were to find the mechanisms behind methotrexate-induced renal damage by studying renal tubular cells (LLC-PK1cells), and to evaluate the usefulness of laboratory tests in assessing glomerular function in pediatric cancer patients by comparing an isotope clearance method with alternative methods. The aim was also to study the long-term effects of bone marrow transplantation (BMT) and high-dose methotrexate (HD-MTX) treatment in renal function. Results: Methotrexate induced time-dependent renal tubular cell swelling and cell death. In patients treated with HD-MTX a significant decrease in GFR was noted after a follow-up time of one to ten years. One year after BMTthe GFR was reduced, especially in patients treated with total body irradiation (TBI). GFR recovered slightly but remained stable thereafter. In glomerular function assessment the serum cystatin C (cysC) concentration showed a significant association with GFR measured by the isotope method. Conclusions: Methotrexate induced acute damage in renal tubular cells. In assessing GFR the isotope method still remains the method of choice, but the assay of cystatin C was the most reliable of other alternatives. Long-term follow-up of renal function is needed in BMT patients and patients treated with HD-MTX.

Hankintaprosessin tehostaminen tietojärjestelmien avulla

Tämän tutkielman tavoitteena oli selvittää miten hankintaprosessia voidaan tehostaa tietojärjestelmien avulla. Toiminnan tehostaminen vaatii yleensä tietojärjestelmien uusimista, joten tavoitteena oli myös selvittää mitä seikkoja pitää ottaa huomioon tietojärjestelmähankintaa suunniteltaessa ja tehtäessä vaatimusmäärittelyä. Tutkielma toteutettiin kvalitatiivisena tapaustutkimuksena, jossa aineisto kerättiin haastatteluiden ja havainnoinnin avulla. Tapaukseksi valittiin kohdeyrityksen osto-osasto ja sen hankintaprosessi. Kohdeyrityksessä tietojärjestelmien uusimisen avulla pystyttäisiin tehostamaan etenkin tiedon hallintaa. Nykyisissä järjestelmissä tiedot ovat hajallaan eri paikoissa ja tarvittavien tietojen hakuun kuluu ylimääräistä aikaa. Kohdeyrityksen hankintaprosessia tehostavia tietojärjestelmiä ovat ennustamiseen, raportointiin, tilausperusteiseen ostotyöhön, laskutukseen ja painatusmateriaalien sähköiseen hallintaan liittyvät järjestelmät. Myös paremmin hankintatoimen tarpeet huomioiva toiminnanohjausjärjestelmä tehostaisi hankintaprosessia.

Pulsatile Blood Flow Simulations in Computed Tomography(CT) Scan-Based and Idealized Geometries of Human Aorta

Blood flow in human aorta is an unsteady and complex phenomenon. The complex patterns are related to the geometrical features like curvature, bends, and branching and pulsatile nature of flow from left ventricle of heart. The aim of this work was to understand the effect of aorta geometry on the flow dynamics. To achieve this, 3D realistic and idealized models of descending aorta were reconstructed from Computed Tomography (CT) images of a female patient. The geometries were reconstructed using medical image processing code. The blood flow in aorta was assumed to be laminar and incompressible and the blood was assumed to be Newtonian fluid. A time dependent pulsatile and parabolic boundary condition was deployed at inlet. Steady and unsteady blood flow simulations were performed in real and idealized geometries of descending aorta using a Finite Volume Method (FVM) code. Analysis of Wall Shear Stress (WSS) distribution, pressure distribution, and axial velocity profiles were carried out in both geometries at steady and unsteady state conditions. The results obtained in thesis work reveal that the idealization of geometry underestimates the values of WSS especially near the region with sudden change of diameter. However, the resultant pressure and velocity in idealized geometry are close to those in real geometry

Numerical Investigation of Flow Past a Circular Cylinder and in a Staggered Tube Bundle Using Various Turbulence Models

Transitional flow past a three-dimensional circular cylinder is a widely studied phenomenon since this problem is of interest with respect to many technical applications. In the present work, the numerical simulation of flow past a circular cylinder, performed by using a commercial CFD code (ANSYS Fluent 12.1) with large eddy simulation (LES) and RANS (κ - ε and Shear-Stress Transport (SST) κ - ω! model) approaches. The turbulent flow for ReD = 1000 & 3900 is simulated to investigate the force coefficient, Strouhal number, flow separation angle, pressure distribution on cylinder and the complex three dimensional vortex shedding of the cylinder wake region. The numerical results extracted from these simulations have good agreement with the experimental data (Zdravkovich, 1997). Moreover, grid refinement and time-step influence have been examined. Numerical calculations of turbulent cross-flow in a staggered tube bundle continues to attract interest due to its importance in the engineering application as well as the fact that this complex flow represents a challenging problem for CFD. In the present work a time dependent simulation using κ – ε, κ - ω! and SST models are performed in two dimensional for a subcritical flow through a staggered tube bundle. The predicted turbulence statistics (mean and r.m.s velocities) have good agreement with the experimental data (S. Balabani, 1996). Turbulent quantities such as turbulent kinetic energy and dissipation rate are predicted using RANS models and compared with each other. The sensitivity of grid and time-step size have been analyzed. Model constants sensitivity study have been carried out by adopting κ – ε model. It has been observed that model constants are very sensitive to turbulence statistics and turbulent quantities.

A study of combustion phenomena in circulating fluidized beds by developing and applying experimental and modelling methods for laboratory-scale reactors

The results shown in this thesis are based on selected publications of the 2000s decade. The work was carried out in several national and EC funded public research projects and in close cooperation with industrial partners. The main objective of the thesis was to study and quantify the most important phenomena of circulating fluidized bed combustors by developing and applying proper experimental and modelling methods using laboratory scale equipments. An understanding of the phenomena plays an essential role in the development of combustion and emission performance, and the availability and controls of CFB boilers. Experimental procedures to study fuel combustion behaviour under CFB conditions are presented in the thesis. Steady state and dynamic measurements under well controlled conditions were carried out to produce the data needed for the development of high efficiency, utility scale CFB technology. The importance of combustion control and furnace dynamics is emphasized when CFB boilers are scaled up with a once through steam cycle. Qualitative information on fuel combustion characteristics was obtained directly by comparing flue gas oxygen responses during the impulse change experiments with fuel feed. A one-dimensional, time dependent model was developed to analyse the measurement data Emission formation was studied combined with fuel combustion behaviour. Correlations were developed for NO, N2O, CO and char loading, as a function of temperature and oxygen concentration in the bed area. An online method to characterize char loading under CFB conditions was developed and validated with the pilot scale CFB tests. Finally, a new method to control air and fuel feeds in CFB combustion was introduced. The method is based on models and an analysis of the fluctuation of the flue gas oxygen concentration. The effect of high oxygen concentrations on fuel combustion behaviour was also studied to evaluate the potential of CFB boilers to apply oxygenfiring technology to CCS. In future studies, it will be necessary to go through the whole scale up chain from laboratory phenomena devices through pilot scale test rigs to large scale, commercial boilers in order to validate the applicability and scalability of the, results. This thesis shows the chain between the laboratory scale phenomena test rig (bench scale) and the CFB process test rig (pilot). CFB technology has been scaled up successfully from an industrial scale to a utility scale during the last decade. The work shown in the thesis, for its part, has supported the development by producing new detailed information on combustion under CFB conditions.

Ion selective electrodes for determination of low and ultra low concentrations of lead (II) in natural waters

Att övervaka förekomsten av giftiga komponenter i naturliga vattendrag är nödvändigt för människans välmående. Eftersom halten av föroreningar i naturens ekosystem bör hållas möjligast låg, pågår en ständig jakt efter kemiska analysmetoder med allt lägre detektionsgränser. I dagens läge görs miljöanalyser med dyr och sofistikerad instrumentering som kräver mycket underhåll. Jonselektiva elektroder har flera goda egenskaper som t.ex. bärbarhet, låg energiförbrukning, och dessutom är de relativt kostnadseffektiva. Att använda jonselektiva elektroder vid miljöanalyser är möjligt om deras känslighetsområde kan utvidgas genom att sänka deras detektionsgränser. För att sänka detektionsgränsen för Pb(II)-selektiva elektroder undersöktes olika typer av jonselektiva membran som baserades på polyakrylat-kopolymerer, PVC och PbS/Ag2S. Fast-fas elektroder med membran av PbS/Ag2S är i allmänhet enklare och mer robusta än konventionella elektroder vid spårämnesanalys av joniska föroreningar. Fast-fas elektrodernas detektionsgräns sänktes i detta arbete med en nyutvecklad galvanostatisk polariseringsmetod och de kunde sedan framgångsrikt användas för kvantitativa bestämningar av bly(II)-halter i miljöprov som hade samlats in i den finska skärgården nära tidigare industriområden. Analysresultaten som erhölls med jonselektiva elektroder bekräftades med andra analytiska metoder. Att sänka detektionsgränsen m.hj.a. den nyutvecklade polariseringsmetoden möjliggör bestämning av låga och ultra-låga blyhalter som inte kunde nås med klassisk potentiometri. Den verkliga fördelen med att använda dessa blyselektiva elektroder är möjligheten att utföra mätningar i obehandlade miljöprov trots närvaron av fasta partiklar vilket inte är möjligt att göra med andra analysmetoder. Jag väntar mig att den nyutvecklade polariseringsmetoden kommer att sätta en trend i spårämnesanalys med jonselektiva elektroder.

Education export as a form of knowledge transfer : The potential of Finnish education export in Chile

Knowledge transfer is a complex process. Knowledge transfer in the form of exporting education products from one system of education to another is particularly complicated, because each system has been developed in a particular context to meet the requirements seen as relevant at each time. National innovation systems are often seen to form an essential framework within which the development of a country, its economy and level of knowledge are considered and promoted. These systems are orientated towards the future, and as such they also provide a framework for the knowledge transfer related to the development of education. In the best of circumstances they are able to facilitate and boost this transfer both from the viewpoint of the provider and the recipient. The leading thought and the idea of the study is that education export is a form of knowledge transfer, which is illustrated by the existing models included. The purpose of this study is to explore, analyze and describe the factors and phenomena related to education export, and more specifically, those related to the experiences and potential of Finnish education export to Chile. For better understanding, of the multiplicity of the issue involved, the current status of education export between Finland and Chile and he existing efforts within the Finnish innovation network will be outlined as well as new forms of co-operation between Finland and Chile in educational matters explored. Several countries have started to commercialize their education system in order to establish themselves as emerging education exporters. Moreover, the demand for education reform is accurate in many developing countries. This offers a good match between Finland and Chile to be the example countries of the research. The main research findings suggest that there are several business areas in education export. These include degrees in education, training services and education technologies for example The factors that influence education export can be divided into four groups, including academic, cultural, political and economic aspects. Challenges to overcome include the lack of product or services to be sold, lack of market and cultural knowledge of the buyer country, financing and lack of suitable pricing model. National innovation systems could be seen as enabling entities for successful education export. The extensive networks that national innovation systems aim to form, could operate as a basis for joining the forces in selling knowledge as well as receiving knowledge in a constructive way.

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.

Modeling of fireside deposit formation in two industrial furnaces

Fireside deposits can be found in many types of utility and industrial furnaces. The deposits in furnaces are problematic because they can reduce heat transfer, block gas paths and cause corrosion. To tackle these problems, it is vital to estimate the influence of deposits on heat transfer, to minimize deposit formation and to optimize deposit removal. It is beneficial to have a good understanding of the mechanisms of fireside deposit formation. Numerical modeling is a powerful tool for investigating the heat transfer in furnaces, and it can provide valuable information for understanding the mechanisms of deposit formation. In addition, a sub-model of deposit formation is generally an essential part of a comprehensive furnace model. This work investigates two specific processes of fireside deposit formation in two industrial furnaces. The first process is the slagging wall found in furnaces with molten deposits running on the wall. A slagging wall model is developed to take into account the two-layer structure of the deposits. With the slagging wall model, the thickness and the surface temperature of the molten deposit layer can be calculated. The slagging wall model is used to predict the surface temperature and the heat transfer to a specific section of a super-heater tube panel with the boundary condition obtained from a Kraft recovery furnace model. The slagging wall model is also incorporated into the computational fluid dynamics (CFD)-based Kraft recovery furnace model and applied on the lower furnace walls. The implementation of the slagging wall model includes a grid simplification scheme. The wall surface temperature calculated with the slagging wall model is used as the heat transfer boundary condition. Simulation of a Kraft recovery furnace is performed, and it is compared with two other cases and measurements. In the two other cases, a uniform wall surface temperature and a wall surface temperature calculated with a char bed burning model are used as the heat transfer boundary conditions. In this particular furnace, the wall surface temperatures from the three cases are similar and are in the correct range of the measurements. Nevertheless, the wall surface temperature profiles with the slagging wall model and the char bed burning model are different because the deposits are represented differently in the two models. In addition, the slagging wall model is proven to be computationally efficient. The second process is deposit formation due to thermophoresis of fine particles to the heat transfer surface. This process is considered in the simulation of a heat recovery boiler of the flash smelting process. In order to determine if the small dust particles stay on the wall, a criterion based on the analysis of forces acting on the particle is applied. Time-dependent simulation of deposit formation in the heat recovery boiler is carried out and the influence of deposits on heat transfer is investigated. The locations prone to deposit formation are also identified in the heat recovery boiler. Modeling of the two processes in the two industrial furnaces enhances the overall understanding of the processes. The sub-models developed in this work can be applied in other similar deposit formation processes with carefully-defined boundary conditions.