Dynamical mode decomposition of Gurney flap wake flow

The present work uses dynamic mode decomposition (DMD) to analyze wake flow of NACA0015 airfoil with Gurney flap. The physics of DMD is first introduced. Then the PIV-measured wake flow velocity field is decomposed into dynamical modes. The vortex shedding pattern behind the trailing edge and its high-order harmonics have been captured with abundant information such as frequency, wavelength and convection speed. It is observed that high-order dynamic modes convect faster than low-order modes; moreover the wavelength of the dynamic modes scales with the corresponding frequency in power law.     

Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

Interfacial instability in turbulent flow over a liquid film in a channel

We revisit the stability of a deformable interface that separates a fully-developed turbulent gas flow from a thin layer of laminar liquid. Although this problem has received considerable attention previously, a model that requires no fitting parameters and that uses a base-state profile that has been validated against experiments is, as yet, unavailable. Furthermore, the significance of wave-induced perturbations in turbulent stresses remains unclear. To address these outstanding issues, we investigate this problem and introduce a turbulent base-state velocity that requires specification of a flow rate or a pressure drop only; no adjustable parameters are necessary. This base state is validated extensively against available experimental data as well as the results of direct numerical simulations. In addition, the effect of perturbations in the turbulent stress distributions is investigated, and demonstrated to be small for cases wherein the liquid layer is thin. The detailed modelling of the liquid layer also elicits two unstable modes, ‘interfacial’ and ‘internal’, with the former being the more dominant of the two. We show that it is possible for interfacial roughness to reduce the growth rate of the interfacial mode in relation to that of the internal one, promoting the latter, to the status of most dangerous mode. Additionally, we introduce an approximate measure to distinguish between ‘slow’ and ‘fast’ waves, the latter being the case for ‘critical-layer’-induced instabilities; we demonstrate that for the parameter ranges studied, the large majority of the waves are ‘slow’. Finally, comparisons of our linear stability predictions are made with experimental data in terms of critical parameters for onset of wave-formation, wave speeds and wavelengths; these yield agreement within the bounds of experimental error.     

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.     

A proper orthogonal decomposition of a simulated supersonic shear layer

The Karhunen–Loève procedure is applied to the analysis of an ensemble of snapshots obtained from a conditionally sampled localized shear layer simulation. The computed set of optimal basis functions is used to economically characterize sampled flow realizations. Pictorially it is seen that the essential features (and roughly 80% of the energy) of typical flows are captured by retaining roughly 10–20 parameters in the expansion. Smaller-scale features are resolved by retaining more terms in the series.     

Examination of wall damping for the k-ε turbulence model using direct simulations of turbulent channel flow

Handler, Hendricks and Leighton have recently reported results for the direct numerical simulation (DNS) of a turbulent channel flow at moderate Reynolds number. These data are used to evaluate the terms in the exact and modelled transport equations for the turbulence kinetic energy k and the isotropic dissipation function ε. Both modelled transport equations show significant imbalances in the high-shear region near the channel walls. The model for the eddy viscosity is found to yield distributions for the production terms which do not agree well with the distributions calculated from the DNS data. The source of the imbalance is attributed to the wall-damping function required in eddy viscosity models for turbulent flows near walls. Several models for the damping function are examined, and it is found that the models do not vary across the channel as does the damping when evaluated from the DNS data. The Lam-Bremhorst model and the standard van Driest model are found to give reasonable agreement with the DNS data. Modification of the van Driest model to include an effective origin yields very good agreement between the modelled production and the production calculated from the DNS data, and the imbalance in the modelled transport equations is significantly reduced.     

Modelling a recirculating density-driven turbulent flow

A mathematical model of turbulent density-driven flows is presented and is solved numerically. A form of the k–? turbulence model is used to characterize the turbulent transport, and both this non-linear model and a sediment transport equation are coupled with the mean-flow fluid motion equations. A partitioned, Newton–Raphson-based solution scheme is used to effect a solution. The model is applied to the study of flow through a circular secondary sedimentation basin.     

Linear stability of channel entrance flow

The spatial stability of two dimensional, steady channel flow is investigated in the downstream entry zone for both exponentially and algebraically growing disturbances. A model based on previous work is presented for the base flow which represents a small deformation of plane Poiseuille flow. The base flow evolution towards the fully developed state comes from the experimental and theoretical study of M. Asai and J.M. Floryan [M. Asai, J.M. Floryan, Certain aspects of channel entrance flow, Phys. Fluids 16 (2004) 1160–1163]. This flow is found to be more stable than the parabolic Poiseuille flow. The most destabilizing base flow defect is then calculated using a variational method. The compromise between the destabilizing effect of the defect, which diffuses downstream, and the instability growth is found to be insufficient to provoke transition in the downstream laminar flow.     

Detached direct numerical simulations of turbulent two-phase bubbly channel flow

DNS simulations of two-phase turbulent bubbly channel flow at Re_τ = 180 (Reynolds number based on friction velocity and channel half-width) were performed using a stabilized finite element method (FEM) and a level set approach to track the air/water interfaces.     

