基于卡门尺度和滤波方法的SST方程改进

基于滤波方法和卡门尺度对原始剪切应力输运(shear stress transport, SST)湍流模型进行了改进,提出了一种卡门尺度修正的滤波SST 方法. 湍流多尺度效应必须在分离流场模拟中给予反映,该方法减弱了雷诺平均(Reynolds averaged Navier-Stokes, RANS)方法时间平均特性对于流场脉动量的压迫作用,在流场中引入了大涡模拟(large eddy simulation, LES)方法的亚格子模型,形成一种新型的脱体涡模拟方法(detached eddy simulation,DES)方法;同时,为了降低原始DES方法在网格加密过程中产生网格诱发的雷诺应力损耗,利用卡门尺度对滤波因子进行修正. 平板边界层算例中,卡门尺度对于RANS方法的跟随性远远强于DES方法,在边界层内的速度型和RANS方法吻合很好,而DES方法在加密过程中速度型的鲁棒性较差,说明卡门尺度在有效地保护了边界层内使用RANS求解,降低速度型偏离对数率现象的产生;HGR-01翼型算例证明BY-SST方法可以有效的避免网格诱导分离现象的产生;证明BY-SST方法在分离流动中的精度高于DES方法.     

侧向多喷口干扰复杂流动数值模拟研究

采用具有高分辨率的NND格式，通过数值求解N-S方程对典型外形多喷口侧向喷流 复杂干扰流动进行了数值模拟. 为了提高计算效率，采用了LU-SGS隐式算法. 采用分块对接 网格技术，生成高质量的贴体计算网格，精确模拟喷口截面. 对比分析了不同计算格式、限 制器形式、网格拓扑及流动形态(层流与湍流)对喷流干扰流场结构和压力分布特性的影响， 研究和分析了喷口附近流场的涡系结构、波系结构和喷流干扰引起的气动力特性. 在上述研 究的基础上，针对典型飞行器外形的侧向喷流干扰特性进行了详细的数值模拟，得到了喷口 参数(喷口位置、数目等)及来流条件对喷流干扰流场结构、气动力特性的影响规律，并对 其流动机理进行了相应的分析. 研究表明，发展的针对多喷口侧喷干扰的数值计算方法 是成功的，可以应用于飞行器侧向喷流干扰的流场结构分析及气动力特性数值预测.     

基于Navier-Stokes方程残差的隐式大涡模拟有限元模型

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.      

A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

90°弯管内流动的理论模型及流动特性的数值研究

从三维不可压缩雷诺时均Navier-Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk-ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。     

基于组合神经网络的雷诺平均湍流模型 多次修正方法

求解雷诺平均(Reynolds-averaged Navier-Stokes, RANS)方程依然是工程应用中有效且实用的方法, 但对雷诺应力建模的不确定性会导致该方法的预测精度具有很大差异. 随着人工智能的发展, 湍流闭合模型结合机器学习元素的数据驱动方法被认为是提高RANS模型预测性能的有效手段, 然而这种数据驱动方法的稳定性和预测精度仍有待进一步提高. 本文通过构建一个全连接神经网络对RANS方程中的涡黏系数进行预测以实现雷诺应力的隐式求解,该神经网络记作涡黏系数神经网络(eddy viscosity neural network, EVNN). 此外, 也使用张量基神经网络(tensor basis neural network, TBNN)预测未封闭量与解析量之间的高阶涡黏关系, 并利用基张量保证伽利略不变性. 最后, 采用多次修正的策略实现修正模型对流场预测的精度闭环. 上述方法使用大涡模拟(large eddy simulation, LES)方法产生的高保真数据, 以及RANS模拟获得的基线数据对由EVNN和TBNN组合的神经网络进行训练, 然后用训练好的模型预测新的RANS模拟的流场. 通过与高保真LES结果进行对比, 结果表明, 相比于原始RANS模型, 修正模型对后验速度场、下壁面平均压力系数和摩擦力系数的预测精度均有较大提升. 可以发现对雷诺应力线性部分的隐式处理可以增强数值求解的稳定性, 对雷诺应力非线性部分的修正可以提升模型对流场各向异性特征预测的性能, 并且多次修正后的模型表现出更高的预测精度. 因此, 该算法在数据驱动湍流建模和工程应用中具有很大的应用潜力.      

