Modeling of nonreacting and reacting turbulent spray jets using a fully stochastic separated flow approach

Numerical modeling of several turbulent nonreacting and reacting spray jets is carried out using a fully stochastic separated flow (FSSF) approach. As is widely used, the carrier-phase is considered in an Eulerian framework, while the dispersed phase is tracked in a Lagrangian framework following the stochastic separated flow (SSF) model. Various interactions between the two phases are taken into account by means of two-way coupling. Spray evaporation is described using a thermal model with an infinite conductivity in the liquid phase. The gas-phase turbulence terms are closed using the k-epsilon model. A novel mixture fraction based approach is used to stochastically model the fluctuating temperature and composition in the gas phase and these are then used to refine the estimates of the heat and mass transfer rates between the droplets and the surrounding gas-phase. In classical SSF (CSSF) methods, stochastic fluctuations of only the gas-phase velocity are modeled. Successful implementation of the FSSF approach to turbulent nonreacting and reacting spray jets is demonstrated. Results are compared against experimental measurements as well as with predictions using the CSSF approach for both nonreacting and reacting spray jets. The FSSF approach shows little difference from the CSSF predictions for nonreacting spray jets but differences are significant for reacting spray jets. In general, the FSSF approach gives good predictions of the flame length and structure but further improvements in modeling may be needed to improve the accuracy of some details of the Predictions. (C) 2011 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Mass Distribution Studies in Effervescent Sprays

An experimental study on liquid mass distribution in effervescent sprays using water and air as working fluids is presented in this paper. Optical patternation and techniques of image processing are employed for analyzing the spray. The flow regime inside the effervescent atomizer largely dictates the mass distribution patterns. The patterns are seen to vary from concentrated, poorly atomized liquid lumps to uniformly distributed, fine droplets as the flow regime changes from bubbly flow to annular flow. A large variety of instantaneous spray patterns are observed in bubbly flow regime indicating a highly unsteady atomization process. However, relatively better consistency in spray patterns is observed at higher gas flow rates. Thus, the degree of unsteadiness gradually diminishes as gas flow rate is increased. The axial evolution of the spray in annular mode shows good mixing of liquid and gas across the interface.

I. The effects of variations in certain parameters upon the efficiency of in vivo hemodialysis with a Kiil dialyzer. II. Interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube.

Part I

Regression analyses are performed on in vivo hemodialysis data for the transfer of creatinine, urea, uric acid and inorganic phosphate to determine the effects of variations in certain parameters on the efficiency of dialysis with a Kiil dialyzer. In calculating the mass transfer rates across the membrane, the effects of cell-plasma mass transfer kinetics are considered. The concept of the effective permeability coefficient for the red cell membrane is introduced to account for these effects. A discussion of the consequences of neglecting cell-plasma kinetics, as has been done to date in the literature, is presented.

A physical model for the Kiil dialyzer is presented in order to calculate the available membrane area for mass transfer, the linear blood and dialysate velocities, and other variables. The equations used to determine the independent variables of the regression analyses are presented. The potential dependent variables in the analyses are discussed.

Regression analyses were carried out considering overall mass-transfer coefficients, dialysances, relative dialysances, and relative permeabilities for each substance as the dependent variables. The independent variables were linear blood velocity, linear dialysate velocity, the pressure difference across the membrane, the elapsed time of dialysis, the blood hematocrit, and the arterial plasma concentrations of each substance transferred. The resulting correlations are tabulated, presented graphically, and discussed. The implications of these correlations are discussed from the viewpoint of a research investigator and from the viewpoint of patient treatment.

Recommendations for further experimental work are presented.

Part II

The interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube in the wavy flow regime has been measured using resistance wave gages. The median water depth, r.m.s. wave height, wave frequency, extrema frequency, and wave velocity have been measured as functions of air and water flow rates. Reynolds numbers, Froude numbers, Weber numbers, and bulk velocities for each phase may be calculated from these measurements. No theory for wave formation and propagation available in the literature was sufficient to describe these results.

