Flow patterns and performance of distillation trays

It is known that distillation tray efficiency depends on the liquid flow pattern, particularly for large diameter trays. Scale·up failures due to liquid channelling have occurred, and it is known that fitting flow control devices to trays sometirr.es improves tray efficiency. Several theoretical models which explain these observations have been published. Further progress in understanding is at present blocked by lack of experimental measurements of the pattern of liquid concentration over the tray. Flow pattern effects are expected to be significant only on commercial size trays of a large diameter and the lack of data is a result of the costs, risks and difficulty of making these measurements on full scale production columns. This work presents a new experiment which simulates distillation by water cooling. and provides a means of testing commercial size trays in the laboratory. Hot water is fed on to the tray and cooled by air forced through the perforations. The analogy between heat and mass transfer shows that the water temperature at any point is analogous to liquid concentration and the enthalpy of the air is analogous to vapour concentration. The effect of the liquid flow pattern on mass transfer is revealed by the temperature field on the tray. The experiment was implemented and evaluated in a column of 1.2 m. dia. The water temperatures were measured by thennocouples interfaced to an electronic computerised data logging system. The "best surface" through the experimental temperature measurements was obtained by the mathematical technique of B. splines, and presented in tenos of lines of constant temperature. The results revealed that in general liquid channelling is more imponant in the bubbly "mixed" regime than in the spray regime. However, it was observed that severe channelling also occurred for intense spray at incipient flood conditions. This is an unexpected result. A computer program was written to calculate point efficiency as well as tray efficiency, and the results were compared with distillation efficiencies for similar loadings. The theoretical model of Porter and Lockett for predicting distillation was modified to predict water cooling and the theoretical predictions were shown to be similar to the experimental temperature profiles. A comparison of the repeatability of the experiments with an errors analysis revealed that accurate tray efficiency measurements require temperature measurements to better than ± 0.1 °c which is achievable with conventional techniques. This was not achieved in this work, and resulted in considerable scatter in the efficiency results. Nevertheless it is concluded that the new experiment is a valuable tool for investigating the effect of the liquid flow pattern on tray mass transfer.

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.

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.

Effect of liquid flow patterns on distillation trays

The thesis describes experimental work on sieve trays in an air-water simulator, 2.44 m in diameter. The liquid flow pattern, for flowrates similar to those used in commercial scale distillation, was observed experimentally by water cooling experiments, in which the temperature of the water is measured at over 100 positions over the tray area. The water is cooled by the rising air which is forced through the tray. A heat and mass transfer analogy is drawn whereby the water temperature is mapped to liquid concentration in mass transfer, and the water temperature profiles reveal how liquid channelling may reduce the tray efficiency. The first experiment was to observe the flow of water only over an unperforated tray. With the exception of very low weir loads, the flow separated at the ends of the inlet downcomer. This caused liquid to flow straight across the tray between the downcomers and large circulating regions to be formed in the side regions of the tray. The effect of the air crossflow on the flow pattern was then observed on a sieve tray of 10% free area with 1 mm diameter holes (such as is used in cryogenic distillation). The flow patterns developed on the tray were similar to those produced with water only on the unperforated tray, but at low weir loads the air crossflow prevented separation of the water flow and the associated circulating regions. At higher weir loads, liquid channelling down the centre of the tray and circulation in the side regions occurred. The percentage of the tray occupied by circulating liquid depended upon the velocity of the liquid entering the tray, which was set by the weir load and size of the gap under the inlet downcomer. The water cooling experiments showed that the temperature of the water in a circulating region is much lower than in other parts of the tray, indicating that the driving force for heat transfer is reduced. In a column section where trays (and circulating areas) are mounted on top of each other, the circulating regions will cause air (or vapour) passing through them to have a reduced change in temperature or concentration leading a loss in tray efficiency.

Flow patterns, performance and scale-up of distillation trays

