Analysis of laminar flow forced convection heat transfer in the entrance region of a circular pipe

A numerical solution, for incompressible, steady-state, laminar flow heat transfer in the combined entrance region of a circular tube is presented for the case of constant wall heat flux and constant wall temperature. The development of velocity profile is obtained from Sparrow's entrance region solution. This velocity distribution is used in solving the energy equation numerically to obtain temperature profiles. Variation of the heat transfer coefficient for these two different boundary conditions for the early stages of boundary layer formation on the pipe wall is obtained. Local Nusselt numbers are calculated and the results are compared with those given byUlrichson andSchmitz. The effect of the thermal boundary conditions is studied by comparing the uniform wall heat flux results with uniform wall temperature.     

Oscillating heat pipe simulation considering dryout phenomena

In heat transport devices such as oscillating heat pipe (OHP), dryout phenomena is very important and avoided in order to give the optimum performance. However, from the previous studies (including our studies), the dryout phenomena in OHP and its mechanism are still unclear. In our studies of OHP (Senjaya and Inoue in Appl Thermal Eng 60:251–255, 2013; Int J Heat Mass Transfer 60:816–824, 2013; Int J Heat Mass Transfer 60:825–835, 2013), we introduced the importance and roles of liquid film in the operating principle of OHP. In our previous simulation, the thickness of liquid film was assumed to be uniform along a vapor plug. Then, dryout never occurred because there was the liquid transfer from the liquid film in the cooling section to that in the heating section. In this research, the liquid film is not treated uniformly but it is meshed similarly with the vapor plugs and liquid slugs. All governing equations are also solved in each control volume of liquid film. The simulation results show that dryout occurs in the simulation without bubble generation and growth. Dryout is started in the middle of vapor plug, because the liquid supply from the left and right liquid slugs cannot reach until the liquid film in the middle of vapor plug, and propagates to the left and right sides of a vapor plug. By inserting the bubble generation and growth phenomena, dryout does not occur because the wall of heating section is always wetted during the bubble growth and the thickness of liquid film is almost constant. The effects of meshing size of liquid film and wall temperature of heating section are also investigated. The results show that the smaller meshing size, the smaller liquid transfer rate and the faster of dryout propagation. In the OHP with higher wall temperature of heating section, dryout and its propagation also occur faster.     

Conjugate Natural Convection in a Porous Enclosure with Non-Uniform Heat Generation

Effects of a conductive wall on natural convection in a square porous enclosure having internal heating at a rate proportional to a power of temperature difference is studied numerically in this article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (0 ???Ra ???1000), the internal heating and the local exponent parameters (0 ????? ???5), (1 ????? ???3), the wall to porous thermal conductivity ratio (0.44 ???Kr ???9.9) and the ratio of wall thickness to its width (0.02 ???D ???0.5). The results are presented to show the effect of these parameters on the fluid flow and heat transfer characteristics. It is found a strong internal heating can generate significant maximum fluid temperature more than the conductive solid wall. Increasing value thermal conductivity ratio and/or decreasing the thickness of solid wall can increase the maximum fluid temperature. It is also found that at very low Rayleigh number, the heat transfer across the porous enclosure remain stable for any values of the thermal conductivity ratio.     

Binary diffusion and heat transfer in laminar mixed convection channel flows with uniform wall heat flux: extremely thin film thickness

A detailed numerical study has been performed to investigate the combined heat and mass transfer in laminar mixed convection channel flows with uniform wall heat flux. In an initial effort the liquid film on the channel wall is assumed to be extremely thin in thickness. Major dimensionless groups governing the present problem areGr _T,Gr _Mx,Pr,Sc, φ andRe. Results are specifically presented for an air-water system under various conditions. The effects of wall heating flux, the Reynolds number and the relative humidity of the moist air in the ambient on the momentum, heat and mass transfer in the flow are investigated in great detail.     

Optimal discrete distribution of heat flux elements for in-tube laminar forced convection

