气固两相平面混合层的直接数值模拟

本文采用直接数值模拟(DNS)方法研究了气固两相混合层的涡量场和颗粒扩散,分析了湍流混合层中涡结构的卷起和配对过程,并讨论了大尺度涡结构的配对过程对平均速度、雷诺应力的影响。同时分析了混合层中不同Stokes数的颗粒在涡结构的作用下的混合和扩散。结果显示Stokes数为1的颗粒主要分布在流场大尺度涡结构的外边界,而Stokes数为0.01在涡结构的作用下,在流场中充分混合。     

高雷诺数气固湍流射流的直接数值模拟

本文对流动雷诺数 Re=5990的空间发展的气固两相湍流射流进行了直接数值模拟。其中对流场的求解采用具有四阶精度的紧致差分格式,对颗粒场的跟踪采用拉格朗日方法。结果表明,湍流拟序结构逐渐由对称模式发展到非对称模式;较小 Stokes 数的颗粒在流场中均匀分布,较大 Stokes 数的颗粒沿横向没有明显的扩散,而 Stokes 数为 1的量级的颗粒则大量聚集在大涡结构的外围。     

基于大涡模拟的小型转子发动机缸内湍流研究

本文应用大涡模拟的方法对小型转子发动机的缸内流场进行了三维数值分析,分别从拟序结构、湍流脉动和湍动能三个方面研究小型转子发动机缸内的瞬态湍流特征。计算结果表明,大涡模拟结合Q准则可以有效的识别小型转子发动机缸内流场的大尺度拟序结构。进气过程中高强度的涡团主要分布于燃烧室的前部和后部,压缩过程中高强度的涡团分布从燃烧室中部向后部转移。湍流在垂直于转子截面上的法向脉动最为剧烈,其中中部区域的剧烈程度远高于两侧。湍流在平行于前后端面的截面上的法向脉动最弱,没有明显的波动。从进气口打开到压缩上止点期间,缸内亚网格湍动能和平均气流速度出现两次明显的峰值,这分别与进气能量和大尺度涡团的破碎有直接的关系。     

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

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

两相平板混合层流动的双向耦合离散涡法数值模拟

提出了一种新的双向耦合离散涡方法,其能方便地将颗粒对流场的作用结合进离散涡模拟中,同时又不破坏离散涡方法固有的优点。采用这种方法对两相平板混合层流动进行了数值模拟。模拟结果表明,在颗粒作用存在的情况下,混合层流动仍能保持良好的自相似性,其动量厚度沿流向仍呈线性发展,说明颗粒不影响大涡结构间配对与合并的随机性。同时,颗粒对混合层流动高速端自由流有阻滞作用,并会使连续相湍流衰减。     

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

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

混合层流场中涡结构对流速度的特性

基于大涡模拟和光线追踪方法, 对光线穿越流场后的光程分布与混合层流场中涡结构之间的关系进行了分析, 提出了一种基于涡核位置提取的涡结构瞬时对流速度定量计算方法, 并使用直接几何测量数据进行了验证. 通过对不同尺寸的涡结构、涡-涡配对及融合过程中的涡结构和强压缩性流场中涡结构瞬时对流速度的定量数值计算, 揭示了混合层流场中涡结构对流速度的特性: 对单个涡结构而言, 其瞬时对流速度具有脉动特性, 且脉动幅度随涡结构尺寸和流场压缩性而变化; 在涡-涡配对及融合过程中, 涡对中各个涡结构的瞬时对流速度都表现出类似正弦波动的特点. 针对混合层流场中涡结构对流速度的特性, 给出了其背后的物理原因.     

空间发展湍流气粒两相平面混合层的大涡模拟与统计

本文对空间发展的湍流气固两相平面混合层流动进行了大涡模拟研究,其中气相亚网格尺度（SGS）使用结构函数模型,气相控制方程组采用SIMPLE方法求解,固体颗粒运动用拉格朗日方法计算。计算结果正确重现了流体涡结构的卷起、合并和破碎过程,以及小尺寸颗粒在涡边缘（低涡度区）的局部富集现象。对直径分别为42μm、72μm和135μm分别进行了模拟,并将统计结果和实验测量结果（Hishida　et　al[1]）比较,表明两者的平均速度吻合很好,但颗粒数密度和脉动速度存在较明显的差异,因此有必要对亚网格应力和颗粒之间的耦合作用以及拟序结构的三维性对颗粒运动的影响开展深入研究。     

