Flow of Colloid Liquid in Horizontal Cell under Heating from Sidewall

Based on numerical simulation, the influence of sedimentation on the convective flow of colloid liquids filling a horizontal cell heated from the sidewall is studied. The set of nonlinear equations is solved by the finite-difference method using explicit schemes. Three convective patterns differing in spatial structure and behavior in time are distinguished. The transition between the patterns is accompanied by a jump in the dimensionless heat flow. Bifurcation diagrams of the convection patterns (the dependences of the heat flow intensity on the Rayleigh number) are given. It is shown that the weak flow of a colloidal suspension exists at a low temperature gradient, the intensity of which is several orders of magnitude lower than the intensity of the flow of a homogeneous liquid under the same parameters. The concentration in the flow with a weak intensity is redistributed in such a way that the density gradient becomes almost vertical, and the heat flow across the layer is absent at the same time. The transition from a weak to a strong one-vortex flow filling the entire cell proceeds abruptly. It is found the threshold of the transition from a weak to intense flow depends on the Boltzmann number characterizing the degree of gravitational stratification. One more flow, namely, a three-vortex flow with an intermediate intensity is generated upon a decrease in the Rayleigh number. Stream-function and concentration fields are manifested for all the observed types of flows.     

Complex surface rays associated with inhomogeneous skimming and Rayleigh waves

The Rayleigh wave, that propagates at the free surface of semi-infinite anisotropic medium, is composed of three inhomogeneous partial waves, each propagating along the surface with a different attenuation along the depth. Since this wave does not exhibit an attenuation on the surface, let us call it the homogeneous Rayleigh wave. The associated slowness corresponds to the real solution of the Rayleigh dispersion equation. Besides this classical solution, an infinite number of complex solutions of the Rayleigh dispersion equation exits. For such particular Rayleigh waves, the slowness vector, i.e. the identical component on the surface of the slowness of each partial waves, is taken to be complex. Thus, these Rayleigh waves are attenuated on the surface and as shown here, their attenuation is normal to the ray direction (or the energy velocity direction). Similarly to the infinite inhomogeneous plane waves which can be associated with complex rays, we call these waves, inhomogeneous Rayleigh waves. We use the inhomogeneous skimming waves, which are inhomogeneous plane waves, and the inhomogeneous Rayleigh waves to explain differently the usual diffraction phenomena on the free surface which cannot be explained by the real ray theory. For example, the arrival time of the wave packet observed beyond the cusp is in perfect accordance with the arrival time of some specific inhomogeneous Rayleigh waves. We show that these results are in agreement with the computation of the Green function. They apply to the theory of surface waves in linear elastodynamics with intrinsic anisotropy as well as to the theory of surface waves in linearised (incremental) elastodynamics with strain-induced anisotropy (also known as small-amplitude waves superimposed on the large static homogeneous deformation of a non-linear solid).     

Effect of a latitudinal temperature gradient on convective motions arising in a rotating spherical layer

The hydrodynamics of planetary atmospheres and Interiors are frequently directly or indirectly connected with convective motions taking place in rotating liquid spherical layers in the field of a central force. Convective stability in a spherical layer at rest, in a central gravity field, was first discussed in [1, 2]. It was shown that the critical Rayleigh number Rao at which convective instability sets in and the wave number of the critical perturbations depend essentially on the thickness of the layer. As in the plane case, the problem of the convective stability of a spherical layer is found to be degenerate, and the form of the critical perturbations cannot be determined from the linear problem. In actuality, minimization of the Rayleigh number permits establishing only the wave numberl for the spherical harmonic Y _l ^m (θ, ?), realized at the limit of stability; the parameter m remains indeterminate and thus 2l+1 independent convective modes correspond to Rao. In [3] a study was made of the convective stability of a liquid in a slowly rotating thin spherical layer. It was shown that the presence of rotation eliminates the degeneracy; at the limit of stability there arise motions corresponding to the Y _l l (θ, ?) -harmonic with a degenerate maximum at the equator, and propagating in a wave manner toward the side opposite to the rotation. In the present work a study is made of the convective stability of a flow of liquid, arising in a rotating spherical layer due to a nonuniform distribution of the temperatures at one of the boundaries of the layer. In such a statement of the problem it is possible to model large-scale motions in the atmospheres of large planets having internal sources of heat and absorbing solar radiation near the cloud cover of the atmosphere. It is established that, depending on the relationships between the parameters imparting the rotation and the inhomogeneous distribution of the temperature, there is either stabilization or destabilization of the layer in comparison with a fixed layer of the same thickness and with the same, but uniformly distributed heat flux supplied to the layer. A study is made of the form of the corresponding critical perturbations.     

