Parameter-free symmetry-preserving regularization modeling of a turbulent differentially heated cavity Trias, F.X.; Verstappen, R.W.C.P.; Gorobets, A.; Soria, M.; Oliva, A. Link to publication in University of Groningen/UMCG research database Citation for published version (APA): Trias, F. X., Verstappen, R. W. C. P., . Parameter-free symmetrypreserving regularization modeling of a turbulent differentially heated cavity. Computers & fluids, 39(10), 1815-1831. DOI: 10.1016/j.compfluid.2010 Copyright Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons).Take-down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal. For technical reasons the number of authors shown on this cover page is limited to 10 maximum.
a b s t r a c tSince direct numerical simulations of buoyancy driven flows cannot be computed at high Rayleigh numbers, a dynamically less complex mathematical formulation is sought. In the quest for such a formulation, we consider regularizations (smooth approximations) of the non-linearity: the convective term is altered to reduce the production of small scales of motion by means of vortex stretching. In doing so, we propose to preserve the symmetry and conservation properties of the convective terms exactly. This requirement yielded a novel class of regularizations [Comput Fluids 2008;37:887] that restrain the convective production of smaller and smaller scales of motion in an unconditionally stable manner, meaning that the velocity cannot blow up in the energy-norm (in 2D also: enstrophy-norm). The numerical algorithm used to solve the governing equations preserves the symmetry and conservation properties too. In the present work, a criterion to determine dynamically the regularization parameter (local filter length) is proposed: it is based on the requirement that the vortex stretching must stop at the scale set by the grid. Therefore, the proposed method constitutes a parameter-free turbulence model. The resulting regularization method is tested for a 3D natural convection flow in an air-filled (Pr = 0.71) differentially heated cavity of height aspect ratio 4. Direct comparison with DNS results at Rayleigh number 6.4 Â 10 8 6 Ra 6 10 11 shows fairly good agreement even for very coarse grids. Finally, the robustness of the method is tested by performing simulations with Ra up to 10
17. A 2/7 scaling law of Nusselt number has been obtained for the investigated range of Ra.