Lattice Boltzmann Modelling for Natural Convection in Power-Law Fluids within a Partially Heated Square Enclosure

2020 ◽  
Author(s):  
Siva Subrahmanyam Mendu ◽  
P.K. Das

Abstract The present paper reports the numerical investigations for steady-state natural convection in power-law fluids inside a square enclosure embedded with bottom discrete heaters. The Lattice Boltzmann Method (LBM) is employed to model the flow and heat transfer phenomenon at different combinations of power-law index, Rayleigh number, and heat source length for a constant Prandtl number. The buoyancy force is incorporated in the collision term of the LBM through Boussinesq approximation. Simulation outcomes are furnished using streamlines and, temperature contours, velocity profiles and variation of heat transfer on the non-adiabatic walls to probe natural convection phenomena. The results indicate that the temperature and the flow fields in the enclosure are symmetric about the vertical centerline. The detailed physical interpretations have been provided for the reported results. Further, the increase in the power-law index means a rise in viscosity and a decrease in thermal energy transport for other constant parameters. The outcomes also specify that the intensity of circulation and heat transfer enhances with the increase of Rayleigh number and size of the localized heater. Finally, though, a rise in the size of the confined heat source enhances the rate of total thermal transport, it does not change the trend of fluid flow and local heat transfer rate.

2018 ◽  
Vol 29 (11) ◽  
pp. 1850105
Author(s):  
Dinggen Li ◽  
Haifeng Zhang ◽  
Peixin Ye ◽  
Zihao Yu

In this paper, natural convection of power-law Al2O3-water nanofluids with temperature-dependent properties in a square enclosure with a circular cylinder is studied. The governing equations of the flow and temperature fields are solved by the lattice Boltzmann method (LBM), and the curved velocity and thermal boundary conditions are treated by immersed boundary method (IBM). The effects of Rayleigh number, power-law index, nanoparticle volume fractions, radius of circular cylinder, nanoparticle diameter and temperature difference on flow and heat transfer characteristics are discussed in detail. The results indicate that the heat transfer rate is increased with the increases of Rayleigh number, radius of circular cylinder and temperature difference, while it generally decreases with an increase in power-law index and nanoparticle diameter. Additionally, it is observed that there is an optimal volume fraction at which the maximum heat transfer enhancement is obtained, and the value of it is found to increase slightly with decreasing the nanoparticle diameter, and to increase remarkably with increasing the temperature difference.


2017 ◽  
Vol 9 (5) ◽  
pp. 1094-1110
Author(s):  
Lei Wang ◽  
Zhenhua Chai ◽  
Baochang Shi

AbstractIn this paper, the magnetic field effects on natural convection of power-law nanofluids in rectangular enclosures are investigated numerically with the lattice Boltzmann method. The fluid in the cavity is a water-based nanofluid containing Cu nanoparticles and the investigations are carried out for different governing parameters including Hartmann number (0.0≤Ha≤20.0), Rayleigh number (104≤Ra≤106), power-law index (0.5≤n≤1.0), nanopartical volume fraction (0.0≤ϕ≤0.1) and aspect ratio (0.125≤AR≤8.0). The results reveal that the flow oscillations can be suppressed effectively by imposing an external magnetic field and the augmentation of Hartmann number and power-law index generally decreases the heat transfer rate. Additionally, it is observed that the average Nusselt number is increased with the increase of Rayleigh number and nanoparticle volume fraction. Moreover, the present results also indicate that there is a critical value for aspect ratio at which the impact on heat transfer is the most pronounced.


Author(s):  
Mustapha Faraji ◽  
El Mehdi Berra

Abstract This paper reported the mathematical modeling and numerical simulation of natural convection flow of Cu/water nanofluid in a square enclosure using the lattice Boltzmann method (LBM). The cavity is heated from below by heat source and cooled by the top wall. The vertical walls are adiabatic. After validating the numerical code against the numerical and experimental data, simulations were performed for different Rayleigh numbers (104–0.5 × 107), nanoparticles volume fractions (0–8%), and cavity inclination angle (0 deg–90 deg). The effects of the studied parameters on the streamlines, on isotherms distributions within the enclosure, and on the local and average Nusselt numbers are investigated. It was found that heat transfer and fluid flow structure depend closely on the nanoparticle concentration. Results show differences in stream separation between a base fluid and the nanofluid. Also, adding small nanoparticles fractions, less than 6%, to the base fluid enhances the heat transfer for higher Rayleigh numbers and cavity inclination angle less than 30 deg. It is concluded that the optimal dilute suspension of copper nanoparticles can be applied as a passive way to enhance heat transfer in natural convection engineering applications.


2021 ◽  
Vol 39 (5) ◽  
pp. 1405-1416
Author(s):  
Hamza Daghab ◽  
Mourad Kaddiri ◽  
Said Raghay ◽  
Ismail Arroub ◽  
Mohamed Lamsaadi ◽  
...  

