LES, T-RANS and hybrid simulations of thermal convection at high Ra numbers

Abstract The paper reports on application of different approaches to the simulations of thermal convection at high Rayleigh ( Ra ) numbers. Based on new well-resolved LES in 10 7  ⩽  Ra  ⩽ 10 9 range, the performance of a T-RANS (using a low- Re three-equation 〈 k 〉–〈 e 〉–〈 θ 2 〉 ASM/AFM subscale model) and a hybrid approach (utilizing the concept of “seamless” RANS/LES merging) have been compared. Targeting accurate predictions of heat transfer at very high Ra numbers, the near-wall behaviour of the subscale turbulence contributions is analyzed in details. Whilst the application of the conventional LES on a coarse grid resulted in huge underprediction of the Nusselt number (50% at Ra  = 10 9 ), thanks to a well-tuned subscale model, the T-RANS results showed excellent agreement for heat transfer with both the fine-resolved LES (for Ra  ⩽ 10 9 ) and experimental data for over a ten-decades range of Ra (10 6  ⩽  Ra  ⩽ 2 × 10 16 ). The visible absence of the fine-scale motion in the T-RANS, however, means that the T-RANS will perform poorly in flows where there are no strong large-scale forcing, as, e.g., in a side-heated vertical channel. In order to sensitize the T-RANS approach to high-frequency instabilities, different ways of hybrid RANS/LES merging based on “seamless” approach have been investigated. It is demonstrated that the hybrid approach is capable of capturing a significantly larger portion of the fine-structure spectrum than is possible with T-RANS, whilst also returning accurate predictions of heat transfer and turbulence statistics.

[1]  Roland Schiestel,et al.  Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations , 2005 .

[2]  D. C. Threlfall,et al.  Free convection in low-temperature gaseous helium , 1975, Journal of Fluid Mechanics.

[3]  Sasa Kenjeres,et al.  Transient analysis of Rayleigh–Bénard convection with a RANS model , 1999 .

[4]  Michael A. Leschziner,et al.  Merging LES and RANS Strategies: Zonal or Seamless Coupling? , 2004 .

[5]  Sasa Kenjeres,et al.  Statistical Structure at the Wall of the High Reynolds Number Turbulent Boundary Layer , 2001 .

[6]  Robert McDougall Kerr,et al.  Rayleigh number scaling in numerical convection , 1996, Journal of Fluid Mechanics.

[7]  M. Markus,et al.  On-off intermittency and intermingledlike basins in a granular medium. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  C. G. Speziale Turbulence modeling for time-dependent RANS and VLES : a review , 1998 .

[9]  B. Geurts,et al.  Direct and large-eddy simulation IV , 2001 .

[10]  Bernardus J. Geurts Modern Simulation Strategies for Turbulent Flow , 2001 .

[11]  J. Ullrich,et al.  Rescattering of ultralow-energy electrons for single ionization of Ne in the tunneling regime. , 2003, Physical review letters.

[12]  Joseph Smagorjnsky,et al.  The Beginnings of Numerical Weather Prediction and General Circulation Modeling: Early Recollections , 1983 .

[13]  Libchaber,et al.  Scaling relations in thermal turbulence: The aspect-ratio dependence. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[14]  K. Hanjalić Synergy of experiments and computer simulations in research of turbulent convection , 2005 .

[15]  S. Zaleski,et al.  Scaling of hard thermal turbulence in Rayleigh-Bénard convection , 1989, Journal of Fluid Mechanics.

[16]  Douglas K. Lilly,et al.  A comparison of two dynamic subgrid closure methods for turbulent thermal convection , 1994 .

[17]  Jackson R. Herring,et al.  Prandtl number dependence of Nusselt number in direct numerical simulations , 2000, Journal of Fluid Mechanics.

[18]  B. E. Launder,et al.  Closure Strategies for Turbulent and Transitional Flows , 2002 .

[19]  J. Domaradzki,et al.  Large eddy simulations of Rayleigh–Bénard convection using subgrid scale estimation model , 2000 .

[20]  G. Grötzbach,et al.  Spatial resolution requirements for direct numerical simulation of the Rayleigh-Bénard convection , 1983 .

[21]  F. Busse,et al.  Large scale structures in Rayleigh-Bénard convection at high Rayleigh numbers. , 2003, Physical review letters.

[22]  Thomas M. Eidson,et al.  Numerical simulation of the turbulent Rayleigh–Bénard problem using subgrid modelling , 1985, Journal of Fluid Mechanics.

[23]  K. Hanjali,et al.  Numerical insight into flow structure in ultraturbulent thermal convection. , 2002 .

[24]  K. R. Sreenivasan,et al.  Turbulent convection at very high Rayleigh numbers , 1999, Nature.

[25]  Richard J Goldstein,et al.  High-Rayleigh-number convection of pressurized gases in a horizontal enclosure , 2002, Journal of Fluid Mechanics.

[26]  M. Wörner Direkte Simulation turbulenter Rayleigh- Benard-Konvektion in flüssigem Natrium , 1994 .

[27]  F. Chillà,et al.  Turbulent Rayleigh–Bénard convection in gaseous and liquid He , 2001 .

[28]  R. Kraichnan Turbulent Thermal Convection at Arbitrary Prandtl Number , 1962 .

[29]  Kenjeres,et al.  Convective rolls and heat transfer in finite-length rayleigh-Benard convection: A two-dimensional numerical study , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[30]  Simulation and Identification of Deterministic Structures in Thermal and Magnetic Convection , 2002, Annals of the New York Academy of Sciences.