A boundary layer analysis of heat transfer by free convection from permeable horizontal cylinders of elliptic cross-section in porous media using a thermal non-equilibrium model

This work uses a thermal non-equilibrium model to study the free convection boundary layer flow driven by temperature gradients near a permeable horizontal cylinder of elliptic cross-section with constant wall temperature in a fluid-saturated porous medium. A coordinate transformation is used to obtain the nonsimilar boundary layer equations. The transformed boundary layer equations are then solved by the cubic spline collocation method. Results for the local Nusselt numbers are presented as functions of the porosity scaled thermal conductivity ratio, the heat transfer coefficient between solid and fluid phases, the transpiration parameter, and the aspect ratio when the major axis of the elliptical cylinder is vertical (slender orientation) and horizontal (blunt orientation). An increase in the porosity scaled thermal conductivity ratio or the heat transfer coefficient between the solid and fluid phases increases the heat transfer rates. Moreover, the use of suction (positive transpiration parameter) tends to increase the heat transfer rates between the porous medium and the surface.     

The thermal modeling of a matrix heat exchanger using a porous medium and the thermal non-equilibrium model

Heat transfer and fluid flow characteristics through a porous medium were investigated using numerical simulations and experiment. For the numerical simulations two models were created: a two-dimensional numerical model and a Fluent™ computational fluid dynamics (CFD) porous media model. The experimental investigation consisted of a flow channel with a porous medium section that was heated from below by a heat source. The results of the numerical models were compared to the experimental data in order to determine the accuracy of the models. The numerical model was then modified to better simulate a matrix heat exchanger. This numerical model then generated temperature profiles that were used to calculate the heat transfer coefficient of the matrix heat exchanger and develop a correlation between the Nusselt number and the Reynolds number.     

Effects of thermal nonequilibrium and non-uniform temperature gradients on the onset of convection in a heterogeneous porous medium

The simultaneous effect of local thermal nonequilibrium (LTNE), vertical heterogeneity of permeability, and non-uniform basic temperature gradient on the criterion for the onset of Darcy-Benard convection is studied. The eigenvalue problem is solved numerically using the Galerkin method. The interaction of various types of permeability heterogeneity and non-uniform basic temperature gradient functions on the stability characteristics of the system is analyzed. It is observed that the linear variation (about the mean) of the permeability and the basic temperature gradient with depth has no added effect on the criterion for the onset of convection. However, the concurrent variation in heterogeneous permeability and non-uniform basic temperature gradient functions has more stabilizing effect on the system, while opposite is the trend when the effect of non-uniform basic temperature gradient alone is present.     

Effect of thermal modulation on the onset of convection in a rotating fluid layer

The stability of a rotating horizontal fluid layer heated from below is examined when, the walls of the layer are subjected to time-periodic temperature modulation. The linear stability analysis is used to study the effect of infinitesimal disturbances. A regular perturbation method based on small amplitude of applied temperature field is used to compute the critical values of Rayleigh number and wavenumber. The shift in critical Rayleigh number is calculated as a function of frequency of modulation, Taylor number and Prandtl number. It is established that the instability can be enhanced by the rotation at low frequency symmetric modulation and with moderate to high frequency lower wall temperature modulation, whereas the stability can be enhanced by the rotation in case of asymmetric modulation. The effect of Taylor number and Prandtl number on the stability of the system is also discussed. We found that by proper tuning of modulation frequency, Taylor number and Prandtl number it is possible to advance or delay the onset of convection.     

The onset of Brinkman ferroconvection using a thermal non-equilibrium model

The onset of ferroconvection in a horizontal layer of ferrofluid saturated Brinkman porous medium heated uniformly from below in the presence of a uniform vertical magnetic field is investigated when the solid and fluid phases of the porous medium are in local thermal non-equilibrium (LTNE) using linear stability theory. The modified Brinkman–Forchheimer-extended Darcy equation is employed to describe the flow in the porous medium and a two-field model for energy equation each representing the solid and fluid phases separately is used. It is established that the principle of exchange of stability is valid. The authenticity of LTNE model over LTE model and also the ferromagnetic effects on the stability of the system are discussed in detail. The system is found to be more stable when the magnetic forces alone are present. Asymptotic analysis for both small and large values of scaled inter-phase heat transfer coefficient H_t is presented and the results are found to be in good agreement with those obtained from the exact formula. The established results in the literature have been reproduced as particular cases from the present study.     

