Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling effects of subgrid-scale mixture fraction variance in LES of a piloted diffusion flame

A posteriori analysis of the statistics of two large-eddy simulation (LES) solutions describing a piloted methane–air (Sandia D) flame is performed on a series of grids with progressively increased resolution reaching about 10.5 million cells. Chemical compositions, density and temperature fields are modelled with a steady flamelet approach and parametrised by the mixture fraction. The difference between the LES solutions arises from a different numerical treatment of the subgrid scale (SGS) mixture fraction variance – an important quantity of interest in non-premixed combustion modelling. In the first case (model I), the variance transport equation is solved directly, while in the second (model II), an equation for the square of the mixture fraction is solved, and the variance is computed from its definition. The comparison of the LES solutions is based on the convergence properties of their statistics with respect to the turbulence resolution length scale. The dependence of the LES statistics is analysed for velocity and the mixture fraction fields, and tested for convergence. For the most part, the statistics converge for the finest grids, but the variance of the mixture fraction shows some residual grid dependence in the high-gradient regions of the jet near field. The SGS variance given by model I exhibits realisability everywhere, whereas in regions of the flame model II is non-realisable, predicting negative variances. Furthermore, the LES statistics of model I exhibit superior convergence behaviour.     

水动力作用对污染物在河流水沙两相中分配的影响研究进展

通过总结近年来有关水动力作用对河流中泥沙吸附/释放重金属、磷等污染物影响研究成果,分析了水流紊动强度和流速对污染物在河流水沙两相中分配的影响规律。当水流紊动强度较低时,泥沙主要以床沙形式存在,泥沙对污染物的吸附/释放速率和强度较低;随着紊动强度增大,水沙界面切应力随之增大并促使泥沙悬浮,污染物与悬浮泥沙颗粒接触面积增大,使得泥沙对污染物的吸附/释放速率和强度显著增加。流速对污染物在水沙两相间分配的影响较为复杂。泥沙静止时,流速的增大会减小边界层厚度,增大水体溶解氧含量、氧化还原电位等参数,并使污染物在水体里的扩散由分子扩散转变为以紊动扩散为主,增强污染物在水-沙界面的交换通量,对污染物在水沙两相间的分配产生影响;随着流速进一步增大,床沙逐渐起动、悬浮,此时除上述因素外,泥沙运动状态、悬浮颗粒间碰撞强度、悬浮物的絮凝等均会对泥沙吸附/释放污染物产生影响,由于影响因素较多,作用机制复杂,目前关于泥沙运动对污染物在水沙两相间分配的影响所得结论仍存在分歧。相应的机理研究及数学模型的建立均做了较大简化。综观当前的研究成果,水沙运动及其化学生物过程等对污染物在河流水沙两相中的分配作用机理和耗氧有机物、有毒有机物等污染物与泥沙的作用机理是这一领域的未来研究重点。     

GPU‐accelerated direct numerical simulations of decaying compressible turbulence employing a GKM‐based solver

Gas Kinetic Method‐based flow solvers have become popular in recent years owing to their robustness in simulating high Mach number compressible flows. We evaluate the performance of the newly developed analytical gas kinetic method (AGKM) by Xuan et al. in performing direct numerical simulation of canonical compressible turbulent flow on graphical processing unit (GPU)s. We find that for a range of turbulent Mach numbers, AGKM results shows excellent agreement with high order accurate results obtained with traditional Navier–Stokes solvers in terms of key turbulence statistics. Further, AGKM is found to be more efficient as compared with the traditional gas kinetic method for GPU implementation. We present a brief overview of the optimizations performed on NVIDIA K20 GPU and show that GPU optimizations boost the speedup up‐to 40x as compared with single core CPU computations. Hence, AGKM can be used as an efficient method for performing fast and accurate direct numerical simulations of compressible turbulent flows on simple GPU‐based workstations. Copyright © 2016 John Wiley & Sons, Ltd.     

某冰水堆积体层流-紊流过渡状态渗流特性试验研究

