scholarly journals Large-eddy simulation of three-dimensional dunes in a steady, unidirectional flow. Part 2. Flow structures

2013 ◽  
Vol 734 ◽  
pp. 509-534 ◽  
Author(s):  
Mohammad Omidyeganeh ◽  
Ugo Piomelli

AbstractWe performed large-eddy simulations of the flow over a series of three-dimensional (3D) dunes at laboratory scale. The bedform three-dimensionality was imposed by shifting a standard two-dimensional (2D) dune shape in the streamwise direction according to a sine wave. The turbulence statistics were discussed in Part 1 of this article (Omidyeganeh & Piomelli, J. Fluid Mech., vol. 721, 2013, pp. 454–483). Coherent flow structures and their statistics are discussed concentrating on two cases with the same crestline amplitudes and wavelengths but different crestline alignments: in-phase and staggered. The present paper shows that the induced large-scale mean streamwise vortices are the primary factor that alters the features of the instantaneous flow structures. Wall turbulence is insensitive to the crestline alignment; alternating high- and low-speed streaks appear in the internal boundary layer developing on the stoss side, whereas over the node plane (the plane normal to the spanwise direction at the node of the crestline), they are inclined towards the lobe plane (the plane normal to the spanwise direction at the most downstream point of the crestline) due to the mean spanwise pressure gradient. Spanwise vortices (rollers) generated by Kelvin–Helmholtz instability in the separated shear layer appear regularly over the lobe with much larger length scale than those over the saddle (the plane normal to the spanwise direction at the most upstream point of the crestline). Rollers over the lobe may extend to the saddle plane and affect the reattachment features; their shedding is more frequent than in 2D geometries. Vortices shed from the separated shear layer in the lobe plane undergo a three-dimensional instability while being advected downstream, and rise toward the free surface. They develop into a horseshoe shape (similar to the 2D case) and affect the whole channel depth, whereas those generated near the saddle are advected downstream and toward the bed. When the tip of such a horseshoe reaches the free surface, the ejection of flow at the surface causes ‘boils’ (upwelling events on the surface). Strong boil events are observed on the surface of the lobe planes of 3D dunes more frequently than in the saddle planes. They also appear more frequently than in the corresponding 2D geometry. The crestline alignment of the dune alters the dynamics of the flow structures, in that they appear in the lobe plane and are advected towards the saddle plane of the next dune, where they are dissipated. Boil events occur at a higher frequency in the staggered alignment, but with less intensity than in the in-phase alignment.

Author(s):  
Chuang Jin ◽  
Giovanni Coco ◽  
Rafael O. Tinoco ◽  
Pallav Ranjan ◽  
Jorge San Juan ◽  
...  

Author(s):  
J. Kouwa ◽  
Y. Iso ◽  
F. Polidoro ◽  
S. Gautier

Convective heat transfer in the cavity between two corotating disks is of great importance for turbomachinery applications. The complex three dimensional and unsteady flow structures induced by the Coriolis forces inside the cavity, and therefore the resulting heat transfer, are challenging to be measured in an experiment or predicted by simulation. In this paper a simplified cavity geometry, characterized experimentally by Long at al., has been chosen. The results obtained with a Very Large Eddy Simulation using Lattice-Boltzmann Method for two operating point with different rotation speeds are compared to the experimental heat transfer coefficients at the wall. The simulation results show the characteristic flow structures and behavior induced by the different regimes. A sensitivity analysis of the results is presented, both for numerical parameters such as grid resolution and for physical parameters, namely the throughflow velocity profile and shroud temperature.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 243
Author(s):  
Abolfazl Asnaghi ◽  
Rickard E. Bensow

The simulation of fully turbulent, three-dimensional, cavitating flow over Delft twisted foil is conducted by an implicit large eddy simulation (LES) approach in both smooth and tripped conditions, the latter by including leading-edge roughness. The analysis investigates the importance of representing the roughness elements on the flow structures in the cavitation prediction. The results include detailed comparisons of cavitation pattern, vorticity distribution, and force predictions with the experimental measurements. It is noted that the presence of roughness generates very small cavitating vortical structures which interact with the main sheet cavity developing over the foil to later form a cloud cavity. Very similar to the experimental observation, these interactions create a streaky sheet cavity interface which cannot be captured in the smooth condition, influencing both the richness of structures in the detached cloudy cavitation as well as the extent and transport of vapour. It is further found to have a direct impact on the pressure distribution, especially in the mid-chord region where the shed cloud cavity collapses.


2015 ◽  
Vol 09 (05) ◽  
pp. 1540002 ◽  
Author(s):  
Chang Lin ◽  
Ming-Jer Kao ◽  
Guang-Wei Tzeng ◽  
Wei-Ying Wong ◽  
James Yang ◽  
...  

