Assessment of turbulence model predictions for a centrifugal compressor simulation

Steady-state computational fluid dynamics (CFD) simulations are an essential tool in the design process of centrifugal compressors. Whilst global parameters, such as pressure ratio and efficiency, can be predicted with reasonable accuracy, the accurate prediction of detailed compressor flow fields is a much more significant challenge. Much of the inaccuracy is associated with the incorrect selection of turbulence model. The need for a quick turnaround in simulations during the design optimisation process, also demands that the turbulence model selected be robust and numerically stable with short simulation times.
In order to assess the accuracy of a number of turbulence model predictions, the current study used an exemplar open CFD test case, the centrifugal compressor ‘Radiver’, to compare the results of three eddy viscosity models and two Reynolds stress type models. The turbulence models investigated in this study were (i) Spalart-Allmaras (SA) model, (ii) the Shear Stress Transport (SST) model, (iii) a modification to the SST model denoted the SST-curvature correction (SST-CC), (iv) Reynolds stress model of Speziale, Sarkar and Gatski (RSM-SSG), and (v) the turbulence frequency formulated Reynolds stress model (RSM-ω). Each was found to be in good agreement with the experiments (below 2% discrepancy), with respect to total-to-total parameters at three different operating conditions. However, for the off-design conditions, local flow field differences were observed between the models, with the SA model showing particularly poor prediction of local flow structures. The SST-CC showed better prediction of curved rotating flows in the impeller. The RSM-ω was better for the wake and separated flow in the diffuser. The SST model showed reasonably stable, robust and time efficient capability to predict global and local flow features.

Mixing performance evaluation of additive manufactured milli-scale reactors

The mixing performance of three passive milli-scale reactors with different geometries was investigated at different Reynolds numbers. The effects of design and operating characteristics such as mixing channel shape and volume flow rate were investigated. The main objective of this work was to demonstrate a process design method that uses on Computational Fluid Dynamics (CFD) for modeling and Additive Manufacturing (AM) technology for manufacture. The reactors were designed and simulated using SolidWorks and Fluent 15.0 software, respectively. Manufacturing of the devices was performed with an EOS M-series AM system. Step response experiments with distilled Millipore water and sodium hydroxide solution provided time-dependent concentration profiles. Villermaux-Dushman reaction experiments were also conducted for additional verification of CFD results and for mixing efficiency evaluation of the different geometries. Time-dependent concentration data and reaction evaluation showed that the performance of the AM-manufactured reactors matched the CFD results reasonably well. The proposed design method allows the implementation of new and innovative solutions, especially in the process design phase, for industrial scale reactor technologies. In addition, rapid implementation is another advantage due to the virtual flow design and due to the fast manufacturing which uses the same geometric file formats.

Investigation expérimentale du décollement dans l’aspirateur d’une turbine bulbe

La présente thèse propose une étude expérimentale du décollement dans le diffuseur d’un modèle de turbine hydroélectrique bulbe. Le décollement se produit quand la turbine est opérée à forte charge et il réduit la section effective de récupération du diffuseur. La diminution de la performance du diffuseur à forte charge engendre une baisse brusque de l’efficacité de la turbine et de la puissance extraite. Le modèle réduit de bulbe est fidèle aux machines modernes avec un diffuseur particulièrement divergent. Les performances de la turbine sont mesurées sur une large gamme de points d’opération pour déterminer les conditions les plus intéressantes pour l’étude du décollement et pour étudier la distribution paramétrique de ce phénomène. La pression est mesurée le long de l’aspirateur par des capteurs dynamiques affleurants alors que les champs de vitesse dans la zone de décollement sont mesurés avec une méthode PIV à deux composantes. Les observations à la paroi sont pour leur part faites à l’aide de brins de laine. Pour un débit suffisant, le gradient de pression adverse induit par la géométrie du diffuseur affaiblit suffisamment la couche limite, entraînant ainsi l’éjection de fluide de la paroi le long d’une large enveloppe tridimensionelle. Le décollement instationnaire tridimensionnel se situe dans la même zone du diffuseur indépendamment du point d’opération. L’augmentation du débit provoque à la fois une extension de la zone de décollement et une augmentation de l’occurrence de ses manifestations. La position et la forme du front de décollement fluctue significativement sans périodicité. L’analyse topologique et celle des tourbillons des champs de vitesse instantanés montrent une topologie du front de décollement complexe qui diffère beaucoup d’une réalisation à l’autre. Bien que l’écoulement soit turbulent, les tourbillons associés aux foyers du front sont clairement plus gros et plus intenses que ceux de la turbulence. Cela suggère que le mécanisme d’enroulement menant aux tourbillons du décollement est clairement distinct des mécanismes de la turbulence.

