Hierarchy of structure functions for passive scalars advected by turbulent flows

A hierarchical model is proposed for the joint moments of the passive scalar dissipation and the velocity dissipation in fluid turbulence. This model predicts that the joint probability density function (PDF) of the dissipations is a bivariate log-Poisson. An analytical calculation of the scaling exponents of structure functions of the passive scalar is carried out for this hierarchical model, showing a good agreement with the results of direct numerical simulations and experiments.

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Simulation of Vortex-dominant and Turbulent Flows in Natural Environment

In the present talk, the simulation of vortex dominant and turbulent flows are primarily addressed. To cope with complicated circumstances in environmental flows we illustrate the strategy of combining simplified physical model and suitable algorithm by a few examples.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Experimental investigation on TBAB clathrate hydrate slurry flows in a horizontal tube: Forced convective heat transfer behaviors

Tetra-n-butyl-ammonium bromide (TBAB) clathrate hydrate slurry (CHS) is one kind of secondary refrigerants, which is promising to be applied into air-conditioning or latent-heat transportation systems as a thermal storage or cold carrying medium for energy saving. It is a solid-liquid two phase mixture which is easy to produce and has high latent heat and good fluidity. In this paper, the heat transfer characteristics of TBAB slurry were investigated in a horizontal stainless steel tube under different solid mass fractions and flow velocities with constant heat flux. One velocity region of weakened heat transfer was found. Moreover, TBAB CHS was treated as a kind of Bingham fluids, and the influences of the solid particles, flow velocity and types of flow on the forced convective heat transfer coefficients of TBAB CHS were investigated. At last, criterial correlations of Nusselt number for laminar and turbulent flows in the form of power function were summarized, and the error with experimental results was within 20%.

Small-scale turbulent fluctuations beyond Taylor's frozen-flow hypothesis

The space-time cross-correlation function C-T(r, tau) of local temperature fluctuations in turbulent Rayleigh-Benard convection is obtained from simultaneous two-point time series measurements. The obtained C-T(r, tau) is found to have the scaling form C-T(r(E), 0) with r(E)=[(r-U tau)(2)+ V-2 tau(2)](1/2), where U and V are two characteristic velocities associated with the mean and rms velocities of the flow. The experiment verifies the theory and demonstrates its applications to a class of turbulent flows in which the requirement of Taylor's frozen flow hypothesis is not met.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.

Variation of drag coefficient with wind wave development

Turbulent air flows over developing wind waves in the air-sea boundary layer are numerically simulated without considering wave breaking. Influences of wind waves on air flows are considered using a model of significant wave and surface roughness, with a formula proposed for calculating the surface roughness, k - epsilon model is adopted to simulate turbulent flows. The results of the drag coefficient and turbulence characteristics agree well with the observations.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Short-and resonant-range interactions between scales in turbulence and their applications

Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence, were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.

A closure theory of intermittency of turbulence

The variational approach to the closure problem of turbulence theory, proposed in an earlier article [Phys. Fluids 26, 2098 (1983); 27, 2229 (1984)], is extended to evaluate the flatness factor, which indicates the degree of intermittency of turbulence. Since the flatness factor is related to the fourth moment of a turbulent velocity field, the corresponding higher-order terms in the perturbation solution of the Liouville equation have to be considered. Most closure methods discard these higher-order terms and fail to explain the intermittency phenomenon. The computed flatness factor of the idealized model of infinite isotropic turbulence ranges from 9 to 15 and has the same order of magnitude as the experimental data of real turbulent flows. The intermittency phenomenon does not necessarily negate the Kolmogorov k−5/3 inertial range spectrum. The Kolmogorov k−5/3 law and the high degree of intermittency can coexist as two consistent consequences of the closure theory of turbulence. The Kolmogorov 1941 theory [J. Fluid Mech. 62, 305 (1974)] cannot be disqualified merely because the energy dissipation rate fluctuates.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.