An Experimental Study of Natural Convection in a Rectangular Enclosure Using PIV Measurement Techniques

Author(s):  
K. M. Akyuzlu ◽  
J. Farkas

An experimental study is conducted to determine the circulation patterns inside a rectangular enclosure due to natural convection using a Particle Image Velocimeter (PIV). Experiments were conducted using two different fluids (air and water) and for rectangular enclosures with aspect ratios 0.5 and 1.0. Natural convection in enclosures has been experimentally studied in the past. Many of these studies cited in the literature use some kind of an optical method like interferograms, shadowgraphs, streak photographs, or multi-exposure photographs to visualize the flow patterns in the enclosure. The present study employs a commercial two-dimensional PIV to capture, instantaneously, the circulation patterns inside the test section. The test cavity in the present setup is of rectangular shape, which is 5 inches (127 mm) wide, where the height of the enclosure can be changed to obtain aspect ratios of 0.5 and 1.0. The depth of the rectangular enclosure measures 12 inches (305 mm) to minimize the effect of walls normal to the two dimensional flow patterns that are expected in this type of arrangement. The walls of the cavity are made of Aluminum plates. These plates are kept at constant but different temperatures during the experiments. In the present study, hollow glass sphere particles with 10 microns in diameter were used as seeding for water experiments and fine particles/flakes of ash generated from burned incense were used as seeding in the air experiments. For each working fluid, the experiments were repeated for different aspect ratios and for different wall temperature differences which corresponded to Rayleigh numbers in the range of 106 and 107. Velocity fields were captured at steady state for each experiment using the two-dimensional PIV system. Numerical studies were also carried out using a commercial CFD software. Comparisons of the numerical and experimental results indicate a good match in terms of circulation patterns and velocity magnitudes in the core of the buoyancy driven flow. Discrepancies in measured and predicted values of velocities are more pronounced near to the boundaries of the enclosure. Separate measurements with finer interrogation areas and different PIV setting were required to improve the accuracy of the measurements near the corners (top and bottom) of the enclosure. The results of these measurements are also presented.

Author(s):  
K. M. Akyuzlu ◽  
Y. Pavri ◽  
A. Antoniou

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 (Pr=1.0). One of the vertical walls of the enclosure is kept at a higher temperature then the opposing vertical wall. The top and the bottom of the enclosure are assumed insulated. The physics based mathematical model for this problem consists of conservation of mass, momentum (two-dimensional Navier-Stokes equations) and energy equations for the enclosed fluid subjected to appropriate boundary conditions. The working fluid is assumed to be compressible through a simple ideal gas relation. The governing equations are discretized using second order accurate central differencing for spatial derivatives and first order forward finite differencing for time derivatives where the computation domain is represented by a uniform orthogonal mesh. The resulting nonlinear equations are then linearized using Newton’s linearization method. The set of algebraic equations that result from this process are then put into a matrix form and solved using a Coupled Modified Strongly Implicit Procedure (CMSIP) for the unknowns (primitive variables) of the problem. A numerical experiment is carried out for a benchmark case (driven cavity flow) to verify the accuracy of the proposed solution procedure. Numerical experiments are then carried out using the proposed compressible flow model to simulate the development of the buoyancy driven circulation patterns for Rayleigh numbers between 103 and 105. Finally, an attempt is made to determine the effect of compressibility of the working fluid by comparing the results of the proposed model to that of models that use incompressible flow assumptions together with Boussinesq approximation.


Author(s):  
K. M. Akyuzlu ◽  
J. Farkas

An experimental study was conducted to observe/visualize, the formation of circulation patterns inside a square cavity due to the movement of a lid at constant velocity. Lid driven cavity flow is one of the benchmark studies used in the verification/improvement of CFD codes for internal flow applications/predictions. Previous work on this topic is primarily focused on improving the steady state predictions of the CFD codes using different numerical schemes and algorithms. Furthermore, almost all of the studies reported in computational fluid mechanics literature relates to steady state predictions of lid or shear driven flows. Experimental work that is reported in these studies is limited in scope and number. This paper reports on the measurements we made using Particle Image Velocimeter (PIV) technique to determine the flow field as it develops from stagnation to steady state inside a square cavity driven by a lid. For this purpose, we employed a 2-D PIV system, which uses a double-cavity, Nd:Yag laser to illuminate the test cavity. Experiments were conducted using water as the working fluid inside a square cavity that is one inch (25.4 mm) high and one inch wide. The depth of the cavity is five inches (127 mm) to ensure two-dimensional circulations patterns. Hollow glass sphere particles with 10 microns in diameter were used as seeding of the working fluid, water. Experiments were repeated for different lid velocities corresponding to lid Reynolds numbers (laminar to beginning of transition of turbulence.) Velocity fields were captured during the development of the circulations patters each being unique for the time of the measurement and value of the lid velocity. The center of the circulation pattern and its path inside the cavity is constructed from the captured images as steady state is attained. Also, the strength of the circulation (as manifested by the increase in the diameter of the circulation) is determined at different times for different Reynolds numbers.


Author(s):  
K. M. Akyuzlu ◽  
S. Nemani ◽  
K. Chakravarthy

