绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

绕振荡水翼流动及其转捩特性的数值计算研究

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明:与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.      

柔性非定常细长升力面理论及涡面作用

本文讨论柔性非定常升力面的细长体理论。对于前缘问题,后缘问题和升力面外存在有旋涡的第二前缘问题的的流动作了分别的讨论,得出了统一的理论和流动的解。在后缘问题中联立求解了Kutta条件涡面积分方程和涡量守恒微分方程,得到了各个不同位置不同时问的旋涡分布,γ(y,x,t),因而比较清楚地阐明了后缘旋涡脱落和运动的机理以及非定常涡面的发展过程。对于第二前缘问题,考虑了涡面离开流动对称面的位移影响,并由此克服了涡面和物面重叠带来的困难。分析表明,当重叠部分的涡面无限接近物体(升力面)表面时,涡面的作用趋于消失。此外,在分析旋涡作用时发现,截面上涡面的流体动力效果可以由一个单一参数(记为W)来表示。W是有速度的量纲,物理上它是在物面上旋涡诱导速度的加权平均值。     

射流对绕水翼云空化流动抑制机理研究

为理解绕水翼云空化流动的发展机理和探究水翼吸力面开孔射流的影响,采用密度修正的RNG k-ε湍流模型和Schnerr-Sauer空化模型对原始NACA66(mod)水翼和采用射流后的水翼的云空化非定常过程进行模拟和对比分析;采用在水翼吸力面近壁区设立监测线的方法对近壁区的流场进行监测,得到近壁区汽相体积分数、回射流速度、压力及压力梯度的时空分布云图;开展了云空化流场特性的涡动力学分析,进而分析水翼云空化的发生机理和射流抑制空化的抑制机理.结果表明:游离型空泡在下游溃灭时产生强烈的局部高压,其向上游传播导致前缘空穴的一次回缩,而空穴的二次回缩受回射流的影响.回射流的发展区域受限于较高的压力梯度,高的压力梯度一直存在,但回射流在一个周期内的首次出现需要时间的积累.在水翼吸力面射流使得射流孔附近压力升高,弥补了由于空化和绕流造成的压降,压力梯度增大,抗逆压能力增强,对回射流起到阻挡作用;另一方面,射流使得回射流区域面积和回射流的强度也有所减小,从而对云空化的发展起到抑制的效果.Q准则的涡结构云图相比于汽相体积分数云图能显示复杂的流动结构,前缘附着型空穴和尾缘游离型空穴内存在旋涡,回射流对空穴存在剪切作用造成空穴脱落.而射流对空穴和回射流的剪切和阻挡使云空化发展得到抑制.     

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.      

主动射流控制水翼空化的数值模拟与分析

为了改善高速流动工况下水翼吸力面上流场的空化特性,提出了水翼表面主动射流对绕水翼周围流动加以控制的方法.基于密度分域滤波的FBDCM混合湍流模型联合Zwart-Gerber-Belamri空化模型,分析了来流空化数为0.83,来流攻角为8°,射流位置距水翼前缘为x=0.19c时,主动射流对于水翼吸力面上流动的空化特性和水动力特性影响.对回射流的强度进行了量化分析,以探究回射流与流场空化特性的关系.数值分析结果表明,在射流水翼吸力面上的时均空泡体积为原始水翼的1/15,使得流场内空化流动由云空化状态转变为较为稳定的片空化状态,显著地削弱了云空化的发展.此外,射流极大地改善了水翼的水动力性能,使得水翼的升阻比较原始水翼提高了22.9%,空泡的脱落频率减少了26.2%,空泡脱落所引起的振幅减小了9.1%.射流大幅降低了水翼吸力面上低压区面积,水翼吸力面上流体的逆向压力减小,回射流强度降低;同时,射流使水翼吸力面上的边界层减薄,增强了流动的抗逆压梯度能力,一定程度上阻挡了回射流向水翼前缘的流动,这也从机理上分析了主动射流抑制空化的原因.     

射流对绕水翼云空化流动抑制机理研究

为理解绕水翼云空化流动的发展机理和探究水翼吸力面开孔射流的影响,采用密度 修正的RNG $k$-$\varepsilon $湍流模型和Schnerr-Sauer空化模型对原始NACA66(mod) 水翼和采用射流后的 水翼的云空化非定常过程进行模拟和对比分析;采用在水翼吸力面近壁区设立监测线的方法对近壁区的流场进行监测,得到 近壁区汽相体积分数、回射流速度、压力及压力梯度的时空分布云图;开展了云空化流场特性的涡动力学分析,进而分析水 翼云空化的发生机理和射流抑制空化的抑制机理. 结果表明:游离型空泡在下游溃灭时产生强烈的局部高压,其向上游传播 导致前缘空穴的一次回缩,而空穴的二次回缩受回射流的影响. 回射流的发展区域受限于较高的压力梯度,高的压力梯度一 直存在,但回射流在一个周期内的首次出现需要时间的积累. 在水翼吸力面射流使得射流孔附近压力升高,弥补了由于空化 和绕流造成的压降,压力梯度增大,抗逆压能力增强,对回射流起到阻挡作用;另一方面,射流使得回射流区域面积和回射 流的强度也有所减小,从而对云空化的发展起到抑制的效果. $Q$准则的涡结构云图相比于汽相体积分数云图能显示复杂的 流动结构,前缘附着型空穴和尾缘游离型空穴内存在旋涡,回射流对空穴存在剪切作用造成空穴脱落. 而射流对空穴和回射 流的剪切和阻挡使云空化发展得到抑制.      

