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.

Characterization of the hairpin vortex solution in plane Couette flow

Quantitative evidence that establishes the existence of the hairpin vortex state (HVS) in plane Couette flow (PCF) is provided in this work. The evidence presented in this paper shows that the HVS can be obtained via homotopy from a flow with a simple geometrical configuration, namely, the laterally heated flow (LHF). Although the early stages of bifurcations of LHF have been previously investigated, our linear stability analysis reveals that the root in the LHF yields multiple branches via symmetry breaking. These branches connect to the PCF manifold as steady nonlinear amplitude solutions. Moreover, we show that the HVS has a direct bifurcation route to the Rayleigh-Bénard convection. © 2010 The American Physical Society.

The flow behaviour of shallow fluidised beds

The flow behaviour of shallow gas-fluidised beds was studied. experimentally using a rotational viscometer, and an inclined open channel. Initially, tests were carried out with the viscometer in order to establish qualitative trends in the flow properties of a variety of materials over a wide range of fluidising conditions. Also, a technique was developed which enabled quantitative viscosity data to be extracted from the experimental results. The flow properties were found to be sensitive to the size, size-range and density of the fluidised material, the type of distributor used, and the moisture content of the fluidising gas. Tests in beds up to 120 mm deep showed that the fluidity of the bed improves with reduction in depth; and indicated a range of flow behaviour from shear-thinning to Newtonian, depending chiefly on fluidising velocity .. Later, an apparatus was built which provided for a steady, continuous flow of fluidised material down an inclined open channel of 3m length x 0.15m square, up to a mass flowrate of 10 kg/s (35 ton/hr). This facility has enabled data to be obtained that is of practical value in industrial applications; which is otherwise difficult in view of the present limited understanding of the true mechanism of fluidised flow. A correlation has been devised, based on analogy with laminar liquid flow, which describes the channel flow behaviour with reasonable accuracy over the whole range of shear-rates used. 1he channeI results indicated that at low fluidiising velocities the flow was adversely affected by settlement of a stagnant layer of particles on to the distributor, which gave rise to increased flow resistance. Conversely, at higher fluidising velocities the resistance at the distributor appeared to be less than at the walls. In view of this, and also because of the disparity in shear-rates between the two types of apparatus, it is not possible as yet to predict exactly the flow behaviour in an open channel from small-scale viscometer tests.

Transition in plane parallel shear flows heated internally

The stability of internally heated inclined plane parallel shear flows is examined numerically for the case of finite value of the Prandtl number, Pr. The transition in a vertical channel has already been studied for 0≤Pr≤100 with or without the application of an external pressure gradient, where the secondary flow takes the form of travelling waves (TWs) that are spanwise-independent (see works of Nagata and Generalis). In this work, in contrast to work already reported (J. Heat Trans. T. ASME 124 (2002) 635-642), we examine transition where the secondary flow takes the form of longitudinal rolls (LRs), which are independent of the steamwise direction, for Pr=7 and for a specific value of the angle of inclination of the fluid layer without the application of an external pressure gradient. We find possible bifurcation points of the secondary flow by performing a linear stability analysis that determines the neutral curve, where the basic flow, which can have two inflection points, loses stability. The linear stability of the secondary flow against three-dimensional perturbations is also examined numerically for the same value of the angle of inclination by employing Floquet theory. We identify possible bifurcation points for the tertiary flow and show that the bifurcation can be either monotone or oscillatory. © 2003 Académie des sciences. Published by Elsevier SAS. All rights reserved.

Concrete beams with stirrups subject to torsion, shear and bending

This thesis examines experimentally and theoretically the behaviour and ultimate strength of rectangular reinforced concrete members under combined torsion, shear and bending. The experimental investigation consists of the test results of 38 longitudinally and transversely reinforced concrete beams subjected to combined loads, ten beams of which were tested under pure torsion and self-weight. The behaviour of each test beam from application of the first increment of load until failure is presented. The effects of concrete strength, spacing of the stirrups, the amount of longitudinal steel and the breadth of the section on the ultimate torsional capacity are investigated. Based on the skew-bending mechanism, compatibility, and linear stress-strain relationship for the concrete and the steel, simple rational equations are derived for the three principal modes of failure for the following four types of failure observed in the tests: TYPE I Yielding the reinforcement, at failure, before crushing the concrete. TYPE II Yielding of the web steel only, at failure, before crushing the concrete. TYPE III Yielding of the longitudinal steel only, at failure, before crushing the concrete. TYPE IV Crushing of the concrete, at failure, before yielding of any of the reinforcement.

Prestressed concrete beams with stirrups subjected to torsion, bending and shear

