Selective catalytic reduction of nitrogen oxides with ammonia in forced unsteady state reactors - Case based reasoning and mathematical model simulation reasoning

The application of forced unsteady-state reactors in case of selective catalytic reduction of nitrogen oxides (NOx) with ammonia (NH3) is sustained by the fact that favorable temperature and composition distributions which cannot be achieved in any steady-state regime can be obtained by means of unsteady-state operations. In a normal way of operation the low exothermicity of the selective catalytic reduction (SCR) reaction (usually carried out in the range of 280-350°C) is not enough to maintain by itself the chemical reaction. A normal mode of operation usually requires supply of supplementary heat increasing in this way the overall process operation cost. Through forced unsteady-state operation, the main advantage that can be obtained when exothermic reactions take place is the possibility of trapping, beside the ammonia, the moving heat wave inside the catalytic bed. The unsteady state-operation enables the exploitation of the thermal storage capacity of the catalyticbed. The catalytic bed acts as a regenerative heat exchanger allowing auto-thermal behaviour when the adiabatic temperature rise is low. Finding the optimum reactor configuration, employing the most suitable operation model and identifying the reactor behavior are highly important steps in order to configure a proper device for industrial applications. The Reverse Flow Reactor (RFR) - a forced unsteady state reactor - corresponds to the above mentioned characteristics and may be employed as an efficient device for the treatment of dilute pollutant mixtures. As a main disadvantage, beside its advantages, the RFR presents the 'wash out' phenomena. This phenomenon represents emissions of unconverted reactants at every switch of the flow direction. As a consequence our attention was focused on finding an alternative reactor configuration for RFR which is not affected by the incontrollable emissions of unconverted reactants. In this respect the Reactor Network (RN) was investigated. Its configuration consists of several reactors connected in a closed sequence, simulating a moving bed by changing the reactants feeding position. In the RN the flow direction is maintained in the same way ensuring uniformcatalyst exploitation and in the same time the 'wash out' phenomena is annulated. The simulated moving bed (SMB) can operate in transient mode giving practically constant exit concentration and high conversion levels. The main advantage of the reactor network operation is emphasizedby the possibility to obtain auto-thermal behavior with nearly uniformcatalyst utilization. However, the reactor network presents only a small range of switching times which allow to reach and to maintain an ignited state. Even so a proper study of the complex behavior of the RN may give the necessary information to overcome all the difficulties that can appear in the RN operation. The unsteady-state reactors complexity arises from the fact that these reactor types are characterized by short contact times and complex interaction between heat and mass transportphenomena. Such complex interactions can give rise to a remarkable complex dynamic behavior characterized by a set of spatial-temporal patterns, chaotic changes in concentration and traveling waves of heat or chemical reactivity. The main efforts of the current research studies concern the improvement of contact modalities between reactants, the possibility of thermal wave storage inside the reactor and the improvement of the kinetic activity of the catalyst used. Paying attention to the above mentioned aspects is important when higher activity even at low feeding temperatures and low emissions of unconverted reactants are the main operation concerns. Also, the prediction of the reactor pseudo or steady-state performance (regarding the conversion, selectivity and thermal behavior) and the dynamicreactor response during exploitation are important aspects in finding the optimal control strategy for the forced unsteady state catalytic tubular reactors. The design of an adapted reactor requires knowledge about the influence of its operating conditions on the overall process performance and a precise evaluation of the operating parameters rage for which a sustained dynamic behavior is obtained. An apriori estimation of the system parameters result in diminution of the computational efforts. Usually the convergence of unsteady state reactor systems requires integration over hundreds of cycles depending on the initial guess of the parameter values. The investigation of various operation models and thermal transfer strategies give reliable means to obtain recuperative and regenerative devices which are capable to maintain an auto-thermal behavior in case of low exothermic reactions. In the present research work a gradual analysis of the SCR of NOx with ammonia process in forced unsteady-state reactors was realized. The investigation covers the presentationof the general problematic related to the effect of noxious emissions in the environment, the analysis of the suitable catalysts types for the process, the mathematical analysis approach for modeling and finding the system solutions and the experimental investigation of the device found to be more suitable for the present process. In order to gain information about the forced unsteady state reactor design, operation, important system parameters and their values, mathematical description, mathematicalmethod for solving systems of partial differential equations and other specific aspects, in a fast and easy way, and a case based reasoning (CBR) approach has been used. This approach, using the experience of past similarproblems and their adapted solutions, may provide a method for gaining informations and solutions for new problems related to the forced unsteady state reactors technology. As a consequence a CBR system was implemented and a corresponding tool was developed. Further on, grooving up the hypothesis of isothermal operation, the investigation by means of numerical simulation of the feasibility of the SCR of NOx with ammonia in the RFRand in the RN with variable feeding position was realized. The hypothesis of non-isothermal operation was taken into account because in our opinion ifa commercial catalyst is considered, is not possible to modify the chemical activity and its adsorptive capacity to improve the operation butis possible to change the operation regime. In order to identify the most suitable device for the unsteady state reduction of NOx with ammonia, considering the perspective of recuperative and regenerative devices, a comparative analysis of the above mentioned two devices performance was realized. The assumption of isothermal conditions in the beginningof the forced unsteadystate investigation allowed the simplification of the analysis enabling to focus on the impact of the conditions and mode of operation on the dynamic features caused by the trapping of one reactant in the reactor, without considering the impact of thermal effect on overall reactor performance. The non-isothermal system approach has been investigated in order to point out the important influence of the thermal effect on overall reactor performance, studying the possibility of RFR and RN utilization as recuperative and regenerative devices and the possibility of achieving a sustained auto-thermal behavior in case of lowexothermic reaction of SCR of NOx with ammonia and low temperature gasfeeding. Beside the influence of the thermal effect, the influence of the principal operating parameters, as switching time, inlet flow rate and initial catalyst temperature have been stressed. This analysis is important not only because it allows a comparison between the two devices and optimisation of the operation, but also the switching time is the main operating parameter. An appropriate choice of this parameter enables the fulfilment of the process constraints. The level of the conversions achieved, the more uniform temperature profiles, the uniformity ofcatalyst exploitation and the much simpler mode of operation imposed the RN as a much more suitable device for SCR of NOx with ammonia, in usual operation and also in the perspective of control strategy implementation. Theoretical simplified models have also been proposed in order to describe the forced unsteady state reactors performance and to estimate their internal temperature and concentration profiles. The general idea was to extend the study of catalytic reactor dynamics taking into account the perspectives that haven't been analyzed yet. The experimental investigation ofRN revealed a good agreement between the data obtained by model simulation and the ones obtained experimentally.

