Hydrodynamic analysis of particle collection efficiency: comparing downflow and upflow filtration

Models of the filtration phenomenon describe the mass balance in bed filtration in terms of particle removal mechanisms, and allow for the determination of global particle removal efficiencies. These models are defined in terms of the geometry and characteristic elements of granule collectors, particles and fluid, and also the composition of the balance of forces that act in the particle collector system. This work analyzes particles collection efficiency comparing downflow and upflow direct filtration, taking into account the contribution of the gravitational factor of the settling removal efficiency in future proposal of initial collection efficiency models for upflow filtration. A qualitative analysis is also made of the proposal for the collection efficiency models for particle removal in direct downflow and upflow filtration using a Computational Fluid Dynamics (CFD) tool. This analysis showed a strong influence of gravitational factor in initial collection efficiency (t = 0) of particles, as well as the reasons of their values to be smaller for upflow filtration in comparison with the downflow filtration.

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.

Impacts of epigeic, anecic and endogeic earthworms on metal and metalloid mobility and availability

The introduction of earthworms into soils contaminated with metals and metalloids has been suggested to aid restoration practices. Eisenia veneta (epigeic), Lumbricus terrestris (anecic) and Allolobophora chlorotica (endogeic) earthworms were cultivated in columns containing 900 g soil with 1130, 345, 113 and 131 mg kg1 of As, Cu, Pb and Zn, respectively, for up to 112 days, in parallel with earthworm-free columns. Leachate was produced by pouring water on the soil surface to saturate the soil and generate downflow. Ryegrass was grown on the top of columns to assess metal uptake into biota. Different ecological groups affected metals in the same way by increasing concentrations and free ion activities in leachate, but anecic L. terrestris had the greatest effect by increasing leachate concentrations of As by 267%, Cu by 393%, Pb by 190%, and Zn by 429% compared to earthworm-free columns. Ryegrass grown in earthworm-bearing soil accumulated more metal and the soil microbial community exhibited greater stress. Results are consistent with earthworm enhanced degradation of organic matter leading to release of organically bound elements. The degradation of organic matter also releases organic acids which decrease the soil pH. The earthworms do not appear to carry out a unique process, but increase the rate of a process that is already occurring. The impact of earthworms on metal mobility and availability should therefore be considered when inoculating earthworms into contaminated soils as new pathways to receptors may be created or the flow of metals and metalloids to receptors may be elevated.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Use of the HDS system (hydric depuration with soils) associated to other constructed wetlands techniques for public water supply

Two projects of water treatment for public water supply were developed and operated by using combined systems of constructed wetlands. One of the projects was carried out in the town of Analandia, Sao Paulo, Brazil and wetlands with floating aquatic plants associated to the HDS system were used. Nearly 6480 inhabitants were supplied. The other conducted project was an experimental station in partnership with SABESP (Sao Paulo State Sanitation Agency/Brazil), for the pretreatment of 1700 l.s-1 of waters from the Cotia River, which is used for the population's supply after conventional treatment at the Lower Cotia Water Treatment Station. For this pilot project, wetlands with emergents and floating plants associated to the HDS system were used. The proposed objectives were achieved in both projects.

Influência do leito de contato na eficiência de remoção de ferro de águas subterrâneas

Este trabalho teve como objetivo avaliar a influência da unidade leito de contato na eficiência de remoção ferro em Estação de Tratamento de Água (ETA) tipo desferrização que trata água subterrânea. Os ensaios foram realizados em uma instalação piloto, composta por um aerador tipo tabuleiro, um leito de contato de fluxo ascendente e um filtro rápido de fluxo descendente, que tem capacidade para tratar cerca de 1,7m3/h. A ETA - piloto foi instalada no Setor Profissional do Campus Guamá da UFPA, sendo a água bruta proveniente de manancial subterrâneo, que contem ferro em teores elevados (em torno de 2mg/L). Na ETA - Piloto foram aplicadas taxas de filtração de 180m3/m2xdia, 270m3/m2xdia e 360m3/m2xdia em triplicata em duas fases experimentais: a primeira com leito de contato e segunda sem leito de contato. Foram monitoradas a carreira de filtração, a perda de carga e as variáveis: pH, Cor Aparente, Ferro (Total e Ferroso) e Turbidez. Para a Fase I as eficiências de remoção de ferro total obtidas para as taxas de filtração de 180m3/m2xdia, 270m3/m2xdia e 360m3/m2xdia foram de 93% a 94%, 85% a 95% e 91% a 92%, respectivamente. Para a Fase II as eficiências de remoção de ferro total ficaram em torno de 91% a 93%, 91% a 93% e 88% a 93%, respectivamente. As carreiras de filtração com leito de contato para as taxas de 180m3/m2xdia, 270m3/m2xdia e 360m3/m2xdia apresentaram tempo de duração média de 41h, 30h e 22,7h respectivamente, enquanto que na ausência de leito de contato a duração média foi de 16h, 12h e 4,7h, respectivamente. Dessa forma pode-se concluir que o leito de contato é uma unidade indispensável em uma ETA tipo desferrização, pois melhora o desempenho de sistema, produzindo efluente de boa qualidade e carreiras de filtração mais prolongadas.