In this paper, numerical study on natural convection heat transfer for confined thermo-dependent power-law fluids is conducted. The geometry of interest is a fluid-filled square enclosure where a uniform flux heating element embedded on its lower wall is cooled from the vertical walls while the remaining parts of the cavity are insulated, without slipping conditions at all the solid boundaries. The governing partial differential equations written in terms of non-dimensional velocities, pressure and temperature formulation with the corresponding boundary conditions are discretized using a finite volume method in a staggered grid system. Coupled equations of conservation are solved through iterative Semi Implicit Method for Pressure Linked Equation (SIMPLE) algorithm. The effects of pertinent parameters, which are Rayleigh number (103 ≤ Ra ≤ 106), power-law index (0.6 ≤ n ≤ 1.4), Pearson number (0 ≤ m ≤ 20) and length of the heat source (0.2 ≤ W ≤ 0.8) on the cooling performance are investigated. The results indicate that the cooling performance of the enclosure is improved with increasing Pearson and Rayleigh numbers as well as with decreasing power-law index and heat source length.


2017 ◽  
Vol 140 (2) ◽  
Author(s):  
Jacob K. Mulamootil ◽  
Sukanta K. Dash

Natural convection heat transfer from an array of horizontal rectangular fins on a vertical flat plate in non-Newtonian power-law fluids has been studied. The underlying physical principles affecting heat transfer were studied using comprehensive solutions obtained from numerical investigations. Heat transfer to the power-law fluid was found to depend on the fluid rheology (power-law index) and significantly on the geometric parameters (interfin spacing, fin length) as well. The dependence was quantified using the Nusselt number (Nu) and fin effectiveness (Q/Q0). The present study shows that compared to a fin analyzed in isolation, the spatial arrangement of multiple fins relative to one another in an array does have a significant effect on the flow field around subsequent fins in power-law fluids. Therefore, the average heat transfer coefficient of the natural convection system is affected significantly. The variation of Nu with the dimensionless fin length (l/L), dimensionless interfin spacing (S/L), and fluid power-law index (n) was plotted. The dependence was found to be counter intuitive to expectations based on studies for natural convection from vertical flat plates to power-law fluids. In the present study involving fins, shear-thinning fluids (n < 1) show a decrease in heat transfer and shear-thickening fluids (n > 1) show an enhancement in heat transfer for higher l/L values. The results of the study may be useful in the design of natural convection systems that employ power-law fluids to enhance or control heat transfer.


2019 ◽  
Vol 141 (6) ◽  
Author(s):  
Noureddine Abouricha ◽  
Mustapha El Alami ◽  
Ayoub Gounni

Turbulent natural convection in a large-scale cavity has taken a great attention thanks to its importance in many engineering applications such as building. In this work, the lattice Boltzmann method (LBM) is used to simulate turbulent natural convection heat transfer in a small room of housing heated from below by means of a heated floor. The ceiling and the four vertical walls of the room are adiabatic except for a portion of one vertical wall. This portion simulates a glass door with a cold temperature θc = 0. The cavity is filled by air (Pr = 0.71) and heated from below with uniformly imposed temperature θh = 1. The effects of the heat source length (Lr) and Rayleigh number (Ra) on the flow structure and heat transfer are studied for ranges of 0.2 ≤ Lr ≤ 0.8 and 5 × 106 ≤Ra ≤ 108. The heat transfer is examined in terms of local and mean Nusselt numbers. The results show that an increase in Rayleigh number or in heat source length increases the temperature in the core of the cavity. The flow structure shows that turbulent natural convection regime is fully developed for Ra = 108. Correlations for mean Nusselt number as a function with Ra for different values of Lr are expressly derived.


Author(s):  
S. M. Dash ◽  
S. Sahoo

In this article, the natural convection process in a two-dimensional cold square enclosure is numerically investigated in the presence of two inline square heat sources. Two different heat source boundary conditions are analyzed, namely, case 1 (when one heat source is hot) and case 2 (when two heat sources are hot), using the in-house developed flexible forcing immersed boundary–thermal lattice Boltzmann model. The isotherms, streamlines, local, and surface-averaged Nusselt number distributions are analyzed at ten different vertical eccentric locations of the heat sources for Rayleigh number between 103 and 106. Distinct flow regimes including primary, secondary, tertiary, quaternary, and Rayleigh–Benard cells are observed when the mode of heat transfer is changed from conduction to convection and heat sources eccentricity is varied. For Rayleigh number up to 104, the heat transfer from the enclosure is symmetric for the upward and downward eccentricity of the heat sources. At Rayleigh number greater than 104, the heat transfer from the enclosure is better for downward eccentricity cases that attain a maximum when the heat sources are near the bottom enclosure wall. Moreover, the heat transfer rate from the enclosure in case 2 is nearly twice that of case 1 at all Rayleigh numbers and eccentric locations. The correlations for heat transfer are developed by relating Nusselt number, Rayleigh number, and eccentricity of the heat sources.