A new technique is proposed to enhance the heat transfer from a discretely heated pipe to a developing laminar fluid flow. Unlike the common heating situation where the fluid is continuously heated along the pipe wall with uniform heat flux, the proposed technique consists in heating the fluid with stepwise distributed heat flux, namely by placing insulated segments between the heated segments. Applying this technique, the effective length of the thermal entrance region is enlarged and as a result, the average heat transfer is invigorated. In order to maximize the heating performance, an optimal placement of the insulated segments between the heated segments is calculated according to constructal design. This serves to describe the optimal stepwise distribution of the heat flux. Owing that the total heat load is considered fixed, the maximization of the heating performance translates into the minimization of the peak temperature (‘hot spot’) of the pipe wall. The analytical results demonstrate that the optimal location of the insulated segments along with the reduction of the peak temperature strongly depend on the Graetz number. It is also shown that for intermediate values of the Graetz number, the peak temperatures are remarkably reduced in response to the optimal placement of the insulated/heated segments.     

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 pressure work and viscous dissipation on Graetz problem for gas flows in parallel-plate channels

The analytical solutions are obtained for the Graetz problem with pressure work and viscous dissipation in the thermal entrance region of the parallel-plate channels for two basic boundary conditions of uniform wall temperature and uniform wall heat flux involving fully developed laminar gas flows. The asymptotic Nusselt number is found to be zero instead of the conventionally accepted value of 7.54 for the uniform wall temperature case and (140/17)/ [1+(27/17) PrEc] for uniform wall heat flux case. The effects of pressure work and viscous dissipation contribute significantly to the asymptotic results for heat transfer and cannot be neglected under any circumstances in the case of uniform wall temperature. Sample results are presented to illustrate the effects of pressure work and viscous dissipation on heat transfer characteristics in the thermal entrance region.     

Viscous dissipation effects on thermal entrance region heat transfer in pipes with uniform wall heat flux

The analytical solution to Graetz problem with uniform wall heat flux is extended by including the viscous dissipation effect in the analysis. The analytical solution obtained reduces to that of Siegel, Sparrow and Hallman neglecting viscous dissipation as a limiting case. The sample developing temperature profiles, wall and bulk temperature distributions and the local Nusselt number variations are presented to illustrate the viscous dissipation effects. It is found that the role of viscous dissipation on thermal entrance region heat transfer is completely different for heating and cooling at wall. In the case of cooling at wall, a critical value of Brinkman number, Br _c=−11/24, exists beyond which (−11/24<Br<0) the fluid bulk temperature will always be less than the uniform entrance temperature indicating the predominance of cooling effect over the viscous heating effect. On the other hand, with Br < Br _c the bulk temperature T _b will approach the wall temperature T _w at some downstream position and from there onward the bulk temperature T _b becomes less than the wall temperature T _w with T _w > B _b > T ₀ indicating overall heating effect for the fluid. The numerical results for the case of cooling at wall Br < 0 are believed to be of some interest in the design of the proposed artctic oil pipeline.     

Unsteady entrance heat transfer in laminar pipe flows with convection from the ambient

The characteristics of unsteady entrance heat transfer in the combined entrance heat transfer region of laminar pipe flows resulting from time-varying inlet temperature are numerically investigated. Three non-dimensional parameters,Nu ₀, a^*, andf are identified in the study. Also, their effects on the non-dimensional duct wall temperature, fluid bulk temperature, and duct wall heat flux are discussed in great detail. Comparisons are made with the zero thermal capacity wall solution.     

Experimental and numerical study of developing turbulent flow and heat transfer in convergent/divergent square ducts

 The work reported in this paper is a systematic experimental and numerical study of friction and heat transfer characteristics of divergent/convergent square ducts with an inclination angle of 1^∘ in the two direction at cross section. The ratio of duct length to average hydraulic diameter is 10. For the comparison purpose, measurement and simulation are also conducted for a square duct with constant cross section area, which equals to the average cross section area of the convergent/divergent duct. In the numerical simulation the flow is modeled as being three-dimensional and fully elliptic by using the body-fitted finite volume method and the kɛ turbulence model. The uniform heat flux boundary condition is specified to simulate the electrical heating used in the experiments. The heat transfer performance of the divergent/convergent ducts is compared with the duct with uniform cross section under three constraints (identical mass flow rate, pumping power and pressure drop). The agreement of the experimental and numerical results is quite good except at the duct inlet. Results show that for the three ducts studied there is a weak secondary flow at the cross section, and the circumference distribution of the local heat transfer coefficient is not uniform, with an appreciable reduction in the four corner regions. In addition, the acceleration/deceleration caused by the cross section variation has a profound effect on the turbulent heat transfer: compared with the duct of constant cross section area, the divergent duct generally shows enhanced heat transfer behavior, while the convergent duct has an appreciable reduction in heat transfer performance. Received on 18 September 2000 / Published online: 29 November 2001     

On the calculation of wall temperatures in the post dryout heat transfer region