Development of a model to predict the life of pneumatic conveyor bends subject to erosive wear

The erosion processes resulting from flow of fluids (gas-solid or liquid-solid) are encountered in nature and many industrial processes. The common feature of these erosion processes is the interaction of the fluid (particle) with its boundary thus resulting in the loss of material from the surface. This type of erosion in detrimental to the equipment used in pneumatic conveying systems. The puncture of pneumatic conveyor bends in industry causes several problems. Some of which are: (1) Escape of the conveyed product causing health and dust hazard; (2) Repairing and cleaning up after punctures necessitates shutting down conveyors, which will affect the operation of the plant, thus reducing profitability. The most common occurrence of process failure in pneumatic conveying systems is when pipe sections at the bends wear away and puncture. The reason for this is particles of varying speed, shape, size and material properties strike the bend wall with greater intensity than in straight sections of the pipe. Currently available models for predicting the lifetime of bends are inaccurate (over predict by 80%. The provision of an accurate predictive method would lead to improvements in the structure of the planned maintenance programmes of processes, thus reducing unplanned shutdowns and ultimately the downtime costs associated with these unplanned shutdowns. This is the main motivation behind the current research. The paper reports on two aspects of the first phases of the study-undertaken for the current project. These are (1) Development and implementation; and (2) Testing of the modelling environment. The model framework encompasses Computational Fluid Dynamics (CFD) related engineering tools, based on Eulerian (gas) and Lagrangian (particle) approaches to represent the two distinct conveyed phases, to predict the lifetime of conveyor bends. The method attempts to account for the effect of erosion on the pipe wall via particle impacts, taking into account the angle of attack, impact velocity, shape/size and material properties of the wall and conveyed material, within a CFD framework. Only a handful of researchers use CFD as the basis of predicting the particle motion, see for example [1-4] . It is hoped that this would lead to more realistic predictions of the wear profile. Results, for two, three-dimensional test cases using the commercially available CFD PHOENICS are presented. These are reported in relation to the impact intensity and sensitivity to the inlet particle distributions.

Análise de métodos para aplicação de ventilação natural em projetos de edificações em Natal-RN

Natural ventilation is an efficient bioclimatic strategy, one that provides thermal comfort, healthful and cooling to the edification. However, the disregard for quality environment, the uncertainties involved in the phenomenon and the popularization of artificial climate systems are held as an excuse for those who neglect the benefits of passive cooling. The unfamiliarity with the concept may be lessened if ventilation is observed in every step of the project, especially in the initial phase in which decisions bear a great impact in the construction process. The tools available in order to quantify the impact of projected decisions consist basically of the renovation rate calculations or computer simulations of fluids, commonly dubbed CFD, which stands for Computational Fluid Dynamics , both somewhat apart from the project s execution and unable to adapt for use in parametric studies. Thus, we chose to verify, through computer simulation, the representativeness of the results with a method of simplified air reconditioning rate calculation, as well as making it more compatible with the questions relevant to the first phases of the project s process. The case object consists of a model resulting from the recommendations of the Código de Obras de Natal/ RN, customized according to the NBR 15220. The study has shown the complexity in aggregating a CFD tool to the process and the need for a method capable of generating data at the compatible rate to the flow of ideas and are discarded during the project s development. At the end of our study, we discuss the necessary concessions for the realization of simulations, the applicability and the limitations of both the tools used and the method adopted, as well as the representativeness of the results obtained

Análise de componentes arquitetônicos para potencialização da ventilação natural com ênfase em captadores de vento

The building envelope is the principal mean of interaction between indoors and environment, with direct influence on thermal and energy performance of the building. By intervening in the envelope, with the proposal of specific architectural elements, it is possible to promote the use of passive strategies of conditioning, such as natural ventilation. The cross ventilation is recommended by the NBR 15220-3 as the bioclimatic main strategy for the hot and humid climate of Natal/RN, offering among other benefits, the thermal comfort of occupants. The analysis tools of natural ventilation, on the other hand, cover a variety of techniques, from the simplified calculation methods to computer fluid dynamics, whose limitations are discussed in several papers, but without detailing the problems encountered. In this sense, the present study aims to evaluate the potential of wind catchers, envelope elements used to increase natural ventilation in the building, through CFD simplified simulation. Moreover, it seeks to quantify the limitations encountered during the analysis. For this, the procedure adopted to evaluate the elements implementation and efficiency was the CFD simulation, abbreviation for Computer Fluid Dynamics, with the software DesignBuilder CFD. It was defined a base case, where wind catchers were added with various settings, to compare them with each other and appreciate the differences in flows and air speeds encountered. Initially there has been done sensitivity tests for familiarization with the software and observe simulation patterns, mapping the settings used and simulation time for each case simulated. The results show the limitations encountered during the simulation process, as well as an overview of the efficiency and potential of wind catchers, with the increase of ventilation with the use of catchers, differences in air flow patterns and significant increase in air speeds indoors, besides changes found due to different element geometries. It is considered that the software used can help designers during preliminary analysis in the early stages of design

Shock Train/Boundary-Layer Interaction in Rectangular Scramjet Isolators

Numerous studies of the dual-mode scramjet isolator, a critical component in preventing inlet unstart and/or vehicle loss by containing a collection of flow disturbances called a shock train, have been performed since the dual-mode propulsion cycle was introduced in the 1960s. Low momentum corner flow and other three-dimensional effects inherent to rectangular isolators have, however, been largely ignored in experimental studies of the boundary layer separation driven isolator shock train dynamics. Furthermore, the use of two dimensional diagnostic techniques in past works, be it single-perspective line-of-sight schlieren/shadowgraphy or single axis wall pressure measurements, have been unable to resolve the three-dimensional flow features inside the rectangular isolator. These flow characteristics need to be thoroughly understood if robust dual-mode scramjet designs are to be fielded. The work presented in this thesis is focused on experimentally analyzing shock train/boundary layer interactions from multiple perspectives in aspect ratio 1.0, 3.0, and 6.0 rectangular isolators with inflow Mach numbers ranging from 2.4 to 2.7. Secondary steady-state Computational Fluid Dynamics studies are performed to compare to the experimental results and to provide additional perspectives of the flow field. Specific issues that remain unresolved after decades of isolator shock train studies that are addressed in this work include the three-dimensional formation of the isolator shock train front, the spatial and temporal low momentum corner flow separation scales, the transient behavior of shock train/boundary layer interaction at specific coordinates along the isolator's lateral axis, and effects of the rectangular geometry on semi-empirical relations for shock train length prediction. A novel multiplane shadowgraph technique is developed to resolve the structure of the shock train along both the minor and major duct axis simultaneously. It is shown that the shock train front is of a hybrid oblique/normal nature. Initial low momentum corner flow separation spawns the formation of oblique shock planes which interact and proceed toward the center flow region, becoming more normal in the process. The hybrid structure becomes more two-dimensional as aspect ratio is increased but corner flow separation precedes center flow separation on the order of 1 duct height for all aspect ratios considered. Additional instantaneous oil flow surface visualization shows the symmetry of the three-dimensional shock train front around the lower wall centerline. Quantitative synthetic schlieren visualization shows the density gradient magnitude approximately double between the corner oblique and center flow normal structures. Fast response pressure measurements acquired near the corner region of the duct show preliminary separation in the outer regions preceding centerline separation on the order of 2 seconds. Non-intrusive Focusing Schlieren Deflectometry Velocimeter measurements reveal that both shock train oscillation frequency and velocity component decrease as measurements are taken away from centerline and towards the side-wall region, along with confirming the more two dimensional shock train front approximation for higher aspect ratios. An updated modification to Waltrup \& Billig's original semi-empirical shock train length relation for circular ducts based on centerline pressure measurements is introduced to account for rectangular isolator aspect ratio, upstream corner separation length scale, and major- and minor-axis boundary layer momentum thickness asymmetry. The latter is derived both experimentally and computationally and it is shown that the major-axis (side-wall) boundary layer has lower momentum thickness compared to the minor-axis (nozzle bounded) boundary layer, making it more separable. Furthermore, it is shown that the updated correlation drastically improves shock train length prediction capabilities in higher aspect ratio isolators. This thesis suggests that performance analysis of rectangular confined supersonic flow fields can no longer be based on observations and measurements obtained along a single axis alone. Knowledge gained by the work performed in this study will allow for the development of more robust shock train leading edge detection techniques and isolator designs which can greatly mitigate the risk of inlet unstart and/or vehicle loss in flight.

Simulação computacional como ferramenta de auxílio ao projeto: aplicação em edifícios naturalmente ventilados no clima de Natal/RN

Natural air ventilation is the most import passive strategy to provide thermal comfort in hot and humid climates and a significant low energy strategy. However, the natural ventilated building requires more attention with the architectural design than a conventional building with air conditioning systems, and the results are less reliable. Therefore, this thesis focuses on softwares and methods to predict the natural ventilation performance from the point of view of the architect, with limited resource and knowledge of fluid mechanics. A typical prefabricated building was modelled due to its simplified geometry, low cost and occurrence at the local campus. Firstly, the study emphasized the use of computational fluid dynamics (CFD) software, to simulate the air flow outside and inside the building. A series of approaches were developed to make the simulations possible, compromising the results fidelity. Secondly, the results of CFD simulations were used as the input of an energy tool, to simulate the thermal performance under different rates of air renew. Thirdly, the results of temperature were assessed in terms of thermal comfort. Complementary simulations were carried out to detail the analyses. The results show the potentialities of these tools. However the discussions concerning the simplifications of the approaches, the limitations of the tools and the level of knowledge of the average architect are the major contribution of this study

Performance assessment of a vertical axis turbine in a marine current flume tank and CFD modelling

As advances in numerical modelling techniques support the increased confidence in predictions from computer simulations, the need remains to have experimental verification built into the design process. This paper outlines the experimental investigation carried out on a shielded vertical axis turbine in a marine environment. The experiments consist of performance measurements and the use of particle image velocimetry on a small scale device in a marine current flume. The results demonstrate that the performance of the device can be modelled numerically; in particular, the results show that the numerical model used can correctly predict the increase in performance with Reynolds number.

Simulations of droplet dispersion in the wakes of cylinders and icing tunnel spray bars

Simulations of droplet dispersion behind cylinder wakes and downstream of icing tunnel spray bars were conducted. In both cases, a range of droplet sizes were investigated numerically with a Lagrangian particle trajectory approach while the turbulent air flow was investigated with a hybrid Reynolds-Averaged Navier-Stokes/Large-Eddy Simulations approach scheme. In the first study, droplets were injected downstream of a cylinder at sub-critical conditions (i.e. with laminar boundary layer separation). A stochastic continuous random walk (CRW) turbulence model was used to capture the effects of sub-grid turbulence. Small inertia droplets (characterized by small Stokes numbers) were affected by both the large-scale and small-scale vortex structures and closely followed the air flow, while exhibiting a dispersion consistent with that of a scalar flow field. Droplets with intermediate Stokes numbers were centrifuged by the vortices to the outer edges of the wake, yielding an increased dispersion. Large Stokes number droplets were found to be less responsive to the vortex structures and exhibited the least dispersion. Particle concentration was also correlated with vorticity distribution which yielded preferential bias effects as a function of different particle sizes. This trend was qualitatively similar to results seen in homogenous isotropic turbulence, though the influence of particle inertia was less pronounced for the cylinder wake case. A similar study was completed for droplet dispersion within the Icing Research Tunnel (IRT) at the NASA Glenn Research Center, where it is important to obtain a nearly uniform liquid water content (LWC) distribution in the test section (to recreate atmospheric icing conditions).. For this goal, droplets are diffused by the mean and turbulent flow generated from the nozzle air jets, from the upstream spray bars, and from the vertical strut wakes. To understand the influence of these three components, a set of simulations was conducted with a sequential inclusion of these components. Firstly, a jet in an otherwise quiescent airflow was simulated to capture the impact of the air jet on flow turbulence and droplet distribution, and the predictions compared well with experimental results. The effects of the spray bar wake and vertical strut wake were then included with two more simulation conditions, for which it was found that the air jets were the primary driving force for droplet dispersion, i.e. that the spray bar and vertical strut wake effects were secondary.

Desenvolvimento de um método para correção de curvas de desempenho em bombas centrífugas submersas operando com fluidos viscosos

Electrical Submersible Pump (ESP) is used as an artificial lift technique. However, pumping viscous oil is generally associated with low Reynolds number flows. This condition leads to a performance degradation respect to the performance expected from the regular operation with water that most of the centrifugal pumps are originally designed for. These issues are considered in this investigation through a numerical study of the flow in two different multistage, semi-axial type ESPs. This investigation is carried out numerically using a Computational Fluid Dynamics (CFD) package, where the transient RANS equations are solved numerically. The turbulence is modeled using the SST model. Head curves for several operating conditions are compared with manufacturer’s curves and experimental data for a three-stage ESP, showing good agreement for a wide range of fluid viscosities and rotational speeds. Dimensionless numbers (n, n, n e Re) are used to investigate performance degradation of the ESPs. In addition, flow phenomena through the impellers of the ESPs are investigated using flow field from numerical results. Results show that performance degradation is directly related to rotational Reynolds number, Re. In addition, it was verified that performance degradation occurs for constant normalized specific speedn, which shows that performance degradation occurs similarly for different centrifugal pumps. Moreover, experimental data and numerical results agreed with a correlation from literature between head and flow correction factors proposed by Stepanoff (1967). A definition of modified Reynolds number was proposed and relates the head correction factor to viscosity. A correlation between head correction factor and the modified Reynolds number was proposed, which agreed well with numerical and experimental data. Then, a method to predict performance degradation based on the previous correlations was proposed. This method was compared with others from literature. In general, results and conclusions from this work can also be useful to bring more information about the flow of highly viscous fluids in pumps, especially in semi-axial, multistage ESPs.

In situ characterization of sediment oxygen demand across the spectrum of non-resuspension to full resuspension with varying bed shear stress

Sediment oxygen demand (SOD) can be a significant oxygen sink in various types of water bodies, particularly slow-moving waters with substantial organic sediment accumulation. In most settings where SOD is a concern, the prevailing hydraulic conditions are such that the impact of sediment resuspension on SOD is not considered. However, in the case of Bubbly Creek in Chicago, Illinois, the prevailing slack water conditions are interrupted by infrequent intervals of very high flow rates associated with pumped combined sewer overflow (CSO) during intense hydrologic events. These events can cause resuspension of the highly organic, nutrient-rich bottom sediments, resulting in precipitous drawdown of dissolved oxygen (DO) in the water column. While many past studies have addressed the dependence of SOD on near-bed velocity and bed shear stress prior to the point of sediment resuspension, there has been limited research that has attempted to characterize the complex and dynamic phenomenon of resuspended-sediment oxygen demand. To address this issue, a new in situ experimental apparatus referred to as the U of I Hydrodynamic SOD Sampler was designed to achieve a broad range of velocities and associated bed shear stresses. This allowed SOD to be analyzed across the spectrum of no sediment resuspension associated with low velocity/ bed shear stress through full sediment resuspension associated with high velocity / bed shear stress. The current study split SOD into two separate components: (1) SODNR is the sediment oxygen demand associated with non-resuspension conditions and is a surface sink calculated using traditional methods to yield a value with units (g/m2/day); and (2) SODR is the oxygen demand associated with resuspension conditions, which is a volumetric sink most accurately characterized using non-traditional methods and units that reflect suspension in the water column (mg/L/day). In the case of resuspension, the suspended sediment concentration was analyzed as a function of bed shear stress, and a formulation was developed to characterize SODR as a function of suspended sediment concentration in a form similar to first-order biochemical oxygen demand (BOD) kinetics with Monod DO term. The results obtained are intended to be implemented into a numerical model containing hydrodynamic, sediment transport, and water quality components to yield oxygen demand varying in both space and time for specific flow events. Such implementation will allow evaluation of proposed Bubbly Creek water quality improvement alternatives which take into account the impact of SOD under various flow conditions. Although the findings were based on experiments specific to the conditions in Bubbly Creek, the techniques and formulations developed in this study should be applicable to similar sites.

Slot Film Cooling: A Comprehensive Experimental Characterization

When components of a propulsion system are exposed to elevated flow temperatures there is a risk for catastrophic failure if the components are not properly protected from the thermal loads. Among several strategies, slot film cooling is one of the most commonly used, yet poorly understood active cooling techniques. Tangential injection of a relatively cool fluid layer protects the surface(s) in question, but the turbulent mixing between the hot mainstream and cooler film along with the presence of the wall presents an inherently complex problem where kinematics, thermal transport and multimodal heat transfer are coupled. Furthermore, new propulsion designs rely heavily on CFD analysis to verify their viability. These CFD models require validation of their results, and the current literature does not provide a comprehensive data set for film cooling that meets all the demands for proper validation, namely a comprehensive (kinematic, thermal and boundary condition data) data set obtained over a wide range of conditions. This body of work aims at solving the fundamental issue of validation by providing high quality comprehensive film cooling data (kinematics, thermal mixing, heat transfer). 3 distinct velocity ratios (VR=uc/u∞) are examined corresponding to wall-wake (VR~0.5), min-shear (VR ~ 1.0), and wall-jet (VR~2.0) type flows at injection, while the temperature ratio TR= T∞/Tc is approximately 1.5 for all cases. Turbulence intensities at injection are 2-4% for the mainstream (urms/u∞, vrms/u∞,), and on the order of 8-10% for the coolant (urms/uc, vrms/uc,). A special emphasis is placed on inlet characterization, since inlet data in the literature is often incomplete or is of relatively low quality for CFD development. The data reveals that min-shear injection provides the best performance, followed by the wall-jet. The wall-wake case is comparably poor in performance. The comprehensive data suggests that this relative performance is due to the mixing strength of each case, as well as the location of regions of strong mixing with respect to the wall. Kinematic and thermal data show that strong mixing occurs in the wall-jet away from the wall (y/s>1), while strong mixing in the wall-wake occurs much closer to the wall (y/s<1). Min-shear cases exhibit noticeably weaker mixing confined to about y/s=1. Additionally to these general observations, the experimental data obtained in this work is analyzed to reveal scaling laws for the inlets, near-wall scaling, detecting and characterizing coherent structures in the flow as well as to provide data reduction strategies for comparison to CFD models (RANS and LES).

Simulação numérica do escoamento monofásico em uma bomba centrífuga radial

Centrifugal pumps are vastly used in many industrial applications. Knowledge of how these components behave in several circumstances is crucial for the development of more efficient and, therefore, less expensive pumping installations. The combination of multiple impellers, vaned diffusers and a volute might introduce several complex flow characteristics that largely deviate from regular inviscid pump flow theory. Computational Fluid Dynamics can be very helpful to extract information about which physical phenomena are involved in such flows. In this sense, this work performs a numerical study of the flow in a two-stage centrifugal pump (Imbil ITAP 65-330/2) with a vaned diffuser and a volute. The flow in the pump is modeled using the software Ansys CFX, by means of a multi-block, transient rotor-stator technique, with structured grids for all pump parts. The simulations were performed using water and a mixture of water and glycerin as work fluids. Several viscosities were considered, in a range between 87 and 720 cP. Comparisons between experimental data obtained by Amaral (2007) and numerical head curves showed a good agreement, with an average deviation of 6.8% for water. The behavior of velocity, pressure and turbulence kinetic energy fields was evaluated for several operational conditions. In general, the results obtained by this work achieved the proposed goals and are a significant contribution to the understanding of the flow studied.

Simulação numérica da convecção mista em cavidade preenchida com meio poroso heterogêneo e homogêneo

In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.