Flow-thermodynamics interactions in decaying anisotropic compressible turbulence with imposed temperature fluctuations

The fundamental nature of the non-linear flow-thermodynamics interactions in a compressible turbulent flow with imposed temperature fluctuations is investigated. Direct numerical simulations (DNS) of decaying anisotropic compressible turbulence (turbulent Mach number 0.06–0.6) with imposed temperature fluctuations are performed to examine: (i) interactions between solenoidal and dilatational kinetic energy; (ii) partition between dilatational kinetic energy and thermodynamic potential energy; and (iii) redistribution of solenoidal and dilatational kinetic energy among the various Reynolds stress components. It is found that solenoidal kinetic energy levels and return-to-isotropy are weakly dependent on Mach number but independent of imposed temperature fluctuations in the parameter range studied. The dilatational kinetic energy generated is proportional to the square of the pressure fluctuations associated with the initial solenoidal and temperature fluctuations and thus a strong function of Mach number and heat release intensity. The energy exchange between dilatational kinetic and potential energy is driven by a strong proclivity toward equipartition. Consequently, the dynamics of pressure-dilatation ( ${\overline{pd}}$ ), which is the mechanism of this energy exchange between dilatational and potential energies, is dictated entirely by the requirement to impose energy equipartition. Based on the results, we provide a physical picture of the solenoidal–dilatational–potential energy interactions and the action of pressure-dilatation. The identification of the fundamental precepts underlying the various interactions is of great utility for turbulence closure model development.     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

DNS of viscoelastic turbulent channel flow with rectangular orifice at low Reynolds number

Direct numerical simulations of turbulent viscoelastic-fluid flow in a channel with a rectangular orifice were performed to investigate the influence of viscoelasticity on turbulence statistics and turbulent structures downstream of the orifice. The geometry considered is periodic rectangular orifices with 1:2 expansion. The constitutive equation follows the Giesekus model, valid for polymer (or surfactant) solutions, which are generally capable of reducing the turbulent frictional drag in a smooth channel. The friction Reynolds number and the Weissenberg number were set to 100 and 20-30, respectively. A drag reduction of about 20% was achieved in the viscoelastic flows. The onset Reynolds number for the transition from a symmetric to an asymmetric state was found to be shifted to higher values than that for the Newtonian flow. In the viscoelastic flow, the turbulent kinetic energy was decreased and fewer turbulent eddies were observed, as the Kelvin-Helmholtz vortices were quickly damped. Away from the orifice, quasi-streamwise vortices in the viscoelastic flow were sustained for a longer period, accompanied by energy exchange from elastic energy of the viscoelastic fluid to kinetic energy.     

A two‐scale low‐Reynolds number turbulence model

In this study, a two‐scale low‐Reynolds number turbulence model is proposed. The Kolmogorov turbulence time scale, based on fluid kinematic viscosity and the dissipation rate of turbulent kinetic energy (ν, ε), is adopted to address the viscous effects and the rapid increasing of dissipation rate in the near‐wall region. As a wall is approached, the turbulence time scale transits smoothly from a turbulent kinetic energy based (k, ε) scale to a (ν, ε) scale. The damping functions of the low‐Reynolds number models can thus be simplified and the near‐wall turbulence characteristics, such as the ε distribution, are correctly reproduced. The proposed two‐scale low‐Reynolds number turbulence model is first examined in detail by predicting a two‐dimensional channel flow, and then it is applied to predict a backward‐facing step flow. Numerical results are compared with the direct numerical simulation (DNS) budgets, experimental data and the model results of Chien, and Lam and Bremhorst respectively. It is proved that the proposed two‐scale model indeed improves the predictions of the turbulent flows considered. Copyright © 2000 John Wiley & Sons, Ltd.     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.     

可压缩均匀剪切湍流的二阶矩模拟

基于对可压缩湍流中脉动压力场和脉动速度场特征的理论分析以及ＤＮＳ结果,建立了可均匀剪切湍流中压力－变形率关联的压缩性修正模式,应用这个模式,加上Ｓａｒｋａｒ等建立的脉动体胀率项（ｄｉｌａｔａｔｉｏｎａｌ ｔｅｒｍｓ）的模式,预测可压缩均匀剪切湍流随时间的发展,所得雷诺应力各是性张量的平衡值与Ｂｌａｉｓｄｅｌｌ等的ＤＮＳ数据非常一致。这个模式准确地预测出均匀剪切湍流中压缩性导致的雷诺应力结构的“流向     

Effect of Roughness on Wall-Bounded Turbulence

