新型大涡数值模拟亚格子模型及应用

基于湍流大小尺度间动量输运的结构函数方程，提出了一种新的湍流大涡模型(LES)亚格子涡粘模式．新亚格子涡粘系数正比于纵向速度增量的扭率，它表征大小尺度湍流间的能量输运和耗散之比．新模式通过各向同性湍流直接数值模拟数据库的检验，并用于槽道湍流的大涡模拟计算，将所得结果与DNS结果进行了比较．     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.     

Modelling of turbulent flows in transonic axial-flow compressor NASA Rotor 37 with local weak-equilibrium damping of eddy viscosity coefficient

A local damping of eddy viscosity depending on the ratio of the production of turbulent energy to turbulence dissipation rate is proposed at the computation of flows in transonic axial compressors. The results of the numerical modeling of flows in compressor NASA Rotor 37 are presented, and the computed distributions of the increase in the total temperature and total pressure are compared with experimental data for different rotation frequencies. An increase in the accuracy of modelling was obtained in all considered regimes.     

A one-equation turbulence model for recirculating flows

A one-equation turbulence model which relies on the turbulent kinetic energy transport equation has been developed to predict the flow properties of the recirculating flows. The turbulent eddy-viscosity coefficient is computed from a recalibrated Bradshaw’s assumption that the constant a₁ = 0.31 is recalibrated to a function based on a set of direct numerical simulation (DNS) data. The values of dissipation of turbulent kinetic energy consist of the near-wall part and isotropic part, and the isotropic part involves the von Karman length scale as the turbulent length scale. The performance of the new model is evaluated by the results from DNS for fully developed turbulence channel flow with a wide range of Reynolds numbers. However, the computed result of the recirculating flow at the separated bubble of NACA4412 demonstrates that an increase is needed on the turbulent dissipation, and this leads to an advanced tuning on the self-adjusted function. The improved model predicts better results in both the non-equilibrium and equilibrium flows, e.g. channel flows, backward-facing step flow and hump in a channel.     

等离子体射流产生与特性的实验研究

本文报道热等离子体射流产生及射流特性的实验研究结果。采用同一个直流等离子体发生器,工作气体流量小时产生出层流等离子体长射流,射流长度随气体流量或弧电流的增加而明显增加;工作气体流量大时则产生出湍流等离子体短射流,此时射流长度几乎与工作气体流量或弧电流无关;在层流与湍流等离子体射流工况之间,存在一个流动状况不稳定的过渡区,此时等离子体射流的平均长度随气流量的加大而减小,但随弧电流的加大而明显加大。层流等离子体长射流有相当好的刚性。     

Temporal evolution of flame stretch due to turbulence and the hydrodynamic instability

The temporal evolution of the strain rate on a turbulent premixed flame was measured experimentally using cinema-stereoscopic particle image velocimetry. Turbulence strains a flame due to velocity gradients associated both directly with the turbulence and those caused by the hydrodynamic instability, which are initiated by the turbulence. The development of flame wrinkles caused by both of these mechanisms was observed. Wrinkles generated by the turbulence formed around vortical structures, which passed through the flame and were attenuated. After the turbulent structures had passed, the hydrodynamic instability flow pattern developed and caused additional strain. The hydrodynamic instability also caused the growth of small flame front perturbations into large wrinkles. In the moderately turbulent flame investigated, it was found that the evolution of the strain rate caused by turbulence–flame interactions followed a common pattern involving three temporal regimes. In the first, the turbulence exerted extensive (positive) strain on the flame, creating a wrinkle that had negative curvature (concave towards the reactants). This was followed by a transition period, leading into the third regime in which the flow pattern and strain rate were dominated by the hydrodynamic instability mechanism. It was also found that the magnitudes of the strain rate in the first and third regimes were similar. Hence, the hydrodynamic instability mechanism caused significant strain on a flame and should be included in turbulent combustion models.     

Cross-scale effects in solar-wind turbulence

