Effects of Various Substrate Openings for Electronic Cooling under Forced and Natural Convection

This study experimentally investigates the heat transfer effects of forced convection and natural convection under different substrate openings design. A computational fluid dynamics (CFD) model was established and implemented to verify and explain the experimental results and heat transfer behavior. It is found that different opening position will destroy the growth of the boundary layer on substrates to alter the cooling ability for both forced under low Reynolds number and natural convection. Nevertheless, having too many opening may reduce heat conduction and affect the overall heat transfer performance. This study provides future researchers with a guideline on designing and electronic package manufacturing.

Simulation of the Flow in a Packed-Bed with and without a Static Mixer by Using CFD Technique

The major focus of this work was to characterize hydrodynamics in a packed-bed with and without static mixer by using Computational Fluid Dynamic (CFD). The commercial software: COMSOL MULTIPHYSICSTM Version 3.3 was used to simulate flow fields of mixed-gas reactants i.e. CO and H2. The packed-bed was a single tube with the inside diameter of 0.8 cm and the length of 1.2 cm. The static mixer was inserted inside the tube. The number of twisting elements was 1 with 0.8 cm in diameter and 1.2 cm in length. The packed-bed with and without static mixer were both packed with approximately 700 spherical structures representing catalyst pellets. Incompressible Navier-Stokes equations were used to model the gas flow inside the beds at steady state condition, in which the inlet Reynolds Number (Re) was 2.31. The results revealed that, with the insertion of static mixer, the gas was forced to flow radially inward and outward between the central portion of the tube and the tube wall. This could help improving the overall performance of the packed-bed, which could be utilized for heterogeneous catalytic reaction such as reforming and Fischer- Tropsch reactions.

Fluid Flow and Heat Transfer Structures of Oscillating Pipe Flows

The RANS method with Saffman-s turbulence model was employed to solve the time-dependent turbulent Navier-Stokes and energy equations for oscillating pipe flows. The method of partial sums of the Fourier series is used to analyze the harmonic velocity and temperature results. The complete structures of the oscillating pipe flows and the averaged Nusselt numbers on the tube wall are provided by numerical simulation over wide ranges of ReA and ReR. Present numerical code is validated by comparing the laminar flow results to analytic solutions and turbulence flow results to published experimental data at lower and higher Reynolds numbers respectively. The effects of ReA and ReR on the velocity, temperature and Nusselt number distributions have been di scussed. The enhancement of the heat transfer due to oscillating flows has also been presented. By the way of analyzing the overall Nusselt number over wide ranges of the Reynolds number Re and Keulegan- Carpenter number KC, the optimal ratio of the tube diameter over the oscillation amplitude is obtained based on the existence of a nearly constant optimal KC number. The potential application of the present results in sea water cooling has also been discussed.

Numerical Simulation of Convection Heat Transfer in a Lid-Driven Cavity with an Open Side

In this manuscript, the LBM is applied for simulating of Mixed Convection in a Lid-Driven cavity with an open side. The cavity horizontal walls are insulated while the west Lid-driven wall is maintained at a uniform temperature higher than the ambient. Prandtl number (Pr) is fixed to 0.71 (air) while Reynolds number (Re) , Richardson number (Ri) and aspect ratio (A) of the cavity are changed in the range of 50-150 , of 0.1-10 and of 1-4 , respectively. The numerical code is validated for the standard square cavity, and then the results of an open ended cavity are presented. Result shows by increasing of aspect ratio, the average Nusselt number (Nu) on lid- driven wall decreases and with same Reynolds number (Re) by increasing of aspect ratio (A), Richardson number plays more important role in heat transfer rate.

An Experimental Study of the Effect of Coil Step on Heat Transfer Coefficient in Shell- Side of Shell-and-Coil Heat Exchanger

In this study the mixed convection heat transfer in a coil-in-shell heat exchanger for various Reynolds numbers and various dimensionless coil pitch was experimentally investigated. The experiments were conducted for both laminar and turbulent flow inside coil and the effects of coil pitch on shell-side heat transfer coefficient of the heat exchanger were studied. The particular difference in this study in comparison with the other similar studies was the boundary conditions for the helical coils. The results indicate that with the increase of coil pitch, shell-side heat transfer coefficient is increased.

