scispace - formally typeset
Search or ask a question

Showing papers in "Journal of Heat Transfer-transactions of The Asme in 2007"


Journal ArticleDOI
TL;DR: In this article, the effects of various parameters such as the ratio of the thermal conductivity of nanoparticles to that of a base fluid, volume fraction, nanoparticle size, and temperature on the effective thermal conductivities of nanofluids were investigated.
Abstract: The addition of a small amount of nanoparticles in heat transfer fluids results in the new thermal phenomena of nanofluids (nanoparticle-fluid suspensions) reported in many investigations. However, traditional conductivity theories such as the Maxwell or other macroscale approaches cannot explain the thermal behavior of nanofluids. Recently, Jang and Choi proposed and modeled for the first time the Brownian-motion-induced nanoconvection as a key nanoscale mechanism governing the thermal behavior of nanofluids, but did not clearly explain this and other new concepts used in the model. This paper explains in detail the new concepts and simplifying assumptions and reports the effects of various parameters such as the ratio of the thermal conductivity of nanoparticles to that of a base fluid, volume fraction, nanoparticle size, and temperature on the effective thermal conductivity of nanofluids. Comparison of model predictions with published experimental data shows good agreement for nanofluids containing oxide, metallic, and carbon nanotubes.

399 citations


Journal ArticleDOI
TL;DR: In this article, the thermal conductivity of water and ethylene glycol-based nanofluids containing alumina, zinc-oxide, and titanium-dioxide nanoparticles is measured using the transient hot-wire method.
Abstract: The thermal conductivity of water- and ethylene glycol-based nanofluids containing alumina, zinc-oxide, and titanium-dioxide nanoparticles is measured using the transient hot-wire method. Measurements are performed by varying the particle size and volume fraction, providing a set of consistent experimental data over a wide range of colloidal conditions. Emphasis is placed on the effect of the suspended particle size on the effective thermal conductivity. Also, the effect of laser-pulse irradiation, i.e., the particle size change by laser ablation, is examined for ZnO nanofluids. The results show that the thermal-conductivity enhancement ratio relative to the base fluid increases linearly with decreasing the particle size but no existing empirical or theoretical correlation can explain the behavior. It is also demonstrated that high-power laser irradiation can lead to substantial enhancement in the effective thermal conductivity although only a small fraction of the particles are fragmented.

356 citations


Journal ArticleDOI
TL;DR: In this article, the thermal conductivity of individual single-wall carbon nanotubes has been investigated and the results indicate that thermal conductivities increases with nanotube length, varying from about 10 W/m to 375W/m K depending on the various simulation conditions.
Abstract: Despite the significant amount of research on carbon nanotubes, the thermal conductivity of individual single-wall carbon nanotubes has not been well established. To date only a few groups have reported experimental data for these molecules. Existing molecular dynamics simulation results range from several hundred to 6600 W/m K and existing theoretical predictions range from several dozens to 9500 W/m K. To clarify the several-order-of-magnitude discrepancy in the literature, this paper utilizes molecular dynamics simulation to systematically examine the thermal conductivity of several individual (10, 10) single-wall carbon nanotubes as a function of length, temperature, boundary conditions and molecular dynamics simulation methodology. Nanotube lengths ranging from 5 nm to 40 nm are investigated. The results indicate that thermal conductivity increases with nanotube length, varying from about 10 W/m to 375 W/m K depending on the various simulation conditions. Phonon decay times on the order of hundreds of fs are computed. These times increase linearly with length, indicating ballistic transport in the nanotubes. A simple estimate of speed of sound, which does not require involved calculation of dispersion relations, is presented based on the heat current autocorrelation decay. Agreement with the majority of theoretical/computational literature thermal conductivity data is achieved for the nanotube lengths treated here. Discrepancies in thermal conductivity magnitude with experimental data are primarily attributed to length effects, although simulation methodology, stress, and intermolecular potential may also play a role. Quantum correction of the calculated results reveals thermal conductivity temperature dependence in qualitative agreement with experimental data.

296 citations


