Heterogeneous mass transfer in fluidized beds by computational fluid dynamics

The main objective of this research is to estimate and characterize heterogeneous mass transfer coefficients in bench- and pilot-scale fluidized bed processes by the means of computational fluid dynamics (CFD). A further objective is to benchmark the heterogeneous mass transfer coefficients predicted by fine-grid Eulerian CFD simulations against empirical data presented in the scientific literature. First, a fine-grid two-dimensional Eulerian CFD model with a solid and gas phase has been designed. The model is applied for transient two-dimensional simulations of char combustion in small-scale bubbling and turbulent fluidized beds. The same approach is used to simulate a novel fluidized bed energy conversion process developed for the carbon capture, chemical looping combustion operated with a gaseous fuel. In order to analyze the results of the CFD simulations, two one-dimensional fluidized bed models have been formulated. The single-phase and bubble-emulsion models were applied to derive the average gas-bed and interphase mass transfer coefficients, respectively. In the analysis, the effects of various fluidized bed operation parameters, such as fluidization, velocity, particle and bubble diameter, reactor size, and chemical kinetics, on the heterogeneous mass transfer coefficients in the lower fluidized bed are evaluated extensively. The analysis shows that the fine-grid Eulerian CFD model can predict the heterogeneous mass transfer coefficients quantitatively with acceptable accuracy. Qualitatively, the CFD-based research of fluidized bed process revealed several new scientific results, such as parametrical relationships. The huge variance of seven orders of magnitude within the bed Sherwood numbers presented in the literature could be explained by the change of controlling mechanisms in the overall heterogeneous mass transfer process with the varied process conditions. The research opens new process-specific insights into the reactive fluidized bed processes, such as a strong mass transfer control over heterogeneous reaction rate, a dominance of interphase mass transfer in the fine-particle fluidized beds and a strong chemical kinetic dependence of the average gas-bed mass transfer. The obtained mass transfer coefficients can be applied in fluidized bed models used for various engineering design, reactor scale-up and process research tasks, and they consequently provide an enhanced prediction accuracy of the performance of fluidized bed processes.

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.

Estimation of Mass Transfer Velocity Based on Measured Turbulence Parameters

The aim of this study is to quantify the mass transfer velocity using turbulence parameters from simultaneous measurements of oxygen concentration fields and velocity fields. The surface divergence model was considered in more detail, using data obtained for the lower range of beta (surface divergence). It is shown that the existing models that use the divergence concept furnish good predictions for the transfer velocity also for low values of beta, in the range of this study. Additionally, traditional conceptual models, such as the film model, the penetration-renewal model, and the large eddy model, were tested using the simultaneous information of concentration and velocity fields. It is shown that the film and the surface divergence models predicted the mass transfer velocity for all the range of the equipment Reynolds number used here. The velocity measurements showed viscosity effects close to the surface, which indicates that the surface was contaminated with some surfactant. Considering the results, this contamination can be considered slight for the mass transfer predictions. (C) 2009 American Institute of Chemical Engineers AIChE J, 56: 2005-2017; 2010

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.

Effects of geological inhomogeneity on high Rayleigh number steady state heat and mass transfer in fluid-saturated porous media heated from below

A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, hut it is significant on the species concentration distribution in the medium.

Effects of bone remodelling on calcium mass transfer during haemodialysis

Background. During haemodialysis, calcium balance can affect, or be affected by, mineral metabolism. However, when dialysate calcium concentration (d[Ca]) is chosen or kinetic models are employed to calculate calcium balance, bone remodelling is rarely considered. In this study, we examined whether bone remodelling affects calcium mass transfer during haemodialysis. Methods. We dialysed 23 patients using a d[Ca] of 1.0, 1.25, 1.5 or 1.75 mmol/L. Calcium mass transfer was measured and associated with remodelling bone factors. Results. Calcium balance varied widely depending on the d[Ca]. Calcium removal was -578 +/- 389, -468 +/- 563, +46 +/- 400 and +405 +/- 413 mg when a d[Ca] of 1.0, 1.25, 1.5 or 1.75 mmol/L was used, respectively (1.0 and 1.25 VS 1.5 and 1.75 mmol/L, P<0.001; 1.5 vs 1.75 mmol/L, P<0.05). Univariate analysis showed that calcium balance correlated with calcium gradient, parathyroid hormone (PTH), osteocalcin and dialysis vintage. Multivariate analysis revealed that calcium balance was dependent on calcium gradient, PTH and osteocalcin. Conclusions. These results suggest that bone remodelling could affect calcium mass transfer during haemodialysis.

Mass transfer correlations for rotating drum bioreactors

