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.

Does the Ingestion of Isoxaflutole Herbicide Affect the Midgut and Salivary Glands of Pentatomidae Predators?

Podisus nigrispinus is a generalist predator naturally occurring in agricultural and forestry systems that effectively contributes to the population balance of phytophagous insects, especially defoliating caterpillars. Histological changes were evaluated in the salivary glands and midgut of P. nigrispinus caused by ingestion of systemic herbicide isoxaflutole. These predator females were fed with leaves of eucalyptus plants, Tenebrio molitor pupae or water, contaminated or not by herbicide. Salivary glands and midguts were dissected, processed and analyzed under a light microscope. Activity level and cell morphology of the salivary glands and midgut showed differences among insects fed on plants, contaminated water or pupae. The epithelia of the salivary gland and midgut of individuals which had no contact with the herbicide showed homogeneous cytoplasm, nucleus with predominance of decondensed chromatin and evident nucleoli, intense cell activity features. As for the insects in contact with contaminated food, they presented undeveloped nucleus and condensed chromatin. The luminal contents of the salivary glands in the contaminated insects had become more acidophilus than in insects without poisoning, as well as having heterogeneous and granular secretion, being more evident in the bioassay in which the insects fed on contaminated water. There was a marked morphological change in the midgut cells in contaminated insects. High degree of apoptosis, disorganization and secretory vacuoles in the epithelial cytoplasm were observed. The apical portion of the midgut cells proved undeveloped, irregular and partially destroyed. It is concluded that isoxaflutole causes morphological changes in the digestive system of the predator P. nigrispinus.

CFD modeling of dispersion water feed in wastewater cleaning application

Fluid particle breakup and coalescence are important phenomena in a number of industrial flow systems. This study deals with a gas-liquid bubbly flow in one wastewater cleaning application. Three-dimensional geometric model of a dispersion water system was created in ANSYS CFD meshing software. Then, numerical study of the system was carried out by means of unsteady simulations performed in ANSYS FLUENT CFD software. Single-phase water flow case was setup to calculate the entire flow field using the RNG k-epsilon turbulence model based on the Reynolds-averaged Navier-Stokes (RANS) equations. Bubbly flow case was based on a computational fluid dynamics - population balance model (CFD-PBM) coupled approach. Bubble breakup and coalescence were considered to determine the evolution of the bubble size distribution. Obtained results are considered as steps toward optimization of the cleaning process and will be analyzed in order to make the process more efficient.

Randzonenveränderung beim Bohren und ihre Auswirkung auf Folgebearbeitungsverfahren

In der Praxis kommt es bei der spanenden Bearbeitung immer wieder zu großen Standwegunterschieden identischer Werkzeuge bei vordergründig identischen Bearbeitungsrandbedingungen. Insbesondere bei Fertigungsschritten, die das Bohren als Vorbearbeitung erfordern, kommt es gelegentlich zu atypischen Verschleißerscheinungen, die auf das Entstehen einer Neuhärtezone an der Werkstückoberfläche beim Bohren zurückgeführt werden. Grundsätzlich sind Randzonenveränderungen eine Folge der mechanischen und thermischen Beanspruchung bei der Bearbeitung. Beim Eindringen des Schneidkeils kommt es zu Kornverzerrungen, welche die Werkstückhärte bis in eine Tiefe von 40 bis 80 µm erhöhen können. Überdies wird die Randzone des Werkstücks durch den Bearbeitungsvorgang und den Spantransport erhitzt und durch den Kühlschmierstoff bzw. die so genannte Selbstabschreckung im Anschluss sehr schnell abgekühlt. So kann es in Abhängigkeit der Randbedingungen zu Gefügeänderungen mit härtesteigernder (Sekundärabschreckung) oder härtemindernder (Anlasseffekte) Wirkung kommen. Nicht zuletzt beeinflussen beide Beanspruchungsarten auch das Ausmaß der Eigenspannungen in der Werkstückoberfläche. In dieser Arbeit werden die beim Kernlochbohren erzeugten Randzonenveränderungen sowie die Standzeit von Folgebearbeitungswerkzeugen, wie Gewindebohrern und Gewindeformern, und deren Abhängigkeit vom Verschleißzustand des Kernlochbohrers untersucht. Weiterhin werden mit Hilfe einer Energiebilanz die Anteile herausgefiltert, die primär die Eigenschaften der Bohrungsrandzone beeinflussen. Dies geschieht mit Hilfe einer mathematischen Modellierung des Bohrprozesses, in der die Einflüsse der Schneidkantengeometrie, der Schneidkantenverrundung, der Schneidkantenfase sowie des Freiflächenverschleißes berücksichtigt werden.

Analysis of energy dissipation in stirred suspension polymerisation reactors using computational fluid dynamics

