Sellutehtaan valkaisusäätöjen evaluointi

Tämän diplomityön tarkoituksena oli arvioida UPM-Kymmene Oyj, Kaukaan sellutehtaan havukuitulinjan uuden valkaisusäätömenetelmän, kompensoidun kappakerroinsäädön, ja vaaleusanalysaattoreiden käyttöönoton vaikutukset klooridioksidikulutukseen ja sellun laatuun. Kirjallisuusosassa esitelläänvalkaisun tärkeimpiä on-line ja in-line mittalaitteita sekä niitä hyödyntäviä säätöratkaisuja. Tiettävästi paras saatavilla oleva säätömenetelmä, kompensoitu kappakerroinsäätö määrittää klooridioksidiannoksen vaiheeseen tulevan massan kappaluvun tai vaaleuden perusteella. Valkaisukemikaaliannos optimoidaan jatkuvatoimisten vaaleus- ja jäännösmittausten sekä vaiheen jälkeisen saavutetun kappaluvuntai vaaleuden takaisinkytkennällä. Kokeellisessa osassa arvioitiin havulinjan valkaisusäätömenetelmän käyttöönoton vaikutukset klooridioksidikulutuksiin ja sellun laatuun vertailemalla uuden ja vanhan säätömenetelmän aikaisia kemikaalikulutuksia sekä välikappaluvun ja loppuvaaleuden itseisarvopoikkeamia tavoitteistaan. Lisäksi selvitettiin annossäätöön ja kompensoituun vaaleussäätöön perustuvan vanhan sekä kompensoitua kappakerroinsäätöä soveltavan uuden säätömenetelmän toimintaperiaatteet. Työssä kartoitettiin myös valkaisulinjan tarpeet uusille mittalaitteille sekä analysaattorin mittausten vertailukelpoisuus laboratorioanalyysien kanssa. Edelleen työssä arvioitiin tehtaan X koivuvalkaisulinjan uuden säätömenetelmän toimivuus. Kompensoidun kappakerroinsäädön käyttöönoton myötä valkaisun klooridioksidikulutukset pienenivät sekä välikappaluku ja loppuvaaleus saatiin pidettyä lähempänä tavoitteittaan. Säätömenetelmän toimivuuttaedisti vaaleusanalysaattoreiden mittausten korreloivuus laboratorioanalyysien kanssa. Vastaavat vaikutukset kemikaalikulutuksiin ja sellun laatuun havaittiin myös tehtaan X valkaisusäädön käyttöönotosta.

Kosketukseton hitsin pintageometrian määrittäminen

Yleisesti tiedetään hitsin pintageometrian vaikuttavan rakenteen väsymislujuuteen. Nopean, edullisen ja luotettavan pintageometrian mittausmenetelmän kehittäminen on askel kohti tarkempaa ja varmempaa rakenteen väsymislujuuden tarkastelua. Tässä työssä on tutkittu hitsejä, joiden pinnan geometria on mitattu norjalaisen SINTEF -yrityksen kehittämällä rakenteellisen valon menetelmällä. Osana työtä kehitettiin MatLab -pohjainen ohjelma, jolla jälkikäsitellään mittauksesta saadut x-y-z -mittapisteet. Mittausdatan jälkikäsittelyssä saadaan mittauksesta määritettyähitsin reunan pyöristys, liittymäkulma, a-mitta, reunahaava ja kateettisuhde. Kehitettyä menetelmää käyttämällä mitattiin lähes 300 voimaakantamatontaristiliitoksen hitsiä. Mittaustuloksia verrattiin vastaavista kappaleista tehtyihin hiemittauksiin. Manuaalisen hieestä tehdyn mittauksen havaittiin olevan tarkempi ja pystyttiin havaitsemaan paikallisempia muotoja. Rakenteellisen valon mittauksissa tapahtunut heijastelu saatiin pienenemään käsittelemällä mitattava pinta mattavalkoisella maalilla. Rakenteellisen valon mittatarkkuudeksi saatiin noin 0,2 mm. Pohjautuen mitattuun hitsin reunan pyöristykseen ja liittymäkulmaan voidaan yksinkertaista kaavaa käyttämällä laskea hitsin jännityskonsentraatio ja näin saada alkuarvaus väsymislujuudelle. Myös muiden tekijöiden tiedetään vaikuttavan hitsin väsymislujuuteen, joten pyöristyksen ja liittymäkulman avulla tehdyt arviot eivät ole absoluuttisen oikeita. Tämä havaittiin väsytyskokeilla, joista yhdessä väsymisvaurio ei syntynyt suurimmankaan jännityskonsentraation alueella.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Direct torque control of permanent magnet synchronous machines - analysis and implementation

