Semi-empirical modeling of pool boiling heat transfer in binary mixtures

Pool boiling heat transfer has been investigated for various binary mixtures, including acetone/isopropanol, water/acetone, water/methanol, water/ethanol, water/isopropanol, water/monoethanolamine, water/diethanolamine and water/triethyleneglycol as test solutions. Many correlations have been developed to predict the pool boiling heat transfer coefficient in mixtures in the past few decades, however the predicted values are not confirming. In addition, the application of many existing correlations requires some individual adjusting parameters that may be not available for every system. In this investigation, a new set of experimental data are presented. These data have been compared to major existing correlations. It is observed that the pool boiling heat transfer coefficients in mixtures are less than the ideal boiling heat transfer coefficient. A new semi-empirical model has been proposed based on the mass transfer resistance to predict the boiling heat transfer coefficient with satisfactory accuracy. The new model does not include any tuning parameter and is applicable to any given binary system. The performance of the proposed model is superior to most existing correlations.     

Process of selective laser sintering of polymer powders: Modeling,simulation, and validation

Selective laser sintering (SLS) of polymer powders involves multiphysical transient phenomena. A numerical tool for simulating such a process is developed on the basis of the reliable modeling of the corresponding thermo-physical transient phenomena and appropriate numerical methods. The present paper addresses modeling, simulation, and validation aspects that are indispensable for studying and optimizing SLS process. The coupled multiphysical models are detailed, and the numerical tool based on the finite volume method is presented, with validations in terms of numerical and physical accuracy, by considering the shrinkage involved in the process and the successive layers deposition. A parametric analysis is finally proposed in order to test the reliability of the model in terms of representing real physical phenomena and thermal history experienced by the material during the process.     

Concepts and realization of microstructure heat exchangers for enhanced heat transfer

Microstructure heat exchangers have unique properties that make them useful for numerous scientific and industrial applications. The power transferred per unit volume is mainly a function of the distance between heat source and heat sink—the smaller this distance, the better the heat transfer. Another parameter governing for the heat transfer is the lateral characteristic dimension of the heat transfer structure; in the case of microchannels, this is the hydraulic diameter. Decreasing this characteristic dimension into the range of several 10s of micrometers leads to very high values for the heat transfer rate.

Another possible way of increasing the heat transfer rate of a heat exchanger is changing the flow regime. Microchannel devices usually operate within the laminar flow regime. By changing from microchannels to three dimensional structures, or to planar geometries with microcolumn arrays, a significant increase of the heat transfer rate can be achieved.

Microheat exchangers in the form of both microchannel devices (with different hydraulic diameters) and microcolumn array devices (with different microcolumn layouts) are presented and compared. Electrically heated microchannel devices are presented, and industrial applications are briefly described.     

Computational modeling of biomagnetic micropolar blood flow and heat transfer in a two-dimensional non-Darcian porous medium

We study theoretically and computationally the incompressible, non-conducting, micropolar, biomagnetic (blood) flow and heat transfer through a two-dimensional square porous medium in an (x,y) coordinate system, bound by impermeable walls. The magnetic field acting on the fluid is generated by an electrical current flowing normal to the x–y plane, at a distance l beneath the base side of the square. The flow regime is affected by the magnetization B ₀ and a linear relation is used to define the relationship between magnetization and magnetic field intensity. The steady governing equations for x-direction translational (linear) momentum, y-direction translational (linear) momentum, angular momentum (micro-rotation) and energy (heat) conservation are presented. The energy equation incorporates a special term designating the thermal power per unit volume due to the magnetocaloric effect. The governing equations are non-dimensionalized into a dimensionless (ξ,η) coordinate system using a set of similarity transformations. The resulting two point boundary value problem is shown to be represented by five dependent non-dimensional variables, f _ξ  (velocity), f _η (velocity), g (micro-rotation), E (magnetic field intensity) and θ (temperature) with appropriate boundary conditions at the walls. The thermophysical parameters controlling the flow are the micropolar parameter (R), biomagnetic parameter (N _H ), Darcy number (Da), Forchheimer (Fs), magnetic field strength parameter (Mn), Eckert number (Ec) and Prandtl number (Pr). Numerical solutions are obtained using the finite element method and also the finite difference method for Ec=2.476×10⁻⁶ and Prandtl number Pr=20, which represent realistic biomagnetic hemodynamic and heat transfer scenarios. Temperatures are shown to be considerably increased with Mn values but depressed by a rise in biomagnetic parameter (N _H ) and also a rise in micropolarity (R). Translational velocity components are found to decrease substantially with micropolarity (R), a trend consistent with Newtonian blood flows. Micro-rotation values are shown to increase considerably with a rise in R values but are reduced with a rise in biomagnetic parameter (N _H ). Both translational velocities are boosted with a rise in Darcy number as is micro-rotation. Forchheimer number is also shown to decrease translational velocities but increase micro-rotation. Excellent agreement is demonstrated between both numerical solutions. The mathematical model finds applications in blood flow control devices, hemodynamics in porous biomaterials and also biomagnetic flows in highly perfused skeletal tissue. Dedicated to Professor Y.C. Fung (1919-), Emeritus Professor of Biomechanics, Bioengineering Department, University of California at San Diego, USA for his seminal contributions to biomechanics and physiological fluid mechanics over four decades and his excellent encouragement to the authors, in particular OAB, with computational biofluid dynamics research.     