Thermal Non-equilibrium Natural Convection in a Square Enclosure with Heat-Generating Porous Layer on Inner Walls

Differentially heated enclosure with heat-generating porous layer on inner walls is studied computationally for non-Darcy flow and thermal non-equilibrium models. In this study, this problem is investigated for different internal and external Rayleigh numbers, Darcy numbers, porosity-scaled thermal conductivity ratio, solid-/fluid-scaled heat transfer coefficient and dimensionless thickness of the porous layer. The results indicate that the dimensionless thickness of the porous layer has an important effect on the heat transfer in the enclosure. It was found that the thermal non-equilibrium model is needed for small values of the porosity-scaled thermal conductivity ratio and the solid-/fluid-scaled heat transfer coefficient. It is shown that the convection of heat due to internal heat generation is increased in the enclosure when the ratio of internal Rayleigh number to external Rayleigh number is larger.     

Transport through baffles in bottom heated top cooled enclosures: parametric studies

This paper presents parametric studies on the heat transfer and fluid exchange through single-hole baffles located at the median height in bottom heated top cooled enclosures. Results indicate that when the baffle area-opening ratio is smaller than 2%, the heat transfer in the enclosure is dominated by the transport through the baffle opening. Even with such small baffle openings, increasing the enclosure aspect ratio still enhances the transport across the baffle. The characteristic length scale of flow in the enclosure is a combination of baffle opening diameter and the chamber height. The Nusselt number that characterize the heat transfer through the baffle-hole is linearly correlated with the Rayleigh number based on baffle opening diameter and the temperature difference between the bulk temperatures in the two chambers, while no effects of Prandtl numbers are observed. The mechanism of transport across the baffle opening varies from conduction dominated, combined conduction and convection, and convection dominated regimes as Rayleigh number increases.     

Computational fluid dynamics applied to internal gas-turbine blade cooling: a review

This paper reviews current capabilities for predicting flow in the cooling passages and cavities of jet engines. Partly because of the need to enhance heat transfer coefficients, these flow domains entail complicated passage shapes where the flow is turbulent, strongly three-dimensional (3-D) and where flow separation and impingement, complicated by strong effects of rotation, pose severe challenges for the modeler. This flow complexity means that more elaborate models of turbulent transport are needed than in other areas of turbine flow analysis. The paper attempts to show that progress is being made, particularly in respect to the flow in serpentine blade-cooling passages. The first essential in modeling such flows is to adopt a low Reynolds number model for the sublayer region. The usual industrial practice of using wall functions cannot give a better than qualitative account of effects of rotation and curvature. It is shown that Rayleigh number effects can modify heat transfer coefficients in the cooling passages by at least 50%. The use of second-moment closure in the modeling is shown to be bringing marked improvements in the quality of predictions. Areas where, at present, more computational fluid dynamics (CFD) applications are encouraged are impingement cooling and pin-fin studies.     

Turbulent thermal convection in wide horizontal fluid layers

Turbulence in thermal convection is investigated for flows in which the production of turbulence energy is due solely to buoyancy, and the statistics of the flow are homogeneous in horizontal planes. New experimental results for high Rayleigh number unsteady turbulent convection in a horizontal layer heated from below and insulated from above are presented and compared to turbulent Rayleigh convection, convection in the planetary boundary layer, and laboratory penetrative convection. Mean temperature fields are correlated in terms of wall layer scales and convection scales. Joint statistics of turbulent temperature and horizontal velocity and vertical velocity through fourth order are presented for the core region of the convection layer.This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Bounds on convective heat transport in a rotating porous layer