Thedirect torque control (DTC) has become an accepted vector control method besidethe current vector control. The DTC was first applied to asynchronous machines,and has later been applied also to synchronous machines. This thesis analyses the application of the DTC to permanent magnet synchronous machines (PMSM). In order to take the full advantage of the DTC, the PMSM has to be properly dimensioned. Therefore the effect of the motor parameters is analysed taking the control principle into account. Based on the analysis, a parameter selection procedure is presented. The analysis and the selection procedure utilize nonlinear optimization methods. The key element of a direct torque controlled drive is the estimation of the stator flux linkage. Different estimation methods - a combination of current and voltage models and improved integration methods - are analysed. The effect of an incorrect measured rotor angle in the current model is analysed andan error detection and compensation method is presented. The dynamic performance of an earlier presented sensorless flux estimation method is made better by improving the dynamic performance of the low-pass filter used and by adapting the correction of the flux linkage to torque changes. A method for the estimation ofthe initial angle of the rotor is presented. The method is based on measuring the inductance of the machine in several directions and fitting the measurements into a model. The model is nonlinear with respect to the rotor angle and therefore a nonlinear least squares optimization method is needed in the procedure. A commonly used current vector control scheme is the minimum current control. In the DTC the stator flux linkage reference is usually kept constant. Achieving the minimum current requires the control of the reference. An on-line method to perform the minimization of the current by controlling the stator flux linkage reference is presented. Also, the control of the reference above the base speed is considered. A new estimation flux linkage is introduced for the estimation of the parameters of the machine model. In order to utilize the flux linkage estimates in off-line parameter estimation, the integration methods are improved. An adaptive correction is used in the same way as in the estimation of the controller stator flux linkage. The presented parameter estimation methods are then used in aself-commissioning scheme. The proposed methods are tested with a laboratory drive, which consists of a commercial inverter hardware with a modified software and several prototype PMSMs.

Dynamic analysis of belt-drives using the absolute nodal coordinate formulation

Belt-drive systems have been and still are the most commonly used power transmission form in various applications of different scale and use. The peculiar features of the dynamics of the belt-drives include highly nonlinear deformation,large rigid body motion, a dynamical contact through a dry friction interface between the belt and pulleys with sticking and slipping zones, cyclic tension of the belt during the operation and creeping of the belt against the pulleys. The life of the belt-drive is critically related on these features, and therefore, amodel which can be used to study the correlations between the initial values and the responses of the belt-drives is a valuable source of information for the development process of the belt-drives. Traditionally, the finite element models of the belt-drives consist of a large number of elements thatmay lead to computational inefficiency. In this research, the beneficial features of the absolute nodal coordinate formulation are utilized in the modeling of the belt-drives in order to fulfill the following requirements for the successful and efficient analysis of the belt-drive systems: the exact modeling of the rigid body inertia during an arbitrary rigid body motion, the consideration of theeffect of the shear deformation, the exact description of the highly nonlinear deformations and a simple and realistic description of the contact. The use of distributed contact forces and high order beam and plate elements based on the absolute nodal coordinate formulation are applied to the modeling of the belt-drives in two- and three-dimensional cases. According to the numerical results, a realistic behavior of the belt-drives can be obtained with a significantly smaller number of elements and degrees of freedom in comparison to the previously published finite element models of belt-drives. The results of theexamples demonstrate the functionality and suitability of the absolute nodal coordinate formulation for the computationally efficient and realistic modeling ofbelt-drives. This study also introduces an approach to avoid the problems related to the use of the continuum mechanics approach in the definition of elastic forces on the absolute nodal coordinate formulation. This approach is applied to a new computationally efficient two-dimensional shear deformable beam element based on the absolute nodal coordinate formulation. The proposed beam element uses a linear displacement field neglecting higher-order terms and a reduced number of nodal coordinates, which leads to fewer degrees of freedom in a finite element.

