High Rayleigh Number Natural Convection Inside 2D Porous Enclosures Using the Lattice Boltzmann Method

2011 ◽  
Vol 133 (6) ◽  
Author(s):  
Ramanathan Vishnampet ◽  
Arunn Narasimhan ◽  
V. Babu

Lattice Boltzmann method (LBM) is employed to investigate natural convection inside porous medium enclosures at high Rayleigh numbers. Volume averaged porous medium model is coupled with the lattice Boltzmann formulation of the momentum and energy equations for fluid flow. A parallel implementation of the single relaxation time LBM is used, which allows the porous medium modified Rayleigh number Ram to be as high as 108. Heat transfer results in the form of the enclosure averaged Nusselt number Nu are obtained for higher modified Rayleigh numbers 104≤Ram≤108. The Nu values are compared with values in the absence of the form drag term. The form drag due to the porous medium is found to influence Nu considerably. The effect of the form drag on Nu is studied by using a form drag modified Rayleigh number RaC (extended from Ram). Utilizing the results for Nu in the high Ram range, a correlation is proposed between Nu and RaC for Darcy numbers 10−6≤Da≤10−2, encompassing the non-Darcy flow regime.

2019 ◽  
Vol 29 (10) ◽  
pp. 3659-3684 ◽  
Author(s):  
Rasul Mohebbi ◽  
Mohsen Izadi ◽  
Nor Azwadi Che Sidik ◽  
Gholamhassan Najafi

Purpose This paper aims to study the natural convection of a nanofluid inside a cavity which contains obstacles using lattice Boltzmann method (LBM). The results have focused mainly on various parameters such as number and aspect ratio of roughness elements and different nanoparticle volume fraction. The isotherms and streamlines are presented to describe the hydrodynamics and thermal behaviors of the nanofluid flow throughout the enclosure. Design/methodology/approach The methodology of this paper consists of mathematical model, statement of the problem, nanofluid thermophysical properties, lattice Boltzmann method, LBM for fluid flow, LBM for heat transfer, numerical strategy, boundary conditions, Nusselt (Nu) number calculation, code validation and grid independence. Findings Natural convection heat transfers of a nanofluid inside cavities with and without rough elements have been studied. Lattice Boltzmann technique has been used as numerical approach. The results showed that at higher Rayleigh number (Ra = 106), there are denser streamlines near the left (source) and right wall (sink) which results in better cooling and enhances convective heat rejection to the heat sink. After a distinctive aspect ratio of rough elements (A = 0.1), change in streamline pattern which arises from increasing of aspect ratio does not have an important effect on isotherms. Results indicate that for lower Rayleigh number (Ra = 103), no variation in average Nu is observed with increasing in number of roughness, while for higher one (Ra = 106) average Nu decreases from N = 0 (smooth cavity) up to N = 4 and then remains constant (N = 6). Originality/value Currently, no argumentative and comprehensive extraction can be concluded without fully understanding the role of different arrangement of roughness. Some geometrical parameters such as aspect ratio, number and position of rough elements have been considered. Also, the effect of nanoparticle concentration was studied at different Ra number. Briefly, using LBM, this paper aims to investigate the natural convection of a nanofluid flow on the thermal and hydrodynamics parameters in the presence of rough element with various arrangements.


2013 ◽  
Vol 2013 ◽  
pp. 1-10 ◽  
Author(s):  
Hassen M. Ouakad

This paper presents an investigation of the density-driven problem that rises during the CO2sequestration into saline aquifer. The lattice Boltzmann method (LBM) is implemented in a way to solve this mixing problem (the brine problem along with the solute transport). The CO2-brine interface was located at the top of the considered domain. Different Rayleigh numbers were used in order to investigate this problem. When Rayleigh number is low, we got steady-state concentration contours describing a Rayleigh-Bénard type of convection. Moreover, when the Rayleigh number was selected to be big enough, we observe that the system is less stable and a convective fingering is initiated. This instability is caused by a higher density difference between the brine and the sequestrated CO2. Note here that the turbulence is not taken into account in the study. After the onset this convective instability, the brine with a high CO2concentration migrates down into the porous medium. This study is based on a statistical LBM theory without assuming periodicity in any directions and without considering any type of disturbances in order to turnon the instability behavior.


2016 ◽  
Vol 26 (4) ◽  
pp. 735-747
Author(s):  
Kai Feldhusen ◽  
Ralf Deiterding ◽  
Claus Wagner