Turbulent Forced Convection Flow in a Channel over Periodic Grooves Using Nanofluids

Turbulent forced convection flow in a 2-dimensional channel over periodic grooves is numerically investigated. Finite volume method is used to study the effect of turbulence model. The range of Reynolds number varied from 10000 to 30000 for the ribheight to channel-height ratio (B/H) of 2. The downstream wall is heated by a uniform heat flux while the upstream wall is insulated. The investigation is analyzed with different types of nanoparticles such as SiO2, Al2O3, and ZnO, with water as a base fluid are used. The volume fraction is varied from 1% to 4% and the nanoparticle diameter is utilized between 20nm to 50nm. The results revealed 114% heat transfer enhancement compared to the water in a grooved channel by using SiO2 nanoparticle with volume fraction and nanoparticle diameter of 4% and 20nm respectively.

Grid Independence Study of Flow Past a Square Cylinder Using the Multi-Relaxation-Time Lattice Boltzmann Method

Numerical calculations of flow around a square cylinder are presented using the multi-relaxation-time lattice Boltzmann method at Reynolds number 150. The effects of upstream locations, downstream locations and blockage are investigated systematically. A detail analysis are given in terms of time-trace analysis of drag and lift coefficients, power spectra analysis of lift coefficient, vorticity contours visualizations and phase diagrams. A number of physical quantities mean drag coefficient, drag coefficient, Strouhal number and root-mean-square values of drag and lift coefficients are calculated and compared with the well resolved experimental data and numerical results available in open literature. The results had shown that the upstream, downstream and height of the computational domain are at least 7.5, 37.5 and 12 diameters of the cylinder, respectively.

Numerical Investigation of the Effect of Flow and Heat Transfer of a Semi-Cylindrical Obstacle Located in a Channel

In this study, a semi-cylinder obstacle placed in a channel is handled to determine the effect of flow and heat transfer around the obstacle. Both faces of the semi-cylinder are used in the numerical analysis. First, the front face of the semi-cylinder is stated perpendicular to flow, than the rear face is placed. The study is carried out numerically, by using commercial software ANSYS 11.0. The well-known κ-ε model is applied as the turbulence model. Reynolds number is in the range of 104 to 105 and air is assumed as the flowing fluid. The results showed that, heat transfer increased approximately 15 % in the front faze case, while it enhanced up to 28 % in the rear face case.

Numerical Study of Flow around Flat Tube between Parallel Walls

Flow around a flat tube is studied numerically. Reynolds number is defined base on equivalent circular tube and it is varied in range of 100 to 300. Equations are solved by using finite volume method and results are presented in form of drag and lift coefficient. Results show that drag coefficient of flat tube is up to 66% lower than circular tube with equivalent diameter. In addition, by increasing l/D from 1 to 2, the drag coefficient of flat tube is decreased about 14-27%.

Heat Transfer and Frictional Characteristics in Rectangular Channel with Inclined Perforated Baffles

A numerical study on the turbulent flow and heat transfer characteristics in the rectangular channel with different types of baffles is carried out. The inclined baffles have the width of 19.8 cm, the square diamond type hole having one side length of 2.55 cm, and the inclination angle of 5o. Reynolds number is varied between 23,000 and 57,000. The SST turbulence model is applied in the calculation. The validity of the numerical results is examined by the experimental data. The numerical results of the flow field depict that the flow patterns around the different baffle type are entirely different and these significantly affect the local heat transfer characteristics. The heat transfer and friction factor characteristics are significantly affected by the perforation density of the baffle plate. It is found that the heat transfer enhancement of baffle type II (3 hole baffle) has the best values.

Vortex Shedding at the End of Parallel-plate Thermoacoustic Stack in the Oscillatory Flow Conditions

This paper investigates vortex shedding processes occurring at the end of a stack of parallel plates, due to an oscillating flow induced by an acoustic standing wave within an acoustic resonator. Here, Particle Image Velocimetry (PIV) is used to quantify the vortex shedding processes within an acoustic cycle phase-by-phase, in particular during the “ejection" of the fluid out of the stack. Standard hot-wire anemometry measurement is also applied to detect the velocity fluctuations near the end of the stack. Combination of these two measurement techniques allowed a detailed analysis of the vortex shedding phenomena. The results obtained show that, as the Reynolds number varies (by varying the plate thickness and drive ratio), different flow patterns of vortex shedding are observed by the PIV measurement. On the other hand, the time-dependent hot-wire measurements allow obtaining detailed frequency spectra of the velocity signal, used for calculating characteristic Strouhal numbers. The impact of the plate thickness and the Reynolds number on the vortex shedding pattern has been discussed. Furthermore, a detailed map of the relationship between the Strouhal number and Reynolds number has been obtained and discussed.