The understanding of the small-scale termination of the turbulent energy cascade in collisionless plasmas is nowadays one of the outstanding problems in space physics. In the absence of collisional viscosity, the dynamics at small scales is presumably kinetic in nature; the identification of the physical mechanism which replaces energy dissipation and establishes the link between macroscopic and microscopic scales would open a new scenario in the study of turbulent heating in space plasmas. We present a numerical analysis of kinetic effects along the turbulent energy cascade in solar-wind plasmas which provides an effective unified interpretation of a wide set of spacecraft observations and shows that, simultaneously with an increase in the ion perpendicular temperature, strong bursts of electrostatic activity in the form of ion-acoustic turbulence are produced together with accelerated beams in the ion distribution function.     

Classification of turbulence modification by dispersed spheres using a novel dimensionless number

A novel dimensionless parameter, the particle moment number Pa, was derived using dimensional analysis of the particle-laden Navier-Stokes equations, in order to understand the underlying physics of turbulence modification by particles. A set of 80 previous experimental measurements where the turbulent kinetic energy was modified by particles was examined, and all results could clearly be divided into three groups in Re-Pa mappings. The turbulence attenuation region was observed between the augmentation regions with two critical particle momentum numbers.     

The effect of turbulence damping on acoustic wave propagation in tubes

The attenuation of sound due to the interaction between a low Mach number turbulent boundary layer and acoustic waves can be significant at low frequencies or in narrow tubes. In a recent publication by the present authors the acoustics of charge air coolers for passenger cars has been identified as an interesting application where turbulence attenuation can be of importance. Favourable low-frequency damping has been observed that could be used for control of the in-duct sound that is created by the engine gas exchange process. Analytical frequency-dependent models for the eddy viscosity that controls the momentum and thermal boundary layers are available but are restricted to thin acoustic boundary layers. For cases with cross-sections of a few millimetres a model based on thin acoustic boundary layers will not be applicable in the frequency range of interest.In the present paper a frequency-dependent axis-symmetric numerical model for interaction between turbulence and acoustic waves is proposed. A finite element scheme is used to formulate the time harmonic linearized convective equations for conservation of mass, momentum and energy into one coupled system of equations. The turbulence is introduced with a linear model for the eddy viscosity that is added to the shear viscosity. The proposed model is validated by comparison with experimental data from the literature.     

On mean flow universality of turbulent wall flows. II. Asymptotic flow analysis

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence statistics, their universality, the convergence of statistics towards their asymptotic profiles, and the overall physical flow organization. In extension of recent studies focusing on the mean flow at moderate and relatively high Re, the latter questions are addressed with respect to three canonical wall-bounded flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Main results reported here are the asymptotic logarithmic law for the mean velocity and corresponding scale-separation laws for bulk flow properties, the Reynolds shear stress, the turbulence production and turbulent viscosity. A scaling analysis indicates that the establishment of a self-similar turbulence state is the condition for the development of a strict logarithmic velocity profile. The resulting overall physical flow structure at extreme Re is discussed.     

A theoretical model of turbulent fiber suspension and its application to the channel flow

A theoretical model of turbulent fiber suspension is developed by deriving the equations of Reynolds averaged Navier-Stokes,turbulence kinetic energy and turbulence dissipation rate with the additional term of fibers.In order to close the above equations,the equation of probability distribution function for mean fiber orientation is also derived.The theoretical model is applied to the turbulent channel flow and the corresponding equations are solved numerically.The numerical results are verified by comparisons with the experimental ones.The effects of Reynolds number,fiber concentration and fiber aspect-ratio on the velocity profile,turbulent kinetic energy and turbulent dissipation rate are analyzed.Based on the numerical data,the expression for the velocity profile in the turbulent fiber suspension channel flow,which includes the effect of Reynolds number,fiber concentration and aspect-ratio,is proposed.     

The foundations of turbulent motion flow between parallel planes

