Simulation of Micro-Scale Jet Impingement Heat Transfer

2000 ◽  
Author(s):  
Paul A. Boeschoten ◽  
Deborah V. Pence ◽  
James A. Liburdy

Abstract The heat transfer performance of a micro-scale, axisymmetric, confined jet impinging on a flat surface at high Mach numbers (0.2 to 0.6) and low Reynolds numbers (419 to 1310) was computationally studied. The flow is characterized by Knudsen numbers, based on the jet radius, large enough (0.0013) to warrant slip-flow boundary conditions at the impinging surface. The effects of Mach number, compressibility, and slip-flow on heat transfer results are presented, along with the local Nusselt number distributions, and velocity and temperature fields near the impingement surface. Results for uniform wall heat flux show that the wall temperature decreases with increasing Mach number, with a local minimum at r/D = 0.7. The slip velocity also increases with Mach number with peak values also near r/D = 0.7. The resulting Nusselt number increases with increasing Mach number, and a local maximum in the Nusselt number is observed at r/D = 0.6, not at the centerline. In general, compressibility improves heat transfer due to increased fluid density near the impinging surface. Also, inclusion of slip-velocity increases the rate of heat transfer. However, the accompanying temperature-jump condition at the wall is found to reduce the local heat transfer rate. The net effect of the slip-flow boundary conditions applied in this study was an overall reduction in heat transfer.

2003 ◽  
Vol 125 (3) ◽  
pp. 447-453 ◽  
Author(s):  
Deborah V. Pence ◽  
Paul A. Boeschoten ◽  
James A. Liburdy

A computational study is presented of the heat transfer performance of a micro-scale, axisymmetric, confined jet impinging on a flat surface with an embedded uniform heat flux disk. The jet flow occurs at large, subsonic Mach numbers (0.2 to 0.8) and low Reynolds numbers (419 to 1782) at two impingement distances. The flow is characterized by a Knudsen number of 0.01, based on the viscous boundary layer thickness, which is large enough to warrant consideration of slip-flow boundary conditions along the impingement surface. The effects of Mach number, compressibility, and slip-flow on heat transfer are presented. The local Nusselt number distributions are shown along with the velocity, pressure, density and temperature fields near the impingement surface. Results show that the wall temperature decreases with increasing Mach number, M, exhibiting a minimum local value at r/R=1.6 for the highest M. The slip velocity also increases with M, showing peak values near r/R=1.4 for all M. The resulting Nusselt number increases with increasing M, and local maxima are observed near r/R=1.20, rather than at the centerline. In general, compressibility improves heat transfer due to increased fluid density near the impinging surface. The inclusion of slip-velocity and the accompanying wall temperature jump increases the predicted rate of heat transfer by as much as 8–10% for M between 0.4 and 0.8.


2015 ◽  
Vol 138 (4) ◽  
Author(s):  
Pamela Vocale ◽  
Gian Luca Morini ◽  
Marco Spiga

In this work, hydrodynamically and thermally fully developed gas flow through elliptical microchannels is numerically investigated. The Navier–Stokes and energy equations are solved by considering the first-order slip flow boundary conditions and by assuming that the wall heat flux is uniform in the axial direction, and the wall temperature is uniform in the peripheral direction (i.e., H1 boundary conditions). To take into account the microfabrication of the elliptical microchannels, different heated perimeter lengths are analyzed along the microchannel wetted perimeter. The influence of the cross section geometry on the convective heat transfer coefficient is also investigated by considering the most common values of the elliptic aspect ratio, from a practical point of view. The numerical results put in evidence that the Nusselt number is a decreasing function of the Knudsen number for all the considered configurations. On the contrary, the role of the cross section geometry in the convective heat transfer depends on the thermal boundary condition and on the rarefaction degree. With the aim to provide a useful tool for the designer, a correlation that allows evaluating the Nusselt number for any value of aspect ratio and for different working gases is proposed.


Author(s):  
H.P. Rani ◽  
V Leela ◽  
Pulla Nagabhushanam ◽  
R Gangadhara Reddy