Journal ArticleDOI
TL;DR: In this article, the authors present a systematic experimental method of studying the heat transfer behavior of buoyancy-driven nanofluids, where the authors estimate the range of volume fractions that results in maximum thermal enhancement and the impact of volume fraction on Nusselt number.
Abstract: This paper presents a systematic experimental method of studying the heat transfer behavior of buoyancy-driven nanofluids. The presence of nanoparticles in buoyancy-driven flows affects the thermophysical properties of the fluid and consequently alters the rate of heat transfer. The focus of this paper is to estimate the range of volume fractions that results in maximum thermal enhancement and the impact of volume fraction on Nusselt number. The test cell for the nanofluid is a two-dimensional rectangular enclosure with differentially heated vertical walls and adiabatic horizontal walls filled with 27 nm Al 2 O 3 -H 2 O nanofluid. Simulations were performed to measure the transient and steady-state thermal response of nanofluid to imposed isothermal condition. The volume fraction is varied between 0% and 8%. It is observed that the trend of the temporal and spatial evolution of temperature profile for the nanofluid mimics that of the carrier fluid. Hence, the behaviors of both fluids are similar. Results shows that for small volume fraction, 0.2 ≤ O≤2% the presence of the nanoparticles does not impede the free convective heat transfer, rather it augments the rate of heat transfer. However, for large volume fraction O>2%, the convective heat transfer coefficient declines due to reduction in the Rayleigh number caused by increase in kinematic viscosity. Also, an empirical correlation for Nu o as a function of O and Ra has been developed, and it is observed that the nanoparticle enhances heat transfer rate even at a small volume fraction.

207 citations


Journal ArticleDOI
TL;DR: In this paper, an atomistic Green function was developed to simulate phonon transport across a strained germanium (or silicon) thin film between two semi-infinite silicon contacts.
Abstract: An atomistic Green 's function method is developed to simulate phonon transport across a strained germanium (or silicon) thin film between two semi-infinite silicon (or germanium) contacts. A plane-wave formulation is employed to handle the translational symmetry in directions parallel to the Interfaces. The phonon transmission function and thermal conductance across the thin film are evaluated for various atomic configurations. The contributions from lattice straining and material heterogeneity are evaluated separately, and their relative magnitudes are characterized. The dependence of thermal conductance on film thickness is also calculated, verifying that the thermal conductance reaches an asymptotic value for very thick films. The thermal boundary resistance of a single Si/Ge interface is computed and agrees well with analytical model predictions. Multiple-interface effects on thermal resistance are investigated, and the results indicate that the first few interfaces have the most significant effect on the overall thermal resistance.

202 citations


Journal ArticleDOI
TL;DR: In this paper, a stacked microchannel heat sink is fabricated using silicon micromachining techniques and thermal performance is characterized through experimental measurements and numerical simulations, where effects of coolant flow direction, flow rate allocation among layers, and nonuniform heating are studied.
Abstract: One of the promising liquid cooling techniques for microelectronics is attaching a microchannel heat sink to, or directly fabricating microchannels on, the inactive side of the chip. A stacked microchannel heat sink integrates many layers of microchannels and manifold layers into one stack. Compared with single-layered microchannels, stacked microchannels provide larger flow passages, so that for a fixed heat load the required pressure drop is significantly reduced. Better temperature uniformity can be achieved by arranging counterflow in adjacent microchannel layers. The dedicated manifolds help to distribute coolant uniformly to microchannels. In the present work, a stacked microchannel heat sink is fabricated using silicon micromachining techniques. Thermal performance of the stacked microchannel heat sink is characterized through experimental measurements and numerical simulations. Effects of coolant flow direction, flow rate allocation among layers, and nonuniform heating are studied. Wall temperature profiles are measured using an array of nine platinum thin-film resistive temperature detectors deposited simultaneously with thin-film platinum heaters on the backside of the stacked structure. Excellent overall cooling performance (0.09 ° C/W cm 2 ) for the stacked microchannel heat sink has been shown in the experiments. It has also been identified that over the tested flow rate range, counterflow arrangement provides better temperature uniformity, while parallel flow has the best performance in reducing the peak temperature. Conjugate heat transfer effects for stacked microchannels for different flow conditions are investigated through numerical simulations. Based on the results, some general design guidelines for stacked microchannel heat sinks are provided.

182 citations


Journal ArticleDOI
TL;DR: In this paper, a series of experimental investigations were designed to systematically examine the effects of the geometric dimensions (i.e., coating thickness, volumetric porosity, and pore size, as well as the surface conditions of the porous coatings) on the pool-boiling performance and characteristics.
Abstract: To better understand the mechanisms that govern the behavior of pool boiling on horizontal highly conductive microporous coated surfaces, a series of experimental investigations were designed to systematically examine the effects of the geometric dimensions (i.e., coating thickness, volumetric porosity, and pore size, as well as the surface conditions of the porous coatings) on the pool-boiling performance and characteristics. The study was conducted using saturated distilled water at atmospheric pressure (101 kPa) and porous surfaces fabricated from sintered isotropic copper wire screens. For nucleate boiling on the microporous coated surfaces, two vapor ventilation modes were observed to exist: (i) upward and (ii) mainly from sideways leakage to the unsealed sides and partially from the center of porous surfaces. The ratio of the heater size to the coating thickness, the friction factor of the two-phase flow to single-phase flow inside the porous coatings, as well as the input heat flux all govern the vapor ventilation mode that occurs. In this investigation, the ratio of heater size to coating thickness varies from 3.5 to 38 in order to identify the effect of heater size on the boiling characteristics. The experimental results indicate that the boiling performance and characteristics are also strongly dependent on the volumetric porosity and mesh size, as well as the surface conditions when the heater size is given. Descriptions and discussion of the typical boiling characteristics; the progressive boiling process, from pool nucleate boiling to film boiling; and the boiling performance curves on conductive microporous coated surfaces are all systematically presented.