Reported in this thesis are test results of 37 eccentrically prestressed beams with stirrups. Single variable parameters were investigated including the prestressing force, the prestressing steel area, the concrete strength, the aspect ratio h/b and the stirrups size and spacing. Interaction of bending, torsion and shear was also investigated by testing a series of beams subjected to varying bending/torsional moment ratios. For the torsional strength an empirical expression of linear format is proposed and can be rearranged in a non-dimensional interaction form: T/To+V/Vo+M/Mo+Ps/Po+Fs/Fo=Pc2/Fsp. This formula which is based on an average experimental steel stress lower than the yield point is compared with 243 prestressed beams containing ' stirrups, including the author's test beams, and good agreement is obtained. For the theoretical analysis of the problem of torsion combined with bending and shear in concrete beams with stirrups, the method of torque-friction is proposed and developed using an average steel stress. A general linear interaction equation for combined torsion with bending and/or shear is proposed in the following format: (fi) T/Tu=1 where (fi) is a combined loading factor to modify the pure ultimate strength for differing cases of torsion with bending and/or shear. From the analysis of 282 reinforced and prestressed concrete beams containing stirrups, including the present investigation, good agreement is obtained between the method and the test results. It is concluded that the proposed method provides a rational and simple basis for predicting the ultimate torisional strength and may also be developed for design purposes.

Bifurcation contrasts between plane Poiseuille flow and plane Magnetohydrodynamic flow

The stability characteristics of an incompressible viscous pressure-driven flow of an electrically conducting fluid between two parallel boundaries in the presence of a transverse magnetic field are compared and contrasted with those of Plane Poiseuille flow (PPF). Assuming that the outer regions adjacent to the fluid layer are perfectly electrically insulating, the appropriate boundary conditions are applied. The eigenvalue problems are then solved numerically to obtain the critical Reynolds number Rec and the critical wave number ac in the limit of small Hartmann number (M) range to produce the curves of marginal stability. The non-linear two-dimensional travelling waves that bifurcate by way of a Hopf bifurcation from the neutral curves are approximated by a truncated Fourier series in the streamwise direction. Two and three dimensional secondary disturbances are applied to both the constant pressure and constant flux equilibrium solutions using Floquet theory as this is believed to be the generic mechanism of instability in shear flows. The change in shape of the undisturbed velocity profile caused by the magnetic field is found to be the dominant factor. Consequently the critical Reynolds number is found to increase rapidly with increasing M so the transverse magnetic field has a powerful stabilising effect on this type of flow.

Cell exclusion in couette flow:evaluation through flow visualisation and mechanical forces

Cell exclusion is the phenomenon whereby the hematocrit and viscosity of blood decrease in areas of high stress. While this is well known in naturally occurring Poiseuille flow in the human body, it has never previously been shown in Couette flow, which occurs in implantable devices including blood pumps. The high-shear stresses that occur in the gap between the boundaries in Couette flow are known to cause hemolysis in erythrocytes. We propose to mitigate this damage by initiating cell exclusion through the use of a spiral-groove bearing (SGB) that will provide escape routes by which the cells may separate themselves from the plasma and the high stresses in the gap. The force between two bearings (one being the SGB) in Couette flow was measured. Stained erythrocytes, along with silver spheres of similar diameter to erythrocytes, were visualized across a transparent SGB at various gap heights. A reduction in the force across the bearing for human blood, compared with fluids of comparable viscosity, was found. This indicates a reduction in the viscosity of the fluid across the bearing due to a lowered hematocrit because of cell exclusion. The corresponding images clearly show both cells and spheres being excluded from the gap by entering the grooves. This is the first time the phenomenon of cell exclusion has been shown in Couette flow. It not only furthers our understanding of how blood responds to different flows but could also lead to improvements in the future design of medical devices.

Effect of oscillatory flow on nucleation kinetics of butyl paraben

More than 165 induction times of butyl paraben-ethanol solution in a batch moving fluid oscillation baffled crystallizer with various amplitudes (1-9 mm) and frequencies (1.0-9.0 Hz) have been determined to study the effect of COBR operating conditions on nucleation. The induction time decreases with increasing amplitude and frequency at power density below about 500 W/m3; however, a further increase of the frequency and amplitude leads to an increase of the induction time. The interfacial energies and pre-exponential factors in both homogeneous and heterogeneous nucleation are determined by classical nucleation theory at oscillatory frequency 2.0 Hz and amplitudes of 3 or 5 mm both with and without net flow. To capture the shear rate conditions in oscillatory flow crystallizers, a large eddy simulation approach in a computational fluid dynamics framework is applied. Under ideal conditions the shear rate distribution shows spatial and temporal periodicity and radial symmetry. The spatial distributions of the shear rate indicate an increase of average and maximum values of the shear rate with increasing amplitude and frequency. In continuous operation, net flow enhances the shear rate at most time points, promoting nucleation. The mechanism of the shear rate influence on nucleation is discussed.

