一种改进后的亚格子特征尺度表达式

针对大涡模拟涡粘性亚格子模式中使用网格尺度为特征尺度存在的不足,提出一种改进的亚格子特征尺度表达式,利用结合Sagaut混合尺度模式获得的改进模式与Smagorinsky模式、Germano模式和Sagaut模式等进行对比研究.时间发展混合层和中性大气边界层的大涡模拟结果表明,改进的有效亚格子特征尺度能合理反映亚格子脉动的时空分布特性,同时扩展了适用范围.在对耗散程度、流动演化过程、可解湍流强度和可解雷诺应力等方面的模拟,改进模式优于基准模式,表明亚格子脉动有效特征尺度在大涡模拟中具有一定的理论基础和广泛的应用价值.     

平行射流大涡模拟及SGS模型的比较

高精度算法和亚网格尺度(SGS)模型是大涡模拟的关键,对声源流场和气动声场的准确预测有着重要的影响。本文采用前人具有部分DNS数据的平板射流为考核算例,使用具有4阶空间精度和3阶时间精度的高精度算法,系统地比较Smagorinsky模型(SM)、选择多尺度模型(SMSM)、以及前人采用动态Smagorinsky模型(DSM)的计算结果。对射流速度、剪切层、湍流脉动以及涡结构的多尺度特性的结果分析表明,Smagorinsky模型的黏性耗散过大,而选择多尺度模型则能揭示转捩过程中的黏性变化,其结果比DSM更接近DNS的数据。同时,SMSM是局部模型,在大规模分布式内存的并行计算方面易于编程,具有良好的应用前景。     

二维槽道湍流拟序结构的大涡模拟

本文采用大涡模拟的方法,对二维槽道湍流流动进行了数值模拟。采用Chorin的分步投影法求解大尺度涡运动的Navier-Stokes方程,小尺度涡采用三种亚格子（SGS）模式分别模拟,给出了不同亚格子涡粘性模式下的模拟结果。对固壁面采用了壁函数。模拟结果再现了二维槽道流动拟序结构的发展演变过程。通过对不同入口速度下的瞬态流场的比较,揭示了入口速度分布对流场的影响。     

双大涡模拟方法模拟提升管内气固两相流动

基于Smagorinsky涡黏模型以及颗粒动理学理论,建立了气固两相流双大涡模拟模型。考虑大涡模拟中过滤尺度的影响,给出颗粒相亚格子压力和热传导系数计算模型。考虑颗粒聚团对两相作用的影响,给出了考虑颗粒聚团作用的气固两相多尺度曳力系数模型。数值模拟了提升管内气固两相流动特性,合理地预测出了提升管内气固两相环-核流动结构。模拟结果与Knowlton等实测结果相吻合。     

湍流边界层拟序结构的大涡模拟研究

采用动力亚格子模型,利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。从瞬时速度和脉动 速度场、脉动速度相关、均方根脉动涡量分布、以及瞬时涡量场等多个方面,对湍流边界层流动的拟序结构进行了分析, 包括近壁区小尺度湍流结构和瞬态过程,如条纹结构、喷射和扫掠过程、以及近壁旋涡结构等。     

均匀入口鼓泡流化床的TFM及DEM模拟

分别运用双流体模型（TFM）及离散元模型（DEM）的软球模型，模拟了二维及三维实验尺度的均匀入口鼓泡流化床内的气固两相流动特性，其中TFM模型结合了气相κ-ε湍流模型，而DEM模型则结合气相κ-ε湍流模型或Smagorinsky亚网格涡黏模型（SGS）。通过对比本文模拟结果与他人实验及计算结果发现，TFM与DEM软球模...     

钝体燃烧室热态流场的大涡模拟

对圆盘凸台钝体燃烧室的热态流场进行了大涡模拟,利用PIV实验数据,显示了大涡模拟对平均速度及均方根速度脉动量的预测与实验中的测量值基本吻合,模拟获得的钝体后回流区、钝体边缘存在的切应力层等湍流特征信息也得到了实验验证。采用两种不同亚网格模型的大涡模拟结果与实验对比表明,在流动情况特别复杂的区域,Germano动态亚网格模型对流场,尤其是湍流脉动的预测准确性要高于标准Smagorinsky亚网格模型。     

壁面在展向作周期运动的槽道湍流的大涡模拟