Using the background field variational method, bounds on convective heat transport in a rotating porous layer heated from below are derived from the primitive equations. The enhancement of heat transport beyond the minimal conduction value (the Nusselt number Nu) is bounded in terms of the dimensionless temperature difference across the layer (the Rayleigh number Ra) according to

This rigorous upper bound shows that rotation has a retarding effect on convective heat transport.     

Natural convection in an inclined enclosure containing internal energy sources and cooled from below

Natural convection in an inclined enclosure from below and containing internally heated fluid has been investigated using a finite difference calculation procedure. Results have been obtained for Rayleigh number values up to 10⁶ and for inclination angles of 30 and 60°. For internal Rayleigh numbers that are much larger than the external Rayleigh number, the flow rises in the interior and moves down both the hot and cold walls. On the other hand, if the external Rayleigh number has a larger magnitude, the flow moves upwards along the hot surface and downwards along the cold surface. For the latter situation, the inner core is multicellular in nature at large external Rayleigh numbers. The average heat flux ratio along the cold surface (convective heat flux/corresponding conduction heat flux) increases with increasing external Rayleigh number and decreasing internal ratio is non-monotonic in nature. The heat flux ratio along both surfaces is observed to be strongly dependent on the inclination angle at high external Rayleigh numbers. A maximum in the local heat flux along the cold surface is obtained in the vicinity of x/L = 1 where hot fluid, either from the interior or directly from the opposite hot wall, meets the surface. Along the hot wall, a maximum in the heat flux ra flo     

Double-diffusive parallel flow induced in a horizontal Brinkman porous layer subjected to constant heat and mass fluxes: analytical and numerical studies

Fluid flow and heat and mass transfer induced by double-diffusive natural convection in a horizontal porous layer subjected to vertical gradients of temperature and concentration are studied analytically and numerically using the Brinkman-extended Darcy model. Both cases of rigid and free horizontal boundaries are examined in the present study. The parameters governing the problem are the Rayleigh number R_T, the Lewis number Le, the buoyancy ratio N, the Darcy number Da and the aspect ratio Ar. The analytical solution is based on the parallel flow approximation. The critical Rayleigh number corresponding to the onset of the parallel flow in this system is determined analytically as a function of Le, N and Da. For sufficiently small Da, both free and rigid boundaries yield results which are identical to those predicted by the Darcy model. The present investigation shows that there exists a region in the plane (N, Le) where the convective flow is not possible in the layer regardless of the Rayleigh and Darcy numbers considered. Received on 21 December 1998     

Stability of Flow with Viscous Dissipation in a Horizontal Porous Layer with an Open Boundary Having a Prescribed Temperature Gradient

A buoyancy-induced stationary flow with viscous dissipation in a horizontal porous layer is investigated. The lower boundary surface is impermeable and subject to a uniform heat flux. The upper open boundary has a prescribed, linearly varying, temperature distribution. The buoyancy-induced basic velocity profile is parallel and non-uniform. The linear stability of this basic solution is analysed numerically by solving the disturbance equations for oblique rolls arbitrarily oriented with respect to the basic velocity field. The onset conditions of thermal instability are governed by the Rayleigh number associated with the prescribed wall heat flux at the lower boundary, by the horizontal Rayleigh number associated with the imposed temperature gradient on the upper open boundary, and by the Gebhart number associated with the effect of viscous dissipation. The critical value of the Rayleigh number for the onset of the thermal instability is evaluated as a function of the horizontal Rayleigh number and of the Gebhart number. It is shown that the longitudinal rolls, having axis parallel to the basic velocity, are the most unstable in all the cases examined. Moreover, the imposed horizontal temperature gradient tends to stabilise the basic flow, while the viscous dissipation turns out to have a destabilising effect.     

Laminar natural convection in a partially divided rectangular cavity at high Rayleigh number

Finite element predictions of two-dimensional laminar natural convection in a partially divided rectangular cavity at high Rayleigh number are presented. The walls are differentially heated, the horizontal surfaces are insulated and the cavity contains a partial vertical divider which is centrally located and whose height is varied. Detailed results are presented for an aluminium half-divider in water for Rayleigh number up to 10¹¹ and compared directly with recent experiments in a cavity of aspect ratio 1/2. The predicted flow and heat transfer are in good agreement with the measurements and confirm the existence of a high Rayleigh number regime with characteristic behaviour that differs significantly from that found at lower Rayleigh number. In addition, the effects of the divider height, the divider conductivity, the fluid Prandtl number and the cavity aspect ratio are studied. The results show that a direct simulation of the complex flow and heat transfer that occurs in partially divided cavities is possible for realistic physical conditions.     