Abstract The heat transfer characteristics of mixed convective two-phase flow in an inclined rotating micro-porous channel kept in a transverse magnetic field are investigated numerically. The counterflow arrangement is assumed within the channel. Slip velocity and asymmetric thermal boundary conditions are assumed. The governing energy equation involves the local thermal non-equilibrium (LTNE) between the two phases. The LTNE implications of control parameters on the flow field variables and the average Nusselt number, Nu, are highlighted and pertinent observations are documented. When confined to a few specific cases, the current results are consistent with previous research work. The effect of inclination angle on fluid velocity is determined by the wall temperature difference ratio. According to the findings, for certain values of the wall temperature differential ratio, the velocity increases with the angle, however it takes on a dual character for other values. The Nusselt number (Nu) is expected to increase with the Biot number, Hartmann number, and rotation parameter, while Nu decreases as the Knudsen number increases. The results show that as the wall temperature ratio increases, the Nu converges to a common minimum value. The database was generated from the validated CFD model covering a range of control parameters arising in the system. The multilayer perceptron (MLP) networks were trained using this CFD dataset to predict Nu. The average relative error in Nu's prediction is found to be ±2%.


Author(s):  
Nian Xiao ◽  
John Elsnab ◽  
Susan Thomas ◽  
Tim Ameel

Two analytical models are presented in which the continuum momentum and energy equations, coupled with second-order slip flow and temperature jump boundary conditions, are solved. An isothermal boundary condition is applied to a microchannel with a circular cross section. The flow is assumed to be hydrodynamically fully developed and thermal field is either fully developed or thermally developing from the tube entrance. A traditional first-order slip boundary condition is found to over predict the slip velocity compared to the second-order model. Heat transfer increases at the upper limit of the slip regime for the second-order model. The maximum second-order correction to the first-order Nusselt number is on the order of 18% for air. The second-order effect is also more significant in the entrance region of the tube. The Nusselt number decreases relative to the no-slip value when slip and temperature jump effects are of the same order or when temperature jump effects dominate. When temperature jump effects are small, the Nusselt number increases relative to the no-slip value. Comparisons to a previously reported model for an isoflux boundary condition indicate that the Nusselt number for the isoflux boundary condition exceeds that for the isothermal case at all axial locations.


Author(s):  
Arman Sadeghi ◽  
Abolhassan Asgarshamsi ◽  
Mohammad Hassan Saidi

Fluid flow and heat transfer at microscale have attracted an important research interest in recent years due to the rapid development of microelectromechanical systems (MEMS). Fluid flow in microdevices has some characteristics which one of them is rarefaction effect related with gas flow. In this research, hydrodynamically and thermally fully developed laminar rarefied gas flow in annular microducts is studied using slip flow boundary conditions. Two different cases of the thermal boundary conditions are considered, namely: uniform temperature at the outer wall and adiabatic inner wall (Case A) and uniform temperature at the inner wall and adiabatic outer wall (Case B). Using the previously obtained velocity distribution, energy conservation equation subjected to relevant boundary conditions is numerically solved using fourth order Runge-Kutta method. The Nusselt number values are presented in graphical form as well as tabular form. It is realized that for the case A increasing aspect ratio results in increasing the Nusselt number, while the opposite is true for the case B. The effect of aspect ratio on Nusselt number is more notable at smaller values of Knudsen number, while its effect becomes slighter at large Knudsen numbers. Also increasing Knudsen number leads to smaller values of Nusselt number for the both cases.


Author(s):  
X. Li ◽  
J. L. Gaddis ◽  
T. Wang

The flow field of a 2-D laminar confined impinging slot jet is investigated. Numerical results indicate that there exist two different solutions in some range of geometric and flow parameters. The two steady flow patterns are obtained under identical boundary conditions but only with different initial flow fields. Three different exit boundary conditions are investigated to eliminate artificial effects. The different flow patterns are observed to significantly affect the heat transfer. A flow visualization experiment is carried out to verify the computational results and both flow patterns are observed. The bifurcation mechanism is interpreted and discussed.


Author(s):  
Karthik Krishna ◽  
Mark Ricklick