This work evaluates the spatial distribution of normalised rates of droplet breakage and droplet coalescence in liquidliquid dispersions maintained in agitated tanks at operation conditions normally used to perform suspension polymerisation reactions. Particularly, simulations are performed with multiphase computational fluid dynamics (CFD) models to represent the flow field in liquidliquid styrene suspension polymerisation reactors for the first time. CFD tools are used first to compute the spatial distribution of the turbulent energy dissipation rates (e) inside the reaction vessel; afterwards, normalised rates of droplet breakage and particle coalescence are computed as functions of e. Surprisingly, multiphase simulations showed that the rates of energy dissipation can be very high near the free vortex surfaces, which has been completely neglected in previous works. The obtained results indicate the existence of extremely large energy dissipation gradients inside the vessel, so that particle breakage occurs primarily in very small regions that surround the impeller and the free vortex surface, while particle coalescence takes place in the liquid bulk. As a consequence, particle breakage should be regarded as an independent source term or a boundary phenomenon. Based on the obtained results, it can be very difficult to justify the use of isotropic assumptions to formulate particle population balances in similar systems, even when multiple compartment models are used to describe the fluid dynamic behaviour of the agitated vessel. (C) 2011 Canadian Society for Chemical Engineering

Sm-Nd isotope and Mn/Ca ratios for cleaned Neogloboquadrina pachyderma from ODP Hole 105-647A (Table 1)