Vortex Shedding on Combined Bodies at Incidence to a Uniform Air Stream

Vortex-shedding phenomenon of the flow around combined two bodies having various geometries and sizes has been investigated experimentally in the Reynolds number range between 4.1x103 and 1.75x104. To see the effect of the rotation of the bodies on the vortex shedding, the combined bodies were rotated from 0° to 180°. The combined models have a cross section composing of a main circular cylinder and an attached circular or square cylinder. Results have shown that Strouhal numbers for two cases were changed considerably with the angle of incidence, while it was found to be largely independent of Reynolds number at 150. Characteristics of the vortex formation region and location of flow attachments, reattachments, and separations were observed by means of the flow visualizations. Depending on the inclination angle the effects of flow attachment, separation and reattachment on vortex-shedding phenomenon have been discussed.

Bifurcations and Chaotic Solutions of Two-dimensional Zonal Jet Flow on a Rotating Sphere

We study bifurcation structure of the zonal jet flow the streamfunction of which is expressed by a single spherical harmonics on a rotating sphere. In the non-rotating case, we find that a steady traveling wave solution arises from the zonal jet flow through Hopf bifurcation. As the Reynolds number increases, several traveling solutions arise only through the pitchfork bifurcations and at high Reynolds number the bifurcating solutions become Hopf unstable. In the rotating case, on the other hand, under the stabilizing effect of rotation, as the absolute value of rotation rate increases, the number of the bifurcating solutions arising from the zonal jet flow decreases monotonically. We also carry out time integration to study unsteady solutions at high Reynolds number and find that in the non-rotating case the unsteady solutions are chaotic, while not in the rotating cases calculated. This result reflects the general tendency that the rotation stabilizes nonlinear solutions of Navier-Stokes equations.

Experimental Investigation of the Transient Cooling Characteristics of an Industrial Glass Tempering Unit

Energy consumption rate during the cooling process of industrial glass tempering process is considerably high. In this experimental study the effect of dimensionless jet to jet distance (S/D) and jet to plate distance (H/D) on the cooling time have been investigated. In the experiments 4 mm thick glass samples have been used. Cooling unit consists of 16 mutually placed seamless aluminum nozzles of 8 mm in diameter and 80 mm in length. Nozzles were in staggered arrangement. Experiments were conducted with circular jets for H/D values between 1 and 10, and for S/D values between 2 and 10. During the experiments Reynolds number has been kept constant at 30000. Experimental results showed that the longest cooling time with 87 seconds has been observed in the experiments for S/D=10 and H/D=10 values, while the shortest cooling time with 42.5 seconds has been measured in the experiments for S/D=2 and H/D=4 values.

Heat Transfer from Two Cam Shaped Cylinders in Side-by-Side Arrangement

Heat transfer from two cam shape cylinder in side-byside arrangement had been studied numerically. The transverse gap between the centers of cylinders (T) is allowed to vary to change the pitch ratio (T/Deq). The equivalent diameter of the cylinder (Deq) is 27.6 mm and pitch ratio varies in range of 1≤T/Deq≤3. The Reynolds numbers based on equivalent circular cylinder are within 50≤ Reeq≤150. Results show that Nusselt number of cylinders increases about 1 to 36 percent when pitch ratio increases from 1 to 3.

Mixed Convection in a Vertical Heated Channel: Influence of the Aspect Ratio

