Biofiltration of Air/Styrene and Air/Styrene/Acetone mixtures in a bubble column reactor

The goal of this work was the treatment of polluted waste gases in a bubble column reactor (BCR), in order to determinate the maximum value of reactor’s efficiency (RE), varying the inlet concentration (C in) of the pollutants. The gaseous mixtures studied were: (i) air with styrene and (ii) air with styrene and acetone. The liquid phase used to contain the biomass in the reactor was a basal salt medium (BSM), fundamental for the microorganisms’ development. The reactor used in this project consists of a glass column of 620mm height and inside diameter 75mm. In all essays there were continually measured: pH, dissolved oxygen and liquid’s temperature. Temperature and pH were controlled (T=24ºC, 7.0 ≤ pH ≤ 7.7). In all experiments the liquid volume (including the biomass) used in the reactor was kept constant (1.5L) as well as the total gas flowrate (1 L/min). Concerning the goal of the work, some parameters were calculated: the organic load (OL), removal efficiency (RE), elimination capacity (EC), biomass concentration (xf) and dry biomass concentration (Xdw). In a first series of experiments, the gas mixture used was air with styrene, varying its concentration from 191 mg.m-3 to 6500 mg.m-3.It was concluded that the RE maximum value (97%) was obtained for C in Sty = 4200 mg.m-3. For the maximum tested value of C in Sty, RE obtained was 20%. In a second step, the gaseous mixture included acetone, varying C in Sty between 225 mg.m-3 and 2659 mg.m-3 and C in Ac between 153mg.m-3 and 1389 mg.m-3. The aim of these tests was the determination of C in Ac for which RE was maximum, obtaining C in Ac = 750 mg.m-3. A third series of experiments was performed, in which C in Ac was maintained equal to that value and C in Sty was varied until higher values (5422 mg.m-3). RE maximum values obtained in this last series were 100% for styrene and 40% for acetone. One important conclusion is the fact that the microorganisms available degrade better styrene than acetone. On the ambit of this study, it was possible to identify the species available in biomass: Xanthobacter antotrophicus py2, Enterobacter aerogenes, Nocardia, Corynebacterium Spp., Rhodococcus rhodochrous e Pseudomonas Sp.

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 simulation of an airlift bubble column reactor