A non-equilibrium post dryout heat transfer model for calculating the wall temperature distribution in vertical upflows is presented in this study. The model is based upon the three path heat transfer formulation developed by MIT researchers (Laverty & Rohsenow 1964, Forslund & Rohsenow 1968, Hynek et al. 1969 and Plummer et al. 1974) that involves heat transfer from wall to vapor, from wall to droplets in contact with the wall and from vapor to liquid droplets in the vapor core. Downstream gradients for the bulk vapor temperature, vapor quality, droplet size and vapor velocities are identical to those used by Hynek et al. (1969) and Plummer et al. (1974). Conditions at the dryout location are calculated using a modified version of a technique developed by Hynek et al. (1969).A procedure for determining an average droplet diameter based on a size distribution is introduced. Migration of droplets through the boundary layer and droplet deposition flux are predicted with the model of Gani? & Rohsenow (1979). Heat transfer from the wall to the impinging liquid droplets is calculated with a correlation by Holman & McGinnis (1969). Mechanisms contributing to wall to droplet heat transfer are identified as (a) droplet-wall contact, (b) intensive droplet evaporation inside the boundary layer, and (c) destruction of the boundary layer due to droplet migration to, and rebound from, the hot surface. The significance of the average droplet size and size distribution is demonstrated through its control over the free stream evaporation and droplet deposition rates.Predicted uniform heat flux wall temperature profiles for water, nitrogen and freon 12 are in good agreement with the data of Era et al. (1966), Bennett et al. (1967), Forslund & Rohsenow (1968), Ling et al. (1971), Groeneveld (1972) and Janssen & Kervinen (1975).     

Heat dispersion in stationary mixed convection flow about horizontal surfaces in porous media

An analysis has been performed to study the influence of velocity dependent dispersion on transverse heat transfer in mixed convection flow above a horizontal wall of prescribed temperature in a saturated porous medium. The Boussinesq approximation and boundary layer analysis were used to numerically obtain gravity affected temperature and velocity distributions within the frames of Darcy's law and a total thermal diffusivity tensor comprising both of constant coefficient heat conduction and velocity proportional mechanical heat dispersion. Dependending on Pe_∞, the molecular Peclét number basing on the effective thermal diffusivity and the velocity of the oncoming flow, density coupling has distinct influences on heat transfer rates between the wall surface and the porous medium flow region. For small Peclét numbers, when heat conduction is the prevailing mechanism, wall heat fluxes are the higher the larger the density difference between the oncoming and the near wall fluid is. The opposite is true for larger Peclét numbers, when mechanical heat dispersion is the main cause of heat spreading. For Pe_∞ tending to infinity these wall heat fluxes approach finite maximum values in the total heat diffusivity model, they grow beyond any limit if only constant coefficient heat conduction is considered. Thus, the inclusion of mechanical heat dispersion effects yields physically more realistic predictions. Received on 18 September 1996     

Analysis of laminar thermal entrance region of elliptical and rectangular channels with Kantorowich method

The Kantorowich method of variational calculus is introduced to solve the problem of laminar forced heat convection in a channel of an arbitrary cross section. Employing this procedure, the development of the temperature field and the heat transfer data in the thermal entrance region of elliptical and rectangular channels are determined, assuming a linear variation of wall temperature in the direction of flow. The uniform and the fully developed velocity profiles are considered. The local Nusselt numbers are tabulated for different aspect ratios.     

Advanced near-wall modeling for engine heat transfer

Recent developments in the engine heat transfer modeling tend to improve existing wall heat transfer models (temperature wall functions) which mostly rely on the standard or low-Re variants of k-ε turbulence model. Presently applied mesh resolutions already allow for first near-wall computational cells reaching the buffer or locally even viscous/conductive sub-layer, thus increasing the importance of more sophisticated modeling approach. As temperature gradient-induced density and fluid property variations become significant, wall heat transfer is strongly influenced by property variations (viscous/conductive sub-layer) and predictive capability of the turbulence model (buffer region), standard wall laws being inadequate anymore, even for attached boundary layers. The present approach relies on the k-ζ-f turbulence model and formulates a compressible wall function of Han and Reitz in the framework of hybrid wall treatment. The model is validated against spark ignition (SI) engine heat transfer measurements. Predicted wall heat flux evolutions on the cylinder head exhibit very good agreement with the experimental data, being superior to similar numerical predictions available in the published literature.     

A heat transfer measurement of jet impingement with high injection temperature