Paperikoneen laatusäätöjärjestelmän uusinnan vaikutukset laatusäätöjen suorituskykyyn

Tutkimuksen kohteena olleen UPM-Kymmene Oyj Kajaanin tehtaan PK3:n laatusäätöjärjestelmä ja mittapalkki uusittiin, jolloin haluttiin selvittää uusinnan vaikutuksia laatusäätöjen suorituskykyyn ja paperin laatuun. Työn kirjallisessa osassa perehdyttiin paperinvalmistusprosessin osiin kyseisen sanomalehtipaperikoneen tapauksessa sekä keskeisimpiin paperin laatuominaisuuksiin liittyviin mittaus- ja säätölaitteisiin sekä niiden toimintaan. Seurattaviksi paperin laatusuureiksi valittiin neliömassa, kuivamassa, kosteus ja paksuus, jotka ovat sanomalehtipaperin tärkeimpiä online-mitattavia ominaisuuksia. Paperin laatusuureiden seurantaan käytetään erilaisia tunnuslukuja ja työkaluja, joita on esitelty tässä työssä. Laatusuureiden konesuuntaisen ja poikkisuuntaisen seurannan tunnusluvuksi valittiin yleisesti käytössä oleva 2σ-keskiarvohajonta. Säätöjen suorituskykyä seurattiin suorituskykykolmion ohjausmatkaindeksien (CTI) ja erosuureen integraalien (IAE) avulla. Kokeellisessa osassa kerättiin mittaustietoja sekä vanhan että uuden laatusäätöjärjestelmän aikana. Seurattavat ajotilanteet paperikoneella jaettiin stabiiliin ajoon ja muutostilanteisiin, jotka käsittävät katkot ja lajinvaihtotilanteet. Stabiilin ajon aikana selvitettiin laatusuureiden hajontojen ja säätöjen suorituskykyindeksien normaaleissa tasoissa tapahtuneet muutokset. Muutostilanteiden osalta haluttiin selvittää, nopeuttaako järjestelmäuusinta katkoista toipumista ja lajinvaihtoaikaa. Stabiilin ajon seurannasta saatujen tulosten perusteella neliömassan ja kuivamassan konesuuntaiset hajonnat kasvoivat järjestelmäuusinnan myötä, mutta kosteuden konesuuntaiset hajonnat pienenivät. Laatusuureiden poikkisuuntaisista hajonnoista neliömassan sekä kuivamassan hajonnat kasvoivat ja kosteuden sekä paksuuden hajonnat pienenivät joidenkin lajien osalta. Poikkisuuntaisten laatusuureiden, etenkin paksuuden, toipuminen katkon jälkeen nopeutui. Myös lajinvaihtoon kuluva aika lyheni poikkisuuntaisilla laatusuureilla. Muutostilanteiden konesuuntaisten hajontojen asettumisajat eivät juuri parantuneet.

Development of beam and plate finite elements based on the absolute nodal coordinate formulation

