基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×106时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

通气空泡尾部气泡流仿真和特性分析

基于欧拉-欧拉双流体模型对通气空泡尾部气泡流进行数值仿真, 并采用基于多尺寸分组模型的总体平衡方法预估气泡尺寸分布. 应用改进后湍流耗散系数计算模型, 考虑了气泡体积含量对在湍流作用下气泡扩散现象的影响. 基于上述模型对两种试验工况下流场进行了数值仿真.结果表明模型对空泡尾部回流区特性进行了准确预示, 在回流区高湍流度作用下气泡迅速破碎成小气泡. 并进一步得到试验体尾流区空泡体积分数和速度分布. 尾流区水流速度分布保持了流体经过非流线型对称体时产生的尾流分布规律. 仿真结果与试验数据相一致, 模型适用性得到验证.      

高速列车近尾流区湍流涡旋的数值计算及动力学讨论

主要针对缩尺比例1:30的高速动车组空气动力学模型的近尾流区域流场进行数值模拟,分析讨论了湍动能和湍流能量的产生,得到以下结论:对应于各展向位置的湍动能沿流向的变化规律与近尾流区涡旋结构携带能量的展向外移现象有关;在尾车鼻端附近,湍流涡旋具有显著的湍动能,并且各方向上的能量分量具有相同的量级,反映出湍流涡旋是高度三维的流动结构;靠近尾车鼻端的近尾流区湍流涡旋具有较强的从平均流动中提取能量的能力,并且结果表明,来自于车体底面和侧面的剪切流动沿垂向分别在一定空间范围内发挥重要的影响作用;大涡特征尺度沿流向增大,其中较小的涡旋对湍流能量的产生有主要贡献,同时,由于受到列车侧面较厚剪切层的影响,对湍流能量有贡献的涡旋所对应的积分尺度范围增大,从而使位于尾车鼻端附近的涡流能够在更大的流向范围内获取用于维持湍流尾流的能量.     

螺旋桨梢涡不稳定性机理与演化模型研究

螺旋桨尾流场的涡流特性是一个基础但又十分复杂的流体力学问题, 它的复杂性源于其蕴含复杂的漩涡系统, 且该漩涡系统会在高速的剪切层流动中不断演化, 其流体动力学行为, 如由稳定态演变为不稳定态的机理以及复杂工况环境中的流动现象, 一直是流体力学领域的难点和备受关注的热点问题. 从工程应用的角度看, 桨后梢涡的演化特性与船舶结构物的宏观特性直接相关, 更好地理解多工况下螺旋桨尾流的动力学特性, 将有助于改善与振动、噪声以及结构问题等相关的推进器性能, 对综合性能优良的下一代螺旋桨的设计和优化有着重要的现实意义. 本文基于延迟分离涡模拟、大涡模拟和无湍流模型模拟方法以及粒子图像测速流场测试分别开展了螺旋桨尾流动力学特性的数值与试验研究, 对螺旋桨尾流不稳定性的触发机理进行了揭示. 基于均匀来流中螺旋桨梢涡的演化机理, 提出了螺旋桨梢涡演化模型. 该模型能够较为准确地模拟螺旋桨梢涡的演化过程, 预测螺旋桨梢涡融合的时间和位置, 对螺旋桨流噪声预报和控制以及性能优良的螺旋桨设计具有重要意义.      

翼梢小翼若干几何参数对翼尖涡流场的影响研究

