非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

涡旋波流场的涡量测试与计算及特征参数的影响

从实验和模拟计算两个方面,对涡旋波流场的速度矢量图和涡量进行了定量测试和综合分析,并对涡旋波流动与其特征影响参数之间的关系进行了研究。结果表明:涡旋波流动只在二维层流状态下产生,涡旋波形成后,Re数对旋涡尺寸的影响较小,涡量随Re数的增大而增大,但涡量不随时间的增加而单调增加。随着St数的减小,旋涡尺寸明显增大,涡量却随之减小。本文也对不同槽道内的涡旋波流场进行了数值模拟计算,进而确定了涡旋波流场的形成条件。     

用PIV测试涡旋波流场的速度和剪应力分布

采用PIV瞬态流场测试技术,对二维槽道中的涡旋波流场不同相位上的速度分布和应力分布进行了测试和计算,本文定量地描述了槽道中涡旋波的形成过程及发展规律。通过调整振荡流的振幅和活塞行程,分析了Re数和Sr数对涡旋波流动的影响,得出了槽道内剪应力的分布状况以及平均剪切应力的周期性变化规律。深入分析了涡旋波流场强化传递现象的本质。     

槽道内涡波流场的POD分析

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

槽道内涡波流场的POD分析

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

分区Lagrangian涡方法及瞬时起动圆柱初期流动的模拟

本文提出一种分区Lagrangian涡方法:将附着流动和分离流动分开处理,在附着区解边界层方层,只在分离区用涡方法解N-S方程。由于将尺度不同的区域分开了,求解分离区流动的涡方法中,每一时间步上物面引出的涡数在较小程度上依赖于Re数。这样,求解高Re数流动时,流场内的涡数,因而计算机内存和时间得以大大减小。用该方法计算了瞬时起动圆柱的初期流动,与实验结果比较相符很好。     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论.通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中.对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制.虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.     

均匀旋转对圆柱水动力及流动结构的影响

基于格子Boltzmann方法 (LBM)对均匀旋转控制下的低雷诺数(Re=100)圆柱绕流问题进行了数值模拟,得到了转速比从0～10变化下,旋转控制对圆柱水动力及流动结构的影响规律.使用动态模态分解(DMD)对流场特征进行提取,并分析了施加旋转控制之后转速比对流场不同模态和增长率的影响.结果表明,随着转速比增大,圆柱下游流动结构依次呈现出卡门涡街、剪切层、反向剪切层、单侧涡和附着涡5种结构;阻力系数时均值先减小,随后在转速进入单侧涡区间后增大,升力系数与力矩系数的时均值均单调增加,同时,在出现涡脱落的两个转速区间内,水动力出现了明显的波动,且二次失稳时波动幅度更大. DMD的结果表明,圆柱下游的流动结构主要受圆柱壁面的旋转影响而发生改变并产生全新流动模态;旋转会对流动稳定性产生影响:在未充分发展阶段,旋转对流动稳定性的影响不显著,而在充分发展后,各转速下的流场不稳定模态数均远少于未充分发展阶段,随着转速比的增大,流动稳定性会产生不同程度的增强或减弱,且无涡脱落时的稳定性高于有涡脱落时,因此,通过旋转控制抑制尾涡脱落可以有效增强流动的稳定性.     

多段翼型缝翼流动速度的定常和非定常特性研究

用热线风速仪研究多段翼型前缘缝翼在不同条件下流动速度的定常性和非定常性。结合多段翼型定常流动Navier-Stokes方程数值模拟的结果,分析了迎角、后缘襟翼参数(偏角、缝道宽度、搭接量)对缝翼定常和非定常流动速度的影响规律。研究结果表明:在缝翼后缘处,流动分为缝道加速流动区、缝翼尾流区、缝翼上表面以上的主流区;缝翼尾流区流动速度非定常性主要表现在中低频率范围(2k Hz以下),而缝道加速流动区和缝翼上表面以上的主流区流动非定常性常表现出高频特性(2k Hz以上);在失速前随迎角增加,或者当襟翼偏角从20°向30°增加时,缝道流动加速,槽区涡减小;缝翼槽区涡形成和振荡是中低频率范围流动非定常性的机理,而缝翼鼻尖脱落涡是缝翼槽区涡振荡的激励因素。     

Some aspects of streamwise vortex behavior during oblique shock wave/vortex interaction

