Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection

G.L. Kooij, M.A. Botchev, E.M.A. Frederix, B.J. Geurts, S. Horn, D. Lohse, E.P. van der Poel, O. Shishkina, R.J.A.M. Stevens, R. Verzicco

Research output: Contribution to journalArticleAcademicpeer-review

10 Citations (Scopus)

Abstract

Computational codes for direct numerical simulations of Rayleigh–Bénard (RB) convection are compared in terms of computational cost and quality of the solution. As a benchmark case, RB convection at Ra=108 and Pr=1 in a periodic domain, in cubic and cylindrical containers is considered. A dedicated second-order finite-difference code (AFID/RBFLOW) and a specialized fourth-order finite-volume code (GOLDFISH) are compared with a general purpose finite-volume approach (OPENFOAM) and a general purpose spectral-element code (NEK5000). Reassuringly, all codes provide predictions of the average heat transfer that converge to the same values. The computational costs, however, are found to differ considerably. The specialized codes AFID/RBFLOW and GOLDFISH are found to excel in efficiency, outperforming the general purpose flow solvers NEK5000 and OPENFOAM by an order of magnitude with an error on the Nusselt number Nu below 5%. However, we find that Nu alone is not sufficient to assess the quality of the numerical results: in fact, instantaneous snapshots of the temperature field from a near wall region obtained for deliberately under-resolved simulations using NEK5000 clearly indicate inadequate flow resolution even when Nu is converged. Overall, dedicated special purpose codes for RB convection are found to be more efficient than general purpose codes.

Original languageEnglish
Pages (from-to)1-8
Number of pages8
JournalComputers & Fluids
Volume166
DOIs
Publication statusPublished - 30 Apr 2018

Fingerprint

Direct numerical simulation
Nusselt number
Containers
Costs
Temperature distribution
Heat transfer
Convection

Keywords

  • Direct numerical simulations
  • Heat transfer
  • Rayleigh–Bénard convection

Cite this

Kooij, G.L. ; Botchev, M.A. ; Frederix, E.M.A. ; Geurts, B.J. ; Horn, S. ; Lohse, D. ; van der Poel, E.P. ; Shishkina, O. ; Stevens, R.J.A.M. ; Verzicco, R. / Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection. In: Computers & Fluids. 2018 ; Vol. 166. pp. 1-8.
@article{d741cf8ff4674fd88bf384f3e6b7a011,
title = "Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–B{\'e}nard convection",
abstract = "Computational codes for direct numerical simulations of Rayleigh–B{\'e}nard (RB) convection are compared in terms of computational cost and quality of the solution. As a benchmark case, RB convection at Ra=108 and Pr=1 in a periodic domain, in cubic and cylindrical containers is considered. A dedicated second-order finite-difference code (AFID/RBFLOW) and a specialized fourth-order finite-volume code (GOLDFISH) are compared with a general purpose finite-volume approach (OPENFOAM) and a general purpose spectral-element code (NEK5000). Reassuringly, all codes provide predictions of the average heat transfer that converge to the same values. The computational costs, however, are found to differ considerably. The specialized codes AFID/RBFLOW and GOLDFISH are found to excel in efficiency, outperforming the general purpose flow solvers NEK5000 and OPENFOAM by an order of magnitude with an error on the Nusselt number Nu below 5{\%}. However, we find that Nu alone is not sufficient to assess the quality of the numerical results: in fact, instantaneous snapshots of the temperature field from a near wall region obtained for deliberately under-resolved simulations using NEK5000 clearly indicate inadequate flow resolution even when Nu is converged. Overall, dedicated special purpose codes for RB convection are found to be more efficient than general purpose codes.",
keywords = "Direct numerical simulations, Heat transfer, Rayleigh–B{\'e}nard convection",
author = "G.L. Kooij and M.A. Botchev and E.M.A. Frederix and B.J. Geurts and S. Horn and D. Lohse and {van der Poel}, E.P. and O. Shishkina and R.J.A.M. Stevens and R. Verzicco",
year = "2018",
month = "4",
day = "30",
doi = "10.1016/j.compfluid.2018.01.010",
language = "English",
volume = "166",
pages = "1--8",
journal = "Computers & Fluids",
issn = "0045-7930",
publisher = "Elsevier",

}