Ion-exchange resins as stationary phase in reactive chromatography

Dynamic behavior of bothisothermal and non-isothermal single-column chromatographic reactors with an ion-exchange resin as the stationary phase was investigated. The reactor performance was interpreted by using results obtained when studying the effect of the resin properties on the equilibrium and kinetic phenomena occurring simultaneously in the reactor. Mathematical models were derived for each phenomenon and combined to simulate the chromatographic reactor. The phenomena studied includes phase equilibria in multicomponent liquid mixture¿ion-exchange resin systems, chemicalequilibrium in the presence of a resin catalyst, diffusion of liquids in gel-type and macroporous resins, and chemical reaction kinetics. Above all, attention was paid to the swelling behavior of the resins and how it affects the kinetic phenomena. Several poly(styrene-co-divinylbenzene) resins with different cross-link densities and internal porosities were used. Esterification of acetic acid with ethanol to produce ethyl acetate and water was used as a model reaction system. Choosing an ion-exchange resin with a low cross-link density is beneficial inthe case of the present reaction system: the amount of ethyl acetate as well the ethyl acetate to water mole ratio in the effluent stream increase with decreasing cross-link density. The enhanced performance of the reactor is mainly attributed to increasing reaction rate, which in turn originates from the phase equilibrium behavior of the system. Also mass transfer considerations favor the use ofresins with low cross-link density. The diffusion coefficients of liquids in the gel-type ion-exchange resins were found to fall rapidly when the extent of swelling became low. Glass transition of the polymer was not found to significantlyretard the diffusion in sulfonated PS¿DVB ion-exchange resins. It was also shown that non-isothermal operation of a chromatographic reactor could be used to significantly enhance the reactor performance. In the case of the exothermic modelreaction system and a near-adiabatic column, a positive thermal wave (higher temperature than in the initial state) was found to travel together with the reactive front. This further increased the conversion of the reactants. Diffusion-induced volume changes of the ion-exchange resins were studied in a flow-through cell. It was shown that describing the swelling and shrinking kinetics of the particles calls for a mass transfer model that explicitly includes the limited expansibility of the polymer network. A good description of the process was obtained by combining the generalized Maxwell-Stefan approach and an activity model that was derived from the thermodynamics of polymer solutions and gels. The swelling pressure in the resin phase was evaluated by using a non-Gaussian expression forthe polymer chain length distribution. Dimensional changes of the resin particles necessitate the use of non-standard mathematical tools for dynamic simulations. A transformed coordinate system, where the mass of the polymer was used as a spatial variable, was applied when simulating the chromatographic reactor columns as well as the swelling and shrinking kinetics of the resin particles. Shrinking of the particles in a column leads to formation of dead volume on top of the resin bed. In ordinary Eulerian coordinates, this results in a moving discontinuity that in turn causes numerical difficulties in the solution of the PDE system. The motion of the discontinuity was eliminated by spanning two calculation grids in the column that overlapped at the top of the resin bed. The reactive and non-reactive phase equilibrium data were correlated with a model derived from thethermodynamics of polymer solution and gels. The thermodynamic approach used inthis work is best suited at high degrees of swelling because the polymer matrixmay be in the glassy state when the extent of swelling is low.

