scholarly journals Natural convection of Casson fluid in a square enclosure

2020 ◽  
Vol 16 (5) ◽  
pp. 1245-1259
Author(s):  
Mohammad Saeid Aghighi ◽  
Christel Metivier ◽  
Hamed Masoumi

PurposeThe purpose of this paper is to analyze the natural convection of a yield stress fluid in a square enclosure with differentially heated side walls. In particular, the Casson model is considered which is a commonly used model.Design/methodology/approachThe coupled conservation equations of mass, momentum and energy related to the two-dimensional steady-state natural convection within square enclosures are solved numerically by using the Galerkin's weighted residual finite element method with quadrilateral, eight nodes elements.FindingsResults highlight a small degree of the shear-thinning in the Casson fluids. It is shown that the yield stress has a stabilizing effect since the convection can stop for yield stress fluids while this is not the case for Newtonian fluids. The heat transfer rate, velocity and Yc obtained with the Casson model have the smallest values compared to other viscoplastic models. Results highlight a weak dependence of Yc with the Rayleigh number: Yc∼Ra0.07. A supercritical bifurcation at the transition between the convective and the conductive regimes is found.Originality/valueThe originality of the present study concerns the comprehensive and detailed solutions of the natural convection of Casson fluids in square enclosures with differentially heated side walls. It is shown that there exists a major difference between the cases of Casson and Bingham models, and hence using the Bingham model for analyzing the viscoplastic behavior of the fluids which follow the Casson model (such as blood) may not be accurate. Finally, a correlation is proposed for the mean Nusselt number Nu¯.

Author(s):  
Sahin Yigit ◽  
Nilanjan Chakraborty

PurposeThis paper aims to numerically analyse natural convection of yield stress fluids in rectangular cross-sectional cylindrical annular enclosures. The laminar steady-state simulations have been conducted for a range of different values of normalised internal radius (ri/L1/8 to 16, whereLis the difference between outer and inner radii); aspect ratio (AR=H/Lfrom 1/8 to 8 whereHis the enclosure height); and nominal Rayleigh number (Rafrom 103to 106) for a single representative value of Prandtl number (Pris 500).Design/methodology/approachThe Bingham model has been used to mimic the yield stress fluid motion, and numerical simulations have been conducted for both constant wall temperature (CWT) and constant wall heat flux (CWHF) boundary conditions for the vertical side walls. The conservation equations of mass, momentum and energy have been solved in a coupled manner using the finite volume method where a second-order central differencing scheme is used for the diffusive terms and a second-order up-wind scheme is used for the convective terms. The well-known semi-implicit method for pressure-linked equations algorithm is used for the coupling of the pressure and velocity.FindingsIt is found that the mean Nusselt number based on the inner peripheryNu¯iincreases (decreases) with an increase inRa(Bn) due to augmented buoyancy (viscous) forces irrespective of the boundary condition. The ratio of convective to diffusive thermal transport increases with increasingri/Lfor both Newtonian (i.e.Bn= 0) and Bingham fluids regardless of the boundary condition. Moreover, the mean Nusselt numberNu¯inormalised by the corresponding Nusselt number due to pure conductive transport (i.e.Nu¯i/(Nu¯i)cond) shows a non-monotonic trend with increasingARin the CWT configuration for a given set of values ofRa,Pr,Lifor both Newtonian (i.e.Bn= 0) and Bingham fluids, whereasNu¯i/(Nu¯i)condincreases monotonically with increasingARin the CWHF configuration. The influences of convective thermal transport strengthen while thermal diffusive transport weakens with increasingAR, and these competing effects are responsible for the non-monotonicNu¯i/(Nu¯i)condvariation withARin the CWT configuration.Originality/valueDetailed scaling analysis is utilised to explain the observed influences ofRa,BN,ri/LandAR, which along with the simulation data has been used to propose correlations forNu¯i.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Nur Husnina Saadun ◽  
Nurul Aini Jaafar ◽  
Md Faisal Md Basir ◽  
Ali Anqi ◽  
Mohammad Reza Safaei

Purpose The purpose of this study is to solve convective diffusion equation analytically by considering appropriate boundary conditions and using the Taylor-Aris method to determine the solute concentration, the effective and relative axial diffusivities. Design/methodology/approach >An analysis has been conducted on how body acceleration affects the dispersion of a solute in blood flow, which is known as a Bingham fluid, within an artery. To solve the system of differential equations analytically while validating the target boundary conditions, the blood velocity is obtained. Findings The blood velocity is impacted by the presence of body acceleration, as well as the yield stress associated with Casson fluid and as such, the process of dispersing the solute is distracted. It graphically illustrates how the blood velocity and the process of solute dispersion are affected by various factors, including the amplitude and lead angle of body acceleration, the yield stress, the gradient of pressure and the Peclet number. Originality/value It is witnessed that the blood velocity, the solute concentration and also the effective and relative axial diffusivities experience a drop when either of the amplitude, lead angle or the yield stress rises.


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.


2015 ◽  
Vol 137 (5) ◽  
Author(s):  
Şahin Yİğİt ◽  
Robert J. Poole ◽  
Nilanjan Chakraborty

The effects of inclination 180deg≥φ≥0deg on steady-state laminar natural convection of yield-stress fluids, modeled assuming a Bingham approach, have been numerically analyzed for nominal values of Rayleigh number Ra ranging from 103 to 105 in a square enclosure of infinite span lying horizontally at φ=0deg, then rotated about its axis for φ>0deg cases. It has been found that the mean Nusselt number Nu¯ increases with increasing values of Rayleigh number but Nu¯ values for yield-stress fluids are smaller than that obtained in the case of Newtonian fluids with the same nominal value of Rayleigh number Ra due to the weakening of convective transport. For large values of Bingham number Bn (i.e., nondimensional yield stress), the mean Nusselt number Nu¯ value settles to unity (Nu¯=1.0) as heat transfer takes place principally due to thermal conduction. The mean Nusselt number Nu¯ for both Newtonian and Bingham fluids decreases with increasing φ until reaching a local minimum at an angle φ* before rising with increasing φ until φ=90deg. For φ>90deg the mean Nusselt number Nu¯ decreases with increasing φ before assuming Nu¯=1.0 at φ=180deg for all values of Ra. The Bingham number above which Nu¯ becomes unity (denoted Bnmax) has been found to decrease with increasing φ until a local minimum is obtained at an angle φ* before rising with increasing φ until φ=90deg. However, Bnmax decreases monotonically with increasing φ for 90deg<φ<180deg. A correlation has been proposed in terms of φ, Ra, and Bn, which has been shown to satisfactorily capture Nu¯ obtained from simulation data for the range of Ra and φ considered here.


Author(s):  
Mohd. Ashique Hassan ◽  
Manabendra Pathak ◽  
Mohd. Kaleem Khan

In this study a computational investigation of two-dimensional, steady-state, natural convection of viscoplastic fluid in a square enclosure has been presented. The enclosure has been locally heated from the bottom wall using a constant heat flux source and symmetrically cooled from both the side walls. The other walls are maintained as insulated surfaces. Finite volume based code has been used in the simulation and Bingham model has been used to model the rheology of the enclosed viscoplastic fluids. Simulations have been made for three different heating lengths of the bottom wall. The flow phenomenon and heat transfer inside the enclosure have been investigated for different properties of viscoplastic fluid, heating conditions and heated length. It has been observed that for a particular thermal condition the heat transfer coefficient or the Nusselt number decrease with the increase in yield stress value of the fluid due to weakening of convective circulation.


Sign in / Sign up

Export Citation Format

Share Document