分别采用3种亚格子模式:传统的Smagorinsky模式、动力Smagorinsky模式和Cui(2004)基于Kolmogorov方程所提出的新模式,对壁面在展向作周期运动的槽道湍流进行了大涡模拟,以考察这3种模式对平均运动为三维、非定常的湍流流动的模拟能力.通过对湍流基本统计量的分析,发现动力模式和新模式都可以较好地预测这种三维非定常的湍流流动;对相位平均的湍流统计量,动力模式的结果略优于新模式;传统的Smagorinsky模式对这种流动的预测结果是最差的.     

尾流的双尺度简化模型及其大涡模拟建模

叶轮机中存在着复杂的湍流流动结构,如何准确地模拟湍流流动一直是工程界关注的核心问题。由于湍流在大部分区域都处于强非均衡态,根据均衡能量传输假设建立起来的RANS湍流模型和大涡模拟亚格子模型都不再适用。非均衡湍流的一个经典的模型流动是尾流,但通常意义上的尾流往往是无数不同尺度相互作用的结果,难以进行简化分析。本文构造了一个极度简化的双尺度尾流模型进行了直接数值模拟(DNS),结果显示在转捩前期就存在反向的非均衡能量传输。分别采用各种传统的RANS湍流模型和大涡模拟亚格子模型进行了数值模拟,发现各种模型都无法准确预测尾流转捩。在考虑非均衡能量传输机制后,通过流向速度梯度扭率构造的理性亚格子应力模型可以使大涡模拟模拟的转捩位置相对于标准Smagorinsky模型有所改善;不采用涡粘假设的改进速度增量(IVI)模型可以得到更加接近DNS的结果。这些结果为湍流模型进一步的改进提供了思路。     

提升管内稠密气粒两相流动大涡模拟

本文采用双流体模型,引入颗粒动力学理论,对提升管内的稠密气粒两相流动进行了大涡模拟。采用改进的分步投影法对滤波后的方程进行显式求解,小尺度量采用Smagorinsky亚格子模式模拟。模拟结果给出的颗粒相速度分布、浓度分布与实验值基本吻合,气固两相存在速度滑移。模拟结果合理预报出了提升管内的环-核流动结构。     

Large eddy simulation of channel flow with and without periodic constrictions using the explicit algebraic subgrid-scale model

We analyse the performance of the explicit algebraic subgrid-scale (SGS) stress model (EASSM) in large eddy simulation (LES) of plane channel flow and the flow in a channel with streamwise periodic hill-shaped constrictions (periodic hill flow) which induce separation. The LESs are performed with the Code_Saturne which is an unstructured collocated finite volume solver with a second-order spatial discretisation suitable for LES of incompressible flow in complex geometries. At first, performance of the EASSM in LES of plane channel flow at two different resolutions using the Code_Saturne and a pseudo-spectral method is analysed. It is observed that the EASSM predictions of the mean velocity and Reynolds stresses are more accurate than the conventional dynamic Smagorinsky model (DSM). The results with the pseudo-spectral method were, in general, more accurate. In the second step, LES with the EASSM of flow separation in the periodic hill flow is compared to LES with the DSM, no SGS model and a highly resolved LES data using the DSM. Results show that the mean velocity profiles, the friction and pressure coefficients, the length and shape of the recirculation bubble, as well as the Reynolds stresses are considerably better predicted by the EASSM than the DSM and the no SGS model simulations. It was also observed that in some parts of the domain, the resolved strain-rate and SGS shear stress have the same sign. The DSM cannot produce a correct SGS stress in this case, in contrast to the EASSM.     

Evaluating the modulated gradient model in large eddy simulation of channel flow with OpenFOAM

Recently, a new family of subgrid-scale (SGS) models, termed as gradient-based models, has been introduced to calculate the SGS stresses in large eddy simulation (LES). In the present work, the modulated gradient model (MGM) was implemented in the OpenFOAM package, and the pimpleFoam solver was improved to be adopted with non-eddy viscosity models. The MGM is a new, nonlinear model that uses the local equilibrium hypothesis to assess the SGS kinetic energy and the velocity gradient tensor to calculate the relative weight of the different components of the SGS stress tensor. To evaluate the accuracy of the MGM along with the modified pimpleFoam solver, a turbulent channel flow was simulated at the three different frictional Reynolds numbers of 180, 395 and 590. Furthermore, the results were compared with direct numerical simulation data, as well as the numerical results obtained by the established SGS models such as the dynamic Smagorinsky model (DSM). A suitable accuracy for the first- and second-order turbulence parameters was reported. Moreover, it was demonstrated that MGM is computationally efficient compared to the DSM in treating channel flow.     