基于DES方法平底水力旋流器内部流场特征的数值分析

采用分离涡模拟方法(Detached Eddy Simulation,DES)针对平底水力旋流器的特殊结构进行了清水流场三维湍流流动的数值模拟.与RNG k-ε模型、雷诺应力模型的计算结果比较后发现,DES方法所预测出的切向速度结果与实验结果吻合最好,从而验证了DES方法在水力旋流器数值模拟方面的可行性.进一步分析了在一定入口流量下平底水力旋流器内部流场的切向速度、轴向速度、径向速度及静压力的分布特征以及局部流动现象,获得了有关平底水力旋流器内流场流动结构的丰富的信息,讨论了圆柱结构内的流动特征对固液分离可能存在的影响.分析表明平底水力旋流器的内流场呈现出兰金涡的特点,并且为非对称分布,其内部的能量损耗主要集中在受强制涡影响的中心轴线附近区域内.本文结果为进一步研究平底水力旋流器的流场特征提供了参考.     

大涡模拟建筑物对近源大气污染物扩散的影响

采用大涡模拟(LES)方法对一个由上风低源排放的烟流在孤立建筑物周围扩散的现象进行了数值模拟。模拟中,离散网格的尺度控制在100倍Kolmogorov湍流尺度的量级以内,亚格子应力由Smagorinsky模型加以模拟,并采用了随机数方法生成符合指定湍流特征的入流边界条件。数值模拟结果表明:建筑物尾流区平均浓度场的分布和大小不仅取决于局地流场的影响,在很大程度上也受到因建筑物阻挡出现的烟流抬升和马蹄涡对烟流的捕获的影响。通过对比风洞试验数据显示:大涡模拟方法能较好地再现气流绕建筑物的流动形态和建筑物周围的浓度分布,因而具有广泛的应用前景。     

可压缩多介质粘性流体和湍流的大涡模拟

在可压缩多介质粘性流体动力学高精度计算方法MVPPM(multi-viscous-fluid piecewise parabolicmethod)基础上,引入Smagorinsky和Vreman亚格子湍流模型,采用大涡数值模拟方法求解可压缩粘性流体NS(Navier-Stokes)方程,给出适用于可压缩多介质流体界面不稳定性发展演化至湍流阶段的计算方法和二维计算程序MVFT(multi-viscosity-fluid and turbulence)。在2种亚格子湍流模型下计算了LANL(Los Ala-mos National Laboratory)激波管单气柱RM不稳定性实验,分析了气柱的形状、流场速度以及涡的特征,通过与LANL实验和计算结果的比较可知,Vreman模型略优于Smagorinsky模型,MVFT方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

LES/RANS方法混合函数特性研究

将两方程k-ω SST湍流模型和Sagaut的混合尺度亚格子模型通过一个混合函数相结合, 构造出一种混合大涡/雷诺平均N-S方程模拟方法(hybird large eddy simulation/reynolds-averaged navier-stokes, Hybrid LES/RANS), 采用这种混合模拟方法结合5阶WENO格式对Ma=2.8平板湍流边界层进行了数值模拟, 并在计算区域上游入口处采用“回收/调节”方法生成湍流脉动边界条件, 通过考查RANS区域向LES区域的过渡参数及网格分辨率对这种混合模拟方法进行了评价. 计算结果表明: 该文采用的混合模拟方法可以捕捉到湍流边界层中的大尺度结构且入口边界层平均参数不会发生漂移, 混合函数应当将RANS区域和LES区域的过渡点设置在对数律层和尾迹律层的交界处, 而过渡应当迅速以获得正确的雷诺剪切应力分布, 在该文采用的模型及数值方法的条件下, 流向及展向的网格小至与Escudier混合长相当时, 能够获得可以接受的脉动速度的单点-二阶统计值.     

Hybrid URANS/LES Turbulence Simulation of Vortex Shedding Behind a Triangular Flameholder