Perussyyanalyysien käyttö ydinvoimalaitosten tapahtumien tutkinnassa

Tämä diplomityö käsittelee perussyyanalyysien käyttöä Loviisan voimalaitoksella tapahtumien tutkinnassa. Diplomityössä käydään läpi erilaisia maailmalla käytössä olevia perussyyanalyysimenetelmiä. Nykyään Loviisan voimalaitoksella ei ole ohjeessa määritelty miten perussyyanalyysit tulee tehdä, joten diplomityön tulosten perusteella on tarkoitus valita Loviisan voimalaitokselle sopiva menetelmä. Diplomityössä tehdään eräs perussyyanalyysi kahdella eri menetelmällä, jotka tuntuivat sopivilta menetelmiltä Loviisan voimalaitoksille. Perussyyanalyysimenetelmät, joita käytetään, ovat Loviisassa käytetty sovellettu HPES-menetelmä ja yhdysvaltalaisten ydinvoimaviranomaisten kehittämä HPIP-menetelmä. Perussyyanalyysien valmistuttua ja tuloksia vertaillessa tultiin siihen lopputulokseen, ettei Loviisan voimalaitoksella ole tarvetta vaihtaa nykyisin käytössä olevaa mutta ohjeistamatonta menetelmää. Loviisan voimalaitokselle tehtävään uuteen perussyyanalyysiohjeeseen liitetään menetelmäkuvaus sovelletusta HPES-menetelmästä. Menetelmäkuvauksessa kuvataan miten HPES-menetelmää käytetään perussyyanalyysien tekoon Loviisan voimalaitoksella.

Pain-evoked alterations on gingival blood flow and gingival crevicular fluid (GCF) neuropeptide SP and collagenase-2 (MMP-8) levels

Previous studies have demonstrated that clinical pulpal pain can induce the expression of pro-inflammatory neuropeptides in the adjacent gingival crevice fluid (GCF). Vasoactive agents such as substance P (SP) are known to contribute to the inflammatory type of pain and are associated with increased blood flow. More recent animal studies have shown that application of capsaicin on alveolar mucosa provokes pain and neurogenic vasodilatation in the adjacent gingiva. Pain-associated inflammatory reactions may initiate expression of several pro- and anti-inflammatory mediators. Collagenase-2 (MMP-8) has been considered to be the major destructive protease, especially in the periodontitis-affected gingival crevice fluid (GCF). MMP-8 originates mostly from neutrophil leukocytes, the first line of defence cells that exist abundantly in GCF, especially in inflammation. With this background, we wished to clarify the spatial extensions and differences between tooth-pain stimulation and capsaicin-induced neurogenic vasodilatation in human gingiva. Experiments were carried out to study whether tooth stimulation and capsaicin stimulation of alveolar mucosa would induce changes in GCF MMP-8 levels and whether tooth stimulation would release neuropeptide SP in GCF. The experiments were carried out on healthy human volunteers. During the experiments, moderate and high intensity painful tooth stimulation was performed by a constant current tooth stimulator. Moderate tooth stimulation activates A-delta fibres, while high stimulation also activates C-fibres. Painful stimulation of the gingiva was achieved by topical application of capsaicin-moistened filter paper on the mucosal surface. Capsaicin is known to activate selectively nociceptive C-fibres of stimulated tissue. Pain-evoked vasoactive changes in gingivomucosal tissues were mapped by laser Doppler imaging (LDI), which is a sophisticated and non-invasive method for studying e.g. spatial and temporal characteristics of pain- and inflammation-evoked blood flow changes in gingivomucosal tissues. Pain-evoked release of MMP-8 in GCF samples was studied by immunofluorometric assay (IFMA) and Western immunoblotting. The SP levels in GCF were analysed by Enzyme immunoassay (EIA). During the experiments, subjective stimulus-evoked pain responses were determined by a visual analogue pain scale. Unilateral stimulation of alveolar mucosa and attached gingiva by capsaicin evoked a distinct neurogenic vasodilatation in the ipsilateral gingiva, which attenuated rapidly at the midline. Capsaicin stimulation of alveolar mucosa provoked clear inflammatory reactions. In contrast to capsaicin stimuli, tooth stimulation produced symmetrical vasodilatations bilaterally in the gingiva. The ipsilateral responses were significantly smaller during tooth stimulation than during capsaicin stimuli. The current finding – that tooth stimulation evokes bilateral vasodilatation while capsaicin stimulation of the gingiva mainly produces unilateral vasodilatation – emphasises the usefulness of LDI in clarifying spatial features of neurogenic vasoactive changes in the intra-oral tissues. Capsaicin stimulation of the alveolar mucosa induced significant elevations in MMP-8 levels and activation in GCF of the adjacent teeth. During the experiments, no marked changes occurred in MMP-8 levels in the GCF of distantly located teeth. Painful stimulation of the upper incisor provoked elevations in GCF MMP-8 and SP levels of the stimulated tooth. The GCF MMP-8 and SP levels of the non-stimulated teeth were not changed. These results suggest that capsaicin-induced inflammatory reactions in gingivomucosal tissues do not cross the midline in the anterior maxilla. The enhanced reaction found during stimulation of alveolar mucosa indicates that alveolar mucosa is more sensitive to chemical irritants than the attached gingiva. Analysis of these data suggests that capsaicin-evoked neurogenic inflammation in the gingiva can trigger the expression and activation of MMP-8 in GCF of the adjacent teeth. In this study, it is concluded that experimental tooth pain at C-fibre intensity can induce local elevations in MMP-8 and SP levels in GCF. Depending on the role of MMP-8 in inflammation, in addition to surrogated tissue destruction, the elevated MMP-8 in GCF may also reflect accelerated local defensive and anti-inflammatory reactions.