Experimental study of combined free and forced laminar convection in an asymmetrically heated two-dimensional aiding flow

An experimental study of mixed convection in an asymmetrically heated two-dimensional flow of water has been made. Experiments in fully developed turbulent flow (Re maximum=9000) were made initially to establish the satisfactory operation of the equipment. In the mixed convection laminar flow tests to determine local Nusselt numbers, the Reynolds number range was from about 100 to 1000 with the Rayleigh number based on heat flux varying from about 10 to 300. The data, which are presented graphically, have been compared with the results of the theoretical analyses of other workers. The effect of buoyancy in the immediate entry of the flow [(x/de)/Re Pr <.05] was found to be negligible, but at greater distances there were significant increases in the heat transfer coefficient with the larger values of Rayleigh number.     

Thermal Instability in a Horizontal Porous Channel with Horizontal Through Flow and Symmetric Wall Heat Fluxes

The linear thermoconvective instability of the basic parallel flow in a plane and horizontal porous channel is investigated. The boundary walls are assumed to be impermeable and subject to symmetric and uniform heat fluxes. The wall heat fluxes produce either a net heating or a net cooling of the fluid saturated porous medium. A horizontal mass flow rate is externally impressed leading to a stationary basic state with a temperature gradient inclined to the vertical. A region of possibly unstable thermal stratification exists either in the lower half-channel (boundary heating), or in the upper half-channel (boundary cooling). The convective instability of the basic flow is governed by the Rayleigh number and by the Péclet number. In the case of boundary heating, the thermal instability arises when the Rayleigh number exceeds its critical value, that depends on the Péclet number. The change of the critical Rayleigh number as a function of the Péclet number is determined numerically for arbitrary normal modes oblique to the basic flow direction. The most dangerous modes are the longitudinal rolls, with a wave vector perpendicular to the basic velocity. There exists a minimum value of the Péclet number, 19.1971, below which no linear instability is detected.     

Effects of Time-Periodic Thermal Boundary Conditions and Internal Heating on Heat Transport in a Porous Medium

The effects of time-periodic boundary temperatures and internal heating on Nusselt number in the Bénard–Darcy convective problem has been considered. The amplitudes of temperature modulation at the lower and upper surfaces are considered to be very small. By performing a weakly non-linear stability analysis, the Nusselt number is obtained in terms of the amplitude of convection, which is governed by the non-autonomous Ginzburg–Landau equation, derived for the stationary mode of convection. The effects of internal Rayleigh number, amplitude and frequency of modulation, thermo-mechanical anisotropies, and Vadasz number on heat transport have been analyzed and depicted graphically. Increasing values of internal Rayleigh number results in the enhancement of heat transport in the system. Further, the study establishes that the heat transport can be controlled effectively by a mechanism that is external to the system.     

Natural convection of SiO_2-water nanofluid in square cavity with thermal square column

A square with a thermal square column is a simple but nontrivial research prototype for nanofluid research. However, until now, the effects of the temperature of the square column on the heat and mass transfer of nanofluids have not been revealed comprehensively, especially on entropy generation. To deepen insight into this important field, the natural convection of the SiO_2-water nanofluid in a square cavity with a square thermal column is studied numerically in this study. The effects of the thermal column temperature(T = 0.0, 0.5, 1.0, 1.5), the Rayleigh number(ranging from 10~3 to 10~6),and the volume fraction of the nanoparticle(varying from 0.01 to 0.04) on the fluid flow,heat transfer, and entropy generation are investigated, respectively. It is found that, no matter at a low or high Rayleigh number, the volume fraction of the nanoparticle shows no considerable effects on the flow field and temperature field for all the temperatures of the thermal column. With an increase in the volume fraction, the mean Nusselt number increases slightly. At the same time, it is found that, with an increase in the temperature of the thermal column, the average Nusselt number gradually decreases at all values of the Rayleigh number. Meanwhile, it is found that, at a high Rayleigh number, the heat transfer mechanism is the main parameter affecting the increase in the total entropy generation rather than the volume fraction. In addition, no matter at a high or low Rayleigh number, when T = 0.5, the total entropy generation is the minimum.     

