Direct and inverse mixed convections in an enclosure with ventilation ports

The numerical study presented in this work describes the direct and inverse mixed convection problems in a slot-ventilated enclosure subjected to an unknown heat flux on one side. Particularly, the interaction of internal natural convection with the cold ventilated flow leads to various flow fields depending on the Richardson number, Reynolds number, and the functional form of the imposed boundary heat flux. Fluid and heat transport structures across the enclosure are visualized by the streamlines and heatlines, respectively. Subsequently, an iterative conjugate gradient method is applied such that the gradient of the cost function is introduced when the appropriate sensitivity and adjoint problems are defined for a domain of arbitrary geometries. In this approach, no a priori information is needed about the unknown boundary heat fluxes to be determined. The accuracy of the heat flux profile solutions is shown to depend strongly on the values of Reynolds number and flux functional forms. Effects of measurement errors on the accuracy of estimation are also investigated. The present work is significant for the flow control simultaneously involving the natural convection and forced convection.     

The lattice Boltzmann modeling on the nanofluid natural convective transport in a cavity filled with a porous foam

In this work, the natural convective transport was numerically investigated for nanofluids in a metal-foam cavity. A lattice Boltzmann (LB) model for the nanofluid natural convection in a porous medium was established by using the volume-averaging method. The velocity and temperature fields were obtained, and flow and thermal characteristics of the nanofluid convection in a porous medium were presented. The effects of the Rayleigh Number, the Darcy Number, the porosity, the solid thermal conductivity of porous medium, the nanoparticle thermal conductivity and the nanoparticle concentration on natural convection were examined. The average velocity was put forward to evaluate the convection effect and the natural convection onset was also discussed. It is shown that the Nusselt number of the natural convection increases with an increase in the Darcy number, the Rayleigh number, the porosity and the effective thermal conductivity. The change from the heat conduction regime to the convection regime is clearly shown from the numerical result, which verifies the onset point of the nanofluid natural convection in a porous medium. The highly conductive porous foam and the nanofluid can promote the thermal performance of the natural convection, which own great potential in practical thermal applications.     

Identification of wall heat flux for turbulent forced convection by inverse analysis

An inverse problem for turbulent forced convection between parallel flat plates is investigated. The space- and time-dependent heat flux at the upper wall is estimated from the temperature measurements taken inside the flow. In the present study, the conjugate gradient method is adopted for the estimation of the unknown wall heat flux. No prior information is needed for the functional form of the wall heat flux in the inverse analysis. The effects of the measurement errors, the functional form of the wall heat flux, and the location of the sensors on the accuracy of the estimation are investigated. The reconstruction of the wall heat flux is satisfactory when simulated exact or noisy data are input to the inverse analysis. The sensitivity coefficients are discussed in this paper. As expected, it is shown that the accuracy of the estimation can be improved when the sensors are located closer to the upper wall.     

Visualization of heat flow using Bejan’s heatline due to natural convection of water near 4 °C in thick walled porous cavity

A numerical study on natural convection heat transfer of cold water near 4 °C in a thick bottom walled cavity filled with a porous medium has been performed. It is assumed that the cavity is isothermally heated from the outside of the thick bottom wall and cooled from ceiling. The finite-difference method has been used to solve the governing partial differential equations of heat and fluid flow. Effects of thermal conductivity ratio, Rayleigh number and bottom wall thickness on heat transfer from the bottom to the ceiling have been studied. The heatline visualization technique has been used to demonstrate the path of heat transport through the enclosure. Moreover, streamlines and isotherms have been used to present fluid flow and temperature distributions. The obtained results show that multiple circulation cells are formed in the cavity and the local Nusselt numbers at the bottom wall and solid–fluid interface are highly affected by formed cells. The increase of Rayleigh number and thermal conductivity ratio increases heat transfer through the cavity. However, the increase of thickness of the bottom wall reduces the mean Nusselt number. Almost one-dimensional conduction heat transfer is observed in the solid bottom wall of the cavity.     

Effect of horizontal pressure gradient on Rayleigh–Bénard convection of a Newtonian nanoliquid in a high porosity medium using a local thermal non-equilibrium model

