Spectral Dynamics in a Channel Flow With Microbubbles

Author(s):  
Jose Alfredo Jimenez Bernal ◽  
Yassin A. Hassan ◽  
Carlos E. Estrada Perez ◽  
Claudia del Carmen Gutierrez-Torres

Measurements were carried out near the upper wall of a turbulent water channel flow at Reynolds number equal to 10256 (considering total height of the channel and the bulk velocity). Particle Image Velocimetry (PIV) measurement technique was used to record instantaneous velocity fields in the x-y plane. The evaluation of the one dimensional spectra for the fluctuating streamwise velocity was performed in the streamwise wavenumber direction (k1) with and without microbubbles in the boundary layer. A redistribution of the energy from larger wavenumber to lower wavenumber was observed for the spectra close to the wall. However, the opposite trend was obtained in the outer region. The energy content with microbubbles was less than the energy without microbubbles in the near wall region of the boundary layer.

2010 ◽  
Vol 664 ◽  
pp. 33-50 ◽  
Author(s):  
M. H. KRANE ◽  
L. M. GREGA ◽  
T. WEI

Measurements of the near-wall velocity field of the flow over cylinders aligned with a uniform flow are presented. The broader objective of this investigation was to quantify and understand the role of transverse curvature in the limit as cylinder diameter approaches zero. The specific goal was to begin with a turbulent boundary layer over a larger radius cylinder and see what happens as the radius is reduced. Spatially and temporally resolved digital particle image velocimetry (DPIV) measurements were made on three different radius cylinders, 0.14 cm ≤ a ≤ 3.05 cm, extending along the length of a large free-surface water tunnel. Mean and fluctuating profiles are presented at a fixed streamwise location and free-stream speed. For the first time, spatially resolved measurements were made very close to the wall, permitting direct determination of wall shear stress, i.e. uτ, from near-wall velocity profiles. The measurements revealed a region close to the wall for small radii where the mean streamwise velocity profile is inflectional. This has significant implications on assumptions regarding what happens in the limit of a vanishing cylinder radius.


2010 ◽  
Vol 132 (4) ◽  
Author(s):  
Bérengère Podvin ◽  
Yann Fraigneau ◽  
Julien Jouanguy ◽  
J.-P. Laval

We use proper orthogonal decomposition (POD) to estimate the flow in the near-wall region based on information from the outer buffer layer. Our goal is to assess how the flow structures in the inner wall region are connected to those further away from the wall, and to investigate the nature of the coupling between the inner and the outer region in the POD framework. Reconstructions are carried out for numerical simulations of a plane channel flow at two different Reynolds numbers. We show that elongated structures with a spanwise wavelength smaller than a critical value tend to be concentrated in the inner layer. The critical wavelength is shown to scale with the inner layer height, and interactions between the inner and the outer layer appear to take place predominantly over a self-similar, height-dependent, range of wavenumbers, in agreement with Townsend’s attached eddy hypothesis. The reconstructed field appears to capture an adequate energy content and to remain correlated with the real field even close to the wall, which reflects the persistence of energetic structures over the extent of the buffer layer.


1996 ◽  
Vol 310 ◽  
pp. 269-292 ◽  
Author(s):  
Hugh M. Blackburn ◽  
Nagi N. Mansour ◽  
Brian J. Cantwell

An investigation of topological features of the velocity gradient field of turbulent channel flow has been carried out using results from a direct numerical simulation for which the Reynolds number based on the channel half-width and the centreline velocity was 7860. Plots of the joint probability density functions of the invariants of the rate of strain and velocity gradient tensors indicated that away from the wall region, the fine-scale motions in the flow have many characteristics in common with a variety of other turbulent and transitional flows: the intermediate principal strain rate tended to be positive at sites of high viscous dissipation of kinetic energy, while the invariants of the velocity gradient tensor showed that a preference existed for stable focus/stretching and unstable node/saddle/saddle topologies. Visualization of regions in the flow with stable focus/stretching topologies revealed arrays of discrete downstream-leaning flow structures which originated near the wall and penetrated into the outer region of the flow. In all regions of the flow, there was a strong preference for the vorticity to be aligned with the intermediate principal strain rate direction, with the effect increasing near the walls in response to boundary conditions.


Author(s):  
Abolfazl Shiri ◽  
William K. George

The turbulence natural convection boundary layer inside a infinite vertical channel with differentially heated walls is analyzed based on a similarity solution methodology. The differences between mean temperature and velocity profiles in a boundary layer along a vertical flat plate and in a channel flow, make it necessary to introduce new sets of scaling parameters. In the limit as H* → ∞, two distinctive parts are considered: an outer region which dominates the core of the flow and inner constant heat flux region close to the walls. The proper inner scaling velocity is showed to be determined by the outer parameters due to momentum integral. The theory is contrasted with the one suggested by George & Capp (1), the deficiencies of which are identified.


