Computational Simulation of Turbulent Natural Convection in a Volumetrically Heated Square Cavity

Author(s):  
Camila Braga Vieira ◽  
Bojan Niceno ◽  
Jian Su

This work aimed to analyze the turbulent natural convection in a volumetrically heated fluid with Prandtl number equal to 0.6, representing the oxide material layer of a corium. Four turbulence models were scrutinized in order to select the most appropriate one for turbulence modeling based on Reynolds Averaged Navier-Stokes equations (RANS) of natural convection in a molten core. The turbulence models scrutinized are the standard k-ε, Shear Stress Transport (SST), low-Reynolds-k-ε (Launder-Sharma) and also an elliptic blending model ν2-f. The simulations were carried out in a square cavity with isothermal walls, for Rayleigh numbers (Ra) ranging from 109 to 1011. The numerical simulations, performed in an open-source of Computational Fluid Dynamics (CFD) - OpenFOAM (Open Field Operation and Manipulation), provided outcomes of average Nusselt number as function of Ra number, which were in a reasonable agreement with an experimental correlation and other authors’ simulations. It was also possible to observe the limitations and robustness of each model analyzed, enabling to conclude that the most adequate turbulence models for the present physical problem were SST and ν2-f.

2021 ◽  
Vol 39 (1) ◽  
pp. 227-234
Author(s):  
Khelifa Hami

This contribution represents a critical view of the advantages and limits of the set of mathematical models of the physical phenomena of turbulence. Turbulence models can be grouped into two categories, depending on how turbulent quantities are calculated: direct numerical simulations (DNS) and RANS (Reynolds Averaged Navier-Stokes Equations) models. The disadvantage of these models is that they require enormous computing power, inaccessible, especially for large and complicated geometries. For this reason, hybrid models (combinations between DNS and RANS methods) have been developed, for example, the LES (“Large Eddy Simulation”) or DES (“Detached Eddy Simulation”) models. They represent a compromise - are less precise than DNS, but more precise than RANS models. The results presented in this contribution will allow and facilitate future research in the field the choice of the model approach necessary for the case studies whatever their difficulty factor.


Author(s):  
Juan B. V. Wanderley ◽  
Gisele H. B. Souza ◽  
Carlos Levi

Author’s previous work Wanderley [1] presented an efficient numerical method to investigate VIV phenomenon on circular cylinders. The numerical model solves the unsteady Reynolds Average Navier–Stokes equations for slightly compressible flows using the Beam–Warming implicit factored scheme. In the present work, the effect of the turbulence model on the results is evaluated for both Baldwin Lomax and k-ε models. To demonstrate the quality of the numerical method, results for the transversal oscillation of a cylinder laterally supported by spring and damper are compared with experimental data. The application of the turbulence models showed the much better agreement of the k-ε model with the experimental results.


Author(s):  
Basma Souayeh ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya ◽  
Taieb Lili

The present work deals with the prediction of a natural convection flow in a square cavity, partially heated by an obstacle placed at the bottom wall. The two transverse walls and the top wall of the cavity are supposed to be cold, the remaining walls are kept insulated. The main parameter of numerical investigations is the Rayleigh number (engine convection) varying from 103 to 105. When Ra is fixed at 107, the flow and thermal fields bifurcate and undergoes an unsteady behavior at critical positions. Flow patterns corresponding to the unsteady state are presented and analyzed in the current study. The simulations were conducted using a numerical approach based on the finite volume method and the projection method, which are implemented in a computer code in order to solve the Navier-Stokes equations.


2001 ◽  
Author(s):  
Masaru Ishizuka ◽  
Guoyi Peng ◽  
Shinji Hayama

Abstract In the present work, an important basic flow phenomena, the natural convection induced flow, is studied numerically. Three-dimensional Navier-Stokes equations along with the temperature equation are solved on the basis of finite difference method. Generalized coordinate system is used so that sufficient grid resolution could be achieved in the body surface boundary layer region. Differential terms with respect to time are approximated by forward differences, diffusions terms are approximated by the implicit Euler form, convection terms in the Navier-Stokes equations are approximated by the third order upwind difference scheme. The heat flux at the body surface of heater is specified. The results of calculation showed a satisfactory agreement with the measured data and led to a good understanding of the overall flow and thermal behavior inside electronic equipment cabinet model which is very difficult, if not impossible, to gather by experiment.


2018 ◽  
Vol 855 ◽  
pp. 43-66 ◽  
Author(s):  
Ke Wu ◽  
Bruno D. Welfert ◽  
Juan M. Lopez

The dynamic response to shear of a fluid-filled square cavity with stable temperature stratification is investigated numerically. The shear is imposed by the constant translation of the top lid, and is quantified by the associated Reynolds number. The stratification, quantified by a Richardson number, is imposed by maintaining the temperature of the top lid at a higher constant temperature than that of the bottom, and the side walls are insulating. The Navier–Stokes equations under the Boussinesq approximation are solved, using a pseudospectral approximation, over a wide range of Reynolds and Richardson numbers. Particular attention is paid to the dynamical mechanisms associated with the onset of instability of steady state solutions, and to the complex and rich dynamics occurring beyond.


2018 ◽  
Vol 240 ◽  
pp. 01006 ◽  
Author(s):  
Nadezhda Bondareva ◽  
Mikhail Sheremet

Present study is devoted to numerical simulation of heat and mass transfer inside a cooper profile filled with paraffin enhanced with Al2O3 nanoparticles. This profile is heated by the heat-generating element of constant volumetric heat flux. Two-dimensional approximation of melting process is described by the Navier-Stokes equations in non-dimensional variables such as stream function, vorticity and temperature. The enthalpy formulation has been used for description of the heat transfer. The influence of volume fraction of nanoparticles and intensity of heat generation on melting process and natural convection in liquid phase has been studied.


Author(s):  
Jiho You ◽  
Jinmo Lee ◽  
Donghyun You

A computational simulation methodology, which combines a computational fluid dynamics technique and a computational structural dynamics technique, is employed to design a deformable foil of which kinematics is inspired by the propulsive motion of a fin or a tail of fish and cetacean. The unsteady incompressible Navier-Stokes equations are solved using a second-order accurate finite-difference method and an immersed-boundary method to effectively impose boundary conditions on complex moving boundaries. A finite-element-based structural dynamics solver is employed to compute the deformation of the foil due to interaction with fluid. A phase angle between pitching and heaving motions as well as the flexibility of the foil, which is represented by the Youngs modulus are varied to find out how these factors affect the propulsion efficiency.


Sign in / Sign up

Export Citation Format

Share Document