后台阶分离流动中大涡结构演变的数值模拟

本文对后台阶分离流动中涡结构的演变进行了大涡模拟，研究了流场结构的变化规律。详细讨论了随着雷诺数的增加流场结构的典型特征的变化规律，指出流场中的涡结构随着雷诺数的增大变得十分复杂和丰富，回流区的数目、大小及其出现的位置也显著地不同。这些结果与已有的一些实验值和流场显示结果是吻合的。在此基础上，进一步研究了高雷诺数时流场中大尺度涡结构的瞬时发展和演化过程，展示了其中大涡的产生、追随、吸引、合并和破碎等过程。对于高雷诺数情况，对大涡模拟得到的数值结果进行了统计，得到的时均速度分布以及台阶后方的回流区长度与现有的其他实验结果符合得很好。本研究是针对后台阶分离流动深入开展湍流控制以及两相流动研究的基础。     

可压缩气固混合层中离散相与连续相的相互作用研究

尽管已有许多文献采用数值模拟方法研究两相流问题,但主要是集中不可压流动方面.本文采用Eul-er-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流动进行了研究.气相流场采用非定常全Navier-Stokes方程描述,并应用具有空间三阶精度的WNND(Weighted Non-Oscillatory, Contai-ning No Free Parameters and Dissipative)格式进行数值高散.固相方程采用二阶单边三点差分离散.在考虑流场对固粒作用的同时,也计及颗粒对流场的反作用.主要研究混合层大尺度涡对颗粒扩散特性的影响及颗粒对流场结构的影响问题.在对流马赫数为0.5时,研究不同Stokes数颗粒在连续流场中的扩散特性,而在对流马赫数为0.8时研究了不同Stokes数颗粒对流场小激波结构的影响.     

稀疏两相湍流的惯性颗粒倾向性弥散

气相采用大涡模拟、颗粒相采用拉格朗日轨道模型的方法对后台阶突扩流、充分发展槽道流和圆湍射流3种典型的稀疏气固两相流动进行了数值模拟,研究了颗粒倾向性弥散的特征和规律. 研究表明颗粒的跟随性和倾向性相联系,颗粒惯性和大涡结构同时决定颗粒的倾向性分布特征. Stokes数量级为1(气相时间参考尺度取为宏观特征时间尺度)左右的颗粒,倾向性分布特征最强烈. 颗粒倾向分布于低涡量(或是低脉动速度)的湍流区域.      

可压缩燃烧反应转捩混合层直接数值模拟

针对三维时间发展可压缩氢/氧非预混燃烧反应平面自由剪切混合层,采用5阶迎风/6阶对称紧致混合差分格式以及3阶显式Runge-Kutta时间推进方法,直接数值模拟了伴随燃烧产物生成和反应能量释放, 流动受扰动激发失稳并转捩的演化过程. 在转捩初期, 获得了${\it\Lambda}$涡、马蹄涡等典型的大尺度拟序结构,观察到了流动失稳后发生双马蹄涡三维对并的现象, 大尺度结构呈较好的对称性.在流动演化后期, 大尺度结构逐次破碎形成小尺度结构, 混合层进入转捩末期,呈明显的不对称性.      

颗粒物在人体上呼吸道模型中沉降规律的数值研究

本文采用RNG k-ε湍流模型对简化的上呼吸道内的流场进行数值仿真,并结合随机涡相互作用模型对上呼吸道内颗粒的局部运动沉积特性进行了数值模拟研究.结果表明,惯性冲击是微米级颗粒物主要的沉降机理,惯性参量能很好地衡量颗粒的沉降规律.口腔中流动较平稳,颗粒受扰动较小.颗粒由于惯性冲击而运动靠近咽部后部并易于沉积在咽部后壁面...     

风口区空气流动的有限元模拟与实验研究

通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.     

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

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

斑马鱼直线加速游动试验研究

海洋生物低噪音、高速、髙效游动能力是任何人造水下航行器所无法比拟的. 借助时间解析粒子图像测速技术对斑马鱼直线加速游动过程进行精细流场测量, 对其运动学行为特性和动力学机理进行分析. 同时应用双正交分解对涡量场进行模态分解, 获取流场的时间演化和空间分布特征. 从流动机理的角度探究斑马鱼游动过程的流动结构特征及旋涡动态演化特性. 试验结果表明: 流动可视化展现了整体涡流尾迹的结构分布,方便探究运动特性与旋涡尾迹之间的耦合关系. 斑马鱼从运动开始时体干保持着鲹科式的运动规律, 游动时的动能主要由前几次大幅的摆尾过程提供, 后续的摆尾主要调整方向及姿态. 两次不同方向的摆尾动作会形成一对方向相反的旋涡, 并在时序下旋涡逐渐脱落. 同时尾流的涡量变化在一定程度上反映鱼体的游向的变化. 基于双正交分解分解后的时间演化结果验证本次试验在时间上涡量场具有合理的恒定幅度, 空间分布表明低阶空间模态表征斑马鱼游动的主要涡流动结构, 高阶空间模态表征涡流动的细节结构. 研究鱼类游动时的摆尾推进机制与动力学特性能够为高效率的仿鱼类推进装置设计提供一定科学参考.      