The water surface level distribution generally is not adequately represented as a stationary Gaussian process. Five types of deviation from the Gaussian process function were noted in this work. The presence of the tube walls and the relatively large interfacial shear stresses precludes the use of simple statistical analyses to describe the interfacial structure. A detailed study of the behavior of individual fluid elements near the interface may be necessary to describe adequately wavy two-phase flow in systems similar to the one used in this work.

Effects of forced solution flow on lysozyme crystal growth

In order to study quantitatively the effects of forced solution on crystal growth, we designed a new set of experimental equipment, in particular, a microchannel mixer was used as crystallization container so that the consumption of protein samples was much reduced and thus an exact syringe pump could be used for precise control of the flow rates. Since the mixer’s section was designed to be rectangular, the solution velocity in its center was steady and constant, and thus repeatable experiments were facilitated. Experimental results showed that the effects of forced solution on protein crystal growth were different under different levels of supersaturation, and new results were obtained for cases of high supersaturation. When the supersaturation is σ = 2.3, with increasing flow rates the growth rates of the lysozyme crystal’s (110) face hardly change when the flow rates are lower than 1300 μm/s, and decrease quickly afterwards. When the flow rate reaches 2000 μm/s, the crystal nearly ceases to grow. When the supersaturation is σ = 2.7, with increasing flow rates the (110) face growth rates increase at the beginning then reach the maximum values at 1700 μm/s – 1900 μm/s and decrease afterwards, approaching zero or so when the flow rate reaches 12000 μm/s. The higher the supersaturation, the larger the flow rate at which the crystal ceases to grow. © 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Advanced polymer membrane development in pervaporation dehydration and lateral flow diagnostics

The work in this thesis concerns the advanced development of polymeric membranes of two types; pervaporation and lateral-flow. The former produced from a solution casting method and the latter from a phase separation. All membranes were produced from casting lacquers. Early research centred on the development of viable membranes. This led to a supported polymer blend pervaporation membrane. Selective layer: plasticized 4:1 mass ratio sodium-alginate: poly(vinyl-alcohol) polymer blend. Using this membrane, pervaporation separation of ethanol/water mixtures was carefully monitored as a function of film thickness and time. Contrary to literature expectations, these films showed increased selectivity and decreased flux as film thickness was reduced. It is argued that morphology and structure of the polymer blend changes with thickness and that these changes define membrane efficiency. Mixed matrix membrane development was done using spherical, discreet, size-monodisperse mesoporous silica particles of 1.8 - 2μm diameter, with pore diameters of ~1.8 nm were incorporated into a poly(vinyl alcohol) [PVA] matrix. Inclusion of silica benefitted pervaporation performance for the dehydration of ethanol, improving flux and selectivity throughout in all but the highest silica content samples. Early lateral-flow membrane research produced a membrane from a basic lacquer composition required for phase inversion; polymer, solvent and non-solvent. Results showed that bringing lacquers to cloud point benefits both the pore structure and skin layers of the membranes. Advancement of this work showed that incorporation of ethanol as a mesosolvent into the lacquer effectively enhances membrane pore structure resulting in an improvement in lateral flow rates of the final membranes. This project details the formation mechanics of pervaporation and lateral-flow membranes and how these can be controlled. The principle methods of control can be applied to the formation of any other flat sheet polymer membranes, opening many avenues of future membrane research and industrial application.

Review of models for predicting the discharge rates of bulk particulates from silos and bins

Various models for predicting discharge rates have been developed over the last four decades by many research workers (notably Beverloo [1], Johanson [2], Brown [3], Carleton [4], Crewdson [5], Nedderman [6], Gu [7].). In many cases these models offer comparable approaches to the prediction of discharge rates of bulk particulates from storage equipment when solely gravity is acting to initiate flow (since they invariably consider the use of mass-flow design equipment). The models that have been developed consider a wide range of bulk particulates (coarse, incompressible, fine, cohesive) and most contemporary works have incorporated validation against test programmes. Research currently underway at The Wolfson Centre for Bulk Solids Handling Technology, University of Greenwich, has considered the relative performance of these models with respect to a range of bulk properties and with particular focus upon the flexibility of the models to cater for different geometrical factors for vessels.