TanDEM-X High Resolution DEMs and Their Applications to Flow Modeling

Lava flow modeling can be a powerful tool in hazard assessments; however, the ability to produce accurate models is usually limited by a lack of high resolution, up-to-date Digital Elevation Models (DEMs). This is especially obvious in places such as Kilauea Volcano (Hawaii), where active lava flows frequently alter the terrain. In this study, we use a new technique to create high resolution DEMs on Kilauea using synthetic aperture radar (SAR) data from the TanDEM-X (TDX) satellite. We convert raw TDX SAR data into a geocoded DEM using GAMMA software [Werner et al., 2000]. This process can be completed in several hours and permits creation of updated DEMs as soon as new TDX data are available. To test the DEMs, we use the Harris and Rowland [2001] FLOWGO lava flow model combined with the Favalli et al. [2005] DOWNFLOW model to simulate the 3-15 August 2011 eruption on Kilauea's East Rift Zone. Results were compared with simulations using the older, lower resolution 2000 SRTM DEM of Hawaii. Effusion rates used in the model are derived from MODIS thermal infrared satellite imagery. FLOWGO simulations using the TDX DEM produced a single flow line that matched the August 2011 flow almost perfectly, but could not recreate the entire flow field due to the relatively high DEM noise level. The issues with short model flow lengths can be resolved by filtering noise from the DEM. Model simulations using the outdated SRTM DEM produced a flow field that followed a different trajectory to that observed. Numerous lava flows have been emplaced at Kilauea since the creation of the SRTM DEM, leading the model to project flow lines in areas that have since been covered by fresh lava flows. These results show that DEMs can quickly become outdated on active volcanoes, but our new technique offers the potential to produce accurate, updated DEMs for modeling lava flow hazards.

Proceedings of the 10th Annual Biochemical Engineering Symposium

This the tenth in a series of symposia devoted to talks by students on their biochemical engineering research. The first, third, fifth, and ninth were at Kansas State University in Manhattan, the second and fourth were at the University of Nebraska–Lincoln, the sixth was in Kansas City in conjunction with the 81st American Institute of Chemical Engineers National Meeting, the seventh was at Iowa State University in Ames, and the eighth was held at the University of Missouri–Columbia. Contents"Combined Autohydrolysis-Organosolv Pretreatment of Lignocellulosic Materials," Robert A. Lewis, Colorado State University "An Investigation of Cellulase Activity Assays," Minhhuong Nguyen, University of Missouri–Columbia "Action Pattern of a Xylobiohydrolase from Aspergillus niger," Mary M. Frederick, Iowa State University "Estimation of Heats of Combustion of Biomass from Elemental Analysis Using Available Electron Concepts," Snehal A. Patel, Kansas State University "Design of a Wheat Straw to Ethanol Conversion Facility," Michael M. Meagher, Colorado State University "Effects of Salt, Heat, and Physical Form on the Fermentation of Bananas," Carl Drewel, University of Missouri–Columbia "Gas Hold-up in the Downflow Section of a Split Cylinder Airlift Column," Vasanti Deshpande, Kansas State University "Measurement of Michaelis Constants for Soluble and Immobilized Glucoamylase," Robert A. Lesch, Iowa State University "Kinetics of Alkaline Oxidation and Degradation of Sugars," Alfred R. Fratzke, Iowa State University "Stability of Cereal Protein During Microbial Growth on Grain Dust," Bamidele O. Solomon, Kansas State University

(Appendix) Petrographic information, trace element and Sr isotope ratios of Pito Deep samples

