Three-dimensional model of bubbling fluidized bed combustion

 Diplomityön tarkoituksena on kehittää kolmiulotteinen malli kerrosleijupoltolle. Työn kirjallisuusosa sisältää seuraavat perusteet kerrosleijupolton tekniikasta: yleistiedot, leijutus- ja palamisilmiöt, kiinteän aineen ja kaasun sekoittuminen, päästöt ja lämmönsiirto. Lisäksi palamissysteemin mallinnuksen perusteet ja ratkaisumenetelmät ovat esitelty. Työn mallinnusosassa kehitetty koodi on ohjelmoitu Fortran-ohjelmointikielellä. Kehitetty malli perustuu olemassa olevaan malliin kiertoleijupoltosta. Yhtälö kiintoainekonsentraatioprofiilille on vaihdettu ja kiertovirta on poistettu koodista. Mallilla on tehty herkkyystarkasteluja polttoaineen ja kaasun sekoittumisen sekä reaktiokertoimen vaikutukselle. Visualisointi on tehty ohjelmassa Tecplot 360 ja mallinnustuloksia on vertailtu mitattuihin tuloksiin. Mallin laskemattulokset vastaavat hyvin mittaustuloksia ja kokemusperäisiä tietoja; monissa tapauksissa malli pystyy kvantitatiivisesti kuvaamaan parametrien variointia ja kaikissa tapauksissa malli antaa ainakin kvalitatiivisesti oikeita tuloksia. Työhön liittyvän kehityksen ja mallinnuskokemuksen perusteella on tehty ehdotukset mallin tulevaa kehitystä ja mittauksia varten.

Model based study of a calcium oxide looping process for post-combustion carbon dioxide capture

Calcium oxide looping is a carbon dioxide sequestration technique that utilizes the partially reversible reaction between limestone and carbon dioxide in two interconnected fluidised beds, carbonator and calciner. Flue gases from a combustor are fed into the carbonator where calcium oxide reacts with carbon dioxide within the gases at a temperature of 650 ºC. Calcium oxide is transformed into calcium carbonate which is circulated into the regenerative calciner, where calcium carbonate is returned into calcium oxide and a stream of pure carbon dioxide at a higher temperature of 950 ºC. Calcium oxide looping has proved to have a low impact on the overall process efficiency and would be easily retrofitted into existing power plants. This master’s thesis is done in participation to an EU funded project CaOling as a part of the Lappeenranta University of Technology deliverable, reactor modelling and scale-up tools. Thesis concentrates in creating the first model frame and finding the physically relevant phenomena governing the process.

Reaction paper : what? why? how?

Abstract

Detection of bovine foetal DNA from amniotic fluid using the polymerase chain reaction

Selostus: Sikiön DNA:n tunnistaminen naudan sikiövedestä polymeraasiketjureaktion avulla

Universally primed polymerase chain reaction analysis of Fusarium avenaceum isolated from wheat and barley in Finland

Selostus: Vehnästä ja ohrasta eristettyjen F. avenaceum -punahomekantojen analysointi UP-PCR-menetelmällä

Optimisation of data analysis method for a fluidized bed combustion test rig

Nowadays the used fuel variety in power boilers is widening and new boiler constructions and running models have to be developed. This research and development is done in small pilot plants where more faster analyse about the boiler mass and heat balance is needed to be able to find and do the right decisions already during the test run. The barrier on determining boiler balance during test runs is the long process of chemical analyses of collected input and outputmatter samples. The present work is concentrating on finding a way to determinethe boiler balance without chemical analyses and optimise the test rig to get the best possible accuracy for heat and mass balance of the boiler. The purpose of this work was to create an automatic boiler balance calculation method for 4 MW CFB/BFB pilot boiler of Kvaerner Pulping Oy located in Messukylä in Tampere. The calculation was created in the data management computer of pilot plants automation system. The calculation is made in Microsoft Excel environment, which gives a good base and functions for handling large databases and calculations without any delicate programming. The automation system in pilot plant was reconstructed und updated by Metso Automation Oy during year 2001 and the new system MetsoDNA has good data management properties, which is necessary for big calculations as boiler balance calculation. Two possible methods for calculating boiler balance during test run were found. Either the fuel flow is determined, which is usedto calculate the boiler's mass balance, or the unburned carbon loss is estimated and the mass balance of the boiler is calculated on the basis of boiler's heat balance. Both of the methods have their own weaknesses, so they were constructed parallel in the calculation and the decision of the used method was left to user. User also needs to define the used fuels and some solid mass flowsthat aren't measured automatically by the automation system. With sensitivity analysis was found that the most essential values for accurate boiler balance determination are flue gas oxygen content, the boiler's measured heat output and lower heating value of the fuel. The theoretical part of this work concentrates in the error management of these measurements and analyses and on measurement accuracy and boiler balance calculation in theory. The empirical part of this work concentrates on the creation of the balance calculation for the boiler in issue and on describing the work environment.

