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.

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

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

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Effect of antifoams on oxygen mass transfer in pulp and paper mill effluents

Vaahdonestoaineiden haitallinen vaikutus hapen liukenemisnopeuteen biologisen puhdistamon jätevesissä on yleisesti tunnettua. Aineiden eri vaikutusmekanismien takia on silti vaikea etukäteen arvioida, miten ja kuinka paljon aineensiirto muuttuu. Työn tavoitteena oli saada tietoa vaahdonestoaineiden ja muiden pinta-aktiivisten aineiden vaikutuksesta kuplakokoon, kaasun tilavuusosuuteen ja kaasu-neste aineensiirtoon. Työn teoriaosassa on kuvailtu vaahdon muodostumiseen vaikuttavia tekijöitä sekä eri vaahdonestoaineiden vaikutusmekanismeja sellu- ja paperitehtaan jätevedessä. Edelleen on esitetty useita hapen siirtoa estäviä ja parantavia aineita. Työn kokeellisessa osassa tutkittiin kahdenkymmenenviiden eri pinta-aktiivisen aineen vaikutusta hapen liukenemisnopeuteen yksivaiheisessa kuplakolonnissa. Kokeet tehtiin kahdella pitoisuudella, kahdella eri kaasunjakolaitteella ja kolmella eri kaasun tyhjäputkinopeudella. Aineensiirtokokeiden rinnalla tutkittiin jätevesien laatu- ja fysikaalisiaominaisuuksia, niiden vaikutusta hapen liukenemisnopeuteen sekä testattavien koeaineiden vaikutusta fysikaalisin ominaisuuksiin. Kokeet osoittavat että pinta-aktiivisten aineiden vaikutus hapen aineensiirtoon vaihtelee riippuen kaasunjakolaitteesta ja aineen pitoisuuksista. Testatuista vaahdonestonaineista pienin negatiivinen vaikutus oli aineella AT 35 ja positiivinen vaikutustodettiin olevan vaahdonestoaineiden komponenteilla: P2, S1, F4 ja T9.

Mittausmenetelmien kehittaminen kuplakolonnille

Kuplakolonnireaktoreiden CFD-mallinnus on talla hetkella voimakkaasti kehittyva tutkimusalue. Kaksifaasivirtauksen luotettava simulointi ja mallintaminen on haastavaa kuplakolonnireaktorissa tapahtuvien ilmioiden monimutkaisuuden vuoksi. Reaktorin kayttaytymiseen vaikuttavat tekijat, kuten kolonnin hydrodynamiikka ja aineensiirto, tulee tuntea hyvin ennen mallien tekoa. Tassa tyossa on kokeellisesti tutkittu erilaisten mittausmenetelmien soveltuvuutta kuplakolonnin hydrodynamiikan tutkimiseen. Mittausmenetelmissa on keskitytty erityisesti CFD-mallien vaatimiin paikallisiin mittauksiin. Lisaksi tyossa on arvioitu mittausmenetelmien soveltuvuutta j a luotettavuutta CFD-mallien validointiin. Tyon kirjallisuusosassa on perehdytty kuplakolonnireaktorin hydrodynaamiseen kayttaytymiseen ja siihen vaikuttaviin tekijoihin. Naita ovat mm. reaktorityypit, kaasun dispergointi, virtaustyypit ja -alueet, kaasun tilavuusosuus, kaasukuplan koko ja kuplan nousunopeus. Mittauksia tehtiin kahdessa erikokoisessa kuplakolonnissa, joista pienemman halkaisija oli 0,078 m ja suuremman 0,182 m. Molempien kolonnien nestepinnan korkeus oli 4,62 m. Mittaukset tehtiin vesijohtovedella ja epaorgaanisella prosessiliuoksella. Hydrodynaamisista ominaisuuksista mitattiin kaasun tilavuusosuus, kaasukuplan koko seka kaasukuplan nousunopeus. Kaasun tilavuusosuusmittaukset tehtiin paaasiassa paine-eromittauksella ja joissakin tapauksissa pinnanmittausmenetelmalla. Kuplakoko- ja kuplan nousunopeusmittaukset tehtiin suumopeusvideokameralla ja laser Doppler-anemometrilla. Mittauksissa kaytettiin kahdeksaa erilaista kaasunjakolaitetta, joilla selvitettiin kaasunjakolaitteen ominaisuuksien vaikutusta kolonnin hydrodynamiikkaan. Tuloksista havaittiin, etta nestefaasin ominaisuuksilla oli suuri vaikutus kolonnin hydrodynaamiseen kayttaytymiseen. En kaasunjakolaitteilla vesijohtovedella mitatut hydrodynaamiset ominaisuudet eivat poikenneet paljoa toisistaan, kun taas prosessiliuoksella kaasunjakolaitteiden valille saatiin huomattavat erot. Mittausmenetelmista laser Doppler-anemometri ei kaytettavissa olleella optiikalla soveltunut kaasukuplien mittaamiseen. Kuplat olivat menetelmalle liian suuria. Suumopeusvideokamerallaja paine-eromittauksella paastiin hyviin tuloksiin.

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.

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.

