Concentration polarization model of spiral-wound membrane modules with application to batch-mode RO desalination of brackish water

Batch-mode reverse osmosis (batch-RO) operation is considered a promising desalination method due to its low energy requirement compared to other RO system arrangements. To improve and predict batch-RO performance, studies on concentration polarization (CP) are carried out. The Kimura-Sourirajan mass-transfer model is applied and validated by experimentation with two different spiral-wound RO elements. Explicit analytical Sherwood correlations are derived based on experimental results. For batch-RO operation, a new genetic algorithm method is developed to estimate the Sherwood correlation parameters, taking into account the effects of variation in operating parameters. Analytical procedures are presented, then the mass transfer coefficient models are developed for different operation processes, i.e., batch-RO and continuous RO. The CP related energy loss in batch-RO operation is quantified based on the resulting relationship between feed flow rates and mass transfer coefficients. It is found that CP increases energy consumption in batch-RO by about 25% compared to the ideal case in which CP is absent. For continuous RO process, the derived Sherwood correlation predicted CP accurately. In addition, we determined the optimum feed flow rate of our batch-RO system.

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.

Local shell-side coefficients in shell and tube exchanges the use of a mass transfer technique

A diffusion-controlled electrochemical mass transfer technique has been employed in making local measurements of shell-side coefficients in segmentally baffled shell and tube heat exchangers. Corresponding heat transfer data are predicted through the Chilton and Colburn heat and mass transfer analogy. Mass transfer coefficients were measured for baffle spacing lengths of individual tubes in an internal baffle compartment. Shell-side pressure measurements were also made. Baffle compartment average coefficients derived from individual tube coefficients are shown to be in good agreement with reported experimental bundle average heat transfer data for a heat exchanger model of similar geometry. Mass transfer coefficients of individual tubes compare favourably with those obtained previously by another mass transfer technique. Experimental data are reported for a variety of segmental baffle configurations over the shell-side Reynolds number range 100 to 42 000. Baffles with zero clearances were studied at three baffle cuts and two baffle spacings. Baffle geometry is shown to have a large effect on the distribution of tube coefficients within the baffle compartment. Fluid "jetting" is identified with some baffle configurations. No simple characteristic velocity is found to correlate zonal or baffle compartment average mass transfer data for the effect of both baffle cut and baffle spacing. Experiments with baffle clearances typical of commercial heat exchangers are also reported. The effect of leakage streams associated with these baffles is identified. Investigations were extended to double segmental baffles for which no data had previously been published. The similarity in the shell-side characteristics of this baffle arrangement and two parallel single segmental baffle arrangements is demonstrated. A general relationship between the shell-side mass transfer performance and pressure drop was indicated by the data for all the baffle configurations examined.

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

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

Mass transfer characteristics of two-aqueous-phase liquid-liquid mixtures

Mass transfer rates were studied using the falling drop method. Cibacron Blue 3 GA dye was the transferring solute from the salt phase to the PEG phase. Measurements were undertaken for several concentrations of the dye and the phase-forming solutes and with a range of different drop sizes, e.g. 2.8, 3.0 and 3.7 mm. The dye was observed to be present in the salt phase as finely dispersed solids but a model confirmed that the mass transfer process could still be described by an equation based upon the Whitman two-film model. The overall mass transfer coefficient increased with increasing concentration of the dye. The apparent mass transfer coefficient ranged from 1 x 10-5 to 2 x 10 -4 m/s. Further experiments suggested that mass transfer was enhanced at high concentration by several mechanisms. The dye was found to change the equilibrium composition of the two phases, leading to transfer of salt between the drop and continuous phases. It also lowered the interfacial tension (i.e. from 1.43 x 10-4 N/m for 0.01% w/w dye concentration to 1.07 x 10-4 N/m for 0.2% w/w dye concentration) between the two phases, which could have caused interfacial instabilities (Marangoni effects). The largest drops were deformable, which resulted in a significant increase in the mass transfer rate. Drop size distribution and Sauter mean drop diameter were studied on-line in a 1 litre agitated vessel using a laser diffraction technique. The effects of phase concentration, dispersed phase hold-up and impeller speed were investigated for the salt-PEG system. An increase in agitation speed in the range 300 rpm to 1000 rpm caused a decrease in mean drop diameter, e.g. from 50 m to 15 m. A characteristic bimodal drop size distribution was established within a very short time. An increase in agitation rate caused a shift of the larger drop size peak to a smaller size.

