scholarly journals Model-based design of transverse wall oscillations for turbulent drag reduction

2012 ◽  
Vol 707 ◽  
pp. 205-240 ◽  
Author(s):  
Rashad Moarref ◽  
Mihailo R. Jovanović

AbstractOver the last two decades, both experiments and simulations have demonstrated that transverse wall oscillations with properly selected amplitude and frequency can reduce turbulent drag by as much as $40\hspace{0.167em} \% $. In this paper, we develop a model-based approach for designing oscillations that suppress turbulence in a channel flow. We utilize eddy-viscosity-enhanced linearization of the turbulent flow with control in conjunction with turbulence modelling to determine skin-friction drag in a simulation-free manner. The Boussinesq eddy viscosity hypothesis is used to quantify the effect of fluctuations on the mean velocity in flow subject to control. In contrast to the traditional approach that relies on numerical simulations, we determine the turbulent viscosity from the second-order statistics of the linearized model driven by white-in-time stochastic forcing. The spatial power spectrum of the forcing is selected to ensure that the linearized model for uncontrolled flow reproduces the turbulent energy spectrum. The resulting correction to the turbulent mean velocity induced by small-amplitude wall movements is then used to identify the optimal frequency of drag-reducing oscillations. In addition, the control net efficiency and the turbulent flow structures that we obtain agree well with the results of numerical simulations and experiments. This demonstrates the predictive power of our model-based approach to controlling turbulent flows and is expected to pave the way for successful flow control at higher Reynolds numbers than currently possible.

Author(s):  
Elie Bou-Zeid ◽  
Charles Meneveau ◽  
Marc B. Parlange

We study turbulent flow over surfaces with varying roughness scales, using large eddy simulation (LES). The goal is to use LES results to formulate effective boundary conditions in terms of effective roughness height and blending height, to be used for RANS. The LES are implemented with the dynamic Smagorinsky model based on the Germano identity. However, as is well-known, when this identity is applied locally, it yields a coefficient with unphysically strong fluctuations and averaging is needed for better realism and numerical stability. The traditional approach consists of averaging over homogeneous directions, for example horizontal planes in channel flow. This requirement for homogeneous directions in the flow field and the concomitant inability to handle complex geometries renders the use of this model questionable in studying the effect of surface heterogeneity. Instead, a new version of the Lagrangian dynamic subgrid-scale (SGS) model [1] is implemented. A systematic set of simulations of flow over patches of differing roughness is performed, covering a wide range of patch length scales and surface roughness values. The simulated mean velocity profiles are analyzed to identify the height of the blending layer and used to measure the effective roughness length. Extending ideas introduced by Miyake [2] and Claussen [3], we have proposed a simple expression for effective surface roughness and blending height knowing local surface patch roughness values and their lengths [4]. Results of the model agreed well with the LES results when the heterogeneous surface consisted of patches of equal sizes. The model is tested here for surfaces with patches of different sizes.


1984 ◽  
Vol 106 (1) ◽  
pp. 134-139 ◽  
Author(s):  
H. G. Elrod ◽  
I. Anwar ◽  
R. Colsher

This paper presents some new equations for the treatment of turbulent lubricating films when the effects of inertia cannot be neglected. The eddy-viscosity concept is used to represent the turbulent stresses in terms of mean-velocity gradient. Transient and steady-state operation are both considered by means of a generalized stream-function-pressure development.


1974 ◽  
Vol 63 (4) ◽  
pp. 673-693 ◽  
Author(s):  
Russ E. Davis

The perturbation of a turbulent flow by an organized wavelike disturbance is examined using a dynamical, rather than phenomenological, approach. On the basis of the assumption that an infinitesimal perturbation results in a linear change in the statistics of the turbulence, and that the turbulence is either weak or that the turbulent perturbations are quasi-Gaussian, a method for predicting the perturbation turbulent Reynolds stresses is developed. The novel aspect of the analysis is that all averaging is delayed until the dynamical equations have been solved rather than attempting to find, apriori, equations for averaged quantities. When applied to long-wave perturbations the analysis indicates that the perturbation shear stress is of primary dynamical importance, and that this stress is determined by the principal component of mean shear through a relation which depends on the spectrum of the turbulent velocity component parallel to the gradient of the undisturbed mean velocity (the component perpendicular to the wall in a turbulent boundary layer). Theoretical arguments and observations are used to estimate the form of this spectrum in a constant stress shear layer. This results in a prediction of the constitutive law relating turbulent stress and the mean flow. The law is visco-elastic in nature, and is in agreement with the known constitutive relation for stress perturbations to a constant stress boundary layer; it resembles the eddy viscosity relation used successfully by others in describing perturbations in turbulent flows. The details of the constitutive law depend on how well the turbulence obeys Taylor's hypothesis that phase velocity equals mean flow velocity, and some insight into this question is given.


2016 ◽  
Vol 790 ◽  
pp. 104-127 ◽  
Author(s):  
S. Russo ◽  
P. Luchini