176 citations


Journal ArticleDOI
TL;DR: A localized radial basis function (RBF) meshless method is developed for coupled viscous fluid flow and convective heat transfer problems and is shown to produce accurate results while requiring a much-reduced effort in problem preparation in comparison to other traditional numerical methods.
Abstract: A localized radial basis function (RBF) meshless method is developed for coupled viscous fluid flow and convective heat transfer problems. The method is based on new localized radial-basis function (RBF) expansions using Hardy Multiquadrics for the sought-after unknowns. An efficient set of formulae are derived to compute the RBF interpolation in terms of vector products thus providing a substantial computational savings over traditional meshless methods. Moreover, the approach developed in this paper is applicable to explicit or implicit time marching schemes as well as steady-state iterative methods. We apply the method to viscous fluid flow and conjugate heat transfer (CHT) modeling. The incompressible Navier‐Stokes are time marched using a Helmholtz potential decomposition for the velocity field. When CHT is considered, the same RBF expansion is used to solve the heat conduction problem in the solid regions enforcing temperature and heat flux continuity of the solid/fluid interfaces. The computation is accelerated by distributing the load over several processors via a domain decomposition along with an interface interpolation tailored to pass information through each of the domain interfaces to ensure conservation of field variables and derivatives. Numerical results are presented for several cases including channel flow, flow in a channel with a square step obstruction, and a jet flow into a square cavity. Results are compared with commercial computational fluid dynamics code predictions. The proposed localized meshless method approach is shown to produce accurate results while requiring a much-reduced effort in problem preparation in comparison to other traditional numerical methods. DOI: 10.1115/1.2402181

151 citations


Journal ArticleDOI
TL;DR: In this article, the thermal and hydraulic performances of silicon-based, low aspect ratio micropin-fin cold plates under cross flow conditions are reported, and the experiments were conducted using water, over a Reynolds number (Re) ranging from 40 to 1000, where Re is the Reynolds number based on maximum velocity and the fin diameter.
Abstract: Experimental results of the thermal and hydraulic performances of silicon-based, low aspect ratio micropin-fin cold plates under cross flow conditions are reported. The pins were both circular and square in shape with dimensions (diameter for circular and sides for square) ranging from 50 μm to 150 μm. The test chip contained 20 integral 75 X 75 μm temperature sensors which were used to determine the thermal resistance (K W -1 ) of the cold plates. The experiments were conducted using water, over a Reynolds number (Re) ranging from 40 to 1000. The data show that the average Nusselt number (Nu) based on the fin diameter varies as Re 0.84 for Re 100, where Re is the Reynolds number based on maximum velocity and the fin diameter. Analysis of the Fanning friction factor (f) data shows that f varies as Re -1.35 for Re 100.

140 citations


Journal ArticleDOI
TL;DR: In this paper, the stagnation point flow toward a stretching vertical sheet is investigated and the features of the flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed.
Abstract: The stagnation point flow toward a stretching vertical sheet is investigated in this study. The temperature and velocity of the sheet as well as the velocity of the external flow are assumed to vary in a power law with the distance from the stagnation point. The governing system of equations is first transformed into a dimensionless form, and then the resulting equations are solved numerically by a finite-difference method. The features of the flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed. Both assisting and opposing flows are considered. It is found that, for opposing flow, dual solutions exist in the neighborhood of the separation region, whereas for assisting flow the solution is unique.

138 citations


Journal ArticleDOI
TL;DR: In this paper, a new superposition model is developed to correlate the heat transfer data in the saturated boiling regime in microchannel flows, and specific features of flow boiling in microchannels are incorporated while deriving analytical solutions for the convection enhancement factor and nucleate boiling suppression factor.
Abstract: Flow boiling heat transfer to water in microchannels is experimentally investigated. The dimensions of the microchannels considered are 275 X 636 and 406X1063 μm 2 . The experiments are conducted at inlet water temperatures in the range of 67-95°C and mass fluxes of 221-1283 kg/m 2 s. The maximum heat flux investigated in the tests is 129 W/cm 2 and the maximum exit quality is 0.2. Convective boiling heat transfer coefficients are measured and compared to predictions from existing correlations for larger channels. While an existing correlation was found to provide satisfactory prediction of the heat transfer coefficient in subcooled boiling in microchannels, saturated boiling was not well predicted by the correlations for macrochannels. A new superposition model is developed to correlate the heat transfer data in the saturated boiling regime in microchannel flows. In this model, specific features of flow boiling in microchannels are incorporated while deriving analytical solutions for the convection enhancement factor and nucleate boiling suppression factor. Good agreement with the experimental measurements indicates that this model is suitable for use in analyzing boiling heat transfer in microchannel flows.

