Mechanisms of drying of particulate slurries

The literature on the evaporation of pure liquid drops and the drying of drops of solutions and slurries has been reviewed with particular reference to spray drying. A 0.1-0.2 mm glass filament-thermocouple was constructed and used to study simultaneously, heat and mass transfer from a single suspended drop placed in a humidity and temperature controlled, 28 mm OD vertical wind tunnel. Heat conduction through the filament was minimised eg at 100¦C it accounted for only 9.3% of the total heat transferred to a drop. Evaporation of single water drops was also studied in a 101 mm OD vertical wind tunnel. The Nusselt number was found to be a function of the Reynolds, Prandtl and Transfer number over an air temperature range of 17¦C to 107¦C. The proposed correlation is: Nu = 2+(-12.96B+0.76)Re¦-5Pr0-33 Experimental drying studies were carried out on single suspended 1 to 2.5 mm diameter drops of aqueous sodium sulphate decahydrate, sodium chloride, potassium sulphate, copper sulphate and sodium acetate solutions and slurries at temperatures of 20¦C to 124¦C. Dried crusts were examined by Scanning Electron Microscopy. The drying history of any material depended upon the nature of the crust formed. Sodium acetate formed a non-rigid skin prior to the formation of a rigid crust. A modified receding evaporation interface model was proposed for the drying of solutions and slurries. This covered both the constant rate period prior to crust formation and the subsequent falling rate period. The model was solved numerically for the variation in core temperature, drop weight and crust thickness. Good agreement was obtained between model predictions and experimental results for materials forming rigid crusts i.e. sodium sulphate decahydrate, sodium chloride, potassium sulphate and copper sulphate. However, the drying histories of drops of 10-20% weight initial concentration sodium acetate were unpredictable since formation of a non-rigid skin deviated from the model assumption of a rigid outer surface. At higher initial concentrations (40% weight) where a rigid crust was formed for sodium acetate, good agreement was obtained between experimental results and model predictions. Single suspended drop studies are concluded to provide a valuable insight into the drying mechanisms of specific solutions and slurries.

A solar powered liquid-desiccant cooling system for greenhouses

Liquid desiccant systems are of potential interest as a means of cooling greenhouses to temperatures below those achieved by conventional means. However, only very little work has been done on this technology with previous workers focussing on the cooling of human dwellings using expensive desiccants such as lithium salts. In this study we are designing a system for greenhouse cooling based on magnesium chloride desiccant which is an abundant and non-toxic substance. Magnesium chloride is found in seawater, for example, and is a by-product from solar salt works. We have carried out a detailed experimental study of the relevant properties of magnesium rich solutions. In addition we have constructed a test rig that includes the main components of the cooling system, namely a dehumidifier and solar regenerator. The dehumidifier is a cross-flow device that consists of a structured packing made of corrugated cellulose paper sheets with different flute angles and embedded cooling tubes. The regenerator is of the open type with insulated backing and fabric covering to spread the flow of desiccant solution. Alongside these experiments we are developing a mathematical model in gPROMS® that combines and simulates the heat and mass transfer processes in these components. The model can be applied to various geographical locations. Here we report predictions for Havana (Cuba) and Manila (Philippines), where we find that average wet-bulb temperatures can be lowered by 2.2 and 3°C, respectively, during the month of May.

Bed expansion:extending correlations for fluidized beds with inserts to open fluidized beds

Bed expansion occurs during the operation of gas-fluidized beds and is influenced by particle properties, gas properties and distributor characteristics. It has a significant bearing on heat and mass transfer phenomena within the bed. A method of predicting bed expansion behavior from other fluidizing parameters would be a useful tool in the design process, dispensing with the need for small-scale trials. This study builds on previous work on fluidized beds with vertical inserts to produce a correlation that links a modified particle terminal velocity, minimum fluidizing velocity and distributor characteristics with bed voidage in the relationship with P as the pitch between holes in the perforated distributor plate. © 2007 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Evaporation maps for ternary non-ideal liquid mixtures

This thesis deals with the evaporation of non-ideal liquid mixtures using a multicomponent mass transfer approach. It develops the concept of evaporation maps as a convenient way of representing the dynamic composition changes of ternary mixtures during an evaporation process. Evaporation maps represent the residual composition of evaporating ternary non-ideal mixtures over the full range of composition, and are analogous to the commonly-used residue curve maps of simple distillation processes. The evaporation process initially considered in this work involves gas-phase limited evaporation from a liquid or wetted-solid surface, over which a gas flows at known conditions. Evaporation may occur into a pure inert gas, or into one pre-loaded with a known fraction of one of the ternary components. To explore multicomponent masstransfer effects, a model is developed that uses an exact solution to the Maxwell-Stefan equations for mass transfer in the gas film, with a lumped approach applied to the liquid phase. Solutions to the evaporation model take the form of trajectories in temperaturecomposition space, which are then projected onto a ternary diagram to form the map. Novel algorithms are developed for computation of pseudo-azeotropes in the evaporating mixture, and for calculation of the multicomponent wet-bulb temperature at a given liquid composition. A numerical continuation method is used to track the bifurcations which occur in the evaporation maps, where the composition of one component of the pre-loaded gas is the bifurcation parameter. The bifurcation diagrams can in principle be used to determine the required gas composition to produce a specific terminal composition in the liquid. A simple homotopy method is developed to track the locations of the various possible pseudo-azeotropes in the mixture. The stability of pseudo-azeotropes in the gas-phase limited case is examined using a linearized analysis of the governing equations. Algorithms for the calculation of separation boundaries in the evaporation maps are developed using an optimization-based method, as well as a method employing eigenvectors derived from the linearized analysis. The flexure of the wet-bulb temperature surface is explored, and it is shown how evaporation trajectories cross ridges and valleys, so that ridges and valleys of the surface do not coincide with separation boundaries. Finally, the assumption of gas-phase limited mass transfer is relaxed, by employing a model that includes diffusion in the liquid phase. A finite-volume method is used to solve the system of partial differential equations that results. The evaporation trajectories for the distributed model reduce to those of the lumped (gas-phase limited) model as the diffusivity in the liquid increases; under the same gas-phase conditions the permissible terminal compositions of the distributed and lumped models are the same.

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.

Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