In this study a detached eddy simulation (DES) model, which belongs to the group of hybrid URANS/LES turbulence models, is used for the simulation of vortex shedding behind a triangular obstacle. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS shear-stress transport (SST) model is used. In the other regions with higher grid resolution a LES model, which uses a transport equation for the turbulent subgrid energy, is applied. The DES model is first investigated for two standard test cases, namely decaying homogeneous isotropic turbulence and the backward facing step, respectively. For the decaying homogeneous isotropic turbulence test case the evolution of the energy spectra in wavenumber space for different times are studied for both the DES and a Smagorinsky type LES model. Different grid resolutions are analyzed with a special emphasis on the modeling constant connecting the filter length scale to the grid size. The results are compared to experimental data. The backward facing step test case is used to study the model behavior for a case with a transition region between a RANS modeling approach close to the wall and LES based modeling in the intense shear flow region. The final application is the simulation of the vortex shedding behind a triangular obstacle. First, the influence of the inlet condition formulation is studied in detail as they can have a significant influence especially for LES based models. Detailed comparisons between simulation and experiment for the flow structure past the obstacle and statistical quantities such as the shedding frequency are shown. Finally the additional temporal and spatial information provided by the DES model is used to show the predicted anisotropy of turbulence.     

Non-zonal detached eddy simulation coupled with a steady RANS solver in the wall region

Xiao and Jenny (2012) proposed an interesting hybrid LES/RANS method in which they use two solvers and solve the RANS and LES equations in the entire computational domain. In the present work this method is simplified and used as a hybrid RANS-LES method, a wall-modeled LES. The two solvers are employed in the entire domain. Near the walls, the flow is governed by the steady RANS solver; drift terms are added to the DES equations to ensure that the time-averaged DES fields agree with the steady RANS field. Away from the walls, the flow is governed by the DES solver; in this region, the RANS field is set to the time-averaged LES field. The disadvantage of traditional DES models is that the RANS models in the near-wall region – which originally were developed and tuned for steady RANS – are used as URANS models where a large part of the turbulence is resolved. In the present method – where steady RANS is used in the near-wall region – the RANS turbulence models are used in a context for which they were developed. In standard DES methods, the near-wall accuracy can be degraded by the unsteady agitation coming from the LES region. It may in the present method be worth while to use an accurate, advanced RANS model. The EARSM model is used in the steady RANS solver. The new method is called NZ S-DES . It is found to substantially improve the predicting capability of the standard DES. A great advantage of the new model is that it is insensitive to the location of the RANS-LES interface.     

The State of the Art of Hybrid RANS/LES Modeling for the Simulation of Turbulent Flows

This review presents the state of the art of hybrid RANS/LES modeling for the simulation of turbulent flows. After recalling the modeling used in RANS and LES methodologies, we propose in a first step a theoretical formalism developed in the spectral space that allows to unify the RANS and LES methods from a physical standpoint. In a second step, we discuss the principle of the hybrid RANS/LES methods capable of representing a RANS-type behavior in the vicinity of a solid boundary and an LES-type behavior far away from the wall boundary. Then, we analyze the principal hybrid RANS/LES methods usually used to perform numerical simulation of turbulent flows encountered in engineering applications. In particular, we investigate the very large eddy simulation (VLES), the detached eddy simulation (DES), the partially integrated transport modeling (PITM) method, the partially averaged Navier-Stokes (PANS) method, and the scale adaptive simulation (SAS) from a physical point of view. Finally, we establish the connection between these methods and more precisely, the link between PITM and PANS as well as DES and PITM showing that these methods that have been built by different ways, practical or theoretical manners have common points of comparison. It is the opinion of the author to consider that the most appropriate method for a particular application will depend on the expectations of the engineer and the computational resources the user is prepared to expend on the problem.     

Simulation of wing-body junction flows with hybrid RANS/LES methods