Evaporative cooling is extremely important for large-scale operation of rotating drum bioreactors (RDBs). Outlet water vapour concentrations were measured for a RDB containing wet wheat bran with the aim of determining the mass transfer coefficient for evaporation from the bran bed to the headspace. Mass transfer was expressed as the mass transfer coefficient times the area for transfer per unit volume of void space in the drum. Values of ka' were determined under combinations of aeration superficial velocities ranging from 0.006 to 0.017 ms(-1) and rotation rates ranging from 0 to 9 rpm. Mass transfer coefficients were evaluated using a variety of residence time distributions (RTDs) for flow in the gas phase including plug flow and well-mixed and a Central Jet RTD based on RTD studies. If plug flow is assumed, the degree of holdup at low effective Peclet (Pe(eff)) numbers gives an apparent under-estimate of ka' compared with empirical correlations. Values of ka' calculated using the Central Jet RTD agree well with values of ka' from literature correlations. There was a linear relationship between ka' and effective Peclet number: ka' = 2.32 x 10(-3) Pe(eff). (C) 2002 Elsevier Science B.V. All rights reserved.

Oxygen mass transfer impact on citric acid production by Yarrowia lipolytica from crude glycerol

Production of citric acid from crude glycerol from biodiesel industry, in batch cultures of Yarrowia lipolytica W29 was performed in a lab-scale stirred tank bioreactor in order to assess the effect of oxygen mass transfer rate in this bioprocess. An empirical correlation was proposed to describe oxygen volumetric mass transfer coefficient (kLa) as a function of operating conditions (stirring speed and specific air flow rate) and cellular density. kLa increased according with a power function with specific power input and superficial gas velocity, and slightly decreased with cellular density. The increase of initial kLa from 7 h-1 to 55 h-1 led to 7.8-fold increase of citric acid final concentration. Experiments were also performed at controlled dissolved oxygen (DO) and citric acid concentration increased with DO up to 60% of saturation. Thus, due to the simpler operation setting an optimal kLa than at controlled DO, it can be concluded that kLa is an adequate parameter for the optimization of citric acid production from crude glycerol by Y. lipolytica and to be considered in bioprocess scale-up. Our empirical correlation, considering the operating conditions and cellular density, will be a valid tool for this purpose.

Mass transfer effects on reactive separation processes - discovery of reactive arrheotropes

Membrane reactor, reactive membrane separation, arrheotrope, azeotrope, dusty gas model, esterification, residue curve map, distillation, kinetics, singular point, bifurcation

Heat and mass transfer in tubular inorganic membranes

Inorganic membranes, permeation, diffusion, heat transfer, mass transfer, axial dispersion

Development of a new method for the characterisation of bioreactors concerning heat and mass transfer

Since the specific heat transfer coefficient (UA) and the volumetric mass transfer coefficient (kLa) play an important role for the design of biotechnological processes, different techniques were developed in the past for the determination of these parameters. However, these approaches often use imprecise dynamic methods for the description of stationary processes and are limited towards scale and geometry of the bioreactor. Therefore, the aim of this thesis was to develop a new method, which overcomes these restrictions. This new approach is based on a permanent production of heat and oxygen by the constant decomposition of hydrogen peroxide in continuous mode. Since the degradation of H2O2 at standard conditions only takes place by the support of a catalyst, different candidates were investigated for their potential (regarding safety issues and reaction kinetic). Manganese-(IV)-oxide was found to be suitable. To compensate the inactivation of MnO2, a continuous process with repeated feeds of fresh MnO2 was established. Subsequently, a scale-up was successfully carried out from 100 mL to a 5 litre glass bioreactor (UniVessel®)To show the applicability of this new method for the characterisation of bioreactors, it was compared with common approaches. With the newly established technique as well as with a conventional procedure, which is based on an electrical heat source, specific heat transfer coefficients were measured in the range of 17.1 – 24.8 W/K for power inputs of about 50 – 70 W/L. However, a first proof of concept regarding the mass transfer showed no constant kLa for different dilution rates up to 0.04 h-1.Based on this, consecutive studies concerning the mass transfer should be made with higher volume flows, due to more even inflow rates. In addition, further experiments are advisable, to analyse the heat transfer in single-use bioreactors and in larger common systems.

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequentin practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of thesample, thus producing a transfer of mass from some components to other ones, butpreserving the total mass present in the system. This evolution is traditionally modelledas a system of ordinary di erential equations of the mass of each component. However,this kind of evolution can be decomposed into a compositional change, expressed interms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despiteof some subcompositions behaving linearly.The goal is to study the characteristics of such simplicial systems of di erential equa-tions such as linearity and stability. This is performed extracting the compositional differential equations from the mass equations. Then, simplicial derivatives are expressedin coordinates of the simplex, thus reducing the problem to the standard theory ofsystems of di erential equations, including stability. The characterisation of stabilityof these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and theassociated behaviour of the orbits are the main tools. For a three component system,these orbits can be plotted both in coordinates of the simplex or in a ternary diagram.A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is aradioactive decay