Extensional behavior influence on viscoelastic turbulent channel flow

In this work we use in the simulation of a viscoelastic turbulent channel flow a modification of the finitely extensible of non-linear elastic dumbbells with the Peterlin approximation (FENE-P) constitutive model for dilute polymer solutions, applicable to high extensional deformations. The new feature introduced by this modification is that the free energy of the polymer (since it is assumed to be entirely entropically driven) remains always bounded (FENE-PB). The characteristics of the model under steady shear flow, pure elongational flow and transient extensional behavior are presented. It is found that the FENE-PB model is more shear thinning than FENE-P. Most importantly, it also shows a higher extensional viscosity than the FENE-P model. Although the steady-state Trouton ratio asymptotically reaches at high extensional rates the same limit as the FENE-P model, the transition from the Newtonian value is sharper and faster. We use the FENE-PB model in direct numerical simulations (DNS) of viscoelastic turbulent channel flow using spectral approximations. The results for various statistics of the flow and the polymer conformation, when compared against those obtained with the original FENE-P model and the same rheological parameters, show an enhanced polymer-induced drag reduction effect and enhanced deformation of the polymer molecules. This indicates that it is not only the asymptotic but also details from the extensional rheological behavior that matter in quantitatively specifying turbulent viscoelastic flow behavior.     

Hierarchical structures in a turbulent pipe flow

A hierarchical structure (HS) analysis (β-test and γ-test) is applied to a fully developed turbulent pipe flow. Velocity signals are measured at two cross sections in the pipe and at a series of radial locations from the pipe wall. Particular attention is paid to the variation of turbulent statistics at wall units 10<y⁺<3000. It is shown that at all locations the velocity fluctuations satisfy the She–Leveque hierarchical symmetry (Phys. Rev. Lett. 72 (1994) 336). The measured HS parameters, β and γ, are interpreted in terms of the variation of fluid structures. Intense anisotropic fluid structures generated near the wall appear to be more singular than the most intermittent structures in isotropic turbulence and appear to be more outstanding compared to the background fluctuations; this yields a more intermittent velocity signal with smaller γ and β. As turbulence migrates into the logarithmic region, small-scale motions are generated by an energy cascade and large-scale organized structures emerge which are also less singular than the most intermittent structures of isotropic turbulence. At the center, turbulence is nearly isotropic, and β and γ are close to the 1994 She–Leveque predictions. A transition is observed from the logarithmic region to the center in which γ drops and the large-scale organized structures break down. We speculate that it is due to the growing eddy viscosity effects of widely spread turbulent fluctuations in a similar way as in the breakdown of the Taylor vortices in a turbulent Couette–Taylor flow at high Reynolds numbers.     

Temporal large-eddy simulation of unstratified and stably stratified turbulent channel flows

Recently, Pruett et al. [Pruett, C.D., Gatski, T.B., Grosch, C.E., Thacker, W.D., 2003. The temporally filtered Navier–Stokes equations: properties of the residual stress. Phys. Fluids 15, 2127–2140] proposed an approach to large-eddy simulation (LES) based on time-domain filtering; their approach was termed temporal large-eddy simulation or TLES. In a continuation of their work, Pruett and collaborators tested their methodology by successfully performing TLES of unstratified turbulent channel flow up to Reynolds number of 590 (based on channel half-height and friction velocity) [Pruett, C.D., Thomas, B.C., Grosch, C.E., Gatski, T.B., 2006. A temporal approximate deconvolution model for LES. Phys. Fluids 18, 028104, 4p]. Here, we carefully analyze the TLES methodology in order to understand the role of its key components and in the process compare TLES to more traditional approaches of spatial LES. Furthermore, we extend the methodology to stably stratified turbulent channel flow.     

Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow

The present study reports detailed statistics for velocity and transfer rates of heavy particles dispersed in turbulent boundary layers. Statistics have been extracted from a homogeneous source of data covering a large target parameter space and are used here to analyze the effects of gravity and lift on particle dispersion and deposition in a systematic way. Datasets were obtained performing Direct Numerical Simulation (DNS) of particle-laden turbulent upward/downward flow in a vertical channel. Six values for the particle timescale (the particle Stokes number, St) ranging three orders of magnitude were considered to analyze the deposition process as the controlling mechanism was shifting from turbulent diffusion to inertia-moderated crossing trajectories. For the particle timescales examined, gravity and lift do not influence the qualitative behavior of particles even though velocity profiles and deposition coefficients are modified in a non-monotonic fashion, reaching an optimum for St ? 15. Physical mechanisms for the different behavior are discussed. Raw data and statistics obtained from the present DNS are made available at http://cfd.cineca.it (mirror site: http://158.110.32.35/download/database) and can be used to test simple models and closure equations for multiphase RANS and Large Eddy simulations.     

Flow structure of microbubble-laden turbulent channel flow measured by PIV combined with the shadow image technique

