MHD natural convection flow from an isothermal horizontal circular cylinder under consideration of temperature dependent viscosity

Magnetohydrodynamic (MHD) natural convection laminar flow from an iso-thermal horizontal circular cylinder immersed in a fluid with viscosity proportional to a linear function of temperature will be discussed with numerical simulations. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equa-tions are reduced to convenient form, which are solved numerically by two very efficient methods, namely, (i) Implicit finite difference method together with Keller box scheme and (ii) Direct numerical scheme. Numerical results are presented by velocity and temperature distributions of the fluid as well as heat transfer characteristics, namely the shearing stress and the local heat transfer rate in terms of the local skin-friction coefficient and the local Nusselt number for a wide range of magnetohydrodynamic parameter, viscosity-variation parameter and viscous dissipation parameter. MHD flow in this geometry with temperature dependent viscosity is absent in the literature. The results obtained from the numerical simulations have been veri-fied by two methodologies.

Scaling for the Prandtl number of the natural convection boundary layer of an inclined flat plate under uniform surface heat flux

An improved scaling analysis and direct numerical simulations are performed for the unsteady natural convection boundary layer adjacent to a downward facing inclined plate with uniform heat flux. The development of the thermal or viscous boundary layers may be classified into three distinct stages: a start-up stage, a transitional stage and a steady stage, which can be clearly identified in the analytical as well as the numerical results. Previous scaling shows that the existing scaling laws of the boundary layer thickness, velocity and steady state time scale for the natural convection flow on a heated plate of uniform heat flux provide a very poor prediction of the Prandtl number dependency of the flow. However, those scalings perform very well with Rayleigh number and aspect ratio dependency. In this study, a modified Prandtl number scaling is developed using a triple layer integral approach for Pr > 1. It is seen that in comparison to the direct numerical simulations, the modified scaling performs considerably better than the previous scaling.

Natural convection in a triangular enclosure due to non-uniform cooling on top

Natural convection in a triangular enclosure subject to non-uniformly cooling at the inclined surfaces and uniformly heating at the base is investigated numerically. The numerical simulations of the unsteady flows over a range of Rayleigh numbers and aspect ratios are carried out using Finite Volume Method. Since the upper surface is cooled and the bottom surface is heated, the air flow in the enclosure is potentially unstable to Rayleigh Benard instability. It is revealed that the transient flow development in the enclosure can be classified into three distinct stages; an early stage, a transitional stage and a steady stage. It is also found that the flow inside the enclosure strongly depends on the governing parameters; Rayleigh number and aspect ratio. The asymmetric behaviour of the flow about the geometric centre line is discussed in detailed. The heat transfer through the roof and the ceiling as a form of Nusselt number is also reported in this study.

Natural convection in the attic-shaped spaces

Heat transfer through an attic space into or out of buildings is an important issue for attic-shaped houses in both hot and cold climates. One of the important objectives for design and construction of houses is to provide thermal comfort for occupants. In the present energy-conscious society, it is also a requirement for houses to be energy efficient, i.e. the energy consumption for heating or air-conditioning houses must be minimized. Relevant to these objectives, research into heat transfer in attics has been conducted for about three decades. The transient behaviour of an attic space is directly relevant to our daily life. Certain periods of the day or night may be considered as having a constant ambient temperature (e.g. during 11am - 2pm or 11pm - 2am). However, at other times during the day or night the ambient temperature changes with time (e.g. between 5am - 9am or 5pm - 9pm). Therefore, the analysis of steady state solution is not sufficient to describe the fluid flow and heat transfer in the attic space. The discussion of the transient development of the boundary is required. A theoretical understanding of the transient behaviour of the flow in the enclosure is performed through scaling analysis for sudden and ramp heating conditions. A proper identification of the timescales, the velocity and the thickness relevant to the flow that develops inside the cavity makes it possible to predict theoretically the basic flow features that will survive once the thermal flow in the enclosure reaches a steady state. Those scaling predictions have been verified by a series of numerical simulations.

Prandtl number scaling of the unsteady natural convection boundary layer adjacent to a vertical flat plate for Pr > 1 subject to ramp surface heat flux

