Numerical investigation for natural convection of a nanofluid in an inclined L-shaped cavity in the presence of an inclined magnetic field

The problem of natural convection in an inclined L-shaped enclosure filled with Cu/water nanofluid that operates under differentially heated walls in the presence of an inclined magnetic field is presented in this paper. The fully implicit finite difference method is used to solve the governing equations. A comparison with previously published results in special case of the present study is performed and a very good agreement is found. Heat transfer and fluid flow are examined for parameters of the Hartmann number (0 ≤ Ha ≤ 100), the nanoparticles volume fraction (0% ≤ ϕ ≤ 20%), the cavity inclination angle (0° ≤ ϑ ≤ 300°), the magnetic field inclination angle (0° ≤ γ ≤ 270°), the cavity aspect ratio (0.25 ≤ AR ≤ 0.6) and the Rayleigh number (10³ ≤ Ra ≤ 10⁶). It is found that, the presence of the magnetic field in the fluid region causes a significant reduction in the fluid flow and heat transfer characteristics. Also, a good enhancement in the heat transfer rate can be obtained by adding the copper nanoparticles to the base fluid.     

Natural convection of water-based nanofluids in an inclined enclosure with a heat source

This study investigates natural convection heat transfer of water-based nanofluids in an inclined square enclosure where the left vertical side is heated with a constant heat flux, the right side is cooled, and the other sides are kept adiabatic. The governing equations are solved using polynomial differential quadrature (PDQ) method. Calculations were performed for inclination angles from 0° to 90°, solid volume fractions ranging from 0% to 20%, constant heat flux heaters of lengths 0.25, 0.50 and 1.0, and a Rayleigh number varying from 10⁴ to 10⁶. The ratio of the nanolayer thickness to the original particle radius is kept at a constant value of 0.1. The heat source is placed at the center of the left wall. Five types of nanoparticles are taken into consideration: Cu, Ag, CuO, Al₂O₃, and TiO₂. The results show that the average heat transfer rate increases significantly as particle volume fraction and Rayleigh number increase. The results also show that the length of the heater is also an important parameter affecting the flow and temperature fields. The average heat transfer decreases with an increase in the length of the heater. As the heater length is increased, the average heat transfer rate starts to decrease for a smaller inclination angle (it starts to decrease with inclination at 90° for ? = 0.25, 60° for ? = 0.50, 45° for ? = 1.0, respectively).     

Natural convection in a square cavity filled with a porous medium: Effects of various thermal boundary conditions

Natural convection flows in a square cavity filled with a porous matrix has been studied numerically using penalty finite element method for uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls. Darcy–Forchheimer model is used to simulate the momentum transfer in the porous medium. The numerical procedure is adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.71 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous thermal boundary conditions. Numerical results are presented in terms of stream functions, temperature profiles and Nusselt numbers. Non-uniform heating of the bottom wall produces greater heat transfer rate at the center of the bottom wall than uniform heating case for all Rayleigh numbers but average Nusselt number shows overall lower heat transfer rate for non-uniform heating case. It has been found that the heat transfer is primarily due to conduction for Da ⩽ 10⁻⁵ irrespective of Ra and Pr. The conductive heat transfer regime as a function of Ra has also been reported for Da ⩾ 10⁻⁴. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes the power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Plume separation from an adiabatic horizontal thin fin placed at different heights on the sidewall of a differentially heated cavity

The separation of plumes from an adiabatic horizontal thin fin attached to a sidewall of a differentially heated cavity at quasi-steady stage is experimentally and numerically studied at three Rayleigh numbers (0.92 × 10⁹, 1.84 × 10⁹ and 3.68 × 10⁹) and over a range of fin positions. Regular plume separation is observed over the present range of parameters during the quasi-steady stage. Both experimental and numerical results reveal that the plume separation frequency increases with the Rayleigh number and decreases with the fin height measured from the leading edge. A higher Rayleigh number leads to a more unstable flow above the horizontal thin fin which in turn leads to a higher plume separation frequency. The decrease of the plume separation frequency with the increasing fin height is mainly due to the reduction of the adverse temperature gradient in the unstable layer above the thin fin as a result of the cavity-wide temperature stratification. It is further revealed that the heat transfer through the sidewall is improved by the presence of the thin fin. An optimum fin height for maximum heat transfer enhancement has been identified for the case with a Rayleigh number of 0.92 × 10⁹. For the other two higher Rayleigh number cases considered in this study, the heat transfer through the sidewall monotonically decreases with the fin height.     

