Bed expansion:extending correlations for fluidized beds with inserts to open fluidized beds

Bed expansion occurs during the operation of gas-fluidized beds and is influenced by particle properties, gas properties and distributor characteristics. It has a significant bearing on heat and mass transfer phenomena within the bed. A method of predicting bed expansion behavior from other fluidizing parameters would be a useful tool in the design process, dispensing with the need for small-scale trials. This study builds on previous work on fluidized beds with vertical inserts to produce a correlation that links a modified particle terminal velocity, minimum fluidizing velocity and distributor characteristics with bed voidage in the relationship with P as the pitch between holes in the perforated distributor plate. © 2007 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Evaporation maps for ternary non-ideal liquid mixtures

This thesis deals with the evaporation of non-ideal liquid mixtures using a multicomponent mass transfer approach. It develops the concept of evaporation maps as a convenient way of representing the dynamic composition changes of ternary mixtures during an evaporation process. Evaporation maps represent the residual composition of evaporating ternary non-ideal mixtures over the full range of composition, and are analogous to the commonly-used residue curve maps of simple distillation processes. The evaporation process initially considered in this work involves gas-phase limited evaporation from a liquid or wetted-solid surface, over which a gas flows at known conditions. Evaporation may occur into a pure inert gas, or into one pre-loaded with a known fraction of one of the ternary components. To explore multicomponent masstransfer effects, a model is developed that uses an exact solution to the Maxwell-Stefan equations for mass transfer in the gas film, with a lumped approach applied to the liquid phase. Solutions to the evaporation model take the form of trajectories in temperaturecomposition space, which are then projected onto a ternary diagram to form the map. Novel algorithms are developed for computation of pseudo-azeotropes in the evaporating mixture, and for calculation of the multicomponent wet-bulb temperature at a given liquid composition. A numerical continuation method is used to track the bifurcations which occur in the evaporation maps, where the composition of one component of the pre-loaded gas is the bifurcation parameter. The bifurcation diagrams can in principle be used to determine the required gas composition to produce a specific terminal composition in the liquid. A simple homotopy method is developed to track the locations of the various possible pseudo-azeotropes in the mixture. The stability of pseudo-azeotropes in the gas-phase limited case is examined using a linearized analysis of the governing equations. Algorithms for the calculation of separation boundaries in the evaporation maps are developed using an optimization-based method, as well as a method employing eigenvectors derived from the linearized analysis. The flexure of the wet-bulb temperature surface is explored, and it is shown how evaporation trajectories cross ridges and valleys, so that ridges and valleys of the surface do not coincide with separation boundaries. Finally, the assumption of gas-phase limited mass transfer is relaxed, by employing a model that includes diffusion in the liquid phase. A finite-volume method is used to solve the system of partial differential equations that results. The evaporation trajectories for the distributed model reduce to those of the lumped (gas-phase limited) model as the diffusivity in the liquid increases; under the same gas-phase conditions the permissible terminal compositions of the distributed and lumped models are the same.

Computational modelling of non-equilibrium condensing steam flows in low-pressure steam turbines

The steam turbines play a significant role in global power generation. Especially, research on low pressure (LP) steam turbine stages is of special importance for steam turbine man- ufactures, vendors, power plant owners and the scientific community due to their lower efficiency than the high pressure steam turbine stages. Because of condensation, the last stages of LP turbine experience irreversible thermodynamic losses, aerodynamic losses and erosion in turbine blades. Additionally, an LP steam turbine requires maintenance due to moisture generation, and therefore, it is also affecting on the turbine reliability. Therefore, the design of energy efficient LP steam turbines requires a comprehensive analysis of condensation phenomena and corresponding losses occurring in the steam tur- bine either by experiments or with numerical simulations. The aim of the present work is to apply computational fluid dynamics (CFD) to enhance the existing knowledge and understanding of condensing steam flows and loss mechanisms that occur due to the irre- versible heat and mass transfer during the condensation process in an LP steam turbine. Throughout this work, two commercial CFD codes were used to model non-equilibrium condensing steam flows. The Eulerian-Eulerian approach was utilised in which the mix- ture of vapour and liquid phases was solved by Reynolds-averaged Navier-Stokes equa- tions. The nucleation process was modelled with the classical nucleation theory, and two different droplet growth models were used to predict the droplet growth rate. The flow turbulence was solved by employing the standard k-ε and the shear stress transport k-ω turbulence models. Further, both models were modified and implemented in the CFD codes. The thermodynamic properties of vapour and liquid phases were evaluated with real gas models. In this thesis, various topics, namely the influence of real gas properties, turbulence mod- elling, unsteadiness and the blade trailing edge shape on wet-steam flows, are studied with different convergent-divergent nozzles, turbine stator cascade and 3D turbine stator-rotor stage. The simulated results of this study were evaluated and discussed together with the available experimental data in the literature. The grid independence study revealed that an adequate grid size is required to capture correct trends of condensation phenomena in LP turbine flows. The study shows that accurate real gas properties are important for the precise modelling of non-equilibrium condensing steam flows. The turbulence modelling revealed that the flow expansion and subsequently the rate of formation of liquid droplet nuclei and its growth process were affected by the turbulence modelling. The losses were rather sensitive to turbulence modelling as well. Based on the presented results, it could be observed that the correct computational prediction of wet-steam flows in the LP turbine requires the turbulence to be modelled accurately. The trailing edge shape of the LP turbine blades influenced the liquid droplet formulation, distribution and sizes, and loss generation. The study shows that the semicircular trailing edge shape predicted the smallest droplet sizes. The square trailing edge shape estimated greater losses. The analysis of steady and unsteady calculations of wet-steam flow exhibited that in unsteady simulations, the interaction of wakes in the rotor blade row affected the flow field. The flow unsteadiness influenced the nucleation and droplet growth processes due to the fluctuation in the Wilson point.