Plasmapolttoprosessin kehittäminen ja optimointi

Ydinvoimaloidenprimaarivesikierron puhdistukseen käytetään ioninvaihtohartsia. Käytönjälkeen ioninvaihtohartsi luokitellaan matalaja keskiaktiivisiin jätteisiin. Plasmakäsittelyllä käytetyn ioninvaihtohartsin tilavuutta voidaan pienentää sekä sen orgaaninen luonne poistaa. Plasmakäsittelyn tarkoituksena on hapettaa orgaaninen aines oksideiksi, jotka poistuvat prosessista savukaasuina. Epäorgaaninen aines, joka sisältää radioaktiivisen aineksen, on tarkoitus hapettaa oksideiksi ja sulfideiksi, jotka voidaan kerätä talteen tuhkana. Tässä diplomityössä käsitellään käytetyn ioninvaihtohartsin käsittelyyn suunnitellun plasmapolttoprosessin kehittämistä ja optimointia. Ioninvaihtohartsin plasmakäsittelyssä syntyvien reaktiotuotteiden selvittäminen suoritettiin tarkastelemalla ainetaseita sekä aihetta käsitteleviä tutkimuksia. Näiden perusteella parannettiin jäähdytystä, suunniteltiin jatkuvatoiminen syöttömenetelmä sekä laadittiin toimintaalueen reunaehdot laitteistolle. Koelaitteistossa 6,5 kW:n rfteho syötetään sovitinpiirin ja kuparisen induktiokelan kautta plasmaan. Plasmakaasuna on käytetty hapenja argonin seoskaasua. Plasmapolttoa on seurattu massaspektrometrilla, optisella emissiospektrometrilla, lämpösekä painemittareilla. Laskennan ja kokeiden pohjalta selvitettiin optimaalinen seossuhde plasmakaasulle, paineen ja tehon noston vaikutus hartsin polttonopeuteen sekä jatkuvatoimisen syöttömenetelmän edut panostoimiseen syöttöön. Rfgeneraattorin teho rajoitti jatkuvatoimisen polttonopeuden 130 g/h ja hetkellisen polttonopeuden 175 g/h. Radioaktiivisten aineiden pidätys oli 93,5 % cesiumin osalta. Tulosten perusteella 4 kg/h ioninvaihtohartsia polttavan laitteiston tehon lähteeksi tarvitaan 65 kW rfgeneraattori. Palamattoman hartsin ja tuhkan kulkeutuminen partikkelisuodattimille sekä reaktiotuotteena syntyvien rikinoksidien käsittely vaatii vielä jatkotutkimusta.

Numerical Modelling of Combustion and Emission Formation in a Diesel Engine with Different Fuel Injection Characteristics

