Hybrid LES-RANS using synthesized turbulent fluctuations for forcing in the interface region

The main bottleneck in using Large Eddy Simulations at high Reynolds number is the requirement of very fine meshes near walls. One of the main reasons why hybrid LES-RANS was invented was to eliminate this limitation. In this method unsteady RANS (URANS) is used near walls and LES is used away from walls. The present paper evaluates a method for improving standard LES-RANS. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane between the LES and URANS regions in order to trigger the equations to resolve turbulence. The turbulent fluctuations are taken from synthesized homogeneous turbulence assuming a modified von Kármán spectrum. Both isotropic and non-isotropic fluctuations are evaluated. The new approach is applied to fully developed channel flow and it is shown that the imposed fluctuations considerably improve the predictions. It is found that increasing the prescribed turbulent length scale of the synthesized turbulence provides excellent agreement with the classical log-law.     

Non-zonal detached eddy simulation coupled with a steady RANS solver in the wall region

Xiao and Jenny (2012) proposed an interesting hybrid LES/RANS method in which they use two solvers and solve the RANS and LES equations in the entire computational domain. In the present work this method is simplified and used as a hybrid RANS-LES method, a wall-modeled LES. The two solvers are employed in the entire domain. Near the walls, the flow is governed by the steady RANS solver; drift terms are added to the DES equations to ensure that the time-averaged DES fields agree with the steady RANS field. Away from the walls, the flow is governed by the DES solver; in this region, the RANS field is set to the time-averaged LES field. The disadvantage of traditional DES models is that the RANS models in the near-wall region – which originally were developed and tuned for steady RANS – are used as URANS models where a large part of the turbulence is resolved. In the present method – where steady RANS is used in the near-wall region – the RANS turbulence models are used in a context for which they were developed. In standard DES methods, the near-wall accuracy can be degraded by the unsteady agitation coming from the LES region. It may in the present method be worth while to use an accurate, advanced RANS model. The EARSM model is used in the steady RANS solver. The new method is called NZ S-DES . It is found to substantially improve the predicting capability of the standard DES. A great advantage of the new model is that it is insensitive to the location of the RANS-LES interface.     

Towards a Unified Turbulence Simulation Approach for Wall-Bounded Flows

A hybrid Reynolds-averaged Navier–Stokes/Large-Eddy Simulation (RANS/LES) methodology has received considerable attention in recent years, especially in its application to wall-bounded flows at high-Reynolds numbers. In the conventional zonal hybrid approach, eddy-viscosity-type RANS and subgrid scale models are applied in the RANS and LES zones, respectively. In contrast, the non-zonal hybrid approach uses only a generalized turbulence model, which provides a unified simulation approach that spans the continuous spectrum of modeling/simulation schemes from RANS to LES. A particular realization of the non-zonal approach, known as partially resolved numerical simulation (PRNS), uses a generalized turbulence model obtained from a rescaling of a conventional RANS model through the introduction of a resolution control function F _R , where F _R is used to characterize the degree of modeling required to represent the unresolved scales of turbulent motion. A new generalized functional form for F _R in PRNS is proposed in this study, and its performance is compared with unsteady RANS (URANS) and LES computations for attached and separated wall-bounded turbulent flows. It is demonstrated that PRNS behaves similarly to LES, but outperforms URANS in general.     

Log-layer mismatch and commutation error in hybrid RANS/LES simulation of channel flow

Hybrid approach combining large eddy simulation (LES) with the Reynolds-averaged Navier–Stokes equation (RANS) is expected to accurately simulate wall-bounded turbulent flows at high Reynolds numbers. As an important issue in developing hybrid methods, it is known that the log layers in the RANS and LES regions are not lined up in hybrid RANS/LES simulations of channel flow. Although several methods including additional filtering near the RANS/LES interface have been proposed to eliminate the log-layer mismatch, there is no obvious physical justification for the methods and some ad hoc tuning is necessary. In this work, the commutation error terms in the filtered velocity equations are investigated to justify the method of additional filtering. It is shown that the additional filtering can be considered as a finite difference approximation to extra terms due to the non-commutivity between the hybrid filter and the spatial derivative. Moreover, an expression determining the filter width and its location for the additional filtering is obtained. To validate the expression, a hybrid simulation of channel flow is carried out. The additional filtering with the filter width derived is shown to be effective in eliminating the log-layer mismatch and improving the mean velocity profile.     