Local heat transfer from an impinging high temperature jet is studied using a method based on the heat thin foil technique and on the infrared thermography. Heat thin foil technique is used to impose several heat fluxes. For each flux, the temperature distribution is recorded using infrared imaging. Local heat transfer coefficients and adiabatic wall temperatures are determined by means of a linear regression method. This procedure is validated for a single round jet impinging on a flat plate for a range of injection temperatures. To cite this article: M. Fénot et al., C. R. Mecanique 333 (2005).     

Flow and heat transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet subject to a transverse magnetic field

An analysis is performed for flow and heat transfer of a steady laminar boundary layer flow of an electrically conducting fluid of second grade in a porous medium subject to a transverse uniform magnetic field past a semi-infinite stretching sheet with power-law surface temperature or power-law surface heat flux. The effects of viscous dissipation, internal heat generation of absorption and work done due to deformation are considered in the energy equation. The variations of surface temperature gradient for the prescribed surface temperature case (PST) and surface temperature for the prescribed heat flux case (PHF) with various parameters are tabulated. The asymptotic expansions of the solutions for large Prandtl number are also given for the two heating conditions. It is shown that, when the Eckert number is large enough, the heat flow may transfer from the fluid to the wall rather than from the wall to the fluid when Eckert number is small. A physical explanation is given for this phenomenon.     

Mixed convection film condensation from downward flowing vapors onto a sphere with variable wall temperature

A model is developed for the study of mixed- convection film condensation from downward flowing vapors onto a sphere with variable wall temperature. The model combined natural convection dominated and forced convection dominated film condensation, concerning effects of pressure gradient (P), interfacial vapor shear drag and non-uniform wall temperature variation (A), has been investigated and solved numerically. The effect of pressure gradient on the dimensionless mean heat transfer, NuˉRe^−1/2 will remain almost uniform with increasing P until for various corresponding available values of F. Meanwhile, the dimensionless mean heat transfer, NuˉRe^−1/2 is increasing significantly with F for its corresponding available values of P. Although the non-uniform wall temperature variation has an appreciable influence on the local film flow and heat transfer; however, the dependence of mean heat transfer on A can be almost negligible. Received on 10 October 1996     

Analysis of viscous dissipation effect on thermal entrance heat transfer in laminar pipe flows with convective boundary conditions

Consideration is given to the influence of viscous dissipation on the thermal entrance region laminar pipe flow heat transfer with convective boundary condition. The Eigenfunction series expansion technique is employed to solve the governing energy equation. The results for axial distributions of dimensionless bulk and wall temperatures, local Nusselt number as well as modified local Nusselt number are presented graphically forNu ₀ =0.1, 2, and 100. The complicated variations of conventional local Nusselt number is due to the inappropriate definition of conventional heat transfer coefficient in this problem. A modified local heat transfer coefficient, based on the difference of bulk fluid temperature and wall temperature, is introduced. Its value can clearly indicate the extent and the direction of heat exchange between the fluid in the pipe and the ambient. The effects of outside Nusselt number are also investigated. Significant viscous dissipation effects have been observed for large Br.     

Single bubble growth in saturated pool boiling of binary mixtures

Nucleate pool boiling experiments for binary mixtures, which are consisted of R11 and R113, were performed with constant wall temperature condition. Results for binary mixtures were also compared with pure fluids. A microscale heater array and Wheatstone bridge circuits were used to maintain the constant temperature of the heating surface and to obtain heat flow rate measurements with high temporal and spatial resolutions. Bubble growth images were captured using a high-speed CCD camera synchronized with the heat flow rate measurements.The departure time for binary mixtures was longer than that for pure fluids, and binary mixtures had a higher onset of nucleate boiling (ONB) temperature than pure fluids. In the asymptotic growth region, the bubble growth rate was proportional to a value between t^1/6 and t^1/4. The bubble growth behavior was analyzed to permit comparisons with binary mixtures and pure fluids at the same scale using dimensionless parameters. There was no discernible difference in the bubble growth behavior between binary mixtures and pure fluids for a given ONB temperature. And the departure radius and time were well predicted within a ±30% error.The minimum heat transfer coefficient of binary mixtures occurred near the maximum |y−x| value, and the average required heat flux during bubble growth did not depend on the mass fraction of R11 as more volatile component in binary mixtures. Finally, the results showed that for binary mixtures, a higher ONB temperature had the greatest effect on reducing the heat transfer coefficient.