Immersed volume method for solving natural convection, conduction and radiation of a hat‐shaped disk inside a 3D enclosure

Author(s):  
E. Hachem ◽  
H. Digonnet ◽  
E. Massoni ◽  
T. Coupez

PurposeThe purpose of this paper is to present an immersed volume method that accounts for solid conductive bodies (hat‐shaped disk) in calculation of time‐dependent, three‐dimensional, conjugate heat transfer and fluid flow.Design/methodology/approachThe incompressible Navier‐Stokes equations and the heat transfer equations are discretized using a stabilized finite element method. The interface of the immersed disk is defined and rendered by the zero isovalues of a level set function. This signed distance function allows turning different thermal properties of each component into homogeneous parameters and it is coupled to a direct anisotropic mesh adaptation process enhancing the interface representation. A monolithic approach is used to solve a single set of equations for both fluid and solid with different thermal properties.FindingsIn the proposed immersion technique, only a single grid for both air and solid is considered, thus, only one equation with different thermal properties is solved. The sharp discontinuity of the material properties was captured by an anisotropic refined solid‐fluid interface. The robustness of the method to compute the flow and heat transfer with large materials properties differences is demonstrated using stabilized finite element formulations. Results are assessed by comparing the predictions with the experimental data.Originality/valueThe proposed method demonstrates the capability of the model to simulate an unsteady three‐dimensional heat transfer flow of natural convection, conduction and radiation in a cubic enclosure with the presence of a conduction body. A previous knowledge of the heat transfer coefficients between the disk and the fluid is no longer required. The heat exchange at the interface is solved and dealt with naturally.

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

Purpose The purpose of this paper is to examine numerically the three natural convection of air induced by temperature difference between a cold outer cubic enclosure and a hot inner cylinder. Simulations have been carried out for Rayleigh numbers ranging from 103 to 107 and titled angle of the enclosure from 0° to 90°. The developed mathematical model is governed by the coupled equations of continuity, momentum and energy, and is solved by finite volume method. The effects of cylinder inclination and Rayleigh number on fluid flow and heat transfer are presented. The distribution of isocontours of temperature and isosurfaces of velocity eventually reaches a steady state in the range of Rayleigh numbers between 103 and 107 for titled inclination of 90°; however, for the remaining inclinations, Rayleigh number must be in the range 103-106 to avoid unsteady state, which is manifested by the division of the area containing the maximum local heat transfer rate into three parts for a Rayleigh number equal to 107 and an inclination of 90°. We mention that instability study is not included in the present paper, which is solely devoted to three-dimensional calculations. Results also indicate that optimal average heat transfer rate is obtained for both high Rayleigh number of 106 and high inclination of 90° for the two cases of the inner cylinder and cubical enclosure. Design/methodology/approach The manuscript deals with prediction of the three-dimensional natural convection phenomena in a cubical cavity induced by an isothermal cylinder at the center with different inclinations by simulating the flow using highly numerical methods such as finite volume method. Findings It is found that the local Nusselt number through active walls for titled inclination set at 90°, the symmetry of the flow is conserved and the area containing the maximum heat transfer is divided into three smaller areas situated near the upper portion of the wall, taking the maximum value. That may be due to the preparation of local occurrence of instabilities and bifurcation phenomena that appear for Ra > 107, which is not included in the present paper to save journal space. It was found also that an optimal heat transfer appears when the cylinder orientation becomes vertical (a = 90°). For this inclination, buoyancy forces act upward, corresponding to an aiding situation. In addition, heat transfer rate is increasing with Rayleigh numbers, so correlations of average Nusselt through the cubical cavity and the cylinder are established as function of two parameters (Ra, a). Comparisons of the numerical results with those obtained from all correlations show good agreements. Originality/value To the author’s knowledge, studies have thus far adressed three-dimensional cuboids enclosures induced by an inner shape which the location is changed. However, no study has examined three-dimensional natural convection between the inner isothermal cylinder and outer cooled cubical enclosure when the outer enclosure is tilted.


Author(s):  
M.A. Mansour ◽  
Sameh Elsayed Ahmed ◽  
Ali J. Chamkha

Purpose This paper aims to investigate the entropy generation due to magnetohydrodynamic natural convection flow and heat transfer in a porous enclosure filled with Cu-water nanofluid in the presence of viscous dissipation effect. The left and right walls of the cavity are thermally insulated. There are heated and cold parts, and these are placed on the bottom and top wall, respectively, whereas the remaining parts are thermally insulated. Design/methodology/approach The finite volume method is used to solve the dimensionless partial differential equations governing the problem. A comparison with previously published woks is presented and is found to be in an excellent agreement. Findings The minimization of entropy generation and local heat transfer according to different values of the governing parameters are presented in details. It is found that the presence of magnetic field has negative effects on the local entropy generation because of heat transfer and the local total entropy generation. Also, the increase in the heated part length leads to a decrease in the local Nusselt number. Originality/value This problem is original, as it has not been considered previously.


