Modelling of combustion and sorbent reactions in three-dimensional flow environment of a circulating fluidized bed furnace

This thesis presents a three-dimensional, semi-empirical, steady state model for simulating the combustion, gasification, and formation of emissions in circulating fluidized bed (CFB) processes. In a large-scale CFB furnace, the local feeding of fuel, air, and other input materials, as well as the limited mixing rate of different reactants produce inhomogeneous process conditions. To simulate the real conditions, the furnace should be modelled three-dimensionally or the three-dimensional effects should be taken into account. The only available methods for simulating the large CFB furnaces three-dimensionally are semi-empirical models, which apply a relatively coarse calculation mesh and a combination of fundamental conservation equations, theoretical models and empirical correlations. The number of such models is extremely small. The main objective of this work was to achieve a model which can be applied to calculating industrial scale CFB boilers and which can simulate all the essential sub-phenomena: fluid dynamics, reactions, the attrition of particles, and heat transfer. The core of the work was to develop the model frame and the required sub-models for determining the combustion and sorbent reactions. The objective was reached, and the developed model was successfully used for studying various industrial scale CFB boilers combusting different types of fuel. The model for sorbent reactions, which includes the main reactions for calcitic limestones, was applied for studying the new possible phenomena occurring in the oxygen-fired combustion. The presented combustion and sorbent models and principles can be utilized in other model approaches as well, including other empirical and semi-empirical model approaches, and CFD based simulations. The main achievement is the overall model frame which can be utilized for the further development and testing of new sub-models and theories, and for concentrating the knowledge gathered from the experimental work carried out at bench scale, pilot scale and industrial scale apparatus, and from the computational work performed by other modelling methods.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

REPRODUCTIVE SYSTEM AND POLLEN FLOW IN PROGENIES OF Qualea grandiflora Mart., A TYPICAL SPECIES OF THE BRAZILIAN CERRADO

This study analyzed the reproductive system and the pollen dispersion pattern of Qualea grandiflora progenies. This is a typical species from the Brazilian Cerrado about which there are not too many studies from the genetics point of view. The study was conducted in an area of 2.2 hectares located in the Conservation Unit managed by the Forest Institute of the state of São Paulo, Brazil (Assis State Forest). Total genomic DNA of 300 seeds from 25 plants (12 seeds from each plant) was extracted and amplified using specific primers to obtain microsatellite markers. Results showed that selfing is frequent among adults and progenies, and the species reproduces by outcrossing between related and unrelated individuals (0.913). The single-locus outcrossing rate was 0.632, which indicates that mating between unrelated individuals is more frequent than between related plants. The selfing rate was low (0.087), that is, the species is allogamous and self-fertilization is reduced. About 35% of the plants in the progenies were full-sibs, and about 57%, half-sibs. Besides, about 8% of the progenies were selfing siblings. The genetic differentiation coefficient within progenies was 0.139, whereas the fixation rate was about 27%. The estimate of the effective size revealed that the genetic representativeness of descent was lower than expected in random mating progenies: The analyzed samples corresponded to only 13.2 individuals of an ideal panmictic population. In environmental recovery programs, seeds, preferably from different fruits, should be collected from 95 trees to preserve the genetic diversity of the species.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Calibration of angles of nozzles and deposition of an axial-flow sprayer on dwarf cashew

The study aimed to determine an optimum angle for the nozzles axial-flow sprayers a deposition for better vertical distribution focused on cashew. In laboratory tests were conducted adjusting the angle of the nozzle axial-flow sprayers. The experimental design was randomized blocks in a 2x3 factorial with four replications. The treatment for this test were two settings (with and without the adjustment of the angles of the nozzles ) and tree application volumes 273, 699 and 954 L ha-¹.The study was conducted in an orchard of dwarf cashew, with eight years of age. It was concluded that the volumetric distribution profile showed better vertical distribution uniformity when the angles of the nozzles were regulated for the canopy, the adjustment of the angles of the nozzles for the canopy provided greater deposition of droplets, the increased volume of application resulted in higher depositions in the leaves of the crop.

Powder flow criteria for design of vertical silo walls

For design of vertical silos walls involving the storage of bulk solids to be safe and reliable, it is important knowing the largest possible number of variables such as: flow properties, silo geometry and pattern of flow desired. In order to validate the theories of flow prediction and design of conical hoppers, the flow properties of two bulk solids were determined, the theories of Jenike's flowability and Enstad and Walker for hopper design were analyzed and the results were compared with those experimentally obtained in a reduced model of a semicircular-section silo. Results show that Enstad theory for the hopper design is adequate to occur mass flow inside the silo, and for the sizing of the discharge outlet, the Walker's theory was closer to the appropriate than Jenike's theory, which was higher around 100% than the experimental hopper outlet.

Air flow humidification for air-assisted boom sprayer

The air-assisted ground spray is fairly widespread. However, due to the unpredictable weather conditions, the operational efficiency is impaired by stops on grounds of low humidity and high temperatures. The aim of this work was to assess an air humidification method and evaluate its impact on temperature and air humidity for the air curtain of the air-assisted sprayer. With respect to relative air humidity, it has increased in 6.59%, being the maximum change when inserting 1.92 L min-1. So, it is concluded that the pipeline humidification might significantly reduce temperature and enhance air humidity. The treatments performed in this study consisted of a varied flow of a humidity device, related to weather conditions. Temperature and relative air humidity were measured at 1.0 m height from right to left of middle point of the machine, corresponding to the end of the spray boom, in the middle and end of right spray boom. The readings were also performed at three different distances from the end of the pipeline and at 0.25 and 0.50 m from that to the soil. The results show that 0.48 L min-1 in the humidification system has promoted a better efficiency in reducing air-temperature, on average 2.52 ºC when compared to the non-humidified one.