This work is concerned with the nature of liquid flow across industrial sieve trays operating in the spray, mixed, and the emulsified flow regimes. In order to overcome the practical difficulties of removing many samples from a commercial tray, the mass transfer process was investigated in an air water simulator column by heat transfer analogy. The temperature of the warm water was measured by many thermocouples as the water flowed across the single pass 1.2 m diameter sieve tray. The thermocouples were linked to a mini computer for the storage of the data. The temperature data were then transferred to a main frame computer to generate temperature profiles - analogous to concentration profiles. A comprehensive study of the existing tray efficiency models was carried out using computerised numerical solutions. The calculated results were compared with experimental results published by the Fractionation Research Incorporation (FRl) and the existing models did not show any agreement with the experimental results. Only the Porter and Lockett model showed a reasonable agreement with the experimental results for cenain tray efficiency values. A rectangular active section tray was constructed and tested to establish the channelling effect and the result of its effect on circular tray designs. The developed flow patterns showed predominantly flat profiles and some indication of significant liquid flow through the central region of the tray. This comfirms that the rectangular tray configuration might not be a satisfactory solution for liquid maldistribution on sieve trays. For a typical industrial tray the flow of liquid as it crosses the tray from the inlet to the outlet weir could be affected by the mixing of liquid by the eddy, momentum and the weir shape in the axial or the transverse direction or both. Conventional U-shape profiles were developed when the operating conditions were such that the froth dispersion was in the mixed regime, with good liquid temperature distribution while in the spray regime. For the 12.5 mm hole diameter tray the constant temperature profiles were found to be in the axial direction while in the spray regime and in the transverse direction for the 4.5 mm hole tray. It was observed that the extent of the liquid stagnant zones at the sides of the tray depended on the tray hole diameter and was larger for the 4.5 mm hole tray. The liquid hold-up results show a high liquid hold-up at the areas of the tray with low liquid temperatures, this supports the doubts about the assumptions of constant point efficiency across an operating tray. Liquid flow over the outlet weir showed more liquid flow at the centre of the tray at high liquid loading with low liquid flow at both ends of the weir. The calculated results of the point and tray efficiency model showed a general increase in the calculated point and tray efficiencies with an increase in the weir loading, as the flow regime changed from the spray to the mixed regime the point and the tray efficiencies increased from approximately 30 to 80%.Through the mixed flow regime the efficiencies were found to remain fairly constant, and as the operating conditions were changed to maintain an emulsified flow regime there was a decrease in the resulting efficiencies. The results of the estimated coefficient of mixing for the small and large hole diameter trays show that the extent of liquid mixing on an operating tray generally increased with increasing capacity factor, but decreased with increasing weir loads. This demonstrates that above certain weir loads, the effect of eddy diffusion mechanism on the process of liquid mixing on an operating tray to be negligible.

Groundwater flow patterns within the Triassic sandstone aquifers, Birmingham, United Kingdom

This paper presents the development of a modelling study for part of the Birmingham area. Restricted access and model resolutions have limited wide applications of some of the previously developed models. The study area covers approximately 221 km2, and is underlain geologically, by a multi-layer setup with varied hydraulic properties. The basal aquifer unit is the Kidderminster sandstone Formation, overlain by the Wildmoor and Bromsgrove sandstone Formations. The presence of the Birmingham fault which acts as low permeability barrier demarcates the eastern and southern boundaries. The western boundary is defined by the presence of crystallised rocks and coal measures, while a groundwater divide defines the northern boundary. The estimated recharge flux is 112 mm/yr. The ranges of calibrated values obtained for horizontal and vertical hydraulic conductivities are 5.787x10-6 - 2.315x10-5  m/s and 5.787x10-8  - 1.157x10-7  m/s, respectively. Corresponding values obtained for the specific yield and specific storage are 0.10 - 0.12, and 1x10 -4 - 5x10 -4. The calculated numerical error is generally much less than 0.1 %. Hydraulic layering within the Permo-Triassic sandstone aquifer is thought to account for the large vertical anisotropy. Although, uncertainties are associated with the use of a simplistic delay approach to characterise the effects of the unsaturated zone, the modelled values are comparable with those obtained in the literature, and the flow pattern predictions appear to be realistic. © Research India Publications.

Variations of local heat transfer coefficient in piped flow of viscous liquids

Measurements were carried out to determine local coefficients of heat transfer in short lengths of horizontal pipe, and in the region of an discontinuity in pipe diameter. Laminar, transitional and turbulent flow regimes were investigated, and mixtures of propylene glycol and water were used in the experiments to give a range of viscous fluids. Theoretical and empirical analyses were implemented to find how the fundamental mechanism of forced convection was modified by the secondary effects of free convection, temperature dependent viscosity, and viscous dissipation. From experiments with the short tube it was possible to determine simple empirical relationships describing the axial distribution of the local 1usselt number and its dependence on the Reynolds and Prandtl numbers. Small corrections were made to account for the secondary effects mentioned above. Two different entrance configurations were investigated to demonstrate how conditions upstream could influence the heat transfer coefficients measured downstream In experiments with a sudden contraction in pipe diameter the distribution of local 1u3se1t number depended on the Prandtl number of the fluid in a complicated way. Graphical data is presented describing this dependence for a range of fluids indicating how the local Nusselt number varied with the diameter-ratio. Ratios up to 3.34:1 were considered. With a sudden divergence in pipe diameter, it was possible to derive the axial distribution of the local Nusse1t number for a range of Reynolds and Prandtl numbers in a similar way to the convergence experiments. Difficulty was encountered in explaining some of the measurements obtained at low Reynolds numbers, and flow visualization techniques wore used to determine the complex flow patterns which could lead to the anomalous results mentioned. Tests were carried out with divergences up to 1:3.34 to find the way in which the local Nusselt number varied with the diameter ratio, and a few experiments were carried out with very large ratios up .to 14.4. A limited amount of theoretical analysis of the 'divergence' system was carried out to substantiate certain explanations of the heat transfer mechanisms postulated.