2016 ◽  
Vol 831 ◽  
pp. 83-91
Author(s):  
Lahoucine Belarche ◽  
Btissam Abourida

The three-dimensional numerical study of natural convection in a cubical enclosure, discretely heated, was carried out in this study. Two heating square sections, similar to the integrated electronic components, are placed on the vertical wall of the enclosure. The imposed heating fluxes vary sinusoidally with time, in phase and in opposition of phase. The temperature of the opposite vertical wall is maintained at a cold uniform temperature and the other walls are adiabatic. The governing equations are solved using Control volume method by SIMPLEC algorithm. The sections dimension ε = D / H and the Rayleigh number Ra were fixed respectively at 0,35 and 106. The average heat transfer and the maximum temperature on the active portions will be examined for a given set of the governing parameters, namely the amplitude of the variable temperatures a and their period τp. The obtained results show significant changes in terms of heat transfer, by proper choice of the heating mode and the governing parameters.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Mahyar Ashouri ◽  
Mohammad Mehdi Zarei ◽  
Ali Moosavi

Purpose The purpose of this paper is to investigate the effects of geometrical parameters, eccentricity and perforated fins on natural convection heat transfer in a finned horizontal annulus using three-dimensional lattice Boltzmann flux solver. Design/methodology/approach Three-dimensional lattice Boltzmann flux solver is used in the present study for simulating conjugate heat transfer within an annulus. D3Q15 and D3Q7 models are used to solve the fluid flow and temperature field, respectively. The finite volume method is used to discretize mass, momentum and energy equations. The Chapman–Enskog expansion analysis is used to establish the connection between the lattice Boltzmann equation local solution and macroscopic fluxes. To improve the accuracy of the lattice Boltzmann method for curved boundaries, lattice Boltzmann equation local solution at each cell interface is considered to be independent of each other. Findings It is found that the maximum heat transfer rate occurs at low fin spacing especially by increasing the fin height and decreasing the internal-cylindrical distance. The effect of inner cylinder eccentricity is not much considerable (up to 5.2% enhancement) while the impact of fin eccentricity is more remarkable. Negative fin eccentricity further enhances the heat transfer rate compared to a positive fin eccentricity and the maximum heat transfer enhancement of 91.7% is obtained. The influence of using perforated fins is more considerable at low fin spacing although some heat transfer enhancements are observed at higher fin spacing. Originality/value The originality of this paper is to study three-dimensional natural convection in a finned-horizontal annulus using three-dimensional lattice Boltzmann flux solver, as well as to apply symmetry and periodic boundary conditions and to analyze the effect of eccentric annular fins (for the first time for air) and perforated annular fins (for the first time so far) on the heat transfer rate.


Author(s):  
Alireza Rahimi ◽  
Aravindhan Surendar ◽  
Aygul Z. Ibatova ◽  
Abbas Kasaeipoor ◽  
Emad Hasani Malekshah

Purpose This paper aims to investigate the three-dimensional natural convection and entropy generation in the rectangular cuboid cavities included by chamfered triangular partition made by polypropylene. Design/methodology/approach The enclosure is filled by multi-walled carbon nanotubes (MWCNTs)-H2O nanofluid and air as two immiscible fluids. The finite volume approach is used for computation. The fluid flow and heat transfer are considered with combination of local entropy generation due to fluid friction and heat transfer. Moreover, a numerical method is developed based on three-dimensional solution of Navier–Stokes equations. Findings Effects of side ratio of triangular partitions (SR = 0.5, 1 and 2), Rayleigh number (103 < Ra < 105) and solid volume fraction (f = 0.002, 0.004 and 0.01 Vol.%) of nanofluid are investigated on both natural convection characteristic and volumetric entropy generation. The results show that the partitions can be a suitable method to control fluid flow and energy consumption, and three-dimensional solutions renders more accurate results. Originality/value The originality of this work is to study the three-dimensional natural convection and entropy generation of a stratified system.


2019 ◽  
Vol 29 (8) ◽  
pp. 2792-2808 ◽  
Author(s):  
Behnam Rafiei ◽  
Hamed Masoumi ◽  
Mohammad Saeid Aghighi ◽  
Amine Ammar

