scholarly journals A new linear forcing method for isotropic turbulence with controlled integral length scale

2021 ◽  
Vol 33 (4) ◽  
pp. 045127
Author(s):  
Jérémie Janin ◽  
Fabien Duval ◽  
Christophe Friess ◽  
Pierre Sagaut
2015 ◽  
Vol 772 ◽  
pp. 361-385 ◽  
Author(s):  
Anikesh Pal ◽  
Sutanu Sarkar

Direct numerical simulations are performed to study the evolution of a towed stratified wake subject to external turbulence in the background. A field of isotropic turbulence is combined with an initial turbulent wake field and the combined wake is simulated in a temporally evolving framework similar to that of Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482). Simulations are performed for external turbulence whose initial level varies between zero and a moderate intensity of up to 7 % relative to the free stream and whose initial integral length scale is of the same order as that of the wake turbulence. A series of simulations are carried out at a Reynolds number of 10 000 and Froude number of 3. Background turbulence, especially at a level of 3 % or above, is found to have substantial quantitative effects in the stratified simulations. Turbulence inside the wake increases due to the entrainment of external turbulence, and the energy transfer through turbulent production from mean to fluctuating velocity also increases, leading to reduced mean velocity. The profiles of normalized mean and turbulence quantities in the stratified wake exhibit little change in the vertical direction but the horizontal spread increases in comparison to the case with undisturbed background. The spatial organization of the internal wave field is disrupted even at the 1 % level of external turbulence. However, key characteristics of stratified wakes such as the formation of coherent pancake vortices and the long lifetime of the mean wake are robust to the presence of fluctuations in the background. A corresponding series of simulations for the unstratified situation is carried out at the same Reynolds number of 10 000 and with similar levels of external turbulence. The change of mean and turbulence statistics is found to be weaker in the unstratified cases compared with the corresponding stratified cases and also weaker relative to that found by Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482) at a similar level of external turbulence relative to the free stream and similar integral length scale. Theoretical arguments and additional simulations are provided to show that the level of external turbulence relative to wake turbulence (dissimilar between the present investigation and Rind & Castro (J. Fluid Mech., vol. 710, 2012a, p. 482)) is a key governing parameter in both stratified and unstratified backgrounds.


2015 ◽  
Vol 766 ◽  
pp. 499-527 ◽  
Author(s):  
Ugo Piomelli ◽  
Amirreza Rouhi ◽  
Bernard J. Geurts

AbstractWe propose a new length scale as a basis for the modelling of subfilter motions in large-eddy simulations (LES) of turbulent flow. Rather than associating the model length scale with the computational grid, we put forward an approximation of the integral length scale to achieve a non-uniform flow coarsening through spatial filtering that reflects the local, instantaneous turbulence activity. Through the introduction of this grid-independent, solution-specific length scale it becomes possible to separate the problem of representing small-scale turbulent motions in a coarsened flow model from that of achieving an accurate numerical resolution of the primary flow scales. The formulation supports the notion of grid-independent LES, in which a prespecified reliability measure is used. We investigate a length-scale definition based on the resolved turbulent kinetic energy (TKE) and its dissipation. The proposed approach, which we call integral length-scale approximation (ILSA) model, is illustrated for turbulent channel flow at high Reynolds numbers and for homogeneous isotropic turbulence (HIT). We employ computational optimization of the model parameter based on various measures of subfilter activity, using the successive inverse polynomial interpolation (SIPI) and establish the efficiency of this route to subfilter modelling.


2017 ◽  
Vol 818 ◽  
pp. 184-204 ◽  
Author(s):  
Alejandro J. Puga ◽  
John C. LaRue

Time-resolved velocity measurements are obtained using a hot-wire in a nearly homogeneous and isotropic flow downstream of an active grid for a range of Taylor Reynolds numbers from$191$to$659$. It is found that the dimensionless dissipation rate,$C_{\unicode[STIX]{x1D716}}$, is nearly a constant for sufficiently high values of Taylor Reynolds number,$R_{\unicode[STIX]{x1D706},u_{q}}$, and is approximately equal to$0.87$. This value is approximately$5\,\%$less than the value reported by Boset al.(Phys. Fluids, vol. 19 (4), 2007, 045101), which is obtained using DNS/LES (direct numerical simulation combined with large eddy simulation) for decaying homogeneous and isotropic turbulence, and is in excellent agreement with the active grid experiment of Thormann & Meneveau (Phys. Fluids, vol. 26 (2), 2014, 025112.). The results presented herein show that deviation from isotropy may cause inconsistencies in the computation of$C_{\unicode[STIX]{x1D716}}$. As a result, it is suggested that the velocity scale be the square root of the turbulence kinetic energy. The integral length scale measurements obtained from the longitudinal velocity correlation are in close agreement with the integral length scale measured from the peak of the energy spectrum,$\unicode[STIX]{x1D705}E_{11}(\unicode[STIX]{x1D705})$, where$\unicode[STIX]{x1D705}$is the wavenumber and$E_{11}(\unicode[STIX]{x1D705})$is the one-dimensional power spectrum of the downstream velocity.


Author(s):  
Mark W. Johnson

A numerical procedure for predicting the receptivity of laminar boundary layers to freestream turbulence consisting of vortex arrays with arbitrary orientation has been developed. Results show that the boundary layer is most receptivity to those vortices which have their axes approximately in the streamwise direction and vortex wavelengths of approximately 1.2 δ. The computed near wall gains for isotropic turbulence are similar in magnitude to previously published experimental values used to predict transition. The new procedure is therefore capable of predicting the development of the fluctuations in the laminar boundary layer from values of the freestream turbulence intensity and length scale and hence determining the start of transition without resorting to any empirical correlation.