Flow field analysis of some mixing and conveying screw element regions, within a closely intermeshing, co-rotating twin-screw extruder.

Grafting of antioxidants and other modifiers onto polymers by reactive extrusion, has been performed successfully by the Polymer Processing and Performance Group at Aston University. Traditionally the optimum conditions for the grafting process have been established within a Brabender internal mixer. Transfer of this batch process to a continuous processor, such as an extruder, has, typically, been empirical. To have more confidence in the success of direct transfer of the process requires knowledge of, and comparison between, residence times, mixing intensities, shear rates and flow regimes in the internal mixer and in the continuous processor.The continuous processor chosen for the current work in the closely intermeshing, co-rotating twin-screw extruder (CICo-TSE). CICo-TSEs contain screw elements that convey material with a self-wiping action and are widely used for polymer compounding and blending. Of the different mixing modules contained within the CICo-TSE, the trilobal elements, which impose intensive mixing, and the mixing discs, which impose extensive mixing, are of importance when establishing the intensity of mixing. In this thesis, the flow patterns within the various regions of the single-flighted conveying screw elements and within both the trilobal element and mixing disc zones of a Betol BTS40 CICo-TSE, have been modelled using the computational fluid dynamics package Polyflow. A major obstacle encountered when solving the flow problem within all of these sets of elements, arises from both the complex geometry and the time-dependent flow boundaries as the elements rotate about their fixed axes. Simulation of the time dependent boundaries was overcome by selecting a number of sequential 2D and 3D geometries, used to represent partial mixing cycles. The flow fields were simulated using the ideal rheological properties of polypropylene and characterised in terms of velocity vectors, shear stresses generated and a parameter known as the mixing efficiency. The majority of the large 3D simulations were performed on the Cray J90 supercomputer situated at the Rutherford-Appleton laboratories, with pre- and postprocessing operations achieved via a Silicon Graphics Indy workstation. A mechanical model was constructed consisting of various CICo-TSE elements rotating within a transparent outer barrel. A technique has been developed using coloured viscous clays whereby the flow patterns and mixing characteristics within the CICo-TSE may be visualised. In order to test and verify the simulated predictions, the patterns observed within the mechanical model were compared with the flow patterns predicted by the computational model. The flow patterns within the single-flighted conveying screw elements in particular, showed good agreement between the experimental and simulated results.

Spray drier simulation and air flow pattern studies

The literature pertaining to the key stages of spray drying has been reviewed in the context of the mathematical modelling of drier performance. A critical review is also presented of previous spray drying models. A new mathematical model has been developed for prediction of spray drier performance. This is applicable to slurries of rigid, porous crust-forming materials to predict trajectories and drying profiles for droplets with a distribution of sizes sprayed from a centrifugal pressure nozzle. The model has been validated by comparing model predictions to experimental data from a pilot-scale counter-current drier and from a full-scale co-current drier. For the latter, the computed product moisture content was within 2%, and the computed air exit temperature within 10oC of experimental data. Air flow patterns have been investigated in a 1.2m diameter transparent countercurrent spray tower by flow visualisation. Smoke was introduced into various zones within the tower to trace the direction, and gauge the intensity, of the air flow. By means of a set of variable-angle air inlet nozzles, a variety of air entry configurations was investigated. The existence of a core of high rotational and axial velocity channelling up the axis of the tower was confirmed. The stability of flow within the core was found to be strongly dependent upon the air entry arrangement. A probe was developed for the measurement of air temperature and humidity profiles. This was employed for studying evaporation of pure water drops in a 1.2m diameter pilot-scale counter-current drier. A rapid approach to the exit air properties was detected within a 1m distance from the air entry ports. Measured radial profiles were found to be virtually flat but, from the axial profiles, the existence of plug-flow, well-mixed-flow and some degree of air short-circuiting can be inferred. The model and conclusions should assist in the improved design and optimum operation of industrial spray driers.

