Improving the Performance of the Micro Jet Array Air Cooling MEMS Device Using CFD

Author(s):  
Yangki Jung ◽  
R. Panneer Selvam

In this work ways to improve the performance of the MJA is investigated using computer modeling. Finite-difference method is used to analyze and understand the heat transfer mechanism in the MEMS based air micro-jet array (MJA) impinging cooling device. The three-dimensional Navier-Stokes (NS) equations using the compressible flow are solved. The computed temperature distribution at the bottom of the MJA is in good agreement with the experimental measurement data. To improve the heat transfer performance, an optimum model is developed. This enhances the heat removal capacity about 100% (0.5 W/cm2°C). The details of computer modeling and the flow visualization to understand the flow and heat transfer in the MJA are presented.   This paper was also originally published as part of the Proceedings of the ASME 2005 Heat Transfer Summer Conference.

Author(s):  
Andojo Ongkodjojo ◽  
Alexis R. Abramson ◽  
Norman C. Tien

This work is a continuation of previous investigations aimed at developing an innovative microfabricated air-cooling technology that employs an electrohydrodynamic corona discharge (i.e. ionic wind pump) [1], [2]. This technology enables the miniaturization of cooling systems for next generation electronics. Our single ionic wind pump element consists of two parallel collecting electrodes between which a single emitting tip is positioned. Two-dimensional (2-D) and three-dimensional (3-D) simulations using COMSOL Multiphysics™ are additionally employed to predict the temperature distribution, the flow field, and the heat removal capacity of the device in operation. One such model utilizes a small gap between collector and emitter electrodes and demonstrates an improvement in the COP (coefficient of performance) of a single device. Comparisons are made with experimental temperature data on an actual device. The purpose of this work is therefore to optimize the performance of a single microfabricated ionic wind pump to enable the development of an array of these elements for use in larger-scale heat transfer applications.


Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.


Author(s):  
Anil K. Tolpadi ◽  
James A. Tallman ◽  
Lamyaa El-Gabry

Conventional heat transfer design methods for turbine airfoils use 2-D boundary layer codes (BLC) combined with empiricism. While such methods may be applicable in the mid span of an airfoil, they would not be very accurate near the end-walls and airfoil tip where the flow is very three-dimensional (3-D) and complex. In order to obtain accurate heat transfer predictions along the entire span of a turbine airfoil, 3-D computational fluid dynamics (CFD) must be used. This paper describes the development of a CFD based design system to make heat transfer predictions. A 3-D, compressible, Reynolds-averaged Navier-Stokes CFD solver with k-ω turbulence modeling was used. A wall integration approach was used for boundary layer prediction. First, the numerical approach was validated against a series of fundamental airfoil cases with available data. The comparisons were very favorable. Subsequently, it was applied to a real engine airfoil at typical design conditions. A discussion of the features of the airfoil heat transfer distribution is included.


Author(s):  
Vijay K. Garg

A multi-block, three-dimensional Navier-Stokes code has been used to compute heat transfer coefficient on the blade, hub and shroud for a rotating high-pressure turbine blade with 172 film-cooling holes in eight rows. Film cooling effectiveness is also computed on the adiabatic blade. Wilcox’s k-ω model is used for modeling the turbulence. Of the eight rows of holes, three are staggered on the shower-head with compound-angled holes. With so many holes on the blade it was somewhat of a challenge to get a good quality grid on and around the blade and in the tip clearance region. The final multi-block grid consists of 4784 elementary blocks which were merged into 276 super blocks. The viscous grid has over 2.2 million cells. Each hole exit, in its true oval shape, has 80 cells within it so that coolant velocity, temperature, k and ω distributions can be specified at these hole exits. It is found that for the given parameters, heat transfer coefficient on the cooled, isothermal blade is highest in the leading edge region and in the tip region. Also, the effectiveness over the cooled, adiabatic blade is the lowest in these regions. Results for an uncooled blade are also shown, providing a direct comparison with those for the cooled blade. Also, the heat transfer coefficient is much higher on the shroud as compared to that on the hub for both the cooled and the uncooled cases.


2000 ◽  
Author(s):  
M. Singh ◽  
P. K. Panigrahi ◽  
G. Biswas

Abstract A numerical study of rib augmented cooling of turbine blades is reported in this paper. The time-dependent velocity field around a pair of symmetrically placed ribs on the walls of a three-dimensional rectangular channel was studied by use of a modified version of Marker-And-Cell algorithm to solve the unsteady incompressible Navier-Stokes and energy equations. The flow structures are presented with the help of instantaneous velocity vector and vorticity fields, FFT and time averaged and rms values of components of velocity. The spanwise averaged Nusselt number is found to increase at the locations of reattachment. The numerical results are compared with available numerical and experimental results. The presence of ribs leads to complex flow fields with regions of flow separation before and after the ribs. Each interruption in the flow field due to the surface mounted rib enables the velocity distribution to be more homogeneous and a new boundary layer starts developing downstream of the rib. The heat transfer is primarily enhanced due to the decrease in the thermal resistance owing to the thinner boundary layers on the interrupted surfaces. Another reason for heat transfer enhancement can be attributed to the mixing induced by large-scale structures present downstream of the separation point.


Author(s):  
Alejandro Herna´ndez Rossette ◽  
Zdzislaw Mazur C. ◽  
Jesu´s Cordero Guridi ◽  
Eric Chumacero Polanco

