No description available.

This thesis presents an experimental and computational fluid dynamics (CFD) study on a rectangular S-bend with straight and diffusing passages with passive effusion cooling. Experimental tests were performed at both cold and hot flow conditions over a range of Reynolds numbers from 2.5e5 to 4.5e5. Hot flow testing was conducted with the primary flow temperature up to 300 °C. Severe backpressure penalties occurred with full-surface passive effusion injection in cold flow tests. Moderate penalties occurred with reduced surface coverage whereby the performance was affected by the S-bend secondary fields with injection at different locations. High surface cooling effectiveness with full-coverage of cooling film was measured; the impacts from the S-bend secondary flow fields were measured to be minimal. The CFD study revealed the importance of using experimental flow boundary conditions for simulations. Using the standard k-ε model with wall functions was confirmed as appropriate for simulating the S-bend flow with effusion cooling. A coarse-grid CFD methodology using a porous wall boundary condition to simulate the effects of effusion cooling was investigated. From a design perspective, this model is preferable for quantifying the injection flow rate since the actual mass flow rate is not known. Comparison to the alternative solution using uniform mass flow boundary conditions showed that both models incorrectly predicted the momentum. The porous wall model, however, is promising for practical design applications of S-duct flow fields with effusion injections. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2013-12-23 14:20:32.38

Maintaining IAQ (Internal Air Quality) and thermal comfort of occupants in buildings have been a challenge to overcome satisfying the two ends: criteria for sustainability and cost effectiveness. Although there was a movement for mechanical ventilation systems in the recent past, in addition to the cost involved, they are found to not deliver the desired air quality, lead to social consequences such as sick building syndrome, contribute to environmental consequences related to ozone-depleting substances with increasing energy consumption, generate noise and having difficulties in cleaning and maintaining. These consequences compelled research on natural ventilation systems, which were used in ancient buildings. Although it has been found that natural ventilation of buildings can become a substantial architectural design tool that leads to “breathing architecture,” fluctuations in indoor temperature and air quality makes depending entirely on natural ventilation less effective. The combination of natural and mechanical ventilation, the hybrid ventilation or mixed-mode ventilation, systems utilizes advantages and eliminates drawbacks from both mechanical and entirely dependent natural ventilation systems. Hybrid ventilation systems, which have been utilized in historical buildings, with less investment cost and reduction of energy usage have been found to be a solution to provide acceptable standards of IAQ and thermal comfort through natural air circulation in buildings. This research study was carried out to verify the effectiveness of a hybrid ventilation system in an auditorium built around 60 years back for its effectiveness as a provider of thermal comfort to its occupants. Computational Fluid Dynamic (CFD) modeling was carried out on a Finite Element (FE) model owing to its capability of offering a wide range of flexible analytical solutions, lower realization time and comparative cost effectiveness to experimental methods of modeling. This verification of the system has revealed that hybrid ventilation systems could provide effective thermal comfort in buildings designed specifically to allow circulation of air through the system. The results of the study were in agreement with measured data and the expected flow of air through the building when the thermal load due to metabolism of occupants was not included in the analysis. In addition, the expected results complied with similar studies on natural/hybrid ventilation systems. With the addition of the thermal load, as a uniform heat flux from the flow of the auditorium, it was observed that the conditioning of the air throughout the space was better than the without thermal load scenario. In the case modeling people as cylinders, with a convective heat flux, it was observed that the air flow direction changes and the seating level of the auditorium do not get sufficient air flow to maintain a comfortable air quality.  Ineffective simulation of the inlet louver was assumed to be the primary reason for this scenario and other reasons such as the seating arrangement modeling too could have effects on the result. As conclusions of the study it was found that the whole building system properties have to be selected, as the control component to produce operating commands, to circulate air through the building in accordance with the air flow: both velocity and patterns, required to maintain thermal comfort of all occupants. Air inflow could be through windows as acquisition components to collect indoor and outdoor climatic parameters and air outflow could be mechanically controlled through exhausted fans turning on or off as the operating component in the system. The result of the study ensures the method of solutions through CFD to be utilized to provide effective and less costly path to verify systems such as natural or hybrid air flow systems through buildings.  The whole system studied could be applied with suitable contextual modifications to any new location, with similar cost effective modeling, to produce less fuel consuming building systems leading to sustainability of built environment.