Thermal stratification studies in a side heated water pool for advanced heavy water reactor applications

Strong heat source at the isolation condenser wall of an Advanced Heavy Water Reactor, results in natural convection in gravity driven water pool, which leads to a thermally stratified pool. Governing equations simulating fluid flow and heat distribution are solved numerically by a general purpose Computational Fluid Dynamics solver developed at Indian Institute of Technology, Kanpur. Incompressible finite volume method with non-staggered grid arrangement is used in this exercise. This algorithm is fully implicit and semi-coupled. Turbulent natural convection in a boundary layer for high Rayleigh numbers is analyzed by the Lam–Bremhorst k − ε turbulence model. Analysis of unsteady laminar natural convection in a side-heated water cavity is also done for different values of Rayleigh number. Results show a warm fluid layer floating on the top of gradually colder layer (along the vertical direction) that indicates the presence of thermal stratification phenomenon. This fact necessitates additional safety features in such a system so that the detrimental effect such as stratification is minimized.     

Turbulent penetrative convection with an internal heat source

An infinite fluid with a vertical cubic temperature profile in the absence of fluid motion is considered as a model for penetrative convection in which a central unstably stratified fluid layer is bounded above and below by stably stratified layers. Turbulence statistics from direct and large eddy numerical simulations for the mean temperature gradient, the velocity and temperature variances and the heat flux are presented for Rayleigh numbers R up to four orders of magnitude above critical. By means of a simplified second-moment closure, analytical scaling laws for the statistics are determined. For high Rayleigh numbers, the mean temperature gradient approaches zero in a central well-mixed layer, a reduced positive (stable) value in upper and lower partially mixed layers, and an unmixed value far above and below. The temperature variance is a factor of R^1/3 larger in the partially mixed layers compared to the well-mixed layer; the velocity variance and heat flux scales the same in both layers. Approximation of the three layers by a two layer model yields an estimate for the height of the mixed layer: the height decreases slowly with increasing Rayleigh number and at the highest Rayleigh number simulated is approximately 30% longer than the unstable layer in the absence of fluid motion.     

Investigation of the heat transfer mechanism in supersonic turbulent boundary layers

A method is described for calculating turbulent Prandtl numbers from Mach number and total temperature profiles in supersonic boundary layers. The calculations are based on boundary layer measurements in the Mach number range from 3.5 to 5. The investigations clearly indicate that in addition to accurate profile measurements reliable values of shear stress and heat flux at the wall must exist, in order to be able to calculate the turbulent Prandtl number in the viscous regime of the boundary layer. For flow conditions with and without heat transfer, the derived turbulent Prandtl numbers indicate that the turbulent transport of heat decreases much faster towards the wall than the turbulent transport of momentum. The results of the analysis show that only the unequivocal qualitative result of increasing turbulent Prandtl numbers in the viscous region of the boundary layer, can be expected. The variation of the turbulent Prandtl number can be described successfully using a simple approximation, based on the mixing length concept, and is applied to the calculation of total temperature distribution using the law of the wall for compressible flow.     

Rayleigh Surface Waves on a Kelvin-Voigt Viscoelastic Half-Space

In this paper we consider the propagation of Rayleigh surface waves in an exponentially graded half-space made of an isotropic Kelvin-Voigt viscoelastic material. Here we take into account the effect of the viscoelastic dissipation energy upon the corresponding wave solutions. As a consequence we introduce the damped in time wave solutions and then we treat the Rayleigh surface wave problem in terms of such solutions. The explicit form of the secular equation is obtained in terms of the wave speed and the viscoelastic inhomogeneous profile. Furthermore, we use numerical methods and computations to solve the secular equation for some special homogeneous materials. The results sustain the idea, existent in literature on the argument, that there is possible to have more than one surface wave for the Rayleigh wave problem.