2019 ◽  
Vol 36 (1) ◽  
pp. 12-21 ◽  
Author(s):  
Xiaotong Cui ◽  
Nan Jiang ◽  
Xiaobo Zheng ◽  
Zhanqi Tang

Abstract This study experimentally investigates the impact of a single piezoelectric (PZT) actuator on a turbulent boundary layer from a statistical viewpoint. The working conditions of the actuator include a range of frequencies and amplitudes. The streamwise velocity signals in the turbulent boundary layer flow are measured downstream of the actuator using a hot-wire anemometer. The mean velocity profiles and other basic parameters are reported. Spectra results obtained by discrete wavelet decomposition indicate that the PZT vibration primarily influences the near-wall region. The turbulent intensities at different scales suggest that the actuator redistributes the near-wall turbulent energy. The skewness and flatness distributions show that the actuator effectively alters the sweep events and reduces intermittency at smaller scales. Moreover, under the impact of the PZT actuator, the symmetry of vibration scales’ velocity signals is promoted and the structural composition appears in an orderly manner. Probability distribution function results indicate that perturbation causes the fluctuations in vibration scales and smaller scales with high intensity and low intermittency. Based on the flatness factor, the bursting process is also detected. The vibrations reduce the relative intensities of the burst events, indicating that the streamwise vortices in the buffer layer experience direct interference due to the PZT control.


2019 ◽  
Vol 872 ◽  
pp. 198-217 ◽  
Author(s):  
Duosi Fan ◽  
Jinglei Xu ◽  
Matthew X. Yao ◽  
Jean-Pierre Hickey

A novel approach to identify internal interfacial layers, or IILs, in wall-bounded turbulent flows is proposed. Using a fuzzy cluster method (FCM) on the streamwise velocity component, a unique and unambiguous grouping of the uniform momentum zones (UMZs) is achieved, thus allowing the identification of the IILs. The approach overcomes some of the key limitations of the histogram-based IIL identification methods. The method is insensitive to the streamwise domain length, can be used on inhomogeneous grids, uses all the available flow field data, is trivially extended to three dimensions and does not need user-defined parameters (e.g. number of bins) other than the number of zones. The number of zones for a given snapshot can be automatically determined by an a priori algorithm based on a kernel density estimation algorithm, or KDE. This automated approach is applied to compute the average number of UMZs as a function of Reynolds number $Re_{\unicode[STIX]{x1D70F}}$ in turbulent channel flows in several numerical simulations. This systematic approach reveals a dependence of the Reynolds number on the average number of UMZs in the channel flow; this supports previously reported observations in the boundary layer. The fuzzy clustering approach is applied to the turbulent boundary layer (experimental, planar particle image velocimetry) and channel flow (numerical, direct numerical simulation) at varying Reynolds numbers. The interfacial layers are characterized by a strong concentration of spanwise vorticity, with the outer-most layer located at the upper edge of the log layer. The three-dimensional interface identification reveals a streak-like organization. The large-scale motion (LSM) at the outer region of the channel flow boundary layer modulates the outer IIL. The corrugations of the outer IIL are aligned with the LSM and the conditional correlation of the inner and outer IIL height shows that extreme near-wall events leave their mark on the outer IIL corrugations.


1992 ◽  
Vol 114 (3) ◽  
pp. 598-606 ◽  
Author(s):  
N. Kasagi ◽  
Y. Tomita ◽  
A. Kuroda

A direct numerical simulation (DNS) of the fully developed thermal field in a two-dimensional turbulent channel flow of air was carried out. The isoflux condition was imposed on the two walls so that the local mean temperature increased linearly in the streamwise direction. With any buoyancy effect neglected, temperature was considered as a passive scalar. The computation was executed on 1,589,248 grid points by using a spectral method. The statistics obtained were root-mean-square temperature fluctuations, turbulent heat fluxes, turbulent Prandtl number, and dissipation time scales. They agreed fairly well with existing experimental and numerical simulation data. Each term in the budget equations of temperature variance, its dissipation rate, and turbulent heat fluxes was also calculated. It was found that the temperature fluctuation θ′ was closely correlated with the streamwise velocity fluctuation u′, particularly in the near-wall region. Hence, the distribution of budget terms for the streamwise and wall-normal heat fluxes, u′θ′ and v′θ′, were very similar to those for the two Reynolds stress components, u′u′ and u′v′.


1976 ◽  
Vol 27 (3) ◽  
pp. 217-228 ◽  
Author(s):  
Noor Afzal ◽  
K P Singh

SummaryIn an axisymmetric turbulent boundary layer along a circular cylinder at constant pressure, measurements have been made of mean velocity profile and turbulence characteristics: longitudinal velocity fluctuations, Reynolds shear stress, transverse correlation and spectrum. It has been found that the qualitative behaviour of an axisymmetric turbulent boundary layer is similar to that of a two-dimensional boundary layer in the wall region, where as in the outer region the effects of transverse curvature are observed.


