PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Numerical simulation of the flow around a simplified vehicle model with active flow control

Large-eddy simulation (LES) was used to study the influence and the resulting flow mechanisms of active flow control applied to a two-dimensional vehicle geometry. The LES results were validated against existing Particle Image Velocimetry (PIV) and force measurement data. This was followed by an exploration of the influence of flow actuation on the near-wake flow and resulting aerodynamic forces. Not only was good agreement found with the previous experimental study, but new knowledge was gained in the form of a complex interaction of the actuation with the coherent flow structures. The resulting time-averaged flow shows a strong influence of the extension of the actuation slots and the lateral solid walls on the near-wake flow structures and thereby on the resulting drag.     

Skin friction reduction by large air bubbles in a horizontal channel flow

Microbubble and air film methods are believed to be applicable to skin friction reduction in ships. Small bubbles are dispersed into the turbulent boundary layer in the former case, and wide air layers cover the wall surface in the latter case. Previous studies did not specifically address the intermediate case between the microbubble and air film conditions. This study is concerned with the possibility and mechanism of drag reduction using relatively large air bubbles compared to the boundary layer thickness in a horizontal turbulent channel flow. The relationship between local skin friction and the bubble’s interfacial structure is investigated by synchronizing the measurement of wall-shear stress with the image acquisition of bubbles. The bubble sizes range from 2 to 90 mm approximately. As a result, a negative correlation between the local skin friction and the local void fraction is confirmed by the time-resolved measurement. A new observation is the fact that the local skin friction decreases drastically in the rear part of individual large bubbles, and rapidly increases after the bubble’s rear interface passes. This characteristic underlies the bubble-size dependency of the average skin friction in the intermediate bubble size condition.     

Numerical prediction of locally forced turbulent boundary layer

An unsteady numerical simulation was performed to analyze flow structure behind a local suction/blowing in a flat-plate turbulent boundary layer. The local forcing was given to the boundary layer flow by means of a sinusoidally oscillating jet. A version of the unsteady k––f_μ model [Fluid Dyn. Res. 26 (6) (2000) 421] was employed. The Reynolds number based on the momentum thickness was about Re_θ=1700. The forcing frequency was varied in the range 0.011f⁺0.044 with a fixed forcing amplitude A_o=0.4. The predicted results were compared and validated with the experimental data. It was shown that the unsteady locally forced boundary layer flow is predicted well by the k––f_μ model. The time-dependent numerical flow visualizations were demonstrated during one period of the local forcing. The effect of the pitch angle of local forcing on the reduction of skin friction was examined.     

Subgrid Scale Modeling in Large-Eddy Simulation of Turbulent Combustion Using Premixed Flamelet Chemistry

Large-eddy simulation (LES) of turbulent combustion with premixed flamelets is investigated in this paper. The approach solves the filtered Navier–Stokes equations supplemented with two transport equations, one for the mixture fraction and another for a progress variable. The LES premixed flamelet approach is tested for two flows: a premixed preheated Bunsen flame and a partially premixed diffusion flame (Sandia Flame D). In the first case, we compare the LES with a direct numerical simulation (DNS). Four non-trivial models for the chemical source term are considered for the Bunsen flame: the standard presumed beta-pdf model, and three new propositions (simpler than the beta-pdf model): the filtered flamelet model, the shift-filter model and the shift-inversion model. A priori and a posteriori tests are performed for these subgrid reaction models. In the present preheated Bunsen flame, the filtered flamelet model gives the best results in a priori tests. The LES tests for the Bunsen flame are limited to a case in which the filter width is only slightly larger than the flame thickness. According to the a posteriori tests the three models (beta-pdf, filtered flamelet and shift-inversion) show more or less the same results as the trivial model, in which subgrid reaction effects are ignored, while the shift-filter model leads to worse results. Since LES needs to resolve the large turbulent eddies, the LES filter width is bounded by a maximum. For the present Bunsen flame this means that the filter width should be of the order of the flame thickness or smaller. In this regime, the effects of subgrid reaction and subgrid flame wrinkling turn out to be quite modest. The LES-results of the second case (Sandia Flame D) are compared to experimental data. Satisfactory agreement is obtained for the main species. Comparison is made between different eddy-viscosity models for the subgrid turbulence, and the Smagorinsky eddy-viscosity is found to give worse results than eddy-viscosities that are not dominated by the mean shear. Paper presented on the Eccomas Thematic Conference Computational Combustion 2007, submitted for a special issue of Flow, Turbulence and Combustion.     

A Similarity Subgrid Model for Premixed Turbulent Combustion