The focus of this dissertation is to develop finite elements based on the absolute nodal coordinate formulation. The absolute nodal coordinate formulation is a nonlinear finite element formulation, which is introduced for special requirements in the field of flexible multibody dynamics. In this formulation, a special definition for the rotation of elements is employed to ensure the formulation will not suffer from singularities due to large rotations. The absolute nodal coordinate formulation can be used for analyzing the dynamics of beam, plate and shell type structures. The improvements of the formulation are mainly concentrated towards the description of transverse shear deformation. Additionally, the formulation is verified by using conventional iso-parametric solid finite element and geometrically exact beam theory. Previous claims about especially high eigenfrequencies are studied by introducing beam elements based on the absolute nodal coordinate formulation in the framework of the large rotation vector approach. Additionally, the same high eigenfrequency problem is studied by using constraints for transverse deformation. It was determined that the improvements for shear deformation in the transverse direction lead to clear improvements in computational efficiency. This was especially true when comparative stress must be defined, for example when using elasto-plastic material. Furthermore, the developed plate element can be used to avoid certain numerical problems, such as shear and curvature lockings. In addition, it was shown that when compared to conventional solid elements, or elements based on nonlinear beam theory, elements based on the absolute nodal coordinate formulation do not lead to an especially stiff system for the equations of motion.

MCMC Analysis Of Classical Time Series Algorithms

Identification of order of an Autoregressive Moving Average Model (ARMA) by the usual graphical method is subjective. Hence, there is a need of developing a technique to identify the order without employing the graphical investigation of series autocorrelations. To avoid subjectivity, this thesis focuses on determining the order of the Autoregressive Moving Average Model using Reversible Jump Markov Chain Monte Carlo (RJMCMC). The RJMCMC selects the model from a set of the models suggested by better fitting, standard deviation errors and the frequency of accepted data. Together with deep analysis of the classical Box-Jenkins modeling methodology the integration with MCMC algorithms has been focused through parameter estimation and model fitting of ARMA models. This helps to verify how well the MCMC algorithms can treat the ARMA models, by comparing the results with graphical method. It has been seen that the MCMC produced better results than the classical time series approach.

Sensitive Troponin Assays: Diagnostic and Prognostic Use in Cardiology

Cardiac troponins (cTns) are the recommended biochemical markers in the diagnosis of myocardial infarction (MI). They are very sensitive and tissue-specific but are limited by their delayed appearance in the circulation. Biochemical markers with more rapid release kinetics, e.g. myoglobin and especially heart-type fatty acid-binding protein (H-FABP), have been used to enhance the early identification of MI. The implementation of cTns into clinical practice has shown that cardiomyocyte injury occurs in many other clinical conditions than MI. The aim of this study was to evaluate the impact of modern and highly sensitive cTnI assays on the early diagnosis of MI. In a patient cohort with suspected MI, such a sensitive cTnI assay enhanced the early diagnostic accuracy when compared to a less sensitive cTnI assay and to myoglobin. When compared to H-FABP during the early hours after symptom onset, the sensitive cTnI assay showed at least similar and, after 6 hours, superior diagnostic accuracy. A positive cTnI test result had superior prognostic value when compared to H-FABP, even among early presenters. The prognostic value of cTn in acute heart failure (AHF) was evaluated in 364 patients who participated in the FINN-AKVA study. The patients presented with AHF but no acute coronary syndrome (ACS). Up to half of the patients had elevated cTn levels which were associated with higher 6-month mortality. The magnitude of cTn elevation was directly proportional to mortality. Finally, the clinical spectrum of cTnI elevations was evaluated in 991 cTnI positive emergency department (ED) patients. 83% of the patients had MI and 17% had cTnI elevation due to other clinical conditions. The latter patient group was characterized by lower absolute cTnI levels and – importantly – higher in-hospital mortality when compared to the MI patients. In conclusion, the use of a highly sensitive cTnI assay enhances the early diagnostic accuracy and risk stratification in suspected MI patients. Cardiac troponin elevations are highly prevalent also in other acute clinical conditions and indicate an adverse outcome of these patients.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

A chart of the southern part of South America : with the track of the centurion from the island of St. Catherines to the island of Juan Fernandes in which is inserted the variation and soundings observ'd on board her, together with her deviation from her estimated course in passing round Cap Horn, occasiond by the force of the currents / R. W. Seale sculp

Nimeketiedot nimiönkehyksissä

Together Everyone Achieves More? Herding Behavior in the Finnish Stock Market

This thesis examines whether or not Finnish stock markets has herding behavior. Sample data is from 2004 to 2013. Including total of 2516 market days. Market wide herding, up and down market herding, extreme price movement herding and turnover volume herding are measured in this thesis. Methods used in this thesis are cross-sectional absolute dispersion and cross-sectional standard deviation. This thesis found no signs of herding in the Finnish stock market.

