Direct Numerical Simulation and RANS Comparison of Turbulent Convective Heat Transfer in a Staggered Ribbed Channel With High Blockage

2016 ◽  
Vol 139 (2) ◽  
Author(s):  
Luca Marocco ◽  
Andrea Franco

A turbulent convective flow of an incompressible fluid inside a staggered ribbed channel with high blockage at ReH ≈ 4200 is simulated with direct numerical simulation (DNS) and Reynolds-averaged Navier–Stokes (RANS) techniques. The DNS results provide the reference solution for comparison of the RANS turbulence models. The k–ε realizable, k–ω SST, and v2¯–f model are accurately analyzed for their strengths and weaknesses in predicting the flow and temperature field for this geometry. These three models have been extensively used in literature to simulate this configuration and boundary conditions but with discordant conclusions upon their performance. The v2¯–f model performs much better than the k–ε realizable while the k–ω SST model results to be inadequate.

2021 ◽  
Vol 9 (3) ◽  
pp. 264
Author(s):  
Shanti Bhushan ◽  
Oumnia El Fajri ◽  
Graham Hubbard ◽  
Bradley Chambers ◽  
Christopher Kees

This study evaluates the capability of Navier–Stokes solvers in predicting forward and backward plunging breaking, including assessment of the effect of grid resolution, turbulence model, and VoF, CLSVoF interface models on predictions. For this purpose, 2D simulations are performed for four test cases: dam break, solitary wave run up on a slope, flow over a submerged bump, and solitary wave over a submerged rectangular obstacle. Plunging wave breaking involves high wave crest, plunger formation, and splash up, followed by second plunger, and chaotic water motions. Coarser grids reasonably predict the wave breaking features, but finer grids are required for accurate prediction of the splash up events. However, instabilities are triggered at the air–water interface (primarily for the air flow) on very fine grids, which induces surface peel-off or kinks and roll-up of the plunger tips. Reynolds averaged Navier–Stokes (RANS) turbulence models result in high eddy-viscosity in the air–water region which decays the fluid momentum and adversely affects the predictions. Both VoF and CLSVoF methods predict the large-scale plunging breaking characteristics well; however, they vary in the prediction of the finer details. The CLSVoF solver predicts the splash-up event and secondary plunger better than the VoF solver; however, the latter predicts the plunger shape better than the former for the solitary wave run-up on a slope case.


Author(s):  
Afaque Shams ◽  
Tomasz Kwiatkowski

Detailed knowledge of a coolant flow in a fuel assembly of a reactor core has always been a major factor in the design of new nuclear systems. In this regard, traditionally adopted subchannel analysis codes cannot take into account local phenomena, which are quite essential. On the other hand, Computational Fluid Dynamic (CFD) is being recognized as a valuable research tool for thermal-hydraulics phenomenon in the fuel assembly geometries. Because of the high Reynolds number and geometric complexities, the practical CFD calculations are mostly limited to pragmatic Reynolds Averaged Navier-Stokes (RANS) type modelling approaches. A good prediction of the flow and heat transport inside the fuel rod bundle is a challenge for such RANS turbulence models and these models need to be validated. Although the measurement techniques are constantly getting improved, however, the CFD-grade experiments of flow mixing and heat transfer in the subchannel scale are often impossible or quite costly to be performed. In addition, lack of experimental databases makes it impossible to validate and/or calibrate the available RANS turbulence models for certain flow situations. In that context, Direct Numerical Simulation (DNS) can serve as a reference for model development and validation. The aim of this work is to design a numerical experiment in order to generate a high quality DNS database for a tight lattice bare rod bundle, which will serve as a reference for the validation purpose. The considered geometric design is based on the well-known Hooper experiment, which contains a bare rod bundle with pitch-to-diameter ratio of P/D = 1.107. Performing a DNS computation corresponding to the Hooper experiment requires a huge computational power. Hence, a wide range of unsteady RANS (URANS) study has been performed to scale-down the Reynolds number such that it is feasible for a DNS computation and at the same time it still preserves the main flow characteristics. In addition to the flow field, a parametric study for three different passive scalars is performed to take into account the heat transfer analysis. These passive scalars correspond to the Prandtl numbers of air, water and liquid metal fluids. The heat transfer of these three fluids has been studied in combination with two different boundary conditions at the walls, i.e. a constant temperature and a constant heat flux. Finally, the obtained URANS results are used to compute the Kolmogorov and Batchelor length scales in order to estimate the overall meshing requirements for the targeted DNS.