Author(s):  
Yiannis Andreopoulos ◽  
Amir H. Danesh-Yazdi ◽  
Oleg Goushcha ◽  
Niell Elvin

Turbulent flows carry mechanical energy distributed over a range of temporal and spatial scales and their interaction with a thin immersed piezoelectric beam results in a strain field which generates electrical charge. This energy harvesting method can be used for developing self-powered electronic devices such as flow sensors. In the present experimental work, various energy harvesters were placed in a turbulent boundary layer or inside a decaying flow field of homogeneous and isotropic turbulence. The role of large instantaneous turbulent structures in this rather complex fluid-structure interaction is discussed in interpreting the electrical output results. The forces acting on the vibrating beams have been measured dynamically and a theory has been developed which incorporates the effects of mean local velocity, turbulence intensity, the relative size of the beam’s length to the integral length scale of turbulence, the structural properties of the beam and the electrical properties of the active piezoelectric layer to provide reasonable estimates of the mean electrical power output. Experiments have been carried out in which these fluidic harvesters are immersed first in inhomogeneous turbulence like that encountered in boundary layers developing over solid walls and homogeneous and isotopic turbulence for which a simplified analytical description exists. It was found that there is a non-linear effect of turbulence length scales on the power output of the fluidic harvesters.


Author(s):  
Marco Vanni

The stresses acting on aggregates smaller than the Kolmogorov length scale in homogeneous isotropic turbulence were estimated by a two-scale numerical simulation. The fluid dynamics at the scales larger than the Kolmogorov length scale was calculated by a Direct Numerical Simulation of the turbulent flow, in which the aggregates were modeled as point particles. Then, we adopted Stokesian Dynamics to evaluate the phenomena governed by the smooth velocity field of the smallest scales. At this level the disordered structure of the aggregates was modeled in detail, in order to take into account the role that the primary particles have in generating and transferring the internal stress. From this result, it was possible to evaluate the internal forces acting at intermonomer contacts and determine the occurrence of breakup as a consequence of the failure of intermonomer bonds. The method was applied to disordered aggregates with isostatic and highly hyperstatic structures, respectively.


2019 ◽  
Vol 871 ◽  
pp. 271-304 ◽  
Author(s):  
Adrian C. H. Lai ◽  
Adrian Wing-Keung Law ◽  
E. Eric Adams

Buoyant jets or forced plumes are discharged into a turbulent ambient in many natural and engineering applications. The background turbulence generally affects the mixing characteristics of the buoyant jet, and the extent of the influence depends on the characteristics of both the jet discharge and ambient. Previous studies focused on the experimental investigation of the problem (for pure jets or plumes), but the findings were difficult to generalize because suitable scales for normalization of results were not known. A model to predict the buoyant jet mixing in the presence of background turbulence, which is essential in many applications, is also hitherto not available even for a background of homogeneous and isotropic turbulence (HIT). We carried out experimental and theoretical investigations of a buoyant jet discharging into background HIT. Buoyant jets were designed to be in the range of $1<z/l_{M}<5$, where $l_{M}=M_{o}^{3/4}/F_{o}^{1/2}$ is the momentum length scale, with $z/l_{M}<\sim 1$ and $z/l_{M}>\sim 6$ representing the asymptotic cases of pure jets and plumes, respectively. The background turbulence was generated using a random synthetic jet array, which produced a region of approximately isotropic and homogeneous field of turbulence to be used in the experiments. The velocity scale of the jet was initially much higher, and the length scale smaller, than that of the background turbulence, which is typical in most applications. Comprehensive measurements of the buoyant jet mixing characteristics were performed up to the distance where jet breakup occurred. Based on the experimental findings, a critical length scale $l_{c}$ was identified to be an appropriate normalizing scale. The momentum flux of the buoyant jet in background HIT was found to be conserved only if the second-order turbulence statistics of the jet were accounted for. A general integral jet model including the background HIT was then proposed based on the conservation of mass (using the entrainment assumption), total momentum and buoyancy fluxes, and the decay function of the jet mean momentum downstream. Predictions of jet mixing characteristics from the new model were compared with experimental observation, and found to be generally in agreement with each other.


2010 ◽  
Vol 650 ◽  
pp. 5-55 ◽  
Author(s):  
FRANCESCO LUCCI ◽  
ANTONINO FERRANTE ◽  
SAID ELGHOBASHI

This study investigates the two-way coupling effects of finite-size solid spherical particles on decaying isotropic turbulence using direct numerical simulation with an immersed boundary method. We fully resolve all the relevant scales of turbulence around freely moving particles of the Taylor length-scale size, 1.2≤d/λ≤2.6. The particle diameter and Stokes number in terms of Kolmogorov length- and time scales are 16≤d/η≤35 and 38≤τp/τk≤178, respectively, at the time the particles are released in the flow. The particles mass fraction range is 0.026≤φm≤1.0, corresponding to a volume fraction of 0.01≤φv≤0.1 and density ratio of 2.56≤ρp/ρf≤10. The maximum number of dispersed particles is 6400 for φv=0.1. The typical particle Reynolds number is of O(10). The effects of the particles on the temporal development of turbulence kinetic energy E(t), its dissipation rate (t), its two-way coupling rate of change Ψp(t) and frequency spectra E(ω) are discussed.In contrast to particles with d < η, the effect of the particles in this study, with d > η, is that E(t) is always smaller than that of the single-phase flow. In addition, Ψp(t) is always positive for particles with d > η, whereas it can be positive or negative for particles with d < η.


Sign in / Sign up

Export Citation Format

Share Document