A velocity-estimation subgrid model constrained by subgrid scale dissipation

Purely dissipative eddy-viscosity subgrid models have proven very successful in large-eddy simulations (LES) at moderate resolution. Simulations at coarse resolutions where the underlying assumption of small-scale universality is not valid, warrant more advanced models. However, non-eddy viscosity models are often unstable due to the lack of sufficient dissipation. This paper proposes a simple modeling approach which incorporates the dissipative nature of existing eddy viscosity models into more physically appealing non-eddy viscosity SGS models. The key idea is to impose the SGS dissipation of the eddy viscosity model as a constraint on the non-eddy viscosity model when determining the coefficients in the non-eddy viscosity model. We propose a new subgrid scale model (RSEM), which is based on estimation of the unresolved velocity field. RSEM is developed in physical space and does not require the use of finer grids to estimate the subgrid velocity field. The model coefficient is determined such that total SGS dissipation matches that from a target SGS model in the mean or least-squares sense. The dynamic Smagorinsky model is used to provide the target dissipation. Results are shown for LES of decaying isotropic turbulence and turbulent channel flow. For isotropic turbulence, RSEM displays some level of backward dissipation, while yielding as good results as the dynamic Smagorinsky model. For channel flow, the results from RSEM are better than those from the dynamic Smagorinsky model for both statistics and instantaneous flow structures.     

横向粗糙元壁面槽道湍流的大涡模拟研究

采用大涡模拟和浸没边界法相结合对不同高度和不同间距横向粗糙元壁面槽道湍流进行了模拟,得到了光滑壁面和粗糙壁面湍流的流向平均速度分布,雷诺剪切应力,脉动速度均方根和近壁区拟序结构。结果发现横向粗糙元降低了流向平均速度,增大了流动阻力,粗糙壁面湍流的雷诺剪切应力大于光滑壁面。粗糙元降低了流向脉动速度,增强了展向和法向脉动速度。粗糙元高度越高,对湍流流动影响越大,而粗糙元间距对湍流统计特性的影响不大。粗糙壁面仍然存在着和光滑壁面类似的条带结构。     

On the effect of an anisotropy-resolving subgrid-scale model on large eddy simulation predictions of turbulent open channel flow with wall roughness

A large eddy simulation (LES) was conducted of turbulent flow in a channel with a rough wall on one side and a free surface on the other by adopting an anisotropy-resolving subgrid-scale (SGS) model. A shear Reynolds number of Re_τ = 395 was used based on the mean friction velocity and channel height. To investigate the grid dependency of the LES results caused by the SGS model, three grid resolutions were tested under the same definition of a roughness shape by using the immersed boundary method. The results obtained were compared with direct numerical simulation data with and without the wall roughness and those without the extra anisotropic term. The primary focus was on how the present anisotropic SGS model with coarser grid resolutions can properly provide the effects of roughness on the mean velocity and turbulent stresses, leading to a considerable reduction of the computational cost of LES.     

Visualisation and analysis of large-scale vortex structures in three-dimensional turbulent lid-driven cavity flow

In this paper, large eddy simulation (LES) of a three-dimensional turbulent lid-driven cavity (LDC) flow at Re = 10,000 has been performed using the multiple relaxation time lattice Boltzmann method. A Smagorinsky eddy viscosity model was used to represent the sub-grid scale stresses with appropriate wall damping. The prediction for the flow field was first validated by comparing the velocity profiles with previous experimental and LES studies, and then subsequently used to investigate the large-scale three-dimensional vortical structures in the LDC flow. The instantaneous three-dimensional coherent structures inside the cavity were visualised using the second invariant (Q), Δ criterion, λ₂ criterion, swirling strength (λ_ci) and streamwise vorticity. The vortex structures obtained using the different criteria in general agree well with each other. However, a cleaner visualisation of the large vortex structures was achieved with the λ_ci criterion and also when the visualisation is based on the vortex identification criteria expressed in terms of the swirling strength parameters. A major objective of the study was to perform a three-dimensional proper orthogonal decomposition (POD) on the fluctuating velocity fields. The higher energy POD modes efficiently extracted the large-scale vortical structures within the flow which were then visualised with the swirling strength criterion. Reconstruction of the instantaneous fluctuating velocity field using a finite number of POD modes indicated that the large-scale vortex structures did effectively approximate the large-scale motion. However, such a reduced order reconstruction of the flow based on the large-scale vortical structures was clearly not as effective in predicting the small-scale details of the fluctuating velocity field which relate to the turbulent transport.     