Purpose The purpose of this paper is to analyze the effects of complex boundary conditions on natural convection of a yield stress fluid in a square enclosure heated from below (uniformly and non-uniformly) and symmetrically cooled from the sides. Design/methodology/approach The governing equations are solved numerically subject to continuous and discontinuous Dirichlet boundary conditions by Galerkin’s weighted residuals scheme of finite element method and using a non-uniform unstructured triangular grid. Findings Results show that the overall heat transfer from the heated wall decreases in the case of non-uniform heating for both Newtonian and yield stress fluids. It is found that the effect of yield stress on heat transfer is almost similar in both uniform and non-uniform heating cases. The yield stress has a stabilizing effect, reducing the convection intensity in both cases. Above a certain value of yield number Y, heat transfer is only due to conduction. It is found that a transition of different modes of stability may occur as Rayleigh number changes; this fact gives rise to a discontinuity in the variation of critical yield number. Originality/value Besides the new numerical method based on the finite element and using a non-uniform unstructured grid for analyzing natural convection of viscoplastic materials with complex boundary conditions, the originality of the present work concerns the treatment of the yield stress fluids under the influence of complex boundary conditions.


Author(s):  
Leo A. Carrilho ◽  
Jamil Khan ◽  
Michael E. Conner ◽  
Abdel Mandour ◽  
Milorad B. Dzodzo

The effects of artificial roughness for the purpose of thermal performance improvement in pressurized water nuclear reactors are investigated. The artificial roughness consists of two-dimensional ribs parallel to the turbulent flow. The fuel rod bundle subchannel is preliminarily modeled as an annulus using the finite element method in ANSYS/FLOTRAN. The Navier-Stokes equations are solved from the SST (Shear Stress Transport) turbulence model for the simulated annulus thermal-flow. The analyses are performed for ribs dimensions and pitch provided by published previous work. It is found that, heat transfer and differential pressure have similar behavior with highest heat transfer occurring at the reattachment point. The finite element model describes well the characteristics of turbulent flow in smooth and rough rod when compared to previous semi-empirical models. Next paper extends the analysis by comparing numerical results with experimental test data and sensitivity analyses for different roughness configurations.


2019 ◽  
Vol 29 (7) ◽  
pp. 2352-2388
Author(s):  
Jory Seguin ◽  
Song Gao ◽  
Wagdi George Habashi ◽  
Dario Isola ◽  
Guido Baruzzi

Purpose This paper aims to describe the physical and numerical modeling of a new computational fluid dynamics solver for hypersonic flows in thermo-chemical non-equilibrium. The code uses a blend of numerical techniques to ensure accuracy and robustness and to provide scalability for advanced hypersonic physics and complex three-dimensional (3D) flows. Design/methodology/approach The solver is based on an edge-based stabilized finite element method (FEM). The chemical and thermal non-equilibrium systems are loosely-coupled to provide flexibility and ease of implementation. Chemical non-equilibrium is modeled using a laminar finite-rate chemical kinetics model while a two-temperature model is used to account for thermodynamic non-equilibrium. The systems are solved implicitly in time to relax numerical stiffness. Investigations are performed on various canonical hypersonic geometries in two-dimensional and 3D. Findings The comparisons with numerical and experimental results demonstrate the suitability of the code for hypersonic non-equilibrium flows. Although convergence is shown to suffer to some extent from the loosely-coupled implementation, trading a fully-coupled system for a number of smaller ones improves computational time. Furthermore, the specialized numerical discretization offers a great deal of flexibility in the implementation of numerical flux functions and boundary conditions. Originality/value The FEM is often disregarded in hypersonics. This paper demonstrates that this method can be used successfully for these types of flows. The present findings will be built upon in a later paper to demonstrate the powerful numerical ability of this type of solver, particularly with respect to robustness on highly stretched unstructured anisotropic grids.


Author(s):  
Alireza Rahimi ◽  
Abbas Kasaeipoor ◽  
Emad Hasani Malekshah ◽  
Mohammad Mehdi Rashidi ◽  
Abimanyu Purusothaman

Purpose This study aims to investigate the three-dimensional natural convection and entropy generation in a cuboid enclosure filled with CuO-water nanofluid. Design/methodology/approach The lattice Boltzmann method is used to solve the problem numerically. Two different multiple relaxation time (MRT) models are used to solve the problem. The D3Q7–MRT model is used to solve the temperature field, and the D3Q19 is used to solve the fluid flow of natural convection within the enclosure. Findings The influences of different Rayleigh numbers (103 < Ra < 106) and solid volume fractions (0 < f < 0.04) on the fluid flow, heat transfer, total entropy generation, local heat transfer irreversibility and local fluid friction irreversibility are presented comprehensively. To predict thermo–physical properties, dynamic viscosity and thermal conductivity, of CuO–water nanofluid, the Koo–Kleinstreuer–Li (KKL) model is applied to consider the effect of Brownian motion on nanofluid properties. Originality/value The originality of this work is to analyze the three-dimensional natural convection and entropy generation using a new numerical approach of dual-MRT-based lattice Boltzmann method.


Sign in / Sign up

Export Citation Format

Share Document