It is found in the literature that the existing scaling results for the boundary layer thickness, velocity and steady state time for the natural convection flow over an evenly heated plate provide a very poor prediction of the Prandtl number dependency of the flow. However, those scalings provide a good prediction of two other governing parameters’ dependency, the Rayleigh number and the aspect ratio. Therefore, an improved scaling analysis using a triple-layer integral approach and direct numerical simulations have been performed for the natural convection boundary layer along a semi-infinite flat plate with uniform surface heat flux. This heat flux is a ramp function of time, where the temperature gradient on the surface increases with time up to some specific time and then remains constant. The growth of the boundary layer strongly depends on the ramp time. If the ramp time is sufficiently long, the boundary layer reaches a quasi steady mode before the growth of the temperature gradient is completed. In this mode, the thermal boundary layer at first grows in thickness and then contracts with increasing time. However, if the ramp time is sufficiently short, the boundary layer develops differently, but after the wall temperature gradient growth is completed, the boundary layer develops as though the startup had been instantaneous.

Prandtl number scaling of natural convection of the flow on an inclined at plate due to uniform surface heat flux

A new scaling analysis has been performed for the unsteady natural convection boundary layer under a downward facing inclined plate with uniform heat flux. The development of the thermal or viscous boundary layers may be classified into three distinct stages including an early stage, a transitional stage and a steady stage, which can be clearly identified in the analytical as well as numerical results. Earlier scaling shows that the existing scaling laws of the boundary layer thickness, velocity and steady state time scales for the natural convection flow on a heated plate of uniform heat flux provide a very poor prediction of the Prandtl number dependency. However, those scalings performed very well with Rayleigh number and aspect ratio dependency. In this study, a modifed Prandtl number scaling has been developed using a triple-layer integral approach for Pr > 1. It is seen that in comparison to the direct numerical simulations, the new scaling performs considerably better than the previous scaling.

Unsteady natural convection within a differentially heated enclosure of sinusoidal corrugated side walls

Numerically investigation of natural convection within a differentially heated modified square enclosure with sinusoidally corrugated side walls has been performed for different values of Rayleigh number. The fluid inside the enclosure considered is air and is quiescent, initially. The top and bottom surfaces are flat and considered as adiabatic. Results reveal three main stages: an initial stage, a transitory or oscillatory stage and a steady stage for the development of natural convection flow inside the corrugated cavity. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. Investigation has been performed for the Rayleigh number, Ra ranging from 105 to 108 with variation of corrugation amplitude and frequency. Constant physical properties for the fluid medium have been assumed. Results have been presented in terms of the isotherms, streamlines, temperature plots, average Nusselt numbers, traveling waves and thermal boundary layer thickness plots, temperature and velocity profiles. The effects of sudden differential heating and its consequent transient behavior on fluid flow and heat transfer characteristics have been observed for the range of governing parameters. The present results show that the transient phenomena are greatly influenced by the variation of the Rayleigh Number with corrugation amplitude and frequency.

Natural convection of nanofluids in open rectangular cavity

An investigation of the effect of nano particles on natural convection of water based nanofluids contained in an open rectangular cavity is carried out numerically. The flow pattern and heat transfer characteristics are studied for different values of volume fraction in the range 0   0.2 , Rayleigh number in the range 9 1 Ra 10 and the nano particles with different thermo physical properties. It was found that for low Rayleigh numbers, heat transfer exhibits a decreasing trend for increasing values of volume fraction of oxide nanofluids, whereas for higher values of Rayleigh numbers, an increasing trend of heat transfer was observed due to increase in the volume fraction of nanofluids.

Finite volume approach of natural convection in a triangular enclosure with localized heating from below

In this study, natural convection heat transfer and buoyancy driven flows have been investigated in a right angled triangular enclosure. The heater located on the bottom wall while the inclined wall is colder and the remaining walls are maintained as adiabatic. Governing equations of natural convection are solved through the finite volume approach, in which buoyancy is modeled via the Boussinesq approximation. Effects of different parameters such as Rayleigh number, aspect ratio, prantdl number and heater location are considered. Results show that heat transfer increases when the heater is moved toward the right corner of the enclosure. It is also revealed that increasing the Rayleigh number, increases the strength of free convection regime and consequently increases the value of heat transfer rate. Moreover, larger aspect ratio enclosure has larger Nusselt number value. In order to have better insight, streamline and isotherms are shown.

Natural convection flow in a strong cross magnetic field with radiation

The problem of MHD natural convection boundary layer flow of an electrically conducting and optically dense gray viscous fluid along a heated vertical plate is analyzed in the presence of strong cross magnetic field with radiative heat transfer. In the analysis radiative heat flux is considered by adopting optically thick radiation limit. Attempt is made to obtain the solutions valid for liquid metals by taking Pr≪1. Boundary layer equations are transformed in to a convenient dimensionless form by using stream function formulation (SFF) and primitive variable formulation (PVF). Non-similar equations obtained from SFF are then simulated by implicit finite difference (Keller-box) method whereas parabolic partial differential equations obtained from PVF are integrated numerically by hiring direct finite difference method over the entire range of local Hartmann parameter, $xi$ . Further, asymptotic solutions are also obtained for large and small values of local Hartmann parameter $xi$ . A favorable agreement is found between the results for small, large and all values of $xi$ . Numerical results are also demonstrated graphically by showing the effect of various physical parameters on shear stress, rate of heat transfer, velocity and temperature.