Solventless continuous flow homogeneous hydroformylation of 1-octene

The hydroformylation of 1-octene under continuous flow conditions is described. The system involves dissolving the catalyst, made in situ from [ Rh(acac)(CO)(2)] (acacH = 2,4- pentanedione) and [RMIM][TPPMS] ( RMIM = 1-propyl (Pr), 1-pentyl (Pn) or 1-octyl (O)-3-methyl imidazolium, TPPMS = Ph2P(3-C6H4SO3)), in a mixture of nonanal and 1-octene and passing the substrate, 1-octene, together with CO and H-2 through the system dissolved in supercritical CO2 (scCO(2)). [PrMIM][TPPMS] is poorly soluble in the medium so heavy rhodium leaching (as complexes not containing phosphine) occurs in the early part of the reaction. [PnMIM][ PPMS] affords good rates at relatively low catalyst loadings and relatively low overall pressure (125 bar) with rhodium losses <1 ppm, but the catalyst precipitates at higher catalyst loadings, leading to lower reaction rates. [OMIM][ TPPMS] is the most soluble ligand and promotes high reaction rates, although preliminary experiments suggested that rhodium leaching was high at 5-10 ppm. Optimisation aimed at balancing flows so that the level within the reactor remained constant involved a reactor set up based around a reactor fitted with a sight glass and sparging stirrer with the CO2 being fed by a cooled head HPLC pump, 1-octene by a standard HPLC pump and CO/H-2 through a mass flow controller. The pressure was controlled by a back pressure regulator. Using this set up, [OMIM][ TPPMS] as the ligand and a total pressure of 140 bar, it was possible to control the level within the reactor and obtain a turnover frequency of ca. 180 h(-1). Rhodium losses in the optimised system were 100 ppb. Transport studies showed that 1-octene is preferentially transported over the aldehydes at all pressures, although the difference in mol fraction in the mobile phase was less at lower pressures. Nonanal in the mobile phase suppresses the extraction of 1-octene to some extent, so it is better to operate at high conversion and low pressure to optimise the extraction of the products relative to the substrate. CO and H2 in the mobile phase also suppress the extraction effciency by as much as 80%.

The empirical metallicity dependence of the mass-loss rate of O- and early B-type stars

We present a comprehensive study of the observational dependence of the mass-loss rate in stationary stellar winds of hot massive stars on the metal content of their atmospheres. The metal content of stars in the Magellanic Clouds is discussed, and a critical assessment is given of state-of-the-art mass-loss determinations of OB stars in these two satellite systems and the Milky-Way. Assuming a power-law dependence of mass loss on metal content,. M. Z(m), and adopting a theoretical relation between the terminal flow velocity and metal content, v(infinity). Z(0.13) (Leitherer et al. 1992, ApJ, 401, 596), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 105.2 L circle dot. Within the errors, this result is in agreement with the prediction m = 0.69 +/- 0.10 by Vink et al. (2001, A& A, 369, 574). Absolute empirical values for the mass loss, based on Ha and ultraviolet (UV) wind lines, are found to be a factor of two higher than predictions in this high luminosity regime. If this difference is attributed to inhomogeneities in the wind, and this clumping does not impact the predictions, this would imply that luminous O and early-B stars have clumping factors in their Ha and UV line forming regions of about a factor of four. For lower luminosity stars, the winds are so weak that their strengths can generally no longer be derived from optical spectral lines (essentially Ha) and one must currently rely on the analysis of UV lines. We confirm that in this low-luminosity domain the observed Galactic WLR is found to be much steeper than expected from theory (although the specific sample is rather small), leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L similar to 10(4.75) L circle dot, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars would turn out to be overestimated as a result of wind clumping, but the degree of clumping would be rather independent of metallicity, the scalings derived in this study are expected to remain correct.