Mechanisms of heat and mass transfer to and from single drops freely-suspended in an air stream

A specially-designed vertical wind tunnel was used to freely suspend individual liquid drops of 5 mm initial diameter to investigate drop dynamics, terminal velocity and heat and mass transfer rates. Droplets of distilled, de-ionised water, n-propanol, iso-butanol, monoethanolamine and heptane were studied over a temperature range of 50oC to 82oC. The effects of substances that may provide drop surface rigidity (e.g. surface active agents, binders and polymers) on mass transfer rates were investigated by doping distilled de-ionised water drops with sodium di-octyl sulfo-succinate surfactant. Mass transfer rates decreased with reduced drop oscillation as a result of surfactant addition, confirming the importance of droplet surface instability. Rigid naphthalene spheres and drops which formed a skin were also studied; the results confirmed the reduced transfer rates in the absence of drop fluidity. Following consideration of fundamental drop dynamics in air and experimental results from this study, a novel dimensionless group, the Oteng-Attakora, (OT), number was included in the mass transfer equation to account for droplet surface behaviour and for prediction of heat and mass transfer rates from single drops which exhibit surface instability at Re>=500. The OT number and the modified mass transfer equation are respectively: OT=(ava2/d).de1.5(d/) Sh = 2 + 0.02OT0.15Re0.88Sc0.33 Under all conditions drop terminal velocity increased linearly with the square root of drop diameter and the drag coefficient was 1. The data were correlated with a modified equation by Finlay as follows: CD=0.237.((Re/P0.13)1.55(1/We.P0.13) The relevance of the new model to practical evaporative spray processes is discussed.

Mass transfer coefficient for drying of moist particulate in a bubbling fluidized bed

Experiments on drying of moist particles by ambient air were carried out to measure the mass transfer coefficient in a bubbling fluidized bed. Fine glass beads of mean diameter 125?µm were used as the bed material. Throughout the drying process, the dynamic material distribution was recorded by electrical capacitance tomography (ECT) and the exit air condition was recorded by a temperature/humidity probe. The ECT data were used to obtain qualitative and quantitative information on the bubble characteristics. The exit air moisture content was used to determine the water content in the bed. The measured overall mass transfer coefficient was in the range of 0.0145–0.021?m/s. A simple model based on the available correlations for bubble-cloud and cloud-dense interchange (two-region model) was used to predict the overall mass transfer coefficient. Comparison between the measured and predicted mass transfer coefficient have shown reasonable agreement. The results were also used to determine the relative importance of the two transfer regions.

Development of novel mass spectrometric methods for identifying HOCl-induced modifications to proteins

Protein oxidation is thought to contribute to a number of inflammatory diseases, hence the development of sensitive and specific analytical techniques to detect oxidative PTMs (oxPTMs) in biological samples is highly desirable. Precursor ion scanning for fragment ions of oxidized amino acid residues was investigated as a label-free MS approach to mapping specific oxPTMs in a complex mixture of proteins. Using HOCl-oxidized lysozyme as a model system, it was found that the immonium ions of oxidized tyrosine and tryptophan formed in MS(2) analysis could not be used as diagnostic ions, owing to the occurrence of isobaric fragment ions from unmodified peptides. Using a double quadrupole linear ion trap mass spectrometer, precursor ion scanning was combined with detection of MS(3) fragment ions from the immonium ions and collisionally-activated decomposition peptide sequencing to achieve selectivity for the oxPTMs. For chlorotyrosine, the immonium ion at 170.1 m/z fragmented to yield diagnostic ions at 153.1, 134.1, and 125.1 m/z, and the hydroxytyrosine immonium ion at 152.1 m/z gave diagnostic ions at 135.1 and 107.1 m/z. Selective MS(3) fragment ions were also identified for 2-hydroxytryptophan and 5-hydroxytryptophan. The method was used successfully to map these oxPTMs in a mixture of nine proteins that had been treated with HOCl, thereby demonstrating its potential for application to complex biological samples.

A dynamic model for the activated sludge treatment of coke oven effluents

This is an Inter-Disciplinary Higher Degree (IHD) thesis about Water Pollution Control in the Iron and Steel Industry. After examining the compositions, and various treatment methods, for the major effluent streams from a typical Integrated Iron and Steel works, it was decided to concentrate investigative work on the activated-sludge treatment of coke-oven effluents. A mathematical model of this process was developed in an attempt to provide a tool for plant management that would enable improved performance, and enhanced control of Works Units. The model differs from conventional models in that allowance is made for the presence of two genera of microorganisms, each of which utilises a particular type of substrate as its energy source. Allowance is also made for the inhibitive effect of phenol on thiocyanate biodegradation, and for the self-toxicity of the bacteria when present in a high substrate concentration environment. The enumeration of the kinetic characteristics of the two groups of micro-organisms was shown to be of major importance. Laboratory experiments were instigated in an attempt to determine accurate values of these coefficients. The use of the Suspended Solids concentration was found to be too insensitive a measure of viable active mass. Other measures were investigated, and Adenosine Triphosphate concentration was chosen as the most effective measure of bacterial populations. Using this measure, a model was developed for phenol biodegradation from experimental results which implicated the possibility of storage of substate prior to metabolism. A model for thiocyanate biodegradation was also developed, although the experimental results indicate that much work is still required in this area.

Vibration analysis of submerged rectangular microplates with distributed mass loading

This paper investigates the vibration characteristics of the coupling system of a microscale fluid-loaded rectangular isotropic plate attached to a uniformly distributed mass. Previous literature has, respectively, studied the changes in the plate vibration induced by an acoustic field or by the attached mass loading. This paper investigates the issue of involving these two types of loading simultaneously. Based on Lamb's assumption of the fluid-loaded structure and the Rayleigh–Ritz energy method, this paper presents an analytical solution for the natural frequencies and mode shapes of the coupling system. Numerical results for microplates with different types of boundary conditions have also been obtained and compared with experimental and numerical results from previous literature. The theoretical model and novel analytical solution are of particular interest in the design of microplate-based biosensing devices.

Proteomic analysis of a noninvasive human model of acute inflammation and its resolution:the twenty-one day gingivitis model

The 21-day experimental gingivitis model, an established noninvasive model of inflammation in response to increasing bacterial accumulation in humans, is designed to enable the study of both the induction and resolution of inflammation. Here, we have analyzed gingival crevicular fluid, an oral fluid comprising a serum transudate and tissue exudates, by LC-MS/MS using Fourier transform ion cyclotron resonance mass spectrometry and iTRAQ isobaric mass tags, to establish meta-proteomic profiles of inflammation-induced changes in proteins in healthy young volunteers. Across the course of experimentally induced gingivitis, we identified 16 bacterial and 186 human proteins. Although abundances of the bacterial proteins identified did not vary temporally, Fusobacterium outer membrane proteins were detected. Fusobacterium species have previously been associated with periodontal health or disease. The human proteins identified spanned a wide range of compartments (both extracellular and intracellular) and functions, including serum proteins, proteins displaying antibacterial properties, and proteins with functions associated with cellular transcription, DNA binding, the cytoskeleton, cell adhesion, and cilia. PolySNAP3 clustering software was used in a multilayered analytical approach. Clusters of proteins that associated with changes to the clinical parameters included neuronal and synapse associated proteins.

Uncatalysed and potassium-catalysed pyrolysis of the cell-wall constituents of biomass and their model compounds

Cell-wall components (cellulose, hemicellulose (oat spelt xylan), lignin (Organosolv)), and model compounds (levoglucosan (an intermediate product of cellulose decomposition) and chlorogenic acid (structurally similar to lignin polymer units)) have been investigated to probe in detail the influence of potassium on their pyrolysis behaviours as well as their uncatalysed decomposition reaction. Cellulose and lignin were pretreated to remove salts and metals by hydrochloric acid, and this dematerialized sample was impregnated with 1% of potassium as potassium acetate. Levoglucosan, xylan and chlorogenic acid were mixed with CHCOOK to introduce 1% K. Characterisation was performed using thermogravimetric analysis (TGA) and differential thermal analysis (DTA). In addition to the TGA pyrolysis, pyrolysis-gas chromatography-mass spectrometry (PY-GC-MS) analysis was introduced to examine reaction products. Potassium-catalysed pyrolysis has a huge influence on the char formation stage and increases the char yields considerably (from 7.7% for raw cellulose to 27.7% for potassium impregnated cellulose; from 5.7% for raw levoglucosan to 20.8% for levoglucosan with CHCOOK added). Major changes in the pyrolytic decomposition pathways were observed for cellulose, levoglucosan and chlorogenic acid. The results for cellulose and levoglucosan are consistent with a base catalysed route in the presence of the potassium salt which promotes complete decomposition of glucosidic units by a heterolytic mechanism and favours its direct depolymerization and fragmentation to low molecular weight components (e.g. acetic acid, formic acid, glyoxal, hydroxyacetaldehyde and acetol). Base catalysed polymerization reactions increase the char yield. Potassium-catalysed lignin pyrolysis is very significant: the temperature of maximum conversion in pyrolysis shifts to lower temperature by 70 K and catalysed polymerization reactions increase the char yield from 37% to 51%. A similar trend is observed for the model compound, chlorogenic acid. The addition of potassium does not produce a dramatic change in the tar product distribution, although its addition to chlorogenic acid promoted the generation of cyclohexane and phenol derivatives. Postulated thermal decomposition schemes for chlorogenic acid are presented. © 2008 Elsevier B.V. All rights reserved.

Heat and mass transfer in membrane distillation used for desalination with slip flow

A theoretical model for the transport phenomena in an air gap membrane distillation is presented. The model is based on the conservation equations for the mass, momentum, energy and species within the feed water solution as well as on the mass and energy balances on the membrane sides. The slip flow occurs due to the hydrophobic properties of the membrane. The slip boundary condition applied on the feed saline solution-membrane interface is taken into consideration showing its effects on process parameters particularly permeate flow, heat transfer coefficient and thermal efficiency. The theoretical model was validated with available experimental data and was found to be in good agreement especially when the slip condition is introduced. Increasing slip length from zero to 200 μm was found to increase the permeate flux and the thermal efficiency by 33% and 1.7% respectively.

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.

An operational research-based integrated approach for mass evacuation planning of a city

Large-scale disasters are constantly occurring around the world, and in many cases evacuation of regions of city is needed. ‘Operational Research/Management Science’ (OR/MS) has been widely used in emergency planning for over five decades. Warning dissemination, evacuee transportation and shelter management are three ‘Evacuation Support Functions’ (ESF) generic to many hazards. This thesis has adopted a case study approach to illustrate the importance of integrated approach of evacuation planning and particularly the role of OR/MS models. In the warning dissemination phase, uncertainty in the household’s behaviour as ‘warning informants’ has been investigated along with uncertainties in the warning system. An agentbased model (ABM) was developed for ESF-1 with households as agents and ‘warning informants’ behaviour as the agent behaviour. The model was used to study warning dissemination effectiveness under various conditions of the official channel. In the transportation phase, uncertainties in the household’s behaviour such as departure time (a function of ESF-1), means of transport and destination have been. Households could evacuate as pedestrians, using car or evacuation buses. An ABM was developed to study the evacuation performance (measured in evacuation travel time). In this thesis, a holistic approach for planning the public evacuation shelters called ‘Shelter Information Management System’ (SIMS) has been developed. A generic allocation framework of was developed to available shelter capacity to the shelter demand by considering the evacuation travel time. This was formulated using integer programming. In the sheltering phase, the uncertainty in household shelter choices (either nearest/allocated/convenient) has been studied for its impact on allocation policies using sensitivity analyses. Using analyses from the models and detailed examination of household states from ‘warning to safety’, it was found that the three ESFs though sequential in time, however have lot of interdependencies from the perspective of evacuation planning. This thesis has illustrated an OR/MS based integrated approach including and beyond single ESF preparedness. The developed approach will help in understanding the inter-linkages of the three evacuation phases and preparing a multi-agency-based evacuation planning evacuation