Journal ArticleDOI
TL;DR: In this paper, two shell-and-tube heat exchangers (STHXs) using continuous helical baffles instead of segmental baffles were proposed, designed, and tested in this study.
Abstract: Two shell-and-tube heat exchangers (STHXs) using continuous helical baffles instead of segmental baffles used in conventional STHXs were proposed, designed, and tested in this study. The two proposed STHXs have the same tube bundle but different shell configurations. The flow pattern in the shell side of the heat exchanger with continuous helical baffles was forced to be rotational and helical due to the geometry of the continuous helical baffles, which results in a significant increase in heat transfer coefficient per unit pressure drop in the heat exchanger. Properly designed continuous helical baffles can reduce fouling in the shell side and prevent the flow-induced vibration as well. The performance of the proposed STHXs was studied experimentally in this work. The heat transfer coefficient and pressure drop in the new STHXs were compared with those in the STHX with segmental baffles. The results indicate that the use of continuous helical baffles results in nearly 10% increase in heat transfer coefficient compared with that of conventional segmental baffles for the same shell-side pressure drop. Based on the experimental data, the nondimensional correlations for heat transfer coefficient and pressure drop were developed for the proposed continuous helical baffle heat exchangers with different shell configurations, which might be useful for industrial applications and further study of continuous helical baffle heat exchangers. This paper also presents a simple and feasible method to fabricate continuous helical baffles used for STHXs.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the thermal conductivity of SWNT poly(methylmethacrylate) (PMMA) nanocomposites with random SWNT orientations and loading up to 9 vol % using the comparative technique.
Abstract: Single-walled carbon nanotubes (SWNTs) are considered as promising filler materials for improving the thermal conductivity of conventional polymers. We carefully investigated the thermal conductivity of SWNT poly(methylmethacrylate) (PMMA) nanocomposites with random SWNT orientations and loading up to 9 vol % using the comparative technique. The composites were prepared by coagulation and exhibit ∼250% improvement in the thermal conductivity at 9 vol %. The experimental results were analyzed using the versatile Nielsen model, which accounts for many important factors, including filler aspect ratio and maximum packing fraction. In this work, the aspect ratio was determined by atomic force microscopy (AFM) and used as an input parameter in the Nielsen model. We obtained good agreement between our results and the predictions of the Nielsen model, which indicates that higher aspect ratio fillers are needed to achieve further enhancement. Our analysis also suggests that improved thermal contact between the SWNT network and the matrix material would be beneficial.

Journal ArticleDOI
TL;DR: In this paper, a new design concept is presented to increase the adiabatic effectiveness of film cooling from a row of film-cooling holes by modifying the geometry upstream of the holes to modify the approaching boundary-layer flow and its interaction with the film cooling jets.
Abstract: A new design concept is presented to increase the adiabatic effectiveness of film cooling from a row of film-cooling holes. Instead of shaping the geometry of each hole; placing tabs, struts, or vortex generators in each hole; or creating a trench about a row of holes, this study proposes a geometry modification upstream of the holes to modify the approaching boundary-layer flow and its interaction with the film-cooling jets. Computations, based on the ensemble-averaged Navier-Stokes equations closed by the realizable k-e turbulence model, were used to examine the usefulness of making the surface just upstream of a row of film-cooling holes into a ramp with a backward-facing step. The effects of the following parameters were investigated: angle of the ramp (8.5 deg, 10 deg, 14 deg), distance between the backward-facing step and the row of film-cooling holes (0.5D,D), blowing ratio (0.36, 0.49, 0.56, 0.98), and "sharpness" of the ramp at the corners. Results obtained show that an upstream ramp with a backward-facing step can greatly increase surface adiabatic effectiveness. The laterally averaged adiabatic effectiveness with a ramp can be two or more times higher than without the ramp by increasing upstream and lateral spreading of the coolant.

Journal ArticleDOI
TL;DR: In this paper, the local heat transfer coefficients induced by piezoelectric fans are determined experimentally for a fan vibrating close to an electrically heated stainless steel foil, and the entire temperature field is observed by means of an infrared camera.
Abstract: Piezoelectric fans have been shown to provide substantial enhancements in heat transfer over natural convection while consuming very little power. These devices consist of a piezoelectric material attached to a flexible cantilever beam. When driven at resonance, large oscillations at the cantilever tip cause fluid motion, which in turn results in improved heat transfer rates. In this study, the local heat transfer coefficients induced by piezoelectric fans are determined experimentally for a fan vibrating close to an electrically heated stainless steel foil, and the entire temperature field is observed by means of an infrared camera. Four vibration amplitudes ranging from 6.35 to 10 mm are considered, with the distance from the heat source to the fan tip chosen to vary from 0.01 to 2.0 times the amplitude. The two-dimensional contours of the local heat transfer coefficient transition from a lobed shape at small gaps to an almost circular shape at intermediate gaps. At larger gaps, the heat transfer coefficient distribution becomes elliptical in shape. Correlations developed with appropriate Reynolds and Nusselt number definitions describe the area-averaged thermal performance with a maximum error of less than 12%.

