Modeling of Turbulent Natural Convection in Enclosed Tall Cavities

It was shown in our previous work (J. Appl. Mech. Tech. Phys 57 (7), 1159–1171 (2016)) that the eddy-resolving parameter-free CABARET scheme as applied to two-and three-dimensional de Vahl Davis benchmark tests (thermal convection in a square cavity) yields numerical results on coarse (20 × 20 and 20 × 20 × 20) grids that agree surprisingly well with experimental data and highly accurate computations for Rayleigh numbers of up to 10¹⁴. In the present paper, the sensitivity of this phenomenon to the cavity shape (varying from cubical to highly elongated) is analyzed. Box-shaped computational domains with aspect ratios of 1: 4, 1: 10, and 1: 28.6 are considered. The results produced by the CABARET scheme are compared with experimental data (aspect ratio of 1: 28.6), DNS results (aspect ratio of 1: 4), and an empirical formula (aspect ratio of 1: 10). In all the cases, the CABARET-based integral parameters of the cavity flow agree well with the other authors’ results. Notably coarse grids with mesh refinement toward the walls are used in the CABARET calculations. It is shown that acceptable numerical accuracy on extremely coarse grids is achieved for an aspect ratio of up to 1: 10. For higher aspect ratios, the number of grid cells required for achieving prescribed accuracy grows significantly.     

Natural convection flow in a square cavity revisited: Laminar and turbulent models with wall functions

Numerical simulations have been undertaken for the benchmark problem of natural convection flow in a square cavity. The control volume method is used to solve the conservation equations for laminar and turbulent flows for a series of Rayleigh numbers (Ra) reaching values up to 10¹⁰. The k-? model has been used for turbulence modelling with and without logarithmic wall functions. Uniform and non-uniform (stretched) grids have been employed with increasing density to guarantee accurate solutions, especially near the walls for high Ra-values. ADI and SIP solvers are implemented to accelerate convergence. Excellent agreement is obtained with previous numerical solutions, while some discrepancies with others for high Ra-values may be due to a possibly different implementation of the wall functions. Comparisons with experimental data for heat transfer (Nusselt number) clearly demonstrates the limitations of the standard k-? model with logarithmic wall functions, which gives significant overpredictions.     

Numerical aspects of calculation of confined swirling flows with internal recirculation

Predictions were performed for two different confined swirling flows with internal recirculation zones. The convection terms in the elliptic governing equations were discretized using three different finite differencing schemes: hybrid, quadratic upwind interpolation and skew upwind differencing. For each flow case, calculations were carried out with these schemes and successively refined grids were employed. For the turbulent flow case the k-ε turbulence model was used. The predicted cases were a laminar swirling flow investigated by Bornstein and Escudier, and a turbulent low-swirl case studied by Roback and Johnson. In both cases an internal recirculation zone was present. The laminar case is well predicted when account is taken of the estimated radial velocity component at the chosen inlet plane. The quadratic upwind interpolation and skew upwind schemes predict the main features of the internal recirculation zone also with a coarse grid. The turbulent case is well predicted with the coarse as well as the finer grids, the skew upwind and quadratic upwind interpolation schemes yielding results very close to the measurements. It is concluded that the skew upwind scheme reaches grid independence slightly before the quadratic upwind scheme, both considerably earlier than the hybrid scheme.     

Finite volume multigrid prediction of laminar natural convection: Bench-mark solutions

A finite volume multigrid procedure for the prediction of laminar natural convection flows is presented, enabling efficient and accurate calculations on very fine grids. The method is fully conservative and uses second-order central differencing for convection and diffusion fluxes. The calculations start on a coarse (typically 10 × 10 control volumes) grid and proceed to finer grids until the desired accuracy or maximum affordable storage is reached. The computing times increase thereby linearly with the number of control volumes. Solutions are presented for the flow in a closed cavity with side walls at different temperatures and insulated top and bottom walls. Rayleigh numbers of 10⁴, 10⁵ and 10⁶ are considered. Grids as fine as 640 × 640 control volumes are used and the results are believed to be accurate to within 0–01%. Second-order monotonic convergence to grid-independent values is observed for all predicted quantities.     