Reconstruction and analysis of surface variation using photometric stereo

In many industrial applications, accurate and fast surface reconstruction is essential for quality control. Variation in surface finishing parameters, such as surface roughness, can reflect defects in a manufacturing process, non-optimal product operational efficiency, and reduced life expectancy of the product. This thesis considers reconstruction and analysis of high-frequency variation, that is roughness, on planar surfaces. Standard roughness measures in industry are calculated from surface topography. A fast and non-contact method to obtain surface topography is to apply photometric stereo in the estimation of surface gradients and to reconstruct the surface by integrating the gradient fields. Alternatively, visual methods, such as statistical measures, fractal dimension and distance transforms, can be used to characterize surface roughness directly from gray-scale images. In this thesis, the accuracy of distance transforms, statistical measures, and fractal dimension are evaluated in the estimation of surface roughness from gray-scale images and topographies. The results are contrasted to standard industry roughness measures. In distance transforms, the key idea is that distance values calculated along a highly varying surface are greater than distances calculated along a smoother surface. Statistical measures and fractal dimension are common surface roughness measures. In the experiments, skewness and variance of brightness distribution, fractal dimension, and distance transforms exhibited strong linear correlations to standard industry roughness measures. One of the key strengths of photometric stereo method is the acquisition of higher frequency variation of surfaces. In this thesis, the reconstruction of planar high-frequency varying surfaces is studied in the presence of imaging noise and blur. Two Wiener filterbased methods are proposed of which one is optimal in the sense of surface power spectral density given the spectral properties of the imaging noise and blur. Experiments show that the proposed methods preserve the inherent high-frequency variation in the reconstructed surfaces, whereas traditional reconstruction methods typically handle incorrect measurements by smoothing, which dampens the high-frequency variation.

Crystallization characteristics of bioactive glasses

