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

在可压缩多介质粘性流体动力学高精度计算方法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方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

基于全局动态Vreman模型的复杂非定常湍流的大涡模拟

在湍流数值模拟方法中,大涡模拟方法可以提供丰富的大涡旋信息,已逐渐成为复杂湍流问题数值研究的重要方法。而大涡模拟中,最重要的一环是尽量准确地构建能反映流场物理本质特征的亚格子应力模型。基于该思想,将一种新型的大涡模拟亚格子应力模型-Vreman亚格子应力模型用于高雷诺数三维后台阶流动的求解,计算结果与实验结果进行对比分析结果较吻合,验证了该模型的可靠性。这是对该模型用于无任何均匀流动方向的高雷诺数复杂湍流非定常流动的首次检验,计算结果优于基于传统的Smagorinsky涡粘性的动态亚格子模型。     

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

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

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组，采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscous-flow and turbulence)，模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF₆界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明，初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前，非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大；反射激波加载后，非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步，定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外，还与平面激波诱导的RM不稳定性进行了简单对比，发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用，使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

流对波边界层湍流运动的影响

采用大涡模拟方法和Smagorinsky亚格子模型，求解三维Navier-Stokes方程，研究了波流边界层中的湍流特性．将大涡模拟结果与相应的直接数值模拟结果和实验数据进行比较，吻合较好．获得了不同波雷诺数，不同波流比情况下的大涡模拟数据库，并由此分析了波流边界层中各种湍流统计量，如速度廓线、剪应力、湍流强度等的变化规律．     

大涡模拟及其在湍流燃烧中的应用

大涡模拟作为一种研究湍流流动和湍流燃烧的有效手段,在国际上已经得到广泛应用。本文在回顾了大涡模拟(LES)的基本思想及其实施方法的基础上着重介绍了前人在大涡模拟的亚格子湍流模式和亚格子燃烧模式中的研究成果,同时给出了采用不同亚格子模式的大涡模拟在湍流燃烧中的应用实例,指出了大涡模拟在湍流燃烧中的重要作用,为大涡模拟的进一步发展和应用提供参考。      

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

速度估计模型在大涡模拟中的应用

介绍了速度估计模型的基本思想及其在物理空间的实现。速度估计模型依靠湍流大尺度的非线性作用估计出小尺度,从而可直接求解亚格子应力项而不需要额外的模型。本文采用该模型对不同雷诺数下的各向同性衰减湍流进行了模拟,并与直接数值模拟、理论分析和其他亚格子模型的大涡模拟结果进行了比较。也初步考察了网格分辨率、不同精度的紧致格式对大涡模拟结果的影响。     

基于人工神经网络的湍流大涡模拟方法