Päästöjen vähentäminen on ollut viime vuosina tärkeässä osassa polttomoottoreita kehitettäessä.Monet viralliset tahot asettavat uusia tiukempia päästörajoituksia. Päästörajatovat tyypillisesti olleet tiukimmat autoteollisuuden valmistamille pienille nopeakäyntisille diesel-moottoreille, mutta viime aikoina paineita on kohdistunut myös suurempiin keskinopeisiin ja hidaskäyntisiin diesel-moottoreihin. Päästörajat ovat erilaisia riippuen moottorin tyypistä, käytetystä polttoaineesta ja paikasta missä moottoria käytetään johtuen erilaisista paikallisista laeista ja asetuksista. Eniten huomiota diesel-moottorin päästöissä täytyy kohdistaa typen oksideihin, savun muodostukseen sekä partikkeleihin. Laskennallisen virtausmekaniikan (CFD) avulla on hyvät mahdollisuudet tutkia diesel-moottorin sylinterissä tapahtuvia ilmiöitä palamisen aikana. CFD on hyödyllinen työkalu arvioitaessa moottorin suorituskykyä ja päästöjen muodostumista. CFD:llä on mahdollista testata erilaisten parametrien ja geometrioiden vaikutusta ilman kalliita moottorinkoeajoja. CFD:tä voidaan käyttää myös opetustarkoituksessa lisäämään paloprosessin tuntemusta. Tulevaisuudessa palamissimuloinnit CFD:llä tulevat epäilemättä olemaan tärkeä osa moottorin kehityksessä. Tässä diplomityössä on tehty palamissimuloinnit kahteen erilaisilla poittoaineenruiskutuslaitteistoilla varustettuun Wärtsilän keskinopeaan diesel-moottoriin. W46 moottorin ruiskutuslaitteisto on perinteinen mekaanisesti ohjattu pumppusuutin ja W46-CR moottorissa on elektronisesti ohjattu 'common rail' ruiskutuslaitteisto. Näiden moottorien ja käytössä olevien ruiskutusprofiilien lisäksi on simuloinneilla testattu erilaisia uusia ruiskutusprofiileja, jotta erityyppisten profiilien hyvät ja huonot ominaisuudet tulisivat selville. Matalalla kuormalla kiinnostuksen kohteena on nokipäästöjen muodostus ja täydellä kuormalla NOx-päästöjen muodostus ja polttoaineen kulutus. Simulointien tulokset osoittivat, että noen muodostusta matalalla kuormalla voidaan selvästi vähentää monivaiheisella ruiskutuksella, jossa yksi ruiskutusjakso jaetaan kahteen tai useampaan jaksoon. Erityisen tehokas noen vähentämisessä vaikuttaa olevan ns. jälkiruiskutus (post injection). Matalat NOx-päästöt ja hyvä polttoaineen kulutus täydellä kuormalla on mahdollista saavuttaaasteittain nostettavalla ruiskutusnopeudella.

Syttyminen ja palamisen eteneminen partikkelikerroksessa

