混合层流动拟序结构的大涡模拟

采用大涡模拟方法对空间发展的二维平面混合层进行了数值模拟 ,动量方程采用分步投影法求解 ,亚格子项采用标准Smagorinsky亚格子模式模拟 ,压力泊松方程采用修正的循环消去法快速求解 ,同时求解了标志物输运方程以实现数值流场显示。模拟结果给出了混合层流动的瞬态发展过程以及流动中拟序结构的发展演变过程 ,成功地模拟了混合层发展中的各种瞬态细节过程 ,如涡的卷起、增长 ,涡与涡之间的配对、合并过程 ,以及大涡破碎为小涡的级联过程 ,为各种以混合层流动为原型流动的射流、尾流等工业流动的控制和优化提供了理论基础。     

次谐波对混合层中拟序涡合并过程影响的数值研究

本文用拟谱方法对混合层二涡、三涡的合并规律进行了数值研究,分析了基波和交谐波之间相位差对涡合并过程的影响,并对此进行了解释。计算结果表明：二涡合并有配地和撕裂两种形式,且涡的完全通过撕裂合并的现象仅在相位差为π/2时出现;三涡合并有三种形式。通过计算可见,通过施加多级次谐波并改变它们之间的相位差可以主动控制拟序涡的发展和演变。     

低密度流体界面不稳定性大涡模拟

采用拉伸涡亚格子尺度应力模型对湍流输运中的亚格子作用项进行模式化处理,发展了适用于可压多介质黏性流动和湍流的大涡模拟方法和代码MVFT（multi-viscous flow and turbulence）。利用MVFT代码对低密度流体界面不稳定性及其诱发的湍流混合问题进行了数值模拟。详细分析了扰动界面的发展,流场中冲击波的传播、相互作用、湍流混合区边界的演化规律,以及流场瞬时密度和湍动能的分布和发展。数值模拟获得的界面演化图像和流场中波系结构与实验结果吻合较好。三维和二维模拟结果的比较显示,两者得到的扰动界面位置、波系及湍流混合区边界基本一致,只是后期的界面构型有所不同,这也正说明湍流具有强三维效应。     

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

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

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

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

用格子Boltzmann模型模拟Lotka-Volterra系统

用格子Boltzmann方法模拟了Lotka-Voltcrra系统的动力学行为．通过使用多尺度技术,我们得到了扩散反应系统。以及零扩散的反应系统。从而得到常微分方程系统，这个理论结果表明该格子Boltzmann模型包含了Lotka-Volterra系统的非线性行为，数值结果表明附加的人为扰动可以改变某些初始分布。     

圆射流中拟序涡配对的数值模拟

用三维涡丝法结合涡丝的增加与合并技术,对不可压圆射流场,数值模拟了在扰动中加入次谐波后流场拟序涡的配对,配对情况与实验结果一致,通过对不同时刻涡量图的比较与分析,给出了次谐波初始强度对涡配对的影响,初始强度越大,涡配对出现越早,涡配对位置越往上游移动;同时还得到了基波与次谐波初始相位差和涡配对的关系,没有相位差时,次谐波增长最快,相位差增加时,涡配对的速度降低。说明控制次谐波的初始强度和基波与次谐波的相位差,能起到控制圆射流剪切层乃至控制整个射流场的目的。     

基于Vreman亚格子模型的有限开口空间内热驱动流数值模拟

利用基于Vreman亚格子模型的大涡模拟技术对有开口的单室和双室房间内热驱动流进行了数值模拟,利用函数分析法定量分析了模拟结果的准确性,并与Smagorinsky亚格子模型的模拟结果进行了比较.结果表明,Vreman和Smagorinsky亚格子模型的计算结果均能够满足工程的需求,但Vreman亚格子模型在开口附近区域的温度和U速度计算结果在整体上比Smagorinsky亚格子模型更接近实验值;Vreman亚格子模型未像Sma-gorinsky亚格子模型那样过高地估算壁面附近高温区域的粘性耗散;对于单室房间内热烟气层高度的预测,采用Vreman模型得到的计算结果准确性比Smagorinsky亚格子模型提高近50％.     

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

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