Journal ArticleDOI
TL;DR: In this article, the potential of winglet-type vortex generator arrays for multirow inline-tube heat exchangers with an emphasis on providing fundamental understanding of the relation between local flow behavior and heat transfer enhancement mechanisms was evaluated.
Abstract: This work is aimed at assessing the potential of winglet-type vortex generator (VG) "arrays" for multirow inline-tube heat exchangers with an emphasis on providing fundamental understanding of the relation between local flow behavior and heat transfer enhancement mechanisms. Three different winglet configurations in common-flow-up arrangement are analyzed in the seven-row compact fin-and-tube heat exchanger: (a) single-VG pair; (b) a 3VG-inline array (alternating tube row); and (c) a 3VG-staggered array. The numerical study involves three-dimensional time-dependent modeling of unsteady laminar flow (330≤ Re < 850) and conjugate heat transfer in the computational domain, which is set up to model the entire fin length in the air flow direction. It was found that the impingement of winglet redirected flow on the downstream tube is an important heat transfer augmentation mechanism for the common-flow-up arrangement of vortex generators in the inline-tube geometry. At Re=850 with a constant tube-wall temperature, the 3VG-inline-array configuration achieves enhancements up to 32% in total heat flux and 74% in j factor over the baseline case, with an associated pressure-drop increase of about 41%. The numerical results for the integral heat transfer quantities agree well with the available experimental measurements.

Journal ArticleDOI
Ahmad Fakheri1
TL;DR: The concept of heat exchanger efficiency was introduced in this article, which is defined as the ratio of the heat transferred in the actual heat exchange to the heat that would be transferred in an ideal heat exchange.
Abstract: This paper provides the solution to the problem of defining thermal efficiency for heat exchangers based on the second law of thermodynamics. It is shown that corresponding to each actual heat exchanger, there is an ideal heat exchanger that is a balanced counter-flow heat exchanger. The ideal heat exchanger has the same UA, the same arithmetic mean temperature difference, and the same cold to hot fluid inlet temperature ratio. The ideal heat exchanger's heat capacity rates are equal to the minimum heat capacity rate of the actual heat exchanger. The ideal heat exchanger transfers the maximum amount of heat, equal to the product of UA and arithmetic mean temperature difference, and generates the minimum amount of entropy, making it the most efficient and least irreversible heat exchanger. The heat exchanger efficiency is defined as the ratio of the heat transferred in the actual heat exchanger to the heat that would be transferred in the ideal heat exchanger. The concept of heat exchanger efficiency provides a new way for the design and analysis of heat exchangers and heat exchanger networks.

Journal ArticleDOI
TL;DR: In this article, critical heat flux (CHF) of R-123 in a silicon-based microchannel heat sink was investigated at exit pressures ranging from 227 kPa to 520 kPa.
Abstract: Critical heat flux (CHF) of R-123 in a silicon-based microchannel heat sink was investigated at exit pressures ranging from 227 kPa to 520 kPa. Critical heat flux data were obtained over effective heat fluxes ranging from 53 W/cm 2 to 196 W/cm 2 and mass fluxes from 291 kg/m 2 s to 1118 kg/m 2 s. Flow images and high exit qualities suggest that dryout is the leading CHF mechanism. The effect of mass velocity, exit quality, and system pressure were also examined, and a new correlation is presented to represent the effect of these parameters.

Journal ArticleDOI
TL;DR: In this paper, a critical review is made of recent studies of the modeling of viscous dissipation in a saturated porous medium, with applications to either forced convection or natural convection.
Abstract: A critical review is made of recent studies of the modeling of viscous dissipation in a saturated porous medium, with applications to either forced convection or natural convection. Alternative forms of the viscous dissipation function are discussed. Limitations to the concept of fully developed convection are noted. Special attention is focused on the roles of viscous dissipation and work done by pressure forces (flow work) in natural convection in a two-dimensional box with either lateral or bottom heating.