Direct numerical simulation of turbulent incompressible plane-channel flow between a smooth wall and one covered with regular three-dimensional roughness elements is performed. While the impact of roughness on the mean-velocity profile of turbulent wall layers is well understood, at least qualitatively, the manner in which other features are affected, especially in the outer layer, has been more controversial. We compare results from the smooth- and rough-wall sides of the channel for three different roughness heights of h ⁺= 5.4, 10.8, and 21.6 for Re _τ of 400, to isolate the effects of the roughness on turbulent statistics and the instantaneous turbulence structure at large and small scales. We focus on the interaction between the near-wall and outer-layer regions, in particular the extent to which the near-wall behavior influences the flow further away from the surface. Roughness tends to increase the intensity of the velocity and vorticity fluctuations in the inner layer. In the outer layer, although the roughness alters the velocity fluctuations, the vorticity fluctuations are relatively unaffected. The higher-order moments and the energy budgets demonstrate significant differences between the smooth-wall and rough-wall sides in the processes associated with the wall-normal fluxes of the Reynolds shear stresses and turbulence kinetic energy. The length scales and flow dynamics in the roughness sublayer, the spatially inhomogeneous layer within which the flow is directly influenced by the individual roughness elements, are also examined. Alternative mechanisms involved in producing and maintaining near-wall turbulence in rough-wall boundary layers are also considered. We find that the strength of the inner/outer-layer interactions are greatly affected by the size of the roughness elements.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

基于雷诺应力模型的平衡大气边界层模拟

基于雷诺应力湍流模型（简称RSM模型）,研究了平衡大气边界层风场数值模拟问题．假设流体不可压,且不计雷诺应力输运方程中的对流项、浮力产生项、系统旋转产生项和扩散项,在准各向同性的条件下,推导出RSM模型湍动能k的表达式是标准k-ε模型k常数表达式的0.893倍．考虑k沿高度变化的修正,根据在标准k-ε模型中满足水平均匀性的湍流来流边界条件,提出在RSM模型中产生平衡大气边界层的湍流来流边界条件．基于空风洞的数值模拟结果表明,与工程上常用的湍流来流边界条件相比,基于本文提出的湍流来流边界条件得到的风场水平均匀性更优,且在整个流域内,得到的雷诺应力剖面更合适．从而验证了该湍流来流边界条件的适用性．     

Effect of Particle Clusters on Carrier Flow Turbulence: A Direct Numerical Simulation Study

Experiments indicate that particle clusters that form in fluidized–bed risers can enhance gas-phase velocity fluctuations. Direct numerical simulations (DNS) of turbulent flow past uniform and clustered configurations of fixed particle assemblies at the same solid volume fraction are performed to gain insight into particle clustering effects on gas-phase turbulence, and to guide model development. The DNS approach is based on a discrete-time, direct-forcing immersed boundary method (IBM) that imposes no-slip and no-penetration boundary conditions on each particle’s surface. Results are reported for mean flow Reynolds number Re _p ?=?50 and the ratio of the particle diameter d _p to Kolmogorov scale is 5.5. The DNS confirm experimental observations that the clustered configurations enhance the level of fluid-phase turbulent kinetic energy (TKE) more than the uniform configurations, and this increase is found to arise from a lower dissipation rate in the clustered particle configuration. The simulations also reveal that the particle-fluid interaction results in significantly anisotropic fluid-phase turbulence, the source of which is traced to the anisotropic nature of the interphase TKE transfer and dissipation tensors. This study indicates that when particles are larger than the Kolmogorov scale (d _p ?>?η), modeling the fluid-phase TKE alone may not be adequate to capture the underlying physics in multiphase turbulence because the Reynolds stress is anisotropic. It also shows that multiphase turbulence models should consider the effect of particle clustering in the dissipation model.     

Large eddy simulation of temporally developing juncture flows

Large eddy simulation (LES) results are reported for temporally developing solid–solid and solid–rigid-lid juncture flows. A MacCormack-type scheme that is second-order in time, and fourth-order in space for the convective terms and second-order in space for the viscous terms, is used. The simulations are obtained for a low subsonic Mach number. The subgrid-scale stresses (SGS) are modeled using the dynamic modeling procedure. The turbulent flow field generated on a flat-plate boundary layer is used to initialize the juncture flow simulations. The results of the flat-plate boundary layer simulations are validated with experimental and direct numerical simulations (DNS) data. In juncture flow simulations, the presence of an adjacent solid-wall/rigid-lid boundary altered the mean and the turbulent field, setting up gradients in the anisotropy of normal Reynolds stresses resulting in the formation of turbulence-induced secondary vortices. The relative size of these secondary vortices and the distribution of mean and turbulent quantities are in qualitative agreement with the experimental observations for the solid–solid juncture. The overall distribution of the mean and turbulence quantities showed close resemblance between the solid–solid and the solid–rigid-lid junctures; except for the absence of a second vortical region near the rigid-lid boundary. In agreement with the experimental observations, it was found that the normalized anisotropy term exhibited similarity when plotted against the distance from the boundary, regardless of the type of boundary, i.e. solid-wall or rigid-lid. The turbulent kinetic energy increased near the rigid-lid boundary. While the surface normal velocity fluctuations decreased to zero at the rigid-lid boundary, the other two velocity components showed an increase in their energy, which is also consistent with the experimental observations. © 1998 John Wiley & Sons, Ltd.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Large-eddy simulation of counter-rotating Taylor–Couette flow: The effects of angular velocity and eccentricity

