Towards a Direct Numerical Simulation of a Simplified Pressurized Thermal Shock

The aim of this research work is to perform high quality direct numerical simulations (DNS) of a simplified single phase pressurized thermal shock (PTS) scenario with and without buoyancy effects. In that context, the objectives of this paper are (i) to present the road towards the DNS of a PTS design without buoyancy effects and (ii) to demonstrate that the code NEK5000 is adequate for true DNS analyses. This DNS of the PTS design will serve as a reference to validate low order CFD approaches. The higher order spectral element code NEK5000 is selected to perform the high quality DNS computations. The capabilities of this code, in order to perform the DNS for PTS like geometries, have been extensively assessed for a well-known turbulent channel flow configuration with Re_τ = 180 (turbulent Reynolds number based on the wall friction velocity). Different numerical parameters of NEK5000 have been thoroughly tested and their influence has been studied to obtain high quality turbulence statistics. This assessment of NEK5000 is further extended for the application of highly skewed hexahedral (non-orthogonal) meshes in a turbulent channel flow. The obtained results have shown that NEK5000 is capable of producing high quality DNS solution for a PTS like complex flow configuration for skewed elements (meshes) up to 60 degrees. Finally, this tested numerical framework is adopted to perform the targeted DNS computations of the simplified PTS design.