Journal ArticleDOI
TL;DR: In this paper, the authors used mathematical models to systematically investigate the transport phenomena leading to the formation of porosity and to find possible solutions to reduce or eliminate porosity formation in laser welding.
Abstract: Porosity has been frequently observed in solidified, deep penetration pulsed laser welds. Porosity is detrimental to weld quality. Our previous study shows that porosity formation in laser welding is associated with the weld pool dynamics, keyhole collapse, and solidification processes. The objective of this paper is to use mathematical models to systematically investigate the transport phenomena leading to the formation of porosity and to find possible solutions to reduce or eliminate porosity formation in laser welding. The results indicate that the formation of porosity in pulsed laser welding is caused by two competing factors: one is the solidification rate of the molten metal and the other is the backfilling speed of the molten metal during the keyhole collapse process. Porosity will be formed in the final weld if the solidification rate of the molten metal exceeds the backfilling speed of liquid metal during the keyhole collapse and solidification processes. Porosity formation was found to be strongly related with the depth-to-width aspect ratio of the keyhole. The larger the ratio, the easier porosity will be formed, and the larger the size of the voids. Based on these studies, controlling the laser pulse profile is proposed to prevent/eliminate porosity formation in laser welding. Its effectiveness and limitations are demonstrated in the current studies. The model predictions are qualitatively consistent with reported experimental results.

Journal ArticleDOI
TL;DR: The basics of the finite-volume method are reviewed, and its extension to unstructured meshes through the use of cell-based and control-volume finite-element schemes is traced.
Abstract: Pressure-based finite-volume techniques have emerged as the methods of choice for a wide variety of industrial applications involving incompressible fluid flow. In this paper, we trace the evolution of this class of solution techniques. We review the basics of the finite-volume method, and trace its extension to unstructured meshes through the use of cell-based and control-volume finite-element schemes. A critical component of the solution of incompressible flows is the issue of pressure-velocity storage and coupling. The development of staggered-mesh schemes and segregated solution techniques such as the SIMPLE algorithm are reviewed. Co-located storage schemes, which seek to replace staggered-mesh approaches, are presented. Coupled multigrid schemes, which promise to replace segregated-solution approaches, are discussed. Extensions of pressure-based techniques to compressible flows are presented. Finally, the shortcomings of existing techniques and directions for future research are discussed.

Journal ArticleDOI
TL;DR: In this paper, an exact shape function finite element model is developed that envisions the foam as randomly oriented cylinders in cross-o w with an axially varying coolant temperature.
Abstract: Convective heat transfer in aluminum metal foam sandwich panels is investigated with potential applications to actively cooled thermal protection systems in hypersonic and reentry vehicles. The size eects of the metal foam core are experimentally investigated and the eects of foam thickness on convective transfer are established. Four metal foam specimens are utilized with a relative density of 0.08 and pore density of 20 ppi in a range of thickness from 6.4 mm to 25.4 mm in increments of approximately 6 mm. An exact-shapefunction nite element model is developed that envisions the foam as randomly oriented cylinders in cross o w with an axially varying coolant temperature eld. Our experimental results indicate that larger foam thicknesses produce increased heat transfer levels in metal foams. Initial FE simulations using a fully developed, turbulent velocity prole show the potential of this numerical tool to model convective heat transfer in metal foams. Metal foam sandwich panels have been proposed as alternative multi-functional materials for structural thermal protection systems in hypersonic and re-entry vehicles 1 . 2 This type of construction oers numerous advantages over other actively cooled concepts because of the unique properties of metal foams. These materials, when brazed between metallic face sheets, are readily suited to allow coolant passage without the addition of alien components that may compromise structural performance. Moreover, the mechanical properties can be varied to suit dieren t structural needs by varying the foam relative density. From a heat transfer point of view, these materials have been shown to be exceptional heat exchangers primarily due to the increased surface area available for heat transfer between the solid and uid phases. The thermo-mechanical response of metal foam sandwich panels has been recently studied and characterized. 2 In particular, it has been shown that using air as coolant at sucien tly high velocities, the strain due to buckling of these structures under thermo-mechanical loads can be virtually eliminated. The implementation of these materials in thermal protection systems, however, requires that a proper heat transfer model exists that allows the coupling between the thermo-mechanical and heat transfer problems to be properly analyzed. In other words, it is necessary to understand how dieren t foam properties such as relative density, pore density, and foam thickness will aect the heat loads that this type of structural component can remove. Heat transfer in metal foams has been a subject of active research in recent years. Lu et al. 3 developed an analytical model envisioning the foam as an array of mutually perpendicular cylinders subjected to cross-o w. In this study, a closed-form expression for the convective coecien t of a foam-lled channel with constant wall temperatures was presented based on foam geometry and material and uid properties. These authors reported that the simplifying assumptions used in their analysis were likely to lead to an over-prediction of the actual heat transfer level. This model has been partially validated by Bastawros and Evans 4 who performed forced convection experiments on aluminum foams adhered to silicon substrate face sheets. These authors reported that the predictions of Lu et al. 3 regarding the dependence of the convective coecien t on coolant velocity and strut diameter were qualitatively consistent with their observations, but that the foam thickness eects were not adequately modeled. In particular, they reported that the heat dissipation rate