A study of linear and weakly nonlinear stability analyses of Darcy–Brinkman convection in a water–alumina, nanoliquid-saturated porous layer for stress-free isothermal boundaries, when the solid and nanoliquid phases are in local thermal nonequilibrium, is conducted. The critical eigenvalue is found using the Galerkin approach. The effect of the pressure gradient, thermal conductivity ratio, interphase heat transfer coefficient, inverse Darcy number, and Brinkman number on the heat transport and onset of convection is examined and represented graphically. The critical values of wavenumber and nanoliquid Rayleigh number are found for different problem parameter values. The effect of increasing the porosity-modified ratio of thermal conductivity advances the onset of convection and increases the amount of heat transport, whereas the remaining parameters have the opposite impact on the onset of convection and amount of heat transport. The classical results of the local thermal equilibrium case and Darcy–Bénard convection in the presence of pressure gradient are obtained as a limiting case of the present problem.     

Inverse determination of building heating profiles from the knowledge of measurements within the turbulent slot-vented enclosure

Slot ventilated enclosure flows have been simulated, respectively in displacement ventilation and mixed ventilation covering from the forced convection dominated flow to the natural convection dominated flow. Direct convection simulation together with the turbulent streamlines and turbulent heatlines demonstrate that the enclosure flow pattern, indoor thermal level and heat transfer potential will depend on the interactions of external forced flow and thermal buoyancy driven flows, i.e., Reynolds number and Grashof number. In subsequent inverse convection modeling, the inverse determination of enclosure wall heat flux profiles was conducted by the use of adjoint methodology, in which the direct, sensitivity and adjoint problems are formulated and solved by finite volume method. The effects of the supplying air flow rate, thermal source strength, ventilation mode, flux functional forms, and the measurement errors on the accuracy of inverse turbulent convection estimation have been investigated. The inverse solutions of turbulent convections are of low level accuracy as the flow becomes thermal-driven turbulent flows, and they deteriorate as the noise levels increase. This work is of fundamental importance for the room air flow design and measurements involving the turbulent thermal convections.     

Vascularization with line-to-line trees in counterflow heat exchange

Here we report the heat and fluid flow characteristics of counterflow heat exchangers with tree-shaped line-to-line flow channels. The flow structures of the hot and cold sides are sequences of point-to-line trees that alternate with upside-down trees. The paper shows under what conditions the tree vascularization offers greater heat flow access than corresponding conventional designs with parallel single-scale channels. The analytical part is based on assuming fully developed laminar flow in every channel and negligible longitudinal conduction in the solid. The numerical part consists of simulations of three-dimensional convection coupled with conduction in the solid. It is shown that tree vascularization offers greater heat flow access (smaller global thermal resistance) than parallel channels when the number of pairing levels increases and the available pumping power or pressure drop is specified. When the solid thermal conductivity increases, the heat transfer effectiveness decreases because of the effect of longitudinal heat conduction. The nonuniformity in fluid outlet temperature becomes more pronounced when the number of pairing levels increases and the pumping power (or pressure drop number) increases. The nonuniformity in outlet fluid temperature decreases when the solid thermal conductivity increases.     

Heat transfer from small objects in microgravity: Experiments and analysis

Heat transfer over a sub-millimeter spheroidal solid is of interest in many engineering processes. One important mechanism of heat transfer in the above processes is natural convection which leads to heat transfer rates many times larger than that of pure conduction. Despite the huge literature devoted to natural convection heat transfer rates over spheres (and to a smaller extent over spheroids) there is not a generally accepted correlation especially for small Rayleigh numbers. Existing correlations for external geometries predict a progressively increasing contribution of natural convection to heat transfer with respect to gravity (starting from zero gravity). To test the validity of these correlations, experiments are performed for the estimation of heat transfer rates at low gravity. Heat pulses are given to a miniature thermistor with a nearly spheroidal shape immersed in a liquid and its thermal response is registered during heating in parabolic flights. The contribution of natural convection to heat transfer is undoubtedly estimated from runs in which acceleration varies from 0 to 1.8 g. Surprisingly enough, the experiments showed that the Rayleigh number must take a minimum value before non-negligible effect of natural convection on heat transfer appears (existence of a threshold Rayleigh number). In the absence of natural convection (below Ra_thr) the experimental thermal response curves can be successfully described by approximating solutions of the transient heat conduction equation for the spheroidal geometry of the thermistor. Apparently, additional research is needed regarding the natural convection around sub-millimeter objects for small Rayleigh numbers.     