Syttymistä ja palamisen etenemistä partikkelikerroksessa tutkitaan paloturvallisuuden parantamista sekä kiinteitä polttoaineita käyttävien polttolaitteiden toiminnan tuntemista ja kehittämistä varten. Tässä tutkimuksessa on tavoitteena kerätä yhteen syttymiseen ja liekkirintaman etenemiseen liittyviä kokeellisia ja teoreettisia tutkimustuloksia, jotka auttavat kiinteäkerrospoltto- ja -kaasutus-laitteiden kehittämisessä ja suunnittelussa. Työ on esitutkimus sitä seuraavalle kokeelliselle ja teoreettiselle osalle. Käsittelyssä keskitytään erityisesti puuperäisiin polttoaineisiin. Hiilidioksidipäästöjen vähentämistavoitteet sekä kiinteiden jätteiden energiakäytön lisääminen ja kaatopaikalle viennin vähentäminen aiheuttavat lähitulevaisuudessa kerrospolton lisääntymistä. Kuljetusmatkojen optimoinnin takia joudutaan rakentamaan melko pieniä polttolaitoksia, joissa kerrospolttotekniikka on edullisin vaihtoehto. Syttymispisteellä tarkoitetaan Semenovin määritelmän mukaan tilaa ja ajankohtaa, jolloin polttoaineen ja hapen reaktioissa muodostuva nettoenergia aikayksikössä on yhtäsuuri kuin ympäristöön siirtyvä nettoenergiavirta. Itsesyttyminen tarkoittaa syttymistä ympäristön lämpötilan tai paineen suurenemisen seurauksena. Pakotettu syttyminen tapahtuu, kun syttymispisteen läheisyydessä on esimerkiksi liekki tai hehkuva kiinteä kappale, joka aiheuttaa paikallisen syttymisen ja syttymisrintaman leviämisen muualle polttoaineeseen. Kokeellinen tutkimus on osoittanut tärkeimmiksi syttymiseen ja syttymisrintaman etenemiseen vaikuttaviksi tekijöiksi polttoaineen kosteuden, haihtuvien aineiden pitoisuuden ja lämpöarvon, partikkelikerroksen huokoisuuden, partikkelien koon ja muodon, polttoaineen pinnalle tulevan säteilylämpövirran tiheyden, kaasun virtausnopeuden kerroksessa, hapen osuuden ympäristössä sekä palamisilman esilämmityksen. Kosteuden lisääntyminen suurentaa syttymisenergiaa ja -lämpötilaa sekä pidentää syttymisaikaa. Mitä enemmän polttoaine sisältää haihtuvia aineita sitä pienemmässä lämpötilassa se syttyy. Syttyminen ja syttymisrintaman eteneminen ovat sitä nopeampia mitä suurempi on polttoaineen lämpöarvo. Kerroksen huokoisuuden kasvun on havaittu suurentavan palamisen etenemisnopeutta. Pienet partikkelit syttyvät yleensä nopeammin ja pienemmässä lämpötilassa kuin suuret. Syttymisrintaman eteneminen nopeutuu partikkelien pinta-ala - tilavuussuhteen kasvaessa. Säteilylämpövirran tiheys on useissa polttosovellutuksissa merkittävin lämmönsiirtotekijä, jonka kasvu luonnollisesti nopeuttaa syttymistä. Ilman ja palamiskaasujen virtausnopeus kerroksessa vaikuttaa konvektiiviseen lämmönsiirtoon ja hapen pitoisuuteen syttymisvyöhykkeellä. Ilmavirtaus voi jäähdyttää ja kuumankaasun virtaus lämmittää kerrosta. Hapen osuuden kasvaminen nopeuttaa syttymistä ja liekkirintaman etenemistä kunnes saavutetaan tila, jota suuremmilla virtauksilla ilma jäähdyttää ja laimentaa reaktiovyöhykettä. Palamisilman esilämmitys nopeuttaa syttymisrintaman etenemistä. Syttymistä ja liekkirintaman etenemistä kuvataan yleensä empiirisillä tai säilyvyysyhtälöihin perustuvilla malleilla. Empiiriset mallit perustuvat mittaustuloksista tehtyihin korrelaatioihin sekä joihinkin tunnettuihin fysikaalisiin lainalaisuuksiin. Säilyvyysyhtälöihin perustuvissa malleissa systeemille määritetään massan, energian, liikemäärän ja alkuaineiden säilymisyhtälöt, joiden nopeutta kuvaavien siirtoyhtälöiden muodostamiseen käytetään teoreettisella ja kokeellisella tutkimuksella saatuja yhtälöitä. Nämä mallinnusluokat ovat osittain päällekkäisiä. Pintojen syttymistä kuvataan usein säilyvyysyhtälöihin perustuvilla malleilla. Partikkelikerrosten mallinnuksessa tukeudutaan enimmäkseen empiirisiin yhtälöihin. Partikkelikerroksia kuvaavista malleista Xien ja Liangin hiilipartikkelikerroksen syttymiseen liittyvä tutkimus ja Gortin puun ja jätteen polttoon liittyvä reaktiorintaman etenemistutkimus ovat lähimpänä säilyvyysyhtälöihin perustuvaa mallintamista. Kaikissa malleissa joudutaan kuitenkin yksinkertaistamaan todellista tapausta esimerkiksi vähentämällä dimensioita, reaktioita ja yhdisteitä sekä eliminoimalla vähemmän merkittävät siirtomekanismit. Suoraan kerrospolttoa ja -kaasutusta palvelevia syttymisen ja palamisen etenemisen tutkimuksia on vähän. Muita tarkoituksia varten tehtyjen tutkimusten polttoaineet, kerrokset ja ympäristöolosuhteet poikkeavat yleensä selvästi polttolaitteiden vastaavista olosuhteista. Erikokoisten polttoainepartikkelien ja ominaisuuksiltaan erilaisten polttoaineiden seospolttoa ei ole tutkittu juuri ollenkaan. Polttoainepartikkelien muodon vaikutuksesta on vain vähän tutkimusta.Ilman kanavoitumisen vaikutuksista ei löytynyt tutkimuksia.

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.

Wet Oxidation of Concentrated Wastewaters: Process Combination and Reaction Kinetic Modeling