In this paper, flows past two wing-body junctions, the Rood at zero angle of attack and NASA TN D-712 at 12.5° angle of attack, are investigated with two Reynolds-Averaged Navier-Stokes (RANS) and large eddy simulation (LES) hybrid methods. One is detached eddy simulation (DES) and the other is delayed-DES, both are based on a weakly nonlinear two-equation k–ω model. While the RANS method can predict the mean flow behaviours reasonably accurately, its performance for the turbulent kinetic energy and shear stress, as compared with available experimental data, is not satisfactory. DES, through introducing a length scale in the dissipation terms of the turbulent kinetic energy equation, delivers flow separation, a vortex or the onset of vortex breakdown too early. DDES, with its delayed effect, shows a great improvement in flow structures and turbulence characteristics, and agrees well with measurements.     

Prediction of separation flows around a 6：1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier–Stokes (RANS) and large-eddy-simulation (LES) methods for the separated flows at high angles of attack around a 6:1 prolate spheroid. The RANS/LES hybrid methods studied in this work include the detached eddy simulation (DES) based on Spalart–Allmaras (S–A), Menter’s k–ω shear-stress-transport (SST) and k–ω with weakly nonlinear eddy viscosity formulation (Wilcox–Durbin+, WD+) models and the zonal-RANS/LES methods based on the SST and WD+ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower–upper symmetric-Gauss–Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows. The project supported by the National Natural Science Foundation of China (10502030 and 90505005).     

Experiences with two detached eddy simulation approaches through use of a high-order finite volume solver

In the present study, two advanced detached eddy simulation (DES) approaches, shear-layer-adapted delayed DES and zonal DES in mode II, which are known to help transition from RANS to LES mode, are employed in various flow problems in conjunction with a high-order finite volume solver. The numerical scheme, being only applicable on structured grids, has low-dissipation and low-dispersion features. Such features benefit mostly in the LES mode, minimizing the interference of numerical diffusion with subgrid eddy viscosity. First, corresponding subgrid models are validated via decaying homogeneous turbulence benchmark case. Then, a channel flow problem is chosen to examine these models in attached flow situations. Finally, flow around an airfoil at low Reynolds number is solved using the shear-layer-adapted delayed DES approach only, in an aim to obtain trailing-edge noise spectrum at an observer location. Despite some log-layer mismatch over turbulent boundary layers, which is typical of most DES methods, the combined application of high-resolution numerical method and advanced DES approaches, which are implemented on a stabilized Spalart-Allmaras turbulence model, shows merit in resolution of turbulence in regions of interest.     

High‐order detached eddy simulation,zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth‐order discretization scheme is used and is validated using the test case of a two‐dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS–NLES approach does not suffer from this trait. For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton–Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS–NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. Copyright © 2013 John Wiley & Sons, Ltd.     

Studies of the unsteady supersonic base flows around three afterbodies

Unsteady supersonic base flows around three afterbodies, cylindrical （Cy）, boattailed （BT） and three-step （MS）, are investigated in this paper. Reynolds-averaged Navier-Stokes （RANS） and two RANS/LES （large-eddy simulation） hybrid methods, detached eddy simulation （DES） and delayed-DES （DDES）, are used to predict the base flow characteristics around the baseline Cy afterbody. All the RANS and hybrid methods are based on the two-equation SST （shear-stress transport） model with compressible corrections （CC）. According to the comparison of measurements, both DES and DDES can produce more satisfactory results than RANS. RANS can only present the ＂stable＂ flow pat- terns, while the hybrid methods can demonstrate unsteady flow structures. DDES and DES results are little different from one another although the latter exhibits better agreement with the experiment. DES is taken to investigate the 5° BT and three-step afterbodies. The mean flow data and the instantaneous turbulent coherent structures are compared against available measurements.     

On the Application of LES to Seal Geometries

Five Large Eddy Simulation (LES) and hybrid RANS-NLES (Reynolds-Averaged Navier-Stokes-Numerical-LES) methods are used to simulate flow through a labyrinth seal geometry and are contrasted with RANS solutions. Results show that LES and RANS-NLES is capable of accurately predicting flow behaviour of two seal flows with a scatter of less than 5 %. RANS solutions show the potential to perform poorly for the turbulence models tested. LES and hybrid RANS-NLES are found to be consistent and in agreement with measurements, providing a flexible numerical platform for design investigations. It also allows greater flow physics insights.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).