二维后向台阶流流动特性的实验研究

利用PIV技术,通过系统对150相似文献   

格栅-空腔自激振荡流场扰动形成机理研究

格栅-空腔流动会引发流场自激振荡现象,产生结构振动及噪声问题.研究此类现象的形成机理对相关设备的减振降噪设计具有重要意义.目前已知这种现象是在流场正逆向两种扰动的耦合作用下形成的,但上述两种扰动的形成机理尚未明确.针对该问题,建立了格栅-空腔流动数值模型,并进行非稳态数值模拟.通过分析压力振荡数据及流场结构演化过程,对两种扰动的形成过程进行了研究.研究结果表明,正逆向两种扰动的形成均与格栅间隔中形成的小尺度涡团有关.上游小尺度涡团从格栅间隔中依次脱离并聚集演化为大尺度涡团,形成该流动系统的正向扰动.下游小尺度涡团脱离格栅间隔激发的逆向扰流构成该流动系统的逆向扰动,并对正向扰动产生持续的正反馈作用.     

混合层中粒子扩散的实验研究

本文应用流动显示和LDA测量,对混合层中粒子的扩散问题进行了实验研究。结果表明：粒子的扩散强烈地依赖于St数,具有较小St的粒子的扩散近似地与流体示踪粒子相同,随着St的增大,扩散角度增大,也就是说粒子扩散加快,并且随着粒子惯性的增大,它对湍流脉动的响应将会减弱。     

Dispersion of fine settling particles from an elevated source in an oscillatory turbulent flow

The present paper examines the stream-wise dispersion of suspended fine particles with settling velocities in an oscillatory turbulent shear flow with or without a non-zero mean over a rough-bed surface when the particles are being released from an elevated continuous source. A finite-difference implicit method is employed to solve the unsteady turbulent convective-diffusion equation. A combined scheme of central and four-point upwind differences is used to solve the steady state equation and the Alternating Direction Implicit (ADI) method is adopted for unsteady equation. It is shown how the mixing of settling particles is influenced by the tidal oscillatory current and the corresponding eddy diffusivity when the initial distribution of concentration regarded as a line-source. The vertical concentration profiles of suspended fine particles with settling velocities are presented for different downstream stations for various values of settling velocity and the frequency of the oscillation in tidal flow. For two-dimensional unsteady dispersion equation, the behaviour of iso-concentration lines for different values of settling velocity, frequency of the oscillation, dispersion time and releasing height is studied in terms of the relative importance of convection and eddy diffusion.     

An LES Turbulent Inflow Generator using A Recycling and Rescaling Method

The present paper describes a recycling and rescaling method for generating turbulent inflow conditions for Large Eddy Simulation. The method is first validated by simulating a turbulent boundary layer and a turbulent mixing layer. It is demonstrated that, with input specification of mean velocities and turbulence rms levels (normal stresses) only, it can produce realistic and self-consistent turbulence structures. Comparison of shear stress and integral length scale indicates the success of the method in generating turbulent 1-point and 2-point correlations not specified in the input data. With the turbulent inlet conditions generated by this method, the growth rate of the turbulent boundary/mixing layer is properly predicted. Furthermore, the method can be used for the more complex inlet boundary flow types commonly found in industrial applications, which is demonstrated by generating non-equilibrium turbulent inflow and spanwise inhomogeneous inflow. As a final illustration of the benefits brought by this approach, a droplet-laden mixing layer is simulated. The dispersion of droplets in the near-field immediately downstream of the splitter plate trailing edge where the turbulent mixing layer begins is accurately reproduced due to the realistic turbulent structures captured by the recycling/rescaling method.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Large Eddy Simulations of the Flow in the Near-Field Region of a Turbulent Buoyant Helium Plume

Large eddy simulations are conducted in the near-field region of a large turbulent buoyant helium plume. The CFD package FireFOAM is applied to that purpose. The transient and mean flow dynamics are discussed as a function of grid resolution, with and without the use of the standard Smagorinsky subgrid scale (SGS) model. Small scale structures, formed at the edge of the plume inlet due to baroclinic and gravitational mechanisms and subject to flow instabilities, interact with large scale features of the flow, resulting in a puffing cycle. In general, the LES calculations reproduce the main features of the turbulent plume, with better agreement when the Smagorinsky type SGS model is applied. In particular, the puffing cycle is recovered in the simulations with correct frequency. The mean and rms values of the velocity components are well predicted with use of the SGS model, even on relatively coarse meshes. Agreement for the species mass fraction (mean and rms values) is less satisfactory, but in line with results found in the literature.     