细颗粒泥沙净冲刷和输移的大涡模拟研究

在传统水沙输移数值模拟研究中一般采用雷诺时均模拟技术(Reynolds-averaged simulation,RANS).与RANS相比,大涡模拟技术(large eddy simulation,LES)能够更加精确反映细部流动结构,计算机的发展使得采用LES探讨水流和泥沙运动规律成为可能.本文尝试给出净冲刷条件下悬沙计算的边界条件,采用动态亚格子模式对循环槽道和长槽道中的水流运动和泥沙输移进行了三维大涡模拟研究.利用直接数值模拟(direct numerical simulation,DNS)结果对LES模型进行了率定,计算结果符合良好,在此基础上初步探讨了泥沙浓度、湍动强度和湍动通量等的分布特征.结果表明,净冲刷条件下输沙平衡时泥沙浓度符合Rouse公式分布,单向流动中泥沙浓度沿着流向逐渐增大.泥沙浓度湍动强度和湍动通量都在近底部达到最大值,沿着垂向迅速减小.湍动黏性系数和扩散系数基本上在水深中间处达到最大.湍动Schmidt数沿着水深方向不是常数,在近底部和自由水面附近较大,在水深中间处较小.     

A Rapid and Accurate Switch from RANS to LES in Boundary Layers Using an Overlap Region

An efficient recycling algorithm is developed for injecting resolved turbulent content in a boundary layer as it switches from a Reynolds Averaged Navier-Stokes (RANS) type treatment to a Large Eddy Simulation (LES) type treatment inside a generalized Detached-Eddy Simulation (DES). The motivation is to use RANS in the thinnest boundary-layer area, following the original argument in favour of DES, and LES in the thicker boundary-layer areas especially approaching separation, to improve accuracy and possibly obtain unsteady outputs. The algorithm relies on an overlap of the RANS and LES domains and, therefore, the availability of both RANS and LES solutions in the recycling region, which is about 5 boundary-layer thicknesses long. This permits a smooth transfer of the turbulent stresses from this section to the LES inflow. The continuity of the skin-friction distribution is very good, reflecting the excellent viability of the resolved turbulence. The approach is validated in a flat-plate boundary layer and an airfoil near stall, with mild pressure gradient near the interface, and then applied to the compressible flow over an idealized airliner windshield wiper. The pressure fluctuations at reattachment are 12dB more intense than under a simple boundary layer at the same speed, and the output contains all the quantities needed to calculate the transmission of sound through the glass.     

URANS prediction of flow fluctuations in rod bundle with split-type spacer grid

Turbulent flow in a rod bundle with split-type spacer grid has been studied using Unsteady Reynolds-Averaged Navier–Stokes (URANS) approach. In the previous studies of turbulent flow in rod bundles URANS (as well as steady-state RANS) simulations predicted mean velocity profiles fairly well. However, they severely underpredicted velocity fluctuations, which is investigated in the present study. Our simulations were performed with the Shear Stress Transport (SST) turbulence model and automatic wall-treatment using OpenFOAM, an open-source CFD code. Results of URANS simulations are compared with the measurements of the MATiS-H experiment, which was performed at Korean Atomic Energy Research Institute (KAERI) in 2011–2012.The URANS predictions of velocity fluctuations have been improved by appropriately summing up fluctuations resolved by the basic URANS model and non-resolved fluctuations, which were modelled with the turbulence model. This treatment of turbulent fluctuations, which are directly measured in high-quality experiments, allows more detailed evaluation of various URANS turbulence models. It was found out that the best agreement is achieved when resolved and modelled fluctuations are assumed to be uncorrelated, which indicates that the large-scale structures in this particular flow are distinct in the spectral space from the rest of turbulence. Turbulent flow in the MATiS-H experiment was reproduced by numerous authors using different approaches and our results are among the most accurate.     

Analysis of Head Losses in a Turbine Draft Tube by Means of 3D Unsteady Simulations

