Dns Of Compressible Turbulent Boundary Layer Around A Sharp Cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20 degrees cone angle (or 10 degrees half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

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.

Statistical Interpolation Method of Turbulent Phase Screen

A relative displacement between the grid points of optical fields and those of phase screens may occur in the simulation of light propagation through the turbulent atmosphere. A statistical interpolator is proposed to solve this problem in this paper. It is evaluated by the phase structure function and numerical experiments of light propagation through atmospheric turbulence with/without adaptive optics (AO) and it is also compared with the well-known linear interpolator under the same condition. Results of the phase structure function show that the statistical interpolator is more accurate in comparison with the linear one, especially in the high frequency region. More importantly, the long-exposure results of light propagation through the turbulent atmosphere with/without AO also show that the statistical interpolator is more accurate and reliable than the linear one. (C) 2009 Optical Society of America.

On separated shear layer of blunt circular cylinder

Separated shear layer of blunt circular cylinder has been experimentally investigated for the Reynolds numbers (based on the diameter) ranging from 2.8 x 10(3) to 1.0 x 10(5), with emphasis on evolution of separated shear layer, its structure and distribution of Reynolds shear stress and turbulence kinetic energy. The results demonstrate that laminar separated shear layer experiences 2 similar to 3 times vortex merging before it reattaches, and turbulence separated shear layer takes 5 similar to 6 times vortex merging. In addition, relationship between dimensionless initial frequencies of K-H instability and Reynolds numbers is identified, and reasons for the decay of turbulence kinetic energy and Reynolds shear stress in reattachment region are discussed.

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Observation and modeling for terrestrial processes in alpine meadow

The water-heat transfer process between land and atmosphere in Haibei alpine meadow area has been systematically observed. A multi-layer coupling model for land-atmosphere interaction was presented with special attention paid to the moisture transfer in leaf stomata under unsaturated condition. A profound investigation on the physical process of turbulent transfer inside the vegetation has been performed with a revised formula of water absorption for root system. The present model facilitates the study of vertically distributed physical variables in detail. Numerical simulation was conducted according to the transfer process of Kinesia humility meadow in the area of Haibei Alpine Meadow Ecosystem Station, CAS. The calculated results agree well with observation.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.