The accuracy of large-eddy simulation (LES) of a turbulent premixed Bunsen flame is investigated in this paper. To distinguish between discretization and modeling errors, multiple LES, using different grid sizes h but the same filterwidth Δ, are compared with the direct numerical simulation (DNS). In addition, LES using various values of Δ but the same ratio Δ/h are compared. The chemistry in the LES and DNS is parametrized with the standard steady premixed flamelet for stochiometric methane-air combustion. The subgrid terms are closed with an eddy-viscosity or eddy-diffusivity approach, with an exception of the dominant subgrid term, which is the subgrid part of the chemical source term. The latter subgrid contribution is modeled by a similarity model based upon 2Δ, which is found to be superior to such a model based upon Δ. Using the 2Δ similarity model for the subgrid chemistry the LES produces good results, certainly in view of the fact that the LES is completely wrong if the subgrid chemistry model is omitted. The grid refinements of the LES show that the results for Δ = h do depend on the numerical scheme, much more than for h = Δ/2 and h = Δ/4. Nevertheless, modeling errors and discretization error may partially cancel each other; occasionally the Δ = h results were more accurate than the h ≤ Δ results. Finally, for this flame LES results obtained with the present similarity model are shown to be slightly better than those obtained with standard β-pdf closure for the subgrid chemistry.     

The turbulent character and pressure loss produced by periodic symmetric ribs in a circular duct

The structural character and steady-state statistics of the turbulence inside a rib-wall circular duct is investigated by the large-eddy simulation (LES) methodology. The impetus of this study is to gain an understanding of the principle physics attributing to minimizing the pressure recovered (or maximizing the pressure loss) within the core flow. For a rib periodicity with height (h) to pitch (p) ratio p/h=5, the computational results show that the majority of turbulence produced due to the rib’s presence is concentrated near the rib crest leading edge. Pairs of counter-rotating streamwise vortices form soon after the leading edge that are quickly convected radially toward the core flow. The turbulent activity within the duct trough region is negligible compared to the turbulence levels of the core flow. At this rib periodicity, the separated shear layers from the trailing edge of each rib nearly reattach to the trough floor before reaching the next rib. The resultant irrecoverable pressure loss in the form a centerline frictional coefficient is verified by an ‘at-sea’ test on board a US Navy submarine. Based on the duct diameter, their Reynolds numbers are Re_D^LES=8×10³ and (Re_D^exp)_avg=4×10⁶, respectively.     

LES prediction of space-time correlations in turbulent shear flows

We compare the space-time correlations calculated from direct numerical simulation(DNS) and large-eddy simulation(LES) of turbulent channel flows.It is found from the comparisons that the LES with an eddy-viscosity subgrid scale(SGS) model over-predicts the space-time correlations than the DNS.The overpredictions are further quantified by the integral scales of directional correlations and convection velocities.A physical argument for the overprediction is provided that the eddy-viscosity SGS model alone does not includes the backscatter effects although it correctly represents the energy dissipations of SGS motions.This argument is confirmed by the recently developed elliptic model for space-time correlations in turbulent shear flows.It suggests that enstrophy is crucial to the LES prediction of spacetime correlations.The random forcing models and stochastic SGS models are proposed to overcome the overpredictions on space-time correlations.     

Frontogenesis in turbulent flow through porous media using non-linear Forchheimer equation

We present here both one- and two-dimensional models for turbulent flow through heterogeneous unbounded fluid saturated porous media using non-linear Forchheimer extended Darcy (DF) equation in the presence of gravitational field. The fluid is initially at rest and sets in motion due to a uniform horizontal density gradient. It is shown that a purely horizontal motion develops satisfying non-linear DF equation. Analytical solutions of this non-linear Initial Value Problem are obtained and limiting solutions valid for the Darcy regime in the case of laminar flow are derived. A measure of the stability of the flow is discussed briefly using Richardson number. The comparison between the nature of the solutions satisfying the non-linear and linear initial value problems are made. We found that even in the case of turbulent flow the vertical density gradient varies continuously both with space z and time t but the horizontal density gradient remains unchanged. The existence and uniqueness theorem of the Initial Value Problem is proved. The stability of these solutions are discussed and it is shown that the solutions are qualitatively and quantitatively different for and in the upper and lower half of the region. In particular, we have shown that the solution which is stable for infinitesimal perturbations is also stable for arbitrary perturbations both in time and space.In the case of two-dimensional motion, a piecewise initial density gradient with continuous distribution of density, stream function formulation is used and the solutions are obtained using time-series analysis. In this case solution shows crowding of the density profiles in the lower-half of the channel reflecting an increase in density gradient and incipient of frontogenesis there, because of the increase in circulation of the flow due to piecewise initial density gradient.     

Multiphase drag reduction: Effect of eliminating slugs

Two-phase pressure drop measurements were taken for air/water mixtures in a 0.052-m diameter horizontal pipe with special focus on the superficial liquid velocity range of 0.03–1.2 m/s at superficial gas velocities of 3.8, 5.2, and 6.6 m/s. It was found that the addition of 400 ppm of sodium dodecyl sulfate (SDS) to the water reduced the pressure drop by 25–40% when compared to equal flow rates without SDS. The pressure drop reduction occurred where the SDS eliminated the occurrence of the intermittent flow present with water. It was also found that the same concentration of SDS had virtually no effect on single phase liquid pressure drop. The pressure drop reduction appears to be due solely to the suppression of intermittent flow patterns.