The contribution of the Finite Volume Point Dilution Method to the determination of solute mass flux in an alluvial aquifer.

Groundwater represents one of the most important resources of the world and it is essential to prevent its pollution and to consider remediation intervention in case of contamination. According to the scientific community the characterization and the management of the contaminated sites have to be performed in terms of contaminant fluxes and considering their spatial and temporal evolution. One of the most suitable approach to determine the spatial distribution of pollutant and to quantify contaminant fluxes in groundwater is using control panels. The determination of contaminant mass flux, requires measurement of contaminant concentration in the moving phase (water) and velocity/flux of the groundwater. In this Master Thesis a new solute flux mass measurement approach, based on an integrated control panel type methodology combined with the Finite Volume Point Dilution Method (FVPDM), for the monitoring of transient groundwater fluxes, is proposed. Moreover a new adsorption passive sampler, which allow to capture the variation of solute concentration with time, is designed. The present work contributes to the development of this approach on three key points. First, the ability of the FVPDM to monitor transient groundwater fluxes was verified during a step drawdown test at the experimental site of Hermalle Sous Argentau (Belgium). The results showed that this method can be used, with optimal results, to follow transient groundwater fluxes. Moreover, it resulted that performing FVPDM, in several piezometers, during a pumping test allows to determine the different flow rates and flow regimes that can occurs in the various parts of an aquifer. The second field test aiming to determine the representativity of a control panel for measuring mass flus in groundwater underlined that wrong evaluations of Darcy fluxes and discharge surfaces can determine an incorrect estimation of mass fluxes and that this technique has to be used with precaution. Thus, a detailed geological and hydrogeological characterization must be conducted, before applying this technique. Finally, the third outcome of this work concerned laboratory experiments. The test conducted on several type of adsorption material (Oasis HLB cartridge, TDS-ORGANOSORB 10 and TDS-ORGANOSORB 10-AA), in order to determine the optimum medium to dimension the passive sampler, highlighted the necessity to find a material with a reversible adsorption tendency to completely satisfy the request of the new passive sampling technique.

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.

Computer simulation of liquid flow patterns on distillation trays

This thesis describes work carried out to improve the fundamental modelling of liquid flows on distillation trays. A mathematical model is presented based on the principles of computerised fluid dynamics. It models the liquid flow in the horizontal directions allowing for the effects of the vapour through the use of an increased liquid turbulence, modelled by an eddy viscosity, and a resistance to liquid flow caused by the vapour being accelerated horizontally by the liquid. The resultant equations are similar to the Navier-Stokes equations with the addition of a resistance term.A mass-transfer model is used to calculate liquid concentration profiles and tray efficiencies. A heat and mass transfer analogy is used to compare theoretical concentration profiles to experimental water-cooling data obtained from a 2.44 metre diameter air-water distillation simulation rig. The ratios of air to water flow rates are varied in order to simulate three pressures: vacuum, atmospheric pressure and moderate pressure.For simulated atmospheric and moderate pressure distillation, the fluid mechanical model constantly over-predicts tray efficiencies with an accuracy of between +1.7% and +11.3%. This compares to -1.8% to -10.9% for the stagnant regions model (Porter et al. 1972) and +12.8% to +34.7% for the plug flow plus back-mixing model (Gerster et al. 1958). The model fails to predict the flow patterns and tray efficiencies for vacuum simulation due to the change in the mechanism of liquid transport, from a liquid continuous layer to a spray as the liquid flow-rate is reduced. This spray is not taken into account in the development of the fluid mechanical model. A sensitivity analysis carried out has shown that the fluid mechanical model is relatively insensitive to the prediction of the average height of clear liquid, and a reduction in the resistance term results in a slight loss of tray efficiency. But these effects are not great. The model is quite sensitive to the prediction of the eddy viscosity term. Variations can produce up to a 15% decrease in tray efficiency. The fluid mechanical model has been incorporated into a column model so that statistical optimisation techniques can be employed to fit a theoretical column concentration profile to experimental data. Through the use of this work mass-transfer data can be obtained.