Flow modulation of a planar free shear layer with large bubbles—direct numerical simulations

The flow of a planar free shear layer with cylindrical bubbles is simulated using a finite difference/front tracking scheme. This approach allows direct numerical simulation of the multiphase flow by wholly incorporating the local bubble flow field in conjunction with the large scale vortical structures of the liquid. The role of large bubbles in modifying low Reynolds number ( 250) shear flow structures is investigated, specifically for bubbles whose diameter approaches the scale of the largest liquid eddies. The results indicate that duration of eddy crossing is the main mechanism for flow modulation, which is typically characterized by decreased vortex coherency and size, modified fluctuation statistics and significant variations in pairing/merging phenomena. The comparison of fluctuating statistics and flow field visualization also allowed qualitative discrimination between the modulation of the non-linear eddy dynamics and fluctuations due simply to the random bubble induced perturbations.     

壁面加热湍流相干结构的子波分析实验研究

通过在平板湍流边界层沿流向固壁表面平行放置若干条通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的大尺度流向涡结构,改变了平板湍流边界层中不同尺度结构及其能量分布。采用对壁湍流多尺度结构的子波分析表明,在湍流边界层近壁区域产生规则的流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,抑制了壁湍流各种尺度湍涡结构脉动,特别抑制了能量最大尺度湍涡结构的脉动,减小由于湍流脉动引起的在湍流边界层法向和展向的动量和能量损耗,从而减小了湍流的阻力。     

基于LBM的壁湍流跨尺度能量传递结构统计

壁湍流不同尺度间能量传输特性存在着明显的各向异性, 了解能量不同尺度间传递的空间分布是进一步构造高保真各向异性大涡模拟亚格子模式的前提. 基于格子Boltzmann数值(lattice Boltzmann method, LBM)模拟方法对雷诺数$Re_{\tau } =180$的槽道湍流进行直接数值模拟. 结果与公开的槽道湍流数据库进行对比, 平均速度剖面、雷诺应力和脉动速度均方根均取得了较好的一致性, 验证了LBM方法模拟槽道湍流的可靠性. 对模拟后的数据采用空间滤波方法得到不同尺度间能量交换量的空间分布场, 结合结构识别捕捉方法——聚类分析法对不同尺度间能量传输结构的空间分布特性进行了分析. 结果表明尺度间能量传输结构在全流场物理空间中主要为小尺寸结构, 结构的体积概率密度呈现出$-$4/5幂律, 按结构距壁面最小距离以结构距壁面距离又可将结构划分为附着结构和分离结构, 其中附着结构以较小的数量占比达到了较高的体积占比, 表明附着结构多为大尺寸结构, 进一步的对附着结构的统计表明结构在尺寸上存在着一定的幂律关系, 表明不同尺度间能量输运结构也具有Townsend提出的附面涡的自相似性, 最后对能量正反传附着结构的成对特性研究发现, 能量正传$\!-\!$反传结构对倾向于沿展向并排排列.      

A comparative study of inflow conditions for two‐ and three‐dimensional spatially developing mixing layers using large eddy simulation

A series of spatially developing mixing layers are simulated using the large eddy simulation (LES) technique. A hyperbolic tangent function and data derived from boundary layer simulations are used to generate the inflow condition, and their effects on the flow are compared. The simulations are performed in both two and three dimensions. In two‐dimensional simulations, both types of inflow conditions produce a layer that grows through successive pairings of Kelvin–Helmholtz (K–H) vortices, but the composition ratio is lower for the hyperbolic tangent inflow simulations. The two‐dimensional simulations do not undergo a transition to turbulence. The three‐dimensional simulations produce a transition to turbulence, and coherent structures are found in the post‐transition region of the flow. The composition ratio of the three‐dimensional layers is reduced in comparison to the counterpart two‐dimensional runs. The mechanisms of growth are investigated in each type of simulation, and amalgamative pairing interactions are found in the pre‐transition region of the three‐dimensional simulations, and throughout the entire computational domain of those carried out in two‐dimensions. The structures beyond the post‐transition region of the three‐dimensional simulations appear to behave in a much different manner to their pre‐transition cousins, with no pairing‐type interactions observed in the turbulent flow. In order to accurately simulate spatially developing mixing layers, it is postulated that the inflow conditions must closely correspond to the conditions present in the reference experiment. Copyright © 2007 John Wiley & Sons, Ltd.