Three-Dimensional Mixed Convection in Plane Symmetric-Sudden Expansion: Bifurcated Flow Regime

2006 ◽  
Vol 129 (7) ◽  
pp. 819-826 ◽  
Author(s):  
M. Thiruvengadam ◽  
B. F. Armaly ◽  
J. A. Drallmeier

Simulations of three-dimensional laminar mixed convection in a vertical duct with plane symmetric sudden expansion are presented to illustrate the effects of the buoyancy-assisting force and the duct’s aspect ratio on flow bifurcation and heat transfer. The stable laminar bifurcated flow regime that develops in this geometry at low buoyancy levels leads to nonsymmetric temperature and heat transfer distributions in the transverse direction, but symmetric distributions with respect to the center width of the duct in the spanwise direction. As the buoyancy force increases, due to increases in wall heat flux, flow bifurcation diminishes and both the flow and the thermal fields become symmetric at a critical wall heat flux. The size of the primary recirculation flow region adjacent to the sudden expansion increases on one of the stepped walls and decreases on the other stepped wall as the wall heat flux increases. The maximum Nusselt number that develops on one of the stepped walls in the bifurcated flow regime is significantly larger than the one that develops on the other stepped wall. The critical wall heat flux increases as the duct’s aspect ratio increases for fixed Reynolds number. The maximum Nusselt number that develops in the bifurcated flow regime increases as the duct’s aspect ratio increases for fixed wall heat flux and Reynolds number.

2017 ◽  
Vol 835 ◽  
pp. 1157-1198 ◽  
Author(s):  
Shingo Motoki ◽  
Genta Kawahara ◽  
Masaki Shimizu

Optimal heat transfer enhancement has been explored theoretically in plane Couette flow. The vector field (referred to as the ‘velocity’) to be optimised is time independent and divergence free, and temperature is determined in terms of the velocity as a solution to an advection-diffusion equation. The Prandtl number is set to unity, and consistent boundary conditions are imposed on the velocity and the temperature fields. The excess of a wall heat flux (or equivalently total scalar dissipation) over total energy dissipation is taken as an objective functional, and by using a variational method the Euler–Lagrange equations are derived, which are solved numerically to obtain the optimal states in the sense of maximisation of the functional. The laminar conductive field is an optimal state at low Reynolds number $Re\sim 10^{0}$. At higher Reynolds number $Re\sim 10^{1}$, however, the optimal state exhibits a streamwise-independent two-dimensional velocity field. The two-dimensional field consists of large-scale circulation rolls that play a role in heat transfer enhancement with respect to the conductive state as in thermal convection. A further increase of the Reynolds number leads to a three-dimensional optimal state at $Re\gtrsim 10^{2}$. In the three-dimensional velocity field there appear smaller-scale hierarchical quasi-streamwise vortex tubes near the walls in addition to the large-scale rolls. The streamwise vortices are tilted in the spanwise direction so that they may produce the anticyclonic vorticity antiparallel to the mean-shear vorticity, bringing about significant three-dimensionality. The isotherms wrapped around the tilted anticyclonic vortices undergo the cross-axial shear of the mean flow, so that the spacing of the wrapped isotherms is narrower and so the temperature gradient is steeper than those around a purely streamwise (two-dimensional) vortex tube, intensifying scalar dissipation and so a wall heat flux. Moreover, the tilted anticyclonic vortices induce the flow towards the wall to push low- (or high-) temperature fluids on the hot (or cold) wall, enhancing a wall heat flux. The optimised three-dimensional velocity fields achieve a much higher wall heat flux and much lower energy dissipation than those of plane Couette turbulence.


Author(s):  
Wai Hing Wong ◽  
Normah Mohd. Ghazali

