scholarly journals Large-eddy simulation of three-dimensional dunes in a steady, unidirectional flow. Part 1. Turbulence statistics

2013 ◽  
Vol 721 ◽  
pp. 454-483 ◽  
Author(s):  
Mohammad Omidyeganeh ◽  
Ugo Piomelli

AbstractWe performed large-eddy simulations of flow over a series of three-dimensional dunes at laboratory scale (Reynolds number based on the average channel depth and streamwise velocity was 18 900) using the Lagrangian dynamic eddy-viscosity subgrid-scale model. The bedform three-dimensionality was imposed by shifting a standard two-dimensional dune shape in the streamwise direction according to a sine wave. The statistics of the flow are discussed in 10 cases with in-phase and staggered crestlines, different deformation amplitudes and wavelengths. The results are validated qualitatively against experiments. The three-dimensional separation of flow at the crestline alters the distribution of wall pressure, which in turn may cause secondary flow across the stream, which directs low-momentum fluid, near the bed, toward the lobe (the most downstream point on the crestline) and high-momentum fluid, near the top surface, toward the saddle (the most upstream point on the crestline). The mean flow is characterized by a pair of counter-rotating streamwise vortices, with core radius of the order of the flow depth. However, for wavelengths smaller than the flow depth, the secondary flow exists only near the bed and the mean flow away from the bed resembles the two-dimensional case. Staggering the crestlines alters the secondary motion; the fastest flow occurs between the lobe and the saddle planes, and two pairs of streamwise vortices appear (a strong one, centred about the lobe, and a weaker one, coming from the previous dune, centred around the saddle). The distribution of the wall stress and the focal points of separation and attachment on the bed are discussed. The sensitivity of the average reattachment length, depends on the induced secondary flow, the streamwise and spanwise components of the channel resistance (the skin friction and the form drag), and the contribution of the form drag to the total resistance are also studied. Three-dimensionality of the bed increases the drag in the channel; the form drag contributes more than in the two-dimensional case to the resistance, except for the staggered-crest case. Turbulent-kinetic energy is increased in the separated shear layer by the introduction of three-dimensionality, but its value normalized by the plane-averaged wall stress is lower than in the corresponding two-dimensional dunes. The upward flow on the stoss side and higher deceleration of flow on the lee side over the lobe plane lift and broaden the separated shear layer, respectively, affecting the turbulent kinetic energy.

2000 ◽  
Vol 122 (4) ◽  
pp. 593-603 ◽  
Author(s):  
Allan G. van de Wall ◽  
Jaikrishnan R. Kadambi ◽  
John J. Adamczyk

The unsteady process resulting from the interaction of upstream vortical structures with a downstream blade row in turbomachines can have a significant impact on the machine efficiency. The upstream vortical structures or disturbances are transported by the mean flow of the downstream blade row, redistributing the time-average unsteady kinetic energy (K) associated with the incoming disturbance. A transport model was developed to take this process into account in the computation of time-averaged multistage turbomachinery flows. The model was applied to compressor and turbine geometry. For compressors, the K associated with upstream two-dimensional wakes and three-dimensional tip clearance flows is reduced as a result of their interaction with a downstream blade row. This reduction results from inviscid effects as well as viscous effects and reduces the loss associated with the upstream disturbance. Any disturbance passing through a compressor blade row results in a smaller loss than if the disturbance was mixed-out prior to entering the blade row. For turbines, the K associated with upstream two-dimensional wakes and three-dimensional tip clearance flows are significantly amplified by inviscid effects as a result of the interaction with a downstream turbine blade row. Viscous effects act to reduce the amplification of the K by inviscid effects but result in a substantial loss. Two-dimensional wakes and three-dimensional tip clearance flows passing through a turbine blade row result in a larger loss than if these disturbances were mixed-out prior to entering the blade row. [S0889-504X(00)01804-3]


1994 ◽  
Vol 270 ◽  
pp. 277-296 ◽  
Author(s):  
Hussein Mansy ◽  
Pan-Mei Yang ◽  
David R. Williams

The fine scale three-dimensional structures usually associated with streamwise vortices in the near wake of a circular cylinder have been studied at Reynolds numbers ranging from 170 to 2200. Spatially continuous velocity measurements along lines parallel to the cylinder axis were obtained with a scanning laser anemometer. To detect the streamwise vortices in the amplitude modulated velocity field, it was necessary to develop a spatial decomposition technique to split the total flow into a primary flow component and a secondary flow component. The primary flow is comprised of the mean flow and Strouhal vortices, while the secondary flow is the result of the three-dimensional streamwise vortices that are the essence of transition to turbulence. The three-dimensional flow amplitude increases in the primary vortex formation region, then saturates shortly after the maximum amplitude in the primary flow is reached. In the near-wake region the wavelength decreases approximately like Re−0.5, but increases with downstream distance. A discontinuous increase in wavelength occurs below Re = 300 suggesting a fundamental change in the character of the three-dimensional flow. At downstream distances (x/D = 10-20), the spanwise wavelength decreases from 1.42D to 1.03D as the Reynolds number increases from 300 to 1200.


