Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

Reaction Kinetics and Viscosity Modelling in the Fusion Syntheses of Ca- and Ca/Mg-resinates

Rosin is a natural product from pine forests and it is used as a raw material in resinate syntheses. Resinates are polyvalent metal salts of rosin acids and especially Ca- and Ca/Mg- resinates find wide application in the printing ink industry. In this thesis, analytical methods were applied to increase general knowledge of resinate chemistry and the reaction kinetics was studied in order to model the non linear solution viscosity increase during resinate syntheses by the fusion method. Solution viscosity in toluene is an important quality factor for resinates to be used in printing inks. The concept of critical resinate concentration, c crit, was introduced to define an abrupt change in viscosity dependence on resinate concentration in the solution. The concept was then used to explain the non-inear solution viscosity increase during resinate syntheses. A semi empirical model with two estimated parameters was derived for the viscosity increase on the basis of apparent reaction kinetics. The model was used to control the viscosity and to predict the total reaction time of the resinate process. The kinetic data from the complex reaction media was obtained by acid value titration and by FTIR spectroscopic analyses using a conventional calibration method to measure the resinate concentration and the concentration of free rosin acids. A multivariate calibration method was successfully applied to make partial least square (PLS) models for monitoring acid value and solution viscosity in both mid-infrared (MIR) and near infrared (NIR) regions during the syntheses. The calibration models can be used for on line resinate process monitoring. In kinetic studies, two main reaction steps were observed during the syntheses. First a fast irreversible resination reaction occurs at 235 °C and then a slow thermal decarboxylation of rosin acids starts to take place at 265 °C. Rosin oil is formed during the decarboxylation reaction step causing significant mass loss as the rosin oil evaporates from the system while the viscosity increases to the target level. The mass balance of the syntheses was determined based on the resinate concentration increase during the decarboxylation reaction step. A mechanistic study of the decarboxylation reaction was based on the observation that resinate molecules are partly solvated by rosin acids during the syntheses. Different decarboxylation mechanisms were proposed for the free and solvating rosin acids. The deduced kinetic model supported the analytical data of the syntheses in a wide resinate concentration region, over a wide range of viscosity values and at different reaction temperatures. In addition, the application of the kinetic model to the modified resinate syntheses gave a good fit. A novel synthesis method with the addition of decarboxylated rosin (i.e. rosin oil) to the reaction mixture was introduced. The conversion of rosin acid to resinate was increased to the level necessary to obtain the target viscosity for the product at 235 °C. Due to a lower reaction temperature than in traditional fusion synthesis at 265 °C, thermal decarboxylation is avoided. As a consequence, the mass yield of the resinate syntheses can be increased from ca. 70% to almost 100% by recycling the added rosin oil.

Mass transfer and reaction kinetics in single droplet systems

As reactive extraction grown more and more popular in a variety of technological applications, optimizing its performance becomes more and more important. The process of complex formation is affected by a great number of both physical and chemical properties of all the components involved, and sometimes their interference with one another makes improving the effectiveness of such processes very difficult. In this Master’s Theses, the processes of complex formation between the aqueous phase - represented by copper sulfate water solution, and organic phase – represented by Acorga M5640 solvent extractor, were studied in order to establish the effect these components have on reactive extraction performance and to determine which step is bottlenecking the process the most.

Kinetics of selective oxidation of lactose to lactobionic acid over Pd-active carbon catalyst