Kertas kerja ini membincangkan simulasi berangka ke atas sinki haba saluran mikro dalam penyejukan alatan mikroelektronik. Model Dinamik Bendalir Berkomputer (CFD) tiga dimensi dibina menggunakan pakej komersil, FLUENT, untuk mengkaji fenomenon aliran bendalir dan pemindahan haba konjugat di dalam suatu sinki haba segi empat yang diperbuat daripada silikon. Model ditentusahkan dengan keputusan daripada uji kaji dan pengkajian berangka yang lepas untuk lingkungan nombor Reynolds kurang daripada 400 berdasarkan diameter hidraulik 86 mm. Kajian ini mengambil kira kesan kelikatan bendalir yang bersandaran dengan suhu dan keadaan aliran pra–membangun dari segi hidrodinamik dan haba. Model memberi maklumat tentang taburan suhu dan fluks haba yang terperinci di dalam sinki haba saluran mikro. Kecerunan suhu yang tinggi dicatat pada kawasan pepejal berdekatan dengan sumber. Fluks haba paling tinggi didapati pada dinding tepi saluran mikro diikuti oleh dinding atas dan bawah. Purata pekali pemindahan haba yang lebih tinggi bagi silikon menjadikan ia bahan binaan sinki haba saluran mikro yang lebih baik berbanding dengan kuprum dan aluminium. Peningkatan nisbah aspek saluran mikro yang bersegi empat memberi kecekapan penyejukan yang lebih tinggi kerana kelebaran saluran yang berkurangan memberi kecerunan halaju yang lebih tinggi dalam saluran. Nisbah aspek yang optimum yang diperoleh adalah dalam lingkungan 3.7 – 4.1. Kata kunci: Saluran mikro, CFD, FLUENT, simulasi berangka, penyejukan mikroelektron The paper discusses the numerical simulation of a micro–channel heat sink in microelectronics cooling. A three–dimensional Computational Fluid Dynamics (CFD) model was built using the commercial package, FLUENT, to investigate the conjugate fluid flow and heat transfer phenomena in a silicon–based rectangular microchannel heatsink. The model was validated with past experimental and numerical work for Reynolds numbers less than 400 based on a hydraulic diameter of 86 mm. The investigation was conducted with consideration of temperaturedependent viscosity and developing flow, both hydrodynamically and thermally. The model provided detailed temperature and heat flux distributions in the microchannel heatsink. The results indicate a large temperature gradient in the solid region near the heat source. The highest heat flux is found at the side walls of the microchannel, followed by top wall and bottom wall due to the wall interaction effects. Silicon is proven to be a better microchannel heatsink material compared to copper and aluminum, indicated by a higher average heat transfer. A higher aspect ratio in a rectangular microchannel gives higher cooling capability due to high velocity gradient around the channel when channel width decreases. Optimum aspect ratio obtained is in the range of 3.7 – 4.1. Key words: Microchannel, CFD, FLUENT, numerical simulation, microeletronics cooling


2010 ◽  
Vol 133 (2) ◽  
Author(s):  
Bi-Li Deng ◽  
Xin-Rong Zhang ◽  
Hiroshi Yamaguchi

This study presents a computational investigation of forced convection of supercritical CO2 flow in plane symmetric sudden expansion duct at an expansion ratio of 2 in flow asymmetric regime. Computations were conducted at various Reynolds numbers in flow asymmetric regime and low wall heat fluxes below 500 W/m2 to examine the Reynolds number and thermal effects on the flow and heat transfer of the bifurcated flow. General flow features and temperature distributions are presented. The transitional Reynolds numbers above, which a third recirculation region will appear at different wall heat flux are presented, and thus thermal effects on the flow stability are discussed. Reynolds number and thermal effects on distributions of wall skin friction, pressure coefficient, and Nusselt number are presented and discussed.


2016 ◽  
Vol 20 (5) ◽  
pp. 1621-1632 ◽  
Author(s):  
Hussein Togun ◽  
Tuqa Abdulrazzaq ◽  
Salim Kazi ◽  
Ahmad Badarudin

This paper presents a study of heat transfer to turbulent air flow in the abrupt axisymmetric expansion of an annular pipe. The experimental investigations were performed in the Reynolds number range from 5000 to 30000, the heat flux varied from 1000 to 4000 W/m2, and the expansion ratio was maintained at D/d=1, 1.25, 1.67 and 2. The sudden expansion was created by changing the inner diameter of the entrance pipe to an annular passage. The outer diameter of the inner pipe and the inner diameter of the outer pipe are 2.5 and 10 cm, respectively, where both of the pipes are subjected to uniform heat flux. The distribution of the surface temperature of the test pipe and the local Nusselt number are presented in this investigation. Due to sudden expansion in the cross section of the annular pipe, a separation flow was created, which enhanced the heat transfer. The reduction of the surface temperature on the outer and inner pipes increased with the increase of the expansion ratio and the Reynolds number, and increased with the decrease of the heat flux to the annular pipe. The peak of the local Nusselt number was between 1.64 and 1.7 of the outer and inner pipes for Reynolds numbers varied from 5000 to 30000, and the increase of the local Nusselt number represented the augmentation of the heat transfer rate in the sudden expansion of the annular pipe. This research also showed a maximum heat transfer enhancement of 63-78% for the outer and inner pipes at an expansion ratio of D/d=2 at a Re=30000 and a heat flux of 4000W/m2.


Author(s):  
J. H. Nie ◽  
B. F. Armaly

Simulations of three-dimensional flow and heat transfer in laminar incompressible forced convection in plane symmetric sudden expansion (backward-facing step in rectangular duct) are presented for different Reynolds numbers. The duct’s downstream (H) and upstream (h) heights are 0.04m and 0.02m, respectively, thus providing a step height (S) of 0.01m and an expansion ratio of 2. The duct’s width (W) is 0.08m, thus resulting in an aspect ratio of 4 before and 2 after the expansion, respectively. The incoming flow is considered to be isothermal, hydro-dynamically steady and fully developed. Uniform and constant heat flux is specified on the stepped walls, while the other walls are treated as adiabatic surfaces. The flow appears to be symmetric for the low Reynolds number range that is considered in this study (Re=150). A “jet-like” flow develops near the sidewall and its impingement on the stepped wall creates a swirling flow inside the primary recirculation region adjacent to the stepped wall, and that is responsible for creating a maximum in the Nusselt number distribution. The results reveal that the location where the streamwise component of wall shear stress is zero on the stepped wall does not coincide with the location of the outer edge of the primary recirculation region, especially in the region near the sidewall. Neither one of these boundary lines represents the reattachment region of the separated flow in the region adjacent to the sidewall. The maximum Nusselt number on the stepped wall is located inside the primary recirculation region and is not identical to the “jet-like” flow impingement point. The maximum friction coefficient on the stepped wall is located inside the primary recirculation region, and it is at the center of the duct for the Reynolds number range considered in this study. The minimum friction coefficient on the stepped wall is located at the impingement of the “jet-like” flow.