后台阶流动的数值模拟

访述了大涡模拟的基本思想，指出大涡模拟的效率主要取决于四个因素，即流动中须有大尺度涡存在、合理的计算格式、合适的滤波器和亚格子应力模型。在深入考虑粘性不可压缩流Navier—Stokes方程各个子项作用的基础上，提出二阶全展开Euler—Taylor—Galerkin有限元方法作为大涡模拟的离散格式，并采用Gauss滤波器，对典型算例——后台阶处的流动进行大涡模拟，计算结果与相关文献符合的很好。从计算结果还可以看出大涡模拟与二阶全展开ETG有限元方法的结合在捕捉涡系及反映涡动时变过程方面具有明显的优势，说明大涡模拟适合于边界几何形状复杂区域流动的模拟。同时应用二阶全展开ETG有限元方法对低雷诺数粘性不可压缩后台阶流动进行了计算，得到与相关文献符合良好的计算结果，即该方法也可独立用于对低雷诺数粘性不可压缩流动的计算。     

Vortex structure simulation for supersonic mixing layers using nonlinear PSE method

The method of nonlinear parabolized stability equations (PSE) is applied in the simulation of vortex structures in compressible mixing layer. The spatially-evolving unstable waves, which dominate the vortex structure, are investigated through spatial marching method. The instantaneous flow field is obtained by adding the harmonic waves to basic flow. The results show that T-S waves do not keep growing exponentially as the linear evolution, the energy transfer to high order harmonic modes, and that finally all harmonic modes get saturated due to nonlinear interaction. The mean flow distortion induced by the nonlinear interaction between the harmonic modes and their conjugate harmonic ones, makes great change of the average flow and increases the thickness of mixing layer. PSE methods can well capture the two- and three-dimensional large scale nonlinear vortex structures in mixing layers such as vortex roll-up, vortex pairing, and Λ vortex.     

超声速平面混合层小激波的形成与演变

为了揭示超声速混合层中小激波形成机理及其与涡相互作用的演变过程,本文基于大涡模拟（LES）方法,结合五阶精度混合TCD／WENO格式,对超声速平面混合层在对流马赫数为Mc=0．65条件下的流场结构进行了数值模拟,数值结果详细描述了超声速混合层中小激波的形成过程。研究了小激波形成后,随涡运动而产生的变形、脱落及发展过程。同时,对混合层双涡合并过程中,小激波与相邻涡相互作用所产生的变形与演变过程进行了讨论。     

LBE simulation of coherent vortex motion and heat transfer in jets of cross flow

Lattice Boltzmann equation method is used to simulate the coherent vortex motions and interactions and the heat transfer characteristics of jets in cross flow (JICFs) via TD2G9 model. After validation, the characteristics of cross flow under different Reynolds numbers are illustrated, including the mean profiles, the Reynolds stress tensor, the vortex and temperature fields, the temperature gradients near the walls, and the coherent correlation of vortex motions. The results show that the velocity profiles in JICF can be characterized by three basic regions, which are mainly caused by the mergence of JICF with the main flow. The temperature gradient near the walls can also be categorized by four basic regions, which are caused mainly by the impulse of JICFs too. Coherent vortex motions are found in JICF for Re = 3000, which are proved by strong periodic correlation of flow variables over a fixed area.     

A three-dimensional numerical study on instability of sinusoidal flame induced by multiple shock waves