气固两相三维圆柱绕流的直接数值模拟

本文采用直接数值模拟(DNS)方法研究了气周两相三维圆柱绕流的涡量场和颗粒扩散,并着重讨论了圆柱绕流中卡门涡街的形成和涡结构的转捩过程。同时分析了圆柱绕流中不同Stokes数的颗粒在涡结构作用下的横向扩散。结果显示Stokes数为1的颗粒主要分布在流场大尺度涡结构的外边界,而Stokes数为0．01的颗粒在涡结构的作用下,在流场中充分混合。     

气固两相混合层流场双向耦合的数值研究

采用双向耦合模型对有涡配对的二维气固两相混合展数值模拟,在考虑颗粒对流场反作用基础上进一步对颗粒间通过流体的相互作用进行分析。流场用拟谱方法求解,颗粒用颗粒轨道模型跟踪。结果发现,流场中大涡卷起和配对仍居主导地位;颗粒Ｓｔ数为Ｏ（０．１）～Ｏ（１）时,颗粒减弱了流场雷诺应力强度,加快涡量扩散; Ｓｔ数为Ｏ（１）时,颗粒分布极不均匀,主要集中在渴的边缘．     

A study on large coherent structures and noise emission in a turbulent round jet

A moderate Reynolds number,and high subsonic turbulent round jet is investigated by large eddy simulation.The detailed results(e.g.mean flow properties,turbulence intensities,etc.)are validated against the experimental data,and special attention is paid to study motions of coherent structures and their contributions to far-field noise.Eulerian methods(e.g.Q-criteria andλ2criteria)are utilized for visualizing coherent structures directly for instantaneous flow fields,and Lagrangian coherent structures accounting for integral effect are shown via calculating fields of finite time Lyapunov exponents based on bidimensional velocity fields.All visualizations demonstrate that intrusion of three-dimensional vortical structures into jet core occurs intermittently at the end of the potential core,resulting from the breakdown of helical vortex rings in the shear layer.Intermittencies in the shear layer and on the centerline are studied quantitatively,and distinctively different distributions of probability density function are observed.Moreover,the physical sound sources are obtained through a filtering operation of defined sources in Lighthill’s analogy,and their distributions verify that intrusion of vortical structures into the core region serves as important sound sources,in particular for noise at aft angles.The facts that intermittent behaviors are caused by motions of coherent structures and correlated with noise generation imply that to establish reasonable sound sources in active noise production region based on intermittent coherent structures is one of the key issues for far-field noise prediction.     

Direct numerical simulation of hypersonic transition induced by an isolated cylindrical roughness element

Hypersonic boundary layer transition induced by an isolated cylindrical roughness element is investigated using direct numerical simulation method based on a finite volume formulation. To simulate the transition procedure by resolving the generation and evolvement of small-scale coherent structures, and capture the shock wave at the same time, high-order minimum dispersion and controllable dissipation scheme is validated and then applied. The results are compared with the available measurements in the quiet wind tunnel, such as the dominated frequency and root mean square of pressure. The computational dominated frequency of 19.23 kHz is very close to the experimental one, 21 kHz. Also, the disturbances of the roughness are mostly generated by the “jet” just before the roughness, and then they travel and develop downstream with the shear layer and vortex shedding. The transition is mainly dominated by the instabilities of both the horseshoe vortex and the shear layer.     

The investigation of coherent structures in the wall region of a supersonic turbulent boundary layer based on DNS database

Through temporal mode direct numerical simulation, flow field database of a fully developed turbulent boundary layer on a flat plate with Mach number 4.5 and Reynolds number Reθ =1094 has been obtained. Commonly used detection meth- ods in experiments are applied to detecting coherent structures in the flow field, and it is found that coherent structures do exist in the wall region of a supersonic turbulent boundary layer. The detected results show that a low-speed streak is de- tected by using the Mu-level method, the rising parts of this streak are detected by using the second quadrant method, and the crossing regions from a low-speed streak to the high-speed one are detected by using the VITA method respectively. Notwithstanding that different regions are detected by different methods, they are all accompanied by quasi-stream-wise vortex structures.     

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.     

Effect of wall temperature on hypersonic turbulent boundary layer