Laktoosi eli maitosokeri on tärkein ainesosa useimpien nisäkkäiden tuottamassa maidossa. Sitä erotetaan herasta, juustosta ja maidosta. Laktoosia käytetään elintarvike- ja lääketeollisuuden raaka-aineena monissaeri tuotteissa. Lääketeollisuudessa laktoosia käytetään esimerkiksi tablettien täyteaineena. Hapettamalla laktoosia voidaan valmistaa laktobionihappoa, 2-keto-laktobionihappoa ja laktuloosia. Laktobionihappoa käytetään biohajoavien pintojen ja kosmetiikkatuotteiden valmistuksessa, sekä sisäelinten säilöntäliuoksissa, joissa laktobionihappo estää happiradikaalien aiheuttamien kudosvaurioiden syntymistä. Tässä työssä laktoosia hapetettiin laktobionihapoksi sekoittimella varustetussa laboratoriomittakaavaisessa panosreaktorissa käyttäenkatalyyttinä palladiumia aktiivihiilellä. Muutamissa kokeissa katalyytin promoottorina käytettiin vismuttia, joka hidastaa katalyytin deaktivoitumista. Työn tarkoituksena oli saada lisää tietoa laktoosin hapettamisen kinetiikasta. Laktoosin hapettumisessa laktobionihapoksi havaittiin selektiivisyyteen vaikuttavan muunmuassa reaktiolämpötila, paine, pH ja käytetyn katalyytin määrä. Katalyyttiä kierrättämällä eri kokeiden välillä saatiin paremmat konversiot, selektiivisyydet ja saannot. Parhaat koetulokset saatiin hapetettaessa synteettisellä ilmalla 60 oC lämpötilassa ja 1 bar paineessa. Tehdyissä kokeissa pH:n säätö tehtiin manuaalisesti, joten pH ei pysynyt koko ajan haluttuna. Laktoosin konversio oli parhaimmillaan 95 %. Laktobionihapon suhteellinen selektiivisyys oli 100% ja suhteellinen saanto 100 %. Kinetiikan matemaattinen mallinnus tehtiin Modest-ohjelmalla käyttäen kokeista saatuja mittaustuloksia.Ohjelman avulla estimoitiin parametreja ja saatiin matemaattinen malli reaktorille. Tässä työssä tehtiin kineettinen mallinnus myös ravistelureaktorissa tehdyille laktoosin hapetuskokeille, missä pH pysyi koko ajan haluttuna 'in-situ' titrauksen avulla. Työn yhteydessä selvitettiin myös mahdollisuutta käyttää monoliittikatalyyttejä laktoosin hapetusreaktiossa.

Photoinhibition of Photosystem II. Kinetics, Photoprotection and Mechanism

Photosystem II (PSII) is susceptible to light-induced damage defined as photoinhibition. In natural conditions, plants are capable of repairing the photoinhibited PSII by on-going degradation and re-synthesis of the D1 reaction centre protein of PSII. Photoinhibition is induced by both visible and ultraviolet light and photoinhibition occurs under all light intensities with the same efficiency per photon. In my thesis work, I studied the reaction kinetics and mechanism of photoinhibition of PSII, as well as photoprotection in leaves of higher plants. Action spectroscopy was used to identify photoreceptors of photoinhibition. I found that the action spectrum of photoinhibition in vivo shows resemblance to the absorption spectra of manganese model compounds of the oxygen evolving complex (OEC) suggesting a role for manganese as a photoreceptor of photoinhibition under UV and visible light. In order to study the protective effect of non-photochemical quenching, the action spectrum was measured from leaves of wild type Arabidopsis thaliana and two mutants impaired in nonphotochemical quenching of chlorophyll a excitations. The findings of action spectroscopy and simulations of chlorophyll-based photoinhibition mechanisms suggested that quenching of antenna excitations protects less efficiently than would be expected if antenna chlorophylls were the only photoreceptors of photoinhibition. The reaction kinetics of prolonged photoinhibition was studied in leaves of Cucurbita maxima and Capsicum annuum. The results indicated that photoinhibitory decrease in both the oxygen evolution activity and ratio of variable to maximum fluorescence follows firstorder kinetics in vivo. The persistence of first-order kinetics suggests that already photoinhibited reaction centres do not protect against photoinhibition and that the mechanism of photoinhibition does not have a reversible intermediate. When Cucurbita maxima leaves were photoinhibited with saturating single-turnover flashes and continuous light, the light response curve of photoinhibition was found to be essentially a straight line with both types of illumination, suggesting that similar photoinhibition mechanisms might function during illumination with continuous light and during illumination with short flashes.

A reaction engineering approach to homogeneously catalyzed glycerol hydrochlorination