Soviet research and development in mass transfer chemical engineering unit operations.

"PAL 60-16"

Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

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.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

Patch and whole of surface mass-transfer rate measurements on a stepped rotating cylinder electrode

A recently developed whole of surface electroplating technique was used to obtain mass-transfer rates in the separated flow region of a stepped rotating cylinder electrode. These data are compared with previously reported mass-transfer rates obtained with a patch electrode. It was found that the two methods yield different results, where at lower Reynolds numbers, the mass-transfer rate enhancement was noticeably higher for the whole of the surface electrode than for the patch electrode. The location of the peak mass transfer behind the step, as measured with a patch electrode, was reported to be independent of the Reynolds number in previous studies, whereas the whole of the surface electrode shows a definite Reynolds number dependence. Large eddy simulation results for the recirculating region behind a step are used in this work to show that this difference in behavior is related to the existence of a much thinner fluid layer at the wall for which the velocity is a linear junction of distance from the wall. Consequently, the diffusion layer no longer lies well within a laminar sublayer. It is concluded that the patch electrode responds to the wall shear stress for smooth wall flow as well as for the disturbed flow region behind the step. When the whole of the surface is electro-active, the response is to mass transfer even when this is not a sole function of wall shear stress. The results demonstrate that the choice of the mass-transfer measurement technique in corrosion studies can have a significant effect on the results obtained from empirical data.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors. Part B: Heat, momentum and mass transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the heat, momentum and mass transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user-defined function (UDF). The study completes the fast pyrolysis modelling in bubbling fluidised bed reactors.

Effect of impeller type and mechanical agitation on the mass transfer and power consumption aspects of ASBR operation treating synthetic wastewater

The effect of flow type and rotor speed was investigated in a round-bottom reactor with 5 L useful volume containing 2.0 L of granular biomass. The reactor treated 2.0 L of synthetic wastewater with a concentration of 800 mgCOD/L in 8-h cycles at 30 degrees C. Five impellers, commonly used in biological processes, have been employed to this end, namely: a turbine and a paddle impeller with six-vertical-flat-blades, a turbine and a paddle impeller with six-45 degrees-inclined-flat-blades and a three-blade-helix impeller. Results showed that altering impeller type and rotor speed did not significantly affect system stability and performance. Average organic matter removal efficiency was about 84% for filtered samples, total volatile acids concentration was below 20 mgHAc/L and bicarbonate alkalinity a little less than 400 mgCaCO(3)/L for most of the investigated conditions. However, analysis of the first-order kinetic model constants showed that alteration in rotor speed resulted in an increase in the values of the kinetic constants (for instance, from 0.57 h(-1) at 50 rpm to 0.84 h(-1) at 75 rpm when the paddle impeller with six-45 degrees-inclined-flat-blades was used) and that axial flow in mechanically stirred reactors is preferable over radial-flow when the vertical-flat-blade impeller is compared to the inclined-flat-blade impeller (for instance at 75 rpm, from 0.52 h(-1) with the six-flat-blade-paddle impeller to 0.84 h(-1) with the six-45 degrees-inclined-flat-blade-paddle impeller), demonstrating that there is a rotor speed and an impeller type that maximize solid-liquid mass transfer in the reaction medium. Furthermore, power consumption studies in this reduced reactor volume showed that no high power transfer is required to improve mass transfer (less than 0.6 kW/10(3) m(3)). (C) 2008 Elsevier Ltd. All rights reserved.

Computer-Assisted Numerical Analysis for Oxygen and Carbon Dioxide Mass Transfer in Blood Oxygenators