Turbulent boundary layers at Mach 4.9 with the ratio of wall temperature to recovery temperature from 0.5 to 1.5 are investigated by means of direct numerical simulation. Various fundamental properties relevant to the influence of wall temperature on Morkovin’s scaling, standard and modified strong Reynolds analogies, and coherent vortical structures have been studied. It is identified that the scaling relations proposed for cool and adiabatic wall conditions, such as Morkovin’s scaling and the modified strong Reynolds analogy, are also applicable for hot wall condition. Moreover, the relation between the density and temperature fluctuations under the second-order approximation is derived and verified to provide a reliable prediction. Based on the analysis of coherent vortical structures, it is found that the orientation of vortex core can be quantitatively determined by means of the vector with its direction and modulus using the local strain direction and the imaginary part of the eigenvalue of velocity gradient tensor, respectively. As the increase of wall temperature, the spanwise distance between the two legs of hairpin vortex increases, and the mean swirling strength and the angle of vortical structure with respect to the wall plane also increase in the inner layer. The statistical properties relevant to vortical structures are nearly insensitive to the wall temperature in the outer layer.     

湍流射流与扩散火焰大涡拟序结构的波动特性研究

１前言湍流射流扩散燃烧方式提供了射流火焰与涡团相互作用的最基本形式,对研究在湍流射流剪切边界层内的反应物的卷吸混合、热量及动量的输运及湍流射流结构方面具有其特殊的意义。Ｋａｔｔａ［１］利用浮力与单步反应机理模型对Ｎ２－Ｈ２射流扩散火焰进行了直接数值模...     

高超声速椭圆锥边界层横流转捩特性大涡模拟

横流效应显著影响高超声速飞行器的三维边界层转捩过程, 深化对该流动机制的认识有助于提升和改善飞行器气动性能及热力学环境. 针对HIFiRE5椭圆锥绕流问题, 采用大涡模拟方法计算分析了超声速边界层横流转捩特性, 并揭示其中的流动机理. 参考HIFiRE5风洞模型试验条件, 数值模拟中椭圆锥来流入口处施加人工速度扰动以激发边界层内不稳定扰动波, 进而预测了高超声速边界层流动横流失稳、转捩过程等基本流动特征, 并基于转捩热流分布形态对比, 获得了与试验数据基本吻合的计算结果. 研究发现, 椭圆锥中心线流动汇聚形成的流向涡结构非常容易失稳, 另外在中心线及侧缘之间的中部区域存在较强的横流不稳定性, 两种机制共同作用影响边界层转捩过程. 此外, 分析了来流扰动幅值对边界层横流失稳转捩的影响, 并发现静来流条件下, 横流区域出现两组独立的定常横流涡结构, 而强噪声来流条件下, 中心线主涡和中部横流涡均发生失稳转捩, 且在椭圆锥表面形成多峰状的转捩阵面. 最后, 深入分析流场的压力脉动动力学特性, 揭示了三维边界层发生失稳转捩的非线性演化机制.      

空间发展平板混和层流动的大涡模拟-3D与2D模拟的比较

本文利用大涡模拟方法研究了空间发展平板混合层流动,分别在三维模型和二维模型下,比较了大尺度结构涡的演化过程,以及流向速度、流向脉动速度、横向脉动速度和Reynolds应力的统计时均结果,并与实验结果对照,指出3D模拟无论在物理真实性还是预报结果准确性方面都要明显优于2D模拟。     

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma = 6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at freestream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibility effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely tO OCCUr.     

Holographic flow visualization as a tool for studying three-dimensional coherent structures and instabilities

Holography is capable of three-dimensional (3D) representation of spatial objects such as fluid interfaces and particle ensembles. Based on this, we adapt it into a 3D flow visualization tool called Holographic Flow Visualization (HFV). This technique provides a novel means of studying spatially and temporally evolving complex fluid flow structures marked by a disperse phase or interfaces of different fluids. This paper demonstrates that HFV is a straightforward technique, especially when the In-line Recording Off-axis Viewing (IROV) configuration is used. The technique can be applied either as a stand-alone experimental tool for studying scalar-based coherent structures, flow instabilities, interactions of different fluids driven by fluid dynamics, interfacial phenomena, or as a precursor to volumetric 3D velocity vector field measurement of complex transient flow dynamics. Experimental results in several complex fluid flows and flames demonstrate the effectiveness of HFV. Different methods are used to mark flow structures undergoing different instabilities: 1) a vortex ring grown out of a drop of polymer suspension falling in water, 2) cascade of a bag-shaped drop of milk in water, and 3) internal flow structures of a jet diffusion flame.