0

# Current Issue

### Research Papers: Bio-Heat and Mass Transfer

J. Heat Transfer 136, 101101 (2014) (12 pages);   Paper No: HT-13-1496;   doi:10.1115/1.4027864

The oscillatory rotational motion of the elephant pinna is considered a key mechanism in metabolic heat dissipation. Limited experimental investigations have revealed that the flapping of the elephant's pinna is responsible for surface heat transfer enhancement. The objective of the present experimental and computational work is to investigate the physics of the flow induced by the pinna's motion and its effects on the heat transfer. This was accomplished by designing, fabricating and testing two full-size laboratory models of elephant pinnae: one rigid and one flexible, both instrumented with small size thermocouples for time-dependent surface temperature measurements. A constant heat flux is applied to both sides of each model which is rotated about a fixed edge with a frequency of 2 rad/s in an infinite domain at ambient conditions. Of interest is the study of the impact of the flexural strength of the model's material on surface heat transfer. Additional computer simulations of the flow and thermal fields revealed a hooked-shape vortex tube around the free edges of the flapping pinna. This result is confirmed by the flow visualization with smoke particles. Both experimental and computational results exhibit local surface temperature profiles characterized by a transient unsteady periodic variation followed by a steady periodic phase. Flow visualization indicated significant interaction between the vortical structures shed off the edge and the flexible model's boundary layer. It has been found that the cooling of the flexible model is enhanced by 30%.

### Research Papers: Conduction

J. Heat Transfer 136, 101301 (2014) (8 pages);   Paper No: HT-13-1438;   doi:10.1115/1.4027882

There are many applications for problems involving thermal conduction in two-dimensional (2D) cylindrical objects. Experiments involving thermal parameter estimation are a prime example, including cylindrical objects suddenly placed in hot or cold environments. In a parameter estimation application, the direct solution must be run iteratively in order to obtain convergence with the measured temperature history by changing the thermal parameters. For this reason, commercial conduction codes are often inconvenient to use. It is often practical to generate numerical solutions for such a test, but verification of custom-made numerical solutions is important in order to assure accuracy. The present work involves the generation of an exact solution using Green's functions where the principle of superposition is employed in combining a one-dimensional (1D) cylindrical case with a 1D Cartesian case to provide a temperature solution for a 2D cylindrical. Green's functions are employed in this solution in order to simplify the process, taking advantage of the modular nature of these superimposed components. The exact solutions involve infinite series of Bessel functions and trigonometric functions but these series sometimes converge using only a few terms. Eigenvalues must be determined using Bessel functions and trigonometric functions. The accuracy of the solutions generated using these series is extremely high, being verifiable to eight or ten significant digits. Two examples of the solutions are shown as part of this work for a family of thermal parameters. The first case involves a uniform initial condition and homogeneous convective boundary conditions on all of the surfaces of the cylinder. The second case involves a nonhomogeneous convective boundary condition on a part of one of the planar faces of the cylinder and homogeneous convective boundary conditions elsewhere with zero initial conditions.

### Research Papers: Evaporation, Boiling, and Condensation

J. Heat Transfer 136, 101501 (2014) (9 pages);   Paper No: HT-13-1346;   doi:10.1115/1.4027881

This study describes the effect of an electric field on nucleate boiling and critical heat flux (CHF) in pool boiling. A dielectric liquid of AE-3000 was used as the working fluid. A heating surface was polished to a surface roughness of 0.05 μm. A microsized electrode, in which slits were provided, was designed to generate a nonuniform electric field and produce electrohydrodynamic (EHD) effects with the application of high dc voltages. The obtained results confirmed CHF enhancement as the EHD effects increased CHF to 86.2 W/cm2 with a voltage of −3000 V, which was four times greater than pool boiling in the absence of the electrode. The usual traveling wave on the bubble interface, induced by the Kelvin–Helmholtz instability, was modified by adding the EHD effects. The traveling wave model exhibits the essential features of the phenomenon and shows good agreement with the experimental data.

J. Heat Transfer 136, 101502 (2014) (9 pages);   Paper No: HT-13-1510;   doi:10.1115/1.4027966

A study was conducted to experimentally characterize the pool boiling performance of hydrofluorocarbon HFC-245fa at pressures ranging from 0.15 MPa to 1.1 MPa (reduced pressure range: 0.04–0.31). Pool boiling experiments were conducted using horizontally oriented 1-cm2 heated surfaces to quantify the effects of pressure and a microporous-enhanced coating on heat transfer coefficients and critical heat flux (CHF) values. Results showed that the coating enhanced heat transfer coefficients and CHF by 430% and 50%, respectively. The boiling heat transfer performance of HFC-245fa was then compared with the boiling performance of HFC-134a and hydrofluoroolefin HFO-1234yf.

### Research Papers: Forced Convection

J. Heat Transfer 136, 101701 (2014) (8 pages);   Paper No: HT-13-1225;   doi:10.1115/1.4027855

Turbulent flow and heat transfer of air with variable properties in a set of regular polygonal ducts and circular tube have been numerically simulated. All the ducts have the same hydraulic diameter as their characteristic lengths in the Reynolds number. The flow is modeled as three-dimensional (3D) and fully elliptic by using the finite volume method and the standard k-ε turbulence model. The results showed that the relatively strong secondary flow could be observed with variable properties fluid. For the regular polygonal ducts, the local heat transfer coefficient along circumferential direction is not uniform; there is an appreciable reduction in the corner region and the smaller the angle of the corner region, the more appreciable deterioration the corner region causes. The use of hydraulic diameter for regular polygonal ducts leads to unacceptably large errors in turbulent heat transfer determined from the circular tube correlations. Based on the simulation results, a correction factor is proposed to predict turbulent heat transfer in regular polygonal ducts.