A two-dimensional numeric simulator is developed to predict the nonlinear, convective-reactive, oxygen mass exchange in a cross-flow hollow fiber blood oxygenator. The numeric simulator also calculates the carbon dioxide mass exchange, as hemoglobin affinity to oxygen is affected by the local pH value, which depends mostly on the local carbon dioxide content in blood. Blood pH calculation inside the oxygenator is made by the simultaneous solution of an equation that takes into account the blood buffering capacity and the classical Henderson-Hasselbach equation. The modeling of the mass transfer conductance in the blood comprises a global factor, which is a function of the Reynolds number, and a local factor, which takes into account the amount of oxygen reacted to hemoglobin. The simulator is calibrated against experimental data for an in-line fiber bundle. The results are: (i) the calibration process allows the precise determination of the mass transfer conductance for both oxygen and carbon dioxide; (ii) very alkaline pH values occur in the blood path at the gas inlet side of the fiber bundle; (iii) the parametric analysis of the effect of the blood base excess (BE) shows that V(CO2) is similar in the case of blood metabolic alkalosis, metabolic acidosis, or normal BE, for a similar blood inlet P(CO2), although the condition of metabolic alkalosis is the worst case, as the pH in the vicinity of the gas inlet is the most alkaline; (iv) the parametric analysis of the effect of the gas flow to blood flow ratio (Q(G)/Q(B)) shows that V(CO2) variation with the gas flow is almost linear up to Q(G)/Q(B) = 2.0. V(O2) is not affected by the gas flow as it was observed that by increasing the gas flow up to eight times, the V(O2) grows only 1%. The mass exchange of carbon dioxide uses the full length of the hollow-fiber only if Q(G)/Q(B) > 2.0, as it was observed that only in this condition does the local variation of pH and blood P(CO2) comprise the whole fiber bundle.

Investigation of the mass transfer processes during the desalination of water containing phenol and sodium chloride by electrodialysis

Oxidation processes are used in wastewater treatment when conventional processes are not effective due to the presence of recalcitrant organic contaminants, like phenol. However, the presence of ionic compounds associated with organic pollutants may retard the oxidation. In this work the transport of species contained in an aqueous solution of phenol containing sodium chloride was evaluated in an electrodialysis (ED) system. An experimental study was carried out in which the influence of the process variables on the phenol loss and sodium chloride removal was investigated. Experiments were also performed without current, in order to determine the phenol transfer due to diffusion. The phenol and salt concentration variations in the ED compartments were measured over time, using dedicated procedures and an experimental design to determine the global characteristic parameters. A phenomenological approach was used to relate the phenol, salt and water fluxes with the driving forces (concentration and electric potential gradients). Under ED conditions, two contributions were pointed out for the phenol transport, i.e. diffusion and convection, this latter coming from the water flux due to electroosmosis related to the migration of salts. The fitting of the parameters of the transport equations resulted in good agreement with the experimental results over the range of conditions investigated. (c) 2008 Elsevier B.V. All rights reserved.

Comparison of oxygen mass transfer coefficient in simple and extractive fermentation systems

Aeration and agitation are important variables to ensure effective oxygen transfer rate during aerobic bioprocesses: therefore, the knowledge of the volumetric mass transfer coefficient (k(L)a) is required. In view of selecting the optimum oxygen requirements for extractive fermentation in aqueous two-phase system (ATPS), the k(L)a values in a typical ATPS medium were compared in this work with those in distilled water and in a simple fermentation medium. in the absence of biomass. Aeration and agitation were selected as the independent variables using a 2(2) full factorial design. Both variables showed statistically significant effects on k(L)a, and the highest values of this parameter in both media for simple fermentation (241 s(-1)) and extractive fermentation with ATPS (70.3 s(-1)) were observed at the highest levels of aeration (5 vvm) and agitation (1200 rpm). The k(L)a values were then used to establish mathematical correlations of this response as a function of the process variables. The exponents of the power number (N(3)D(2)) and superficial gas velocity (V(s)) determined in distilled water (alpha = 0.39 and beta = 0.47, respectively) were in reasonable agreement with the ones reported in the literature for several aqueous systems and close to those determined for a simple fermentation medium (alpha=0.38 and beta=0.41). On the other hand, as expected by the increased viscosity in the presence of polyethylene glycol, their values were remarkably higher in a typical medium for extractive fermentation (alpha=0.50 and beta=1.0). A reasonable agreement was found between the experimental data of k(L)a for the three selected systems and the values predicted by the theoretical models, under a wide range of operational conditions. (C) 2009 Elsevier B.V. All rights reserved.

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.

Modeling of non-isothermalvapor membrane separation with thermodynamic models and generalized mass transfer equations

A rigorous unit operation model is developed for vapor membrane separation. The new model is able to describe temperature, pressure, and concentration dependent permeation as wellreal fluid effects in vapor and gas separation with hydrocarbon selective rubbery polymeric membranes. The permeation through the membrane is described by a separate treatment of sorption and diffusion within the membrane. The chemical engineering thermodynamics is used to describe the equilibrium sorption of vapors and gases in rubbery membranes with equation of state models for polymeric systems. Also a new modification of the UNIFAC model is proposed for this purpose. Various thermodynamic models are extensively compared in order to verify the models' ability to predict and correlate experimental vapor-liquid equilibrium data. The penetrant transport through the selective layer of the membrane is described with the generalized Maxwell-Stefan equations, which are able to account for thebulk flux contribution as well as the diffusive coupling effect. A method is described to compute and correlate binary penetrant¿membrane diffusion coefficients from the experimental permeability coefficients at different temperatures and pressures. A fluid flow model for spiral-wound modules is derived from the conservation equation of mass, momentum, and energy. The conservation equations are presented in a discretized form by using the control volume approach. A combination of the permeation model and the fluid flow model yields the desired rigorous model for vapor membrane separation. The model is implemented into an inhouse process simulator and so vapor membrane separation may be evaluated as an integralpart of a process flowsheet.