Thermal and dynamic optimization of the convective flow in Trombe Wall shaped channels by numerical investigation

Two-dimensional, laminar, transitional and turbulent simulations were obtained by solving the fully-elliptic governing equations of the motion established by natural convection in channels, with Trombe Wall configuration, for different geometrical parameters and symmetrical heating. In transitional and turbulent cases, the low-Re k−ω turbulence model has been employed. To validate the numerical results, some comparisons with experimental results taken from literature have been carried out. Numerical results for the average Nusselt number and the non-dimensional induced mass-flow rate have been obtained for a wide and not yet covered range of the Rayleigh number varying from 10⁵ to 10¹². Correlations for the thermal and the mass-flow optimum wall-to-wall spacing have been presented. Finally, additional configurations including discrete heat sources have been studied, in order to obtain thermal and dynamic improvements. These intermediate devices were tested as turbulence generators, in the transitional range of Rayleigh numbers.     

Spectral entropy as a flow state indicator

In many practical applications, the flow state (laminar, transitional, turbulent) might vary in space and/or in time for a given configuration. The aim of the current study is to show that the spectral entropy S_d, obtained from solving the eigenvalue problem for the temporal autocorrelation function, can be used in order to uniquely quantify the flow state and differentiate between laminar, transitional, or turbulent regimes; as such, it delivers a direct measure of turbulence level. Therefore, this quantity might support hybrid numerical simulations by determining the local flow state, identifying in this way the most suitable computational model and switching, e.g., from RANS to LES. The first test of the suggested approach relies on Direct Numerical Simulations (DNS) for decaying Homogeneous Isotropic Turbulence (HIT) performed for ten different Taylor Reynolds numbers. Results obtained by analyzing DNS indicate that S_d is an excellent candidate to quantify turbulence level and transition. To check the robustness of the corresponding analysis, the impact of different resolutions has been investigated, revealing that a correct state estimate is still obtained with a coarser spatial or temporal resolution. Finally, to check the generality of the approach, the entropy thresholds obtained from the DNS analysis have been used with the same algorithm to analyze 1) DNS results obtained for the Taylor-Green vortex benchmark at Re=1600 as well as 2) results obtained through Large Eddy Simulations in a blood nozzle, revealing in both cases a perfect agreement with a traditional, user-based analysis of the flow conditions. Hence, S_d appears to be an excellent quantitative indicator of laminar, transitional, or turbulent flow, allowing an automatic, user-independent analysis of the flow state for a variety of conditions. In principle, it could be used without modification to analyze experimental measurements as well.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Natural convection heat transfer in a uniformly heated horizontal pipe

Natural convection heat transfers inside horizontal pipes were measured. The Rayleigh numbers were varied from 6.8 × 10⁸ to 1.5 × 10¹², while the Prandtl number was fixed at 2,094. Based on the analogy concept, a copper sulfate electroplating system was adopted to measure mass transfer rates in place of heat transfer rates. Test results using single-piece electrodes were in good agreement with the work of Sarac and Korkut. The angle-dependent mass transfer rates, measured using piecewise electrodes, were compared with the results of studies on natural convection in concentric annuli, and showed similar trends. The experiments were expanded to the turbulent region, and a transition criterion was proposed. Angle-dependent natural convection heat transfer correlations for the laminar and turbulent regions were derived.     

Experiments on turbulent natural convection in an enclosed tall cavity

Experiments have been undertaken to investigate the natural convection of air in a tall differentially heated rectangular cavity (2.18 m high by 0.076 m wide by 0.52 m in depth). They were performed with temperature differentials between the vertical plates of 19.6°C and 39.9°C, giving Rayleigh numbers based on the width of 0.86×10⁶ and 1.43×10⁶. Under these conditions the flow in the core of the cavity is fully turbulent and property variations with temperature are comparatively small. A previously used experimental rig has been modified, by fitting partially conducting top and bottom walls and outer guard channels, to provide boundary conditions which avoid the inadequately defined sharp changes in temperature gradient and other problems associated with insufficient insulation on nominally adiabatic walls. Mean and turbulent temperature and velocity variations within the cavity have been measured, together with heat fluxes and turbulent shear stresses. The temperature and flow fields were found to be closely two-dimensional, except close to the front and back walls, and anti-symmetric across the diagonal of the cavity. The partially conducting roof and floor provide locally unstable thermal stratification in the wall jet flows there, which enhances the turbulence as the flow moves towards the temperature controlled plates. The results provide a greatly improved benchmark for the testing of turbulence models in this low turbulence Reynolds number flow.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Direct numerical simulation of turbulent heat transfer in a wall-normal rotating channel flow

