LES Simulations of Rotating Two-Pass Ribbed Duct With Coriolis and Centrifugal Buoyancy Forces at Re=100,000

Author(s):  
Cody Dowd ◽  
Danesh Tafti

The focus of this research is to predict the flow and heat transfer in a rotating two-pass duct geometry with staggered ribs using Large-Eddy Simulations (LES). The geometry consists of a U-Bend with 17 ribs in each pass. The ribs are staggered with an e/Dh = 0.1 and P/e = 10. LES is performed at a Reynolds number of 100,000, a rotation number of 0.2 and buoyancy parameters (Bo) of 0.5 and 1.0. The effects of Coriolis forces and centrifugal buoyancy are isolated and studied individually. In all cases it is found that increasing Bo from 0.5 to 1.0 at Ro = 0.2 has little impact on heat transfer. It is found that in the first pass, the heat transfer is quite receptive to Coriolis forces which augment and attenuate heat transfer at the trailing and leading walls, respectively. Centrifugal buoyancy, on the other hand has a bigger effect in augmenting heat transfer at the trailing wall than in attenuating heat transfer at the leading wall. On contrary, it aids heat transfer in the second half of the first pass at the leading wall by energizing the flow near the wall. The heat transfer in the second pass is dominated by the highly turbulent flow exiting the bend. Coriolis forces have no impact on the augmentation of heat transfer on the leading wall till the second half of the passage whereas it attenuates heat transfer at the trailing wall as soon as the flow exits the bend. Contrary to phenomenological arguments, inclusion of centrifugal buoyancy augments heat transfer over Coriolis forces alone on both the leading and trailing walls of the second pass.

2004 ◽  
Vol 127 (2) ◽  
pp. 306-320 ◽  
Author(s):  
A. K. Saha ◽  
Sumanta Acharya

Large eddy simulations (LES) and unsteady Reynolds averaged Navier-Stokes (URANS) simulations have been performed for flow and heat transfer in a rotating ribbed duct. The ribs are oriented normal to the flow and arranged in a staggered configuration on the leading and trailing surfaces. The LES results are based on a higher-order accurate finite difference scheme with a dynamic Smagorinsky model for the subgrid stresses. The URANS procedure utilizes a two equation k-ε model for the turbulent stresses. Both Coriolis and centrifugal buoyancy effects are included in the simulations. The URANS computations have been carried out for a wide range of Reynolds number (Re=12,500-100,000), rotation number (Ro=0-0.5) and density ratio (Δρ∕ρ=0-0.5), while LES results are reported for a single Reynolds number of 12,500 without and with rotation (Ro=0.12,Δρ∕ρ=0.13). Comparison is made between the LES and URANS results, and the effects of various parameters on the flow field and surface heat transfer are explored. The LES results clearly reflect the importance of coherent structures in the flow, and the unsteady dynamics associated with these structures. The heat transfer results from both LES and URANS are found to be in reasonable agreement with measurements. LES is found to give higher heat transfer predictions (5–10% higher) than URANS. The Nusselt number ratio (Nu∕Nu0) is found to decrease with increasing Reynolds number on all walls, while they increase with the density ratio along the leading and trailing walls. The Nusselt number ratio on the trailing and sidewalls also increases with rotation. However, the leading wall Nusselt number ratio shows an initial decrease with rotation (till Ro=0.12) due to the stabilizing effect of rotation on the leading wall. However, beyond Ro=0.12, the Nusselt number ratio increases with rotation due to the importance of centrifugal-buoyancy at high rotation.


Author(s):  
Mohammad A. Elyyan ◽  
Danesh K. Tafti

LES calculations are conducted for flow in a channel with dimples and protrusions on opposite walls with both surfaces heated at three Reynolds numbers, ReH = 220, 940, and 9300 ranging from laminar, weakly turbulent to fully turbulent, respectively. Turbulence generated by the separated shear layer in the dimple and along the downstream rim of the dimple is primarily responsible for heat transfer augmentation on the dimple surface. On the other hand, augmentation on the protrusion surface is mostly driven by flow impingement and flow acceleration between protrusions, while the turbulence generated in the wake has a secondary effect. Heat transfer augmentation ratios of 0.99 at ReH = 220, 2.9 at ReH = 940, and 2.5 at ReH = 9300 are obtained. Both skin friction and form losses contribute to pressure drop in the channel, with form losses increasing from 45% to 80% with an increase in the Reynolds number. Friction coefficient augmentation ratios of 1.67, 4.82 and 6.37 are obtained at ReH = 220, 940, and 9300, respectively. Based on the geometry studied, it is found that dimples and protrusions may not be viable heat transfer augmentation surfaces when the flow is steady and laminar.


2004 ◽  
Vol 126 (4) ◽  
pp. 627-636 ◽  
Author(s):  
Samer Abdel-Wahab ◽  
Danesh K. Tafti

Results from large eddy simulations (LES) of fully developed flow in a 90 deg ribbed duct are presented with rib pitch-to-height ratio P/e=10 and a rib height-to-hydraulic-diameter ratio e/Dh=0.1. Three rotation numbers Ro=0.18, 0.36, and 0.68 are studied at a nominal Reynolds number based on bulk velocity of 20 000. Centrifugal buoyancy effects are included at two Richardson numbers of Ri=12, 28 (Buoyancy parameter, Bo=0.12 and 0.30) for each rotation case. Heat transfer augmentation on the trailing side of the duct due to the action of Coriolis forces alone asymptotes to a value of 3.7±5% by Ro=0.2. On the other hand, augmentation ratios on the leading surface keep decreasing with an increase in rotation number with values ranging from 1.7 at Ro=0.18 to 1.2 at Ro=0.67. Secondary flow cells augment the heat transfer coefficient on the smooth walls by 20% to 30% over a stationary duct. Centrifugal buoyancy further strengthens the secondary flow cells in the duct cross-section which leads to an additional increase of 10% to 15%. Buoyancy also accentuates the augmentation of turbulence near the trailing wall of the duct and increases the heat transfer augmentation ratio 10% to 20% over the action of Coriolis forces alone. However, it does not have any significant effect at the leading side of the duct. The overall effect of buoyancy on heat transfer augmentation for the ribbed duct is found to be less than 10% over the effect of Coriolis forces alone. Friction on the other hand is augmented 15% to 20% at the highest buoyancy number studied. Comparison with available experiments in the literature show excellent agreement.


