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.

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.

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

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.

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.

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.

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.

The modelling of buoyancy driven flow in bubble columns

Using the analogy between lateral convection of heat and the two-phase flow in bubble columns, alternative turbulence modelling methods were analysed. The k-ε turbulence and Reynolds stress models were used to predict the buoyant motion of fluids where a density difference arises due to the introduction of heat or a discrete phase. A large height to width aspect ratio cavity was employed in the transport of heat and it was shown that the Reynolds stress model with the use of velocity profiles including the laminar flow solution resulted in turbulent vortices developing. The turbulence models were then applied to the simulation of gas-liquid flow for a 5:1 height to width aspect ratio bubble column. In the case of a gas superficial velocity of 0.02 m s-1 it was determined that employing the Reynolds stress model yielded the most realistic simulation results. © 2003 Elsevier B.V. All rights reserved.

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.

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.