The turbulence structure of a horizontal channel flow with microbubbles is experimentally investigated using combined particle image velocimetry (PIV) in order to clarify the mechanism of drag reduction caused by microbubbles. A new system which simultaneously measures the liquid phase and the dispersed bubbles is proposed, based on a combination of particle tracking velocimetry (PTV), laser-induced fluorescence (LIF) and the shadow image technique (SIT). To accurately obtain the velocity of the liquid phase, tracer particles which overlap with the bubble shadow images are almost entirely eliminated in the post-processing. Finally, the turbulence characteristics of the flow field are presented, including measurements for both phases, and the bubble effect on the turbulence is quantified.     

Effects of blowing through a porous strip in a turbulent channel flow

An experimental investigation is performed on a fully developed turbulent channel flow with local injection through a porous strip. The Reynolds number based on the channel half-width was set to 5000. In addition to the no blowing case, measurements are made for three different blowing rates σ = 0.22, 0.36 and 0.58 (where σ is the ratio of momentum flux gain due to the blowing and momentum flux of the incoming channel flow). Measurements carried out with hot-wire anemometry reveal that injection strongly affects both the velocity profiles and the turbulence characteristics. The injection decreases the skin friction coefficient and increases all the Reynolds stresses downstream the blowing strip. The turbulence structure and the bursting phenomena were examined using space-time correlations measurements and conditional analysis. It is found that the injection increases the frequency of occurrence of the bursts.     

Modeling fluid-particle interaction in dilute-phase turbulent liquid-particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid-particle pipe flows in which the inter-stitial fluid effect on the particle fluctuating motion is significant.The impacts of employing different drag correlations and turbulence closure models to describe the fluid-particle interactions(i.e.drag force and long-range interaction)are examined at both the mean and fluctuating velocity levels.The model pre-dictions are validated using experimental data of turbulent liquid-particle flows in a vertical pipe at different particle Reynolds numbers(ReP > 400 and ReP < 400),which characterize the importance of the vortex shedding phenomenon in the fluid-phase turbulence modulation.The results indicate that(1)the fluctuating velocity level predictions at different ReP are highly sensitive to the drag correlation selec-tion and(2)different turbulence closure models must be employed to accurately describe the long-range fluid-particle interaction in each phase.In general,good agreement is found between the model predic-tions and the experimental data at both the mean and fluctuating velocity levels provided that appropriate combinations of the drag correlation and the turbulence closure model are selected depending on Rep.     

Wake response of a stationary finite-sized particle in a turbulent channel flow

Here we consider the effect of a finite-sized stationary particle in a channel flow of modest turbulence at _Reτ=178.12${\mathit{Re}}_{\tau }=178.12$. The size of particle is varied such that the particle Reynolds number ranges from about 40 to 450. The location of the particle is chosen to be either in the buffer layer (_yp+=17.81)$(y_{p+}=17.81)$ or at the channel center. Fully resolved direct numerical simulations of the turbulent channel flow around the particles is performed. Here the ambient turbulence intensity relative to the mean velocity seen by the particle is large (I=23.16%)$(I=23.16\%)$ in the buffer region, while it is substantially lower (I=4.09%)$(I=4.09\%)$ at the channel center. We present results on turbulence modulation due to the particle in terms of wake dynamics and vortex shedding.     

Orthogonal wavelet analysis of counter gradient transport phenomena in turbulent asymmetric channel flow

In this paper four families of orthogonal wavelets are applied to analyze the turbulent counter gradient transport phenomena in fully developed asymmetric channel flows. The results show that: (1) In the instance of counter gradient transport, the principal scale of the coherent structure is responsible for the strong local counter gradient transport; (2) Counter gradient transport phenomena have a strong effect on the intermittency of turbulence; (3) Non-Gaussian part of the principal coherent structure is essential for counter gradient transport phenomena.The project supported by the National Natural Science Foundation of China (10272071, 10472063)     

Open channel turbulent flow over hemispherical ribs

This paper reports an experimental investigation of open channel turbulent flow over hemispherical ribs. A row of ribs consists of hemispheres closely placed to one another in the spanwise direction and cover the entire span of the channel. The pitch-to-height ratio is varied to achieve the so-called d-type, intermediate and k-type roughness. The Reynolds numbers based on water depth, h, and momentum thickness, θ, of the approach flow are respectively, Re_h = 28,100 and Re_θ = 1800. A particle image velocimetry is used to obtain detailed velocity measurements in and above the cavity. Streamlines, mean velocity and time-averaged turbulent statistics are used to study the effects of pitch-to-height ratio on the flow characteristics and also to document similarities and differences between the present work and prior studies over two-dimensional transverse rods. It was observed that interaction between the outer flow and the shear layers generated by ribs is strongest for k-type and least for d-type ribs. The results also show that hemispherical ribs are less effective in augmenting flow resistance compared to two-dimensional transverse ribs. The levels of the Reynolds stresses and budget terms increase with increasing pitch-to-height ratio inside the roughness sublayer.