Ceramic Matrix Composite is a woven material characterized by a significant level of surface waviness of 35–60μm and surface roughness of 5–6μm. To be implemented in a future gas turbine engine they will be cooled traditionally to increase power and efficiency. To analyze the CMC surface effects on heat transfer rate, an impinging circular jet on a simulated CMC surface is studied experimentally and the CMC surface is represented by a high resolution CNC machined surface. The test parameters are jet to plate distance of 7 jet diameters, oblique impingement angles of 45° and 90° and Reynolds numbers of 11,000 to 35,000. The test surface is broken down into constant temperature segments, and individual segment Nusselt number is determined and plotted for the various impingement cases studied. Area-Averaged results show negligible changes in average Nusselt number as compared to the hydrodynamically smooth surface. The impact of the CMC surface feature is negligible compared to the uncertainty in heat transfer coefficient, and therefore traditional design tools can be utilized.


2020 ◽  
Vol 24 (Suppl. 1) ◽  
pp. 13-23
Author(s):  
Zia Ullah ◽  
Muammad Ashraf ◽  
Saqib Zia ◽  
Ishtiaq Ali

The present phenomena address the slip velocity effects on mixed convection flow of electrically conducting fluid with surface temperature and free stream velocity oscillation over a non-conducting horizontal cylinder. To remove the difficulties in illustrating the coupled PDE, the primitive variable formulation for finite dif?ference technique is proposed to transform dimensionless equations into primitive form. The numerical simulations of coupled non-dimensional equations are exam?ined in terms of fluid slip velocity, temperature, and magnetic velocity which are used to calculate the oscillating components of skin friction, heat transfer, and cur?rent density for various emerging parameters magnetic force parameter, ?, mixed convection parameter, ?, magnetic Prandtl number, ?, Prandtl number, and slip factor, SL. It is observed that the effect of slip flow on the non-conducting cylinder is reduced the fluid motion. A minimum oscillating behavior is noted in skin friction at each position but maximum amplitude of oscillation in heat transfer is observed at each position ? = ?/4 and 2?/3. It is further noticed that a fluid velocity increas?es sharply with the impact of slip factor on the fluid-flow mechanism. Moreover, due to frictional forces with lower magnitude between viscous layers, the rise in Prandtl number leads to decrease in skin fiction and heat transfer which is physi?cally in good agreement.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Anuj Kumar Shukla ◽  
Anupam Dewan

Purpose Convective heat transfer features of a turbulent slot jet impingement are comprehensively studied using two different computational approaches, namely, URANS (unsteady Reynolds-averaged Navier–Stokes equations) and SAS (scale-adaptive simulation). Turbulent slot jet impingement heat transfer is used where a considerable heat transfer enhancement is required, and computationally, it is a quite challenging flow configuration. Design/methodology/approach Customized OpenFOAM 4.1, an open-access computational fluid dynamics (CFD) code, is used for SAS (SST-SAS k-ω) and URANS (standard k-ε and SST k-ω) computations. A low-Re version of the standard k-ε model is used, and other models are formulated for good wall-refined calculations. Three turbulence models are formulated in OpenFOAM 4.1 with second-order accurate discretization schemes. Findings It is observed that the profiles of the streamwise turbulence are under-predicted at all the streamwise locations by SST k-ω and SST SAS k-ω models, but follow similar trends as in the reported results. The standard k-ε model shows improvements in the predictions of the streamwise turbulence and mean streamwise velocity profiles in the zone of outer wall jet. Computed profiles of Nusselt number by SST k-ω and SST-SAS k-ω models are nearly identical and match well with the reported experimental results. However, the standard k-ε model does not provide a reasonable profile or quantification of the local Nusselt number. Originality/value Hybrid turbulence model is suitable for efficient CFD computations for the complex flow problems. This paper deals with a detailed comparison of the SAS model with URANS and LES for the first time in the literature. A thorough assessment of the computations is performed against the results reported using experimental and large eddy simulations techniques followed by a detailed discussion on flow physics. The present results are beneficial for scientists working with hybrid turbulence models and in industries working with high-efficiency cooling/heating system computations.


Sign in / Sign up

Export Citation Format

Share Document