Flow patterns on distillation trays

Studies into the two-phase flow patterns produced on a sieve tray were carried out using an air-water simulator of 2.44 m in diameter. The flow patterns were investigated by a number of methods, direct observation using directional flow pointers; by water-cooling to simulate mass transfer; and by measurement of the height of clear liquid across the tray with manometers. The flow rates used were designed to show how the flow pattern changed with the change in the gas and liquid rates. The results from water-only studies on an un-perforated tray were compared with those produced on a sieve tray with holes of 12.7 mm diameter. The presence of regions on the sides of the tray where the liquid was circulating was noted from the water-only experiments. The presence and magnitude of the circulations was reduced when the air was passed through the liquid. These were similar to the findings of Hine (1990) and Chambers (1993). When circulation occurred, the flow separated at the ends of the inlet downcomer and circulations of up to 30% of the tray area were observed. Water-cooling and the manometer measurements were used to show the effect of the flow pattern on the tray efficiency and the height of clear liquid respectively. The efficiency was severely reduced by the presence of circulations. The height of clear liquid tended to rise in these areas. A comparison of data collected on trays with different hole diameters showed that the larger hole diameter inhibited the on-set of separation to a greater extent than small hole diameters. The tray efficiency was affected by a combination of the better mixing on smaller hole trays and detrimental effect of greater circulation on these trays. Work on a rectangular tray geometry was carried out to assess the effect of hole size on the height of clear liquid. It was found that the gradient on the outlet half of the tray was very small and that the highest clear liquid height was given by the highest hole size. Overall, the experiments helped to clarify the effect that the flow pattern had on the operation of the tray. It is hoped that the work can be of use in the development of models to predict the flow pattern and hence the tray efficiency.

The use of flow control devices to improve the flow pattern and throughput of sieve trays

Compared to packings trays are more cost effective column internals because they create a large interfacial area for mass transfer by the interaction of the vapour on the liquid. The tray supports a mass of froth or spray which on most trays (including the most widely used sieve trays) is not in any way controlled. The two important results of the gas/liquid interaction are the tray efficiency and the tray throughput or capacity. After many years of practical experience, both may be predicted by empirical correlations, despite the lack of understanding. It is known that the tray efficiency is in part determined by the liquid flow pattern and the throughput by the liquid froth height which in turn depends on the liquid hold-up and vapour velocity. This thesis describes experimental work on sieve trays in an air-water simulator, 2.44 m in diameter. The liquid flow pattern, for flow rates similar to those used in commercial scale distillation, was observed experimentally by direct observation; by water-cooling, to simulate mass transfer; use of potassium permanganate dye to observe areas of longer residence time; and by height of clear liquid measurements across the tray and in the downcomer using manometers. This work presents experiments designed to evaluate flow control devices proposed to improve the gas liquid interaction and hence improve the tray efficiency and throughput. These are (a) the use of intermediate weirs to redirect liquid to the sides of the tray so as to remove slow moving/stagnant liquid and (b) the use of vapour-directing slots designed to use the vapour to cause liquid to be directed towards the outlet weir thus reducing the liquid hold-up at a given rate i.e. increased throughput. This method also has the advantage of removing slow moving/stagnant liquid. In the experiments using intermediate weirs, which were placed in the centre of the tray. it was found that in general the effect of an intermediate weir depends on the depth of liquid downstream of the weir. If the weir is deeper than the downstream depth it will cause the upstream liquid to be deeper than the downstream liquid. If the weir is not as deep as deep as the downstream depth it may have little or no effect on the upstream depth. An intermediate weir placed at an angle to the direction of flow of liquid increases the liquid towards the sides of the tray without causing an increase in liquid hold-up/ froth height. The maximum proportion of liquid caused to flow sideways by the weir is between 5% and 10%. Experimental work using vapour-directing slots on a rectangular sieve tray has shown that the horizontal momentum that is imparted to the liquid is dependent upon the size of the slot. If too much momentum is transferred to the liquid it causes hydraulic jumps to occur at the mouth of the slot coupled with liquid being entrained, The use of slots also helps to eliminate the hydraulic gradient across sieve trays and provides a more uniform froth height on the tray. By comparing the results obtained of the tray and point efficiencies, it is shown that a slotted tray reduces both values by approximately 10%. This reduction is due to the fact that with a slotted tray the liquid has a reduced residence time Ion the tray coupled also with the fact that large size bubbles are passing through the slots. The effectiveness of using vapour-directing slots on a full circular tray was investigated by using dye to completely colour the biphase. The removal of the dye by clear liquid entering the tray was monitored using an overhead camera. Results obtained show that the slots are successful in their aim of reducing slow moving liquid from the sides of the tray, The net effect of this is an increase in tray efficiency. Measurements of slot vapour-velocity found it to be approximately equal to the hole velocity.