2020 ◽  
Author(s):  
P. Srinivasa Rao

The governing equations of viscous fluid flow are generally represented by Navier–Stokes (NS) equations. The output of Navier Stokes equations is in essence velocity vector from which rest of the flow parameters can be calculated. It is essentially a riotous task, sometimes it becomes so unmanageable that fluid flow over simplest topologies under low Reynold’s numbers also needs the most powerful supercomputing facility to solve, if needed to model the fluid and its behavior under the turbulent conditions the best way out is to solve the averaged NS equations. However in the process of averaging Reynolds introduced certain new terms such as Reynolds Stresses. Therefore it is required to close the system of equations by relating the unknown variables with known ones. Hence we have turbulence models. Direct Numerical Simulation (DNS) is a method of solving NS equations directly that is by forfeiting the need of turbulence models as the equations are not averaged. However originally direct numerical simulation procedure does not need of additional closure equations, it is essential to have very fine grid elements and should be estimated for exceptionally small time steps to achieve precise solutions. In the present chapter an interesting flow through nano-channel problem has been discussed using the indispensable mathematical technique of computational fluid dynamics (CFD) which is DNS.


Author(s):  
Alessandro Chiarini ◽  
Maurizio Quadrio

AbstractA direct numerical simulation (DNS) of the incompressible flow around a rectangular cylinder with chord-to-thickness ratio 5:1 (also known as the BARC benchmark) is presented. The work replicates the first DNS of this kind recently presented by Cimarelli et al. (J Wind Eng Ind Aerodyn 174:39–495, 2018), and intends to contribute to a solid numerical benchmark, albeit at a relatively low value of the Reynolds number. The study differentiates from previous work by using an in-house finite-differences solver instead of the finite-volumes toolbox OpenFOAM, and by employing finer spatial discretization and longer temporal average. The main features of the flow are described, and quantitative differences with the existing results are highlighted. The complete set of terms appearing in the budget equation for the components of the Reynolds stress tensor is provided for the first time. The different regions of the flow where production, redistribution and dissipation of each component take place are identified, and the anisotropic and inhomogeneous nature of the flow is discussed. Such information is valuable for the verification and fine-tuning of turbulence models in this complex separating and reattaching flow.


2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.


2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.


2018 ◽  
Vol 74 (6) ◽  
pp. 840-860 ◽  
Author(s):  
Sagar Dave ◽  
Chetankumar Anghan ◽  
Shaswat Saincher ◽  
Jyotirmay Banerjee

2013 ◽  
Vol 135 (3) ◽  
Author(s):  
C. Chin ◽  
M. Li ◽  
C. Harkin ◽  
T. Rochwerger ◽  
L. Chan ◽  
...  

A numerical study of compressible jet flows is carried out using Reynolds averaged Navier–Stokes (RANS) turbulence models such as k-ɛ and k-ω-SST. An experimental investigation is performed concurrently using high-speed optical methods such as Schlieren photography and shadowgraphy. Numerical and experimental studies are carried out for the compressible impinging at various impinging angles and nozzle-to-wall distances. The results from both investigations converge remarkably well and agree with experimental data from the open literature. From the flow visualizations of the velocity fields, the RANS simulations accurately model the shock structures within the core jet region. The first shock cell is found to be constraint due to the interaction with the bow-shock structure for nozzle-to-wall distance less than 1.5 nozzle diameter. The results from the current study show that the RANS models utilized are suitable to simulate compressible free jets and impinging jet flows with varying impinging angles.


Sign in / Sign up

Export Citation Format

Share Document