2001 ◽  
Vol 426 ◽  
pp. 297-326 ◽  
Author(s):  
MAGNE LYGREN ◽  
HELGE I. ANDERSSON

Turbulent flow between a rotating and a stationary disk is studied. Besides its fundamental importance as a three-dimensional prototype flow, such flow fields are frequently encountered in rotor–stator configurations in turbomachinery applications. A direct numerical simulation is therefore performed by integrating the time-dependent Navier–Stokes equations until a statistically steady state is reached and with the aim of providing both long-time statistics and an exposition of coherent structures obtained by conditional sampling. The simulated flow has local Reynolds number r2ω/v = 4 × 105 and local gap ratio s/r = 0.02, where ω is the angular velocity of the rotating disk, r the radial distance from the axis of rotation, v the kinematic viscosity of the fluid, and s the gap width.The three components of the mean velocity vector and the six independent Reynolds stresses are compared with experimental measurements in a rotor–stator flow configuration. In the numerically generated flow field, the structural parameter a1 (i.e. the ratio of the magnitude of the shear stress vector to twice the mean turbulent kinetic energy) is lower near the two disks than in two-dimensional boundary layers. This characteristic feature is typical for three-dimensional boundary layers, and so are the misalignment between the shear stress vector and the mean velocity gradient vector, although the degree of misalignment turns out to be smaller in the present flow than in unsteady three-dimensional boundary layer flow. It is also observed that the wall friction at the rotating disk is substantially higher than at the stationary disk.Coherent structures near the disks are identified by means of the λ2 vortex criterion in order to provide sufficient information to resolve a controversy regarding the roles played by sweeps and ejections in shear stress production. An ensemble average of the detected structures reveals that the coherent structures in the rotor–stator flow are similar to the ones found in two-dimensional flows. It is shown, however, that the three-dimensionality of the mean flow reduces the inter-vortical alignment and the tendency of structures of opposite sense of rotation to overlap. The coherent structures near the disks generate weaker sweeps (i.e. quadrant 4 events) than structures in conventional two-dimensional boundary layers. This reduction in the quadrant 4 contribution from the coherent structures is believed to explain the reduced efficiency of the mean flow in producing Reynolds shear stress.


2007 ◽  
Vol 589 ◽  
pp. 479-507 ◽  
Author(s):  
E. KIT ◽  
I. WYGNANSKI ◽  
D. FRIEDMAN ◽  
O. KRIVONOSOVA ◽  
D. ZHILENKO

The flow in a turbulent mixing layer resulting from two parallel different velocity streams, that were brought together downstream of a jagged partition was investigated experimentally. The trailing edge of the partition had a short triangular ‘chevron’ shape that could also oscillate up and down at a prescribed frequency, because it was hinged to the stationary part of the partition to form a flap (fliperon). The results obtained from this excitation were compared to the traditional results obtained by oscillating a two-dimensional fliperon. Detailed measurements of the mean flow and the coherent structures, in the periodically excited and spatially developing mixing layer, and its random constituents were carried out using hot-wire anemometry and stereo particle image velocimetry.The prescribed spanwise wavelength of the chevron trailing edge generated coherent streamwise vortices while the periodic oscillation of this fliperon locked in-phase the large spanwise Kelvin–Helmholtz (K-H) rolls, therefore enabling the study of the inter- action between the two. The two-dimensional periodic excitation increases the strength of the spanwise rolls by increasing their size and their circulation, which depends on the input amplitude and frequency. The streamwise vortices generated by the jagged trailing edge distort and bend the primary K-H rolls. The present investigation endeavours to study the distortions of each mode as a consequence of their mutual interaction. Even the mean flow provides evidence for the local bulging of the large spanwise rolls because the integral width (the momentum thickness, θ), undulates along the span. The lateral location of the centre of the ensuing mixing layer (the location where the mean velocity is the arithmetic average of the two streams,y0), also suggests that these vortices are bent. Phase-locked and ensemble-averaged measurements provide more detailed information about the bending and bulging of the large eddies that ensue downstream of the oscillating chevron fliperon. The experiments were carried out at low speeds, but at sufficiently high Reynolds number to ensure naturally turbulent flow.


1999 ◽  
Vol 387 ◽  
pp. 353-396 ◽  
Author(s):  
J. ZHOU ◽  
R. J. ADRIAN ◽  
S. BALACHANDAR ◽  
T. M. KENDALL