Mathematical simulation of radial heat transfer in packed beds by pseudohomogeneous modeling

Uniform flow regime and constant effective thermal conductivity inside packed beds are commonly accepted in the evaluation of the fluid dynamics and heat transfer in such systems.However,several authors have confirmed the presence of an oscillatory velocity profile caused by the effective contribution of porosity profile in the fluid dynamic behavior of packed beds,which directly influences the heat transfer inside the beds.This paper describes the application of a pseudo-homogeneous mathematical ...     

Numerical modeling of fluid flow and heat transfer in a narrow Taylor-Couette-Poiseuille system

We consider turbulent flows in a differentially heated Taylor-Couette system with an axial Poiseuille flow. The numerical approach is based on the Reynolds Stress Modeling (RSM) of [Elena and Schiestel, 1996] and [Schiestel and Elena, 1997] widely validated in various rotor-stator cavities with throughflow ( [Poncet, 2005], [Poncet et al., 2005] and [Haddadi and Poncet, 2008]) and heat transfer (Poncet and Schiestel, 2007). To show the capability of the present code, our numerical predictions are compared very favorably to the velocity measurements of Escudier and Gouldson (1995) in the isothermal case, for both the mean and turbulent fields. The RSM model improves, in particular, the predictions of the k-ε model of Naser (1997). Then, the second order model is applied for a large range of rotational Reynolds (3744 ? Re_i ? 37,443) and Prandtl numbers (0.01 ? Pr ? 12), flow rate coefficient (0 ? C_w ? 30,000) in a very narrow cavity of radius ratio s = R_i/R_o = 0.961 and aspect ratio L = (R_o − R_i)/h = 0.013, where R_i and R_o are the radii of the inner and outer cylinders respectively and h is the cavity height. Temperature gradients are imposed between the incoming fluid and the inner and outer cylinders. The mean hydrodynamic and thermal fields reveal three distinct regions across the radial gap with a central region of almost constant axial and tangential mean velocities and constant mean temperature. Turbulence, which is weakly anisotropic, is mainly concentrated in that region and vanishes towards the cylinders. The mean velocity distributions are not clearly affected by the rotational Reynolds number and the flow rate coefficient. The effects of the flow parameters on the thermal field are more noticeable and considered in details. Correlations for the averaged Nusselt numbers along both cylinders are finally provided according to the flow control parameters Re_i, C_w, and Pr.     

Experimental investigation in pool boiling heat transfer of ammonia/water mixture and heat transfer correlations

The nucleate pool boiling heat transfer coefficient of ammonia/water mixture was investigated on a cylindrical heated surface at low pressure of 4-8 bar and at low mass fraction of 0 < x_NH3 < 0.3 and at different heat flux. The effect of mass fraction, heat flux and pressure on boiling heat transfer coefficient was studied. The results indicate that the heat transfer coefficient in the mixture decreases with increase in ammonia mass fraction, increases with increase in heat flux and pressure in the investigated range. The measured heat transfer coefficient was compared with existing correlations. The experimental data were predicted with an accuracy of ±20% by the correlation of Calus&Rice, correlation of Stephan-Koorner and Inoue-Monde correlation for ammonia/water mixture in the investigated range of low ammonia mass fraction. The empirical constant of the first two correlations is modified by fitting the correlation to the present experimental data. The modified Calus&Rice correlation predicts the present experimental data with an accuracy of ±18% and the modified Stephan-Koorner correlation with an accuracy of ±16%.     

Numerical study of flow and heat transfer during a high-speed micro-drop impact on thin liquid films