The objective of this study is to investigate, experimentally, using Particle Image Velocimeter (PIV), the effect of different heat transfer boundary conditions on natural convection (specifically on flow patterns) inside a storage tank (a rectangular enclosure in this case with an aspect ratio of 0.5.) Purified water is used as the working fluid and it is seeded with 10 micron hollow glass sphere particles. The results of the first set of experiments are used to verify the flow patterns expected for the constant wall temperature boundary condition, the benchmark case. Similar experiments are conducted for the constant heat flux boundary condition. The Rayleigh number for all the cases studied lie between 106 and 107. The time averaged velocity field is determined using standard cross correlation techniques. Streamlines, and velocity contour plots are generated using this velocity field. Finally, the results of both cases are compared to identify the differences in circulation patterns and thermal stratification in the fluid.


Author(s):  
K. M. Akyuzlu ◽  
M. Chidurala

A two-dimensional, mathematical model is adopted to investigate the development of buoyancy driven circulation patterns and temperature stratification inside a rectangular enclosure. One of the vertical walls of the enclosure is kept at a higher temperature then the opposing vertical wall. The top and the bottom of the enclosure are assumed insulated. The physics based mathematical model for this problem consists of conservation of mass, momentum (two-dimensional, unsteady Navier-Stokes equations for turbulent compressible flows), and energy equations for the enclosed fluid subjected to appropriate boundary conditions. A standard two equation turbulence model is used to model the turbulent flow in the enclosure. The compressibility of the working fluid is represented by an ideal gas relation. The conservation equations are discretized using an implicit finite difference technique which employs second order accurate central differencing for spatial derivatives and second order (based on Taylor expansion) finite differencing for time derivatives. The linearized finite difference equations are solved using a Coupled Modified Strongly Implicit Procedure (CMSIP) for the unknowns of the problem. Numerical experiments were then carried out to simulate the development of the buoyancy driven circulation patterns inside rectangular enclosures (with aspects ratios 0.5, 1 and 1.5) filled with a compressible fluid (Pr = 0.72). Experiments were repeated for various wall temperature differences which corresponded to Rayleigh numbers between 104 and 106. Changes in unsteady circulation patterns, temperature contours, and vertical and horizontal velocity profiles were predicted while the flow inside the enclosure transferred from laminar to turbulent flow due to the sudden temperature change imposed on the vertical walls of the enclosure. Only the results of the enclosure with aspect ratio one is presented in this paper. These results indicate that this transition is characterized by unicellular circulation patterns breaking up in to multicellular formations and increase in the values of the predicted wall heat fluxes and Nusselt number as flow becomes turbulent.


1991 ◽  
Vol 113 (1) ◽  
pp. 116-121 ◽  
Author(s):  
H. L. Zhang ◽  
Q. J. Wu ◽  
W. Q. Tao

In this paper, the results of an experimental study of laminar natural convection heat transfer and fluid flow in horizontal annuli between a cylindrical envelope and its inner concentric octagonal heated cylinder are presented. Two octagonal cylinders are investigated: one with a complete surface and the other with two horizontal slots on the top and bottom surfaces. The ratio of the slot width W to H is 0.072. Air is used as the working fluid. The range of Rayleigh number is 2.1×102–1.58×106 for the unslotted case and 1.2×102–1.5×106 for the slotted case. The average heat transfer correlations for the two cases are provided. The results show that the heat transfer intensity of the unslotted octagon is slightly weaker than that in a cylindrical annulus, while for the slotted case, the overall heat transfer enhancement may be as high as 74 percent. The smoke technique is used to visualize the flow patterns. A series of photographs of the flow patterns are provided, which enhances our understanding of the mechanism of heat transfer enhancement for the slotted octagonal case.


1994 ◽  
Vol 116 (3) ◽  
pp. 627-632 ◽  
Author(s):  
Y. Kamotani ◽  
F.-B. Weng ◽  
S. Ostrach ◽  
J. Platt

An experimental study is made of natural convection oscillations in gallium melts enclosed by right circular cylinders with differentially heated end walls. Cases heated from below are examined for angles of inclination (φ) ranging from 0 deg (vertical) to 75 deg with aspect ratios Ar (height/diameter) of 2, 3, and 4. Temperature measurements are made along the circumference of the cylinder to detect the oscillations, from which the oscillatory flow structures are inferred. The critical Rayleigh numbers and oscillation frequencies are determined. For Ar=3 and φ = 0 deg, 30 deg the supercritical flow structures are discussed in detail.


2016 ◽  
Vol 138 (9) ◽  
Author(s):  
Amin Bouraoui ◽  
Rachid Bessaïh

In this paper, a numerical study of three-dimensional (3D) natural convection air-cooling of two identical heat sources, simulating electronic components, mounted in a rectangular enclosure was carried out. The governing equations were solved by using the finite-volume method based on the SIMPLER algorithm. The effects of Rayleigh number Ra, spacing between heat sources d, and aspect ratios Ax in x-direction (horizontal) and Az in z-direction (transversal) of the enclosure on heat transfer were investigated. In steady state, when d is increased, the heat transfer is more important than when the aspect ratios Ax and Az are reduced. In oscillatory state, the critical Rayleigh numbers Racr for different values of spacing between heat sources and their aspect ratios, at which the flow becomes time dependent, were obtained. Results show a strong relation between heat transfers, buoyant flow, and boundary layer. In addition, the heat transfer is more important at the edge of each face of heat sources than at the center region.


Sign in / Sign up

Export Citation Format

Share Document