In the present work, turbulent flow in the annulus of a counter-rotating Taylor-Couette (CRTC) system is studied using large-eddy simulation. The numerical methodology employed is validated, for both the mean and second-order statistics, with the direct numerical simulation (DNS) data available in the literature, for a range of Reynolds numbers from 500 to 4000. Thereafter, turbulent flow occurring in this system at Reynolds numbers of 8000 and 16000 are studied, and the results obtained are analyzed using mean and second-order statistics, vortical structures, velocity vector plots and power energy spectra. Further, the spatio-temporal variation of azimuthal velocity, extracted near the inner cylinder, shows the existence of herringbone like patterns similar to that observed in the previous studies. The effect of eccentricity of the inner cylinder with respect to the outer cylinder is studied, on the turbulent flow in the CRTC system, for two different eccentricity ratios of 0.2 and 0.5 and for two different Reynolds numbers of 1500 and 4000. The results of the eccentric CRTC are analyzed using contours of pressure, mean and second-order statistics, velocity vectors, vortical structures, and turbulence anisotropy maps. It is observed from the eccentric CRTC simulations that the smaller-gap region seems to contain higher amplitude fluctuations and more vortical structures when compared with the larger-gap region. The mean turbulent kinetic energy contours do not change qualitatively with the Reynolds number, however, quantitatively a higher turbulent kinetic energy is observed in the higher Reynolds number case of 4000.     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

可压缩各向同性衰减湍流直接数值模拟研究

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动.     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

An Analysis of Turbulent Motions In and Around a Differentially Forced Density Interface

The Rapid-Distortion-Theory-based analysis proposed by Fernando and Hunt [1] is extended to study the nature of turbulence in and around a density interface sandwiched between turbulent layers with dissimilar properties. It is shown that interfacial motions consist of low-frequency, resonantly excited, nonlinear internal waves and high-frequency, linear internal waves driven by background turbulence. Based on the assumptions that (i) all resonant waves and some nonresonant waves having frequencies close to the resonant frequencies grow rapidly, break, and cause interfacial mixing, (ii) the spectral amplitude of the vertical velocity in the wave-breaking regime is constant, and (iii) kinetic energy is equipartitioned between linear and nonlinear breaking wave regimes, the r.m.s. vertical velocity at the interface and the turbulent kinetic energy flux into the interface are calculated. The migration velocity of the interface is calculated using the additional assumption that the buoyancy flux into a given turbulent layer is a fixed fraction of the turbulent kinetic energy flux supplied to the interface by the same layer. The calculations are found to be in good agreement with the entrainment data obtained in previous laboratory experiments in the parameter regime where the interface is dominated by internal wave dynamics. Received 23 July 1997 and accepted 8 January 1999     

Predicting Turbulent Spectra in Drag-reduced Flows

In the present work we describe how turbulent skin-friction drag reduction obtained through near-wall turbulence manipulation modifies the spectral content of turbulent fluctuations and Reynolds shear stress with focus on the largest scales. Direct Numerical Simulations (DNS) of turbulent channels up to Re _τ = 1000 are performed in which drag reduction is achieved either via artificially removing wall-normal turbulent fluctuations in the vicinity of the wall or via streamwise-travelling waves of spanwise wall velocity. This near-wall turbulence manipulation is shown to modify turbulent spectra in a broad range of scales throughout the whole channel. Above the buffer layer, the observed changes can be predicted, exploiting the vertical shift of the logarithmic portion of the mean streamwise velocity profile, which is a classic performance measure for wall roughness or drag-reducing riblets. A simple model is developed for predicting the large-scale contribution to turbulent fluctuation and Reynolds shear stress spectra in drag-reduced turbulent channels in which a flow control acts at the wall. Any drag-reducing control that successfully interacts with large scales should deviate from the predictions of the present model, making it a useful benchmark for assessing the capability of a control to affect large scales directly.     

Direct numerical simulation of a freely decaying turbulent interfacial flow

Whereas Large Eddy Simulation (LES) of single-phase flows is already widely used in the CFD world, even for industrial applications, LES of two-phase interfacial flows, i.e. two-phase flows where an interface separates liquid and gas phases, still remains a challenging task. The main issue is the development of subgrid scale models well suited for two-phase interfacial flows. The aim of this work is to generate a detailed data base from direct numerical simulation (DNS) of two-phase interfacial flows in order to clearly understand interactions between small turbulent scales and the interface separating the two phases. This work is a first contribution in the study of the interface/turbulence interaction in the configuration where the interface is widely deformed and where both phases are resolved by DNS. To do this, the interaction between an initially plane interface and a freely decaying homogeneous isotropic turbulence (HIT) is studied. The densities and viscosities are the same for both phases in order to focus on the effect of the surface tension coefficient. Comparisons with existing theories built on wall-bounded or free-surface turbulence are carried out. To understand energy transfers between the interfacial energy and the turbulent one, PDFs of the droplet sizes distribution are calculated. An energy budget is carried out and turbulent statistics are performed including the distance to the interface as a parameter. A spectral analysis is achieved to highlight the energy transfer between turbulent scales of different sizes. The originality of this work is the study of the interface/turbulence interactions in the case of a widely deformed interface evolving in a turbulent flow.