The accumulation of aqueous pollutants is becoming a global problem. The search for suitable methods and/or combinations of water treatment processes is a task that can slow down and stop the process of water pollution. In this work, the method of wet oxidation was considered as an appropriate technique for the elimination of the impurities present in paper mill process waters. It has been shown that, when combined with traditional wastewater treatment processes, wet oxidation offers many advantages. The combination of coagulation and wet oxidation offers a new opportunity for the improvement of the quality of wastewater designated for discharge or recycling. First of all, the utilization of coagulated sludge via wet oxidation provides a conditioning process for the sludge, i.e. dewatering, which is rather difficult to carry out with untreated waste. Secondly, Fe2(SO4)3, which is employed earlier as a coagulant, transforms the conventional wet oxidation process into a catalytic one. The use of coagulation as the post-treatment for wet oxidation can offer the possibility of the brown hue that usually accompanies the partial oxidation to be reduced. As a result, the supernatant is less colored and also contains a rather low amount of Fe ions to beconsidered for recycling inside mills. The thickened part that consists of metal ions is then recycled back to the wet oxidation system. It was also observed that wet oxidation is favorable for the degradation of pitch substances (LWEs) and lignin that are present in the process waters of paper mills. Rather low operating temperatures are needed for wet oxidation in order to destruct LWEs. The oxidation in the alkaline media provides not only the faster elimination of pitch and lignin but also significantly improves the biodegradable characteristics of wastewater that contains lignin and pitch substances. During the course of the kinetic studies, a model, which can predict the enhancements of the biodegradability of wastewater, was elaborated. The model includes lumped concentrations suchas the chemical oxygen demand and biochemical oxygen demand and reflects a generalized reaction network of oxidative transformations. Later developments incorporated a new lump, the immediately available biochemical oxygen demand, which increased the fidelity of the predictions made by the model. Since changes in biodegradability occur simultaneously with the destruction of LWEs, an attempt was made to combine these two facts for modeling purposes.

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ä.

Wet Oxidation of TMP Concentrated Paper Mill Process Water. Kinetics of the reaction

Työn tarkoituksena oli tutkia lämpötilan, paineen, pH:n ja katalyytin vaikutusta paperitehtaan TMP-konsentroidun prosessiveden märkähapetuksessa. Teoriaosio sisältää katsauksen sellu- ja paperiteollisuuteen, jätevesien käsittelyyn, nanosuodatuksen ja märkähapetusprosessin toimintaperiaatteet ja sovellukset hybriditeknologialle nanosuodatus-märkähapetuksessa. Empiirinen osa koostuu märkähapetuskokeista eri lämpötiloissa, paineissa, pH:ssa ja eri katalyyseillä. Työssä tutkittiin näiden vaikutusta kemialliseen hapenkulutukseen (COD), Biologiseen hapenkulutukseen (BOD), Välittömästi saatavana olevan biologisen hapenkulutukseen (IABOD), ligniiniin, täysin orgaanisen hiileen (TOC) ja rasvaliukoisten uuteaineiden (LWEs) pitoisuuteen. Tuloksina kokeellisesta työstä saatiin korkeimmat COD:n alenemat ja BOD/COD (biohajoavuus) suurimmilla lämpötilaolosuhteilla (COD:n alenema 70 % ja BOD/COD 97 % 200 °C:ssa ja hapen 10 bar osapaineella). Tutkimuksessa, jossa selvitettiin hapen osapaineen vaikutusta saatiin tuloksena, että hapen osapaineen kasvu parantaa orgaanisen kuormituksen poistoa: COD poisto oli olosuhteilla130°C, 5bar 5 %, olosuhteilla 130 °C, 15bar 15 %, olosuhteilla 170 °C, 5bar 20 % ja olosuhteilla 170 °C, 15bar 50 %. Lähes täydellinen LWEs –poisto saavutettiin 150 °C ja 10bar olosuhteilla, vaikka tässä lämpötilassa ei saavutettu korkeata orgaanisen kuormituksen poistoa. Emäksinen pH vaikutti suosivan hapettavia reaktioita, koska korkein COD:n poisto saavutettiin näissä olosuhteilla; kuitenkin alkalisen väliaineen tehokkuudelle löydettiin tärkeä lämpötilariippuvuus.

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.