2013 ◽  
Vol 135 (3) ◽  
Author(s):  
Sunil Patil ◽  
Danesh Tafti

Large eddy simulations of flow and heat transfer in a square ribbed duct with rib height to hydraulic diameter of 0.1 and 0.05 and rib pitch to rib height ratio of 10 and 20 are carried out with the near wall region being modeled with a zonal two layer model. A novel formulation is used for solving the turbulent boundary layer equation for the effective tangential velocity in a generalized co-ordinate system in the near wall zonal treatment. A methodology to model the heat transfer in the zonal near wall layer in the large eddy simulations (LES) framework is presented. This general approach is explained for both Dirichlet and Neumann wall boundary conditions. Reynolds numbers of 20,000 and 60,000 are investigated. Predictions with wall modeled LES are compared with the hydrodynamic and heat transfer experimental data of (Rau et al. 1998, “The Effect of Periodic Ribs on the Local Aerodynamic and Heat Transfer Performance of a Straight Cooling Channel,”ASME J. Turbomach., 120, pp. 368–375). and (Han et al. 1986, “Measurement of Heat Transfer and Pressure Drop in Rectangular Channels With Turbulence Promoters,” NASA Report No. 4015), and wall resolved LES data of Tafti (Tafti, 2004, “Evaluating the Role of Subgrid Stress Modeling in a Ribbed Duct for the Internal Cooling of Turbine Blades,” Int. J. Heat Fluid Flow 26, pp. 92–104). Friction factor, heat transfer coefficient, mean flow as well as turbulent statistics match available data closely with very good accuracy. Wall modeled LES at high Reynolds numbers as presented in this paper reduces the overall computational complexity by factors of 60–140 compared to resolved LES, without any significant loss in accuracy.


2014 ◽  
Vol 65 (2) ◽  
pp. 103-128 ◽  
Author(s):  
L. Paniagua ◽  
O. Lehmkuhl ◽  
C. Oliet ◽  
C. D. Pérez-Segarra

Author(s):  
Mohamed Abdelhady ◽  
David H. Wood

Abstract Stranded overhead conductor cables are used to transfer electric power, often over large distances. Conductor geometry, as well as environmental conditions, affect the power carrying capacity. This paper studies the flow dynamics and heat transfer for one stranded conductor geometry in air at Reynolds number of 1,000, determined using dynamic Smagorinsky Large Eddy Simulations. Proper Orthogonal Decomposition was used to identify coherent structures. In comparison to a smooth circular cylinder, the conductor strands noticeably affect the flow dynamics and heat transfer, locally and globally.


Author(s):  
Sung-Eun Kim ◽  
Hajime Nakamura

Large eddy simulation has been carried out of turbulent flow and heat transfer around a circular cylinder in crossflow at three subcritical Reynolds numbers (Re = 3,900, 10,000, 18,900) where the flow and heat transfer characteristics change rapidly with the Reynolds number. The computations were carried out using a second-order-accurate finite-volume Navier-Stokes solver that permits use of arbitrary unstructured meshes. A fully implicit, non-iterative fractional-step method was employed to advance the solution in time. The subgrid-scale (SGS) turbulent stresses and heat fluxes were modeled using the dynamic Smagorinsky model. The LES predictions were found to be in good agreement with the experimental data of Hajime and Igarashi (2004). The salient features of turbulent heat transfer in subcritical regime such as the laminar thermal boundary layer and the rapid increase with Reynolds number both in the mean and the r.m.s. Nusselt number in the separated region are closely reproduced by the predictions. The numerical results confirmed that the heat transfer characteristics are closely correlated with the structural change in the underlying flow with the Reynolds number.


Author(s):  
Samer Abdel-Wahab ◽  
Danesh K. Tafti

Results from large eddy simulations (LES) of fully developed flow in a 90° ribbed duct are presented with rib pitch-to-height ratio P/e = 10 and a rib height-to-hydraulic-diameter ratio e/Dh = 0.1. Three rotation numbers Ro = 0.18, 0.36 and 0.68 are studied at a nominal Reynolds number based on bulk velocity of 20,000. Centrifugal buoyancy effects are included at two Richardson numbers of Ri = 12, 28 (Buoyancy number, Bo = 0.12 and 0.30) for each rotation case. Buoyancy strengthens the secondary flow cells in the duct cross-section which leads to an increase of 20% to 30% in heat transfer augmentation at the smooth walls over and above the effect of Coriolis forces. Buoyancy also accentuates the augmentation of turbulence near the trailing wall of the duct and increases the heat transfer augmentation ratio 10% to 20% over the action of Coriolis forces alone. However, it does not have any significant effect at the leading side of the duct. The overall effect of buoyancy on heat transfer augmentation for the ribbed duct is found to be less than 10% over the effect of Coriolis forces alone. Friction on the other hand is augmented 15% to 20% at the highest buoyancy number studied. Comparison with available experiments in the literature show excellent agreement.


Sign in / Sign up

Export Citation Format

Share Document