The simulation of two-phase flow for an experimental airlift reactor (32-l volume) using commercially available software from Fluent Incorporated is presented here (http://www.fluent.co.uk). Data from the simulation is compared with the experimental data obtained by the tracking of a magnetic particle and analysis of the pressure drop to determine the gas hold-up. Comparisons between vertical velocity and gas hold-up were made for a series of experiments where the superficial gas velocity in the riser was adjusted between 0.01 and 0.075 m s-1. © 2003 Elsevier B.V. All rights reserved.

Three-dimensional gas-liquid simulation of an airlift bubble column reactor

Basic hydrodynamic parameters of an airlift reactor with internal loop were estimated experimentally and simulated using commercially available CFD software from Fluent. Circulation velocity in a 32-dm(3)-airlift reactor was measured using the magnetic tracer method, meanwhile the gas hold-up was obtained by analysis of the pressure drop using the method of inverted U-tube manometers. Comparison of simulated (in two and three dimensions) and experimental data was performed at different superficial gas velocities in the riser.

Flow regime, hydrodynamics, floc size distribution and sludge properties in activated sludge bubble column, air-lift and aerated stirred reactors

This paper presents a comparative study how reactor configuration, sludge loading and air flowrate affect flow regimes, hydrodynamics, floc size distribution and sludge solids-liquid separation properties. Three reactor configurations were studied in bench scale activated sludge bubble column reactor (BCR), air-lift reactor (ALR) and aerated stirred reactor (ASR). The ASR demonstrated the highest capacity of gas holdup and resistance, and homogeneity in flow regimes and shearing forces, resulting in producing large numbers of small and compact floes. The fluid dynamics in the ALR created regularly directed recirculation forces to enhance the gas holdup and sludge flocculation. The BCR distributed a high turbulent flow regime and non-homogeneity in gas holdup and mixing, and generated large numbers of larger and looser floes. The sludge size distributions, compressibility and settleability were significantly influenced by the reactor configurations associated with the flow regimes and hydrodynamics.

Use of computational fluid dynamics in the design of bubble column reactors

Investigations into the modelling techniques that depict the transport of discrete phases (gas bubbles or solid particles) and model biochemical reactions in a bubble column reactor are discussed here. The mixture model was used to calculate gas-liquid, solid-liquid and gasliquid-solid interactions. Multiphase flow is a difficult phenomenon to capture, particularly in bubble columns where the major driving force is caused by the injection of gas bubbles. The gas bubbles cause a large density difference to occur that results in transient multi-dimensional fluid motion. Standard design procedures do not account for the transient motion, due to the simplifying assumptions of steady plug flow. Computational fluid dynamics (CFD) can assist in expanding the understanding of complex flows in bubble columns by characterising the flow phenomena for many geometrical configurations. Therefore, CFD has a role in the education of chemical and biochemical engineers, providing the examples of flow phenomena that many engineers may not experience, even through experimentation. The performance of the mixture model was investigated for three domains (plane, rectangular and cylindrical) and three flow models (laminar, k-e turbulence and the Reynolds stresses). mThis investigation raised many questions about how gas-liquid interactions are captured numerically. To answer some of these questions the analogy between thermal convection in a cavity and gas-liquid flow in bubble columns was invoked. This involved modelling the buoyant motion of air in a narrow cavity for a number of turbulence schemes. The difference in density was caused by a temperature gradient that acted across the width of the cavity. Multiple vortices were obtained when the Reynolds stresses were utilised with the addition of a basic flow profile after each time step. To implement the three-phase models an alternative mixture model was developed and compared against a commercially available mixture model for three turbulence schemes. The scheme where just the Reynolds stresses model was employed, predicted the transient motion of the fluids quite well for both mixture models. Solid-liquid and then alternative formulations of gas-liquid-solid model were compared against one another. The alternative form of the mixture model was found to perform particularly well for both gas and solid phase transport when calculating two and three-phase flow. The improvement in the solutions obtained was a result of the inclusion of the Reynolds stresses model and differences in the mixture models employed. The differences between the alternative mixture models were found in the volume fraction equation (flux and deviatoric stress tensor terms) and the viscosity formulation for the mixture phase.

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.

Hydrodynamics and mass transfer coefficient in activated sludge aerated stirred column reactor: Experimental analysis and modeling

The aerated stirred reactor (ASR) has been widely used in biochemical and wastewater treatment processes. The information describing how the activated sludge properties and operation conditions affect the hydrodynamics and mass transfer coefficient is missing in the literature. The aim of this study was to investigate the influence of flow regime, superficial gas velocity (U-G), power consumption unit (P/V-L), sludge loading, and apparent viscosity (pap) of activated sludge fluid on the mixing time (t(m)), gas hold-up (epsilon), and volumetric mass transfer coefficient (kLa) in an activated sludge aerated stirred column reactor (ASCR). The activated sludge fluid performed a non-Newtonian rheological behavior. The sludge loading significantly affected the fluid hydrodynamics and mass transfer. With an increase in the UG and P/V-L, the epsilon and k(L)a increased, and the t(m), decreased. The E, kLa, and tm,were influenced dramatically as the flow regime changed from homogeneous to heterogeneous patterns. The proposed mathematical models predicted the experimental results well under experimental conditions, indicating that the U-G, P/V-L, and mu(ap) had significant impact on the t(m) epsilon, and k(L)a. These models were able to give the tm, F, and kLa values with an error around +/- 8%, and always less than +/- 10%. (c) 2005 Wiley Periodicals, Inc.

Gas-liquid-solid flow modelling in a bubble column

An alternative approach to the modelling of solid-liquid and gas-liquid-solid flows for a 5:1 height to width aspect ratio bubble column is presented here. A modified transport equation for the volume fraction of a dispersed phase has been developed for the investigation of turbulent buoyancy driven flows (Chem. Eng. Proc., in press). In this study, a modified transport equation has been employed for discrete phase motion considering both solid-liquid and gas-liquid-solid flows. The modelling of the three-phase flow in a bubble column was achieved in the following case: injecting a slug of solid particles into the column for 10 s at a velocity of 0.1 m s-1 and then the gas phase flow was initiated with a superficial gas velocity of 0.02 cm s-1. © 2003 Elsevier B.V. All rights reserved.

CFD and bubble column reactors:simulation and experiment

The simulation of two-phase flow in bubble columns using commercially available software fromFluent Incorporated is presented here. Data from a bubble column with a ratio of height to thecolumn diameter of 5 : 1 are compared with simulations and experimental results for time-averaged velocity and Reynolds stress proles are used to validate transient, two-dimensional simulations.The models are based on multiphase biological reactors with applications in the food industry. An example case of the mass transfer of oxygen through the liquid phase is also presented.

Influence of temperature on continuous high gravity brewing with yeasts immobilized on spent grains

Flavor compounds` formation and fermentative parameters of continuous high gravity brewing with yeasts immobilized on spent grains were evaluated at three different temperatures (7, 10 and 15 degrees C). The assays were performed in a bubble column reactor at constant dilution rate (0.05 h(-1)) and total gas flow rate (240 ml/min of CO(2) and 10 ml/min of air), with high-gravity all-malt wort (15 degrees Plato). The results revealed that as the fermentation temperature was increased from 7 to 15 degrees C, the apparent and real degrees of fermentation, rate of extract consumption, ethanol volumetric productivity and consumption of free amino nitrogen (FAN) increased. In addition, beer produced at 15 degrees C presented a higher alcohols to esters ratio (2.2-2.4:1) similar to the optimum values described in the literature. It was thus concluded that primary high-gravity (15 degrees Plato) all-malt wort fermentation by continuous process with yeasts immobilized on spent grains, can be carried out with a good performance at 15 degrees C.

Retamycin production by immobilized cells of Streptomyces olindensis ICB20 in repeated-batch cultures

Repeated-batch cultures of Ca-alginate immobilized cells of Streptomyces olindensis ICB20 for retamycin production were carried out in two different bioreactors: a basket-type stirred tank reactor (BSTR) and a bubble column reactor (BCR). Higher average values of retamycin content (R) and productivity (P-R) were achieved in the BSTR cultures (about 1.7 AU and 0.031 AU h(-1), respectively) compared to those obtained in the BCR cultures (about 0.6 AU and 0.012 AU h(-1), respectively). The BCR, on the other hand, presented significantly better operation stability than the BSTR, which makes the former much more promising regarding future industrial applications. (C) 2008 Elsevier Ltd. All rights reserved.

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.