Author(s):  
Sylvain Aguinaga ◽  
Olivier Simonin ◽  
Jacques Bore´e ◽  
Vincent Herbert

The deposition rate of droplets is strongly linked to their interaction with the boundary layer turbulence. In “industrial simulations”, droplets dispersion is usually modeled using Lagrangian stochastic simulations based on Reynold Average Navier Stokes (RANS) fluid calculations. Wall functions are also used to bound the number of mesh cells in the near wall region. But they also reduce the description of the boundary layer and lead to bad predictions of the droplets deposition rate. This study presents channel flow simulations using wall functions and run with the CFD code Fluent. In such configurations, the stochastic model of Fluent failed to represent the so-called “diffusion-impaction” regime of deposition. The “Concentration Wall Boundary Layer” model presented in this paper has been developed to predict deposition in simulations using industrial meshes with refinement such as y* > 20. This model calculates the deposition rate using only the intrinsic properties of the particles and the turbulent kinetic energy of the fluid expressed at the top of the boundary layer. The data provided by wall functions are then sufficient to calculate the deposition rate. This model is turned into a “Lagrangian stochastic wall boundary condition model” for the commercial CFD code Fluent. Various simulations have shown that this model improves remarkably the deposition predictions in channel flow. The dependence on the boundary cell size and the channel flow mean velocity has been tested. This model draws interesting perspectives to model deposition in complex configurations without requiring prohibiting mesh sizes.


Author(s):  
Yanguang Long ◽  
Jinjun Wang ◽  
Chong Pan

The sharp but irregular interface that separates the instantaneous turbulent and irrotational flows is termed as the turbulent/non-turbulent interface (TNTI). TNTI can be widely observed in various types of flow, such as turbulent boundary layers, jets and combustion flame fronts. Due to its importance on the intermittency and entrainment process, TNTI has been widely explored in its geometry and dynamic properties (da Silva et al., 2014). Most of the studies focus on the TNTIs in smooth plane boundary layer, while few investigate the effects of wall shapes. However, the wall conditions in many engineering applications are complex and heterogeneous, which will induce large-scale heterogeneity (Barros and Christensen, 2014) and require further investigations. To shed new light on the intermittency and entrainment above complex surfaces, the TNTI over spanwise heterogeneity are investigated here with time-resolved stereoscopic PIV (TR-SPIV). The model and TR-SPIV experimental set-up are shown in Fig. 1. The experiments are conducted in the low-speed water channel at Beijing University of Aeronautics and Astronautics. The spanwise distance S between two adjacent ridges is S/(δ) = 1.35, where (δ) is the spanwise-averaged boundary layer thickness. This spanwise distance is selected to induced strong secondary vortices (Vanderwel and Ganapathisubramani, 2015; Wangsawijaya et al., 2020). The Reynolds number based on the streamwise location x is Rex = 7.2×105. The field of view is around 2S×1.8S, and is captured by two CMOS cameras (2048×2048 pixel) with sampling rate as 500Hz. The averaged resolution is about 8 pixels per Kolmogorov scale (calculated at y/(δ) = 0.6), which is high enough for TNTI-related research (Borrell and Jimenez, 2016). The ´TNTI is detected by the magnitude of local enstrophy ω2/2, and the threshold is selected to be the value where changing the threshold has the smallest influence on the TNTI-mean-height (Watanabe et al., 2018). The time-mean velocity and TNTI location are present in Fig.2(a). A pair of counter-rotating largescale secondary vortices (SVs) are induced over the ridge-type roughness. At the position where SVs induce upwash flow, a low-momentum pathway (LMP) can be observed, while the time-mean height of TNTI (yI) is brought higher. As a contrast, where downwash flow induces high-momentum pathway (HMP), (yI) is lower. TNTI properties are further discussed from two aspect. The geometry properties are firstly investigated. The fractal dimension of the TNTI keeps as 2.3 along the spanwise direction. This value is consistent with the result over smooth plate (Borrell and Jimenez, 2016; Wu et al., 2020) and riblets plates(Cui et al., 2019),´ which indicates that the wall shapes do not influence the multiscale properties of the TNTI. The streamwise wavelength of the TNTI (λI) is further obtained by calculating the streamwise pre-multiplied spectrum of the TNTI. It is found that at each spanwise location, λI is identical to the wavelength of streamwise velocity fluctuation at the TNTI mean height. This shows that the large-scale fluctuation of TNTI is controlled by the large-scale streamwise velocity fluctuation structures. Secondly, the p.d.f. of TNTI instantaneous height is investigated, as shown in Fig. 2(b). It can be observed that the p.d.f. of TNTI height above LMP shows a negative skewness, while the p.d.f. above HMP skews positively. A closer look at instantaneous structures shows that the skewness is attributed to the different probability of Q2/Q4 events in LMP and HMP.


Sign in / Sign up

Export Citation Format

Share Document