The neodymium (Nd) isotope composition of ancient seawater is a potentially useful tracer of changes in continental inputs and ocean circulation on timescales of a few ka. Here we present the first Nd isotope record for seawater using sedimentary foraminifera cleaned using standard oxidative-reductive techniques. The data, along with Mn/Ca ratios, suggest that cleaned foraminifera provide a reliable record of Nd in seawater and hold out the prospect of using Nd in foraminifera to examine changes in seawater that accompany glacial-interglacial climatic cycles. The principal potential problem to be overcome with the use of forams as records of trace elements in ancient seawater is their diagenetic Fe-Mn coatings. These contain large amounts of Nd and other trace elements but can be cleaned off using highly reducing reagents. Mn(Ca ratios for the majority of the cleaned sedimentary foraminifera analysed here lie within the range (10-100 µmol/mol) that has yielded success in studies of transition elements in forams. Mass-balance modelling suggests that for residual Mn/Ca ratios <100 µmol/mol, Nd added to the foram in the coating will never shift the measured Nd isotope composition significantly away from the seawater value acquired by the foram test in the water column. Additionally, Nd concentrations measured in cleaned sedimentary foraminifera are comparable with those for a modern sample that has never encountered diagenetic fluids. Finally, core-top planktonic foraminifera for two sites have Nd isotope compositions that are identical to local surface seawater. The data we present here for Labrador Sea forams over the past 2.5 m.y. are interpreted in terms of changes in the seawater isotopic composition. The data show a pronounced shift from epsilon-Nd values of ~-12 to ~-19 in the period 2.5-1.5 Ma. This change is interpreted to result from the initiation of Northern Hemisphere glaciation and the increased derivation of Labrador Sea Nd via ice-rafting from Archaean terranes in central Canada. In combination with stable isotope and foraminiferal relative species abundance data, the new Nd data are consistent with the surface hydrography of the Labrador Sea being dominated by a fluctuating balance between cold, polar waters containing unradiogenic Nd and warm, subtropical waters containing more radiogenic Nd. The major change in Labrador Sea Nd that is observed in the past 2.5 Ma can, on its own, account for the change in the Nd isotope composition of North Atlantic Deep Water over the same time period.

Kinetic study of the thermal degradation of high density polyethylene

The thermal degradation of high density polyethylene has been modelled by the random breakage of polymer bonds, using a set of population balance equations. A model was proposed in which the population balances were lumped into representative sizes so that the experimentally determined molecular weight distribution of the original polymer could be used as the initial condition. This model was then compared to two different cases of the unlumped population balance which assumed unimolecular initial distributions of 100 and 500 monomer units, respectively. The model that utilised the experimentally determined molecular weight distribution was found to best describe the experimental data. The model fits suggested a second mechanism in addition to random breakage at slow reaction rates. (c) 2005 Elsevier Ltd. All rights reserved.

A priori prediction of aggregation efficiency and rate constant for fluidized bed melt granulation

This paper presents a predictive aggregation rate model for spray fluidized bed melt granulation. The aggregation rate constant was derived from probability analysis of particle–droplet contact combined with time scale analysis of droplet solidification and granule–granule collision rates. The latter was obtained using the principles of kinetic theory of granular flow (KTGF). The predicted aggregation rate constants were validated by comparison with reported experimental data for a range of binder spray rate, binder droplet size and operating granulator temperature. The developed model is particularly useful for predicting particle size distributions and growth using population balance equations (PBEs).

Representing spray zone with cross flow as a well-mixed compartment in a high shear granulator

The spray zone is an important region to control nucleation of granules in a high shear granulator. In this study, a spray zone with cross flow is quantified as a well-mixed compartment in a high shear granulator. Granulation kinetics is quantitatively derived at both particle-scale and spray zone-scale. Two spatial decay rates, DGSDR (droplet-granule spatial decay rate) ζDG and DPSDR (droplet-primary particle spatial decay rate) ζDP, which are functions of volume fraction and diameter of particulate species within the powder bed, are defined to simplify the deduction. It is concluded that in cross flow, explicit analytical results show that the droplet concentration is subject to exponential decay with depth which produces a numerically infinite depth of spray zone in a real penetration process. In a well-mixed spray zone, the depth of the spray zone is 4/(ζDG + ζDP) and π2/3(ζDG + ζDP) in cuboid and cylinder shape, respectively. The first-order droplet-based collision rates of, nucleation rate B0 and rewetting rate RW0 are uncorrelated with the flow pattern and shape of the spray zone. The second-order droplet-based collision rate, nucleated granule-granule collision rate RGG, is correlated with the mixing pattern. Finally, a real formulation case of a high shear granulation process is used to estimate the size of the spray zone. The results show that the spray zone is a thin layer at the powder bed surface. We present, for the first time, the spray zone as a well-mixed compartment. The granulation kinetics of a well-mixed spray zone could be integrated into a Population Balance Model (PBM), particularly to aid development of a distributed model for product quality prediction.

A compartmental CFD-PBM model of high shear wet granulation

The conventional, geometrically lumped description of the physical processes inside a high shear granulator is not reliable for process design and scale-up. In this study, a compartmental Population Balance Model (PBM) with spatial dependence is developed and validated in two lab-scale high shear granulation processes using a 1.9L MiPro granulator and 4L DIOSNA granulator. The compartmental structure is built using a heuristic approach based on computational fluid dynamics (CFD) analysis, which includes the overall flow pattern, velocity and solids concentration. The constant volume Monte Carlo approach is implemented to solve the multi-compartment population balance equations. Different spatial dependent mechanisms are included in the compartmental PBM to describe granule growth. It is concluded that for both cases (low and high liquid content), the adjustment of parameters (e.g. layering, coalescence and breakage rate) can provide a quantitative prediction of the granulation process.

Standardized geo-referenced catch, fishing effort and length-frequency data for the Indian Ocean Bigeye Tuna, Thunnus obesus (1952-2014)

Geo-referenced catch and fishing effort data of the bigeye tuna fisheries in the Indian Ocean over 1952-2014 were analysed and standardized to facilitate population dynamics modelling studies. During this sixty-two years historical period of exploitation, many changes occurred both in the fishing techniques and the monitoring of activity. This study includes a series of processing steps used for standardization of spatial resolution, conversion and standardization of catch and effort units, raising of geo-referenced catch into nominal catch level, screening and correction of outliers, and detection of major catchability changes over long time series of fishing data, i.e., the Japanese longline fleet operating in the tropical Indian Ocean. A total of thirty fisheries were finally determined from longline, purse seine and other-gears data sets, from which 10 longline and four purse seine fisheries represented 96% of the whole historical catch. The geo-referenced records consists of catch, fishing effort and associated length frequency samples of all fisheries.

Multi-impurity adsorption model for modeling crystal purity and shape evolution during crystallization processes in impure media

Peer reviewed

Modelling of spectral properties and population kinetics studies of inertial fusión and laboratory astrophysical plasmas

Fundamental research and modelling in plasma atomic physics continue to be essential for providing basic understanding of many different topics relevant to high-energy-density plasmas. The Atomic Physics Group at the Institute of Nuclear Fusion has accumulated experience over the years in developing a collection of computational models and tools for determining the atomic energy structure, ionization balance and radiative properties of, mainly, inertial fusion and laser-produced plasmas in a variety of conditions. In this work, we discuss some of the latest advances and results of our research, with emphasis on inertial fusion and laboratory-astrophysical applications.

On the use of correlated beta random variables with animal population modelling

We give reasons why demographic parameters such as survival and reproduction rates are often modelled well in stochastic population simulation using beta distributions. In practice, it is frequently expected that these parameters will be correlated, for example with survival rates for all age classes tending to be high or low in the same year. We therefore discuss a method for producing correlated beta random variables by transforming correlated normal random variables, and show how it can be applied in practice by means of a simple example. We also note how the same approach can be used to produce correlated uniform triangular, and exponential random variables. (C) 2008 Elsevier B.V. All rights reserved.

Design, modelling and analysis of a six component force balance for hypervelocity wind tunnel testing

A combination of modelling and analysis techniques was used to design a six component force balance. The balance was designed specifically for the measurement of impulsive aerodynamic forces and moments characteristic of hypervelocity shock tunnel testing using the stress wave force measurement technique. Aerodynamic modelling was used to estimate the magnitude and distribution of forces and finite element modelling to determine the mechanical response of proposed balance designs. Simulation of balance performance was based on aerodynamic loads and mechanical responses using convolution techniques. Deconvolution was then used to assess balance performance and to guide further design modifications leading to the final balance design. (C) 2001 Elsevier Science Ltd. All rights reserved.