Designing Multi-Shape Dual-Gas Initial Conditions for the Study of Hydrodynamic Instabilities

Author(s):  
Marta Maria Rasteiro dos Santos ◽  
Yannick Bury ◽  
Stephane Jamme

Abstract The flow resulting from the rotation of a series of thin plates that initially separate two gases of different densities is analysed using Direct Numerical Simulations. The ninety degrees plates' rotation forms a vorticity shear layer and a density interface in between the tips of two neighbouring plates. Results of this study show that the shape of these layers strongly depends on the plate tip-based Reynolds number that can be varied thanks to a parametrisation of the plates' opening law. Different regimes are identified corresponding to single- or multi-mode initial interfaces, with or without the occurrence of starting vortices during the formation of the shear layer. The density interfaces resulting from this procedure are particularly well-suited to serve as initial conditions for the study of the Richtmyer-Meshkov instability-induced mixing. Results of this study also provide a description of vortex formation in stratified flows.

1988 ◽  
Vol 190 ◽  
pp. 491-512 ◽  
Author(s):  
M. F. Unal ◽  
D. Rockwell

Vortex shedding from a circular cylinder is examined over a tenfold range of Reynolds number, 440 ≤ Re ≤ 5040. The shear layer separating from the cylinder shows, to varying degrees, an exponential variation of fluctuating kinetic energy with distance downstream of the cylinder. The characteristics of this unsteady shear layer are interpreted within the context of an absolute instability of the near wake. At the trailing-end of the cylinder, the fluctuation amplitude of the instability correlates well with previously measured values of mean base pressure. Moreover, this amplitude follows the visualized vortex formation length as Reynolds number varies. There is a drastic decrease in this near-wake fluctuation amplitude in the lower range of Reynolds number and a rapid increase at higher Reynolds number. These trends are addressed relative to the present, as well as previous, observations.


2014 ◽  
Vol 759 ◽  
pp. 579-611 ◽  
Author(s):  
Jae Hwan Jung ◽  
Hyun Sik Yoon

AbstractWe consider a twisted cylinder that was designed by rotating the elliptic cross-section along the spanwise direction, resulting in a passive control. The flow over the twisted cylinder is investigated at a subcritical Reynolds number (Re) of 3000 using large eddy simulation based on the finite volume method. For comparison, the flow past smooth and wavy cylinders is also calculated. The twisted cylinder achieves reductions of approximately 13 and 5 % in mean drag compared with smooth and wavy cylinders, respectively. In particular, the root mean square (r.m.s.) value of the lift fluctuation of the twisted cylinder shows a substantial decrease of approximately 96 % compared with the smooth cylinder. The shear layer of the twisted cylinder covering the recirculation region is more elongated than those of the smooth and wavy cylinders, and vortex shedding from the twisted cylinder is considerably suppressed. Consequently, the elongation of the shear layer from the body and the near disappearance of vortex shedding in the near wake with weak vortical strength contributes directly to the reduction of drag and lift oscillation. Various fundamental mechanisms that affect the flow phenomena, three-dimensional separation, pressure coefficient, vortex formation length and turbulent kinetic energy are examined systematically to demonstrate the effect of the twisted cylinder surface. In addition, for the twisted cylinder at $\mathit{Re}=3000$, the effect of the cross-sectional aspect ratio is investigated from 1.25 to 2.25 to find an optimal value that can reduce the drag and lift forces. Moreover, the effect of the Reynolds number on the aerodynamic characteristics is investigated in the range of $3\times 10^{3}\leqslant \mathit{Re}\leqslant 1\times 10^{4}$. We find that as Re increases, the mean drag and the r.m.s. lift coefficient of the twisted cylinder increase, and the vortex formation length decreases.


2007 ◽  
Vol 582 ◽  
pp. 103-132 ◽  
Author(s):  
J. R. CARPENTER ◽  
G. A. LAWRENCE ◽  
W. D. SMYTH

When a stably stratified density interface is embedded in a region of strong velocity shear, hydrodynamic instabilities result. Here we generalize the stratified shear layer to allow an offset between the centre of the shear layer and the density interface. By including this asymmetry, and keeping the density interface thin with respect to the shear layer, the asymmetric Holmboe (AH) instability emerges. This study examines the evolution and mixing behaviour of AH instabilities, and compares the results to the well-known Kelvin–Helmholtz (KH) and Holmboe instabilities. This is done by performing a series of direct numerical simulations (DNS). The simulation results show that there are two different mixing mechanisms present. The first is a feature of KH instabilities and leads to the mixing and production of intermediate density fluid. The second mixing mechanism is found in AH and Holmboe instabilities and consists of regions of mixing and turbulence production that are located on one or both sides of the density interface. Since the Holmboe-type instabilities do not generate a large-scale overturning of the central isopycnal, the density interface is able to retain its identity throughout the mixing event. The amount of mixing that takes place is found to be strongly dependent on the degree of asymmetry in the flow.