Fluidization of solid particles using gas flow is an important process in chemical and pharmaceutical industries. The dynamics of fluidisation are intricately related to particle scale physics. Fluid-particle interactions dominate gas-solid fluidization behaviour for particles with average size and density greater than 10-4 m and 103 kg/m3, respectively, classified as Geldart B and D particles. Inter-particle forces, such as cohesion, play an increasingly important role in the fluidization dynamics of smaller particles, which are classified as Geldart A and C. In particular, interesting fluidization regimes have been noticed for weakly cohesive Geldart A particles, exhibiting a window of uniform fluidization before the onset of bubbling behaviour. Despite widespread industrial interests, the fundamental understanding of the mechanisms that underlie these fluidization regimes is poor. The present study aims to improve the understanding of fluidization dynamics of Geldart A regimes using numerical simulations. A DEM-CFD model was employed to capture the widely separated spatial and temporal scales associated with fluidization behaviour. The model couples the locally averaged Navier-Stokes equation for fluid with a discrete description of the particles. The methodology and its computer implementation are verified and validated to assess the extent of fluidization physics that it is able to capture. Verification cases check the implementation of the inter-phase momentum transfer term, drag model implementation and pressure-velocity coupling. The test cases are employed in order to cover a wide range of flow conditions. Robust validation tests for complex fluidization phenomena such as bubbling, spouting and bidisperse beds have been conducted to assess the predictive capabilities of the DEM-CFD solver. The simulation results for time and spatially averaged fluidziation behaviour are compared to experimental measurements obtained from the literature, and are shown to have capture fluidization physics qualitatively. Robust features of bubbling fluidization, such as minimum fluidization velocity, frequency of pressure drop fluctuations, segregation rates and solid circulation patterns were captured. Furthermore, the DEM-CFD model is critically assessed in terms of model conceptualization and parameter estimation, including those for drag closures, particle-wall boundary conditions, bed height and particle shape effects. The validation studies establish modelling best-practice guidelines and the level of discrepancy against the analytical solutions or experimental measurements. Having developed the model and established its predictive capability, it is used to probe the hydrodynamics of weakly cohesive particles. Cohesive interactions are captured by employing a pair-wise van derWaals force model. The cohesive strength of the granular bed is quantified by the ratio of the maximum van der Waals force to the particle gravitational force, defined as the granular Bond number. The Bond number of the bed is increased systematically from 0-10 to examine the role of cohesion in the fluidization behaviour of fine powders while keeping the particle size and density constant across all the simulations. The idea was to segregate the hydrodynamics associated with size and density of the particles from the inter-particle interactions. The size and density of the particles are carefully chosen at a scale where inter-particle forces are present but minimal [Seville et al., 2000]. The Geldart A fluidization behaviour is captured for granular beds with Bond numbers ranging from 1 to 3. Many robust features of Geldart A fluidization, such as pressure drop overshoot, delay in the onset of bubbling, macroscopic Umf predictions and uniform bed expansion are captured in the DEM-CFD framework. The expanded bed was characterized according to criteria that the particles are highly immobile in this regime and the expanded porosity is related to inlet velocity by Richardson–Zaki correlations. Sudden jumps in the magnitudes of global granular temperature were found near the regime transitions. This observation was used an indicator of the onset of bubbling and quantification of minimum bubbling velocity (Umb). The window of the expanded bed regime (quantified as Umb - Umf) was shown to be an increasing function of cohesive strength of the bed. Furthermore, the stability of the expanded bed was probed by studying the response of the expanded bed to sudden inertial and voidage shocks. A kinematic wave, generated as a response to the voidage shock, was shown to slow down with increasing cohesion and decreasing hydrodynamic forces. Furthermore, predictions of Umb by DEM-CFD simulations for weakly cohesive beds were compared against empirical correlations by Valverde [2013] with an excellent match. Stress analysis of the expanded bed revealed the presence of tensile stresses. As the inlet velocity is increased beyond the minimum fluidization velocity, a longitudinal shift of these negative stresses is observed until they reach the top of the bed. Negative stresses were seen at the bed surface at the onset of bubbling. The role of cohesion stresses in the formation of expanded bed and suppression of bubbling was highlighted. Finally, the microstructure of the expanded bed was probed at different local micro and mescoscopic length scales. Evidence of clustering, agglomeration and cavities were presented in the expanded bed. Expanded bed expansion was shown to have mesostructural inhomogeneities present, which is contrary to the belief of homogeneous expansion.