In mechanical and environmental engineering, mixed convection is a frequently encountered thermal fluid phenomenon which exists in atmospheric environment, urban canopy flows, ocean currents, gas turbines, heat exchangers, and computer chip cooling systems etc... . This paper deals with a numerical investigation of mixed convection in a vertical heated channel. This flow results from the mixing of the up-going fluid along walls of the channel with the one issued from a flat nozzle located in its entry section. The fluiddynamic and heat-transfer characteristics of vented vertical channels are investigated for constant heat-flux boundary conditions, a Rayleigh number equal to 2.57 1010, for two jet Reynolds number Re=3 103 and 2104 and the aspect ratio in the 8-20 range. The system of governing equations is solved with a finite volumes method and an implicit scheme. The obtained results show that the turbulence and the jet-wall interaction activate the heat transfer, as does the drive of ambient air by the jet. For low Reynolds number Re=3 103, the increase of the aspect Ratio enhances the heat transfer of about 3%, however; for Re=2 104, the heat transfer enhancement is of about 12%. The numerical velocity, pressure and temperature fields are post-processed to compute the quantities of engineering interest such as the induced mass flow rate, and average Nusselt number, in terms of Rayleigh, Reynolds numbers and dimensionless geometric parameters are presented.

Steady State Transpiration Cooling System in Ni-Cr Open-Cellular Porous Plate

The steady-state temperature for one-dimensional transpiration cooling system has been conducted experimentally and numerically to investigate the heat transfer characteristics of combined convection and radiation. The Nickel –Chrome (Ni-Cr) open-cellular porous material having porosity of 0.93 and pores per inch (PPI) of 21.5 was examined. The upper surface of porous plate was heated by the heat flux of incoming radiation varying from 7.7 - 16.6 kW/m2 whereas air injection velocity fed into the lower surface was varied from 0.36 - 1.27 m/s, and was then rearranged as Reynolds number (Re). For the report of the results in the present study, two efficiencies including of temperature and conversion efficiency were presented. Temperature efficiency indicating how close the mean temperature of a porous heat plate to that of inlet air, and increased rapidly with the air injection velocity (Re). It was then saturated and had a constant value at Re higher than 10. The conversion efficiency, which was regarded as the ability of porous material in transferring energy by convection after absorbed from heat radiation, decreased with increasing of the heat flux and air injection velocity. In addition, it was then asymptotic to a constant value at the Re higher than 10. The numerical predictions also agreed with experimental data very well.

Force Statistics and Wake Structure Mechanism of Flow around a Square Cylinder at Low Reynolds Numbers

Numerical investigation of flow around a square cylinder are presented using the multi-relaxation-time lattice Boltzmann methods at different Reynolds numbers. A detail analysis are given in terms of time-trace analysis of drag and lift coefficients, power spectra analysis of lift coefficient, vorticity contours visualizations, streamlines and phase diagrams. A number of physical quantities mean drag coefficient, drag coefficient, Strouhal number and root-mean-square values of drag and lift coefficients are calculated and compared with the well resolved experimental data and numerical results available in open literature. The Reynolds numbers affected the physical quantities.

Natural and Mixed Convection Heat Transfer Cooling of Discrete Heat Sources Placed Near the Bottom on a PCB

Steady state experiments have been conducted for natural and mixed convection heat transfer, from five different sized protruding discrete heat sources, placed at the bottom position on a PCB and mounted on a vertical channel. The characteristic length ( Lh ) of heat sources vary from 0.005 to 0.011 m. The study has been done for different range of Reynolds number and modified Grashof number. From the experiment, the surface temperature distribution and the Nusselt number of discrete heat sources have been obtained and the effects of Reynold number and Richardson number on them have been discussed. The objective is to find the rate of heat dissipation from heat sources, by placing them at the bottom position on a PCB and to compare both modes of cooling of heat sources.

Numerical Simulation on Heat Transfer Enhancement in Channel by Triangular Ribs

Turbulent heat transfer to fluid flow through channel with triangular ribs of different angles are presented in this paper. Ansys 14 ICEM and Ansys 14 Fluent are used for meshing process and solving Navier stokes equations respectively. In this investigation three angles of triangular ribs with the range of Reynolds number varied from 20000 to 60000 at constant surface temperature are considered. The results show that the Nusselt number increases with the increase of Reynolds number for all cases at constant surface temperature. According to the profile of local Nusselt number on ribs walled of channel, the peak is at the midpoint between the two ribs. The maximum value of average Nusselt number is obtained for triangular ribs of angel 60°and at Reynolds number of 60000 compared to the Nusselt number for the ribs of angel 90° and 45° and at same Reynolds number. The recirculation regions generated by the ribs corresponding to the velocity streamline show the largest recirculation region at triangular ribs of angle 60° which also provides the highest enhancement of heat transfer.