The evolution of a single hairpin vortex-like structure in the mean turbulent field of a low-Reynolds-number channel flow is studied by direct numerical simulation. The structure of the initial three-dimensional vortex is extracted from the two-point spatial correlation of the velocity field by linear stochastic estimation given a second-quadrant ejection event vector. Initial vortices having vorticity that is weak relative to the mean vorticity evolve gradually into omega-shaped vortices that persist for long times and decay slowly. As reported in Zhou, Adrian & Balachandar (1996), initial vortices that exceed a threshold strength relative to the mean flow generate new hairpin vortices upstream of the primary vortex. The detailed mechanisms for this upstream process are determined, and they are generally similar to the mechanisms proposed by Smith et al. (1991), with some notable differences in the details. It has also been found that new hairpins generate downstream of the primary hairpin, thereby forming, together with the upstream hairpins, a coherent packet of hairpins that propagate coherently. This is consistent with the experimental observations of Meinhart & Adrian (1995). The possibility of autogeneration above a critical threshold implies that hairpin vortices in fully turbulent fields may occur singly, but they more often occur in packets. The hairpins also generate quasi-streamwise vortices to the side of the primary hairpin legs. This mechanism bears many similarities to the mechanisms found by Brooke & Hanratty (1993) and Bernard, Thomas & Handler (1993). It provides a means by which new quasi-streamwise vortices, and, subsequently, new hairpin vortices can populate the near-wall layer.


1972 ◽  
Vol 23 (3) ◽  
pp. 211-228 ◽  
Author(s):  
H P Horton

SummaryMeasurements are presented of the mean flow properties of some three-dimensional turbulent boundary layers re-developing after reattachment behind short separation bubbles yawed at 26.5° to the main stream. For these measurements, Rθ11varied from about 550 to 1450. It was found that, where the pressure gradient parameter (ν/ρu3τ1)∂p/∂s was not greater than about 0.05, the flow in the local external streamline direction conformed well with empirical laws for fully-attached two-dimensional layers with regard to the mean velocity profiles, shape parameter relationships and skin friction laws, giving support to the usual assumption that these two-dimensional relationships may be applied to the streamwise flow in three-dimensional layers, subject to the limitation on the pressure gradient parameter. The cross-flow profiles, on the other hand, were not generally fitted well by the often-used representations of Mager and Johnston. The variations of the kinetic energy dissipation coefficient and the entrainment rate were deduced for one of the layers, both quantities being found to be higher than those predicted by empirical relationships for conventionally-developing two-dimensional layers. However, the energy dissipation is in fair agreement with that in a similarly re-developing two-dimensional flow.


Author(s):  
Young-Woo Yi ◽  
◽  
Bhupendra Singh Chauhan ◽  
Hee-Chang Lim ◽  
◽  
...  

Large Eddy Simulations (LES) has been widely applied and used in several decades to simulate a turbulent boundary layer in the numerical domain. In this study, we aimed to make a synthetic inflow generator (SIG) yielding an appropriate property of turbulent boundary layer in the inlet section and making quick development in the downstream of a three-dimensional domain. In order to achieve turbulent boundary layer quickly in a limited domain, the oscillating term was implemented in the well-defined boundary layer, which was expected to make faster convergence in the calculation. Cholesky decomposition was also applied to possess turbulent statistics such as the randomness and correlation of turbulent flow. In a result, the oscillating inflow did not show the faster convergence, but it indicated a possibility to improve statistical quantities in the downstream. In addition, regarding the mean flow characteristics were very close to the calculation without the oscillating flow. On the other hand, the turbulent statistics were improved depending on the oscillating magnitude.


2008 ◽  
Vol 603 ◽  
pp. 189-206 ◽  
Author(s):  
S. MASUDA ◽  
S. FUKUDA ◽  
M. NAGATA

We analyse the stability of plane Poiseuille flow with a streamwise system rotation. It is found that the instability due to two-dimensional perturbations, which sets in at the well-known critical Reynolds number, Rc = 5772.2, for the non-rotating case, is delayed as the rotation is increased from zero, showing a stabilizing effect of rotation. As the rotation is increased further, however, the laminar flow becomes most unstable to perturbations which are three-dimensional. The critical Reynolds number due to three-dimensional perturbations at this higher rotation case is many orders of magnitude less than the corresponding value due to two-dimensional perturbations. We also perform a nonlinear analysis on a bifurcating three-dimensional secondary flow. The secondary flow exhibits a spiral vortex structure propagating in the streamwise direction. It is confirmed that an antisymmetric mean flow in the spanwise direction is generated in the secondary flow.


1996 ◽  
Vol 118 (2) ◽  
pp. 347-352 ◽  
Author(s):  
R. G. Dominy ◽  
D. A. Kirkham

Interturbine diffusers provide continuity between HP and LP turbines while diffusing the flow upstream of the LP turbine. Increasing the mean turbine diameter offers the potential advantage of reducing the flow factor in the following stages, leading to increased efficiency. The flows associated with these interturbine diffusers differ from those in simple annular diffusers both as a consequence of their high-curvature S-shaped geometry and of the presence of wakes created by the upstream turbine. It is shown that even the simplest two-dimensional wakes result in significantly modified flows through such ducts. These introduce strong secondary flows demonstrating that fully three-dimensional, viscous analysis methods are essential for correct performance modeling.


Sign in / Sign up

Export Citation Format

Share Document