Investigations into the characteristics of turbulent heat transfer and coherent flow structures in a plane-channel subjected to wall-normal system rotation are conducted using direct numerical simulation (DNS). In order to investigate the influence of system rotation on the temperature field, a wide range of rotation numbers are tested, with the flow pattern transitioning from being fully turbulent to being quasilaminar, and eventually, fully laminar. In response to the Coriolis force, secondary flows appear as large vortical structures, which interact intensely with the wall shear layers and have a significant impact on the distribution of turbulence kinetic energy (TKE), turbulence scalar energy (TSE), temperature statistics, and turbulent heat fluxes. The characteristic length scales of turbulence structures responsible for the transport of TSE are the largest at the quasilaminar state, which demands a very large computational domain in order to capture the two-dimensional spectra of temperature fluctuations. The effects of the Coriolis force on the turbulent transport processes of the temperature variance and turbulent heat fluxes are thoroughly examined in terms of their respective budget balances.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

Natural convection heat transfer along vertical rectangular ducts

Experimental investigations have been reported on steady state natural convection from the outer surface of vertical rectangular and square ducts in air. Seven ducts have been used; three of them have a rectangular cross section and the rest have square cross section. The ducts are heated using internal constant heat flux heating elements. The temperatures along the vertical surface and the peripheral directions of the duct wall are measured. Axial (perimeter averaged) heat transfer coefficients along the side of each duct are obtained for laminar and transition to turbulent regimes of natural convection heat transfer. Axial (perimeter averaged) Nusselt numbers are evaluated and correlated using the modified Rayleigh numbers for laminar and transition regime using the vertical axial distance as a characteristic length. Critical values of the modified Rayleigh numbers are obtained for transition to turbulent. Furthermore, total overall averaged Nusselt numbers are correlated with the modified Rayleigh numbers and the area ratio for the laminar regimes. The local axial (perimeter averaged) heat transfer coefficients are observed to decrease in the laminar region and increase in the transition region. Laminar regimes are obtained at the lower half of the ducts and its chance to appear decreases as the heat flux increases.     

Interaction of turbulent natural convection and surface thermal radiation in inclined square enclosures

Two-dimensional numerical studies of flow and temperature fields for turbulent natural convection and surface radiation in inclined differentially heated enclosures are performed. Investigations are carried out over a wide range of Rayleigh numbers from 10⁸ to 10¹², with the angle of inclination varying between 0° and 90°. Turbulence is modeled with a novel variant of the k–ε closure model. The predicted results are validated against experimental and numerical results reported in literature. The effect of the inclination of the enclosure on pure turbulent natural convection and the latter’s interaction with surface radiation are brought out. Profiles of turbulent kinetic energy and effective viscosity are studied to observe the net effect on the intensity of turbulence caused by the interaction of natural convection and surface radiation. The variations of local Nusselt number and average Nusselt number are presented for various inclination angles. Marked change in the convective Nusselt number is found with the orientation of enclosure. Also analyzed is the influence of change in emissivity on the flow and heat transfer. A correlation relevant to practical applications in the form of average Nusselt number, as a function of Rayleigh number, Ra, radiation convection parameter, N _RC and inclination angle of the enclosure, φ is proposed.     

DNS of Fractal-Generated Turbulence