2019 ◽  
Vol 877 ◽  
pp. 35-81 ◽  
Author(s):  
Nek Sharan ◽  
Georgios Matheou ◽  
Paul E. Dimotakis

Aspects of turbulent shear-layer mixing are investigated over a range of shear-layer Reynolds numbers, $Re_{\unicode[STIX]{x1D6FF}}=\unicode[STIX]{x0394}U\unicode[STIX]{x1D6FF}/\unicode[STIX]{x1D708}$, based on the shear-layer free-stream velocity difference, $\unicode[STIX]{x0394}U$, and mixing-zone thickness, $\unicode[STIX]{x1D6FF}$, to probe the role of initial conditions in mixing stages and the evolution of the scalar-field probability density function (p.d.f.) and variance. Scalar transport is calculated for unity Schmidt numbers, approximating gas-phase diffusion. The study is based on direct-numerical simulation (DNS) and large-eddy simulation (LES), comparing different subgrid-scale (SGS) models for incompressible, uniform-density, temporally evolving forced shear-layer flows. Moderate-Reynolds-number DNS results help assess and validate LES SGS models in terms of scalar-spectrum and mixing estimates, as well as other metrics, to $Re_{\unicode[STIX]{x1D6FF}}\lesssim 3.3\times 10^{4}$. High-Reynolds-number LES investigations to $Re_{\unicode[STIX]{x1D6FF}}\lesssim 5\times 10^{5}$ help identify flow parameters and conditions that influence the evolution of scalar variance and p.d.f., e.g. marching versus non-marching. Initial conditions that generate shear flows with different mixing behaviour elucidate flow characteristics in each flow regime and identify elements that induce p.d.f. transition and scalar-variance behaviour. P.d.f. transition is found to be largely insensitive to local flow parameters, such as $Re_{\unicode[STIX]{x1D6FF}}$, or a previously proposed vortex-pairing parameter based on downstream distance, or other equivalent criteria. The present study also allows a quantitative comparison of LES SGS models in moderate- and high-$Re_{\unicode[STIX]{x1D6FF}}$ forced shear-layer flows.


2016 ◽  
Vol 801 ◽  
pp. 475-507 ◽  
Author(s):  
Md. Mahbub Alam ◽  
Honglei Bai ◽  
Yu Zhou

This work aims to provide a systematic experimental study of the wake behind two staggered square cylinders at a Reynolds number $Re=1.3\times 10^{4}$. Four distinct flow regimes, i.e. two single-street modes S-I and S-II and two double-street modes T-I and T-II, are identified based on extensive data, including Strouhal numbers $(St)$, flow structures and their downstream evolution. S-I, S-II and T-II are each further subdivided into two types. The flow characteristics in each regime are presented in terms of shear layer reattachment and impingement, vortex impingement, gap flow behaviour, interaction between cylinders and downstream evolution of the wake. A detailed discussion is made regarding how the physical aspects of the flow are connected to the initial conditions and the $St$ number.


2019 ◽  
Vol 871 ◽  
pp. 595-635 ◽  
Author(s):  
Mohammad Mohaghar ◽  
John Carter ◽  
Gokul Pathikonda ◽  
Devesh Ranjan

The effects of incident shock strength on the mixing transition in the Richtmyer–Meshkov instability (RMI) are experimentally investigated using simultaneous density–velocity measurements. This effort uses a shock with an incident Mach number of 1.9, in concert with previous work at Mach 1.55 (Mohaghar et al., J. Fluid Mech., vol. 831, 2017 pp. 779–825) where each case is followed by a reshock wave. Single- and multi-mode interfaces are used to quantify the effect of initial conditions on the evolution of the RMI. The interface between light and heavy gases ($\text{N}_{2}/\text{CO}_{2}$, Atwood number, $A\approx 0.22$; amplitude to wavelength ratio of 0.088) is created in an inclined shock tube at $80^{\circ }$ relative to the horizontal, resulting in a predominantly single-mode perturbation. To investigate the effects of initial perturbations on the mixing transition, a multi-mode inclined interface is also created via shear and buoyancy superposed on the dominant inclined perturbation. The evolution of mixing is investigated via the density fields by computing mixed mass and mixed-mass thickness, along with mixing width, mixedness and the density self-correlation (DSC). It is shown that the amount of mixing is dependent on both initial conditions and incident shock Mach number. Evolution of the density self-correlation is discussed and the relative importance of different DSC terms is shown through fields and spanwise-averaged profiles. The localized distribution of vorticity and the development of roll-up features in the flow are studied through the evolution of interface wrinkling and length of the interface edge, which indicate that the vorticity concentration shows a strong dependence on the Mach number. The contribution of different terms in the Favre-averaged Reynolds stress is shown, and while the mean density-velocity fluctuation correlation term, $\langle \unicode[STIX]{x1D70C}\rangle \langle u_{i}^{\prime }u_{j}^{\prime }\rangle$, is dominant, a high dependency on the initial condition and reshock is observed for the turbulent mass-flux term. Mixing transition is analysed through two criteria: the Reynolds number (Dimotakis, J. Fluid Mech., vol. 409, 2000, pp. 69–98) for mixing transition and Zhou (Phys. Plasmas, vol. 14 (8), 2007, 082701 for minimum state) and the time-dependent length scales (Robey et al., Phys. Plasmas, vol. 10 (3), 2003, 614622; Zhou et al., Phys. Rev. E, vol. 67 (5), 2003, 056305). The Reynolds number threshold is surpassed in all cases after reshock. In addition, the Reynolds number is around the threshold range for the multi-mode, high Mach number case ($M\sim 1.9$) before reshock. However, the time-dependent length-scale threshold is surpassed by all cases only at the latest time after reshock, while all cases at early times after reshock and the high Mach number case at the latest time before reshock fall around the threshold. The scaling analysis of the turbulent kinetic energy spectra after reshock at the latest time, at which mixing transition analysis suggests that an inertial range has formed, indicates power scaling of $-1.8\pm 0.05$ for the low Mach number case and $-2.1\pm 0.1$ for the higher Mach number case. This could possibly be related to the high anisotropy observed in this flow resulting from strong, large-scale streamwise fluctuations produced by large-scale shear.