J. Heat Transfer 136, 101702 (2014) (10 pages);   Paper No: HT-13-1009;   doi:10.1115/1.4027908

This paper presents a large eddy simulation (LES) of flow and heat transfer in a tandem configuration of two square cylinders at moderate Reynolds number ($Re=16,000$). Compressible LES on a hybrid mesh is used to predict the flow structure and the heat transfer at the wall. The goals of this work are to analyze the flow and the heat transfer around a tandem arrangement of two inline square cylinders as well as to propose a LES approach that can be applied to convective heat transfer problems in industrial configurations. The meshing strategy allows to resolve the flow field until the viscous sublayer with y+ of the order unity. The wall adapting linear eddy model is chosen to model the subgrid turbulent viscosity. Aerodynamics results are validated versus experimental measurements performed on isolated cylinders and on tandem configurations. The main flow structures responsible for heat transfer are analyzed. Finally, heat transfer around both cylinders of the tandem is described.

### Research Papers: Heat Exchangers

J. Heat Transfer 136, 101801 (2014) (10 pages);   Paper No: HT-13-1483;   doi:10.1115/1.4027911

Heat removal capacity, coolant pumping power requirement, and surface temperature nonuniformity are three major challenges facing single-phase flow microchannel compact heat exchangers. In this paper multi-objective optimization has been performed to increase heat removal capacity, and decrease pumping power and temperature nonuniformity in complex networks of microchannels. Three-dimensional (3D) four-floor configurations of counterflow branching networks of microchannels were optimized to increase heat removal capacity from surrounding silicon substrate (15 × 15 × 2 mm). Each floor has four different branching subnetworks with opposite flow direction with respect to the next one. Each branching subnetwork has four inlets and one outlet. Branching patterns of each of these subnetworks could be different from the others. Quasi-3D conjugate heat transfer analysis has been performed by developing a software package which uses quasi-1D thermofluid analysis and a 3D steady heat conduction analysis. These two solvers were coupled through their common boundaries representing surfaces of the cooling microchannels. Using quasi-3D conjugate analysis was found to require one order of magnitude less computing time than a fully 3D conjugate heat transfer analysis while offering comparable accuracy for these types of application. The analysis package is capable of generating 3D branching networks with random topologies. Multi-objective optimization using modeFRONTIER software was performed using response surface approximation and genetic algorithm. Diameters and branching pattern of each subnetwork and coolant flow direction on each floor were design variables of multi-objective optimization. Maximizing heat removal capacity, while minimizing coolant pumping power requirement and temperature nonuniformity on the hot surface, were three simultaneous objectives of the optimization. Pareto-optimal solutions demonstrate that thermal loads of up to 500 W/cm2 can be managed with four-floor microchannel cooling networks. A fully 3D thermofluid analysis was performed for one of the optimal designs to confirm the accuracy of results obtained by the quasi-3D simulation package used in this paper.

### Research Papers: Micro/Nanoscale Heat Transfer

J. Heat Transfer 136, 102401 (2014) (5 pages);   Paper No: HT-13-1369;   doi:10.1115/1.4027910

We performed molecular dynamics (MD) simulations of the interface which is comprised of self-assembled monolayer (SAM) and water solvent to investigate heat transfer characteristics. In particular, local thermal boundary conductance (TBC), which is an inverse of so-called Kapitza resistance, at the SAM–solvent interface was evaluated by using the nonequilibrium MD (NEMD) technique in which the one-dimensional thermal energy flux was imposed across the interface. By using two kinds of SAM terminal with hydrophobic and hydrophilic properties, the local TBCs of these interfaces with water solvent were evaluated, and the result showed a critical difference due to an affinity between SAM and solvent. In order to elucidate the molecular-scale mechanism that makes this difference, microscopic components contributing to thermal energy flux across the interface of hydrophilic SAM and water were evaluated in detail, i.e., the total thermal energy flux is decomposed into the heat transfer modes such as the contribution of molecular transport and that of energy exchange by molecular interactions. These heat transfer modes were also compared with those in the bulk water.

### Technical Brief

J. Heat Transfer 136, 104501 (2014) (6 pages);   Paper No: HT-12-1230;   doi:10.1115/1.4024592

In this paper, the problem of steady two-dimensional magnetohydrodynamic (MHD) stagnation-point flow and heat transfer of an incompressible viscous fluid over a stretching/shrinking sheet is investigated in the presence of velocity and thermal slips. With the help of similarity transformations, the governing Navier–Stokes and the energy equations are reduced to ordinary differential equations, which are then solved numerically using a shooting technique. Interesting solution behavior is observed for the similarity equations with multiple solution branches for certain parameter domain. Fluid velocity increases due to the increasing value of the velocity slip parameter resulting in a decrease in the temperature field. Temperature at a point increases with increase in the thermal slip parameter. The effects of the slips, stretching/shrinking, and the magnetic parameters on the skin friction or the wall shear stress, heat flux from the surface of the sheet, velocity, and temperature profiles are computed and discussed.