Journal ArticleDOI
TL;DR: In this paper, a second-order accurate finite volume code with nonstaggered arrangement of variables is developed employing momentum interpolation for the pressure-velocity coupling for two-dimensional laminar flow and heat transfer past a triangular cylinder placed in a horizontal channel for the range 80≤Re≤200 and blockage ratio 1/12≤β≤1/3.
Abstract: Computational study of two-dimensional laminar flow and heat transfer past a triangular cylinder placed in a horizontal channel is presented for the range 80≤Re≤200 and blockage ratio 1/12≤β≤1/3. A second-order accurate finite volume code with nonstaggered arrangement of variables is developed employing momentum interpolation for the pressure-velocity coupling. Global mode of cross-stream velocity oscillations predict the first bifurcation point increases linearly with blockage ratio with no second bifurcation found in the range of Re studied. The Strouhal number and rms of lift coefficient increase significantly with blockage ratio and Reynolds number while overall Nusselt number remains almost unchanged for different blockage ratios. At lower blockage ratios, flow is found to be similar to the unconfined flow and is more prone to wake instability. Instantaneous streak lines provide an excellent means of visualizing the von Karman vortex street.

Journal ArticleDOI
TL;DR: In this paper, the role of viscous dissipation in liquids flowing through heated microchannels was analyzed by using the conventional theory, and it was shown that the problem of heat transfer enhancement in microdevices cannot be solved by indefinitely reducing the microchannel dimensions because the viscous effects shall offset the gains of high heat transfer coefficients associated with a reduction in the channel size.
Abstract: Many experimental works appeared in the last decade in the open literature dealing with forced convection through microchannels. The earliest experimental results on single-phase flows in microchannels evidenced that for channels having a hydraulic diameter less than 1mm the conventional continuum models can no longer be considered as able to accurately predict pressure drop and convective heat transfer coefficients. This conclusion seemed to be valid for both gas and liquid flows. Sometimes the authors justified this conclusion by invoking new micro-effects, e.g., electrostatic interaction between the fluid and the walls or scaling effects (axial heat conduction, viscous forces, conjugate heat transfer, wall roughness, and so on). In this paper the role of the viscous dissipation in liquids flowing through heated microchannels will be analyzed by using the conventional theory. We will present a correlation between the Brinkman number and the Nusselt number for silicon ⟨100⟩ and ⟨110⟩ microchannels. It will be demonstrated that the fluid is of importance in establishing the exact limit of significance of the viscous dissipation in microchannels; a criterion to analyze the significance of the viscous effects will be presented. The role of the cross-section aspect ratio on the viscous dissipation will be highlighted. The main goal of this work is to demonstrate that the problem of heat transfer enhancement in microdevices cannot be solved by indefinitely reducing the microchannel dimensions because the viscous dissipation effects shall offset the gains of high heat transfer coefficients associated with a reduction in the channel size.

Journal ArticleDOI
TL;DR: In this paper, a single-phase Si microchannel coolers were designed and characterized in single chip modules in a laboratory environment using either water at 22°C or a fluorinated fluid at temperatures between 20 and -40°C as the coolant.
Abstract: High performance single-phase Si microchannel coolers have been designed and characterized in single chip modules in a laboratory environment using either water at 22°C or a fluorinated fluid at temperatures between 20 and -40° C as the coolant. Compared to our previous work, key performance improvements were achieved through reduced channel pitch (from 75 to 60 microns), thinned channel bases (from 425 to 200 microns of Si), improved thermal interface materials, and a thinned thermal test chip (from 725 to 400 microns of Si). With multiple heat exchanger zones and 60 micron pitch microchannels with a water flow rate of 1.25 1pm, an average unit thermal resistance of 15.9 C mm 2 /W between the chip surface and the inlet cooling water was demonstrated for a Si microchannel cooler attached to a chip with Ag epoxy. Replacing the Ag epoxy layer with an In solder layer reduced the unit thermal resistance to 12.0 C mm 2 /W. Using a fluorinated fluid with an inlet temperature of -30°C and 60 micron pitch microchannels with an Ag epoxy thermal interface layer, the average unit thermal resistance was 25.6 C mm 2 /W. This fell to 22.6 C mm 2 /W with an In thermal interface layer. Cooling >500 W/cm 2 was demonstrated with water. Using a fluorinated fluid with an inlet temperature of -30° C, a chip with a power density of 270 W/cm 2 was cooled to an average chip surface temperature of 35°C. Results using both water and a fluorinated fluid are presented for a range of Si microchannel designs with a channel pitch from 60 to 100 microns.