Partikkelikokojakaumamenetelmän testaus INCA Feature

Tässä työssä testattiin partikkelikokojakaumien analysoinnissa käytettävää kuvankäsittelyohjelmaa INCA Feature. Partikkelikokojakaumat määritettiin elektronimikroskooppikuvista INCA Feature ohjelmaa käyttäen partikkeleiden projektiokuvista päällystyspigmenttinä käytettävälle talkille ja kahdelle eri karbonaattilaadulle. Lisäksi määritettiin partikkelikokojakaumat suodatuksessa ja puhdistuksessa apuaineina käytettäville piidioksidi- ja alumiinioksidihiukkasille. Kuvankäsittelyohjelmalla määritettyjä partikkelikokojakaumia verrattiin partikkelin laskeutumisnopeuteen eli sedimentaatioon perustuvalla SediGraph 5100 analysaattorilla ja laserdiffraktioon perustuvalla Coulter LS 230 menetelmällä analysoituihin partikkelikokojakaumiin. SediGraph 5100 ja kuva-analyysiohjelma antoivat talkkipartikkelien kokojakaumalle hyvin samankaltaisen keskiarvon. Sen sijaan Coulter LS 230 laitteen antama kokojakauman keskiarvo poikkesi edellisistä. Kaikki vertailussa olleet partikkelikokojakaumamenetelmät asettivat eri näytteiden partikkelit samaan kokojärjestykseen. Kuitenkaan menetelmien tuloksia ei voida numeerisesti verrata toisiinsa, sillä kaikissa käytetyissä analyysimenetelmissä partikkelikoon mittaus perustuu partikkelin eri ominaisuuteen. Työn perusteella kaikki testatut analyysimenetelmät soveltuvat paperipigmenttien partikkelikokojakaumien määrittämiseen. Tässä työssä selvitettiin myös kuva-analyysiin tarvittava partikkelien lukumäärä, jolla analyysitulos on luotettava. Työssä todettiin, että analysoitavien partikkelien lukumäärän tulee olla vähintään 300 partikkelia. Liian suuri näytemäärä lisää kokojakauman hajontaa ja pidentää analyysiin käytettyä aikaa useaan tuntiin. Näytteenkäsittely vaatii vielä lisää tutkimuksia, sillä se on tärkein ja kriittisin vaihe SEM ja kuva-analyysiohjelmalla tehtävää partikkelikokoanalyysiä. Automaattisten mikroskooppien yleistyminen helpottaa ja nopeuttaa analyysien tekoa, jolloin menetelmän suosio tulee kasvamaan myös paperipigmenttien tutkimuksessa. Laitteiden korkea hinta ja käyttäjältä vaadittava eritysosaaminen tulevat rajaamaan käytön ainakin toistaiseksi tutkimuslaitoksiin.

Heat-transfer improvements in an axial-flux permanent-magnet synchronous machine

Axial-flux machines tend to have cooling difficulties since it is difficult to arrange continuous heat path between the stator stack and the frame. One important reason for this is that no shrink fitting of the stator is possible in an axial-flux machine. Using of liquid-cooled end shields does not alone solve this issue. Cooling of the rotor and the end windings may also be difficult at least in case of two-stator-single-rotor construction where air circulation in the rotor and in the end-winding areas may be difficult to arrange. If the rotor has significant losses air circulation via the rotor and behind the stator yokes should be arranged which, again, weakens the stator cooling. In this paper we study a novel way of using copper bars as extra heat transfer paths between the stator teeth and liquid cooling pools in the end shields. After this the end windings still suffer of low thermal conductivity and means for improving this by high-heat-conductance material was also studied. The design principle of each cooling system is presented in details. Thermal models based on Computational Fluid Dynamics (CFD) are used to analyse the temperature distribution in the machine. Measurement results are provided from different versions of the machine. The results show that significant improvements in the cooling can be gained by these steps.