Numerical simulation of high-speed micro-droplet impingement on thin liquid film covering a heated solid surface has been carried out. Effect of droplet Weber number and liquid film thickness on the characteristics of flow and heat transfer has been investigated using the coupled level set and volume of fluid method. The code is validated against both the experimental and numerical results from the literature. Results show that the crown dynamics is mostly affected by variations in the initial film thickness but is weakly influenced by changes in the Weber number. The liquid within the film can be categorized as three regions based on the heat transfer distribution: the static film region, the transition region, and the impact region. The transient local wall temperature shows three stages: first stage when the temperature decreases rapidly, followed by a second stage in which the temperature starts to rise and then becomes almost constant in the third stage. After drop impact, the local Nusselt number continuously increases until reaching a maximum value, and then decreases approaching the initial impact stage. Our analysis of the change in Weber number shows that larger Weber number contributes to intense temperature variation at the crater core relative to other radial locations. Lastly, the results reveal that the thinner liquid film leads to lower wall temperature and hence, higher average Nusselt number.     

The influence of fluid properties on electrohydrodynamic heat transfer enhancement in liquids under viscous and electrically dominated flow conditions

Fluid property effects on electrohydrodynamic (EHD) heat transfer enhancement were investigated. Heat transfer, pressure drop, electrical power requirements, and the transition between the viscous dominated and electrically dominated flow regimes as a function of fluid properties were examined using three cooling oils having widely varying physical properties. Low viscosity and low electrical conductivity gave the greatest heat transfer enhancement for a given electrical power input. The required electrical power to achieve a specified heat transfer enhancement was greater for working fluids that had a small charge relaxation time, defined as the ratio of the electrical permittivity to the electrical conductivity. These results correlate well with available experimental and analytical data. A theoretical prediction of the effect of fluid properties and forced flow rate on the onset of EHD enhancement was experimentally verified. The onset of significant EHD heat transfer enhancement occurs most readily in low viscosity liquids at low Reynolds number flows for a given electrical power input.     

Use of experiment and an inverse method to study interface heat transfer during solidification in the investment casting process

A technique to determine the thermal boundary conditions existing during the solidification of metallic alloys in the investment casting process is presented. Quantitative information about these conditions is needed so that numerical models of heat transfer in this process produce accurate results. In particular, the variation of the boundary conditions both spatially and temporally must be known. The method used involves the application of a new inverse heat conduction method to thermal data recorded during laboratory experiments of aluminium alloy solidification in investment casting shell moulds. The resultant heat transfer coefficient for the alloy/mould interface is calculated. An experimental programme to determine requisite mould thermal properties was also undertaken. It was observed that there is significant variation of the alloy/mould heat transfer coefficient during solidification. It is found to be highly dependent on the alloy type and on the vertical position below the initial free surface of the liquid metal. The aluminium casting alloys used in this study were 413, A356, 319 (Aluminum Association designations), and commercially pure aluminium. These alloys have significantly different freezing ranges. In particular, it was found that alloys with a high freezing range solidify with rates of heat transfer to the mould which are very sensitive to metallostatic head.     

CFD modeling and experimental validation of heat and mass transfer in wood poles subjected to high temperatures: a conjugate approach

In this article, a coupling method is presented in the case of high thermal treatment of a wood pole and a three-dimensional numerical simulation is proposed. The conservation equations for the wood sample are obtained using diffusion equation with variables diffusion coefficients and the incompressible Reynolds averaged Navier–Stokes equations have been solved for the flow field. The connection between the two problems is achieved by expressing the continuity of the state variables and their respective fluxes through the interface. Turbulence closure is obtained by the use of the standard k–ɛ model with the usual wall function treatment. The model equations are solved numerically by the commercial package ANSYS-CFX10. The wood pole was subjected to high temperature treatment under different operating conditions. The model validation is carried out via a comparison between the predicted values with those obtained experimentally. The comparison of the numerical and experimental results shows good agreement, implying that the proposed numerical algorithm can be used as a useful tool in designing high-temperature wood treatment processes. A parametric study was also carried out to determine the effects of several parameters such as initial moisture content, wood aspect ratio and final gas temperature on temperature and moisture content distributions within the samples during heat treatment.     

A simultaneous PIV and heat transfer study of bubble interaction with free convection flow

Whole field velocity and point temperature and surface heat flux measurements were performed to characterise the interaction of a single rising ellipsoidal air bubble with the free convection flow from a heated flat surface immersed in water at different angles of inclination. Two thermocouples and a hot film sensor were used to characterise heat transfer from the surface, while a time-resolved digital particle image velocimetry technique was used to map the bubble induced flow in a plane parallel to the surface. Heat flux fluctuations, preceding and following the bubble passage, were shown to correlate with the variation in both local flow velocities and fluid temperatures. The largest increases in heat transfer were recorded when both flow and temperature effects combined to enhance the convective cooling simultaneously. Such conditions were shown to be most likely met when the block was inclined at 45°, thus forcing the bubble to slide closer to the heated surface and hence to the thermal boundary layer.     