In this paper, Unsteady RANS (URANS) simulations and Large Eddy Simulations (LES) in the draft tube of a bulb turbine are presented with the objective to understand and locate the head losses in this turbine component. Three operating points of the turbine are considered. Numerical results are compared with experimental velocity measurements for validation. Thanks to a detailed analysis of the energy balance in the draft tube, the physical and hydrodynamic phenomena responsible for head losses in the draft tube are identified. Head losses are due to transfer of mean kinetic energy to the turbulent flow and viscous dissipation of kinetic energy. This occurs mainly in the central vortex structure and next to the walls in the draft tube. Head losses prediction is found to be dependent on the turbulence model used in the simulations, especially in URANS simulations. Using this analysis, the evolution of head losses between the three operating points is understood.     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.     

Hybrid LES–RANS technique based on a one-equation near-wall model

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES–RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196–205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES–RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES–RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.      

Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

High-order Large Eddy Simulations of Confined Rotor-Stator Flows

In many engineering and industrial applications, the investigation of rotating turbulent flow is of great interest. In rotor-stator cavities, the centrifugal and Coriolis forces have a strong influence on the turbulence by producing a secondary flow in the meridian plane composed of two thin boundary layers along the disks separated by a non-viscous geostrophic core. Most numerical simulations have been performed using RANS and URANS modelling, and very few investigations have been performed using LES. This paper reports on quantitative comparisons of two high-order LES methods to predict a turbulent rotor-stator flow at the rotational Reynolds number Re(=?Ωb ²/ν)?=4 × 10⁵. The classical dynamic Smagorinsky model for the subgrid-scale stress (Germano et al., Phys Fluids A 3(7):1760–1765, 1991) is compared to a spectral vanishing viscosity technique (Séverac & Serre, J Comp Phys 226(2):1234–1255, 2007). Numerical results include both instantaneous data and post-processed statistics. The results show that both LES methods are able to accurately describe the unsteady flow structures and to satisfactorily predict mean velocities as well as Reynolds stress tensor components. A slight advantage is given to the spectral SVV approach in terms of accuracy and CPU cost. The strong improvements obtained in the present results with respect to RANS results confirm that LES is the appropriate level of modelling for flows in which fully turbulent and transition regimes are involved.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

The artificial buffer layer and the effects of forcing in hybrid LES/RANS

The present study is focused on large eddy simulations (LES) that use a statistical (RANS) turbulence model near solid walls, and on the artificial buffer layer that is formed at the interface between these two modeling regions. Additional forcing is used to trigger resolved motions in the LES region more quickly, and leads to improved results in several ways. The study investigates the artificial buffer layer and how it changes with the use of forcing in an in-depth manner, with the purpose of increased understanding of the increasingly popular hybrid LES/RANS group of methods.

The artificial buffer layer is shown to extend from below the modeling interface to well above it, in fact up to 20% of the boundary layer thickness for the cases studied here. The artificial buffer layer is found to be similar to the true buffer layer in many aspects, including a high correlation between the streamwise and wall normal velocity components in the ‘superstreaks’. This indicates that while the superstreaks are highly anisotropic and have unphysical length scales, they still contribute to the resolved shear stress. The forcing does not remove the artificial buffer layer, but it does reduce its extent and increases the resolved shear stress. This increase is mainly associated with increased fluctuations of the wall normal velocity.

A simple, low-dimensional forcing model is proposed and tested, with favorable results. The model is simple to implement and easily generalized to more complex geometries.     

A Zonal RANS/LES Approach for Noise Sources Prediction

Hybrid CFD/CAA methods have generally to be used for the numerical simulation of trailing-edge noise (see [9, 20] for instance). This study focuses on the first step of such hybrid methods, which is to predict the unsteady aerodynamic sources by the mean of a 3D unsteady simulation of the flow. Such a simulation is however generally still away from the numerical capabilities of ‘usual’ supercomputers. This paper investigates the use of a zonal LES method (based on the NLDE – Non-Linear Disturbance Equations – technique) for the numerical prediction of the aerodynamic noise sources. This method makes it possible to perform only zonal LES close to the main elements responsible of sound generation, while the overall configuration is only treated by a RANS approach. Attention will be paid to the specific boundary treatment at the interface between the RANS and LES regions. More precisely, the problem of the generation of turbulent inflow conditions for the LES region will be carefully addressed. The method is first assessed in the simulation of a flat plate ended by a blunted trailing-edge, and then applied to the simulation of the flow over a NACA0012 airfoil with blunted trailing-edge.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