为解决使用CFD方法研究翼尖涡流场时存在的传统湍流模式难以准确计算翼尖涡的强旋流场的问题,本文使用添加旋转修正的SST-RC模式计算了光滑机翼的翼尖涡流场,并通过与实验结果进行对比,验证了该模式对于翼尖涡流场的有效性。然后使用该模式计算了翼梢小翼后的翼尖涡流场,发现翼梢小翼作为一种涡扩散器能够将原来强度集中的翼尖涡破碎分解成两个旋涡,这两个旋涡互相作用,加速了翼尖涡的粘性耗散。通过改变小翼的三个几何参数:根弦长Cr、倾斜角θ、外撇角,计算了不同外形翼梢小翼后的翼尖涡流场,同时引入一个翼尖涡的强度尺度S,考察S在尾缘下游的分布,发现添加小翼后S普遍下降了10%~20%。改变小翼根弦长对S的衰减影响不太明显,不同根弦长对应的S值相差2%~4%,Cr=0.7c时效果最佳;改变小翼倾斜角θ对S的影响在初始时差别较大,不同倾斜角对应的S相差约9%,在下游后期S差别较小,相差约3%,θ=20°时效果最佳;与倾斜角θ类似,改变小翼外撇角对S的影响在初始时差别较大,不同外撇角对应的S相差约10%,在下游后期S差别较小,相差约1%,=5°时效果最佳。     

舵几何特征对桨-舵系统尾流场演化的影响

发生在桨和舵之间的干扰会影响螺旋桨尾流的演化,导致尾流场中的湍流在下游增强,恶化船舶的振动和噪声性能,深入分析舵几何参数对桨-舵系统尾流场演化的影响能够为推进器尾流场的调节和减振降噪提供新思路.因此,从弦长、剖面和梯形舵入手分析不同的舵几何参数对螺旋桨尾流场演化特性的影响,使用大漩涡模拟方法模拟流场中的湍流结构,对不同舵弦长、剖面下的螺旋桨尾涡结构演化进行了分析,在舵弦长、剖面影响螺旋桨尾流场演化的研究的基础上分析了梯形舵对螺旋桨尾涡结构的影响,进一步分析了梯形舵影响下的螺旋桨尾流场中湍动能的分布.结果表明舵的弦长和剖面均会影响螺旋桨尾流场的演化,这种影响表现为更大的弦长和更厚的剖面会促进螺旋桨梢涡在舵压力面上的偏移,更薄的舵剖面会带来更强烈的螺旋桨毂涡偏移;涡管轮廓和舵表面脉动压力的对比均表明梯形舵会促进螺旋桨尾流场沿逆舵梯度方向偏移,从而导致螺旋桨的尾涡结构在舵两侧及下游呈现不对称分布,桨-舵系统下游的湍流结构与螺旋桨尾涡-舵碰撞过程、螺旋桨尾涡-舵随边涡干扰过程、螺旋桨梢涡-螺旋桨毂涡干扰有关,偏移更大的螺旋桨尾涡结构会在尾流场中更早地引起湍动能增强.     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

弯度翼型大迎角时分离涡流场的多平衡态现象

基于kω的SST两方程湍流模型,在时间域求解雷诺平均Navier-Stokes方程,模拟弯度翼型大迎角时的分离流动。通过给翼型施加一定形式的扰动,重点关注了翼型弯度对大迎角分离涡流场平衡态转移的影响。研究结果表明:与相同厚度20%以上的对称翼型相比,2%弯度的翼型出现分离涡流场平衡态转移的起始迎角变小2°左右,迎角区间变宽约1°;在厚度相对较小的NACA2416翼型上也发现上述分离涡平衡态转移现象。由此说明翼型弯度在一定程度上促使了分离涡平衡态的转移。     

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

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

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

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

Investigation of turbulent flow past a yawed wavy cylinder