It is suggested that turbulence is a kind of macroscopic viscosity; the usual turbulent motion being a combination of purely dynamical motion and (proper) turbulence. A flow in a channel is considered, the equation for the spectrum of (proper) turbulence and the expression for the spectrum are derived.     

Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows

The localized artificial diffusivity method is investigated in the context of large-eddy simulation of compressible turbulent flows. The performance of different artificial bulk viscosity models are evaluated through detailed results from the evolution of decaying compressible isotropic turbulence with eddy shocklets and supersonic turbulent boundary layer. Effects of subgrid-scale (SGS) models and implicit time-integration scheme/time-step size are also investigated within the framework of the numerical scheme used. The use of a shock sensor along with artificial bulk viscosity significantly improves the scheme for simulating turbulent flows involving shocks while retaining the shock-capturing capability. The proposed combination of Ducros-type sensor with a negative dilatation sensor removes unnecessary bulk viscosity within expansion and weakly compressible turbulence regions without shocks and allows it to localize near the shocks. It also eliminates the need for a wall-damping function for the bulk viscosity while simulating wall-bounded turbulent flows. For the numerical schemes used, better results are obtained without adding an explicit SGS model than with SGS model at moderate Reynolds number. Inclusion of a SGS model in addition to the low-pass filtering and artificial bulk viscosity results in additional damping of the resolved turbulence. However, investigations at higher Reynolds numbers suggest the need for an explicit SGS model. The flow statistics obtained using the second-order implicit time-integration scheme with three sub-iterations closely agrees with the explicit scheme if the maximum Courant–Friedrichs–Lewy is kept near unity.     

Contribution to the nonlinear theory of sound and hydrodynamic turbulence of a compressible liquid

The interaction of sound with hydrodynamic turbulence has been studied in detail. The sound absorption decrement, the correlation time and length and the frequency diffusion coefficient for the acoustic wave packet are calculated. The spectral composition of the sound radiated by a unit, turbulent volume and the spectral energy density of sound in equilibrium with the turbulence are studied. The region of applicability of the kinetic equation for sound with a linear dispersion low is found.     

Validity of Taylor's dissipation-viscosity independence postulate in variable-viscosity turbulent fluid mixtures

G. I. Taylor's postulate [Proc. R. Soc. A 151, 421 (1935)] that dissipation is independent of viscosity at high Reynolds numbers is the foundation of many single-fluid turbulence theories and closure models. The validity of this key postulate in an important class of flows, turbulent mixtures, is not yet clearly established. We devise a simple numerical experiment of decaying turbulence in a mixture of two fluids of vastly different viscosities to examine dissipation scaling. Initially, the two fluids are segregated, and dissipation is directly proportional to viscosity. As turbulence evolves and fluids mix, the velocity gradients rapidly adapt to the viscosity field, and within one-half eddy turnover time, dissipation-viscosity independence is established. Viscosity-weighted velocity-gradient skewness is shown to be constant, leading to the validity of Taylor's postulate in turbulent mixtures.     

Turbulent vortex rings in weakly turbulent and in turbulent ambient flow

Turbulent vortex rings were investigated in weakly turbulent flow and in three different grid generated turbulent flows to clarify the reciprocal action of the vortex ring with defined external turbulence. Assuming self-similarity and turbulent viscosity as proportional to V₀D₀ the equations for the ring diameter D(t) and the velocity of propagation V(t) were derived. The time difference Δt between the virtual origins of 1/V(t) respectively D²(t) led to an invariant term. The equation of momentum is fulfilled. – Position and diameter of the vortex rings were determined till their decay by means of an optical system, which did not disturb the vortex rings. The experimental results in weakly turbulent ambient flow obtained by the author and by others confirm the theory very well. The ambient turbulence was nearly constant in the measuring region; its effect could be described by simply adding its viscosity to the vortex ring’s internal turbulent viscosity. The results could be represented in unified non-dimensional diagrams. Moreover, an explanation was found as to why the mean internal turbulent viscosity is constant.