An innovative approach which combines high order compact schemes, Immersed Boundary Method and an efficient domain decomposition method is used to perform high fidelity Direct Numerical Simulations (DNS) of four spatially evolving turbulent flows, one generated by a regular grid and three generated by fractal square grids. The main results which we have been able to obtain from these simulations are the following: the vorticity field appears more clustered when generated by fractal square grids compared to a regular grid; fractal square grids generate higher vorticities and turbulence intensities than a regular grid; the flow holds clear geometrical imprints of the fractal grids far downstream, a property which could be used in the future for flow design, management and passive control; the DNS obtained with fractal grids confirmed the existence of two turbulent regions, one where the turbulence progressively amplifies closer to the grid (the production region) followed by one where the turbulence decays; the energy spectra of fluctuating turbulent velocities at various locations in the production region of the flow provide some information on how the turbulence is generated at the smallest scales first near the grid where the smallest wakes are dominant, followed by progressively smaller turbulent frequencies further downstream where progressively larger wakes interact.     

Hysteresis phenomenon in turbulent convection

Coherent large-scale circulations of turbulent thermal convection in air have been studied experimentally in a rectangular box heated from below and cooled from above using Particle Image Velocimetry. The hysteresis phenomenon in turbulent convection was found by varying the temperature difference between the bottom and the top walls of the chamber (the Rayleigh number was changed within the range of 10⁷–10⁸). The hysteresis loop comprises the one-cell and two-cells flow patterns while the aspect ratio is kept constant (A=2–2.23). We found that the change of the sign of the degree of the anisotropy of turbulence was accompanied by the change of the flow pattern. The developed theory of coherent structures in turbulent convection (Phys Rev E 66:1–15, 2002, Boundary-Layer Meteorol, 2005) is in agreement with the experimental observations. The observed coherent structures are superimposed on a small-scale turbulent convection. The redistribution of the turbulent heat flux plays a crucial role in the formation of coherent large-scale circulations in turbulent convection.     

Direct Numerical Simulation of Subsonic Round Turbulent Jet

Direct numerical simulation(DNS) of spatially developing round turbulent jet flow with Reynolds number 4,700 was carried out. Over 20 million grid points were used in this simulation. Fully compressible three-dimensional Navier–Stokes equations were solved. High order explicit spatial difference schemes and Runge–Kutta time integration scheme were used to calculate derivatives and time marching, respectively. Non-reflecting boundary conditions and exit zone techniques were adopted. Some refined computational grids were used in order to capture the smallest turbulent structures near the centerline of the jet. Low level disturbance were imposed on the jet inflow velocity to trigger the developing of turbulence. Turbulent statistics such as mean velocity, Reynolds stresses, third order velocity moments were obtained and compared with experimental data. One-dimensional velocity autospectra was also calculated. The inertial region where the spectra decays according to the k^− 5/3 was observed. The quantitative profiles of mean velocity and all of the third order velocity moments which were difficult to measure via experimental techniques were presented here in detail. The jet flow was proven to be close to fully self-similar around 19 jet diameters downstream of jet exit. The statistic data and revealed flow feature obtained in this paper can provide valuable reference for round turbulent jet research.     

Convective heat transfer and temperature stratification in a sphere completely filled with a liquid,with a given heat flux

A number of articles have been devoted to the theoretical and experimental investigation of natural convection in spherical vessels completely filled with a liquid [1–6]. Analytical solutions are known, obtained by the expansion of the sought function in series in powers of the Rayleigh number (see, for example, [1]), valid for very small values of this number. A numerical solution of the nonlinear Boussinesq equations can be used to obtain solutions with larger Rayleigh numbers, but the existing data for spherical regions [2, 3] embrace a relatively narrow range of Rayleigh numbers. The experimental data with a given heat flux, published in [4–6], were obtained with relatively large Rayleigh numbers (Ra*=10⁹?10¹¹) and Prandtl numbers (P= 3?1500). Data on the characteristics of convection in spherical vessels are still not very numerous and, in a number of cases, contradictory. This relates, in particular, to the boundaries of unsteady-state conditions. The present article, continuing [7–9], expounds a method and gives the results of a calculation of convection in a sphere with a thinwalled shell, in a range of Rayleigh and Fourier numbers embracing the principal conditions of unsteady-state laminar convection with a given heat flux.     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.