Design and manufacture of a microdistillation column

Aim of this thesis was to design and manufacture a microdistillation column. The literature review part of this thesis covers stainless steels, material processing and basics about engineering design and distillation. The main focus, however, is on the experimental part. Experimental part is divided into five distinct sections: First part is where the device is introduced and separated into three parts. Secondly the device is designed part by part. It consists mostly of detail problem solving, since the first drawings had already been drawn and the critical dimensions decided. Third part is the manufacture, which was not fully completed since the final assembly was left out of this thesis. Fourth part is the test welding for the device, and its analysis. Finally some ideas for further studies are presented. The main goal of this thesis was accomplished. The device only lacks some final assembly but otherwise it is complete. One thing that became clear during the process was how difficult it is to produce small and precise steel parts with conventional manufacturing methods. Internal stresses within steel plates and thermal distortions can easily ruin small steel structures. Designing appropriate welding jigs is an important task for even simple devices. Laser material processing is a promising tool for this kind of steel processing because of the flexibility, good cutting quality and also precise and low heat input when welding. Next step in this project is the final assembly and the actual distillation tests. The tests will be carried out at Helsinki University of Technology.

Evaluation of Bubble Formation and Break Up in Suppression Pools by Using Pattern Recognition Methods

During a possible loss of coolant accident in BWRs, a large amount of steam will be released from the reactor pressure vessel to the suppression pool. Steam will be condensed into the suppression pool causing dynamic and structural loads to the pool. The formation and break up of bubbles can be measured by visual observation using a suitable pattern recognition algorithm. The aim of this study was to improve the preliminary pattern recognition algorithm, developed by Vesa Tanskanen in his doctoral dissertation, by using MATLAB. Video material from the PPOOLEX test facility, recorded during thermal stratification and mixing experiments, was used as a reference in the development of the algorithm. The developed algorithm consists of two parts: the pattern recognition of the bubbles and the analysis of recognized bubble images. The bubble recognition works well, but some errors will appear due to the complex structure of the pool. The results of the image analysis were reasonable. The volume and the surface area of the bubbles were not evaluated. Chugging frequencies calculated by using FFT fitted well into the results of oscillation frequencies measured in the experiments. The pattern recognition algorithm works in the conditions it is designed for. If the measurement configuration will be changed, some modifications have to be done. Numerous improvements are proposed for the future 3D equipment.

Treatment of Fe-containing Acid Mine Waters in Fixed Bed Adsorption Column with Calcite Side-Stone

The acid mining drainage is considered the most significant environmental pollution problem around the world for the extensive formation acidic leachates containing heavy metals. Adsorption is widely used methods in water treatment due to it easy operation and the availability of a wide variety of commercial adsorbent (low cost). The primary goal of this thesis was to investigate the efficiency of neutralizing agents, CaCO3 and CaSiO3, and metal adsorption materials with unmodified limestone from Company Nordkalk Oy. In addition to this, the side materials of limestone mining were tested for iron adsorption from acidic model solution. This study was executed at Lappeenranta University of Technology, Finland. The work utilised fixed-bed adsorption column as the main equipment and large fluidized column. Atomic absorption spectroscopy (AAS) and x-ray diffraction (XRD) was used to determine ferric removal and the composition of material respectively. The results suggest a high potential for the studied materials to be used a low cost adsorbents in acid mine drainage treatment. From the two studied adsorbents, the FS material was more suitable than the Gotland material. Based on the findings, it is recommended that further studies might include detailed analysis of Gotland materials.

Value-based selling of the column connection solution as part of project sales

The objective of this study was to find out how third party influencers can facilitate value-based selling in a network and how suppliers should aim to impact on these third party influencers to facilitate value-based selling. The study considers construction industry, selling the column connection solution and third party influencers. Third party influencers examined in this study were structural designers. The study also aims to find out structural designers’ value drivers and the differences between the market areas that this study related to. The theoretical part of the study focuses on two separate areas. The first part of the theory focuses on a value-based selling concept: what it is, what it requires and what are the main barriers for value-based selling. The second part of the theory examines value creation in networks. The present knowledge over value creation in networks and different network actors are presented. Project marketing is also discussed briefly because this study’s topic, which is highly related to project business. The results reveal structural designers’ value drivers considering the usage of the column connection solution and present ways how suppliers should aim to impact structural designers to facilitate value-based selling. The main result of the study indicates that third party influencers can have a positive impact on facilitating value-based selling. Structural designers are communicating more or less with all the salient actors in different project phases and they can act as sponsors to support the sales of Peikko’s column connection solution and promote solution to other actors involved to the project. This requires that structural designers can understand the actual benefits of how the solution can improve their and their customers’ business.