大涡模拟方法(LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中相对误差较大, 在后验验证中耗散过强. 近几年来,机器学习方法在湍流建模问题中得到了越来越多的应用.本文介绍了基于人工神经网络(ANN)的湍流亚格子模型的最新进展.详细地讨论了人工神经网络混合模型、空间人工神经网络模型和反卷积人工神经网络模型的构造方法.借助于人工神经网络强大的数据插值能力,新的亚格子模型的先验精度和后验精度均有显著提升. 在先验验证中,新模型所预测的亚格子应力的相关系数超过了0.99,在预测精度上远高于传统的大涡模拟模型. 在后验验证中,新模型对各类湍流统计量和瞬态流动结构的预测都优于隐式大涡模拟方法、动态Smagorinsky模型、动态混合模型等传统模型.因此, 人工神经网络方法在发展复杂湍流的先进大涡模拟模型中具有很大的潜力.      

基于人工神经网络的湍流大涡模拟方法

大涡模拟方法 (LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中相对误差较大,在后验验证中耗散过强.近几年来,机器学习方法在湍流建模问题中得到了越来越多的应用.本文介绍了基于人工神经网络(ANN)的湍流亚格子模型的最新进展.详细地讨论了人工神经网络混合模型、空间人工神经网络模型和反卷积人工神经网络模型的构造方法.借助于人工神经网络强大的数据插值能力,新的亚格子模型的先验精度和后验精度均有显著提升.在先验验证中,新模型所预测的亚格子应力的相关系数超过了0.99,在预测精度上远高于传统的大涡模拟模型.在后验验证中,新模型对各类湍流统计量和瞬态流动结构的预测都优于隐式大涡模拟方法、动态Smagorinsky模型、动态混合模型等传统模型.因此,人工神经网络方法在发展复杂湍流的先进大涡模拟模型中具有很大的潜力.     

Experimental and numerical investigation of inclined air/SF6 interface instability under shock wave

The shock tube experiments of inclined air/SF6 interface instability under the shock wave with the Mach numbers 1.23 and 1.41 are conducted. The numerical simulation is done with the parallel algorithm and the multi-viscous-fluid and turbulence (MVFT) code of the large-eddy simulation (LES). The developing process of the interface accelerated by the shock wave is reproduced by the simulations. The complex wave structures, e.g., the propagation, refraction, and reflection of the shock wave, are clearly revealed in the flows. The simulated evolving images of the interface are consistent with the experimental ones. The simulated width of the turbulent mixing zone (TMZ) and the displacements of the bubble and the spike also agree well with the experimental data. Also, the reliability and effectiveness of the MVFT in simulating the problem of interface instability are validated. The more energies are injected into the TMZ when the shock wave has a larger Mach number. Therefore, the perturbed interface develops faster.     

Shock-induced mixing zone characterization: an attempt by hot-wire diagnostic

An attempt to extract velocity and molar fraction from a single hot-wire trace within a turbulent mixing zone induced by a shock accelerated gaseous interface has been proposed. Experiments have been conducted for negative and positive density jumps across the interface. The hot-wire signals clearly show interfaces between mixed and unmixed regions and the locations of incident and reflected shocks. With some hypotheses on the temperature, velocity and molar fraction profiles within the turbulent mixing zone have been obtained solving an inverse problem. Results show that if the molar fraction profiles follow physically coherent evolutions, those of the local velocity are strongly correlated with the choice of its variation range. So, we reasonably think that the results obtained from single wire have to remain limited to interface and shock locations. And it is only by coupling the present technique with the laser Doppler velocimetry, which we will be able to possibly obtain reliable estimates of the variations of quantities in the turbulent mixing zone.      

Numerical modeling of interaction between a turbulent mixing zone and a local perturbation of the density field in a pycnocline

A numerical model is developed, and interaction between a turbulent mixing zone and a local perturbation of the density field in a pycnocline is studied. It is demonstrated that the flow induced by a local perturbation of the density field can generate turbulence energy in the turbulent mixing zone and increase its lifetime.     

Numerical simulations of the process of multiple shock–flame interactions

Based on a weighted essentially nonoscillatory scheme, the multiple interactions of a flame interface with an incident shock wave and its reshock waves are numerically simulated by solving the compressible reactive Navier–Stokes equations with a single-step Arrhenius chemical reaction. The two-dimensional sinusoidally perturbed premixed flames with different initial perturbed amplitudes are used to investigate the effect of the initial perturbation on the flame evolutions. The results show that the development of the flame interface is directly affected by the initial perturbed amplitudes before the passages of reshock waves, and the perturbation development is mainly controlled by the Richtmyer–Meshkov instability(RMI). After the successive impacts of multiple reshock waves, the chemical reaction accelerates the consumption of reactants and leads to a gradual disappearance of the initial perturbed information. The perturbation developments in frozen flows with the same initial interface as those in reactive flows are also demonstrated.Comparisons of results between the reactive and frozen flows show that a chemical reaction changes the perturbation pattern of the flame interface by decreasing the density gradient,thereby weakening the baroclinic torque in the flame mixing region, and therefore plays a dominant role after the passage of reshock waves.     

Detonation interaction with a diffuse interface and subsequent chemical reaction

We have investigated the interaction of a detonation with an interface separating a combustible from an oxidizing mixture. The ethylene-oxygen combustible mixture had a fuel-rich composition to promote secondary combustion with the oxidizer in the turbulent mixing zone that resulted from the interaction. Diffuse interfaces were created by the formation of a gravity current using a sliding valve that initially separated the test gas and combustible mixture. Opening the valve allowed a gravity current to develop before the detonation was initiated. By varying the delay between opening the valve and initiating the detonation it was possible to achieve a wide range of interface conditions. The interface orientation and thickness with respect to the detonation wave have a profound effect on the outcome of the interaction. Diffuse interfaces result in curved detonation waves with a transmitted shock and following turbulent mixing zone. The impulse was measured to quantify the degree of secondary combustion, which accounted for 1–5% of the total impulse. A model was developed that estimated the volume expansion of a fluid element due to combustion in the turbulent mixing zone and predicted the resulting impulse increment.      

Towards large eddy simulation of isothermal two-phase flows: Governing equations and a priori tests

This article reports on the potential of application of LES in the calculation of turbulent two-phase flows, in the case where each phase is resolved and interfaces remain much larger than the mesh size. In comparison with single-phase flow, successful application of LES to resolve two-phase flow problems should account for the complex interaction between turbulence and interfaces. Non-linear transfers of turbulent energy across the interface have to be accurately modeled. The derivation of the complete filtered two-phase flow governing equations has been formulated to deal with turbulence at the interface in a comprehensive and practical way. Explicit filtering of 2D direct numerical simulations has been employed to evaluate the order of magnitude of the new subgrid contributions. A parametric study on the academic test case of two counter-rotative vortices and a more complex test case of phase inversion in a closed box have been utilized to perform an order of magnitude analysis of different transport mechanisms. Important features of turbulent energy transfer across the interface have been discussed. Analyses of the numerical results have been conducted to derive conclusions on the relative importance of the different subgrid scale contributions, and modeling issues and solutions are provided.     

Numerical simulation of a developed shear turbulence

Results of two-dimensional numerical studies of turbulence that arises at the interface of two flows of poorly compressible gases are described. The results were obtained using a MAKh software system. The interrelation between spatial and time problems on the development of a turbulent zone induced by shear instability is analyzed. A constant that characterizes the degree of turbulent shear mixing is calculated. The effect of the density difference of the mixing fluids on the growth rate of the turbulence zone is studied. For all density differences considered, the coefficient of heterogeneity of the resultant mixture is evaluated. Institute for Technical Physics, Snezhinsk 456770. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 1, pp. 77–83, January–February, 2000.     

青藏铁路冻土区环境问题对工程安全可靠性影响

当大涡模拟用于研究化学反应流动时，传统的滤波方法会导致化学反应项不封闭. 为克服这 个困难，发展了条件滤波大涡模拟方法. 在选择适当的条件变量后，条件滤波的化学反 应项可以表达为一个封闭项. 但同时也带来了新的问题：条件滤波耗散或条件滤波扩散项的 不封闭. 为解决这一问题，采用了直接数值模拟方法研究了它们在大小尺度上的统计特 性. 研究结果表明：条件滤波耗散和扩散对于大尺度的依赖主要体现在大尺度标量场中扩散 层结构的影响，同时小尺度脉动的变化几乎与条件滤波扩散无关，而它对条件滤波耗散却显 现出明显的作用. 在构造条件滤波耗散的亚格子模型时，小尺度脉动的作用不容忽视.