The characteristics of flow fields for a complete evolution of the non-breaking solitary wave, having a wave-height to water-depth ratio of 0.363 and propagating over a 1:5 sloping bottom, are investigated experimentally. This study mainly focuses on the occurrences of both flow separation on the boundary layer under an adverse pressure gradient and subsequent hydraulic jump with the abrupt rising of free surface during rundown motion of the shoaling wave, together with emphasis on the evolution of vortex structures underlying the separated shear layer and hydraulic jump. A flow visualization technique with particle trajectory method and a high-speed particle image velocimetry (HSPIV) system with a high-speed digital camera were used. Based on the instantaneous flow images visualized and/or the ensemble-averaged velocity fields measured, the following interesting features, which are unknown up-to-date, are presented and discussed in this study: (1) Flow bifurcation occurring on both offshore and onshore sides of the explicit demarcation curve and the stagnation point during runup motion; (2) The dependence of the diffuser-like flow field, being changed from the supercritical flow in the shallower region to the subcritical flow in the deeper counterpart, on the Froude number during the early and middle stages of rundown motion; (3) The positions and times for the occurrences of the incipient flow separation and the sudden rising of free surface of the hydraulic jump; (4) The associated movement and evolution of vortex structures under the separated shear layer, the hydraulic jump and/or the high-speed external main stream of the retreated flow; and (5) The entrainment of air bubbles from the free surface into the external main stream of the retreated flow.


2003 ◽  
Vol 478 ◽  
pp. 197-226 ◽  
Author(s):  
M. J. VOGEL ◽  
A. H. HIRSA ◽  
J. M. LOPEZ

The flow in a rectangular cavity driven by the sinusoidal motion of the floor in its own plane has been studied both experimentally and computationally over a broad range of parameters. The stability limits of the time-periodic two-dimensional base state are of primary interest in the present study, as it is within these limits that the flow can be used as a viable surface viscometer (as outlined theoretically in Lopez & Hirsa 2001). Three flow regimes have been found experimentally in the parameter space considered: an essentially two-dimensional time-periodic flow, a time-periodic three-dimensional flow with a cellular structure in the spanwise direction, and a three-dimensional irregular (in both space and time) flow. The system poses a space–time symmetry that consists of a reflection about the vertical mid-plane together with a half-period translation in time (RT symmetry); the two-dimensional base state is invariant to this symmetry. Computations of the two-dimensional Navier–Stokes equations agree with experimentally measured velocity and vorticity to within experimental uncertainty in parameter regimes where the flow is essentially uniform in the spanwise direction, indicating that in this cavity with large spanwise aspect ratio, endwall effects are small and localized for these cases. Two classes of flows have been investigated, one with a rigid no-slip top and the other with a free surface. The basic states of these two cases are quite similar, but the free-surface case breaks RT symmetry at lower forcing amplitudes, and the structure of the three-dimensional states also differs significantly between the two classes.


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.


2010 ◽  
Vol 665 ◽  
pp. 238-273 ◽  
Author(s):  
CHANG-YUE XU ◽  
LI-WEI CHEN ◽  
XI-YUN LU

Numerical investigation of the compressible flow past a wavy cylinder was carried out using large-eddy simulation for a free-stream Mach number M∞ = 0.75 and a Reynolds number based on the mean diameter Re = 2 × 105. The flow past a corresponding circular cylinder was also calculated for comparison and validation against experimental data. Various fundamental mechanisms dictating the intricate flow phenomena, including drag reduction and fluctuating force suppression, shock and shocklet elimination, and three-dimensional separation and separated shear-layer instability, have been studied systematically. Because of the passive control of the flow over a wavy cylinder, the mean drag coefficient of the wavy cylinder is less than that of the circular cylinder with a drag reduction up to 26%, and the fluctuating force coefficients are significantly suppressed to be nearly zero. The vortical structures near the base region of the wavy cylinder are much less vigorous than those of the circular cylinder. The three-dimensional shear-layer shed from the wavy cylinder is more stable than that from the circular cylinder. The vortex roll up of the shear layer from the wavy cylinder is delayed to a further downstream location, leading to a higher-base-pressure distribution. The spanwise pressure gradient and the baroclinic effect play an important role in generating an oblique vortical perturbation at the separated shear layer, which may moderate the increase of the fluctuations at the shear layer and reduce the growth rate of the shear layer. The analysis of the convective Mach number indicates that the instability processes in the shear-layer evolution are derived from oblique modes and bi-dimensional instability modes and their competition. The two-layer structures of the shear layer are captured using the instantaneous Lamb vector divergence, and the underlying dynamical processes associated with the drag reduction are clarified. Moreover, some phenomena relevant to the compressible effect, such as shock waves, shocklets and shock/turbulence interaction, are analysed. It is found that the shocks and shocklets which exist in the circular cylinder flow are eliminated for the wavy cylinder flow and the wavy surface provides an effective way of shock control. As the shock/turbulence interaction is avoided, a significant drop of the turbulent fluctuations around the wavy cylinder occurs. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the passive control of the compressible flow past a wavy surface.


Author(s):  
Qingfang Jiang

AbstractThe impact of Kelvin-Helmholtz billows (KHBs) in an elevated shear layer (ESL) on the underlying atmospheric boundary layer (BL) is examined utilizing a group of large-eddy simulations. In these simulations, KHBs develop in the ESL and experience exponential growth, saturation, and exponential decay stages. In response, strong wavy motion occurs in the BL, inducing rotor circulations near the surface when the BL is stable. During the saturation stage, secondary instability develops in the ESL and the wavy BL almost simultaneously, followed by the breakdown of the quasi-two-dimensional KH billows and BL waves into three-dimensional turbulence. Consequently, during and after a KH event, the underlying BL becomes more turbulent with its depth increased and stratification weakened substantially, suggestive of significant lasting impact of elevated KH billows on the atmospheric BL. The eventual impact of KHBs on the BL is found to be sensitive to both the ESL and BL characteristics.


Sign in / Sign up

Export Citation Format

Share Document