Numerical validation of a κ-ω-κθ-ωθ heat transfer turbulence model for heavy liquid metals

33rd UIT (Italian Union of Thermo-fluid dynamics) Heat Transfer Conference

物理学;力学

Cerroni, D.^1 ; Da Viá, R.^1 ; Manservisi, S.^1 ; Menghini, F.^1 ; Pozzetti, G.^2 ; Scardovelli, R.^1

DIN - Lab. of Montecuccolino, University of Bologna, Via dei Colli 16, Bologna

40136, Italy^1

RUES, University of Luxemburg, 6, rue Richard Coudenhove-Kalergi, Luxembourg

L-1359, Luxembourg^2

关键词: Computational Fluid Dynamics codes;    Experimental correlation;    Heat transfer model;    Molecular diffusion;    Numerical instability;    Numerical validations;    Stability and convergence;    Turbulent Prandtl number;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/655/1/012046/pdf
DOI  :  10.1088/1742-6596/655/1/012046

学科分类：力学，机械学

来源: IOP

【 摘 要 】

The correct prediction of heat transfer in turbulent flows is relevant in almost all industrial applications but many of the heat transfer models available in literature are validated only for ordinary fluids with Pr1. In commercial Computational Fluid Dynamics codes only turbulence models with a constant turbulent Prandtl number of 0.85 - 0.9 are usually implemented but in heavy liquid metals with low Prandtl numbers it is well known that these models fail to reproduce correlations based on experimental data. In these fluids heat transfer is mainly due to molecular diffusion and the time scales of temperature and velocity fields are rather different, so simple turbulence models based on similarity between temperature and velocity cannot reproduce experimental correlations. In order to reproduce experimental results and Direct Numerical Simulation data obtained for fluids with Pr0.025 we introduce a κ-Ε-κθ-Εθturbulence model. This model, however, shows some numerical instabilities mainly due to the strong coupling between κ and Ε on the walls. In order to fix this problem we reformulate the model into a new four parameter κ-ω-κθ-ωθwhere the dissipation rate on the wall is completely independent on the fluctuations. The model improves numerical stability and convergence. Numerical simulations in plane and channel geometries are reported and compared with experimental, Direct Numerical Simulation results and with results obtained with the κ-Ε formulation, in order to show the model capabilities and validate the improved κ-ω model.