Fluid flow through the axial hydrothermal system at fast spreading ridges is investigated using the Sr-isotopic composition of upper crustal samples recovered from a tectonic window at Pito Deep (NE Easter microplate). Samples from the sheeted dike complex collected away from macroscopic evidence of channelized fluid flow, such as faults and centimeter-scale hydrothermal veins, show a range of 87Sr/86Sr from 0.7025 to 0.7030 averaging 0.70276 relative to a protolith with 87Sr/86Sr of ~0.7024. There is no systematic variation in 87Sr/86Sr with depth in the sheeted dike complex. Comparison of these new data with the two other localities that similar data sets exist for (ODP Hole 504B and the Hess Deep tectonic window) reveals that the extent of Sr-isotope exchange is similar in all of these locations. Models that assume that fluid-rock reaction occurs during one-dimensional (recharge) flow lead to significant decreases in the predicted extent of isotopic modification of the rock with depth in the crust. These model results show systematic misfits when compared with the data that can only be avoided if the fluid flow is assumed to be focused in isolated channels with very slow fluid-rock exchange. In this scenario the fluid at the base of the crust is little modified in 87Sr/86Sr from seawater and thus unlike vent fluids. Additionally, this model predicts that some rocks should show no change from the fresh-rock 87Sr/86Sr, but this is not observed. Alternatively, models in which fluid-rock reaction occurs during upflow (discharge) as well as downflow, or in which fluids are recirculated within the hydrothermal system, can reproduce the observed lack of variation in 87Sr/86Sr with depth in the crust. Minimum time-integrated fluid fluxes, calculated from mass balance, are between 1.5 and 2.6 * 10**6 kg/m**2 for all areas studied to date. However, new evidence from both the rocks and a compilation of vent fluid compositions demonstrates that some Sr is leached from the crust. Because this leaching lowers the fluid 87Sr/86Sr without changing the rock 87Sr/86Sr, these mass balance models must underestimate the time-integrated fluid flux. Additionally, these values do not account for fluid flow that is channelized within the crust.

Pore-water chemistry of ODP sites 111-677 and 111-678

Sites 677 and 678 were drilled on ODP Leg 111 to test hypotheses about the nature and pattern of hydrothermal circulation on a mid-ocean ridge flank. Together with earlier results from DSDP Site 501/504 and several heatflow and piston coring surveys covering a 100-km**2 area surrounding the three drill sites, they confirm that hydrothermal circulation persists in this 5.9-m.y.-old crust, both in basement and through the overlying sediments (Langseth et al., 1988, doi:10.2973/odp.proc.ir.111.102.1988). Profiles of sediment pore-water composition with depth at the three drill sites show both vertical and horizontal gradients. The shapes of the profiles and their variation from one site to another result from a combination of vertical and horizontal diffusion, convection, and reaction in the sediments and basement. Chemical species that are highly reactive in the siliceous-calcareous biogenic sediments include bicarbonate (alkalinity), ammonium, sulfate, manganese, calcium, strontium, lithium, silica, and possibly potassium. Reactions include bacterial sulfate reduction, mobilization of Mn2+, precipitation of CaCO3, and recrystallization of calcareous and siliceous oozes to chalk, limestone, and chert. Species with profiles more affected by reaction in basaltic basement than in the sediments include Mg, Ca, Na, K, and oxygen isotopes. Reaction in basement at 60?C and at higher temperatures has produced a highly altered basement formation water that is uniform in composition over distances of several kilometers. As inferred from the composition of the basal sediment pore water at the three sites, this uniformity extends from up flow zone to downflow zone in basement and the sediments. It exists in spite of large variations in heat flow and depth to basement, apparently as a result of homogenization by hydrothermal circulation in basement. Profiles for chlorinity, Na, Mg, and other species in the sediment pore waters confirm that Site 678, drilled on a localized heatflow high identified by Langseth et al. (1988), is a site of long-lived upwelling of warm water from basement through the sediments at velocities of 1 to 2 mm/yr. The upflow through the anomalously thin sediments is apparently localized above an uplifted fault block in basement. This site and other similar sites in the survey area give rise to lateral diffusion and possibly flow through the sediments, which produces lateral gradients in sediment pore-water composition at sites such as 501/504. The complementary pore-water profiles at the low-heatflow Site 677 2 km to the south indicate that downflow is occurring through the sediments there, at comparable rates of 1 to 2 mm/yr.

Local liquid holdups and hysteresis in a 2-D packed bed using X-ray radiography

An X-ray visualization technique has been used for the quantitative determination of local liquid holdups distribution and liquid holdup hysteresis in a nonwetting two-dimensional (2-D) packed bed. A medical diagnostic X-ray unit has been used to image the local holdups in a 2-D cold model having a random packing of expanded polystyrene beads. An aqueous barium chloride solution was used as a fluid to achieve good contrast on X-ray images. To quantify the local liquid holdup, a simple calibration technique has been developed that can be used for most of the radiological methods such as gamma ray and neutron radiography. The global value of total liquid holdup, obtained by X-ray method, has been compared with two conventional methods: drainage and tracer response. The X-ray technique, after validation, has been used to visualize and quantify, the liquid hysteresis phenomena in a packed bed. The liquid flows in preferred paths or channels that carry droplets/rivulets of increasing size and number as the liquid flow rate is increased. When the flow is reduced, these paths are retained and the higher liquid holdup that persists in these regions leads to the holdup hysteresis effect. Holdup in some regions of the packed bed may be an order of magnitude higher than average at a particular flow rate. (c) 2005 American Institute of Chemical Engineers