The production of biodiesel through transesterification has created a surplus of glycerol on the international market. In few years, glycerol has become an inexpensive and abundant raw material, subject to numerous plausible valorisation strategies. Glycerol hydrochlorination stands out as an economically attractive alternative to the production of biobased epichlorohydrin, an important raw material for the manufacturing of epoxy resins and plasticizers. Glycerol hydrochlorination using gaseous hydrogen chloride (HCl) was studied from a reaction engineering viewpoint. Firstly, a more general and rigorous kinetic model was derived based on a consistent reaction mechanism proposed in the literature. The model was validated with experimental data reported in the literature as well as with new data of our own. Semi-batch experiments were conducted in which the influence of the stirring speed, HCl partial pressure, catalyst concentration and temperature were thoroughly analysed and discussed. Acetic acid was used as a homogeneous catalyst for the experiments. For the first time, it was demonstrated that the liquid-phase volume undergoes a significant increase due to the accumulation of HCl in the liquid phase. Novel and relevant features concerning hydrochlorination kinetics, HCl solubility and mass transfer were investigated. An extended reaction mechanism was proposed and a new kinetic model was derived. The model was tested with the experimental data by means of regression analysis, in which kinetic and mass transfer parameters were successfully estimated. A dimensionless number, called Catalyst Modulus, was proposed as a tool for corroborating the kinetic model. Reactive flash distillation experiments were conducted to check the commonly accepted hypothesis that removal of water should enhance the glycerol hydrochlorination kinetics. The performance of the reactive flash distillation experiments were compared to the semi-batch data previously obtained. An unforeseen effect was observed once the water was let to be stripped out from the liquid phase, exposing a strong correlation between the HCl liquid uptake and the presence of water in the system. Water has revealed to play an important role also in the HCl dissociation: as water was removed, the dissociation of HCl was diminished, which had a retarding effect on the reaction kinetics. In order to obtain a further insight on the influence of water on the hydrochlorination reaction, extra semi-batch experiments were conducted in which initial amounts of water and the desired product were added. This study revealed the possibility to use the desired product as an ideal “solvent” for the glycerol hydrochlorination process. A co-current bubble column was used to investigate the glycerol hydrochlorination process under continuous operation. The influence of liquid flow rate, gas flow rate, temperature and catalyst concentration on the glycerol conversion and product distribution was studied. The fluid dynamics of the system showed a remarkable behaviour, which was carefully investigated and described. Highspeed camera images and residence time distribution experiments were conducted to collect relevant information about the flow conditions inside the tube. A model based on the axial dispersion concept was proposed and confronted with the experimental data. The kinetic and solubility parameters estimated from the semi-batch experiments were successfully used in the description of mass transfer and fluid dynamics of the bubble column reactor. In light of the results brought by the present work, the glycerol hydrochlorination reaction mechanism has been finally clarified. It has been demonstrated that the reactive distillation technology may cause drawbacks to the glycerol hydrochlorination reaction rate under certain conditions. Furthermore, continuous reactor technology showed a high selectivity towards monochlorohydrins, whilst semibatch technology was demonstrated to be more efficient towards the production of dichlorohydrins. Based on the novel and revealing discoveries brought by the present work, many insightful suggestions are made towards the improvement of the production of αγ-dichlorohydrin on an industrial scale.

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.

Kuplivan leijukerroksen poltto- ja lämpötekninen tarkastelu