应用热非平衡模型模拟梯度磁场作用下多孔介质方腔内空气热磁对流

梯度磁场可用来控制多孔介质内空气的自然对流传热过程.利用局部热非平衡模型对圆形截流线圈水平放置时三维多孔介质方腔内的空气热磁对流进行了数值研究.控制方程基本变量采用控制容积法离散,求解采用SIMPLE算法.计算过程中Ra为103～105,磁场力数γ为0～150、Da为10-7～100.获得了空气热磁对流的流场和温度场....     

Effect of rotation on the onset of thermal convection in a sparsely packed porous layer using a thermal non-equilibrium model

Linear and nonlinear stability of a rotating fluid-saturated sparsely packed porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium. The extended Darcy–Brinkman model that includes the time derivative and Coriolis terms is employed as a momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for both stationary and oscillatory convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of rotation and thermal diffusion that causes the convection to set in through oscillatory mode rather than stationary. The rotation inhibits the onset of convection in both stationary and oscillatory mode. The Darcy number stabilizes the system towards the oscillatory mode, while it has dual effect on stationary convection. Besides, the effect of porosity modified conductivity ratio, Darcy–Prandtl number and the ratio of diffusivities on the stability of the system is investigated. The nonlinear theory is based on the truncated representation of Fourier series method. The effect of thermal non-equilibrium on heat transfer is brought out. The transient behavior of the Nusselt number is investigated by using the Runge–Kutta method. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

Onset of Darcy-Bénard convection using a thermal non-equilibrium model

In this paper we use a two-field model for the separate modelling of the solid and fluid phase temperature fields in a fluid-saturated porous medium, and, in particular, we consider how the onset criterion for convection in a horizontal layer is affected by the adoption of such a model. In general we find that both the critical Rayleigh number and wavenumber are modified by the presence of thermal non-equilibrium effects. It is shown that the well-known result of Lapwood [Proc. Cambridge Philos. Soc. 44 (1948) 508] which corresponds to local thermal equilibrium (LTE), is recovered when taking the thermal equilibrium limit of the non-equilibrium analysis.We also present asymptotic solutions for both small and large values of H the inter-phase heat transfer coefficient, H, and compare this with the numerical solutions. For intermediate values of H we find that the critical wavenumber is always larger than π, the critical value for the LTE case. In some cases this critical wavenumber may be very large compared with π.     

Entropy generation for forced convection in a porous saturated circular tube with uniform wall temperature

A numerical study is reported to investigate both the First and the Second Law of Thermodynamics for thermally developing forced convection in a circular tube filled by a saturated porous medium, with uniform wall temperature, and with the effects of viscous dissipation included. A theoretical analysis is also presented to study the problem for the asymptotic region applying the perturbation solution of the Brinkman momentum equation reported by (K. Hooman, K., A.A. Ranjbar-Kani, A perturbation based analysis to investigate forced convection in a porous saturated tube, Journal of Computational and Applied Mathematics 162 (2) (2004) 411–419.). Expressions are reported for the temperature profile, the Nusselt number, the Bejan number, and the dimensionless entropy generation rate in the asymptotic region. Numerical results are found to be in good agreement with theoretical counterparts.     

Stability of a fluid-saturated porous medium heated from below by forced convection

A linear stability analysis determining the onset of convection in a bounded rectangular cavity containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top wall, and bottom wall heated by forced convection. The nature of the bottom wall heating necessarily involves the Biot number, Bi. Numerical calculations of the critical Rayleigh number, R_c made over the range of Biot numbers 10⁻⁴?Bi?10⁴ for cavity aspect ratios 0?(a,b)?5 cover all effective bottom heating conditions from the constant heat flux global limit, R_c=27.096 found as Bi→0 to the isothermal global limit, R_c=4π² found as Bi→∞. Marginal stability boundaries, preferred cellular modes and disturbance temperature contours are displayed graphically.     