Author(s):  
Naushad Hasin Khan ◽  
M. A. Hassan

The numerical investigation of laminar natural convection of viscoplastic fluid in a two dimensional square enclosure has been reported in this work. The enclosed fluid is subjected to partial bi-heating from the bottom wall and symmetrical cooling from the sides under steady condition. Yield stress fluid has been heated through two heaters symmetrically placed on the either side of the centre of the bottom wall of the square enclosure. The viscoplastic fluid is the one which requires a minimum critical stress called yield stress to flow otherwise behave as a solid, have been modeled with Herschel–Bulkley model. Such fluids have significant technological relevance due to its wide application ranging from cosmetics products, food processing industries, pharmaceuticals to natural occurring like flow of debris and lava. The solution of governing partial differential equations has been approached using finite volume based formulation. Non uniform set of grid has been used. The effects of yield stress, heat flux, and power law index on the flow and thermal characteristics of the free convection of Herschel-Bulkley fluids have been studied for a particular value of Prandtl number. The flow and thermal fields have been investigated for the following ranges of conditions: Rayleigh number varies between 103 and 106 whereas power law index ranges from 0 to 1. The heat transfer characteristic has been depicted with the help of isotherms and the flow field has been illustrated by streamlines. The onset of convection is substantially delayed due to presence of yield stress of the fluid. This results in enhanced critical Rayleigh number for onset of convection. With increase in the Yield number in turn yield stress, results in the weakening of heat transfer through convection and at a particular relatively higher value of Yield number the heat transfer is solely taken place by conduction mode. Due to the symmetry in both heating and boundary conditions, the obtained isotherms and streamlines of the right half are symmetrical to the left half of the square enclosure. The conductive mode of heat transfer becomes dominated by increasing yield stress and reducing Ra and vice versa. The simultaneous presence of yielded and unyielded region presents an interesting pattern in the convection zone. Furthermore, it can be seen that rise in heat flux, in turn Ra, promotes the buoyancy driven circulation of viscoplastic fluid i.e. enhances natural convective heat transfer. In addition, the effect of power law index has been investigated. Power law index has little effect on thermal distribution and flow field.


2021 ◽  
Vol 3 (3) ◽  
Author(s):  
Mahbuba Tasmin ◽  
Preetom Nag ◽  
Zarin T. Hoque ◽  
Md. Mamun Molla

AbstractA numerical study on heat transfer and entropy generation in natural convection of non-Newtonian nanofluid flow has been explored within a differentially heated two-dimensional wavy porous cavity. In the present study, copper (Cu)–water nanofluid is considered for the investigation where the specific behavior of Cu nanoparticles in water is considered to behave as non-Newtonian based on previously established experimental results. The power-law model and the Brinkman-extended Darcy model has been used to characterize the non-Newtonian porous medium. The governing equations of the flow are solved using the finite volume method with the collocated grid arrangement. Numerical results are presented through streamlines, isotherms, local Nusselt number and entropy generation rate to study the effects of a range of Darcy number (Da), volume fractions (ϕ) of nanofluids, Rayleigh numbers (Ra), and the power-law index (n). Results show that the rate of heat transfer from the wavy wall to the medium becomes enhanced by decreasing the power-law index but increasing the volume fraction of nanoparticles. Increase of porosity level and buoyancy forces of the medium augments flow strength and results in a thinner boundary layer within the cavity. At negligible porosity level of the enclosure, effect of volume fraction of nanoparticles over thermal conductivity of the nanofluids is imperceptible. Interestingly, when the Darcy–Rayleigh number $$Ra^*\gg 10$$ R a ∗ ≫ 10 , the power-law effect becomes more significant than the volume fraction effect in the augmentation of the convective heat transfer process. The local entropy generation is highly dominated by heat transfer irreversibility within the porous enclosure for all conditions of the flow medium. The particular wavy shape of the cavity strongly influences the heat transfer flow pattern and local entropy generation. Interestingly, contour graphs of local entropy generation and local Bejan number show a rotationally symmetric pattern of order two about the center of the wavy cavity.


2017 ◽  
Vol 2017 ◽  
pp. 1-11 ◽  
Author(s):  
L. El Moutaouakil ◽  
Z. Zrikem ◽  
A. Abdelbaki

Laminar natural convection in a water filled square enclosure containing at its center a horizontal hexagonal cylinder is studied by the lattice Boltzmann method. The hexagonal cylinder is heated while the walls of the cavity are maintained at the same cold temperature. Two orientations are treated, corresponding to two opposite sides of the hexagonal cross-section which are horizontal (case I) or vertical (case II). For each case, the results are presented in terms of streamlines, isotherms, local and average convective heat transfers as a function of the dimensionless size of the hexagonal cylinder cross-section (0.1≤B≤0.4), and the Rayleigh number (103≤Ra≤106).


Sign in / Sign up

Export Citation Format

Share Document