非定常流动的快速拉格朗日涡方法数值模拟

采用快速拉格朗日涡方法数值模拟有复杂旋涡运动的非定常流动. 利用离散涡元模拟旋涡的产生、聚集和输送过程. 拉格朗日描述法用来计算离散涡元的移动,而移动速度则利用广义毕奥-萨伐尔公式结合快速多极子展开法计算,修正的涡半径扩散模型用来模拟离散涡元的黏性扩散. 突然起动圆柱和大攻角下突然起动翼型的非定常有涡流动的数值模拟,及其与试验结果的对比验证了方法的有效性. 另外,大攻角下突然起动翼型的计算结果给出了翼型起动后吸力面旋涡的产生、发展,周期性非定常流动的形成,以及尾流旋涡结构等一些重要的流动特征.[关键词] 非定常流有涡流动快速涡方法      

非定常空化流动涡旋运动及其流体动力特性

采用大涡模拟方法对绕水翼云状空化的水动力特性和非定常流场结构进行研究. 基于实验结果对数值方法进行验证,分析空化与流场内部涡旋结构之间的相互作用以及对水翼动力特性的影响. 研究结果表明:大涡模拟方法可以准确模拟绕水翼流动的非定常过程. 在无空化条件下,升阻力系数存在斯特劳哈数St = 0.85 的主频波动,这是由水翼尾部涡旋结构的发展脱落引起的;在云状空化条件下,升阻力系数存在St = 0.34 的高能量密度低频波动,这是由大规模云状空泡团的发展和脱落引起的;云状空化阶段的升阻力系数在St = 0.5~1.5 的范围内都存在较高的波动,这是由于空化现象对水翼尾缘涡旋结构的发展和脱落产生影响,在不同发展阶段,空化现象不同程度地降低尾缘涡旋结构脱落频率.      

旋翼翼型动态失速流场特性PIV试验研究及L-B模型修正

为测量翼型动态失速的非定常涡流场特性,采用3D-PIV 技术,对典型直升机旋翼翼型SC1095 的动态失速流场特性进行测量,发现涡在不同位置处的输运速度不同:位于翼型表面的涡的无量纲速度为0.39,位于尾迹区的涡的无量纲速度为0.55. 利用前缘涡输运速度变化这一特征,改进了经典的翼型动态失速利什曼-贝多斯（Leishman-Beddoes,L-B）模型,将该模型中固定的涡时间常数修正为可以随涡位置变化的时变函数,修正后的模型计算得到翼型法向力峰值相对原L-B 模型提升5%,力矩系数负峰值相对原L-B 模型提升13%,与试验值相比更加吻合,表明修正后的翼型动态失速模型更好地体现了翼型前缘涡的物理特征.      

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

Large aerodynamic forces on a sweeping wing at low Reynolds number

The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically. After an initial start from rest, the wing is made to execute an azimuthal rotation (sweeping) at a large angle of attack and constant angular velocity. The Reynolds number (Re) considered in the present note is 480 (Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root). During the constant-speed sweeping motion, the stall is absent and large and approximately constant lift and drag coefficients can be maintained. The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows. Soon after the initial start, a vortex ring, which consists of the leading-edge vortex (LEV), the starting vortex, and the two wing-tip vortices, is formed in the wake of the wing. During the subsequent motion of the wing, a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength. This prevents the LEV from shedding. As a result, the size of the vortex ring increases approximately linearly with time, resulting in an approximately constant time rate of the first moment of vorticity, or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wingspan. The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force. The project supported by the National Natural Science Foundation of China (10232010)     

Asymmetric vortex shedding flow past an inclined flat plate at high incidence

