用基于M-SST模型的DES数值模拟喷流流场

脱体涡数值模拟方法(dettached eddy simulation,DES)是把雷诺平均Navier-Stokes方程(RANS)方法及大涡模拟方法(LES)结合起来模拟有脱体涡的湍流流场的数值模拟方法,其主要思想是在物面附近解雷诺平均Navier-Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文在剪切应力传输(SST)湍流模型的基础上用DES及混合非结构网格数值模拟具有横向喷流的湍流流场,算法采用Osher逆风格式,利用该套程序(包括网格生成及算法),对导弹在不同马赫数下的喷流流场进行了数值模拟,并与同时开展的实验研究的结果进行了对比,结果表明用该方法处理这类问题是较准确的。     

基于卡门尺度和滤波方法的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方法.     

超声速来流横向狭缝喷流标量输运的混合RANS/LES模拟

采用基于时均湍流模型的混合RANS/LES模拟方法对超声速来流狭缝横向喷流标量输运过程进行大涡模拟,控制方程对流项采用五阶精度的WENO格式求解,过滤后的组分方程中的亚格子组分对流通量采用梯度扩散模化.模拟得到了超声速来流狭缝横向喷流大尺度涡和组分的演化过程,研究表明喷流具有明显的非定常周期性特点,并且组分随时间的空间分布特征决定于流动的大尺度结构.模拟得到的压力的统计时均结果和实验结果定性一致,但喷流前后的分离区较实验值过大,这表明混合RANS/LES需要进一步改进.     

高速流场气动光学RANS/DSMC混合数值模拟算法研究

提出了适用于高速流场气动光学数值模拟的RANS/DSMC混合算法.通过RANS对全局时均流场进行数值模拟,再对其中局部流场的脉动量采用DSMC进行数值模拟,以Maxwellian速度分布实现RANS宏观量信息向DSMC微观量信息的传递.采用超声速环境下尖劈模型对该混合算法进行校验,对比试验结果论证了算法的正确性.在计算...     

高阶矩湍流模型研究进展及挑战

高阶矩模型是湍流模式理论研究中的难点和前沿. 自周培源先生首次建立一般湍流的雷诺应力输运方程起, 为了更精确的预测复杂流动, 人们从未间断过对高阶矩模型的研究. 尤其进入新世纪以来, 随着计算机硬件水平的飞跃和高精度数值算法的突破, 湍流模拟方法正由RANS向LES转变. 而无论对于RANS框架、LES框架还是两者混合, 高阶矩模式都是其中先进封闭模式的代表. 基于此, 本文对高阶矩模型的发展情况进行了总结, 重点包括高阶矩模型中各项的建模方式、尺度提供方程的演化历程和数值求解技术的关键需求. 然后, 通过几类典型湍流问题展示了其相对于传统涡黏模型的优势, 并且给出了部分CFD软件对高阶矩模型的集成情况. 最后对高阶矩湍流模型未来面临的挑战和发展的方向进行了展望.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

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