This study compared experimental and Computational Fluid Dynamics (CFD) results for heating and cooling in a packed bed (N=5.33). The experimental data was compared between heating and cooling, and was also used to validate the CFD model. The validated models were used to compare theoretical heat transfer parameters. For the experiments, it was found that the effective thermal conductivity was comparable for heating and cooling, and the wall Nusselt number for heating was higher. For the CFD results, it was found that both the wall Nusselt number and effective thermal conductivity were comparable for heating and cooling. The wall Nusselt number was slightly higher for cooling, however this difference decreased as the Reynolds number increased.

This study compared experimental and Computational Fluid Dynamics (CFD) results for heating and cooling in a packed bed (N=5.33). The experimental data was compared between heating and cooling, and was also used to validate the CFD model. The validated models were used to compare theoretical heat transfer parameters. For the experiments, it was found that the effective thermal conductivity was comparable for heating and cooling, and the wall Nusselt number for heating was higher. For the CFD results, it was found that both the wall Nusselt number and effective thermal conductivity were comparable for heating and cooling. The wall Nusselt number was slightly higher for cooling, however this difference decreased as the Reynolds number increased.

The ingestion of hot mainstream gas into the wheel-space between the rotor and staler discs is one of the most important internal cooling problems for gas turbine designers. To solve this problem, engineers design a rim seal at the periphery of wheel-space and direct a sealing flow from the internal cooling system to prevent ingress. The main aim of this thesis is to build a simple computational model to predict the scaling effectiveness of externally-induced ingress for engine designers. The axisymmetric model represents a gas turbine wheel-space and provides useful information related to the fluid dynamics and heat transfer in the wheel-space. At the same time, this model saves much computation time and cost for engine designers who currently use complex and time-consuming 3D models. The- computational model in this -thesis is called the prescribed ingestion model. Steady simulations are carried out using the commercial CFD code, ANSYS CFX with meshes built using ICEM CFD. Boundary conditions are applied at the ingress inlet of the model using experimental measurements and a mass-based averaging procedure. Computational parameters such as rotational Reynolds number, non-dimensional sealing flow rate and thermal conditions on the rotor are selected to investigate the fluid dynamics and heat transfer at typical experimental rig operating conditions. Different rim seal geometries arc investigated and results are compared with experimental data. In addition to the prescribed ingestion model, two typical axisymmetric rotor-stator system models without ingress arc established. The aim of these rotor-stator models is to investigate the fluid dynamics and heat transfer of the wheel-space in the situation without ingress. The effects of geometry and turbulence model also arc studied in these simulations. Most results from these simulations are in good agreement with experimental data from the literature, which enhances confidence in the prescribed Ingestion model.

