Flow and Heat Transfer in an Internally Ribbed Duct With Rotation: An Assessment of Large Eddy Simulations and Unsteady Reynolds-Averaged Navier-Stokes Simulations

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):  
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 side walls 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):  
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.


2021 ◽  
Author(s):  
Matthew Searle ◽  
Arnab Roy ◽  
James Black ◽  
Doug Straub ◽  
Sridharan Ramesh

Abstract In this paper, experimental and numerical investigations of three variants of internal cooling configurations — dimples only, ribs only and ribs with dimples have been explored at process conditions (96°C and 207bar) with sCO2 as the coolant. The designs were chosen based on a review of advanced internal cooling features typically used for air-breathing gas turbines. The experimental study described in this paper utilizes additively manufactured square channels with the cooling features over a range of Reynolds number from 80,000 to 250,000. Nusselt number is calculated in the experiments utilizing the Wilson Plot method and three heat transfer characteristics — augmentation in Nusselt number, friction factor and overall Thermal Performance Factor (TPF) are reported. To explore the effect of surface roughness introduced due to additive manufacturing, two baseline channel flow cases are considered — a conventional smooth tube and an additively manufactured square tube. A companion computational fluid dynamics (CFD) simulation is also performed for the corresponding cooling configurations reported in the experiments using the Reynolds Averaged Navier Stokes (RANS) based turbulence model. Both experimental and computational results show increasing Nusselt number augmentation as higher Reynolds numbers are approached, whereas prior work on internal cooling of air-breathing gas turbines predict a decay in the heat transfer enhancement as Reynolds number increases. Comparing cooling features, it is observed that the “ribs only” and “ribs with dimples” configurations exhibit higher Nusselt number augmentation at all Reynolds numbers compared to the “dimples only” and the “no features” configurations. However, the frictional losses are almost an order of magnitude higher in presence of ribs.


Author(s):  
A. K. Saha ◽  
Sumanta Acharya

The flow and heat transfer in ribbed coolant passages of aspect ratios (AR) of 1:1, 4:1, and 1:4 are numerically studied through the solution of the Unsteady Reynolds Averaged Navier-Stokes (URANS) equations. The ribs are oriented normal to the flow and arranged in a staggered configuration on the leading and trailing surfaces. The URANS procedure can resolve large-scale bulk unsteadiness, and utilizes a two equation k-ε model for the turbulent stresses. Both Coriolis and centrifugal buoyancy effects are included in the simulations. The computations are carried out for a fixed Reynolds number of 25000 and density ratio of 0.13 while the Rotation number has been varied between 0.12–0.50. The average duct heat transfer is the highest for the 4:1 AR case. For this case, the secondary flow structures consist of multiple roll cells that direct flow both to the trailing and leading surfaces. The 1:4 AR duct shows flow reversal along the leading surface at high rotation numbers with multiple rolls in the secondary flow structures near the leading wall. For this AR, the potential for conduction-limited heat transfer along the leading surface is identified. At high rotation number, both the 1:1 and 4:1 AR cases exhibit loss of axial periodicity over one inter-rib module. The friction factor reveals an increase with the rotation number for all aspect ratio ducts, and shows a sudden jump in its value at a critical rotation number because of either loss of spatial periodicity or the onset of backflow.


2019 ◽  
Vol 2019 ◽  
pp. 1-9
Author(s):  
Xuan Zhang ◽  
Taocheng Zhao ◽  
Suchen Wu ◽  
Feng Yao

Although roughness is negligible for laminar flow through tubes in classic fluid mechanics, the surface roughness may play an important role in microscale fluid flow due to the large ratio of surface area to volume. To further verify the influence of rough surfaces on microscale liquid flow and heat transfer, a performance test system of heat transfer and liquid flow was designed and built, and a series of experimental examinations are conducted, in which the microchannel material is stainless steel and the working medium is methanol. The results indicate that the surface roughness plays a significant role in the process of laminar flow and heat transfer in microchannels. In microchannels with roughness characteristics, the Poiseuille number of liquid laminar flow relies not only on the cross section shape of the rough microchannels but also on the Reynolds number of liquid flow. The Poiseuille number of liquid laminar flow in rough microchannels increases with increasing Reynolds number. In addition, the Nusselt number of liquid laminar heat transfer is related not only to the cross section shape of a rough microchannel but also to the Reynolds number of liquid flow, and the Nusselt number increases with increasing Reynolds number.


2006 ◽  
Vol 129 (6) ◽  
pp. 685-696 ◽  
Author(s):  
Guoguang Su ◽  
Hamn-Ching Chen ◽  
Je-Chin Han

Computations with multi-block chimera grids were performed to study the three-dimensional turbulent flow and heat transfer in a rotating rectangular channel with staggered arrays of pin-fins. The channel aspect ratio (AR) is 4:1, the pin length to diameter ratio (H∕D) is 2.0, and the pin spacing to diameter ratio is 2.0 in both the stream-wise (S1∕D) and span-wise (S2∕D) directions. A total of six calculations have been performed with various combinations of rotation number, Reynolds number, and coolant-to-wall density ratio. The rotation number and inlet coolant-to-wall density ratio varied from 0.0 to 0.28 and from 0.122 to 0.20, respectively, while the Reynolds number varied from 10,000 to 100,000. For the rotating cases, the rectangular channel was oriented at 150deg with respect to the plane of rotation to be consistent with the configuration of the gas turbine blade. A Reynolds-averaged Navier-Stokes (RANS) method was employed in conjunction with a near-wall second-moment turbulence closure for detailed predictions of mean velocity, mean temperature, and heat transfer coefficient distributions.