The effect of subgrid-scale models on the entrainment of a passive scalar in a turbulent planar jet

Classical large-eddy simulation (LES) modelling assumes that the passive subgrid-scale (SGS) models do not influence large-scale quantities, even though there is now ample evidence of this in many flows. In this work, direct numerical simulation (DNS) and large-eddy simulations of turbulent planar jets at Reynolds number Re_H = 6000 including a passive scalar with Schmidt number Sc = 0.7 are used to study the effect of several SGS models on the flow integral quantities e.g. velocity and scalar jet spreading rates. The models analysed are theSmagorinsky, dynamic Smagorinsky, shear-improved Smagorinsky and the Vreman. Detailed analysis of the thin layer bounding the turbulent and non-turbulent regions – the so-called turbulent/non-turbulent interface (TNTI) – shows that this region raises new challenges for classical SGS models. The small scales are far from equilibrium and contain a high fraction of the total kinetic energy and scalar variance, but the situation is worse for the scalar than for the velocity field. Both a-priori and a-posteriori (LES) tests show that the dynamic Smagorinsky and shear-improved models give the best results because they are able to accurately capture the correct statistics of the velocity and passive scalar fluctuations near the TNTI. The results also suggest the existence of a critical resolution Δ_x, of the order of the Taylor scale λ, which is needed for the scalar field. Coarser passive scalar LES i.e. Δ_x ≥ λ results in dramatic changes in the integral quantities. This fact is explained by the dynamics of the small scales near the jet interface.     

A priori evaluation of large eddy simulation subgrid-scale scalar flux models in isotropic passive-scalar and anisotropic buoyancy-driven homogeneous turbulence

Direct numerical simulation (DNS) of passive (non-buoyant) and active (buoyant) scalar homogeneous turbulence is carried out using a standard pseudo-spectral numerical method. The flow settings simulated include stationary forced and decaying passive-scalar turbulence, as well as decaying anisotropic active-scalar turbulence. The Schmidt number is unity in all cases. The results are compared with, and are found to be in very good agreement with, previous similar DNS studies. The well-validated DNS data are divided into 19 sets, and are employed to study different large eddy simulation (LES) subgrid-scale (SGS) models for the SGS scalar flux. The models examined include three eddy-viscosity-type models (Smagorinsky, Vreman and Sigma with a constant SGS Schmidt number), a Dynamic Structure model and two versions of the Gradient (Gradient and Modulated Gradient) model. The models are investigated with respect to their ability to predict the orientation, and the magnitude, of the SGS scalar flux. Eddy-viscosity models are found to predict the magnitude of the SGS scalar flux accurately, but are poor at predicting the orientation of the SGS scalar flux. The Dynamic Structure and Gradient models are better than eddy-viscosity models at predicting both the magnitude and direction. However, neither of them can be realised in an actual LES, without carrying additional transport equations. Based on these observations, four new models are proposed – combining directions from Dynamic Structure and Gradient models, and magnitudes from Smagorinsky and Vreman eddy-viscosity models. These models are expected to be better than eddy-viscosity and Modulated Gradient models, and this is confirmed by preliminary a posteriori tests.     

Large Eddy Simulation of Turbulent Channel Flow with 3D Roughness Using a Roughness Element Model

Large eddy simulation of turbulent channel flow with dense and small 3D roughness elements is carried out using a roughness element model Profiles of mean Reynolds stress, mean velocity and rms velocity as well as turbulent structures near the wall are obtained. The shear stress in the rough wall is larger than that in the smooth wall side and the rough wall has a larger influence on the channel flow. Profiles of mean streamwise velocity near the wall have logarithmic velocity distributions for both smooth and roughness walls, while there is a velocity decrease for the rough wall due to larger fractional drag. All the three components of rms velocities in the rough wall region are larger than that in the smooth wall region, and the roughness elements on the wall increase turbulent intensity in all directions. The s~reak spacing and average diameter of near wall quasi-s~reamwise vortices increase with the presence of roughness elements on the wall and it is shown that the rough wall induces complex and strong streamwise vortices. Results of dense and small 3D roughness elements in both turbulent statistics and structure, obtained with a relatively simple method, are found to be comparable to related experiments.