Inverse Estimation of the Thermal Conductivity in a One-Dimensional Domain by Taylor Series Approach

The Taylor series approximation is developed for the inverse estimation of thermal conductivity in a one-dimensional domain. The differential governing equation of heat conduction is converted to a discrete system of linear equations in matrix form using the temperature measurement and heat generation at the grid points as well as surface heat flux. The unknown thermal conductivity is estimated by solving the linear algebraic equations directly without iterations. The features of the present method are that no prior information about the functional form of the thermal conductivity is required, nor are any initial guesses or iterations in the calculation process needed. The accuracy and robustness of the present method are verified by comparing the results with the analytical solutions for constant, spatial- and temperature-dependent thermal conductivities. The results show that the inverse solutions are in good agreement with the exact solutions.     

One-step GPS for the estimation of temperature-dependent thermal conductivity

In this paper we are concerned with the estimation of temperature-dependent thermal conductivity of a one-dimensional inverse heat conduction problem. First, we construct a one-step group-preserving scheme (GPS) for the semi-discretization of quasilinear heat conduction equation, and then derive a quasilinear algebraic equation to determine the unknown thermal conductivity under a given initial temperature and a measured temperature perturbed by noise at time T. The new method does not require any prior information on the functional form of thermal conductivity. Several examples are examined to show that the new approach has high accuracy and efficiency, and the number of iterations spent in solving the quasilinear algebraic equation is smaller than five even in a large temperature range.     

Double diffusive LTNE porous convection with Cattaneo effects in the solid

The impact of Cattaneo heat flux law in the solid on the onset of double‐diffusive Darcy porous convection with local thermal nonequilibrium temperatures is investigated. The Fourier law of heat transfer is invoked for the fluid, whereas the Cattaneo heat flux law used to transfer heat in solid skeleton alters the temperature equation from parabolic to hyperbolic. The results are obtained for porous skeletons of aluminum and copper oxides. Both Cattaneo and solute concentration effects reinforce in controlling the onset of oscillatory convection and some novel consequences are observed. Compared with the results perceived in the absence of solute concentration, a manifestation of oscillatory convection with scaled‐interphase heat transfer coefficient as well as solid thermal relaxation time parameter initiates earlier in its presence. The effect of increasing interphase heat transfer coefficient and the Lewis number is to delay and hasten the onset of stationary and oscillatory convection. Besides, the increase in the value of solid thermal relaxation time parameter advances the oscillatory onset. Although the increase in the solute Darcy–Rayleigh number is to delay the stationary onset, it shows a twofold behavior on the onset of oscillatory convection. Before the onset of oscillatory convection, the size of the convection cell gets narrower and after which it becomes much wider. The existing results are retrieved as limiting cases from the current study.     

Non-iteration estimation of thermal conductivity using finite volume method

The finite volume approach is developed for the inverse estimation of thermal conductivity in one-dimensional domain. The differential governing equation of heat conduction is converted to a system of linear equations in matrix form using the temperature data and heat generation at the discrete grid points as well as surface heat flux. The unknown thermal conductivities are obtained by solving the system equations directly. The features of the present method are that no prior information about the functional form of the thermal conductivity is required and no iterations in the calculation process are needed. The accuracy and robust of the present method are verified by comparing examples of inverse estimation of spatially and temperature-dependent thermal conductivities with the exact solutions.     

Conjugate Wall Conduction-Fluid Natural Convection in a Three-Dimensional Inclined Enclosure