Effects of duct inclination angle on thermal entrance region of laminar and transition mixed convection

Experimental investigation was performed on the mixed convection heat transfer of thermal entrance region in an inclined rectangular duct for laminar and transition flow. Air flowed upwardly and downwardly with inclination angles from ?90° to 90°. The duct was made of duralumin plate and heated with uniform heat flux axially. The experiment was designed for determining the effects of inclination angles on the heat transfer coefficients and friction factors at seven orientations (θ = ? 90°, ?60°, ?30°, 0°, 30°, 60° and 90°), six Reynolds numbers (Re ≈ 420, 840, 1290, 1720, 2190 and 2630) within the range of Grashof numbers from 6.8 × 10³ to 4.1 × 10⁴. The optimum inclination angles that yielded the maximum heat transfer coefficients decreased from 30° to ?30° with the increase of Reynolds numbers from 420 to 1720. The heat transfer coefficients first increased with inclination angles up to a maximum value and then decreased. With further increase in Reynolds numbers, the heat transfer coefficients were nearly independent of inclination angles. The friction factors decreased with the increase of inclination angles from ?90° to 90° when Reynolds numbers ranged from 420 to 1290, and independent of inclination angles with higher Reynolds numbers.     

Visualization of natural convection heat transport using heatline method in porous non-isothermally heated triangular cavity

Natural convection heat transfer in a porous media filled and non-isothermally heated from the bottom wall of triangular enclosure is analyzed using finite difference technique. Darcy law was used to write equations of porous media. Dimensionless heatfunction was used to visualize the heat transport due to buoyancy forces. Three different boundary conditions were applied for the vertical and inclined boundaries of triangular enclosures as Case I; both vertical and inclined walls were isothermal, Case II; vertical wall was adiabatic and inclined one was isothermal, Case III; vertical wall is isothermal and inclined one is adiabatic. A cosine function was utilized to get non-isothermal wall condition. The study was performed for different aspect ratios (0.25 ? AR ? 1.0) and Darcy-modified Rayleigh numbers (100 ? Ra ? 1000). It was observed that heat transfer enhancement was formed when vertical and inclined walls were isothermal while bottom wall was at non-uniform temperature. Heat transfer from bottom wall did not vary when the value of aspect ratio was higher than 0.50. In addition, heatline visualization technique was a useful technique for non-isothermally heated and porous media filled triangular enclosures.     

Effects of thermal boundary conditions on natural convection flows within a square cavity

A numerical study to investigate the steady laminar natural convection flow in a square cavity with uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls has been performed. A penalty finite element method with bi-quadratic rectangular elements has been used to solve the governing mass, momentum and energy equations. The numerical procedure adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁵ and Prandtl number Pr, 0.7 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous Dirichlet boundary conditions. Non-uniform heating of the bottom wall produces greater heat transfer rates at the center of the bottom wall than the uniform heating case for all Rayleigh numbers; however, average Nusselt numbers show overall lower heat transfer rates for the non-uniform heating case. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes, power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Lattice Boltzmann simulation of natural convection heat transfer in an elliptical-triangular annulus

A numerical study for steady-state, laminar natural convection in a horizontal annulus between a heated triangular inner cylinder and cold elliptical outer cylinder was investigated using lattice Boltzmann method. Both inner and outer surfaces are maintained at the constant temperature and air is the working fluid. Study is carried out for Rayleigh numbers ranging from 1.0 × 10³ to 5.0 × 10⁵. The effects of different aspect ratios and elliptical cylinder orientation were studied at different Rayleigh numbers. The local and average Nusselt numbers and percentage of increment heat transfer rate were presented. The average Nusselt number was correlated. The results show that by decreasing the value of aspect ratio and/or increasing the Rayleigh number, the Nusselt number increases. Also the heat transfer rate increases when the ellipse positioned vertically.     

Heat transfer by free convection from the inside surface of the vertical and inclined elliptic tube

Free convection from the inside surface of vertical and inclined elliptic tubes of axis ratio (a:b) 2:1 with a uniformly heated surface (constant heat flux) is investigated experimentally. The effects of orientation angle (α) and inclination angle (ϕ) on the heat transfer coefficient were studied. The orientation angle (α) is varied from 0° (when the major axis is horizontal) to 90° (when the major axis is vertical) with steps of 15°. The inclination angle (ϕ) is measured from the horizontal and varied from 15° to 75° with steps of 15°. The vertical position is considered as a special case of the inclined case when ϕ = 90. The experiments covered a range of Rayleigh number, Ra from 2.6 × 10⁶ to 3.6 × 10⁷. The local and average Nusselt numbers are estimated for different orientation angles and inclination angles at different Rayleigh numbers. The results obtained showed that the local Nu increased with the increase of axial distance from the lower end of the elliptic tube until a maximum value near the upper end, and then, it gradually decreased. The average Nu increases with the increase of α or ϕ at the same Ra. The results obtained are correlated by dimensionless groups and with the available data of the inclined and vertical elliptic tubes.     