Natural convection flow with surface radiation along a vertical wavy surface

In this study, natural convection boundary layer flow of thermally radiating fluid along a heated vertical wavy surface is analyzed. Here, the radiative component of heat flux emulates the surface temperature. Governing equations are reduced to dimensionless form, subject to the appropriate transformation. Resulting dimensionless equations are transformed to a set of parabolic partial differential equations by using primitive variable formulation, which are then integrated numerically via iterative finite difference scheme. Emphasis has been given to low Prandtl number fluid. The numerical results obtained for the physical parameters, such as, surface radiation parameter, R, and radiative length parameter, ξ, are discussed in terms of local skin friction and Nusselt number coefficients. Comprehensive interpretation of velocity distribution is also given in the form of streamlines.

Conduction-radiation effect on natural convection flow in a fluid-saturated non-Darcy porous medium enclosed by non-isothermal walls

The effect of conduction-convection-radiation on natural convection flow of Newtonian optically thick gray fluid, confined in a non-Darcian porous media square cavity is numerically studied. For the gray fluid consideration is given to Rosseland diffusion approximation. Further assuming that (i) the temperature of the left vertical wall is varying linearly with height, (ii) cooled right vertical and top walls and (iii) the bottom wall is uniformly-heated. The governing equations are solved using the Alternate Direct Implicit method together with the Successive Over Relaxation technique. The investigation of the effect of governing parameters namely the Forschheimer resistance (Γ), the Planck constant (Rd), and the temperature difference (Δ), on flow pattern and heat transfer characteristics has been carried out. It was seen that the reduction of flow and heat transfer occurs as the Forschheimer resistance is increased. On the other hand both the strength of flow and heat transfer increases as the temperature ratio, Δ, is increased.

Non-Newtonian natural convection flow along an isothermal horizontal circular cylinder using modified power-law model

Laminar two-dimensional natural convection boundary-layer flow of non-Newtonian fluids along an isothermal horizontal circular cylinder has been studied using a modified power-law viscosity model. In this model, there are no unrealistic limits of zero or infinite viscosity. Therefore, the boundary-layer equations can be solved numerically by using marching order implicit finite difference method with double sweep technique. Numerical results are presented for the case of shear-thinning as well as shear thickening fluids in terms of the fluid velocity and temperature distributions, shear stresses and rate of heat transfer in terms of the local skin-friction and local Nusselt number respectively.

Steady natural convection of non-Newtonian power-law fluid in a trapezoidal enclosure

Numerically investigation of free convection heat transfer in a differentially heated trapezoidal cavity filled with non-Newtonian Power-law fluid has been performed in this study. The left inclined surface is uniformly heated whereas the right inclined surface is maintained as uniformly cooled. The top and bottom surfaces are kept adiabatic with initially quiescent fluid inside the enclosure. Finite volume based commercial software FLUENT 14.5 is used to solve the governing equations. Dependency of various flow parameters of fluid flow and heat transfer is analyzed including Rayleigh number, Ra ranging from 10^5 to 10^7, Prandtl number, Pr of 100 to 10,000 and power index, n of 0.6 to 1.4. Outcomes have been reported in terms of isotherms, streamline, and local Nusselt number for various Ra, Pr, n and inclined angles. Grid sensitivity analysis is performed and numerically obtained results have been compared with those results available in the literature and found good agreement.

Effect of an adiabatic fin on natural convection heat transfer in a triangular enclosure

Natural convection thermal boundary layer adjacent to the heated inclined wall of a right angled triangle with an adiabatic fin attached to that surface is investigated by numerical simulations. The finite volume based unsteady numerical model is adopted for the simulation. It is revealed from the numerical results that the development of the boundary layer along the inclined surface is characterized by three distinct stages, i.e. a start-up stage, a transitional stage and a steady stage. These three stages can be clearly identified from the numerical simulations. Moreover, in presence of adiabatic fin, the thermal boundary layer adjacent to the inclined wall breaks initially. However, it is reattached with the downstream boundary layer next to the fin. More attention has been given to the boundary layer development near the fin area.