The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier–Stokes(NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later.The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, η, that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time.The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter η can reasonably character the features of flame interface development induced by the multiple shock waves.     

Dynamics of large-scale eddies in turbulent flows

A model equation based on the equipartition of the turbulent dissipation is proposed for describing the dynamics of large-scale eddies in turbulent flows. The equation is reducible to the equation of motion of an inviscid fluid, so that the motion of the large-scale eddies can be described in terms of inviscid fluid dynamics. It is found that the large-scale eddies are always weakened by the background turbulence and their evolution is slowed down compared with the corresponding inviscid motion. In the case of turbulent mixing layer, its linear growth in downstream direction is accounted for by the exponential growth in time of the perturbation in an inviscid plane vortex sheet.     

An experimental study of the large scale coherent structures in a forced free shear layer

The dynamic characteristics of the large scale coherent structures in a forced free shear layer are experimentally studied by means of flow visualization. The quantitative measurements are acquired by the use of a LDV. It is shown that the development of the coherent structures can be greatly influenced by upstream artificial perturbations and as a result the mixing in the layer can be controlled. Like vortex merging, vortex splitting is also a common evolution pattern in the development of the coherent structures.     

数值研究平板方舵激波-湍流边界层干扰

数值研究了平板方舵激波-湍流边界层干扰流场。模拟出了分离激波与弓型激波砬撞后形成的“λ”激波结构;消晰地显示了分离区中的旋涡结构,发现流场中会出现二次分离涡,并从理论上分析了流场对称面涡心形态与非定常的关系,得到了涡心为不稳定螺旋点或出现极限环是非定常流动特征的新结论。     

Development of LBGK and incompressible LBGK‐based lattice Boltzmann flux solvers for simulation of incompressible flows

This paper presents lattice Boltzmann Bhatnagar–Gross–Krook (LBGK) model and incompressible LBGK model‐based lattice Boltzmann flux solvers (LBFS) for simulation of incompressible flows. LBFS applies the finite volume method to directly discretize the governing differential equations recovered by lattice Boltzmann equations. The fluxes of LBFS at each cell interface are evaluated by local reconstruction of lattice Boltzmann solution. Because LBFS is applied locally at each cell interface independently, it removes the major drawbacks of conventional lattice Boltzmann method such as lattice uniformity, coupling between mesh spacing, and time interval. With LBGK and incompressible LBGK models, LBFS are examined by simulating decaying vortex flow, polar cavity flow, plane Poiseuille flow, Womersley flow, and double shear flows. The obtained numerical results show that both the LBGK and incompressible LBGK‐based LBFS have the second order of accuracy and high computational efficiency on nonuniform grids. Furthermore, LBFS with both LBGK models are also stable for the double shear flows at a high Reynolds number of 10⁵. However, for the pressure‐driven plane Poiseuille flow, when the pressure gradient is increased, the relative error associated with LBGK model grows faster than that associated with incompressible LBGK model. It seems that the incompressible LBGK‐based LBFS is more suitable for simulating incompressible flows with large pressure gradients. Copyright © 2014 John Wiley & Sons, Ltd.     

Generation of the second,fourth, eighth,and subsequent harmonics by a quadratic nonlinear hyperelastic longitudinal plane wave

The perturbation (small-parameter) method is used to obtain the first three approximations for the problem of a harmonic longitudinal plane wave propagating in a quadratic nonlinear hyperelastic material described by the classical Murnaghan model. The subsequent approximations are discussed. The contribution of each approximation to the overall wave pattern is analyzed. It is shown that the third approximation corrects the prediction of the evolution of the initial wave profile. Which of the harmonics dominates depends on the distance traveled by the wave: the second harmonic is generated first, then it transforms into the fourth harmonic, and finally, as the distance increases, the eighth harmonic shows     

槽道内涡波流场的POD分析

对槽道内涡波流场的瞬态速度矢量场进行了2DPIV测量实验，将2DPIV测量的矢量场数据进行POD分析，根据POD分解的各阶模态的能量比确定了表征涡波流场主导结构的前15阶模态。结果表明，POD分解的前15阶模态发现槽道内涡波流场是由槽道壁面剪切层诱导的涡列以及伴随的波状主流组成；流场中大尺度的涡旋发展为涡对，对波状主流的脉动频率产生影响；根据涡波流场中的驻点和鞍点，获取了流场的大尺度涡对、平均流场以及Helmholtz涡环等明显特征；最后根据POD分解的前15阶模态对槽道内涡波流场进行重组，重组流场表征了槽道内层流状态下波状主流的形态和涡旋共存的涡波结构以及驻点和鞍点的位置处涡旋的变化等主要特征，有效地剔除了PIV测量流场中的随机信息，保留了PIV测量流场的主导特征。