Kooij, GL, Botchev, MA, Frederix, EMA, Geurts, BJ, Horn, S, Lohse, D, van der Poel, EP, Shishkina, O, Stevens, RJAM & Verzicco, R 2018, 'Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection', Computers & Fluids, vol. 166, pp. 1-8. https://doi.org/10.1016/j.compfluid.2018.01.010

Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection. / Kooij, G.L.; Botchev, M.A.; Frederix, E.M.A.; Geurts, B.J.; Horn, S.; Lohse, D.; van der Poel, E.P.; Shishkina, O.; Stevens, R.J.A.M.; Verzicco, R.

In: Computers & Fluids, Vol. 166, 30.04.2018, p. 1-8.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection

AU - Kooij, G.L.

AU - Botchev, M.A.

AU - Frederix, E.M.A.

AU - Geurts, B.J.

AU - Horn, S.

AU - Lohse, D.

AU - van der Poel, E.P.

AU - Shishkina, O.

AU - Stevens, R.J.A.M.

AU - Verzicco, R.

PY - 2018/4/30

Y1 - 2018/4/30

N2 - Computational codes for direct numerical simulations of Rayleigh–Bénard (RB) convection are compared in terms of computational cost and quality of the solution. As a benchmark case, RB convection at Ra=108 and Pr=1 in a periodic domain, in cubic and cylindrical containers is considered. A dedicated second-order finite-difference code (AFID/RBFLOW) and a specialized fourth-order finite-volume code (GOLDFISH) are compared with a general purpose finite-volume approach (OPENFOAM) and a general purpose spectral-element code (NEK5000). Reassuringly, all codes provide predictions of the average heat transfer that converge to the same values. The computational costs, however, are found to differ considerably. The specialized codes AFID/RBFLOW and GOLDFISH are found to excel in efficiency, outperforming the general purpose flow solvers NEK5000 and OPENFOAM by an order of magnitude with an error on the Nusselt number Nu below 5%. However, we find that Nu alone is not sufficient to assess the quality of the numerical results: in fact, instantaneous snapshots of the temperature field from a near wall region obtained for deliberately under-resolved simulations using NEK5000 clearly indicate inadequate flow resolution even when Nu is converged. Overall, dedicated special purpose codes for RB convection are found to be more efficient than general purpose codes.

AB - Computational codes for direct numerical simulations of Rayleigh–Bénard (RB) convection are compared in terms of computational cost and quality of the solution. As a benchmark case, RB convection at Ra=108 and Pr=1 in a periodic domain, in cubic and cylindrical containers is considered. A dedicated second-order finite-difference code (AFID/RBFLOW) and a specialized fourth-order finite-volume code (GOLDFISH) are compared with a general purpose finite-volume approach (OPENFOAM) and a general purpose spectral-element code (NEK5000). Reassuringly, all codes provide predictions of the average heat transfer that converge to the same values. The computational costs, however, are found to differ considerably. The specialized codes AFID/RBFLOW and GOLDFISH are found to excel in efficiency, outperforming the general purpose flow solvers NEK5000 and OPENFOAM by an order of magnitude with an error on the Nusselt number Nu below 5%. However, we find that Nu alone is not sufficient to assess the quality of the numerical results: in fact, instantaneous snapshots of the temperature field from a near wall region obtained for deliberately under-resolved simulations using NEK5000 clearly indicate inadequate flow resolution even when Nu is converged. Overall, dedicated special purpose codes for RB convection are found to be more efficient than general purpose codes.

KW - Direct numerical simulations

KW - Heat transfer

KW - Rayleigh–Bénard convection

UR - http://www.scopus.com/inward/record.url?scp=85041475622&partnerID=8YFLogxK

U2 - 10.1016/j.compfluid.2018.01.010

DO - 10.1016/j.compfluid.2018.01.010

M3 - Article

AN - SCOPUS:85041475622

VL - 166

SP - 1

EP - 8

JO - Computers & Fluids

JF - Computers & Fluids

SN - 0045-7930

ER -