Lie group analysis of viscoelastic MHD aligned flow and heat transfer

Exact solutions for an incompressible, viscoelastic, electrically conducting MHD aligned fluid are obtained for velocity components and temperature profiles. Lie Group method is applied to obtain the solution and the symmetries used are of translational type.The English text was polished by Keren Wang and Yunming Chen.     

Peristaltic transport of fourth grade fluid with heat transfer and induced magnetic field

We investigate the influence of an induced magnetic field on the peristaltic flow of an incompressible fourth grade fluid in a symmetric channel with heat transfer. Adopting long wavelength, low Reynolds number and small Deborah number assumptions we derive the solutions for stream function, pressure gradient, temperature, magnetic force function, induced magnetic field and current density. Qualitative agreement is demonstrated between the graphs and expected observations.     

Flow and heat transfer performance of a mini-channel radiator with cylinder disturbed flow

Mini-channel heat sinks have relatively low Nusselt number due to small Reynolds number. For heat transfer enhancement purpose, a mini-channel radiator with cylinder disturbed flow was proposed. The disturbed flow was created by a circular cylinder placed horizontally in front of channels entrance. The performance of heat transfer and pressure drop with/without disturbed flow was studied experimentally. It was found that the friction factor of mini-channel flow was larger than that of the macro-channel flow due to larger surface roughness, and the pressure drop caused by cylinder disturbed flow was less than 5%. It also concluded that the average Nusselt number increases with augment of Reynolds and Prandtl number. The Nusselt number correlations as the function of the Reynolds and Prandtl number were given for evaluation the heat removal performance of similar heat radiators. There is an inflexion point in the empirical formulas when the channel length equals to the thermal entrance length. For the mini-channels heat radiators with disturbed flow, the inflexion Reynolds number is larger than that of without disturbed flow. Due to the flow pulsing caused by circular cylinder placed in front of channels entrance, the thermal entrance length increases. On the other hand, for both mini-channels with or without disturbed flow, the thermal resistance increases with the decrease of pressure drop.     

Direct simulation of turbulent flow and heat transfer in a channel. Part I: Smooth walls

Interest in the use of supercomputers for the direct numerical calculation of turbulence prompts the development of efficient numerical techniques so that calculation at higher Reynolds numbers might be made. This paper presents an efficient pseudo-spectral technique, similar to but different from others that have recently appeared, for the calculation of momentum and heat transfer to a constant-property, turbulent fluid in a two-dimensional channel with walls at different, uniform temperature. The code uses no empiricism, although periodic boundary conditions are used for fluctuating quantities in the streamwise and spanwise directions. Calculations were made for a Prandtl number of 0·72 and Reynolds number based on friction velocity and channel half-height of 180 or 2800 based on channel half-height and average velocity. Calculations of mean velocity profile, turbulence intensities, skewness, flatness, Reynolds stress and eddy diffusivity of heat near a wall compare favourably with experimental results. Representative contour plots of the temperature field near the wall and of the spanwise and streamwise two-point velocity correlations are given. Deficiencies are that the calculation requires many hours on a fast computer with a large high-speed memory and that the grid size in each direction for appropriate resolution is approximately proportional to the square of the Reynolds number and to the Prandtl number raised to some power greater than one.     

Flow of a micropolar fluid due to a porous stretching sheet and heat transfer

The present work investigates the micropolar fluid flow due to a permeable stretching sheet and the resulting heat transfer. Unlike the existing numerical works on the flow phenomenon in the literature, the prime interest here is to analytically work out shape of the solutions and identify whether they are unique. Indeed, unique solutions are detected and presented in the exact formulas for the associated boundary layer equations. Temperature field influenced by the microrotation is also mathematically resolved in the cases of constant wall temperature, constant heat flux and Newtonian heating. To discover the salient physical features of many mechanisms acting on the considered problem, it is adequate to have the analytical velocity and temperature fields and also closed-form skin friction/couple stress/heat transfer coefficients, all as given in the current paper. For instance, the practically significant rate of heat transfer is represented by a single formula valid for all three temperature cases.     

One-dimensional finite element heat transfer solution of a fin with triangular perforations of bases parallel and towered its base

Heat transfer dissipation from a horizontal rectangular fin embedded with equilateral triangular perforations is computed numerically using one-dimensional finite element technique. The bases of the triangles are parallel and toward the fin base. The body of the fin is discretized into a number of subdivisions (finite elements). The number of these elements can be altered as required according to the automatic mesh generation. The heat dissipation of the perforated fin is computed and compared with that of the solid one of the same dimensions and same thermal properties. The comparison refers to acceptable results and heat dissipation enhancement due to certain perforation.