This paper reports an experimental investigation of the vortex shedding wake behind a long flat plate inclined at a small angle of attack to a main flow stream. Detailed velocity fields are obtained with particle-image velocimetry (PIV) at successive phases in a vortex shedding cycle at three angles of attack, α=20°, 25° and 30°, at a Reynolds number Re≈5,300. Coherent patterns and dynamics of the vortices in the wake are revealed by the phase-averaged PIV vectors and derived turbulent properties. A vortex street pattern comprising a train of leading edge vortices alternating with a train of trailing edge vortices is found in the wake. The trailing edge vortex is shed directly from the sharp trailing edge while there are evidences that the formation and shedding of the leading edge vortex involve a more complicated mechanism. The leading edge vortex seems to be shed into the wake from an axial location near the trailing edge. After shedding, the vortices are convected downstream in the wake with a convection speed roughly equal to 0.8 the free-stream velocity. On reaching the same axial location, the trailing edge vortex, as compared to the leading edge vortex, is found to possess a higher peak vorticity level at its centre and induce more intense fluid circulation and Reynolds stresses production around it. It is found that the results at the three angles of attack can be collapsed into similar trends by using the projected plate width as the characteristic length of the flow.     

Parameter dependence of vortex interactions on a two-dimensional plunging plate

The structure and dynamics of the flow field created by a plunging flat-plate airfoil are investigated at a chord Reynolds number of 10,000 while varying plunge amplitude and Strouhal number. Digital particle image velocimetry measurements are used to characterize the shedding patterns and the interactions between the leading- and trailing-edge vortex structures (LEV and TEV), resulting in the development of a wake classification system based on the nature and timing of interactions between the leading- and trailing-edge vortices. The streamwise advancement of the LEV during a plunge cycle and its resulting interaction with the TEV is primarily dependent on reduced frequency; however, for Strouhal numbers above approximately 0.4, significant changes are observed in the formation of vortices shed from the leading and trailing edges, as well as the circulation of the leading-edge vortex. The functional form of the relationship between leading-edge vortex circulation and Strouhal number suggests that the Strouhal number dependence is more specifically a manifestation of the effective angle of attack. Comparison with low-Reynolds-number studies of plunging airfoil aerodynamics reveals a high degree of consistency and suggests applicability of the classification system beyond the range examined in the present work.     

应用PIV对角区非定常马蹄涡结构的实验研究

利用PIV技术研究了柱体与平板层流边界层角区的非定常流动结构,流动显示和PIV测量均 表明角区存在3种非定常的马蹄涡模态,即绕合模态、脱落-绕合模态以及脱落-耗散模态, 一定$Re$数下主涡脱落后既可能表现为脱落-绕合模态,也可能表现为脱落-耗散模态. 这主 要取决于模型头部形状对涡轴造成的拉伸以及耗散和扩散程度. PIV测量表明,随雷诺数增 加主涡下方从壁面喷发的反向二次涡逐步增大形成强度和尺度较大的``涡舌', 该``涡舌' 将突入整个涡系所在的边界层,最终将主涡与上游涡系隔离并使其从旋涡生成区涡系脱落. 马蹄涡非定常摆动时具有较复杂的奇点形态组合和演化,反映涡轴受到了交替的拉伸和压缩 作用.     

Numerical investigations of the vortex interactions for a flow over a pitching foil at different stages

The fluid–structure interaction is investigated numerically for a two-dimensional flow (Re=2.5·10⁶) over a sinusoid-pitching foil by the SST (Shear Stress Transport) k–ω model. Although discrepancies in the downstroke phase, which are also documented in other numerical studies, are observed by comparing with experimental results, our current numerical results are sufficient to predict the mean features and qualitative tendencies of the dynamic stall phenomenon. These discrepancies are evaluated carefully from the numerical and experimental viewpoints.In this study, we have utilized Λ, which is the normalized second invariant of the velocity gradient tensor, to present the evolution of the Leading Edge Vortex (LEV) and Trailing Edge Vortex (TEV). The convective, pressure, and diffusion terms during the dynamic stall process are discussed based on the transport equation of Λ. It is found that the pressure term dominates the rate of the change of the rotation strength inside the LEV. This trend can hardly be observed directly by using the vorticity transport equation due to the zero baroclinic term for the incompressible flow.The mechanisms to delay the stall are categorized based on the formation of the LEV. At the first stage before the formation of the LEV in the upper surface, the pitching foil provides extra momentum into the fluid flows to resist the flow separation, and hence the stall is delayed. At the second stage, a low-pressure area travels with the evolution of the LEV such that the lift still can be maintained. Three short periods at the second stage corresponds to different flow patterns during the dynamic stall, and these short periods can be distinguished according to the trend of the pressure variation inside the LEV. The lift stall occurs when a reverse flow from the lower surface is triggered during the shedding of the LEV. For a reduced frequency k_f=0.15, the formation of the TEV happens right after the lift stall, and the lift can drop dramatically. With a faster reduced frequency k_f=0.25, the shedding of the LEV is postponed into the downstroke, and the interaction between the LEV and TEV becomes weaker correspondingly. Thus, the lift drops more gently after the stall. In order to acquire more reliable numerical results within the downstroke phase, the Large Eddy Simulation (LES), which is capable of better predictions for the laminar-to-turbulent transition and flow reattachment process, will be considered as the future work.     

