Transport equations for the normalized nth-order moments of velocity derivatives in grid turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
S.L. Tang ◽  
R.A. Antonia ◽  
L. Djenidi

Transport equations for the normalized moments of the longitudinal velocity derivative ${F_{n + 1}}$ (here, $n$ is $1, 2, 3\ldots$ ) are derived from the Navier–Stokes (N–S) equations for shearless grid turbulence. The effect of the (large-scale) streamwise advection of ${F_{n + 1}}$ by the mean velocity on the normalized moments of the velocity derivatives can be expressed as $C_1 {F_{n + 1}}/Re_\lambda$ , where $C_1$ is a constant and $Re_\lambda$ is the Taylor microscale Reynolds number. Transport equations for the normalized odd moments of the transverse velocity derivatives ${F_{y,n + 1}}$ (here, $n$ is 2, 4, 6), which should be zero if local isotropy is satisfied, are also derived and discussed in sheared and shearless grid turbulence. The effect of the (large-scale) streamwise advection term on the normalized moments of the velocity derivatives can also be expressed in the form $C_2 {F_{y,n + 1}}/Re_\lambda$ , where $C_2$ is a constant. Finally, the contribution of the mean shear in the transport equation for ${F_{n + 1}}$ can be modelled as $15 B/Re_\lambda$ , where $B$ ( $=S^*{S_{s,n + 1}}$ ) is the product of the non-dimensional shear parameter $S^*$ and the normalized mixed longitudinal-transverse velocity derivatives ${{S_{s,n + 1}}}$ ; if local isotropy is satisfied, $S_{s,n + 1}$ should be zero. These results indicate that if ${F_{n + 1}}$ , ${F_{y,n + 1}}$ and $B$ do not increase as rapidly as $Re_\lambda$ , then the effect of the large-scale structures on small-scale turbulence will disappear when $Re_\lambda$ becomes sufficiently large.

1998 ◽  
Vol 374 ◽  
pp. 29-57 ◽  
Author(s):  
R. A. ANTONIA ◽  
T. ZHOU ◽  
Y. ZHU

All components of the fluctuating vorticity vector have been measured in decaying grid turbulence using a vorticity probe of relatively simple geometry (four X-probes, i.e. a total of eight hot wires). The data indicate that local isotropy is more closely satisfied than global isotropy, the r.m.s. vorticities being more nearly equal than the r.m.s. velocities. Two checks indicate that the performance of the probe is satisfactory. Firstly, the fully measured mean energy dissipation rate 〈ε〉 is in good agreement with the value inferred from the rate of decay of the mean turbulent energy 〈q2〉 in the quasi-homogeneous region; the isotropic mean energy dissipation rate 〈εiso〉 agrees closely with this value even though individual elements of 〈ε〉 indicate departures from isotropy. Secondly, the measured decay rate of the mean-square vorticity 〈ω2〉 is consistent with that of 〈q2〉 and in reasonable agreement with the isotropic form of the transport equation for 〈ω2〉. Although 〈ε〉≃〈εiso〉, there are discernible differences between the statistics of ε and εiso; in particular, εiso is poorly correlated with either ε or ω2. The behaviour of velocity increments has been examined over a narrow range of separations for which the third-order longitudinal velocity structure function is approximately linear. In this range, transverse velocity increments show larger departures than longitudinal increments from predictions of Kolmogorov (1941). The data indicate that this discrepancy is only partly associated with differences between statistics of locally averaged ε and ω2, the latter remaining more intermittent than the former across this range. It is more likely caused by a departure from isotropy due to the small value of Rλ, the Taylor microscale Reynolds number, in this experiment.


2014 ◽  
Vol 748 ◽  
Author(s):  
F. Thiesset ◽  
R. A. Antonia ◽  
L. Djenidi