A study of thermodynamic characteristics and predictive computer modelling of an air to water heat pump system

A study on heat pump thermodynamic characteristics has been made in the laboratory on a specially designed and instrumented air to water heat pump system. The design, using refrigerant R12, was based on the requirement to produce domestic hot water at a temperature of about 50 °C and was assembled in the laboratory. All the experimental data were fed to a microcomputer and stored on disk automatically from appropriate transducers via amplifier and 16 channel analogue to digital converters. The measurements taken were R12 pressures and temperatures, water and R12 mass flow rates, air speed, fan and compressor input powers, water and air inlet and outlet temperatures, wet and dry bulb temperatures. The time interval between the observations could be varied. The results showed, as expected, that the COP was higher at higher air inlet temperatures and at lower hot water output temperatures. The optimum air speed was found to be at a speed when the fan input power was about 4% of the condenser heat output. It was also found that the hot water can be produced at a temperature higher than the appropriate R12 condensing temperature corresponding to condensing pressure. This was achieved by condenser design to take advantage of discharge superheat and by further heating the water using heat recovery from the compressor. Of the input power to the compressor, typically about 85% was transferred to the refrigerant, 50 % by the compression work and 35% due to the heating of the refrigerant by the cylinder wall, and the remaining 15% (of the input power) was rejected to the cooling medium. The evaporator effectiveness was found to be about 75% and sensitive to the air speed. Using the data collected, a steady state computer model was developed. For given input conditions s air inlet temperature, air speed, the degree of suction superheat , water inlet and outlet temperatures; the model is capable of predicting the refrigerant cycle, compressor efficiency, evaporator effectiveness, condenser water flow rate and system Cop.

Effects of temperature-dependent viscosity variation on entropy generation, heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford, is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Péclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case.

Modelling the exothermic heat and the demand of cooling in paper and pulp industry biological wastewater treatment

This master thesis presents a study on the requisite cooling of an activated sludge process in paper and pulp industry. The energy consumption of paper and pulp industry and it’s wastewater treatment plant in particular is relatively high. It is therefore useful to understand the wastewater treatment process of such industries. The activated sludge process is a biological mechanism which degrades carbonaceous compounds that are present in waste. The modified activated sludge model constructed here aims to imitate the bio-kinetics of an activated sludge process. However, due to the complicated non-linear behavior of the biological process, modelling this system is laborious and intriguing. We attempt to find a system solution first using steady-state modelling of Activated Sludge Model number 1 (ASM1), approached by Euler’s method and an ordinary differential equation solver. Furthermore, an enthalpy study of paper and pulp industry’s vital pollutants was carried out and applied to revise the temperature shift over a period of time to formulate the operation of cooling water. This finding will lead to a forecast of the plant process execution in a cost-effective manner and management of effluent efficiency. The final stage of the thesis was achieved by optimizing the steady state of ASM1.

Comparison of heat-treated rapeseed expeller and solvent-extracted soya-bean meal as protein supplements for dairy cows given grass silage-based diets

Sixteen early to mid lactation Finnish Ayrshire dairy cows were used in a cyclic change-over experiment with four 21-day experimental periods and a 4 5 2 factorial arrangement of treatments to evaluate the effects of heat-treated rapeseed expeller and solvent-extracted soya-bean meal protein supplements on animal performance. Dietary treatments consisted of grass silage offered ad libitum supplemented with a fixed amount of a cereal based concentrate (10 kg/day on a fresh weight basis) containing 120, 150, 180 or 210 g crude protein (CP) per kg dry matter (DM). Concentrate CP content was manipulated by replacement of basal ingredients (g/kg) with either rapeseed expeller (R; 120, 240 and 360) or soya-bean meal (S; 80, 160 and 240). Increases in concentrate CP stimulated linear increases (P < 0.05) in silage intake (mean 22.5 and 23.8 g DM per g/kg increase in dietary CP content, for R and S, respectively) and milk production. Concentrate inclusion of rapeseed expeller elicited higher (P < 0.01) milk yield and milk protein output responses (mean 108 and 3.71 g/day per g/kg DM increase in dietary CP content) than soya-bean meal (corresponding values 62 and 2.57). Improvements in the apparent utilization of dietary nitrogen for milk protein synthesis (mean 0.282 and 0.274, for R and S, respectively) were associated with higher (P < 0.05) plasma concentrations of histidine, branched-chain, essential and total amino acids (35, 482, 902 and 2240 and 26, 410, 800 and 2119 mu mol/l, respectively) and lower (P < 0.01) concentrations of urea (corresponding values 4.11 and 4.52 mmol/l). Heat-treated rapeseed expeller proved to be a more effective protein supplement than solvent-extracted soya-bean meal for cows offered grass silage-based diets.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

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.