求解雷诺平均(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模型, 修正模型对后验速度场、下壁面平均压力系数和摩擦力系数的预测精度均有较大提升. 可以发现对雷诺应力线性部分的隐式处理可以增强数值求解的稳定性, 对雷诺应力非线性部分的修正可以提升模型对流场各向异性特征预测的性能, 并且多次修正后的模型表现出更高的预测精度. 因此, 该算法在数据驱动湍流建模和工程应用中具有很大的应用潜力.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全.针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见.开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据.采用直接数值模拟方法对来流马赫数6.0,34?压缩拐角内激波与湍流边界层的干扰问题进行了研究.基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性.通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制.采用动态模态分解方法,探讨了干扰流场的非定常运动历程.研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性.干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力-膨胀项,而对膨胀-耗散项影响较小.高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.     

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程，对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法，亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动，数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明，数值模拟结果 与理论吻合，得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致，湍流流谱定性合理.     

近壁湍流统计性质研究

采用大涡模拟方法,模拟了槽道湍流,得到了不同雷诺数下槽道湍流的结果. 在此基础上,研究了平均速度、雷诺应力、脉动动能和脉动速度均方根的分布;讨论了平均速度的壁面律问题;给出了雷诺应力、脉动动能和脉动速度均方根随雷诺数的变化规律,其中雷诺应力、脉动动能给出了定量公式.      

Simulation of the flow past a circular cylinder in the supercritical regime by blending RANS and variational-multiscale LES models

A strategy which blends a variational multiscale large eddy simulation (VMS-LES) model and a RANS model in a hybrid approach is investigated. A smooth blending function, which is based on the value of a blending parameter, is used for switching from VMS-LES to RANS. Different definitions of the blending parameter are investigated. The capabilities of the novel hybrid approach are appraised in the simulation of the flow around a circular cylinder at a Reynolds number 1.4×10⁵, based on the freestream velocity and on the cylinder diameter, in the presence of turbulent boundary-layer due to turbulent inflow conditions. A second study at Reynolds numbers from Re=6.7×10⁵ to 1.25×10⁶ is also presented. The effect of using the VMS-LES approach in the hybrid model is evaluated. Results are compared to those of other RANS, LES and hybrid simulations in the literature and with experimental data     

A hybrid RANS/LES approach with scale-adaptive capabilities for highly separated flows

A new hybrid RANS/LES approach with scale-adaptive capabilities is developed. The blending function in the SST model is adopted to prevent the invasion of the von Karman length scale to the RANS region, and the compressibility correction proposed by Wilcox is incorporated to produce a realistic shear layer development in compressible flows. The new model is validated for a subcritical flow past a circular cylinder and a supersonic base flow. Time-averaged turbulent statistics predicted by the new model show fairly good agreement with the experimental data, slight improvements over DES simulations, and are much better than SAS results. The main advantage of the new model over the DES method is that the distribution of the blending function reflects local vortex structures instead of grid spacing in the turbulent wake. The sequence of the effect intensity of the compressibility correction from strong to weak is SAS, the new model and DES.     

Hybrid LES-RANS using synthesized turbulent fluctuations for forcing in the interface region

The main bottleneck in using Large Eddy Simulations at high Reynolds number is the requirement of very fine meshes near walls. One of the main reasons why hybrid LES-RANS was invented was to eliminate this limitation. In this method unsteady RANS (URANS) is used near walls and LES is used away from walls. The present paper evaluates a method for improving standard LES-RANS. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane between the LES and URANS regions in order to trigger the equations to resolve turbulence. The turbulent fluctuations are taken from synthesized homogeneous turbulence assuming a modified von Kármán spectrum. Both isotropic and non-isotropic fluctuations are evaluated. The new approach is applied to fully developed channel flow and it is shown that the imposed fluctuations considerably improve the predictions. It is found that increasing the prescribed turbulent length scale of the synthesized turbulence provides excellent agreement with the classical log-law.     

Simulation of Transverse Combustion Instability in a Multi-Injector Combustor using the Time-Domain Impedance Boundary Conditions

The time-domain impedance boundary condition (TDIBC) is used as a reduced-order model (ROM) in large-eddy simulation (LES) to study self-sustained transverse oscillations in an experimentally studied high-pressure, shear coaxial multi-injector combustor. This work is an extension of the recent study using ROM-LES to simulate a single-element combustor that exhibited longitudinal instability. Here, we focus on transverse instability in a seven-injector combustor. The fuel and oxidizer inlets are truncated and the conventional inflow boundary conditions at the original inlet are replaced by an impedance describing function (IDF) in the form of a reflection coefficient that couples with LES through characteristic based boundary conditions at the truncated inlet. The impedance model is also generalized to include the effects of entropy fluctuations at the inflow. The hybrid ROM-LES simulations are compared with LES simulations with the full combustor geometry. Results show very good agreement and confirm that the use of TDIBC within LES is a viable tool to account for complex acoustic/boundary interaction in a physical way without explicitly solving the full geometry at LES level. Some simplifications and approximations have to be invoked and these constraints are also discussed.     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Hybrid LES-RANS: An approach to make LES applicable at high Reynolds number

The main bottle neck for using large eddy simulations (LES) at high Reynolds number is the requirement of very fine meshes near walls. Hybrid LES-Reynolds-averaged Navier-Stokes (RANS) was invented to get rid of this limitation. In this method, unsteady RANS (URANS) is used near walls and away from walls LES is used. The matching between URANS and LES takes place in the inner log-region. In the present paper, a method to improve standard LES-RANS is evaluated. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane in order to provide the equations in the LES region with relevant turbulent structures. The fluctuations are taken from a DNS of a generic boundary layer. Simulations of fully developed channel flow and plane asymmetric diffuser flow are presented. Hybrid LES-RANS is used both with and without forcing conditions.     

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.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

Hybrid Unsteady RANS and PDF Method for Turbulent Non-Reactive and Reactive Flows

A hybrid unsteady Reynolds-averaged numerical simulation (U-RANS) and probability density function (PDF) method is developed for turbulent non-reactive and reactive flows. The resulting modeled equations are solved by a consistent hybrid finite volume and Lagrangian Monte-Carlo particle method. Both turbulent non-reactive and reactive flows in a rectangular channel containing a triangular-shaped bluff-body are simulated. One-step and two-step mechanisms for propane/air combustion are used for the reactive case. The time-averaged results are compared with both experimental data and numerical results from the literature using large eddy simulation (LES) and steady RANS. The results of the present method are in good agreement with the experimental data, and they improve the numerical results available in the literature.     

Simulation of supersonic turbulent flow in the vicinity of an inclined backward-facing step

Large eddy simulation (LES) is a viable and powerful tool to analyse unsteady three-dimensional turbulent flows. In this article, the method of LES is used to compute a plane turbulent supersonic boundary layer subjected to different pressure gradients. The pressure gradients are generated by allowing the flow to pass in the vicinity of an expansion–compression ramp (inclined backward-facing step with leeward-face angle of 25°) for an upstream Mach number of 2.9. The inflow boundary condition is the main problem for all turbulent wall-bounded flows. An approach to solve this problem is to extract instantaneous velocities, temperature and density data from an auxiliary simulation (inflow generator). To generate an appropriate realistic inflow condition to the inflow generator itself the rescaling technique for compressible flows is used. In this method, Morkovin's hypothesis, in which the total temperature fluctuations are neglected compared with the static temperature fluctuations, is applied to rescale and generate the temperature profile at inlet. This technique was successfully developed and applied by the present author for an LES of subsonic three-dimensional boundary layer of a smooth curved ramp. The present LES results are compared with the available experimental data as well as numerical data. The positive impact of the rescaling formulation of the temperature is proven by the convincing agreement of the obtained results with the experimental data compared with published numerical work and sheds light on the quality of the developed compressible inflow generator.