Bioactive glasses are excellent candidates for implant materials, because they can form a chemical bond to bone or guide bone growth, depending on the glass composition. Some compositions have even shown soft tissue attachment and antimicrobial effects. So far, most clinical applications are based on monoliths, plates and particulates of different grain sizes. There is a growing interest in special products such as porous implants sintered from microspheres and fibers drawn from preforms or glass melts. The viscosity range at which these are formed coincides with the crystallization temperature range for most bioactive glasses, thus complicating the manufacturing process. In this work, the crystallization tendency and its kinetics for a series of glasses with their compositions within the range of bioactivity were investigated. The factors affecting crystallization and how it is related to composition were studied by means of thermal analysis and hot stage microscopy. The crystal compositions formed during isothermal and non-isothermal heat treatments were analyzed with SEM-EDXA and X-ray diffraction analysis. The temperatures at which sintering and fiber drawing can take place without interfering with crystallization were determined and glass compositions which are suitable for these purposes were established. The bioactivity of glass fibers and partly crystallized glass plates was studied by soaking them in simulated body fluid (SBF). The thickness of silica, calcium and phosphate rich reaction layers on the glass surface after soaking was used as an indication of the bioactivity. The results indicated that the crystallization tendencies of the experimental glasses are strongly dependent on composition. The main factor affecting the crystallization was found to be the alkali oxide content: the higher the alkali oxide content the lower the crystallization temperature. The primary crystalline phase formed at low temperatures in these glasses was sodium calcium silicate. The crystals were found to form through internal nucleation, leading to bulk crystallization. These glasses had high bioactivity in vitro. Even when partially crystalline, they formed typical reaction layers, indicating bioactivity. In fact, sodium calcium silicate crystals were shown to transform in vitro into hydroxyapatite during soaking. However, crystallization should be avoided because it was shown to retard dissolution, bioactivity reactions and complicate fiber drawing process. Glass compositions having low alkali oxide content showed formation of wollastonite crystals on the surface, at about 300°C above the glass transition temperature. The wide range between glass transition and crystallization allowed viscous flow sintering of these compositions. These glasses also withstood the thermal treatments required for fiber drawing processing. Precipitation of calcium and phosphate on fibers of these glasses in SBF suggested that they were osteoconductive. Glasses showing bioactivity crystallize easily, making their hot working challenging. Undesired crystallization can be avoided by choosing suitable compositions and heat treatment parameters, allowing desired product forms to be attained. Small changes in the oxide composition of the glass can have large effects and therefore a thorough understanding of glass crystallization behavior is a necessity for a successful outcome, when designing and manufacturing implants containing bioactive glasses.

Detection of the Vulnerable Atherosclerotic Plaque. Pre-Clinical Evaluation of Novel PET Imaging Probes With Experimental Models

Atherosclerosis is a life-long vascular inflammatory disease and the leading cause of death in Finland and in other western societies. The development of atherosclerotic plaques is progressive and they form when lipids begin to accumulate in the vessel wall. This accumulation triggers the migration of inflammatory cells that is a hallmark of vascular inflammation. Often, this plaque will become unstable and form vulnerable plaque which may rupture causing thrombosis and in the worst case, causing myocardial infarction or stroke. Identification of these vulnerable plaques before they rupture could save lives. At present, in the clinic, there exists no appropriated, non-invasive method for their identification. The aim of this thesis was to evaluate novel positron emission tomography (PET) probes for the detection of vulnerable atherosclerotic plaques and to characterize, two mouse models of atherosclerosis. These studies were performed by using ex vivo and in vivo imaging modalities. The vulnerability of atherosclerotic plaques was evaluated as expression of active inflammatory cells, namely macrophages. Age and the duration of high-fat diet had a drastic impact on the development of atherosclerotic plaques in mice. In imaging of atherosclerosis, 6-month-old mice, kept on high-fat diet for 4 months, showed matured, metabolically active, atherosclerotic plaques. [18F]FDG and 68Ga were accumulated in the areas representative of vulnerable plaques. However, the slow clearance of 68Ga limits its use for the plaque imaging. The novel synthesized [68Ga]DOTA-RGD and [18F]EF5 tracers demonstrated efficient uptake in plaques as compared to the healthy vessel wall, but the pharmacokinetic properties of these tracers were not optimal in used models. In conclusion, these studies resulted in the identification of new strategies for the assessment of plaque stability and mouse models of atherosclerosis which could be used for plaque imaging. In the used probe panel, [18F]FDG was the best tracer for plaque imaging. However, further studies are warranted to clarify the applicability of [18F]EF5 and [68Ga]DOTA-RGD for imaging of atherosclerosis with other experimental models.

Characterisation of Laser Beam and Paper Material Interaction