An experimental study of the flowfield generated by the interaction of a streamwise vortex having a strong wake-type axial Mach number profile and a two-dimensional oblique shock wave was conducted in a Mach 2.49 flow. The experiments were aimed at investigating the dynamics of supersonic vortex distortion and to study downstream behavior of a streamwise vortex during a strong shock wave/vortex encounter. The experiments involved positioning an oblique shock generator in the form of a two-dimensional wedge downstream of a semi-span, vortex generator wing section so that the wing-tip vortex interacted with the otherwise planar oblique shock wave. Planar laser sheet visualizations of the flowfield indicated an expansion of the vortex core in crossing a spherically blunt-nose shock front. The maximum vortex core diameter occurred at a distance of 12.7 mm downstream of the wedge leading edge where the vortex had a core diameter of more than double its undisturbed value. At distances further downstream the vortex core diameter remained nearly constant, while it appeared to become more diffused at distances far from the wedge leading edge. Measurements of vortex trajectory revealed that the vortex convected in the freestream direction immediately downstream of the bulged-forward shock structure, while it traveled parallel to the wedge surface at distances further downstream. The turbulent distorted vortex structure which formed as a result of the interaction, was found to be sensitive to downstream disturbances in a manner consistent with incompressible vortex breakdown. Physical arguments are presented to relate behavior of streamwise vortices during oblique and normal shock wave interactions. Received 7 September 1996 / Accepted 10 February 1998     

Free surface waves in equilibrium with a vortex

Finite amplitude solitary waves of uniform depth which interact with a stationary point vortex are considered. Waves both with and without a submerged obstacle are computed. The method of solution is collocation of Bernoulli's equation at a finite number of points on the free surface coupled with equations for equilibrium of a point vortex. The stream function and vortex location are found by computing a conformal map of the flow domain to an infinite strip. For a given obstacle the solutions are parametrized with respect to Froude number and vortex circulation. When no obstacle is present there are two families of solutions, in one of which the amplitude of the wave increases by increasing the circulation, while in the other amplitude increases by decreasing the circulation. Beyond a certain critical Froude number the maximum amplitude wave has a sharp crest with an angle of 120 degrees. Similar behavior is observed for the flow past a submerged obstacle except that there is a critical Froude number below which there is no solution at all.     

Effect of Volume Energy Supply on the Stability of a Subsonic Vortex Flow

Calculations of the stability of an axisymmetric vortex flow of viscous heat-conducting gas with volume energy supply are presented. The unperturbed axisymmetric vortex flow was found numerically using a quasi-cylindrical approximation of the Navier-Stokes equations under the assumption of constant peripheral-velocity circulation in the ambient co-current flow. The volume energy supply in the viscous vortex core was modeled by an additional source term in the energy equation. The stability characteristics of the viscous vortex flow in a longitudinal vortex with respect to both axisymmetric and non-axisymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wave number were found using the time-dependent formulation of the linear stability theory for compressible three-dimensional plane-parallel flows.     

Self-intensification in shock wave and vortex interaction

Computer simulation has been performed for the interaction between a shock wave and a vortex ring moving toward the wave. The computed density contours are compared with the pattern of shadowgraphs. A remarkable property found in the simulation is that, during the passage of the shock wave over the vortex ring, the part of the wave propagating through the inside of the ring-vortex is intensified spontaneously at a localized region. Maximum pressure occurs inside the vortex and attains a high value, about several times that of the impinging shock for incident Mach numbers of around 1.2 with the vortex translation Mach number 0.60. This is due to a double-step mechanism of intensification within the flow field by the shock-vortex interaction.     

Application of PIV to velocity measurements in a liquid film flowing down an inclined cylinder

PIV technique is applied for measurements of instant velocity distributions in a liquid film flowing down an inclined tube in the form of a wavy rivulet. An application of special optical calibration is applied to correct distortion effects caused by the curvature of the interface. A vortex flow of liquid is observed inside a wave hump in the reference system moving with wave phase velocity. Conditionally averaged profiles of longitudinal and transverse components of liquid velocity are obtained for different cross-sections of developed non-linear waves. It is shown that the increase in wave amplitude slightly changes the location of the vortex center. The analysis of modification of vortex motion character due to wavy flow conditions, such as tube inclination angle, film Reynolds number, wave excitation frequency, is fulfilled.     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。