AbstractOn the basis of a two-point similarity analysis, the well-known power-law variations for the mean kinetic energy dissipation rate $\overline{\epsilon }$ and the longitudinal velocity variance $\overline{u^2}$ on the axis of a round jet are derived. In particular, the prefactor for $\overline{\epsilon } \propto (x-x_0)^{-4}$, where $x_0$ is a virtual origin, follows immediately from the variation of the mean velocity, the constancy of the local turbulent intensity and the ratio between the axial and transverse velocity variance. Second, the limit at small separations of the two-point budget equation yields an exact relation illustrating the equilibrium between the skewness of the longitudinal velocity derivative $S$ and the destruction coefficient $G$ of enstrophy. By comparing the latter relation with that for homogeneous isotropic decaying turbulence, it is shown that the approach towards the asymptotic state at infinite Reynolds number of $S+2G/R_{\lambda }$ in the jet differs from that in purely decaying turbulence, although $S+2G/R_{\lambda } \propto R_{\lambda }^{-1}$ in each case. This suggests that, at finite Reynolds numbers, the transport equation for $\overline{\epsilon }$ imposes a fundamental constraint on the balance between $S$ and $G$ that depends on the type of large-scale forcing and may thus differ from flow to flow. This questions the conjecture that $S$ and $G$ follow a universal evolution with $R_{\lambda }$; instead, $S$ and $G$ must be tested separately in each flow. The implication for the constant $C_{\epsilon 2}$ in the $k-\overline{\epsilon }$ model is also discussed.


2010 ◽  
Vol 661 ◽  
pp. 341-364 ◽  
Author(s):  
D. CHUNG ◽  
B. J. McKEON

We investigate statistics of large-scale structures from large-eddy simulation (LES) of turbulent channel flow at friction Reynolds numbers Reτ = 2K and 200K (where K denotes 1000). In order to capture the behaviour of large-scale structures properly, the channel length is chosen to be 96 times the channel half-height. In agreement with experiments, these large-scale structures are found to give rise to an apparent amplitude modulation of the underlying small-scale fluctuations. This effect is explained in terms of the phase relationship between the large- and small-scale activity. The shape of the dominant large-scale structure is investigated by conditional averages based on the large-scale velocity, determined using a filter width equal to the channel half-height. The conditioned field demonstrates coherence on a scale of several times the filter width, and the small-scale–large-scale relative phase difference increases away from the wall, passing through π/2 in the overlap region of the mean velocity before approaching π further from the wall. We also found that, near the wall, the convection velocity of the large scales departs slightly, but unequivocally, from the mean velocity.


2009 ◽  
Vol 631 ◽  
pp. 199-230 ◽  
Author(s):  
J. WESTERWEEL ◽  
C. FUKUSHIMA ◽  
J. M. PEDERSEN ◽  
J. C. R. HUNT

Conditionally sampled measurements with particle image velocimetry (PIV) of a turbulent round submerged liquid jet in a laboratory have been taken at Re = 2 × 103 between 60 and 100 nozzle diameters from the nozzle in order to investigate the dynamics and transport processes at the continuous and well-defined bounding interface between the turbulent and non-turbulent regions of flow. The jet carries a fluorescent dye measured with planar laser-induced fluorescence (LIF), and the surface discontinuity in the scalar concentration is identified as the fluctuating turbulent jet interface. Thence the mean outward ‘boundary entrainment’ velocity is derived and shown to be a constant fraction (about 0.07) of the the mean jet velocity on the centreline. Profiles of the conditional mean velocity, mean scalar and momentum flux show that at the interface there are clear discontinuities in the mean axial velocity and mean scalar and a tendency towards a singularity in mean vorticity. These actual or asymptotic discontinuities are consistent with the conditional mean momentum and scalar transport equations integrated across the interface. Measurements of the fluxes of turbulent kinetic energy and enstrophy are consistent with computations by Mathew & Basu (Phys. Fluids, vol. 14, 2002, pp. 2065–2072) in showing that for a jet flow (without forcing) the entrainment process is dominated by small-scale eddying at the highly sheared interface (‘nibbling’), with large-scale engulfing making a small (less than 10%) contribution consistent with concentration measurements showing that the interior of the jet is well mixed. (Turbulent jets differ greatly from the free shear layer in this respect.) To explain the difference between velocity and scalar profiles, their conditional mean gradients are defined in terms of a local eddy viscosity and eddy diffusivity and the momentum and scalar fluxes inside the interface. Since the eddy diffusivity is larger than the eddy viscosity, the scalar profile is flatter inside the interface so that the scalar discontinuity is relatively greater than the mean velocity discontinuity. Theoretical arguments, following Hunt, Eames & Westerweel (in Proc. of the IUTAM Symp. on Computational Physics and New Perspectives in Turbulence, ed. Y. Kaneda, vol. 4, 2008, pp. 331–338, Springer), are proposed for how the vortex sheet develops, how the internal structure of the interface layer relates to the inhomogeneous rotational and irrotational motions on each side and why the dominant entrainment process of jets and wakes differs from that of free shear layers.