Time-dependent evolution of tissue markers by MALDI-MS imaging

We have used MALDI-MS imaging (MALDI-MSI) to monitor the time dependent appearance and loss of signals when tissue slices are brought rapidly to room temperature for short to medium periods of time. Sections from mouse brain were cut in a cryostat microtome, placed on a MALDI target and allowed to warm to room temperature for 30 s to 3 h. Sections were then refrozen, fixed by ethanol treatment and analysed by MALDI-MSI. The intensity of a range of markers were seen to vary across the time course, both increasing and decreasing, with the intensity of some markers changing significantly within 30 s and markers also showed tissue location specific evolution. The markers resulting from this autolysis were compared directly to those that evolved in a comparable 16 h on-tissue trypsin digest, and the markers that evolved in the two studies were seen to be substantially different. These changes offer an important additional level of location-dependent information for mapping changes and seeking disease-dependent biomarkers in the tissue. They also indicate that considerable care is required to allow comparison of biomarkers between MALDI-MSI experiments and also has implications for the standard practice of thaw-mounting multiple tissue sections onto MALDI-MS targets.

Further studies on flow distribution on distillation sieve trays

The thesis describes experimental work on the possibility of using deflection baffles in conventional distillation trays as flow straightening devices, with the view of enhancing tray efficiency. The mode of operation is based on deflecting part of the liquid momentum from the centre of the tray to the segment regions in order to drive stagnating liquid at the edges forward. The first part of the work was a detailed investigation into the two-phase flow patterns produced on a conventional sieve tray having 1 mm hole size perforations. The data provide a check on some earlier work and extend the range of the existing databank, particularly to conditions more typical of industrial operation. A critical survey of data collected on trays with different hole sizes (Hine, 1990; Chambers, 1993; Fenwick, 1996; this work) showed that the hole diameter has a significant influence on the flow regime, the size of the stagnant regions and the hydraulic and mass transfer performance. Five modified tray topologies were created with different configurations of baffles and tested extensively in the 2.44 m diameter air-water pilot distillation simulator for their efficacy in achieving uniform flow across the tray and for their impact on tray loading capacity and mass transfer efficiency. Special attention was given to the calibration of the over 100 temperature probes used in measuring the water temperature across the tray on which the heat and mass transfer analogy is based. In addition to normal tray capacity experiments, higher weir load experiments were conducted using a 'half-tray' mode in order to extend the range of data to conditions more typical of industrial operation. The modified trays show superior flow characteristics compared to the conventional tray in terms of the ability to replenish the zones of exceptionally low temperatures and high residence times at the edges of the tray, to lower the bulk liquid gradient and to achieve a more uniform flow across the tray. These superior flow abilities, however, tend to diminish with increasing weir load because of the increasing tendency for the liquid to jump over the barriers instead of flowing over them. The modified tray topologies showed no tendency to cause undue limitation to tray loading capacity. Although the improvement in the efficiency of a single tray over that of the conventional tray was moderate and in some cases marginal, the multiplier effect in a multiple tray column situation would be significant (Porter et al., 1972). These results are in good agreement with an associated CFD studies (Fischer, 1999) carried out by partners in the Advanced Studies in Distillation consortium. It is concluded that deflection baffles can be used in a conventional distillation sieve tray to achieve better liquid flow distribution and obtain enhanced mass transfer efficiency, without undermining the tray loading capacity. Unlike any other controlled-flow tray whose mechanical complexity impose stringent manufacturing and installation tolerances, the baffled-tray models are simple to design, manufacture and install and thus provide an economic method of retrofitting badly performing sieve trays both in terms of downtime and fabrication. NOTE APPENDICES 2-5 ARE ON A SEPARATE FLOPPY DISK ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY WITH PRIOR ARRANGEMENT

Flow and pressure drop on the shellside of cylindrical heat exchangers