Työn tavoitteena oli luoda menetelmä leijukerroskattilan petin lämpötilojen ja toiminta-alueen arvioimiseksi annetuilla lähtötiedoilla. Toimintatapana oli tutkia polttoaineen ominaisuuksien ja kattilan geometrian vaikutuksia petin energiataseeseen ja tätä kautta petin lämpötilaan kattilan ajotapa huomioiden. Työn teoriaosassa esitellään aluksi biopolttoaineiden ominaisuuksia ja leijukerroskattilan toimintaa petin toimintaympäristön selvittämiseksi. Teoriaosan lopussa käsitellään perusteellisesti petin hydrodynamiikkaa ja reaktiokinettikkaa sekä petin alueen mekaanista mitoitusta. Lisäksi polttoaineen syöttötapa ja palaminen petissä tuodaan esille. Tutkimusosassa selvitetään petin energiatasemallin määrittämistä teoriaan ja diplomityön aikana eri kattiloilla suoritettuihin mittauksiin perustuen. Mittauksista on erillinen selvitys seitsemännessä kappaleessa. Lopuksi kuvaillaan petin energiatasetta mallintavan menetelmän toimintaa sekä siihen syötettäviä lähtötietoja. Tuloksena on saatu suppealla alueella toimiva laajennettavissa oleva petin energiataseen malli. Menetelmän testaamiseen käytetään mittaustietoa kahdelta kattilalta, joita ei olla käytetty petin energiatasemallin virityksessä.

Kineettinen tarkastelu keskisakeusalueella toimivalle happivaiheelle

Työn tavoitteena oli tarkastella keskisakeudessa toimivan happidelignifioinnin kinetiikkaa ja vertailla sitä olemassa olevien kineettisten tutkimuksien kanssa. Lisäksi tutkittiin kuinka happivaiheen suodos vaikuttaa delignifiointiasteeseen sekä selektiivisyyteen. Työn kirjallisuusosassa perehdyttiin happidelignifioinnin kineettisiin malleihin, jotka on kehitetty ennustamaan alkalin, paineen, lämpötilan ja ajan vaikutusta ligniinin poistumiselle sekä hiilihydraattiketjujen katkeilulle. Delignifioituminen on nopeaa ensimmäisen kymmenen minuutin aikana ja hidastuu jäännösdelignifiointivaiheessa. Laboratoriotutkimuksissa selvitettiin kuinka prosessimuuttujat vaikuttavat kappaluvun ja viskositeetin muutokseen. Tuloksia tarkasteltiin lineaarisella regressioanalyysillä, jonka avulla muuttujien vaikutukset saadaan selkeästi esille. Happivaiheen ensimmäisessä vaiheessa suurin vaikutus delignifiointiasteeseen sekä viskositeetin alenemiseen on alkaliannoksella. Pienillä hapen annostuksella saavutettiin hyvä selektiivisyys, suurillakin alkaliannoksilla. Lämpötilan vaikutus on lähes olematon, kun paine on pieni. Kun paine kasvaa, kasvaa myös lämpötilan vaikutus delignifiointiasteeseen. Hiilihydraattien depolymerisoituminen lisääntyy lämpötilan kasvaessa lähes lineaarisesti, riippumatta paineesta. Hyvä loppuselektiivisyys saavutetaan kun ensimmäisen vaiheen selektiivisyys on hyvä. Ensimmäisessä vaiheessa saavutettiin hyvä selektiivisyys, kun alkali ei pääse kulumaan loppuun eikä delignifioituminen jatku liian pitkälle. Toisessa vaiheessa selektiivisyys säilyy parhaiten, kun lisättävä alkaliannos on pieni. Lämpötila ja paine vaikuttavat alkaliannokseen verrattuna hyvin vähän loppuselektiivisyyteen.

Implisiittisen kalibroinnin soveltamismahdollisuudet kemiantekniikassa