2021 ◽  
Vol 929 ◽  
Author(s):  
N. Agastya Balantrapu ◽  
Christopher Hickling ◽  
W. Nathan Alexander ◽  
William Devenport

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.


1949 ◽  
Vol 2 (4) ◽  
pp. 451 ◽  
Author(s):  
AA Townsend

Extending previous work on turbulent diffusion in the wake of a circular-cylinder, a series of measurements have been made of the turbulent transport of mean stream momentum, turbulent energy, and heat in the wake of a cylinder of 0.169 cm. diameter, placed in an air-stream of velocity 1280 cm. sec.-1. It has been possible to extend the measurements to 960 diameters down-stream from the cylinder, and it 1s found that, at distances in excess of 600 diameters, the requirements of dynamical similarity are very nearly satisfied. To account for the observed rates of transport of turbulent energy and heat, it is necessary that only part of this transport be due to bulk convection by the slow large-scale motion of the jets of turbulent fluid emitted by the central, fully turbulent core of the wake, which had been supposed previously to perform most of the transport. The remainder of the transport is carried out by the small-scale diffusive motion of the turbulent eddies within the jets, and may be described by assigning diffusion coefficients to the turbulent fluid. It is found that the diffusion coefficients for momentum and heat are approximately equal, but that for turbulent energy is considerably smaller. On the basis of these hypotheses, it is possible to calculate $he form of the mean velocity distribution in good agreement with experiment, and to give a qualitative explanation of the apparently more rapid diffusion of heat.


2018 ◽  
Vol 857 ◽  
pp. 907-936 ◽  
Author(s):  
A. Cimarelli ◽  
A. Leonforte ◽  
D. Angeli

The separating and reattaching flows and the wake of a finite rectangular plate are studied by means of direct numerical simulation data. The large amount of information provided by the numerical approach is exploited here to address the multi-scale features of the flow and to assess the self-sustaining mechanisms that form the basis of the main unsteadinesses of the flows. We first analyse the statistically dominant flow structures by means of three-dimensional spatial correlation functions. The developed flow is found to be statistically dominated by quasi-streamwise vortices and streamwise velocity streaks as a result of flow motions induced by hairpin-like structures. On the other hand, the reverse flow within the separated region is found to be characterized by spanwise vortices. We then study the spectral properties of the flow. Given the strongly inhomogeneous nature of the flow, the spectral analysis has been conducted along two selected streamtraces of the mean velocity field. This approach allows us to study the spectral evolution of the flow along its paths. Two well-separated characteristic scales are identified in the near-wall reverse flow and in the leading-edge shear layer. The first is recognized to represent trains of small-scale structures triggering the leading-edge shear layer, whereas the second is found to be related to a very large-scale phenomenon that embraces the entire flow field. A picture of the self-sustaining mechanisms of the flow is then derived. It is shown that very-large-scale fluctuations of the pressure field alternate between promoting and suppressing the reverse flow within the separation region. Driven by these large-scale dynamics, packages of small-scale motions trigger the leading-edge shear layers, which in turn created them, alternating in the top and bottom sides of the rectangular plate with a relatively long period of inversion, thus closing the self-sustaining cycle.