Flow patterns on sieve trays

Studies into gas-liquid flow patterns were carried out on commercial scale sieve trays where the ratio of froth depth to flow path length is typical of that found in practice. Experiments were conducted on a 2.44 m diameter air-water distillation simulator, in which flow patterns were investigated by direct observation, using directional flow pointers; by water cooling, to simulate mass transfer; and by height of clear liquid measurements across the tray. The flow rates used are typical of those found in practice. The approach adopted was to investigate the effect of the gas flow on the liquid flow by comparing water only flow patterns across an unperforated tray with air-water flow patterns on perforated trays. Initial gas-liquid contacting experiments on the 6.35 mm hole tray showed that, under certain conditions, the gas flow pattern beneath the test tray can have a significant effect on the tray liquid flow pattern such that gas-driven liquid circulation was produced. This was found to be a function of this particular air-water simulator design, and as far as is known this is the first time that this phenomenon has been observed. Consequently non-uniform gas flow effects were removed by modification of the gas distribution system. By eliminating gas circulation effects, the effect of the gas flow on the separation of liquid flow was similar to that obtained on the 1.0 mm hole tray (Hine, 1990). That is, flow separation occurred at the ends of the inlet downcomer which produced large circulating zones along the tray segments both on the non-perforated and perforated trays. The air when forced through the liquid, inhibited circulating flow such that it only occurred at high water inlet velocities. With the 6.35 mm hole tray, the growth and velocity of circulating flow was reduced at high superficial air velocities, and in the experiments to simulate distillation, liquid was in forward flow over most of the tray.

Hydrodynamics and mass transfer studies in a scheibel extractor

A 10 cm diameter four-stage Scheibel column with dispersed phase wetted packing sections has been constructed to study the hydrodynamics and mass transfer using the system toluene-acetone-water. The literature pertaining to the above extractor has been examined and the important phenomena such as droplet break-up and coalescence, mass transfer and backmixing have been reviewed. A critical analysis of the backmixing or axial mixing models and the corresponding techniques for parameter estimation was applied and an optimization technique based on Marquardt's algorithm was implemented. A single phase sampling technique was developed to estimate the acetone concentration profile in both phases along the column. Column flooding characteristics were investigated under various operating conditions and it was found that, when the impellers were located at about DI/5cm from the upper surface of the pads, the limiting flow rates increased with impeller speed. This unusual behaviour was explained in terms of the pumping effect created by the turbine impellers. Correlations were developed to predict Sauter mean drop diameters. A five-cell with backflow model was used to estimate the column performance (stage efficiency) and phases non-ideality (backflow parameters). Overall mass transfer coefficients were computed using the above model and compared with those calculated using the correlations based on single drop mechanism.