2004 ◽  
Vol 126 (6) ◽  
pp. 997-1005 ◽  
Author(s):  
R. Saleh ◽  
A. B. Rahimi

The unsteady viscous flow and heat transfer in the vicinity of an axisymmetric stagnation point of an infinite moving cylinder with time-dependent axial velocity and with uniform normal transpiration Uo are investigated. The impinging free stream is steady and with a constant strain rate k¯. An exact solution of the Navier–Stokes equations and energy equation is derived in this problem. A reduction of these equations is obtained by use of appropriate transformations for the most general case when the transpiration rate is also time-dependent but results are presented only for uniform values of this quantity. The general self-similar solution is obtained when the axial velocity of the cylinder and its wall temperature or its wall heat flux vary as specified time-dependent functions. In particular, the cylinder may move with constant speed, with exponentially increasing–decreasing axial velocity, with harmonically varying axial speed, or with accelerating–decelerating oscillatory axial speed. For self-similar flow, the surface temperature or its surface heat flux must have the same types of behavior as the cylinder motion. For completeness, sample semisimilar solutions of the unsteady Navier–Stokes and energy equations have been obtained numerically using a finite-difference scheme. Some of these solutions are presented for special cases when the time-dependent axial velocity of the cylinder is a step-function, and a ramp function. All the solutions above are presented for Reynolds numbers, Re=ka¯2/2υ, ranging from 0.1 to 100 for different values of dimensionless transpiration rate, S=Uo/ka¯, where a is cylinder radius and υ is kinematic viscosity of the fluid. Absolute value of the shear-stresses corresponding to all the cases increase with the increase of Reynolds number and suction rate. The maximum value of the shear- stress increases with increasing oscillation frequency and amplitude. An interesting result is obtained in which a cylinder moving with certain exponential axial velocity function at any particular value of Reynolds number and suction rate is axially stress-free. The heat transfer coefficient increases with the increasing suction rate, Reynolds number, Prandtl number, oscillation frequency and amplitude. Interesting means of cooling and heating processes of cylinder surface are obtained using different rates of transpiration. It is shown that a cylinder with certain type of exponential wall temperature exposed to a temperature difference has no heat transfer.


2005 ◽  
Vol 127 (5) ◽  
pp. 486-498 ◽  
Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy forces are included in this study. A direct approach is presented for the unsteady calculation of the nondimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12,500, a rotation number (Ro) of 0.12, and an inlet coolant-to-wall density ratio Δρ/ρ of 0.13. The predicted time and space-averaged Nusselt numbers are shown to compare satisfactorily with the published experimental data. Time sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer behavior. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of the flowfield indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.


Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy force are included in this study. A direct approach is presented for the unsteady calculation of the non-dimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12, 500, a Rotation number (Ro) of 0.12 and an inlet coolant-to-wall density ratio (Δρ/ρ) of 0.13. The time-averaged Nusselt numbers compare satisfactorily with the data of Wagner et al. (J. Turbomachinery, Vol. 114, pp. 847–857). Time-sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer processes. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of 200 snapshots indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.


2006 ◽  
Vol 129 (1) ◽  
pp. 106-115 ◽  
Author(s):  
A. B. Rahimi ◽  
R. Saleh

The unsteady viscous flow and heat transfer in the vicinity of an axisymmetric stagnation point of an infinite rotating circular cylinder with transpiration U0 are investigated when the angular velocity and wall temperature or wall heat flux all vary arbitrarily with time. The free stream is steady and with a strain rate of Γ. An exact solution of the Navier-Stokes equations and energy equation is derived in this problem. A reduction of these equations is obtained by the use of appropriate transformations for the most general case when the transpiration rate is also time-dependent but results are presented only for uniform values of this quantity. The general self-similar solution is obtained when the angular velocity of the cylinder and its wall temperature or its wall heat flux vary as specified time-dependent functions. In particular, the cylinder may rotate with constant speed, with exponentially increasing/decreasing angular velocity, with harmonically varying rotation speed, or with accelerating/decelerating oscillatory angular speed. For self-similar flow, the surface temperature or its surface heat flux must have the same types of behavior as the cylinder motion. For completeness, sample semi-similar solutions of the unsteady Navier-Stokes equations have been obtained numerically using a finite-difference scheme. Some of these solutions are presented for special cases when the time-dependent rotation velocity of the cylinder is, for example, a step-function. All the solutions above are presented for Reynolds numbers, Re=Γa2∕2υ, ranging from 0.1 to 1000 for different values of Prandtl number and for selected values of dimensionless transpiration rate, S=U0∕Γa, where a is cylinder radius and υ is kinematic viscosity of the fluid. Dimensionless shear stresses corresponding to all the cases increase with the increase of Reynolds number and suction rate. The maximum value of the shear stress increases with increasing oscillation frequency and amplitude. An interesting result is obtained in which a cylinder rotating with certain exponential angular velocity function and at particular value of Reynolds number is azimuthally stress-free. Heat transfer is independent of cylinder rotation and its coefficient increases with the increasing suction rate, Reynolds number, and Prandtl number. Interesting means of cooling and heating processes of cylinder surface are obtained using different rates of transpiration.


Sign in / Sign up

Export Citation Format

Share Document