1997 ◽  
Vol 119 (4) ◽  
pp. 810-817 ◽  
Author(s):  
C. Gau ◽  
W. Y. Sheu ◽  
C. H. Shen

Experiments are performed to study (a) slot air jet impingement cooling flow and (b) the heat transfer under acoustic excitations. Both flow visualization and spectral energy evolution measurements along the shear layer are made. The acoustic excitation at either inherent or noninherent frequencies can make the upstream shift for both the most unstable waves and the resulting vortex formation and its subsequent pairing processes. At inherent frequencies the most unstable wave can be amplified, which increases the turbulence intensity in both the shear layer and the core and enhances the heat transfer. Both the turbulence intensity and the heat transfer increase with increasing excitation pressure levels Spl until partial breakdown of the vortex occurs. At noninherent frequencies, however, the most unstable wave can be suppressed, which reduces the turbulence intensity and decreases the heat transfer. Both the turbulence intensity and the heat transfer decreases with increasing Spl, but increases with increasing Spl when the excitation frequency becomes dominant. For excitation at high Reynolds number with either inherent or noninherent frequency, a greater excitation pressure level is needed to cause the enhancement or the reduction in heat transfer. During the experiments, the inherent frequencies selected for excitation are Fo/2 and Fo/4, the noninherent frequencies are 0.71 Fo, 0.75 Fo, and 0.8 Fo, the acoustic pressure level varies from 70 dB to 100 dB, and the Reynolds number varies from 5500 to 22,000.


1994 ◽  
Vol 274 ◽  
pp. 1-21 ◽  
Author(s):  
R. Gopalkrishnan ◽  
M. S. Triantafyllou ◽  
G. S. Triantafyllou ◽  
D. Barrett

It is shown experimentally that free shear flows can be substantially altered through direct control of the large coherent vortices present in the flow.First, flow-visualization experiments are conducted in Kalliroscope fluid at Reynolds number 550. A foil is placed in the wake of a D-section cylinder, sufficiently far behind the cylinder so that it does not interfere with the vortex formation process. The foil performs a heaving and pitching oscillation at a frequency close to the Strouhal frequency of the cylinder, while cylinder and foil also move forward at constant speed. By varying the phase of the foil oscillation, three basic interaction modes are identified. (i) Formation of a street of pairs of counter-rotating vortices, each pair consisting of one vortex from the initial street of the cylinder and one vortex shed by the foil. The width of the wake is then substantially increased. (ii) Formation of a street of vortices with reduced or even reverse circulation compared to that of oncoming cylinder vortices, through repositioning of cylinder vortices by the foil and interaction with vorticity of the opposite sign shed from the trailing edge of the foil. (iii) Formation of a street of vortices with circulation increased through merging of cylinder vortices with vortices of the same sign shed by the foil. In modes (ii) and (iii) considerable repositioning of the cylinder vortices takes place immediately behind the foil, resulting in a regular or reverse Kármán street. The formation of these three interaction patterns is achieved only for specific parametric values; for different values of the parameters no dominant stable pattern emerges.Subsequently, the experiments are repeated in a different facility at larger scale, resulting in Reynolds number 20000, in order to obtain force and torque measurements. The purpose of the second set of experiments is to assess the impact of flow control on the efficiency of the oscillating foil, and hence investigate the possibility of energy extraction. It is found that the efficiency of the foil depends strongly on the phase difference between the oscillation of the foil and the arrival of cylinder vortices. Peaks in foil efficiency are associated with the formation of a street of weakened vortices and energy extraction by the foil from the vortices of the vortex street.


Sign in / Sign up

Export Citation Format

Share Document