It is known already from 1970´s that laser beam is suitable for processing paper materials. In this thesis, term paper materials mean all wood-fibre based materials, like dried pulp, copy paper, newspaper, cardboard, corrugated board, tissue paper etc. Accordingly, laser processing in this thesis means all laser treatments resulting material removal, like cutting, partial cutting, marking, creasing, perforation etc. that can be used to process paper materials. Laser technology provides many advantages for processing of paper materials: non-contact method, freedom of processing geometry, reliable technology for non-stop production etc. Especially packaging industry is very promising area for laser processing applications. However, there are only few industrial laser processing applications worldwide even in beginning of 2010´s. One reason for small-scale use of lasers in paper material manufacturing is that there is a shortage of published research and scientific articles. Another problem, restraining the use of laser for processing of paper materials, is colouration of paper material i.e. the yellowish and/or greyish colour of cut edge appearing during cutting or after cutting. These are the main reasons for selecting the topic of this thesis to concern characterization of interaction of laser beam and paper materials. This study was carried out in Laboratory of Laser Processing at Lappeenranta University of Technology (Finland). Laser equipment used in this study was TRUMPF TLF 2700 carbon dioxide laser that produces a beam with wavelength of 10.6 μm with power range of 190-2500 W (laser power on work piece). Study of laser beam and paper material interaction was carried out by treating dried kraft pulp (grammage of 67 g m-2) with different laser power levels, focal plane postion settings and interaction times. Interaction between laser beam and dried kraft pulp was detected with different monitoring devices, i.e. spectrometer, pyrometer and active illumination imaging system. This way it was possible to create an input and output parameter diagram and to study the effects of input and output parameters in this thesis. When interaction phenomena are understood also process development can be carried out and even new innovations developed. Fulfilling the lack of information on interaction phenomena can assist in the way of lasers for wider use of technology in paper making and converting industry. It was concluded in this thesis that interaction of laser beam and paper material has two mechanisms that are dependent on focal plane position range. Assumed interaction mechanism B appears in range of average focal plane position of 3.4 mm and 2.4 mm and assumed interaction mechanism A in range of average focal plane position of 0.4 mm and -0.6 mm both in used experimental set up. Focal plane position 1.4 mm represents midzone of these two mechanisms. Holes during laser beam and paper material interaction are formed gradually: first small hole is formed to interaction area in the centre of laser beam cross-section and after that, as function of interaction time, hole expands, until interaction between laser beam and dried kraft pulp is ended. By the image analysis it can be seen that in beginning of laser beam and dried kraft pulp material interaction small holes off very good quality are formed. It is obvious that black colour and heat affected zone appear as function of interaction time. This reveals that there still are different interaction phases within interaction mechanisms A and B. These interaction phases appear as function of time and also as function of peak intensity of laser beam. Limit peak intensity is the value that divides interaction mechanism A and B from one-phase interaction into dual-phase interaction. So all peak intensity values under limit peak intensity belong to MAOM (interaction mechanism A one-phase mode) or to MBOM (interaction mechanism B onephase mode) and values over that belong to MADM (interaction mechanism A dual-phase mode) or to MBDM (interaction mechanism B dual-phase mode). Decomposition process of cellulose is evolution of hydrocarbons when temperature is between 380- 500°C. This means that long cellulose molecule is split into smaller volatile hydrocarbons in this temperature range. As temperature increases, decomposition process of cellulose molecule changes. In range of 700-900°C, cellulose molecule is mainly decomposed into H2 gas; this is why this range is called evolution of hydrogen. Interaction in this range starts (as in range of MAOM and MBOM), when a small good quality hole is formed. This is due to “direct evaporation” of pulp via decomposition process of evolution of hydrogen. And this can be seen can be seen in spectrometer as high intensity peak of yellow light (in range of 588-589 nm) which refers to temperature of ~1750ºC. Pyrometer does not detect this high intensity peak since it is not able to detect physical phase change from solid kraft pulp to gaseous compounds. As interaction time between laser beam and dried kraft pulp continues, hypothesis is that three auto ignition processes occurs. Auto ignition of substance is the lowest temperature in which it will spontaneously ignite in a normal atmosphere without an external source of ignition, such as a flame or spark. Three auto ignition processes appears in range of MADM and MBDM, namely: 1. temperature of auto ignition of hydrogen atom (H2) is 500ºC, 2. temperature of auto ignition of carbon monoxide molecule (CO) is 609ºC and 3. temperature of auto ignition of carbon atom (C) is 700ºC. These three auto ignition processes leads to formation of plasma plume which has strong emission of radiation in range of visible light. Formation of this plasma plume can be seen as increase of intensity in wavelength range of ~475-652 nm. Pyrometer shows maximum temperature just after this ignition. This plasma plume is assumed to scatter laser beam so that it interacts with larger area of dried kraft pulp than what is actual area of beam cross-section. This assumed scattering reduces also peak intensity. So result shows that assumably scattered light with low peak intensity is interacting with large area of hole edges and due to low peak intensity this interaction happens in low temperature. So interaction between laser beam and dried kraft pulp turns from evolution of hydrogen to evolution of hydrocarbons. This leads to black colour of hole edges.

Thermal modelling of commercial lithium-ion batteries