澜沧江上游两岸河谷地带广泛分布一套规模巨大、形成于第四纪更新世中晚期的冰水堆积体,这类冰水堆积体往往具有较特别的渗流特性。通过对澜沧江某冰水堆积体的结构特征调查和颗粒组成试验分析,从沉积相的角度,认为该冰水堆积体可分为两层:表层厚约40 m的混合巨粒土,和底部厚约17 m的粗颗粒土。通过对这两层冰水堆积体的渗流特点试验研究,结果表明地下水在表层混合巨粒土中的渗流方式主要表现为层流-紊流过渡阶段的特点,渗流速度与水力梯度成幂函数关系;而在底部粗颗粒土中的运移方式主要表现为层流特点。根据混合巨粒土的渗流特点,文中提出了两种计算其渗透系数的方法,即根据渗透速度和水力梯度关系曲线估算渗透系数,或采用文中提出的公式计算渗透系数。     

An efficient low‐dissipative WENO filter scheme

This paper presents a simple and efficient procedure developed for tracing discontinuities in flow fields. Numerical experiments are carried out to test the new sensor coupled with the associated nonlinear WENO dissipation filter developed to suppress the numerical dissipation. The tests show that, for a problem containing shocks and vortices, the implementation of the new sensor and an optimized WENO scheme can obtain a stable and highly resolved solution. The numerical experiments demonstrated that the new filter scheme performs efficiently both in parallel and serial running for one‐dimensional inviscid flow problems. Direct numerical simulation of a Mach 5 turbulent boundary layer over a flat plate was carried out to demonstrate the applicability of the scheme to the DNS practices. Copyright © 2011 John Wiley & Sons, Ltd.     

Calculation of Turbulent Flow with Pseudo-spectral Method

This paper deals with the numerical simulation of incompressible turbulent boundary flow of a flat plate with the pseudo-spectral matrix method. In order to appear more than 10 nodes in the turbulent base-stratum and transition of 43×43 computational grids,a coordinate transformation is put up from physical panel to computational panel. Several zero turbulent models are computed comparatively. The results are credible when comparing with the previous methods.     

A Novel Scheme for Calculating Pressure Fluctuations Used in Lagrangian-Eulerian Particulate Dispersion Codes

A novel scheme for pressure fluctuations in turbulent flows is developed. The pressure fluctuations are sensitive parameter in some of the fluid phenomena. In the computational methods and modeling turbulence flow, the pressure fluctuations are eliminated after averaging of the Navier-Stokes equations, and only average pressure could be calculated. In this research, the Reynolds-averaged Navier-Stokes equations are computed using SIMPLE method. The Reynolds stress transport model (RSTM) is used to determine the Reynolds stresses and the flow details. The velocity fluctuations are simulated using the Kraichnan model. The Poisson equation for the pressure fluctuations is obtained by taking the divergence of the incompressible momentum equation and algebraic operations, and this equation is numerically solved by finite difference method. The effects of Reynolds number on the pressure fluctuations are studied.     

Effect of pressure gradients on the different stages of roughness induced boundary layer transition

The physical mechanisms of roughness-induced transition (RIT) in pressure gradient boundary layers are studied using direct numerical simulations. Recent investigations have examined RIT processes in zero-pressure-gradient boundary layers (Suryanarayanan et al., 2019). The present study uses a vorticity dynamics point of view to examine how these processes are altered by a locally accelerating or decelerating flow that strains the vorticity field and creates a net vorticity flux at the wall. Flow acceleration is imposed on specific streamwise extents of the flow. This provides an understanding about how the fundamental mechanisms in different stages of RIT are affected by pressure gradients. The present results suggest that both lift-up and subsequent amplification of the unsteady perturbations are mitigated by flow acceleration. The effect on lift-up is explained by the compression (i.e. large negative value of the stretching term) of the wall-normal vorticity by negative dv/dy. Consistent with earlier experimental observations on spots and wedges, favorable pressure gradients reduce turbulent wedge spreading and nearly arrest the spreading when sufficiently strong. This result is also explained in terms of vorticity dynamics.