Työn taustalla oli tavoite parantaa erään teollisuusprosessin toimintaa ja sen tuottoa mallintamalla reaktiovaiheen alussa tapahtuvan välituotteen muodostumisen reaktiokinetiikka sekä perinteisellä tavalla että implisiittisellä kalibroinnilla. Toisena tavoitteena oli selvittää, kuinka implisiittistä kalibrointia voidaan yleisemmin hyödyntää kemiantekniikassa. Implisiittinen kalibrointi on menetelmä, jolla voidaan ratkaista jonkin teoreettisen mallin parametrit suoraan epäsuorasta mittausdatasta (esimerkiksi spektreistä) lähes kokonaan ilman off-line analyysejä. Tämän työn kirjallisuusosassa on esitetty implisiittisen kalibroinnin toimintaperiaate sekä lyhyesti FTIR-spektrometrian perusteita. Työn kokeellisessa osassa on estimoitu tutkitun välituotteen muodostumisen kineettiset parametrit sekä tavanomaisella parametriestimoinnilla että implisiittisellä kalibroinnilla. Lisäksi kokeellisessa osassa on selvitetty lyhyesti tutkitun prosessin FTIR-spektrien lämpötilariippuvuuksia ja esitetty neljä mahdollista uutta sovelluskohdetta implisiittiselle kalibroinnille. Tavanomaisella parametriestimoinnilla saatiin estimoitua varsin yksiselitteiset arvot kineettisille parametreille. Myös mallin sovitus koedataan on hyvä kolmessa kokeessa viidestä. Parametriestimointi implisiittisellä kalibroinnilla onnistui lupaavasti vaikka tulokset eivät ole aivan niin hyviä kuin tavanomaisessa parametriestimoinnissa. Parhaat tulokset implisiittisessä kalibroinnissa saavutettiin suoralla kalibrointitavalla GRR (Generalized Ridge Regression)-kalibrointimenetelmää käyttämällä.

Jäännöshiilen ja kalkkikiven reaktiokinetiikan mallinnus kiertoleijukattilan tulipesässä

Työn tavoitteena oli kehittää mallit jäännöshiilen ja kalkkikiven reaktiokinetiikan sekä jäännöshiilen jauhautumisen ennustamiseksi kiertoleijukattilan tulipesässä. Kehitetyt mallit toimivat tulipesämallin osamalleina. Perustuen mallinnettuihin reaktionopeuksiin ja jauhautumiskäyttäytymiseen tulipesämalli ennustaa kalkkikivihiukkasten rikinsidonnan ja jäännöshiilen jakautumisen erikokoisiksi hiukkasiksi tulipesässä ja tuhkissa. Työssä kehitetyt mallit perustuvat olemassa oleviin kalkkikiven ja polttoaineen reaktiivisuustesteihin laboratorio-kokoluokan leijukerrosreaktorissa. Mallit huomioivat myös tulipesän olosuhteet. Menetelmät kelpoistettiin onnistuneesti kaupallisen kokoluokan kiertoleijukattiloista mitattujen ja tulipesämallilla laskettujen taseiden avulla. Mallien kehittämistä tullaan jatkamaan.

Gas-phase photocatalytic oxidation of volatile organic compounds