The accelerating adoption of electrical technologies in vehicles over the recent years has led to an increase in the research on electrochemical energy storage systems, which are among the key elements in these technologies. The application of electrochemical energy storage systems for instance in hybrid electrical vehicles (HEVs) or hybrid mobile working machines allows tolerating high power peaks, leading to an opportunity to downsize the internal combustion engine and reduce fuel consumption, and therefore, CO2 and other emissions. Further, the application of electrochemical energy storage systems provides an option of kinetic and potential energy recuperation. Presently, the lithium-ion (Li-ion) battery is considered the most suitable electrochemical energy storage type in HEVs and hybrid mobile working machines. However, the intensive operating cycle produces high heat losses in the Li-ion battery, which increase its operating temperature. The Li-ion battery operation at high temperatures accelerates the ageing of the battery, and in the worst case, may lead to a thermal runaway and fire. Therefore, an appropriate Li-ion battery cooling system should be provided for the temperature control in applications such as HEVs and mobile working machines. In this doctoral dissertation, methods are presented to set up a thermal model of a single Li-ion cell and a more complex battery module, which can be used if full information about the battery chemistry is not available. In addition, a non-destructive method is developed for the cell thermal characterization, which allows to measure the thermal parameters at different states of charge and in different points of cell surface. The proposed models and the cell thermal characterization method have been verified by experimental measurements. The minimization of high thermal non-uniformity, which was detected in the pouch cell during its operation with a high C-rate current, was analysed by applying a simplified pouch cell 3D thermal model. In the analysis, heat pipes were incorporated into the pouch cell cooling system, and an optimization algorithm was generated for the estimation of the optimalplacement of heat pipes in the pouch cell cooling system. An analysis of the application of heat pipes to the pouch cell cooling system shows that heat pipes significantly decrease the temperature non-uniformity on the cell surface, and therefore, heat pipes were recommended for the enhancement of the pouch cell cooling system.

Ultrasound-assisted electrochemical treatment of wastewaters containing organic pollutants by using novel Ti/Ta2O5-SnO2 electrodes

Advanced oxidation processes (AOPs) are modern methods using reactive hydroxyl radicals for the mineralization of organic pollutants into simple inorganic compounds, such as CO2 and H2O. Among AOPs electrochemical oxidation (EO) is a method suitable for coloured and turbid wastewaters. The degradation of pollutants occurs on electrocatalytic electrodes. The majority of electrodes contain in their structure either expensive materials (diamond and Pt-group metals) or are toxic for the environment compounds (Sb or Pb). One of the main disadvantages of electrochemical method is the polarization and contamination of electrodes due to the deposition of reaction products on their surface, which results in diminishing of the process efficiency. Ultrasound combined with the electrochemical degradation process eliminates electrode contamination because of the continuous mechanical cleaning effect produced by the formation and collapse of acoustic cavitation bubbles near to the electrode surface. Moreover, high frequency ultrasound generates hydroxyl radicals at water sonolysis. Ultrasound-assisted EO is a non-selective method for oxidation of different organic compounds with high degradation efficiencies. The aim of this research was to develop novel sustainable and cost-effective electrodes working as electrocatalysts and test their activity in electrocatalytic oxidation of organic compounds such as dyes and organic acids. Moreover, the goal of the research was to enhance the efficiency of electrocatalytic degradation processes by assisting it with ultrasound in order to eliminate the main drawbacks of a single electrochemical oxidation such as electrodes polarization and passivation. Novel Ti/Ta2O5-SnO2 electrodes were developed and found to be electrocatalytically active towards water (with 5% Ta content, 10 oxide film layers) and organic compounds oxidation (with 7.5% Ta content, 8 oxide film layers) and therefore these electrodes can be applicable in both environmental and energy fields. The synergetic effect of combined electrolysis and sonication was shown while conducting sonoelectrochemical (EO/US) degradation of methylene blue (MB) and formic acid (FA). Complete degradation of MB and FA was achieved after 45 and 120 min of EO/US process respectively in neutral media. Mineralization efficiency of FA over 95% was obtained after 2 h of degradation using high frequency ultrasound (381, 863, 1176 kHz) combined with 9.1 mA/cm2 current density. EO/US degradation of MB provided over 75% mineralization in 8 h. High degradation kinetic rates and mineralization efficiencies of model pollutants obtained in EO/US experiments provide the preconditions for further extrapolation of this treatment method to pilot scale studies with industrial wastewaters.