Experimental investigation of mixed convection heat transfer from longitudinal fins in a horizontal rectangular channel

Mixed convection heat transfer from longitudinal fins inside a horizontal channel has been investigated for a wide range of modified Rayleigh numbers and different fin heights and spacings. An experimental parametric study was made to investigate effects of fin spacing, fin height and magnitude of heat flux on mixed convection heat transfer from rectangular fin arrays heated from below in a horizontal channel. The optimum fin spacing to obtain maximum heat transfer has also been investigated. During the experiments constant heat flux boundary condition was realized and air was used as the working fluid. The velocity of fluid entering channel was kept nearly constant (0.15 ? w_in ? 0.16 m/s) using a flow rate control valve so that Reynolds number was always about Re = 1500. Experiments were conducted for modified Rayleigh numbers 3 × 10⁷ < Ra¹ < 8 × 10⁸ and Richardson number 0.4 < Ri < 5. Dimensionless fin spacing was varied from S/H = 0.04 to S/H = 0.018 and fin height was varied from H_f/H = 0.25 to H_f/H = 0.80. For mixed convection heat transfer, the results obtained from experimental study show that the optimum fin spacing which yields the maximum heat transfer is S = 8–9 mm and optimum fin spacing depends on the value of Ra¹.     

Experimental study of evaporation heat transfer of R-134a inside a corrugated tube with different tube inclinations

Experimental heat transfer studies during evaporation of R-134a inside a corrugated tube have been carried out. The corrugated tube has been provided with different tube inclination angles of the direction of fluid flow from horizontal, α. The experiments were performed for seven different tube inclinations, α, in a range of − 90° to + 90° and four mass velocities of 46, 81, 110 and 136 kg m^− 2 s^− 1 for each tube inclination angle during evaporation of R-134a. Data analysis demonstrate that the tube inclination angle, α, affects the boiling heat transfer coefficient in a significant manner. The effect of tube inclination angle, α, on heat transfer coefficient, h, is more prominent at low vapor quality and mass velocity. In the low vapor quality region, the heat transfer coefficient, h, for the + 90° inclined tube is about 62% more than that of the − 90° inclined tube. The results also showed that at all mass velocities, the highest average heat transfer coefficient were achieved for α = + 90°. An empirical correlation has also been developed to predict the heat transfer coefficient during flow boiling inside a corrugated tube with different tube inclinations.     

Laminar natural convection in a square cavity: Low Prandtl numbers and large density differences

Steady natural convection at low Prandtl numbers caused by large density differences in a square cavity heated through the side walls is investigated numerically and theoretically. An appropriate dimensionless parameter characterizing the density differences of the working fluid is identified by the Gay-Lussac number. The Boussinesq assumption is achieved when the Gay-Lussac number tends to zero. The Nusselt number is derived for the ranges in Rayleigh number 10 ? Ra ? 10⁸, in Prandtl number 0.0071 ? Pr ? 7.1 and in Gay-Lussac number 0 ? Ga < 2. The effects of the Rayleigh, Prandtl and Gay-Lussac numbers on the Nusselt number are discussed on physical grounds by means of a scale analysis. Finally, based on physical arguments, a heat transfer correlation is proposed, valid for all Prandtl and Gay-Lussac number ranges addressed.     

Cooling of a sphere by natural convection – The applicability of the lumped capacitance method

Two approaches to predicting the sphere cooling process by laminar natural convection were compared in terms of the accuracy of the volume averaged sphere temperature and the heat transfer rate between the sphere and the surrounding fluid. The first approach is based on the formulation of conjugate heat transfer (heat conduction in the sphere and laminar natural convection in the fluid). The second approach includes the lumped capacity method based on the assumption that the temperature in the sphere is spatially uniform and on the Churchill correlation function. The solution to the problem depends on the Rayleigh number (Ra), the Biot number (Bi), the Prandtl number (Pr), and the sphere-to-fluid thermal diffusivity ratio (A). The lumped capacitance method gives fairly accurate results with respect to the conjugate heat transfer method (discrepancy in the volume averaged sphere temperature less than 5%) when A · Bi/Ra^0.452 < 0.05, for Bi < 0.15 and Pr > 1.     