Flowback of solids through distribution plates of gas fluidized beds and associated phenomena

This work is concerned with a study of certain phenomena related to the performance and design of distributors in gas fluidized beds with particular regard to flowback of solid particles. The work to be described is divided into two parts. I. In Part one, a review of published material pertaining to distribution plates, including details from the patent specifications, has been prepared. After a chapter on the determination of the incipient fluidizing velocity, the following aspects of multi-orifice distributor plates in gas fluidized beds have been studied: (i) The effect of the distributor on bubble formation related to the way in which even distribution of bubbles on the top surface of the fluidized bed is obtained, e.g. the desirable pressure drop ratio ?PD/?PB for the even distribution of gas across the bed. Ratios of distributor pressure drop ?PD to bed pressure drop at which stable fluidization occurs show reasonable agreement with industrial practice. There is evidence that larger diameter beds tend to be less stable than smaller diameter beds when these are operated with shallow beds. Experiments show that in the presence of the bed the distributor pressure drop is reduced relative to the pressure drop without the bed, and this pressure drop in the former condition is regarded as the appropriate parameter for the design of the distributor. (ii) Experimental measurements of bubble distribution at the surface has been used to indicate maldistribution within the bed. Maldistribution is more likely at low gas flow rates and with distributors having large fractional free area characteristics (i.e. with distributors having low pressure drops). Bubble sizes obtained from this study, as well as those of others, have been successfully correlated. The correlation produced implies the existence of a bubble at the surface of an orifice and its growth by the addition of excess gas from the fluidized bed. (iii) For a given solid system, the amount of defluidized particles stagnating on the distributor plate is influenced by the orifice spacing, bed diameter and gas flow rate, but independent of the initial bed height and the way the orifices are arranged on the distributor plate. II. In Part two, solids flowback through single and multi-orifice distributors in two-dimensional and cylindrical beds of solids fluidized with air has been investigated. Distributors equipped with long cylindrical nozzles have also been included in the study. An equation for the prediction of free flowback of solids through multi-orifice distributors has been derived. Under fluidized conditions two regimes of flowback have been differentiated, namely Jumping and weeping. Data in the weeping regime have been successfully correlated. The limiting gas velocity through the distributor orifices at which flowback is completely excluded is found to be indepnndent of bed height, but a function of distributor design and physical properties of gas and solid used. A criterion for the prediction of this velocity has been established. The decisive advantage of increasing the distributor thickness or using nozzles to minimize solids flowback in fluidized beds has been observed and the opportunity taken to explore this poorly studied subject area. It has been noted, probably for the first time, that with long nozzles, there exists a critical nozzle length above which uncontrollable downflow of solids occurs. A theoretical model for predicting the critical length of a bundle of nozzles in terms of gas velocity through the nozzles has been set up. Theoretical calculations compared favourably with experiments.

Nitrogen removal and sustainability of vertical flow constructed wetlands for small scale wastewater treatment

The challenge for wastewater professionals is to design and operate treatment processes that support human well being and are environmentally sensitive throughout the life-cycle. This research focuses on one technology for small-scale wastewater treatment: the vertical flow constructed wetland (VFCW), which is herein investigated for the capacity to remove ammonium and nitrate nitrogen from wastewater. Hydraulic regime and presence/absence of vegetation are the basis for a three-phase bench scale experiment to determine oxygen transfer and nitrogen fate in VFCWs. Results show that 90% NH4+-N removal is achieved in aerobic downflow columns, 60% NO3--N removal occurs in anaerobic upflow columns, and 60% removal of total nitrogen can be achieved in downflow-upflow in-series. The experimental results are studied further using a variably saturated flow and reactive transport model, which allows a mechanistic explanation of the fate and transport of oxygen and nitrogen. The model clarifies the mechanisms of oxygen transport and nitrogen consumption, and clarifies the need for readily biodegradable COD for denitrification. A VFCW is then compared to a horizontal flow constructed wetland (HFCW) for life cycle environmental impacts. High areal emissions of greenhouse gases from VFCWs compared to HFCWs are the driver for the study. The assessment shows that because a VFCW is only 25% of the volume of an HFCW designed for the same treatment quality, the VFCW has only 25-30% of HFCW impacts over 12 impact categories and 3 damage categories. Results show that impacts could be reduced by design improvements. Design recommendations are downflow wetlands for nitrification, upflow wetlands for denitrification, series wetlands for total nitrogen removal, hydraulic load of 142 L/m2d, 30 cm downflow wetland depth, 1.0 m upflow wetland depth, recycle, vegetation and medium-grained sand. These improvements will optimize nitrogen removal, minimize gaseous emissions, and reduce wetland material requirements, thus reducing environmental impact without sacrificing wastewater treatment quality.