Laminar conjugate conduction-natural convection heat transfer in a 3-D inclined cubic enclosure comprised of finite thickness conductive walls and central cavity is numerically investigated. The dimensionless governing equations describing the convective flow and wall heat conduction are solved by the high accuracy multidomain pseudospectral method. Computations are performed for different Rayleigh numbers (10³ ≤ Ra* ≤ 10⁶), thermal conductivity ratios (1 ≤ k ≤ 100), dimensionless wall thickness (0 ≤ s ≤ 0.25), and enclosure inclinations (?30° ≤ α ₁ ≤ 30°, 0° ≤ α ₂ ≤ 45°). The effects of the above controlling parameters on the heat transfer performances of the enclosure system are investigated in detail, with emphases on the variations of wall conduction and fluid convection heat transfer, and the interactive heat transfer conditions between solid walls and fluid in the central cavity. Numerical results reveal that the existence of enclosure walls reduces the temperature gradient across the cavity and alters the temperature distribution within the solid walls; thus, the fluid convection is complexly determined by the combined effects of k and s, and is greatly affected by enclosure inclinations at high Rayleigh numbers. Moreover, the temperature distributions and solid-fluid interactive heat transfer conditions are provided for further interpretation and demonstration of the effects of the solid walls.     

Effects of variable properties and non-uniform heating on natural convection flows in vertical channels

The effects of variable properties and non-uniform heating on laminar air flows induced by natural convection in vertical channels are investigated numerically. A full-elliptic model which accounts for variations in viscosity and thermal conductivity with temperature and which determines the density from the state equation has been applied to cases in which variable property effects cannot be neglected (including conditions for which flow reversals may occur). The influence of the Rayleigh number and non-uniformity of the wall heat flux distribution on the critical heat per unit time transferred from the walls in symmetric and asymmetric channels is analyzed for a wide range of Rayleigh numbers. It is shown that the maximum wall temperature can be reduced substantially by selecting an appropriate wall heat flux distribution, although at the cost of slightly more restrictive critical conditions.     

Numerical simulation of natural convection in a horizontal enclosure with a heat-generating conducting body

The physical model considered here is a horizontal layer of fluid heated below and cold above with heat-generating conducting body placed at the center of the layer. The dimensionless thermal conductivities of body considered in the present study are 0.1, 1 and 50. The dimensionless temperature difference ratios considered are 0.0, 0.25, 2.5 and 25. Two-dimensional solution for unsteady natural convection is obtained using an accurate and efficient Chebyshev spectral methodology for variety of Rayleigh number from 10³ to 10⁶. Multi-domain technique is used to handle square-shaped heat-generating conducting body. The fluid flow, heat transfer and time- and surface-averaged Nusselt number are investigated for various ranges of Rayleigh number, thermal conductivity ratio and dimensionless temperature difference ratio. The results for the case of conducting body with heat generation are also compared to those without heat generation to see the effects of heat generation from the conducting body on the fluid flow and heat transfer in the enclosure.     

Numerical study on coupled fluid flow and heat transfer process in parabolic trough solar collector tube

A unified two-dimensional numerical model was developed for the coupled heat transfer process in parabolic solar collector tube, which includes nature convection, forced convection, heat conduction and fluid-solid conjugate problem. The effects of Rayleigh number (Ra), tube diameter ratio and thermal conductivity of the tube wall on the heat transfer and fluid flow performance were numerically analyzed. The distributions of flow field, temperature field, local Nu and local temperature gradient were examined. The results show that when Ra is larger than 10⁵, the effects of nature convection must be taken into account. With the increase of tube diameter ratio, the Nusselt number in inner tube (Nu₁) increases and the Nusselt number in annuli space (Nu₂) decreases. With the increase of tube wall thermal conductivity, Nu₁ decreases and Nu₂ increases. When thermal conductivity is larger than 200 W/(m K), it would have little effects on Nu and average temperatures. Due to the effect of the nature convection, along the circumferential direction (from top to down), the temperature in the cross-section decreases and the temperature gradient on inner tube surface increases at first. Then, the temperature and temperature gradients would present a converse variation at θ near π. The local Nu on inner tube outer surface increases along circumferential direction until it reaches a maximum value then it decreases again.     