This work presents pressure distributions and fluid flow patterns on the shellside of a cylindrical shell-and-tube heat exchanger. The apparatus used was constructed from glass enabling direct observation of the flow using a dye release technique and had ten traversable pressure instrumented tubes permitting detailed pressure distributions to be obtained. The `exchanger' had a large tube bundle (278 tubes) and main flow areas typical of practical designs. Six geometries were studied: three baffle spacings both with and without baffle leakage. Results are also presented of three-dimensional modelling of shellside flows using the Harwell Laboratory's FLOW3D code. Flow visualisation provided flow patterns in the central plane of the bundle and adjacent to the shell wall. Comparison of these high-lighted significant radial flow variations. In particular, separated regions, originating from the baffle tips, were observed. The size of these regions was small in the bundle central plane but large adjacent to the shell wall and extended into the bypass lane. This appeared to reduce the bypass flow area and hence the bypass flow fraction. The three-dimensional flow modelling results were presented as velocity vector and isobar maps. The vector maps illustrated regions of high and low velocity which could be prone to tube vibration and fouling. Separated regions were also in evidence. A non-uniform crossflow was discovered with, in general, higher velocities in the central plane of the bundle than near the shell wall._The form of the isobar maps calculated by FLOW3D was in good agreement with experimental results. In particular, larger pressure drops occurred across the inlet than outlet of a crossflow region and were higher near the upstream than downstream baffle face. The effect of baffle spacing and baffle leakage on crossflow and window pressure drop measurements was identified. Agreement between the current measurements, previously obtained data and commonly used design correlations/models was, in general, poor. This was explained in terms of the increased understanding of shellside flow. The bulk of previous data, which dervies from small-scale rigs with few tubes, have been shown to be unrepresentative of typical commerical units. The Heat Transfer and Fluid Flow Service design program TASC provided the best predictions of the current pressure drop results. However, a number of simple one-dimensional models in TASC are, individually, questionable. Some revised models have been proposed.

Studies in three phase gas-liquid-solid fluidised systems

The work is a logical continuation of research started at Aston some years ago when studies were conducted on fermentations in bubble columns. The present work highlights typical design and operating problems that could arise in such systems as waste water, chemical, biochemical and petroleum operations involving three-phase, gas-liquid-solid fluidisation; such systems are in increasing use. It is believed that this is one of few studies concerned with `true' three-phase, gas-liquid-solid fluidised systems, and that this work will contribute significantly to closing some of the gaps in knowledge in this area. The research work was mainly experimentally based and involved studies of the hydrodynamic parameters, phase holdups (gas and solid), particle mixing and segregation, and phase flow dynamics (flow regime and circulation patterns). The studies have focused particularly on the solid behaviour and the influence of properties of solids present on the above parameters in three-phase, gas-liquid-solid fluidised systems containing single particle components and those containing binary and ternary mixtures of particles. All particles were near spherical in shape and two particle sizes and total concentration levels were used. Experiments were carried out in two- and three-dimensional bubble columns. Quantitative results are presented in graphical form and are supported by qualitative results from visual studies which are also shown as schematic diagrams and in photographic form. Gas and solid holdup results are compared for air-water containing single, binary and ternary component particle mixtures. It should be noted that the criteria for selection of the materials used are very important if true three-phase fluidisation is to be achieved: this is very evident when comparing the results with those in the literature. The fluid flow and circulation patterns observed were assessed for validation of the generally accepted patterns, and the author believes that the present work provides more accurate insight into the modelling of liquid circulation in bubble columns. The characteristic bubbly flow at low gas velocity in a two-phase system is suppressed in the three-phase system. The degree of mixing within the system is found to be dependent on flow regime, liquid circulation and the ratio of solid phase physical properties. Evidence of strong `trade-off' of properties is shown; the overall solid holdup is believed to be a major parameter influencing the gas holdup structure.

Effect of aspect ratio on thermal convection in open and closed systems

Internally heated fluids are found across the nuclear fuel cycle. In certain situations the motion of the fluid is driven by the decay heat (i.e. corium melt pools in severe accidents, the shutdown of liquid metal reactors, molten salt and the passive control of light water reactors) as well as normal operation (i.e. intermediate waste storage and generation IV reactor designs). This can in the long-term affect reactor vessel integrity or lead to localized hot spots and accumulation of solid wastes that may prompt local increases in activity. Two approaches to the modeling of internally heated convection are presented here. These are based on numerical analysis using codes developed in-house and simulations using widely available computational fluid dynamics solvers. Open and closed fluid layers at around the transition between conduction and convection of various aspect ratios are considered. We determine optimum domain aspect ratio (1:7:7 up to 1:24:24 for open systems and 5:5:1, 1:10:10 and 1:20:20 for closed systems), mesh resolutions and turbulence models required to accurately and efficiently capture the convection structures that evolve when perturbing the conductive state of the fluid layer. Note that the open and closed fluid layers we study here are bounded by a conducting surface over an insulating surface. Conclusions will be drawn on the influence of the periodic boundary conditions on the flow patterns observed. We have also examined the stability of the nonlinear solutions that we found with the aim of identifying the bifurcation sequence of these solutions en route to turbulence.