Hybrid URANS/LES Turbulence Simulation of Vortex Shedding Behind a Triangular Flameholder

In this study a detached eddy simulation (DES) model, which belongs to the group of hybrid URANS/LES turbulence models, is used for the simulation of vortex shedding behind a triangular obstacle. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS shear-stress transport (SST) model is used. In the other regions with higher grid resolution a LES model, which uses a transport equation for the turbulent subgrid energy, is applied. The DES model is first investigated for two standard test cases, namely decaying homogeneous isotropic turbulence and the backward facing step, respectively. For the decaying homogeneous isotropic turbulence test case the evolution of the energy spectra in wavenumber space for different times are studied for both the DES and a Smagorinsky type LES model. Different grid resolutions are analyzed with a special emphasis on the modeling constant connecting the filter length scale to the grid size. The results are compared to experimental data. The backward facing step test case is used to study the model behavior for a case with a transition region between a RANS modeling approach close to the wall and LES based modeling in the intense shear flow region. The final application is the simulation of the vortex shedding behind a triangular obstacle. First, the influence of the inlet condition formulation is studied in detail as they can have a significant influence especially for LES based models. Detailed comparisons between simulation and experiment for the flow structure past the obstacle and statistical quantities such as the shedding frequency are shown. Finally the additional temporal and spatial information provided by the DES model is used to show the predicted anisotropy of turbulence.     

Large Eddy and Reynolds-Averaged Navier - Stokes Simulations of Turbulent Flow Over an Airfoil

A Large Eddy Simulation (LES) of turbulent flow over an airfoil near stall is performed. Results of the LES are compared with those of Reynolds-Averaged Navier-Stokes (RANS) simulations using two well-known turbulence models, namely the Baldwin-Lomax model and the Spalart-Allmaras model. The subgrid scale model used for the LES is the filtered structure function model. All simulations are performed using the same structured multi-block code. In order to reduce the CPU time, an implicit time stepping method is used for the LES. The purpose of this study is to show the possibilities and limitations of LES of complex flows associated with aeronautical applications using state of the art simulation techniques. Typical flow features are captured by the LES such as the adverse-pressure gradient and flow retardation. Visualization of instantaneous flow fields shows the typical streaky structures in the near-wall region.     

Highly Resolved LES and URANS of Turbulent Buoyancy-Driven Flow Within Inclined Differentially-Heated Enclosures

In the present study, both LES and unsteady RANS computations are presented, for turbulent natural convection of air inside differentially-heated rectangular tilted cavities using a finite volume code (Code_Saturne), for an aspect ratio of H/L?=?28.6 and Rayleigh number of 0.86×10⁶. Attention is focused on two angles of inclination: 15° to the horizontal with hot lower and cold upper wall, the 15° unstable case, and the mirror image of this case where the angle is the same but with a hot upper and cold lower wall, the 15° stable case. In accordance with recent experimental data, the LES computations for both the stable and unstable tilted cavities returned three-dimensional time-averaged flow fields. In the case of the unstably stratified enclosure, the flow is highly unsteady with coherent turbulent structures in the core of the enclosure. Time-averaged temperature, velocity and resolved turbulence intensities resulting from LES computations show close agreement to measured data. Subsequent comparisons of different URANS schemes with LES are used in order to explore to what extent these models are able to reproduce the large-scale unsteady flow structures. All URANS schemes have been found to be able to reproduce the 3-D unsteady flow features present in the 15° unstable cavity. However, the low-Reynolds-number model tested, as well as requiring a high resolution near-wall grid, also needed a finer grid in the core region than the high-Reynolds-number models, thus making it computationally very expensive. Flow within the 15° stable cavity also shows some 3-D features, although it is significantly less unsteady, and the URANS models tested here have been less successful in reproducing this flow pattern. The overall heat transfer is presented here for both differentially heated enclosures.