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, Detlef; van der Poel, E. P.; Shishkina, Olga; Stevens, R. J. A. M.; Verzicco, R.

doi:10.1016/j.compfluid.2018.01.010

Item

ITEM ACTIONSEXPORT

Add to Basket

Please note that a newer version of this item is available:
https://pure.mpg.de/pubman/item/item_2582893_5

DetailsSummary

Released

Journal Article

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

MPS-Authors

/persons/resource/persons192998

Lohse, Detlef
Max Planck Institute for Dynamics and Self-Organization, Max Planck Society;

/persons/resource/persons180814

Shishkina, Olga
Laboratory for Fluid Dynamics, Pattern Formation and Biocomplexity, Max Planck Institute for Dynamics and Self-Organization, Max Planck Society;

External Resource

No external resources are shared

Fulltext (restricted access)

There are currently no full texts shared for your IP range.

Fulltext (public)

There are no public fulltexts stored in PuRe

Supplementary Material (public)

There is no public supplementary material available

Citation

Kooij, G. L., Botchev, M. A., Frederix, E. M. A., Geurts, B. J., Horn, S., Lohse, D., et al. (2018). Comparison of computational codes for direct numerical simulations of turbulent Rayleigh–Bénard convection. Computers and Fluids, 166, 1-8. doi:10.1016/j.compfluid.2018.01.010.

Cite as: https://hdl.handle.net/21.11116/0000-0001-3621-D

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 and 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.