Substances emitted into the atmosphere by human activities in urban and industrial areas cause environmental problems such as air quality degradation, respiratory diseases, climate change, global warming, and stratospheric ozone depletion. Volatile organic compounds (VOCs) are major air pollutants, emitted largely by industry, transportation and households. Many VOCs are toxic, and some are considered to be carcinogenic, mutagenic, or teratogenic. A wide spectrum of VOCs is readily oxidized photocatalytically. Photocatalytic oxidation (PCO) over titanium dioxide may present a potential alternative to air treatment strategies currently in use, such as adsorption and thermal treatment, due to its advantageous activity under ambient conditions, although higher but still mild temperatures may also be applied. The objective of the present research was to disclose routes of chemical reactions, estimate the kinetics and the sensitivity of gas-phase PCO to reaction conditions in respect of air pollutants containing heteroatoms in their molecules. Deactivation of the photocatalyst and restoration of its activity was also taken under consideration to assess the practical possibility of the application of PCO to the treatment of air polluted with VOCs. UV-irradiated titanium dioxide was selected as a photocatalyst for its chemical inertness, non-toxic character and low cost. In the present work Degussa P25 TiO2 photocatalyst was mostly used. In transient studies platinized TiO2 was also studied. The experimental research into PCO of following VOCs was undertaken: - methyl tert-butyl ether (MTBE) as the basic oxygenated motor fuel additive and, thus, a major non-biodegradable pollutant of groundwater; - tert-butyl alcohol (TBA) as the primary product of MTBE hydrolysis and PCO; - ethyl mercaptan (ethanethiol) as one of the reduced sulphur pungent air pollutants in the pulp-and-paper industry; - methylamine (MA) and dimethylamine (DMA) as the amino compounds often emitted by various industries. The PCO of VOCs was studied using a continuous-flow mode. The PCO of MTBE and TBA was also studied by transient mode, in which carbon dioxide, water, and acetone were identified as the main gas-phase products. The volatile products of thermal catalytic oxidation (TCO) of MTBE included 2-methyl-1-propene (2-MP), carbon monoxide, carbon dioxide and water; TBA decomposed to 2-MP and water. Continuous PCO of 4 TBA proceeded faster in humid air than dry air. MTBE oxidation, however, was less sensitive to humidity. The TiO2 catalyst was stable during continuous PCO of MTBE and TBA above 373 K, but gradually lost activity below 373 K; the catalyst could be regenerated by UV irradiation in the absence of gas-phase VOCs. Sulphur dioxide, carbon monoxide, carbon dioxide and water were identified as ultimate products of PCO of ethanethiol. Acetic acid was identified as a photocatalytic oxidation by-product. The limits of ethanethiol concentration and temperature, at which the reactor performance was stable for indefinite time, were established. The apparent reaction kinetics appeared to be independent of the reaction temperature within the studied limits, 373 to 453 K. The catalyst was completely and irreversibly deactivated with ethanethiol TCO. Volatile PCO products of MA included ammonia, nitrogen dioxide, nitrous oxide, carbon dioxide and water. Formamide was observed among DMA PCO products together with others similar to the ones of MA. TCO for both substances resulted in the formation of ammonia, hydrogen cyanide, carbon monoxide, carbon dioxide and water. No deactivation of the photocatalyst during the multiple long-run experiments was observed at the concentrations and temperatures used in the study. PCO of MA was also studied in the aqueous phase. Maximum efficiency was achieved in an alkaline media, where MA exhibited high fugitivity. Two mechanisms of aqueous PCO – decomposition to formate and ammonia, and oxidation of organic nitrogen directly to nitrite - lead ultimately to carbon dioxide, water, ammonia and nitrate: formate and nitrite were observed as intermediates. A part of the ammonia formed in the reaction was oxidized to nitrite and nitrate. This finding helped in better understanding of the gasphase PCO pathways. The PCO kinetic data for VOCs fitted well to the monomolecular Langmuir- Hinshelwood (L-H) model, whereas TCO kinetic behaviour matched the first order process for volatile amines and the L-H model for others. It should be noted that both LH and the first order equations were only the data fit, not the real description of the reaction kinetics. The dependence of the kinetic constants on temperature was established in the form of an Arrhenius equation.

Fischer-Tropsch Slurry Bubble Column Reactor Modeling

The literature part of the work reviews overall Fischer-Tropsch process, Fischer-Tropsch reactors and catalysts. Fundamentals of Fischer-Tropsch modeling are also presented. The emphasis is on the reactor unit. Comparison of the reactors and the catalysts is carried out to choose the suitable reactor setup for the modeling work. The effects of the operation conditions are also investigated. Slurry bubble column reactor model operating with cobalt catalyst is developed by taking into account the mass transfer of the reacting components (CO and H2) and the consumption of the reactants in the liquid phase. The effect of hydrostatic pressure and the change in total mole flow rate in gas phase are taken into account in calculation of the solubilities. The hydrodynamics, reaction kinetics and product composition are determined according to literature. The cooling system and furthermore the required heat transfer area and number of cooling tubes are also determined. The model is implemented in Matlab software. Commercial scale reactor setup is modeled and the behavior of the model is investigated. The possible inaccuraries are evaluated and the suggestions for the future work are presented. The model is also integrated to Aspen Plus process simulation software, which enables the usage of the model in more extensive Fischer-Tropsch process simulations. Commercial scale reactor of diameter of 7 m and height of 30 m was modeled. The capacity of the reactor was calculated to be about 9 800 barrels/day with CO conversion of 75 %. The behavior of the model was realistic and results were in the right range. The highest uncertainty to model was estimated to be caused by the determination of the kinetic rate.

Modeling of Gold Cyanidation