We identify a benchmark problem simple enough that it can be solved both by an eddy-viscosity model and by direct numerical simulation: this is the linear response of a turbulent flow’s mean-velocity profile to an external volume force. An example of such a force was found in a study of the perturbation induced by bottom topography by Luchini & Charru (J. Fluid Mech., vol. 656, 2010, pp. 337–341). On the other hand, a modification of the method by Quadrio & Luchini (Proceedings of the IX European Turbulence Conference, Southampton, UK, 2002, pp. 715–718) and Luchini et al. (Phys. Fluids, vol. 18, 2006, 121702) to compute the linear impulse response of a wall-bounded turbulent flow allows the response to a volume force to be computed directly. The comparison exhibits significant differences and suggests that there might be fundamental obstacles to designing an eddy-viscosity model that provides the correct result.


1982 ◽  
Vol 49 (4) ◽  
pp. 697-703 ◽  
Author(s):  
L. C. Thomas ◽  
H. M. Kadry

An analysis is presented for fully developed transitional turbulent flow between parallel plates that features the use of a turbulent burst model for the important wall region. Model closure is accomplished by the specification of the mean burst frequency and, for moderate to high Reynolds numbers, by matching with a classical eddy viscosity model for the turbulent core. Predictions obtained for friction factor and mean velocity distribution are compared with experimental data for fully developed transitional turbulent flow in a channel with large aspect ratio. Predictions are also developed for the eddy viscosity within the wall region for transitional turbulent conditions.


2020 ◽  
Vol 640 ◽  
pp. A53
Author(s):  
L. Löhnert ◽  
S. Krätschmer ◽  
A. G. Peeters

Here, we address the turbulent dynamics of the gravitational instability in accretion disks, retaining both radiative cooling and irradiation. Due to radiative cooling, the disk is unstable for all values of the Toomre parameter, and an accurate estimate of the maximum growth rate is derived analytically. A detailed study of the turbulent spectra shows a rapid decay with an azimuthal wave number stronger than ky−3, whereas the spectrum is more broad in the radial direction and shows a scaling in the range kx−3 to kx−2. The radial component of the radial velocity profile consists of a superposition of shocks of different heights, and is similar to that found in Burgers’ turbulence. Assuming saturation occurs through nonlinear wave steepening leading to shock formation, we developed a mixing-length model in which the typical length scale is related to the average radial distance between shocks. Furthermore, since the numerical simulations show that linear drive is necessary in order to sustain turbulence, we used the growth rate of the most unstable mode to estimate the typical timescale. The mixing-length model that was obtained agrees well with numerical simulations. The model gives an analytic expression for the turbulent viscosity as a function of the Toomre parameter and cooling time. It predicts that relevant values of α = 10−3 can be obtained in disks that have a Toomre parameter as high as Q ≈ 10.


Author(s):  
J. Szydlowski ◽  
M. Costes

This paper presents numerical simulations of the flow around a NACA 0015 airfoil at static and dynamic stall. The treatment of these configurations is a very challenging task for CFD applications. The turbulent flow around the static and in pitch oscillation airfoil is computed using different approaches: 2D RANS, 3D RANS and DES methodologies and with finer and finer meshes in order to try to reach a space converged solution. The main conclusion of the paper is that the prediction of static and all the more dynamic stall is not mature with present modeling capabilities.


1973 ◽  
Vol 95 (3) ◽  
pp. 415-421 ◽  
Author(s):  
A. J. Wheeler ◽  
J. P. Johnston

Predictions have been made for a variety of experimental three-dimensional boundary layer flows with a single finite difference method which was used with three different turbulent stress models: (i) an eddy viscosity model, (ii) the “Nash” model, and (iii) the “Bradshaw” model. For many purposes, even the simplest stress model (eddy viscosity) was adequate to predict the mean velocity field. On the other hand, the profile of shear stress direction was not correctly predicted in one case by any model tested. The high sensitivity of the predicted results to free stream pressure gradient in separating flow cases is demonstrated.


Author(s):  
N Kharoua ◽  
L Khezzar

Large eddy simulation of turbulent flow around smooth and rough hemispherical domes was conducted. The roughness of the rough dome was generated by a special approach using quadrilateral solid blocks placed alternately on the dome surface. It was shown that this approach is capable of generating the roughness effect with a relative success. The subgrid-scale model based on the transport of the subgrid turbulent kinetic energy was used to account for the small scales effect not resolved by large eddy simulation. The turbulent flow was simulated at a subcritical Reynolds number based on the approach free stream velocity, air properties, and dome diameter of 1.4 × 105. Profiles of mean pressure coefficient, mean velocity, and its root mean square were predicted with good accuracy. The comparison between the two domes showed different flow behavior around them. A flattened horseshoe vortex was observed to develop around the rough dome at larger distance compared with the smooth dome. The separation phenomenon occurs before the apex of the rough dome while for the smooth dome it is shifted forward. The turbulence-affected region in the wake was larger for the rough dome.


Sign in / Sign up

Export Citation Format

Share Document