Measurement of parallel blade–vortex interaction at low Reynolds numbers

In this study parallel blade–vortex interaction for a Schmidt-propeller configuration has been examined using particle image velocimetry (PIV). This tandem configuration consists of a leading airfoil (forefoil), used to generate a vortical wake of leading-edge vortices (LEVs) and trailing-edge vortices (TEVs) through a pitching or plunging motion, and a trailing airfoil (hindfoil), held fixed with a specified angle of attack and vertical spacing in its wake. The hindfoil incidence (loading) and not the vertical spacing to the incoming vortical wake has been found to dictate the nature of the interaction (inviscid vs. viscous). For cases where the vortex–blade offset is small and the hindfoil is loaded, vortex distortion and vortex-induced separations are observed. By tracking the circulation of the LEV and TEV, it has been found that the vortices are strengthened for the tandem arrangement and in certain cases dissipate quicker in the wake when interacting with the hindfoil. Time-averaged forces obtained using a standard control-volume analysis are then obtained and used to evaluate these vortex-interaction cases. A subsequent analysis of the varying pressure distribution over the suction side of the hindfoil is performed by integrating the Navier–Stokes equations through the velocity field. This allows for a direct comparison of the vortex-induced loading for the various configurations.     

Unsteady flow structure and vorticity convection over the airfoil oscillating at high reduced frequency

Unsteady vortex structures and vorticity convection over the airfoil (NACA 0012), oscillating in the uniform inflow, are studied by flow visualization and velocity measurements. The airfoil, pivoting at one-third of the chord, oscillates periodically near the static stalling angle of attack (AOA) at high reduced-frequency. The phase-triggering and modified phase-averaged techniques are employed to reconstruct the pseudo instantaneous velocity field over the airfoil. During the down stroke cycle, the leading-edge separation vortex is growing and the vortex near the trailing edge begins to shed into the wake. During the upstroke cycle, the leading-edge separation vortex is matured and moves downstream, and the counter clockwise vortex is forming near the trailing edge. Convection speeds and wavelength of the unsteady vortex structure over the airfoil equal to that of the counter clockwise vortex shed into the wake. This kind of vortex structure is termed as “synchronized shedding” type. The wavelength of unsteady vortex structure over the airfoil is significantly different from that at low reduced-frequency. Consistent convection speeds of the leading-edge separation vortex are acquired from the spatial-temporal variations of local circulation and local surface vorticity generation, and equals that predicted from flow visualization. Spatial-temporal variations of the local surface vorticity generation clearly reveal the formation and passage of the leading-edge separation vortex only in the region where the flow does not separate completely from the surface. Significant amounts of the surface vorticity are generated within the leading-edge region of the airfoil during the upstroke cycle. Only negligible amount of surface vorticity is produced within the region of complete flow separation. During the down stroke cycle, the surface vorticity generation is mild along the airfoil surface, except the leading-edge region where a small scale leading-edge separation vortex is forming and growing.     

The influence of surface roughness on cloud cavitation flow around hydrofoils

The aim of this study is to investigate experimentally the effect of surface roughness on cloud cavitation around Clark-Y hydrofoils. High-speed video and particle image velocimetry(PIV) were used to obtain cavitation patterns images(Prog. Aerosp. Sci. 37: 551–581, 2001), as well as velocity and vorticity fields. Results are presented for cloud cavitating conditions around a Clark-Y hydrofoil fixed at angle of attack of α = 8?for moderate Reynolds number of Re = 5.6 × 10~5. The results show that roughness had a great influence on the pattern, velocity and vorticity distribution of cloud cavitation. For cavitating flow around a smooth hydrofoil(A) and a rough hydrofoil(B), cloud cavitation occurred in the form of finger-like cavities and attached subulate cavities, respectively. The period of cloud cavitation around hydrofoil A was shorter than for hydrofoil B.Surface roughness had a great influence on the process of cloud cavitation. The development of cloud cavitation around hydrofoil A consisted of two stages:(1) Attached cavities developed along the surface to the trailing edge;(2) A reentrant jet developed, resulting in shedding and collapse of cluster bubbles or vortex structure. Meanwhile, its development for hydrofoil B included three stages:(1) Attached cavities developed along the surface to the trailing edge, with accumulation and rotation of bubbles at the trailing edge of the hydrofoil affecting the flow field;(2) Development of a reentrant jet resulted in the first shedding of cavities. Interaction and movement of flows from the pressure side and suction side brought liquid water from the pressure side to the suction side of the hydrofoil, finally forming a reentrant jet. The jet kept moving along the surface to the leading edge of the hydrofoil, resulting in large-scale shedding of cloud bubbles. Several vortices appeared and dissipated during the process;(3) Cavities grew and shed again.