Effective notch method assessment of misalignment for fatigue loaded welds

 The objective of this thesis work was to assess axial misalignment in fatigue loaded welds using the effective notch method. As a result, the fatigue behaviour of non-load carrying cruciform fillet welded joint under cyclic tensile loading has been studied. Various degrees of axial misalignment have been found in one series of non-load carrying cruciform fillet welded joints used in a laboratory investigation. As a result, it was important to carry out a comprehensive investigation since axial misalignment forms part of thequality of fatigue loaded structure and can reduce the fatigue strength. To extend the study, the correlation between fatigue strength and stress ratio, as well as stress concentration factor, were also studied. Moreover, a closer investigation of place of crack initiation and its dependence on weld sequence and imperfections of test specimen (angular distortion) was studied. For the fatigue class calculations, FEM (finite element method) and the effectivenotch approach are used. The addressed variable is the axial misalignment whichis introduce by modeling the entire joint. Fracture mechanics based calculations are also used and quantitatively compared with effective notch and experimental results.

Quantile Estimation in Non-Stationary Streaming Data

We provide an incremental quantile estimator for Non-stationary Streaming Data. We propose a method for simultaneous estimation of multiple quantiles corresponding to the given probability levels from streaming data. Due to the limitations of the memory, it is not feasible to compute the quantiles by storing the data. So estimating the quantiles as the data pass by is the only possibility. This can be effective in network measurement. To provide the minimum of the mean-squared error of the estimation, we use parabolic approximation and for comparison we simulate the results for different number of runs and using both linear and parabolic approximations.

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.

Modeling of non-isothermalvapor membrane separation with thermodynamic models and generalized mass transfer equations

A rigorous unit operation model is developed for vapor membrane separation. The new model is able to describe temperature, pressure, and concentration dependent permeation as wellreal fluid effects in vapor and gas separation with hydrocarbon selective rubbery polymeric membranes. The permeation through the membrane is described by a separate treatment of sorption and diffusion within the membrane. The chemical engineering thermodynamics is used to describe the equilibrium sorption of vapors and gases in rubbery membranes with equation of state models for polymeric systems. Also a new modification of the UNIFAC model is proposed for this purpose. Various thermodynamic models are extensively compared in order to verify the models' ability to predict and correlate experimental vapor-liquid equilibrium data. The penetrant transport through the selective layer of the membrane is described with the generalized Maxwell-Stefan equations, which are able to account for thebulk flux contribution as well as the diffusive coupling effect. A method is described to compute and correlate binary penetrant¿membrane diffusion coefficients from the experimental permeability coefficients at different temperatures and pressures. A fluid flow model for spiral-wound modules is derived from the conservation equation of mass, momentum, and energy. The conservation equations are presented in a discretized form by using the control volume approach. A combination of the permeation model and the fluid flow model yields the desired rigorous model for vapor membrane separation. The model is implemented into an inhouse process simulator and so vapor membrane separation may be evaluated as an integralpart of a process flowsheet.

Analysis of torque and speed ripple producing non-idealities of frequency converters in electric drives

Electric motors driven by adjustable-frequency converters may produce periodic excitation forces that can cause torque and speed ripple. Interaction with the driven mechanical system may cause undesirable vibrations that affect the system performance and lifetime. Direct drives in sensitive applications, such as elevators or paper machines, emphasize the importance of smooth torque production. This thesis analyses the non-idealities of frequencyconverters that produce speed and torque ripple in electric drives. The origin of low order harmonics in speed and torque is examined. It is shown how different current measurement error types affect the torque. As the application environment, direct torque control (DTC) method is applied to permanent magnet synchronous machines (PMSM). A simulation model to analyse the effect of the frequency converter non-idealities on the performance of the electric drives is created. Themodel enables to identify potential problems causing torque vibrations and possibly damaging oscillations in electrically driven machine systems. The model is capable of coupling with separate simulation software of complex mechanical loads. Furthermore, the simulation model of the frequency converter's control algorithm can be applied to control a real frequency converter. A commercial frequencyconverter with standard software, a permanent magnet axial flux synchronous motor and a DC motor as the load are used to detect the effect of current measurement errors on load torque. A method to reduce the speed and torque ripple by compensating the current measurement errors is introduced. The method is based on analysing the amplitude of a selected harmonic component of speed as a function oftime and selecting a suitable compensation alternative for the current error. The speed can be either measured or estimated, so the compensation method is applicable also for speed sensorless drives. The proposed compensation method is tested with a laboratory drive, which consists of commercial frequency converter hardware with self-made software and a prototype PMSM. The speed and torque rippleof the test drive are reduced by applying the compensation method. In addition to the direct torque controlled PMSM drives, the compensation method can also beapplied to other motor types and control methods.