Laminar natural convection of non-Newtonian power-law fluids between concentric circular cylinders

The two-dimensional steady-state natural convection of power-law fluids is studied numerically between two concentric horizontal cylinders with different constant temperatures. The governing equations are discretized using finite volume technique based on second order upwind and are solved using the SIMPLE algorithm. The effects of Rayleigh number (10³ ≤ Ra ≤ 10⁵) and Prandtl number (10 ≤ Pr ≤ 10³) on the dimensionless velocity and temperature are investigated for both pseudoplastic and dilatant fluids. Also the mean Nusselt number for various values of governing parameters is obtained and discussed. The results indicate that with increasing the power-law index from 0.6 to 1.4, the mean Nusselt number decreases. In the best case among the range of parameters considered here the heat transfer rate for pseudo-plastic fluid (n = 0.6) is 170% higher than the Newtonian one and for dilatant fluid (n = 1.4) the heat transfer rate is 43% lower than the Newtonian fluid. So the pseudoplastic and dilatant fluids are more efficient than Newtonian fluids for cooling and insulating purposes, respectively. It is shown that as the Rayleigh number increases the cooling effect of pseudoplastic fluid and the insulating effect of dilatant fluid become more pronounced.     

Simulation of heat transfer enhancement in nanofluids using dissipative particle dynamics

The current paper applied dissipative particle dynamics (DPD) approach to investigate heat transfer within nanofluids. The DPD approach was applied to study natural convection in a differential heated enclosure by considering the viscosity and the thermal conductivity of the nanofluid to be dual function of temperature and volume fraction of nanoparticles. Experimental data for viscosity and thermal conductivity are incorporated in the current DPD model to mimic energy transport within nanofluids. This incorporation is done through the modification of the dissipative weighting function that appears in the dissipative force vector and the dissipative heat flux. For the entire range of Rayleigh number considered in this study, it was found that the DPD results show a deterioration in heat transfer in the enclosure due to the presence of nanoparticles for φ > 4%. However, some slight enhancement is shown to take place for small volume fraction of nanoparticles, φ ≤ 4%. The DPD results experienced some degree of compressibility at high values of Rayleigh number Ra ≥ 10⁵.     

Investigation of natural circulation in cavities with uniform heat generation for different Prandtl number fluids

Natural convection in enclosures with uniform heat generation and isothermal side walls is studied here. For the rectangular enclosure, two-dimensional conservation equations are solved using SIMPLE algorithm. Parametric studies are conducted to examine the effects of orientation of the cavity, fluid properties (Pr number), and aspect ratio for Rayleigh numbers up to 10⁶. For a horizontal square cavity, the flow becomes periodically oscillating at Ra = 5 × 10⁴ and chaotic at Ra = 8 × 10⁵. With a slight increase in the inclination angle, the oscillations die and for inclination angles greater than 15⁰, the flow attain a steady state over a range of Ra. It is found that for tall cavities (aspect ratio > 1), the steady-state solution is obtained for all values of Ra considered here. However, for wide cavities (aspect ratio < 1), an oscillatory flow regime is observed. The maximum temperature within the cavity is calculated for the range of Ra, aspect ratio and Pr number. Correlations for the maximum cavity temperature is presented here. The values of critical Rayleigh number at which the convection sets in the rectangular cavity are also studied and two distinct criteria are determined to evaluate the critical Rayleigh number. Further, a three-dimensional simulation is performed for a cubic cavity. It is found that the steady state solutions are obtained for all Rayleigh number, except at Ra = 10⁶. This is in contrast to the predictions for a two-dimensional square cavity, which has an oscillatory zone from Ra = 5 × 10⁴ onwards.     

Natural convection and fluid flow in inclined enclosure with a corner heater

The phenomena of natural convection in an inclined square enclosure heated via corner heater have been studied numerically. Finite difference method is used for solving momentum and energy equations in the form of stream function–vorticity. One wall of the enclosure is isothermal but its temperature is colder than that of heaters while the remaining walls are adiabatic. The numerical procedure adopted in this analysis yields consistent performance over a wide range of parameters; Rayleigh number, Ra (10³ ? Ra ? 10⁶); Prandtl number, Pr (0.07 ? Pr ? 70); dimensionless lengths of heater in x and y directions (0.25 ? hx ? 0.75, 0.25 ? hy ? 0.75); and inclination angle, ? (0° ? ? ? 270°). It is observed that heat transfer is maximum or minimum depending on the inclination angle and depending on the length of the corner heaters. The effect of Prandtl number on mean Nusselt number is more significant for Pr < 1.     