As a gas turbine entry temperature (TET) increases, thermal loading on first stage blades increases too and therefore, a variety of cooling techniques and thermal barrier coatings (TBCs) are used to maintain the blade temperature within the acceptable limits. In this work a multi-block three dimensional Navier-Stokes commercial turbomachinery oriented CFD-code has been used to compute steady state conjugated heat transfer (CHT) on the blade suction and pressure coated sides of a rotating first inter-stage (nozzle and bucket) with cooling holes of a 60 MW Gas turbine. A Spallart Allmaras model was used for modeling the turbulence. Convection and radiation were modeled for a super alloy blade with and without TBC. The CFD simulations were configured with a mesh domain of nozzle and bucket inter-stage in order to predict the fluid parameters at inlet and outlet of bucket for validate with turbine inter-stage parameter data test of gas turbine manufacturer. The effects of blade surface temperature changes were simulated with both configurations coated and uncoated blades.


Author(s):  
Yeon-Sik Kim ◽  
Yoon-Sub Sim ◽  
Eui-Kwang Kim

A new design concept for a decay heat removal system in a liquid metal reactor is proposed. The new design utilizes a thermosyphon to enhance the heat removal capacity and its heat transfer characteristics are analyzed against the current PSDRS (Passive Safety Decay heat Removal System) in the KALIMER (Korea Advanced LIquid MEtal Reactor) design. The preliminary analysis results show that the new design with a thermosyphon yields substantial increase of 20∼40% in the decay heat removal capacity compared to the current design that do not have the thermosyphon. The new design reduces the temperature rise in the cooling air of the system and helps the surrounding structure in maintaining its mechanical integrity for long term operation at an accident. Also the analysis revealed the characteristics of the interactions among various heat transfer modes in the new design.


Metals ◽  
2020 ◽  
Vol 10 (8) ◽  
pp. 1065
Author(s):  
Chunyu Liu ◽  
Xiaodong Li ◽  
Run Luo ◽  
Rafael Macian-Juan

The Small Modular Dual Fluid Reactor (SMDFR) is a novel molten salt reactor based on the dual fluid reactor concept, which employs molten salt as fuel and liquid lead/lead-bismuth eutectic (LBE) as coolant. A unique design of this reactor is the distribution zone, which locates under the core and joins the core region with the inlet pipes of molten salt and coolant. Since the distribution zone has a major influence on the heat removal capacity in the core region, the thermal hydraulics characteristics of the distribution zone have to be investigated. This paper focuses on the thermal hydraulics analysis of the distribution zone, which is conducted by the numerical simulation using COMSOL Multiphysics with the CFD (Computational Fluid Dynamics) module and the Heat Transfer module. The energy loss and heat exchange in the distribution zone are also quantitatively analyzed. The velocity and temperature distributions of both molten salt and coolant at the outlet of the distribution zone, as inlet of the core region, are produced. It can be observed that the outlet velocity profiles are proportional in magnitude to the inlet velocity ones with a similar shape. In addition, the results show that the heat transfer in the center region is enhanced due to the velocity distribution, which could compensate the power peak and flatten the temperature distribution for a higher power density.


2004 ◽  
Vol 126 (4) ◽  
pp. 519-523 ◽  
Author(s):  
Siddharth Bhopte ◽  
Musa S. Alshuqairi ◽  
Dereje Agonafer ◽  
Gamal Refai-Ahmed

The current numerical investigation will examine the effect of an impinging mixed convection air jet on the heat transfer rate of a parallel flat plate heat sink. A three-dimensional numerical model was developed to evaluate the effects of the nozzle diameter d, nozzle-to-target vertical placement H/d, Rayleigh number, and the jet Reynolds number on the heat transfer rates from a discrete heat source. Simulations were performed for a Prandtl number of 0.7 and for Reynolds numbers ranging from 100 to 5000. The governing equations were solved in the dimensionless form using a commercial finite-volume package. Average Nusselt numbers were obtained, at H/d=3 and two jet diameters, for the bare heat source, for the heat source with a base heat sink, and for the heat source with the finned heat sink. The heat transfer rates from the bare heat source surface have been compared with the ones obtained with the heat sink in order to determine the overall performance of the heat sink in an impingement configuration.


2016 ◽  
Vol 138 (9) ◽  
Author(s):  
Amin Bouraoui ◽  
Rachid Bessaïh

In this paper, a numerical study of three-dimensional (3D) natural convection air-cooling of two identical heat sources, simulating electronic components, mounted in a rectangular enclosure was carried out. The governing equations were solved by using the finite-volume method based on the SIMPLER algorithm. The effects of Rayleigh number Ra, spacing between heat sources d, and aspect ratios Ax in x-direction (horizontal) and Az in z-direction (transversal) of the enclosure on heat transfer were investigated. In steady state, when d is increased, the heat transfer is more important than when the aspect ratios Ax and Az are reduced. In oscillatory state, the critical Rayleigh numbers Racr for different values of spacing between heat sources and their aspect ratios, at which the flow becomes time dependent, were obtained. Results show a strong relation between heat transfers, buoyant flow, and boundary layer. In addition, the heat transfer is more important at the edge of each face of heat sources than at the center region.


Sign in / Sign up

Export Citation Format

Share Document