Large Eddy Simulations of Flow and Heat Transfer in Rotating Ribbed Duct Flows

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.


Author(s):  
Angela Wu ◽  
Seunghwan Keum ◽  
Volker Sick

In this study, the effects of the thermal boundary conditions at the engine walls on the predictions of Large-Eddy Simulations (LES) of a motored Internal Combustion Engine (ICE) were examined. Two thermal boundary condition cases were simulated. One case used a fixed, uniform wall temperature, which is typically used in conventional LES modeling of ICEs. The second case utilized a Conjugate Heat Transfer (CHT) modeling approach to obtain temporally and spatially varying wall temperature. The CHT approach solves the coupled heat transfer problem between fluid and solid domains. The CHT case included the solid valves, piston, cylinder head, cylinder liner, valve seats, and spark plug geometries. The simulations were validated with measured bulk flow, near-wall flow, surface temperature, and surface heat flux. The LES quality of both simulations was also discussed. The CHT results show substantial spatial, temporal, and cyclic variability of the wall heat transfer. The surface temperature dynamics obtained from the CHT model compared well with measurements during the compression stroke, but the absolute magnitude was 5 K (or 1.4%) off and the prediction of the drop in temperature after top dead center suffered from temporal resolution limitations. Differences in the predicted flow and temperature fields between the uniform surface temperature and CHT simulations show the impact of the surface temperature on bulk behavior.


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):  
Yasuo Hattori ◽  
Hitoshi Suto ◽  
Keisuke Nakao ◽  
Yuma Hasebe ◽  
Shuji Ishihara

Abstract We performed large-eddy simulations by using an open source CFD code, OpenFOAM to examine the performance of large-eddy simulations for examining turbulence heat-transfer processes of forced-convection in a pipe with non-homogeneous thermal boundary conditions; an accurate description of such processes is of practical interest of nuclear engineering. Special attention was paid to the domain size in the streamwise direction, which must be closely related to the turbulence processes with super structures. Three domain sizes were used: the size was varied from 5R to 100R, where R is the radius of the pipe. The turbulence intensities of temperature fluctuations near the heated surface strongly depended on the domain size. This was because that the turbulence intensities were closely related to large-scale fluid motions, the scale of which is much larger than 25R, and such large-scale fluid motions interfered in the dynamics of turbulence heat transfer processes near the heated surface. This indicates that the large-eddy simulations for estimating the turbulence heat-transfer rate in the pipe must require the large domain size in the streamwise direction with well resolved grids to capture turbulence eddies near the heated surface.


Author(s):  
J S Lee ◽  
R H Pletcher

Turbulent combined flow of forced and natural convection was investigated using large eddy simulations for horizontal and vertical channels with significant heat transfer. The walls were maintained at constant temperatures, one heated and the other cooled, at temperature ratios of 1.01, 1.99, and 3.00, respectively. Results showed that with increasing the Grashof number, large-scale turbulent motions emerged near the wall, resulting in significant changes in turbulent intensities for the horizontal channel flow case. Aiding and opposing flows, however, for the vertical channel, emerge near the heated and cooled walls, respectively, while the pressure gradient drives the mean flow upwards. Buoyancy effects on the mean velocity, temperature, and turbulent intensities were observed near the walls. In the aiding flow, the turbulent intensities and heat transfer were suppressed and the flow was relaminarized at large values of the Grash of number. In the opposing flow, however, turbulence was enhanced with increasing velocity fluctuations.


Sign in / Sign up

Export Citation Format

Share Document