Low-flow estimates in regions of extrapolation of the regionalization equations: a new concept

ABSTRACT Knowledge of natural water availability, which is characterized by low flows, is essential for planning and management of water resources. One of the most widely used hydrological techniques to determine streamflow is regionalization, but the extrapolation of regionalization equations beyond the limits of sample data is not recommended. This paper proposes a new method for reducing overestimation errors associated with the extrapolation of regionalization equations for low flows. The method is based on the use of a threshold value for the maximum specific low flow discharge estimated at the gauging sites that are used in the regionalization. When a specific low flow, which has been estimated using the regionalization equation, exceeds the threshold value, the low flow can be obtained by multiplying the drainage area by the threshold value. This restriction imposes a physical limit to the low flow, which reduces the error of overestimating flows in regions of extrapolation. A case study was done in the Urucuia river basin, in Brazil, and the results showed the regionalization equation to perform positively in reducing the risk of extrapolation.

Portable flow board for storage of fruits and vegetables in mini-chambers with controlled atmosphere

ABSTRACT A portable flow board system was developed in the present study with the aim to facilitate lab-scale experiments of controlled atmosphere (CA) with fruits and vegetables. This sturdy flow board combines ease fabrication, low cost and gas economy. Its functionality is provided by manifolds and gas mixers. Each gaseous component is supplied by a gas cylinder through a differential valve of adjusted pressure control, generally at 6 kPa, and forced through 13 standardized restrictors coupled to each manifold output. Controlled atmospheres are then formed with one, two or three gases in 13 gas mixers affixed to the flow board base, which are further conducted through flexible tubes to storage mini-chambers that can also be used to study metabolic consumption and production of gaseous components. The restrictors used in the flow gaseous components were manufactured from microhematocrit test-type capillary glass tubes following the hot forming method under continuous air flow. The portable flow board showed to be low cost and simple post-harvest equipment that allows preparing controlled atmospheres in open systems with stable composition and flow, in a manner similar to traditional flow boards with control of gas escape by barostats.

Pressure compensating microsprinklers using microtube as a flow controller

ABSTRACT Microsprinkler non-pressure compensating nozzles usually show water flow variation along the lateral line. This study aimed at adapting microtubes into non-compensating system of microsprinklers previous installed in the field, as a self-compensated nozzle, to improve the flow uniformity along the lateral line. Microtubes were adapted to three types of commercial microsprinklers. Tests were conducted, both in the laboratory and in field, to evaluate the microsprinkler performance at four different flows (40, 50, 60 and 70 L h-1) under pressure head range from 75 to 245 kPa. Nozzles presented coefficient of flow-rate variation (CVq) lower than 5.5% and distribution uniformity (DU) greater than 95%, which are classified as excellent. The original spatial water distribution of the microsprinkler did not change by using microtube as a nozzle. This device adapted to non-pressure compensating microsprinklers are functional and operate effectively with flows ranging up to 70 L h-1. Small variations at microsprinkler flows along the lateral line can occur, however, at random manner, which is common for pressure-compensating nozzles. Therefore, the microtube technique is able to control pressure variation in microsprinklers.

Backpressure effects on the flow-pressure relation of driplines

ABSTRACT For drip irrigation design and management, it is necessary to know the relation between flow and pressure acting on emitters. In the case of subsurface drip irrigation, the backpressure phenomenon may change the hydraulic characteristics of emitters. Thus, this study aimed at determining such relationship between flow and pressure of different driplines in surface and subsurface conditions; aiming to find possible differences in hydraulic behavior. We tested four emitter types; two pressure compensating (D5000 and Hydro PCND) and two non-pressure compensating (TalDrip and Jardiline). Emitter flow rates were attained in atmospheric conditions and submerged in water, in which submergence levels represented backpressure. Assays were performed using inlet pressures of 80, 100, 120, and 150 kPa for the Hydro PCND dripline and 25, 50, 100, and 150 kPa for the other ones; the backpressures were of 0.49, 1.47, 2.45, 4.41 and 6.37 kPa with four replications. The emitters had their proportionality constants and discharge exponents changed in submerged applications, representing backpressure effect. Non-pressure compensating emitters had their discharge exponent decreased, while in pressure compensating ones, it was increased. Backpressure reduced emitter flow rates at all evaluated pressures.

On the development of an unstructured grid solver for inert and reactive high speed flow simulations

An unstructured grid Euler solver for reactive compressible flow applications is presented. The method is implemented in a cell centered, finite volume context for unstructured triangular grids. Three different schemes for spatial discretization are implemented and analyzed. Time march is implemented in a time-split fashion with independent integrators for the flow and chemistry equations. The capability implemented is tested for inert flows in a hypersonic inlet and for inert and reactive supersonic flows over a 2-D wedge. The results of the different schemes are compared with each other and with independent calculations using a structured grid code. The strengths and the possible weaknesses of the proposed methods are discussed.

Experimental analysis of the flow between stay and guide vanes of a pump-turbine in pumping mode

An experimental study was conducted in a pump-turbine model in pumping mode, in order to characterize the flow field structure in the region between stay and guide vanes, using mainly the laser-Doppler anemometry in a two-color and back-scattered light-based system. The structure of the steady and unsteady flow was analyzed. The measurements were performed at three operation points. The obtained data provide appropriate boundary conditions and a good base of validation for numerical codes, and for the understanding of main loss mechanisms of this complex flow.