Analysis of free convection about a horizontal cylinder in a porous media using a thermal non-equilibrium model

The thermally non-equilibrium model is used to study the free convection from a horizontal cylinder immersed in porous media. The governing equations are transformed to dimensionless form by introducing the boundary layer dimensionless variables. The resultant parabolic system of differential equations is solved by using an implicit finite difference method based on Keller box algorithm. The results of the developed code are validated with different mesh sizes, which can be used as benchmark results for thermally equilibrium condition. Numerical results are obtained for non-equilibrium model to analyze the effect of the governing parameters, which are the heat transfer coefficient between the solid and fluid phases H and the porosity scaled thermal conductivity ratio K_r. The results show that increasing H or K_r leads to increase in the total average Nusselt number. The value of the average Nusselt number for both fluid and solid phases as well as the total average Nusselt number have approached the corresponding values of the thermally equilibrium model at high value of H × K_r.     

Thermosolutal convection in a fluid‐porous composite medium

In the present study analysis has been performed for thermosolutal convection in a fluid‐porous composite medium, consisting of a fluid‐saturated porous medium followed by an overlaying clean medium. The fluid‐porous composite medium is subjected to both a horizontal solutal and a thermal gradient. Top and bottom walls of the fluid‐porous composite medium are assumed to be impermeable and adiabatic. The Darcy‐Brinkman‐Forchheimer model is used to study the flow through the fluid‐porous composite medium. A single domain approach is taken into consideration for numerical simulation. The solution is done by control volume integration. A comprehensive analysis has been performed for various pertinent parameters to delineate their behavior. Results of the transport phenomenon have been provided in graphical and tabular form, for the complete understanding of the complex phenomenon. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library (wileyonlinelibrary.com/journal/htj). DOI 10.1002/htj.21048     

Combined radiation and convection heat transfer in a porous channel bounded by isothermal parallel plates

Combined radiation and convection heat transfer in a porous medium confined between gray isothermal parallel plates is investigated. The medium is absorbing, emitting and scattering. Cases of boundaries at temperatures higher or lower than the medium are considered. In the porous medium, the boundary effect on the fully developed laminar velocity field as proposed by Kaviany is accounted for. For various values of the extinction coefficient, the scattering albedo, the conduction-radiation parameter and the boundary emissivity, Nusselt number, temperature and heat flux distributions are found for the range of values including the extreme limits of the porous medium shape parameter (PMSP), γ=(W²φ/K)^1/2, where W is the channel width, φ the porosity and K the permeability. For the lower limiting value of the PMSP γ, the effect of the porous medium is negligible and the situation approaches that of Poiseuille flow. For this limiting case, results from the present work are compared with those available in the literature. For medium to high values of the PMSP γ, for the purpose of comparison, some results are presented in tabular form. Radiation is found to have a significant effect on various parameters studied. The discrete transfer method was used for the solution of the radiative part of the energy equation. An iterative finite difference scheme was used to solve the energy equation.     

Heat generation effects in natural convection inside a porous annulus

The interplay between internal heat generation and externally driven natural convection inside a porous medium annulus is studied in detail using numerical methods. The axisymmetric domain is bounded with adiabatic top and bottom walls and differentially heated side walls sustaining steady natural convection of a fluid with Prandtl number, Pr = 5, through a porous matrix of volumetric porosity, ? = 0.4. The generalized momentum equation with Brinkman–Darcy–Forchheimer terms and the local thermal non-equilibrium based two-energy equation model are solved to determine the flow and the temperature distribution. Beyond a critical heat generation value defined using an internal Rayleigh number, Ra_I,cr^?, the convection transits from unicellular to bicellular mode, as the annulus T_max becomes higher than the fixed hot-wall temperature. The Ra_I,cr^? increases proportionately when the permeability based external Rayleigh number Ra_E^? and the solid–fluid thermal conductivity ratio γ are independently increased. A correlation is proposed to predict the overall annulus Nu as a function of Ra_E^?, Ra_I^?, Da and γ. It predicts the results within ± 20% accuracy.     