Numerical studies on laminar natural convection inside inclined cylinders of unity aspect ratio

The effect of cylinder inclination on thermal buoyancy induced flows and internal natural convective heat transfer is explored using CFD simulations. The cylinder’s top and bottom surfaces were maintained at different temperatures while the curved surface was adiabatic. The aspect ratio (length/diameter) of the cylinder was unity and the Prandtl number of the fluid was fixed at 0.71. The Rayleigh number of the confined fluid was varied from 10³ to 3.1 × 10⁴ by changing the specified end wall temperatures. The critical Rayleigh number was estimated to be 3800 for the vertical cylinder. Relaxing the convergence criterion caused false hysteresis in the converged results for the vertical cylinder. Typical natural convective fluid flow and temperature patterns obtained under laminar flow conditions are illustrated for various inclinations ranging from 0° to 180°. Flow visualization studies revealed complex three-dimensional patterns. Different thermal–hydrodynamic regimes were identified and were classified in terms of Rayleigh number and angle of inclination. Empirical correlations for the Nusselt number and maximum velocities in the domain as a function of the inclination angle and Rayleigh number are developed.     

Experimental investigation of mixed convection heat transfer for thermally developing flow in a horizontal circular cylinder

An experimental investigation is presented on mixed (free and forced) convection to study the local and average heat transfer for hydrodynamically fully developed, thermally developing and thermally fully developed laminar air flow in a horizontal circular cylinder. The experimental setup consists of aluminum cylinder as test section with 30 mm inside diameter and 900 mm heated length (L/D = 30), is subjected to a constant wall heat flux boundary condition. The investigation covers Reynolds number range from 400 to 1600, the heat flux varied from 60 W/m² to 400 W/m² and with cylinder inclination angle of θ = 0° (horizontal). The hydrodynamically fully developed condition is achieved by using an aluminum entrance section pipes (calming sections) having the same inside diameter as test section pipe but with variable lengths. The entrance sections included two long calming sections, one with length of 180 cm (L/D = 60), another one with length of 240 cm (L/D = 80) and two short calming sections with lengths 60 cm (L/D = 20), 120 cm (L/D = 40). The surface temperature variation along the cylinder surface, the local and average Nusselt number variation with the dimensionless axial distance Z⁺ were presented. For all entrance sections, it was found an increase in the Nusselt number values as the heat flux increases. It was concluded that the free convection effects tended to decrease the heat transfer results at low Re while to increase the heat transfer results for high Re. The combined convection regime could be bounded by a suitable selection of Re number ranges and the heat flux ranges. The obtained Richardson numbers (Ri) range varied approximately from 0.13 to 7.125. The average Nusselt numbers were correlated with the (Rayleigh numbers/Reynolds numbers). The proposed correlation has been compared with available literature and showed satisfactory agreement.     

Numerical study of assisting and opposing mixed convective nanofluid flows in an inclined circular pipe

A numerical investigation of mixed convection is carried out to study the heat transfer and fluid flow characteristics in an inclined circular pipe using the finite volume method. The pipe has L/D of 500 and it was subjected to a uniform heat flux boundary condition. Four types of nanofluids (Al₂O₃, CuO, SiO₂, and TiO₂ with H₂O) with nanoparticles concentration in the range of 0 ≤ φ ≤ 5% and nanoparticles diameter in the range of 20 ≤ dp ≤ 60 nm were used. The pipe inclination angle was in the range of 30^∘ ≤ θ ≤ 75^∘ using assisting and opposing flow. The influences of Reynolds number in the range of 100 ≤ Re ≤ 2000, and Grashof numbers in the range of 6.3 × 10² ≤ Gr ≤ 8.37 × 10³ were examined. It is found that the velocity and wall shear stress are increased as Re number increases, while the surface temperature decreases. There is no significant effect of increasing Gr number on thermal and flow fields. The velocity and wall shear stress are increased and the surface temperature is decreased as φ and dp are decreased. It is concluded that the surface temperature is increased as the pipe inclination angle increases from the horizontal position (θ = 0^°) to the inclined position (θ = 75^°). In addition, it is inferred that the heat transfer is enhanced using SiO₂ nanofluid compared with other nanofluids types. Furtheremore, it is enhanced using assisting flow compared to opposing flow.