NATURAL CONVECTION IN VERTICAL PARALLEL PLATES WITH AN UNHEATED ENTRY OR UN HEATED EXIT

This study aims to investigate the effects of the unhealed entry or unheated exit section on the free convection heat transfer in airflow in vertical parallel plate channels resulting from the thermal boundary conditions of uniform heat flux (VHF) and uniform wall temperature (UWT). Results of average Nusselt number and dimensionless volume flow rate are presented in terms of the ratio of the length of heated section to the full channel length and a Rayleigh number, ranging from the limit for the fully developed flow to that for single-plate behavior. Analytical equations for dimensionless volume flow rate and average Nusselt number for both unheated restrictions and both thermal boundary conditions have been developed for the fully developed flow limit. The numerical solutions are shown to approach asymptotically the approximate solution for fully developed flow as the Rayleigh number approaches 1 or less. An important finding of the study is that an unheated exit characterizes greater total heat transfer and volume flow rate than an unheated entry does. The presence of the unheated entry or unheated exit severely affects the convection process, especially at low Rayleigh number. A notable effect of an unheated exit on convection characteristics was found for the case of UHF at high Rayleigh number.     

Inverse Identification of Thermal Properties of Charring Ablators

The modified Levenberg-Marquardt method is used for simultaneous estimation of decomposition kinetic coefficients and temperature-dependent thermophysical properties of charring ablators with a moving boundary over a wide temperature range. No prior information is used for the functional forms of the unknown thermal conductivity and specific heat. The procedure used differs from the traditional one in that it does not require prescribed time-dependent surface heat flux, recession rate, and pyrolysis gas mass flow rate. These time-dependent quantities may recover during an iterative procedure. The measured temperatures are simulated numerically by the Charring material ablation code, which accounts for unsteady ablation. The method can determine unknown parameters in an efficient manner with reasonable accuracy, without exact advance knowledge about the net surface heat flux, surface recession, and gas flux through the material.     

WALL HEAT CONDUCTION EFFECT ON NATURAL CONVECTION IN AN ENCLOSURE FILLED WITH A NON-DARCIAN POROUS MEDIUM

ABSTRACT

A numerical analysis has been made of the conjugate natural convection in a rectangular enclosure filled with a fluid-saturated porous medium and surrounded with four solid walls. The conductance of the walls is assumed to be much greater than that of the cavity filled with a porous medium. The main objective was to investigate the influences of the ratio of thermal conductivity of the wall to that of the fluid-porous matrix composite, the Darcy-modified Rayleigh number, the Prandtl number, and the aspect ratio. The streamlines and isotherms are presented; also, the local and average Nusselt numbers are presented along the interface between walls and cavity. A non-Darcian model was employed and the numerical method was SIMPLE-C. The numerical results indicate that the wall heat conduction effects decrease the heat transfer rate. When the wall heat conduction is considered, the greater the conductance of the solid walls surrounding the cavity, the greater is the rate of heat transfer.     

Non-uniform porosity and thermal dispersion effects on natural convection about a heated horizontal cylinder in an enclosed porous medium

A numerical solution has been obtained for transient two-dimensional natural convection from a heated horizontal cylinder embedded in an enclosed porous medium. Non-Darcian effects are taken into consideration in the momentum equation, while the thermal dispersion effect is taken into consideration in the energy equation. The wall effect on porosity is approximated by an exponential function and its effect on thermal dispersion is modeled by a dispersive length. The governing equations in terms of the stream function, vorticity, and temperature are expressed in a body-fitted coordinate system, which were solved numerically by the finite difference method. Results are presented for the streamlines and isotherms, tangential velocity and temperature distributions, as well as the average Nusselt numbers at different values of Rayleigh number, dimensionless particle diameter, and Prandtl number. The non-uniform porosity effect tends to increase the temperature gradient near the wall while the thermal dispersion effect increases the effective thermal conductivity, both resulting in an increase in surface heat flux. The effect of thermal dispersion on natural convection in porous media at low to moderate Rayleigh number is small. With nonuniform porosity and thermal dispersion effects taken into consideration, the predicted average Nusselt numbers are found to be in better agreement with experimental data.