The chemistry of gold dissolution in alkaline cyanide solution has continually received attention and new rate equations expressing the gold leaching are still developed. The effect of leaching parameters on gold gold cyanidation is studied in this work in order to optimize the leaching process. A gold leaching model, based on the well-known shrinking-core model, is presented in this work. It is proposed that the reaction takes place at the reacting particle surface which is continuously reduced as the reaction proceeds. The model parameters are estimated by comparing experimental data and simulations. The experimental data used in this work was obtained from Ling et al. (1996) and de Andrade Lima and Hodouin (2005). Two different rate equations, where the unreacted amount of gold is considered in one equation, are investigated. In this work, it is presented that the reaction at the surface is the rate controlling step since there is no internal diffusion limitation. The model considering the effect of non-reacting gold shows that the reaction orders are consistent with the experimental observations reported by Ling et al. (1996) and de Andrade Lima and Hodouin (2005). However, it should be noted that the model obtained in this work is based on assumptions of no side reactions, no solid-liquid mass transfer resistances and no effect from temperature.

Sinkkirikasteen liukenemisnopeuteen vaikuttavat tekijät ja sinkkirikasteen suoraliuotuksen mallintaminen

Työn kirjallisuusosassa selvitettiin eri tekijöiden (lämpötila, paine, mineralogia, partikkelikoko, sekoitus, kiintoainepitoisuus, liuoksen happo-, rauta- ja happipitoisuus) vaikutusta sinkkirikasteen suoraliuotusprosessin tärkeimpiin ilmiöihin (diffuusio, aineensiirto, reaktiokinetiikka). Kirjallisuusosassa kartoitettiin myös tämän hetkistä tietämystä niistä tekijöistä, jotka ovat oleellisia, kun sinkkirikasteen suoraliuotusprosessia mallinnetaan. Näitä tekijöitä ovat: sinkkirikasteiden liuotuksen kemia, sinkkirikasteiden liuotuksen kinetiikka ja mekanismit, kaasuneste aineensiirto ja kiinteäneste aineensiirto. Lisäksi selvitettiin millä tavoin aikaisemmissa tutkimuksissa sinkkirikasteiden suoraliuotusta on mallinnettu. Mallinnusosassa käsiteltiin atmosfääristä sinkkirikasteen suoraliuotusta, jossa hapettimena toimi ferri(III)rauta. Mallintamisessa käytettiin kirjallisuudessa esitettyjä mittaustuloksia ja mallintaminen tehtiin Modest tietokoneohjelmistolla. Työssä tehty atmosfäärisen suoraliuotuksen mallintaminen labo-ratoriomittakaavassa (laimeat liuokset ja pienet kiintoainepitoisuudet) antoi lupaavia tuloksia. Ongelmia mallin ennustuksen kanssa esiintyi pienissä happopitoisuuksissa, alhaisissa lämpötiloissa sekä pienillä ja suurilla partikkeleilla. Työn kirjallisuusosassa tunnistettiin ne ongelmakohdat, jotka vaativat lisätutkimuksia, jotta sinkkirikasteen atmosfääriselle suoraliuotukselle pystytään kehittämään simulointimalli. Näitä ovat: 1. Hapen liukoisuus ja aineensiirto teollisuuden käyttämissä sinkkirikastelietteissä, 2. Sopivien mittaustulosten puuttuminen, jotta atmosfääristä suoraliuotusta, jossa O2 ja Fe3+ toimivat hapettimena voitaisiin mallintaa, 3. Kiinteäneste aineensiirron merkitys sinkkirikasteen suoraliuotuksessa. Mallinnusosassa osoitettiin, että mekanistisella mallintamisella voidaan simuloida sinkkirikasteen atmosfääristä suoraliuotusta ainakin laboratoriomittakaavassa. Työn perusteella voidaan todeta, että sinkkirikasteen atmosfääriselle suoraliuotukselle voidaan työssä ehdotettujen jatkotutkimusten avulla kehittää numeerinen mekanistinen malli, jolla atmosfääristä sinkkirikasteiden suoraliuotusprosessia voidaan simuloida eri olosuhteissa.