Abstract Utilizing the Boussinesq approximation, a double-population incompressible thermal lattice Boltzmann method (LBM) for forced and natural convection in two and three space dimensions is developed and validated. A block-structured dynamic adaptive mesh refinement (AMR) procedure tailored for the LBM is applied to enable computationally efficient simulations of moderate to high Rayleigh number flows which are characterized by a large scale disparity in boundary layers and free stream flow. As test cases, the analytically accessible problem of a two-dimensional (2D) forced convection flow through two porous plates and the non-Cartesian configuration of a heated rotating cylinder are considered. The objective of the latter is to advance the boundary conditions for an accurate treatment of curved boundaries and to demonstrate the effect on the solution. The effectiveness of the overall approach is demonstrated for the natural convection benchmark of a 2D cavity with differentially heated walls at Rayleigh numbers from 103 up to 108. To demonstrate the benefit of the employed AMR procedure for three-dimensional (3D) problems, results from the natural convection in a cubic cavity at Rayleigh numbers from 103 up to 105 are compared with benchmark results.


Author(s):  
Ayoub Msaddak ◽  
Mohieddine Ben Salah ◽  
Ezeddine Sediki

Lattice Boltzmann method (LBM) is performed to study numerically combined natural convection and surface radiation inside an inclined two-dimensional open square cavity. The cavity is heated by a constant temperature at the wall facing the opening. The walls normal to the heated surface are assumed to be adiabatic, diffuse, gray, and opaque while the open boundary is assumed to be black at ambient temperature. A Bathnagar, Gross and Krook (BGK) collision model with double distribution function (D2Q9-D2Q4) is adopted. Effects of surface radiation, inclination angle, and Rayleigh number on the heat transfer are analyzed and discussed. Results are presented in terms of isotherms, streamlines, and Nusselt number. It was found that the presence of surface radiation enhances the heat transfer. The convective Nusselt number decreases with increasing surface emissivity as well as with Rayleigh number, while the total Nusselt number increases with increasing surface emissivity and Rayleigh number. The inclination angle has also a significant effect on flow and heat transfer inside the cavity. However, the magnitude of total heat transfer decreases considerably when open cavity is tilted downward.


Author(s):  
Mohsen Izadi ◽  
Rasul Mohebbi ◽  
A. Chamkha ◽  
Ioan Pop

PurposeThe purpose of this paper is to consider natural convection of a nanofluid inside of a C-shaped cavity using Lattice Boltzmann method (LBM).Design/methodology/approachEffects of some geometry and flow parameters consisting of the aspect ratio of the cavity, aspect ratio of the heat source; Rayleigh number (Ra = 103− 106) have been investigated. The validity of the method is checked by comparing the present results with ones from the previously published work.FindingsThe results demonstrate that for Ra = 103, the aspect ratio of the heat source has more influence on the average Nusselt number in contrast to the case of Ra = 106. Contrary to the fact that the average Nusselt number increases non-linearly more than twice because of the increase of the aspect ratio of the enclosure at Ra = 103, the average Nusselt number has a linear relation with the aspect ratio for of Ra = 106. Therefore, upon increasing the Rayleigh number, the efficiency of the aspect ratio of the cavity on the thermal convection, gradually diminishes.Originality/valueThe authors believe that all the results, both numerical and asymptotic, are original and have not been published elsewhere.


Author(s):  
Mo Yang ◽  
Yuwei Zhou ◽  
Yuwen Zhang ◽  
Enjie Bian ◽  
Zheng Li

The Lattice Boltzmann method is used to simulate the flow patterns of natural convection in horizontal cylindrical annulus for aspect ratios in the range of 0.4≤A≤10 and for Prandtl numbers varying from 0.1 to 0.7. At Pr = 0.3 and A = 2, flow patterns on the whole range of Rayleigh number are mapped indicating the lower and upper critical values for transitions. At Pr = 0.7, the influence of aspect ratio on flow pattern is analyzed acquiring the result that the oscillation flow never happens at A≤3. At A = 2, several Prandtl numbers are calculated at certain Rayleigh number and the conclusion is that the steady upward flow keeps when 0.5≤Pr. The results are found in good agreement with existed studies.


2015 ◽  
Vol 17 (5) ◽  
pp. 1201-1224 ◽  
Author(s):  
Wei Gong ◽  
Kévyn Johannes ◽  
Frédéric Kuznik

AbstractA new solver is developed to numerically simulate the melting phase change with natural convection. This solver was implemented on a single Nvidia GPU based on the CUDA technology in order to simulate the melting phase change in a 2D rectangular enclosure. The Rayleigh number is of the order of magnitude of 108and Prandlt is 50. The hybrid thermal lattice Boltzmann method (HTLBM) is employed to simulate the natural convection in the liquid phase, and the enthalpy formulation is used to simulate the phase change aspect. The model is validated by experimental data and published analytic results. The simulation results manifest a strong convection in the melted phase and a different flow pattern from the reference results with low Rayleigh number. In addition, the computational performance is estimated for single precision arithmetic, and this solver yields 703.31MLUPS and 61.89GB/s device to device data throughput on a Nvidia Tesla C2050 GPU.


Sign in / Sign up

Export Citation Format

Share Document