A large eddy simulation (LES) study was conducted to investigate the three-dimensional characteristics of the turbulent flow past wavy cylinders with yaw angles from 0° to 60° at a subcritical Reynolds number of 3900. The relationships between force coefficients and vortex shedding frequency with yaw angles for both wavy cylinders and circular cylinders were investigated. Experimental measurements were also performed for the validation of the present LES results. Comparing with corresponding yawed circular cylinders at similar Reynolds number, significant differences in wake vortex patterns between wavy cylinder and circular cylinder were observed at small yaw angles. The difference in wake pattern becomes insignificant at large yaw angles. The mean drag coefficient and the Strouhal number obey the independence principle for circular cylinders at yaw angle less than 45°, while the independence principle was found to be unsuitable for yawed wavy cylinders. In general, the mean drag coefficients and the fluctuating lift coefficients of a yawed wavy cylinder are less than those of a corresponding yawed circular cylinder at the same flow condition. However, with the increase of the yaw angle, the advantageous effect of wavy cylinder on force and vibration control becomes insignificant.     

Experimental study on aerodynamic coefficients of yawed cylinders

This paper presents wind tunnel tests on a stationary cylinder inclined with the flow. The cylinder was positioned at different sets of yaw and vertical angles. The flow regime of the tests remained in the subcritical state. Two load cells were designed and installed to measure the aerodynamic forces, with enough sensitivity to measure vortex shedding frequencies. In this paper, the three aerodynamic force coefficients are normalized using the free stream velocity instead of its normal component. The results show that the drag coefficient and the resultant of the lift and side forces coefficients can be described by an empirical function of the incidence angle. The lift and side force coefficients remain however functions of both the horizontal yaw and vertical angles and cannot be expressed as functions of the incidence angle only. The Independence Principle was observed to become inaccurate for yaw angles larger than 40°. However, the measured Strouhal numbers indicate that the vortex shedding frequencies of a yawed cylinder can be predicted using the Independence Principle.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

复杂工况下的大型风力机气动性能和尾迹研究Investigation of aerodynamic performance and wake of the large scale wind turbine in complicated conditions

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

Large Eddy Simulation of Flow Control Around a Cube Subjected to Momentum Injection

The concept of Momentum Injection (MI) through Moving Surface Boundary layer Control (MSBC) applied to a cubic structure is numerically studied using Large Eddy Simulation at a Reynolds number of 6.7×10⁴. Two small rotating cylinders are used to add the momentum at the front vertical edges of the cube. Two configurations are studied with the yaw angle of 0° and 30°, respectively, with ratio of the rotation velocity of cylinders and the freestream velocity of 2. The results suggest that MI delays the boundary layer separation and reattachment, and thus reduces the drag. A drag reduction of about 6.2 % is observed in the 0° yaw angle case and about 44.1 % reduction in the 30° yaw angle case. In the case of 0° yaw angle, the main change of the flow field is the disappearance of the separation regions near the rotating cylinders and the wake region is slightly changed due to MI. In the 30° yaw angle case, the flow field is changed a lot. Large flow separations near one rotating cylinder and in the wake is significantly reduced, which results in the large drag reduction. Meanwhile, the yaw moment is increased about 50.5 %.     

Traversing field of view and AR-PIV for mid-field wake vortex investigation in a towing tank

Wake vortex flow experiments are performed in a water tank where a 1:48 scaled model of a large transport aircraft A340-300 is towed at the speed of 3 and 5 ms^-1 with values of the angle of attack !={2°, 4°, 8°}. Particle image velocimetry (PIV) measurements are performed in a plane perpendicular to the towing direction describing the streamwise component of the wake vorticity. The instantaneous field of view (I-FOV) is traversed vertically with an underwater moving-camera device tracking the vortex core during the downward motion. An adaptive resolution (AR) image-processing technique is introduced that enhances the PIV interrogation in terms of spatial resolution and accuracy. The main objectives of the investigation are to demonstrate the applicability of PIV diagnostics in wake vortex research with towing-tank facilities. The specific implementation of the traversing field-of-view (T-FOV) technique and the AR image processing are driven by the need to characterize the vortex wake global properties as well as the vortex decay phenomenon in the mid- and far-field. Relevant aerodynamic information is obtained in the mid-field where the time evolution of the vortex structure (core radius and tangential velocity) and of the overall vortex wake (vortex trajectory, descent velocity, circulation) are discussed.     