Journal ArticleDOI
TL;DR: In this article, a multilayer structure consisting of a SiC film coated atop a dielectric photonic crystal (PC) was proposed for coherent emission from both p and s polarizations.
Abstract: Enhancement of thermal emission and control of its direction are important for applications in optoelectronics and energy conversion. A number of structures have been proposed as coherent emission sources, which exhibit a large emissivity peak within a narrow wavelength band and at a well-defined direction. A commonly used structure is the grating, in which the excited surface polaritons or surface waves are coupled with propagating waves in air, resulting in coherent emission for p polarization only. One-dimensional photonic crystals can also support surface waves and may be modified to construct coherent emission sources. The present study investigates coherent emission from a multilayer structure consisting of a SiC film coated atop a dielectric photonic crystal (PC). By exciting surface waves at the interface between SiC and the PC, coherent emission is predicted for both p and s polarizations. In addition to the excitation of surface waves, the emission from the proposed multilayer structure can be greatly enhanced by the cavity resonance mode and the Brewster mode.

Journal ArticleDOI
TL;DR: In this article, a review of phase-change heat transfer in microchannels is presented, including flow boiling and flow condensation, with applications to microchannel heat sinks and microheat exchangers, as well as bubble growth and collapse on microheaters under pulse heating.
Abstract: Recent work on miscroscale phase-change heat transfer, including flow boiling and flow condensation in microchannnels (with applications to microchannel heat sinks and microheat exchangers) as well as bubble growth and collapse on microheaters under pulse heating (with applications to micropumps and thermal inkjet printerheads), is reviewed

Journal ArticleDOI
TL;DR: In this paper, a numerical and experimental study of heat transfer and fluid flow in a single pass counter flow plate heat exchanger with chevron plates has been presented, where the complete geometry of the heat transfer surface and more realistic hydrodynamic and thermal boundary conditions were taken into account.
Abstract: A numerical and experimental study of heat transfer and fluid flow in a single pass counter flow plate heat exchanger with chevron plates has been presented in this paper. CFD analysis of small sized plate heat exchanger was carried out by taking the complete geometry of the heat transfer surface and more realistic hydrodynamic and thermal boundary conditions. A cold channel with two chevron plates and two halves of hot channels on either side having flat periodic boundaries was selected as the computational domain. The numerical model was validated with data from experiments and empirical correlations from literature. Heat transfer and pressure drop data were obtained experimentally with water as the working fluid, in the Reynolds number range 400-1300 and the Prandtl number range 4.4-6.3.

Journal ArticleDOI
TL;DR: In this article, the effects of the mixed convection parameter, Gr/Re, the Knudsen number, Kn, and the ratio of wall temperature difference, r T, on the microchannel hydrodynamic and thermal behaviors are determined.
Abstract: In this study, fully developed mixed convective heat transfer of a Newtonian fluid in an open-ended vertical parallel plate microchannel is analytically investigated by taking the velocity slip and the temperature jump at the wall into account. The effects of the mixed convection parameter, Gr/Re, the Knudsen number, Kn, and the ratio of wall temperature difference, r T , on the microchannel hydrodynamic and thermal behaviors are determined. Finally, a Nu=f(Gr/Re,Kn,r T ) expression is developed.

Journal ArticleDOI
TL;DR: In this paper, the effect of various modeling schemes on simulation results is investigated, such as the effects of turbulence models, dispersed phase modeling, inclusion of different forces (Saffman, thermophoresis, and Brownian), trajectory tracking, and mist injection scheme.
Abstract: Numerical simulation is performed in this study to explore film-cooling enhancement by injecting mist into the cooling air with a focus on investigating the effect of various modeling schemes on simulation results. The effect of turbulence models, dispersed-phase modeling, inclusion of different forces (Saffman, thermophoresis, and Brownian), trajectory tracking, and mist injection scheme is studied. The effect of flow inlet boundary conditions (with/without air supply plenum), inlet turbulence intensity, and the near-wall grid density on simulation results is also included. Simulation of a two-dimensional (2D) slot film cooling with a fixed blowing angle and blowing ratio shows a 2% mist (by mass) injected into the cooling air can increase the cooling effectiveness about 45%. The renormalization group (RNG) k-e model, Reynolds stress model, and the standard k-e turbulence model with an enhanced wall treatment produce consistent and reasonable results while the turbulence dispersion has a significant effect on mist film cooling through the stochastic trajectory calculation. The thermophoretic force slightly increases the cooling effectiveness, but the effect of Brownian force and Saffman lift is imperceptible. The cooling performance deteriorates when the plenum is included in the calculation due to the altered velocity profile and turbulence intensity at the jet exit plane. The results of this paper can provide guidance for corresponding experiments and serve as the qualification reference for future more complicated studies with 3D cooling holes, different blowing ratios, various density ratios, and rotational effect.