A two-dimensional mathematical model is adopted to investigate the development of buoyancy driven circulation patterns and temperature contours inside a rectangular enclosure filled with a compressible fluid where one of the vertical walls of the enclosure is kept at a higher temperature than the opposite one. Fluid thermodynamic and transport properties are assumed to be functions of temperature. The governing equations are discretized using second order accurate differencing for spatial and temporal derivatives and then linearized using Newton's linearization method. The resulting set of algebraic equations is solved using the Coupled Modified Strongly Implicit Procedure for the unknowns of the problem. The results of this study show that the variable property model predicts lower values for wall heat fluxes and Nu number than the constant property one for Rayleigh numbers between 104 and 105.

Natural ventilation is an effective method to save energy required to condition buildings and to improve indoor air quality. Computational fluid dynamics (CFD) was used to model single-sided buoyancy-driven natural ventilation in a single room with a heater and door. The velocity and temperature profiles at the doorway agreed fairly well with published literature that includes Mahajan's experimental [2] and Schaelin et al's numerical studies [1]. The 2D and 3D models predicted the neutral level with a difference of 5.6 % and 0.08 % compared to the experimental results, respectively. Using solutions at the doorway, heat transfer rates were calculated. More realistic situations were studied considering conduction, various ambient conditions, wind speeds, and additional heat sources and furniture in the room. The heat loss through the wall was modeled and the airflow and temperature within the room showed no significant changes despite modeling conduction through the walls. Various ambient temperatures and wind speeds were tested, and the neutral level height and total heat transfer rate through the doorway increased with decreasing ambient temperatures. However, the neutral level did not significantly change as wind speeds varied. Total heat transfer rate at the doorway became positive, that is heat transferred into the room, with wind speed. Lastly, the effect of additional heat sources (mini-refrigerator, monitor and computer) and furniture (bookshelf, desk, chair and box) on airflow and heat transfer in the room was analyzed by comparing with a simple case of a room with a heater. Large velocities and high temperatures were predicted in the vicinity of the heat sources. However, the spatially averaged velocity and temperature did not change significantly despite additional heat sources. The room with furniture was modeled at lower ambient temperature, where the spatially averaged velocities were larger and temperatures were lower than the simple case. The room heated up and reached its thermal comfort level, but the velocities exceeded the maximum acceptable level set by ASHRAE guidelines [8]. Wind was considered simultaneously with the lower temperature, and the room was cooled faster with wind. However, the room was never able to achieve the comfortable level both in velocity and temperature. / Master of Science

No description available.

Over the last two centuries, there has been a continual evolution of how occupied rooms are heated, with inventors competing to design new heating devices. In particular, there is a wide range of convector types, which vary in shape, size, design, material, operating medium and application. With approximately 190 million convectors installed in the UK alone, the question arises regarding the dependencies on the efficiency of heat distribution through convector-induced filling flows. A standard approach to evaluate convector performance is based on the convector strength only, the implication being the stronger the convector the better the performance. This work has gone beyond the limits of a stereotypical assessment in pursuit of answers regarding the physics of convector-induced filling and a new objective method to evaluate the efficiency of this transient process. The ultimate goal has been to provide a deep understanding of filling and stratification induced by a convector, in order to heat rooms rapidly and effectively. An experimental facility has been designed that approximates dynamic similarity between the experimental set-up and a real-life room with a convector. In the experiments, a rectangular sectioned water tank represents a room and a saline source rectangular sectioned panel with sintered side walls provides a convector representation. Experiments have been performed in water with a saline solution to ensure high Rayleigh numbers. Diagnostic techniques involve a combination of a shadowgraph method, a dye-attenuation method, direct salinity measurements and a new application of Particle Image Velocimetry (PIV). Interesting insight into convector-induced buoyancy-driven flows has been gained. As a result, new guidelines aimed at heating rooms more rapidly and effectively have been proposed. The key outcome that can be immediately applied is that, for a given convector strength, heat distribution with height can be improved by adjusting the convector position. For instance, faster filling leading to more uniform heat distribution occurs in rooms with convectors detached from side walls, due to large-scale mixing flows in the early period of filling. Also shorter convectors relative to the room height, positioned close to the floor level, promote faster and more uniform filling. An attempt to describe the transient filling has been made and to do so statistical methods, application specific, have been developed. As a result, the empirical equations describing both the filling rates in different stages of filling and the development of stratification have been derived, which rank the governing parameters, based on their importance, as either dominant or subordinate. Two dominant parameters governing filling flows are the non-dimensional accumulation parameter B and the Rayleigh number ΔRa, which are related to the convector strength. The impact of these two parameters is constant throughout the process. The parameters accounting for the system geometry and filling time (T) are subordinate parameters. Their impact, visible in the early period, decreases as filling continues.