Effect of vortex shedding in unsteady aerodynamic forces for a low Reynolds number stationary wing at low angle of attack

This study elucidates the relation between wake vortex shedding and aerodynamic force fluctuations for a low Reynolds number wing from time resolved particle image velocimetry (TR-PIV) experimental measurements. The results reveal a periodic lift and drag variation within the shedding cycle and resolve the frequencies of those fluctuations from a proper orthogonal decomposition (POD) and power spectral density (PSD) analysis. To show the effect of vortex shedding on the body force fluctuations, the evolution of instantaneous aerodynamic forces is compared to the pressure field of the fluid flow and to the vortical structures in the wake of the airfoil. A six step model describing the vortex-force relation is proposed. It shows that changes in lift such as maximum lift and minimum lift are associated with the detachment of a vortex. It also shows that the minimum or local minimum drag value is obtained at the onset formation of a vortex on the airfoil wake. Similarly, the maximum or local maximum drag is obtained at the onset formation of the saddle on the airfoil wake. The model further explains the asymmetry observed in the unsteady drag force evolution. The model can be used to optimize flow control and fluid-structure interaction applications.     

A nonlinear model for aerodynamic configuration of wake behind horizontal-axis wind turbine

Determination of the aerodynamic configuration of wake is the key to analysis and evaluation of the rotor aerodynamic characteristics of a horizontal-axis wind turbine.According to the aerodynamic configuration, the real magnitude and direction of the onflow velocity at the rotor blade can be determined, and subsequently, the aerodynamic force on the rotor can be determined. The commonly employed wake aerodynamic models are of the cylindrical form instead of the actual expanding one. This is because the influence of the radial component of the induced velocity on the wake configuration is neglected. Therefore, this model should be called a "linear model". Using this model means that the induced velocities at the rotor blades and aerodynamic loads on them would be inexact. An approximately accurate approach is proposed in this paper to determine the so-called "nonlinear" wake aerodynamic configuration by means of the potential theory,where the influence of all three coordinate components of the induced velocity on wake aerodynamic configuration is taken into account to obtain a kind of expanding wake that approximately looks like an actual one. First, the rotor aerodynamic model composed of axial(central), bound, and trailing vortexes is established with the help of the finite aspect wing theory. Then, the Biot-Savart formula for the potential flow theory is used to derive a set of integral equations to evaluate the three components of the induced velocity at any point within the wake. The numerical solution to the integral equations is found,and the loci of all elementary trailing vortex filaments behind the rotor are determined thereafter. Finally, to formulate an actual wind turbine rotor, using the nonlinear wake model, the induced velocity everywhere in the wake, especially that at the rotor blade,is obtained in the case of various tip speed ratios and compared with the wake boundary in a neutral atmospheric boundary layer. Hereby, some useful and referential conclusions are offered for the aerodynamic computation and design of the rotor of the horizontal-axis wind turbine.     

Experimental investigation of a blunt trailing edge flow field with application to sound generation

The unsteady lift generated by turbulence at the trailing edge of an airfoil is a source of radiated sound. The objective of the present research was to measure the velocity field in the near wake region of an asymmetric beveled trailing edge in order to determine the flow mechanisms responsible for the generation of trailing edge noise. Two component velocity measurements were acquired using particle image velocimetry. The chord Reynolds number was 1.9 × 10⁶. The data show velocity field realizations that were typical of a wake flow containing an asymmetric periodic vortex shedding. A phase average decomposition of the velocity field with respect to this shedding process was utilized to separate the large scale turbulent motions that occurred at the vortex shedding frequency (i.e., those responsible for the production of tonal noise) from the smaller scale turbulent motions, which were interpreted to be responsible for the production of broadband sound. The small scale turbulence was found to be dependent on the phase of the vortex shedding process implying a dependence of the broadband sound generated by the trailing edge on the phase of the vortex shedding process.