Feasibility of a solar panel-powered liquid desiccant cooling system for greenhouses

To investigate the technical feasibility of a novel cooling system for commercial greenhouses, knowledge of the state of the art in greenhouse cooling is required. An extensive literature review was carried out that highlighted the physical processes of greenhouse cooling and showed the limitations of the conventional technology. The proposed cooling system utilises liquid desiccant technology; hence knowledge of liquid desiccant cooling is also a prerequisite before designing such a system. Extensive literature reviews on solar liquid desiccant regenerators and desiccators, which are essential parts of liquid desiccant cooling systems, were carried out to identify their advantages and disadvantages. In response to the findings, a regenerator and a desiccator were designed and constructed in lab. An important factor of liquid desiccant cooling is the choice of liquid desiccant itself. The hygroscopicity of the liquid desiccant affects the performance of the system. Bitterns, which are magnesium-rich brines derived from seawater, are proposed as an alternative liquid desiccant for cooling greenhouses. A thorough experimental and theoretical study was carried out in order to determine the properties of concentrated bitterns. It was concluded that their properties resemble pure magnesium chloride solutions. Therefore, magnesium chloride solution was used in laboratory experiments to assess the performance of the regenerator and the desiccator. To predict the whole system performance, the physical processes of heat and mass transfer were modelled using gPROMS® advanced process modelling software. The model was validated against the experimental results. Consequently it was used to model a commercials-scale greenhouse in several hot coastal areas in the tropics and sub-tropics. These case studies show that the system, when compared to evaporative cooling, achieves 3oC-5.6oC temperature drop inside the greenhouse in hot and humid places (RH>70%) and 2oC-4oC temperature drop in hot and dry places (50%

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.

Mechanisms of convective drying in thick beds of solids

The literature on heat and mass transfer mechanisms in the convective drying of thick beds of solids has been critically reviewed. Related mathematical models of heat transfer are also considered. Experimental and theoretical studies were made of the temperature distribution within beds, and of drying rates, with various materials undergoing convective drying. The experimental work covered thick beds of hygroscopic and non-hygroscopic materials (glass beads of different diameters, polystyrene pellets, activated alumina and wood powder) at air temperatures of 54°C to 84°C. Tests were carried out in a laboratory drying apparatus comprising a wind tunnel through which the air, of controlled temperature and humidity, was passed over a sample suspended from a balance. Thermocouples were inserted at different depths within the sample bed. The temperature distribution profiles for both hygroscopic and non-hygroscopic beds exhibited a clear difference between the temperatures at the surface and bottom during the constant rate period. An effective method was introduced for predicting the critical moisture content. During the falling rate the profiles showed the existence of a receding evaporation plane; this divided the system into a hotter dry zone in the upper section and a wet zone near the bottom. A graphical procedure was established to predict accurately the position of the receding evaporation front at any time. A new mathematical model, based on the receding evaporation front phenomenon, was proposed to predict temperature distributions throughout a bed during drying. Good agreement was obtained when the model was validated by comparing its predictions with experimental data. The model was also able to predict the duration of each drying stage. In experiments using sample trays of different diameters, the drying rate was found to increase with a decrease in the effective length of the bed surface. During the constant rate period with trays of a small effective length, i.e. less than 0.08 m, an 'inversion' in temperature distribution occurred in the bed; the bottom temperature increased and became greater than that of the surface. Experimental measurements were verified in several ways to ensure this phenomenon was real. Theoretical explanations are given for both the effective length and temperature inversion phenomena.

Technical and economic study of plantain food products

Plantain (Banana-Musa AAB) is a widely growing but commercially underexploited tropical fruit. This study demonstrates the processing of plantain to flour and extends its use and convenience as a constituent of bread, cake and biscuit. Plantain was peeled, dried and milled to produce flour. Proximate analysis was carried out on the flour to determine the food composition. Drying at temperatures below 70ºC produced light coloured plantain flour. Experiments were carried out to determine the mechanism of drying, the heat and mass transfer coefficients, effect of air velocity, temperature and cube size on the rate of drying of plantain cubes. The drying was diffusion controlled. Pilot scale drying of plantain cubes in a cabinet dryer showed no significant increase of drying rate above 70ºC. In the temperature range found most suitable for plantain drying (ie 60 to 70ºC) the total drying time was adequately predicted using a modified equation based on Fick's Law provided the cube temperature was taken to be about 5ºC below the actual drying air temperature. Studies of baking properties of plantain flour revealed that plantain flour can be substituted for strong wheat flour up to 15% for bread making and up to 50% for madeira cake. A shortcake biscuit was produced using 100% plantain flour and test-marketed. Detailed economic studies showed that the production of plantain fruit and its processing into flour would be economically viable in Nigeria when the flour is sold at the wholesale price of NO.65 per kilogram provided a minimum sale of 25% plantain suckers. There is need for government subsidy if plantain flour is to compete with imported wheat flour. The broader economic benefits accruing from the processing of plantain fruit into flour and its use in bakery products include employment opportunity, savings in foreign exchange and stimulus to home agriculture.

Solar stills:a comprehensive review of designs, performance and material advances

The demand for fresh water production is growing day by day with the increase in world population and with industrial growth. Use of desalination technology is increasing to meet this demand. Among desalination technologies, solar stills require low maintenance and are readily affordable; however their productivity is limited. This paper aims to give a detailed review about the various types of solar stills, covering passive and active designs, single- and multi-effect types, and the various modifications for improved productivity including reflectors, heat storage, fins, collectors, condensers, and mechanisms for enhancing heat and mass transfer. Photovoltaic-thermal and greenhouse type solar stills are also covered. Material advances in the area of phase change materials and nanocomposites are very promising to enhance further performance; future research should be carried out in these and other areas for the greater uptake of solar still technology.

Asymmetric hydrogen flame in a heated micro-channel: role of Darrieus-Landau and thermal-diffusive instabilities

Present work examines numerically the asymmetric behavior of hydrogen/air flame in a micro-channel subjected to a non-uniform wall temperature distribution. A high resolution (with cell size of 25 μm × 25 μm) of two-dimensional transient Navier–Stokes simulation is conducted in the low-Mach number formulation using detailed chemistry evolving 9 chemical species and 21 elementary reactions. Firstly, effects of hydrodynamic and diffusive-thermal instabilities are studied by performing the computations for different Lewis numbers. Then, the effects of preferential diffusion of heat and mass transfer on the asymmetric behavior of the hydrogen flame are analyzed for different inlet velocities and equivalence ratios. Results show that for the flames in micro-channels, interactions between thermal diffusion and molecular diffusion play major role in evolution of a symmetric flame into an asymmetric one. Furthermore, the role of Darrieus–Landau instability found to be minor. It is also found that in symmetric flames, the Lewis number decreases behind the flame front. This is related to the curvature of flame which leads to the inclination of thermal and mass fluxes. The mass diffusion vectors point toward the walls and the thermal diffusion vectors point toward the centerline. Asymmetric flame is observed when the length of flame front is about 1.1–1.15 times of the channel width.

Effect of drying on the physical, chemical and sensorial properties of kiwi

Kiwi fruit is a highly nutritional fruit due to the high level of vitamin C and its strong antioxidant capacity due to a wide number of phytonutrients including carotenoids, lutein, phenolics, flavonoids and chlorophyll [1]. Drying consists of a complex process in which simultaneous heat and mass transfer occur. Several alterations occur during the drying of foods at many levels (physical, chemical, nutritional or sensorial) which are influenced by a number of factors, including processing conditions [2]. Temperature is particularly important because of the effects it produces at the chemical and also at the physical level, particularly colour and texture [3]. In the present work were evaluated the changes in sliced kiwi when exposed to air drying at different temperatures (50, 60, 70, 80 ºC), namely in terms of some chemical properties like ascorbic acid or phenolic compounds, physical characteristics like colour and texture and also at the sensorial level. All experiments followed standard established procedures and several replicates were done to assess each property. The results obtained indicated that moisture was reduced with drying by 74 to 87%, depending on the temperature. Also ascorbic acid decreased with drying, being 7% for 50 ºC and increasing up to 28% for the highest temperature (80 ºC). The phenolic compounds and antioxidant activity were also very much affected by the drying temperature. The water activity of the dried samples varied from 0.658 to 0.753, being compatible with a good preservation. Regarding colour, the total colour difference between the dried samples and the fresh sample was found to vary in the range 9.45 – 17.17. The textural parameters were also much affected by drying, namely hardness which decreased by 45 to 72 %, and all other parameters increased: cohesiveness (approximately doubled), springiness (increased 2 to 3 times) and chewiness which increased up to 2.5 times that off the fresh sample. Adhesiveness, which was observed for the fresh samples (-4.02 N.s) disappeared in all the dried samples. The sensorial analysis made to the dried samples allowed establishing the sensorial profiles as shown in Figure 1.

A wavelet-based adaptive technique for adsorption problems involving steep gradients

A general, fast wavelet-based adaptive collocation method is formulated for heat and mass transfer problems involving a steep moving profile of the dependent variable. The technique of grid adaptation is based on sparse point representation (SPR). The method is applied and tested for the case of a gas–solid non-catalytic reaction in a porous solid at high Thiele modulus. Accurate and convergent steep profiles are obtained for Thiele modulus as large as 100 for the case of slab and found to match the analytical solution.

Experimental Analysis of Froude Number Effect on Air Entrainment in the Hydraulic Jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in rectangular horizontal flume with partially-developed inflow conditions. The vertical distributions of void fraction and air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was provided.

Horizontal rotating cylinder - A compact apparatus for studying the effect of water wetting on carbon dioxide corrosion of mild steel

Water wetting is a crucial issue in carbon dioxide (CO.) corrosion of multiphase flow pipelines made from mild steel. This study demonstrates the use of a novel benchtop apparatus, a horizontal rotating cylinder, to study the effect of water wetting on CO2 corrosion of mild steel in two-phase flow. The setup is similar to a standard rotating cylinder except for its horizontal orientation and the presence of two phases-typically water and oil. The apparatus has been tested by using mass-transfer measurements and CO2 corrosion measurements in single-phase water flow. CO2 corrosion measurements were subsequently performed using a water/hexane mixture with water cuts varying between 5% and 50%. While the metal surface was primarily hydrophilic under stagnant. conditions, a variety of dynamic water wetting situations was encountered as the water cut and fluid velocity were altered. Threshold velocities were identified at various water cuts when the surface became oil-wet and corrosion stopped.

A mechanistic model for carbon dioxide corrosion of mild steel in the presence of protective iron carbonate films - Part 3: Film growth model

A model of iron carbonate (FeCO3) film growth is proposed, which is an extension of the recent mechanistic model of carbon dioxide (CO2) corrosion by Nesic, et al. In the present model, the film growth occurs by precipitation of iron carbonate once saturation is exceeded. The kinetics of precipitation is dependent on temperature and local species concentrations that are calculated by solving the coupled species transport equations. Precipitation tends to build up a layer of FeCO3 on the surface of the steel and reduce the corrosion rate. On the other hand, the corrosion process induces voids under the precipitated film, thus increasing the porosity and leading to a higher corrosion rate. Depending on the environmental parameters such as temperature, pH, CO2 partial pressure, velocity, etc., the balance of the two processes can lead to a variety of outcomes. Very protective films and low corrosion rates are predicted at high pH, temperature, CO2 partial pressure, and Fe2+ ion concentration due to formation of dense protective films as expected. The model has been successfully calibrated against limited experimental data. Parametric testing of the model has been done to gain insight into the effect of various environmental parameters on iron carbonate film formation. The trends shown in the predictions agreed well with the general understanding of the CO2 corrosion process in the presence of iron carbonate films. The present model confirms that the concept of scaling tendency is a good tool for predicting the likelihood of protective iron carbonate film formation.

Estudo comparativo entre o sistema convencional de arrefecimento e o arrefecedor evaporativo para os edifícios industriais portugueses

Em Portugal existem muitos espaços comerciais e industriais em que as necessidades térmicas de arrefecimento são muito superiores às necessidades de aquecimento devido aos ganhos internos que advêm da existência de equipamentos e da iluminação dos edifícios, assim como, da presença das pessoas. A instalação de sistemas convencionais de ar condicionado para espaços comerciais e industriais de grande dimensão está geralmente associada ao transporte de grandes caudais de ar, e consequentemente, a elevados consumos de energia primária, e também, elevados custos de investimento, de manutenção e de operação. O arrefecedor evaporativo é uma solução de climatização com elevada eficiência energética, cujo princípio de funcionamento promove a redução do consumo de energia primária nos edifícios. A metodologia utilizada baseou-se na criação de uma ferramenta informática de simulação do funcionamento de um protótipo de um arrefecedor evaporativo. Foi efetuada a modelação matemática das variáveis dinâmicas envolvidas, dos processos de transferência de calor e de massa, assim como dos balanços de energia que ocorrem no arrefecedor evaporativo. A ferramenta informática desenvolvida permite o dimensionamento do protótipo do arrefecedor evaporativo, sendo determinadas as caraterísticas técnicas (potência térmica, caudal, eficiência energética, consumo energético e consumo e água) de acordo com o tipo de edifício e com as condições climatéricas do ar exterior. Foram selecionados três dimensionamentos de arrefecedores evaporativos, representativos de condições reais de uma gama baixa, média e elevada de caudais de ar. Os resultados obtidos nas simulações mostram que a potência de arrefecimento (5,6 kW, 16,0 kW e 32,8 kW) e o consumo de água (8 l/h, 23,9 l/h e 48,96 l/h) aumentam com o caudal de ar do arrefecedor, 5.000 m3/h, 15.000 m3/h e 30.000 m3/h, respetivamente. A eficácia de permuta destes arrefecedores evaporativos, foi de 69%, 66% e 67%, respetivamente. Verificou-se que a alteração de zona climática de V1 para V2 implicou um aumento de 39% na potência de arrefecimento e de 20% no consumo de água, e que, a alteração de zona climática de V2 para V3 implicou um aumento de 39% na potência de arrefecimento e de 39% no consumo de água. O arrefecedor evaporativo apresenta valores de consumo de energia elétrica entre 40% a 80% inferiores aos dos sistemas de arrefecimento convencionais, sendo este efeito mais intenso quando a zona climática de verão se torna mais severa.

Modelação, simulação e implementação de agitação acústica para aplicações em microfluídica

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

Adaptive numerical simulation of reaction - diffusion systems

Fluidized beds, granulation, heat and mass transfer, calcium dynamics, stochastic process, finite element methods, Rosenbrock methods, multigrid methods, parallelization

Particle transport method for convection-diffusion-reaction problems

Convective transport, both pure and combined with diffusion and reaction, can be observed in a wide range of physical and industrial applications, such as heat and mass transfer, crystal growth or biomechanics. The numerical approximation of this class of problemscan present substantial difficulties clue to regions of high gradients (steep fronts) of the solution, where generation of spurious oscillations or smearing should be precluded. This work is devoted to the development of an efficient numerical technique to deal with pure linear convection and convection-dominated problems in the frame-work of convection-diffusion-reaction systems. The particle transport method, developed in this study, is based on using rneshless numerical particles which carry out the solution along the characteristics defining the convective transport. The resolution of steep fronts of the solution is controlled by a special spacial adaptivity procedure. The serni-Lagrangian particle transport method uses an Eulerian fixed grid to represent the solution. In the case of convection-diffusion-reaction problems, the method is combined with diffusion and reaction solvers within an operator splitting approach. To transfer the solution from the particle set onto the grid, a fast monotone projection technique is designed. Our numerical results confirm that the method has a spacial accuracy of the second order and can be faster than typical grid-based methods of the same order; for pure linear convection problems the method demonstrates optimal linear complexity. The method works on structured and unstructured meshes, demonstrating a high-resolution property in the regions of steep fronts of the solution. Moreover, the particle transport method can be successfully used for the numerical simulation of the real-life problems in, for example, chemical engineering.