Aspect ratio effect on natural convection in a square enclosure with a sinusoidal active thermal wall using a thermal non-equilibrium model

ABSTRACT

A numerical investigation of the aspect ratio effect on natural convection in a square enclosure is carried out by adopting the local thermal non-equilibrium model. The top and bottom walls of the enclosure are adiabatic, the left vertical wall is partially heated and cooled by the sinusoidal thermal boundary condition, and the right vertical wall is maintained at uniform thermal boundary condition. The results show the value of periodicity parameter increasing. The streamlines vary in different patterns, rotating clockwise and counterclockwise simultaneously when N > 1, and the number of clockwise and counterclockwise rotating cells increases with the increase of N and equals the value of N. The sinusoidal local Nusselt number profiles are observed and the wave amplitude of local Nusselt number decreases with the increase of aspect ratio, and the absolute values of average Nusselt number at left wall of porous cavity reach maximum when Ar = 1. The absolute value of solid-to-fluid temperature differences decreases as the inter-phase heat transfer coefficient (H) increases and it increases as the value of aspect ratio increases. The total heat transfer of porous cavity can be enhanced by increasing the aspect ratio and the thermal conductivity ratio.     

Linear stability of convection in a vertical channel filled with nanofluid saturated porous medium

The convection in a vertical channel filled with a porous medium saturated by a nanofluid is studied numerically. The effects of Brownian motion and thermophoresis are incorporated in the model used for nanofluid. Also, the flow within the porous region is governed by Brinkman's equation. The generalized eigenvalue problem for the perturbed state is obtained from a normal mode analysis and solved using the Chebyshev spectral collocation method. The Rayleigh number is expressed as an implicit function of the wavenumber with other parameters. The critical wavenumber and the critical Rayleigh number are calculated for different parameters. The preferred modes under critical conditions are detected.     

Analysis of mixed convection in a vertical porous layer using non-equilibrium model

The problem of two-dimensional steady mixed convection in a vertical porous layer is investigated numerically in the present paper using the thermally non-equilibrium model. The vertical porous layer is assumed to have a finite isothermally heated segment on one vertical wall which is otherwise adiabatic and the other vertical wall is cooled to a constant temperature. Non-dimensionalization of the governing equations results in four parameters for both aiding and opposing flows: (1) Ra, Rayleigh number (2) Pe, Péclet number (3) K_r, thermal conductivity ratio parameter, and (4) H, heat transfer coefficient parameter. The numerical results are presented for 0.01  H  100, 0.01  K_r  100, 0.01  Pe  100 and Ra = 10, 50 and 100. The results show that, the thermal equilibrium model cannot predict the average Nusselt number correctly for small values of H × K_r. In both the aiding and opposing flows, the total average Nusselt number is decreasing with increasing the heat transfer coefficient parameter at low values of Pe, while for high values of Pe, higher H will enhance the total heat transfer rate. Increasing the thermal conductivity ratio leads to increase in the total average Nusselt number. It is found also that the total average Nusselt number depends strongly on the thermal conductivity ratio parameter and depends slightly on the heat transfer coefficient parameter.     

Analytical investigation of heat transfer enhancement in a channel partially filled with a porous material under local thermal non-equilibrium condition

Forced convection through a channel partially filled with a porous medium is investigated analytically in the present work. Thermally developed condition is considered and the local thermal non-equilibrium model is utilized to obtain the exact solutions of both fluid and solid temperature fields for flow inside the porous material as well as for flow in the clear region. Nusselt number is obtained in terms of the porous insert thickness (S), porosity (?) as well as pertinent parameters such as thermal conductivity ratio (k), Biot number (Bi), and Darcy number (Da). The values of S by which the temperature difference between the two phases approach to zero, for different values of Bi, k, and Da number are obtained. It is found that three mechanisms affect the Nu number i: clear fluid conduction ii: internal heat exchange in the porous medium iii: channeling effect in the clear flow. The value of S, which yields the highest Nu number is found to vary linearly from 0.8 to 0.97 as the value of Da decreases from 10⁻³ to 10⁻⁷. At the expense of reasonable pressure drop the optimum thickness of porous material in order to enhance the heat transfer rate is found S = 0.8.