Author(s):  
Shawn Siroka ◽  
Melissa Shallcross ◽  
Stephen Lynch

Cylindrical pins, often called pin fins, are used to create turbulence and promote convective heat transfer within many devices, ranging from computer heat sinks to the trailing edge of jet engine turbine blades. Previous experiments have measured the time-averaged heat transfer over a single pin as well as the flow fields around the pin. However, in this study, focus is placed on the instantaneous heat flux around the centerline of a low aspect-ratio pin within an array. Time-mean and unsteady convective heat flux are measured around the circumference of an isothermal heated test pin via a microsensor located at the surface. The pin is positioned at various locations within a staggered array in a large-scale wind tunnel. Reynolds numbers from 3,000 to 50,000, based on pin diameter and maximum velocity between pins, are tested with a streamwise spacing of 1.73 diameters between rows, a spanwise spacing of 2 diameters, and a pin height of 1 diameter. The time-averaged and standard deviation of convective heat flux around the pin is higher over most of the pin surface for pins in downstream row positions of an array relative to the first row pin, except in the wake which has similar levels for all rows. For a given pin position in the array, as the Reynolds number increases, the point of minimum heat transfer moves circumferentially upstream on the pin fin, corresponding to earlier transition of the pin boundary layer. Also, for a given Reynolds number, the minimum heat transfer point on the pin circumference moves upstream for pins further into the array, due to the high turbulence levels within the array which cause early transition. For a single pin row with no downstream pins, heat transfer fluctuations are very high on the backside of the pin due to the significant unsteadiness in the pin wake, but heat transfer fluctuations are suppressed for a pin with downstream rows due to the confining effects of the close spacing. The results from this study can be used to design pin-fin arrays that take advantage of unsteadiness and increase overall convective heat transfer for various industry components.


1974 ◽  
Vol 64 (1) ◽  
pp. 111-128 ◽  
Author(s):  
F. Durst ◽  
A. Melling ◽  
J. H. Whitelaw

Flow visualization and laser-anemometry measurements are reported in the flow downstream of a plane 3: 1 symmetric expansion in a duct with an aspect ratio of 9·2: 1 downstream of the expansion. The flow was found to be markedly dependent on Reynolds number, and strongly three-dimensional even well away from the channel corners except at the lowest measurable velocities. The measurements at a Reynolds number of 56 indicated that the separation regions behind each step were of equal length. Symmetric velocity profiles existed from the expansion to a fully developed, parabolic profile far downstream, although there were substantial three-dimensional effects in the vicinity of the separation regions. The velocity profiles were in good agreement with those obtained by solving the two-dimensional momentum equation. At a Reynolds number of 114, the two separation regions were of different lengths, leading to asymmetric velocity profiles; three dimensional effects were much more pronounced. At a Reynolds number of 252, a third separation zone was found on one wall, downstream of the smaller of the two separation zones adjacent to the steps. As at the lower Reynolds numbers, the flow was very stable. At higher Reynolds numbers the flow became less stable and periodicity became increasingly important in the main stream; this was accompanied by a highly disturbed fluid motion in the separation zones, as the flow tended towards turbulence.


2013 ◽  
Vol 597 ◽  
pp. 3-8
Author(s):  
Lahoucine Belarche ◽  
Btissam Abourida ◽  
Slawomir Smolen ◽  
Touria Mediouni

Natural convection in inclined cubic cavity, discretely heated, is studied numerically using a three-dimensional finite volume formulation. Two heating square portions are placed on the vertical wall of the enclosure, while the rest of the considered wall is adiabatic. These sections, similar to the integrated electronic components, generate a heat flux q". The opposite vertical wall is maintained at a cold uniform temperature Tc and the other walls are adiabatic. The fluid flow and heat transfer in the cavity are studied for different sets of the governing parameters, namely the Rayleigh number Ra (103 ≤ Ra ≤ 107), the cavity inclination γ (- 45° ≤ γ ≤ 45°) and the position of the heating sections λ (0.3 ≤ λ ≤ 0.7). The dimensions of the heater sections, ε = D / H and the longitudinal aspect ratio of the cavity Ax = H / L are respectively fixed to 0.35 and 1.


Sign in / Sign up

Export Citation Format

Share Document