Exploration of a simple, low cost, micro gas turbine recuperator solution for a domestic combined heat and power unit

For micro gas turbines (MGT) of around 1 kW or less, a commercially suitable recuperator must be used to produce a thermal efficiency suitable for use in UK Domestic Combined Heat and Power (DCHP). This paper uses computational fluid dynamics (CFD) to investigate a recuperator design based on a helically coiled pipe-in-pipe heat exchanger which utilises industry standard stock materials and manufacturing techniques. A suitable mesh strategy was established by geometrically modelling separate boundary layer volumes to satisfy y + near wall conditions. A higher mesh density was then used to resolve the core flow. A coiled pipe-in-pipe recuperator solution for a 1 kW MGT DCHP unit was established within the volume envelope suitable for a domestic wall-hung boiler. Using a low MGT pressure ratio (necessitated by using a turbocharger oil cooled journal bearing platform) meant unit size was larger than anticipated. Raising MGT pressure ratio from 2.15 to 2.5 could significantly reduce recuperator volume. Dimensional reasoning confirmed the existence of optimum pipe diameter combinations for minimum pressure drop. Maximum heat exchanger effectiveness was achieved using an optimum or minimum pressure drop pipe combination with large pipe length as opposed to a large pressure drop pipe combination with shorter pipe length. © 2011 Elsevier Ltd. All rights reserved.

A compartmental CFD-PBM model of high shear wet granulation

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

Practical secure information flow in programming languages

If we classify variables in a program into various security levels, then a secure information flow analysis aims to verify statically that information in a program can flow only in ways consistent with the specified security levels. One well-studied approach is to formulate the rules of the secure information flow analysis as a type system. A major trend of recent research focuses on how to accommodate various sophisticated modern language features. However, this approach often leads to overly complicated and restrictive type systems, making them unfit for practical use. Also, problems essential to practical use, such as type inference and error reporting, have received little attention. This dissertation identified and solved major theoretical and practical hurdles to the application of secure information flow. ^ We adopted a minimalist approach to designing our language to ensure a simple lenient type system. We started out with a small simple imperative language and only added features that we deemed most important for practical use. One language feature we addressed is arrays. Due to the various leaking channels associated with array operations, arrays have received complicated and restrictive typing rules in other secure languages. We presented a novel approach for lenient array operations, which lead to simple and lenient typing of arrays. ^ Type inference is necessary because usually a user is only concerned with the security types for input/output variables of a program and would like to have all types for auxiliary variables inferred automatically. We presented a type inference algorithm B and proved its soundness and completeness. Moreover, algorithm B stays close to the program and the type system and therefore facilitates informative error reporting that is generated in a cascading fashion. Algorithm B and error reporting have been implemented and tested. ^ Lastly, we presented a novel framework for developing applications that ensure user information privacy. In this framework, core computations are defined as code modules that involve input/output data from multiple parties. Incrementally, secure flow policies are refined based on feedback from the type checking/inference. Core computations only interact with code modules from involved parties through well-defined interfaces. All code modules are digitally signed to ensure their authenticity and integrity. ^

Eulerian-Lagrangian two phase debris flow model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier-Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Secure information flow via stripping and fast simulation

Type systems for secure information flow aim to prevent a program from leaking information from H (high) to L (low) variables. Traditionally, bisimulation has been the prevalent technique for proving the soundness of such systems. This work introduces a new proof technique based on stripping and fast simulation, and shows that it can be applied in a number of cases where bisimulation fails. We present a progressive development of this technique over a representative sample of languages including a simple imperative language (core theory), a multiprocessing nondeterministic language, a probabilistic language, and a language with cryptographic primitives. In the core theory we illustrate the key concepts of this technique in a basic setting. A fast low simulation in the context of transition systems is a binary relation where simulating states can match the moves of simulated states while maintaining the equivalence of low variables; stripping is a function that removes high commands from programs. We show that we can prove secure information flow by arguing that the stripping relation is a fast low simulation. We then extend the core theory to an abstract distributed language under a nondeterministic scheduler. Next, we extend to a probabilistic language with a random assignment command; we generalize fast simulation to the setting of discrete time Markov Chains, and prove approximate probabilistic noninterference. Finally, we introduce cryptographic primitives into the probabilistic language and prove computational noninterference, provided that the underling encryption scheme is secure.