2017 ◽  
Vol 826 ◽  
pp. 888-917 ◽  
Author(s):  
Valentin Resseguier ◽  
Etienne Mémin ◽  
Dominique Heitz ◽  
Bertrand Chapron

We present here a new stochastic modelling approach in the constitution of fluid flow reduced-order models. This framework introduces a spatially inhomogeneous random field to represent the unresolved small-scale velocity component. Such a decomposition of the velocity in terms of a smooth large-scale velocity component and a rough, highly oscillating component gives rise, without any supplementary assumption, to a large-scale flow dynamics that includes a modified advection term together with an inhomogeneous diffusion term. Both of those terms, related respectively to turbophoresis and mixing effects, depend on the variance of the unresolved small-scale velocity component. They bring an explicit subgrid term to the reduced system which enables us to take into account the action of the truncated modes. Besides, a decomposition of the variance tensor in terms of diffusion modes provides a meaningful statistical representation of the stationary or non-stationary structuration of the small-scale velocity and of its action on the resolved modes. This supplies a useful tool for turbulent fluid flow data analysis. We apply this methodology to circular cylinder wake flow at Reynolds numbers $Re=100$ and $Re=3900$. The finite-dimensional models of the wake flows reveal the energy and the anisotropy distributions of the small-scale diffusion modes. These distributions identify critical regions where corrective advection effects, as well as structured energy dissipation effects, take place. In providing rigorously derived subgrid terms, the proposed approach yields accurate and robust temporal reconstruction of the low-dimensional models.


Author(s):  
Walter Boscheri ◽  
Giacomo Dimarco ◽  
Lorenzo Pareschi

In this paper, we propose a novel space-dependent multiscale model for the spread of infectious diseases in a two-dimensional spatial context on realistic geographical scenarios. The model couples a system of kinetic transport equations describing a population of commuters moving on a large scale (extra-urban) with a system of diffusion equations characterizing the non-commuting population acting over a small scale (urban). The modeling approach permits to avoid unrealistic effects of traditional diffusion models in epidemiology, like infinite propagation speed on large scales and mass migration dynamics. A construction based on the transport formalism of kinetic theory allows to give a clear model interpretation to the interactions between infected and susceptible in compartmental space-dependent models. In addition, in a suitable scaling limit, our approach permits to couple the two populations through a consistent diffusion model acting at the urban scale. A discretization of the system based on finite volumes on unstructured grids, combined with an asymptotic preserving method in time, shows that the model is able to describe correctly the main features of the spatial expansion of an epidemic. An application to the initial spread of COVID-19 is finally presented.


2016 ◽  
Vol 31 (4) ◽  
pp. 1397-1405
Author(s):  
Weihong Qian ◽  
Ning Jiang ◽  
Jun Du

Abstract Mathematical derivation, meteorological justification, and comparison to model direct precipitation forecasts are the three main concerns recently raised by Schultz and Spengler about moist divergence (MD) and moist vorticity (MV), which were introduced in earlier work by Qian et al. That previous work demonstrated that MD (MV) can in principle be derived mathematically with a value-added empirical modification. MD (MV) has a solid meteorological basis. It combines ascent motion and high moisture: the two elements necessary for rainfall. However, precipitation efficiency is not considered in MD (MV). Given the omission of an advection term in the mathematical derivation and the lack of precipitation efficiency, MD (MV) might be suitable mainly for heavy rain events with large areal coverage and long duration caused by large-scale quasi-stationary weather systems, but not for local intense heavy rain events caused by small-scale convection. In addition, MD (MV) is not capable of describing precipitation intensity. MD (MV) worked reasonably well in predicting heavy rain locations from short to medium ranges as compared with the ECMWF model precipitation forecasts. MD (MV) was generally worse than (though sometimes similar to) the model heavy rain forecast at shorter ranges (about a week) but became comparable or even better at longer ranges (around 10 days). It should be reiterated that MD (MV) is not intended to be a primary tool for predicting heavy rain areas, especially in the short range, but is a useful parameter for calibrating model heavy precipitation forecasts, as stated in the original paper.


Sign in / Sign up

Export Citation Format

Share Document