"Liquid piston gas compression utilizes a liquid to directly compress gas. The benefit of this approach is that liquid can conform to irregular compression chamber volume. The compression chamber is divided into many small little bores in order to increases the surface area to volume ratio. The heat transfer rate increases with increasing surface area to volume ratio. However, as the bore diameter becomes smaller, the viscous force increases. In order to maximize the heat transfer rate and to minimize the viscous force, computational fluid dynamics is used. ANSYS Fluent is used to simulate the liquid piston gas compression cycle. Having created the model in Fluent, different factors, including diameter, length, liquid temperature, and the acceleration are varied in order to understand how each factor affects the heat transfer and viscous energy loss. The results show that both viscous force and heat transfer rate increase as the diameter decreases. The viscous force increases and the heat transfer decreases as the length increases. Both the viscous force and heat transfer increase as the acceleration increases. The viscous force decreases as the liquid temperature increases. Results show that the highest compression efficiency of 86.4% is found with a 3mm bore radius and a short cylinder. The piston acceleration is advised to be below 0.5g in order to avoid surface instability problem."

Pressure screens are used for contaminant removal and fibre length fractionation in the production of pulp and paper products. Axial variations and entry effects in the screen are known to occur and these variations have not been adequately quantified. This thesis describes a fundamental study of the axial variations of several factors that occur within an industrial pressure screen; namely, pulp consistency, fibre length distribution, rotor pressure pulse, and feed annulus tangential velocity. Axial variations of pulp consistency in the screen annulus and the accept chamber of the screen were studied using an internal radial sampling method. Localised pulp samples were taken and evaluated and common measures of screen performance such as fibre passage ratio and fractionation efficiency were calculated along the screen. Consistency generally increased along the length of the screen although under certain conditions the consistency toward the front of the screen was lower than the feed consistency. A two passage ratio model that incorporated forward and reverse passage ratio was derived to elucidate the flow of both fibre and fluid through the screen and their effects on overall screen performance. The passage of fibre through the screen decreased with screen length which generally had a positive effect on the fractionation efficiency toward the back of the screen. The passage of individual fibre length fractions was also studied and it was found that long fibre had a much lower passage than short fibre which caused the average fibre length in the annulus to increase. Rotor induced pressure pulse variations along the screen length were also investigated. The magnitude of the pressure pulse was significantly lower (up to 40 %) at the rear of the screen. The variation in pressure caused by the rotor is due to a Venturi effect and the shape of the rotor. The relative velocity of the fluid and the rotor, called the slip factor, also directly affects the size of the pressure pulse in the annulus. The slip factor decreases along the length of the screen due to the increase in tangential velocity of the fluid. Pressure pulse data was also used to estimate the instantaneous aperture velocity and back-flush ratio. The instantaneous aperture velocity was calculated to vary considerably from the superficial aperture velocity by up to 5 m/s in the forward direction and 10 m/s in the reverse direction. Computational Fluid Dynamics (CFD) was used to model tangential velocity changes in simplified screen annuli with axial through flow. For a smooth screen rotor the mean tangential velocity increased over the entire length of the annulus without reaching a maximum value. A step and bump rotor were modelled and the shape of the pressure pulses showed good agreement with experimentally measured pulses. The mean tangential velocity and the entrance length were found to be heavily dependant on the screen rotor used.

Flow velocity field, vorticity and circulation and wall shear stresses were simulated by FSI approach under conditions of pulsatile flow in a scale model of the rabbit elastin-induced aneurysm. The flow pattern inside the aneurysm sac confirmed the in vitro experimental findings that in diastole time period the flow inside the aneurysm sac is a stable circular clock-wise flow, while in systole time period higher velocity enters into the aneurysm sac and during systole and diastole time period an anti-clock circular flow pattern emerged near the distal neck; in the 3-D aneurysm sac, the kinetic energy per point is about 0.0002 (m2/s2); while in the symmetrical plane of the aneurysm sac, the kinetic energy per point is about 0.00024 (m2/s2). In one cycle, the shape of the intraaneurysmal energy profile is in agreement with the experimental data; The shear stress near the proximal neck experienced higher shear stress (peak value 0.35 Pa) than the distal neck (peak value 0.2 Pa), while in the aneurysm dome, the shear stress is always the lowest (0.0065 Pa). The ratio of shear stresses in the proximal neck vs. distal neck is around 1.75, similar to the experimental findings that the wall shear rate ratio of proximal neck vs. distal neck is 1.5 to 2.