用鲁棒Riemann求解器和运动重叠网格计算旋翼粘性绕流

发展了一种基于鲁棒Riemann求解器和运动重叠网格技术计算直升机悬停旋翼流场的方法。基于惯性坐标系,悬停旋翼流场是非定常流场,控制方程为可压缩Reynolds平均Navier-Stoke方程,其对流项采用Roe近似Reimann求解器离散,使用改进的五阶加权基本无振荡格式进行高阶重构,非定常时间推进采用含牛顿型LUSGS子迭代的全隐式双时间步方法。为实施旋转运动和便于捕捉尾迹,计算采用运动重叠网格技术。计算得到的桨叶表面压力分布及桨尖涡涡核位置都与实验结果吻合较好。数值结果表明:所发展方法对桨尖涡具有较高的分辨率,对激波具有较好的捕捉能力,该方法可进一步推广到前飞旋翼粘性绕流的计算。     

基于Navier-Stokes／Kirchhoff方程的悬停旋翼跨声速气动声学计算

研究了Kirchhoff积分面是否有盖有底,以及是否计及旋翼网格上的流场值,这两个因素对噪声预测结果的影响.发展了一种基于重叠网格的计算悬停旋翼远场噪声的数值方法.数值计算过程分为流场模拟和声场模拟两部分.悬停旋翼流场的数值模拟是在两个相互重叠的网格上进行的:在高质量的旋翼网格上求解Navier-Stokes方程,用于模拟旋翼附近的粘性流动和近场尾涡的捕捉;在远离粘性区域处布置符合悬停流场物理特征的圆柱形背景网格,控制方程为Euler方程,用于远场尾涡的捕捉.计算得到的流场信息插值到用于声场计算的Kirchhoff积分面上.观测点处的噪声可以认为是由这个完全包含桨叶的Kirchhoff积分面上的面元(声源)发声得到.远场声波的传播由Kirchhoff积分公式描述.计算结果表明:采用有盖有底的Kirchhoff积分面并且同时计及旋翼网格流场值时,计算得到的HSI噪声与实验值吻合最好.     

微型四旋翼飞行器最佳旋翼间距布局研究

为研究悬停状态下不同翼间距对微型四旋翼飞行器气动性能的影响,结合整机试验和数值模拟,分析了不同旋翼间距下微型四旋翼飞行器拉力和功耗的变化规律。在样机试验中,通过搭建试验平台对间距比l/d范围1.1～2.0的微型四旋翼飞行器进行了拉力和功耗的测量,确定了相同功耗条件下具有较大拉力的最佳旋翼间距范围。为更直观地得到旋翼间气动干扰对整机气动性能的影响,通过CFD方法对微型四旋翼飞行器流场进行了仿真,得到了不同间距下的压力、流线和涡量分布情况,进而对四旋翼飞行器在不同旋翼间距下表现出的不同气动特性进行对比。结果表明,与无干扰状态下的孤立单旋翼相比,四旋翼间存在的气动干扰在合理的旋翼间距下可以保持涡流完整,并有助于提升四旋翼系统的拉力。最后,通过试验和仿真对比发现,在旋翼间距为1.8d时,四旋翼飞行器具有较大的功率载荷和良好的气动特性,是该四旋翼飞行器整机的最佳气动布局。     

基于动态面搭接技术的直升机旋翼流场分析

为了进一步提高直升机旋翼流场计算的精度与效率,本文根据直升机旋翼的运动特点,发展了一种基于结构化网格动态面搭接技术及低速预处理技术的直升机旋翼流场的非定常数值分析方法。采用该方法分别对Caradonna-Tung、UH-60A旋翼模型的悬停流场及NASA某旋翼无升力前飞流场进行了数值分析;将计算所得的旋翼表面压力分布与风洞试验结果进行了对比,对旋翼的尾涡结构进行了相关分析。分析结果表明:本文所发展的方法具有较高的精度且计算速度快,并能够比较准确地捕捉直升机旋翼的尾涡结构,可以用于直升机旋翼流场的分析研究,从而为开展考虑挥舞及变距运动的直升机旋翼前飞状态非定常气动特性的研究奠定了基础。     

基于运动嵌套网格的前飞旋翼桨叶气动干扰数值模拟

采用运动嵌套网格和求解Navier-Stokes方程的方法,数值模拟了前飞状态的旋翼流场,结果同试验数据吻合得很好;在此基础上,分别对单双桨叶旋翼在两种桨距情况下的前飞流场,进行了对比计算,模拟了各种因素对旋翼桨叶间气动干扰作用的影响,并对其机理进行了探讨,得到一些与工程实际吻合的现象和结论.     

旋翼尾流与地面干扰时地面涡现象的研究

用Ｎ－Ｓ方程对近地飞行时旋翼尾流与地面干扰时产生的地面涡现象进行了数值计算旋翼对流场的作用由分布在特定区域内的动量源项模拟结果表明，旋翼尾流撞到地面后的卷起和轴向流动的拉伸作用是形成地面涡的原因；地面边界层形成的二次分离涡向地面涡内输入（与尾流所携带的涡量）相反的涡量，而使地面涡保持平衡；地面涡的存在和运动使旋翼附近流场大大改变     

桨-舵间距和雷诺数对潜艇尾部伴流场均匀性的影响研究

采用Realized k-ε湍流模型求解Reynolds平均Navier-Stokes方程,对SUBOFF全附体模型在低雷诺数和高雷诺数条件下的三维粘性流场分别进行了数值模拟,研究了桨-舵间距以及雷诺数对螺旋桨盘面伴流场的影响.计算结果表明:适当的加大桨-舵间距对提高潜艇尾流场的均匀性有显著的效果;桨-舵间距对潜艇桨盘面伴流场均匀性的影响在潜艇的低雷诺数状态和高雷诺数状态具有同样的规律.     

侧向多喷口干扰复杂流动数值模拟研究

采用具有高分辨率的NND格式,通过数值求解N-S方程对典型外形多喷口侧向喷流复杂干扰流动进行了数值模拟. 为了提高计算效率,采用了LU-SGS隐式算法. 采用分块对接网格技术,生成高质量的贴体计算网格,精确模拟喷口截面. 对比分析了不同计算格式、限制器形式、网格拓扑及流动形态(层流与湍流)对喷流干扰流场结构和压力分布特性的影响,研究和分析了喷口附近流场的涡系结构、波系结构和喷流干扰引起的气动力特性. 在上述研究的基础上,针对典型飞行器外形的侧向喷流干扰特性进行了详细的数值模拟,得到了喷口参数(喷口位置、数目等)及来流条件对喷流干扰流场结构、气动力特性的影响规律,并对其流动机理进行了相应的分析. 研究表明,发展的针对多喷口侧喷干扰的数值计算方法是成功的,可以应用于飞行器侧向喷流干扰的流场结构分析及气动力特性数值预测.      

全附体潜艇粘性流场的RANS模拟及场量和涡量的校验分析

采用相同拓扑结构和相近网格质量的4套网格和5种湍流模型,对全附体Suboff潜艇粘性流场进行RANS模拟,分析了网格密度、节点空间分布规律和湍流模型对计算精度的影响,详细校验了其力积分量、速度场量和涡量特征。结果表明:网格密度最大的G4网格(140万)计算精度最高,总阻力较实验值误差为0.723%,其采用SST湍流模型时最优。计算得到的压力系数和剪切应力系数分布均与实验值吻合很好;桨盘面速度等值线分布计算精度与文献相当,轴向相对速度0.9以上的计算半径稍大于实验值,其余半径与实验吻合较好;桨盘面上0.25倍半径处速度分量沿周向分布计算精度较文献高,轴向分量与实验值吻合较好,径向分量峰值稍小于实验值,但峰值所处周向位置与实验值一致。成功捕捉到了附体端面绕流诱导对旋涡、附体叶根截面下游处项链形涡对、尾翼端面尾缘上方附着涡蹄、附体马蹄涡系、尾翼截面通道流体挤压作用诱导涡以及桨盘面涡量汇集的潜艇涡量场特征,且围壳端面绕流诱导对旋涡沿流动方向持续稳定,不影响桨盘面涡量场,均与文献中由大涡模拟模拟得到的定性结论一致。研究表明,在网格密度较大、节点分布合理、网格质量较高、湍流模型选取适当和壁面函数使用有效的条件下,RANS模拟潜艇粘性流场的场量和涡量特征同样具有很高的计算精度,能够在工程应用中有力支撑新型艇型设计与性能分析。     

90°弯管内流动的理论模型及流动特性的数值研究

从三维不可压缩雷诺时均Navier-Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk-ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。     

Euler solutions for self-generated rotor blade-vortex interactions

A finite-difference procedure has been developed for the prediction of three-dimensional rotor blade-vortex interactions. The interaction velocity field was obtained through a non-linear superposition of the rotor flow field, computed using the unsteady three-dimensional Euler equations, and the embedded vortex wake flow field, computed using the law of Biot-Savart. In the Euler model, near wake rotational effects were simulated using the surface velocity ‘transpiration’ approach. As a result, a modified surface boundary condition was prescribed and enforced at each time step of the computations to satisfy the tangency boundary condition. For supercritical interactions using an upstream-generated vortex, accuracy of the numerical results were found to rely on the user-specified vortex core radius and vortex strength. For the more general self-generated subcritical interactions, vortex wake trajectories were computed using the lifting-line helicopter/rotor trim code CAMRAD. For these interactions, accuracy of the results were found to rely heavily on the CAMRAD-predicted vortex strength, vortex orientation with respect to the blade, and to a large extent on the user-specified vortex core radius. Results for the one-seventh scale model OLS rotor and for a non-lifting rectangular blade having a NACA0012 section are presented. Comparisons with the experimental windtunnel data are also made.     

鸭式旋翼/机翼飞机悬停及小速度前飞气动干扰实验研究

鸭式旋翼/机翼飞机是一种新概念可垂直起降高速飞行器,为了解该飞机在悬停及小速度前飞时的全机气动干扰特性,在南京航空航天大学开口风洞中进行了飞机全机气动力实验,实验采用多台测力天平分别测量主机翼和机身的气动力.结果表明,悬停时受主机翼高速旋转产生的下洗尾流影响,机身产生了较大的法向力和低头力矩;前飞时下洗尾流对机身的法向力和俯仰力矩有比较严重的干扰,对滚转力矩和偏航力矩干扰较小,对侧向力有一定影响.实验结果为飞机的飞行动力学特性研究以及控制律设计提供了参考.＃     

Effects of camber angle on aerodynamic performance of flapping-wing micro air vehicle

This study investigates the unsteady aerodynamic characteristics of the cambered wings of a flapping-wing micro air vehicle (FW-MAV) in hover. A three-dimensional fluid–structure interaction solver is developed for a realistic modeling of large-deforming wing structure and geometry. Cross-validation is conducted against the experimental results obtained also in the present study to establish more accurate analyses of cambered wings. An investigation is carried out on the unsteady vortex structures around the wings caused by the passive twisting motion. A parametric study is then conducted to evaluate the aerodynamic performance with respect to the camber angle at three different flapping frequencies including normal operating conditions. The camber angles producing the largest thrust and highest propulsive efficiency are estimated at each flapping frequency, and their effects on aerodynamic performance are analyzed in terms of the stroke phase. The timing and magnitude of the passive twisting motion, which are dependent on the camber angle at the operating frequency, greatly affects the unsteady vortex structure. Consequently, the camber angle designed at the operating frequency plays a key role in enhancing the aerodynamic performance of FW-MAVs.     

基于涡尾迹方法的风力机气动性能分析

基于(势流)涡尾迹方法开发了水平轴风力机叶片气动性能分析程序,采用固定尾迹涡模型和自由尾迹模型分别对气动设计性能进行计算分析,得到风力机设计工况下的涡位置、诱导因子、功率系数及扭矩系数等气动性能参数,并与设计结果对比。结果表明,涡尾迹方法能够快速准确地计算风力机叶片气动性能参数,对风力机叶片气动分析,固定尾迹涡模型较自由尾迹模型计算时间短,具有较好的实用性。     

Wind pressures on tapered and set-back tall buildings

Recent tall buildings tend to have irregular and unconventional shapes as a prevailing but unavoidable trend, which is very effective for suppressing across-wind responses. Suppression of across-wind responses is a major factor in safety and habitability design of tall buildings, and the so-called aerodynamic modification method is comprehensively used. While the effectiveness of aerodynamic modification in reducing wind loads has been widely reported, there have been few detailed investigations of pressure fluctuations. The purpose of the present work is to investigate the spatio-temporal characteristics of pressure fluctuations applied to height-modified tall buildings comprehensively, including differences of vortex formation and shedding mechanism based on the previously reported mechanism of a conventional square tall building. The results show that taper and set-back affect on the bandwidth of power spectra and position of peak frequencies. And through taper and set-back, the height at which the vortex begins to form moves up, and due to the small building dimension, the vortex component formed at that height sheds from the building more frequently before an inverted conical vortex is formed over the whole height.     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

蜻蜓滑翔时柔性褶皱前翅气动特性分析

蜻蜓是自然界优秀的飞行家,滑翔是其常见且有效的飞行模式.蜻蜓优异的飞行能力来源于其翅膀的巧妙结构,褶皱是蜻蜓翅膀上最为显著的结构之一,不仅提高了翅膀的刚度,还改变了其气动特性,而飞行过程中柔性翅膀会产生变形是蜻蜓翅膀的另一特性.为揭示蜻蜓在滑翔时,柔性褶皱前翅的变形,探究褶皱和柔性的共同作用对其气动特性的影响,基于逆向工程,依据前人的测量数据和研究成果,通过三维建模软件建立了蜻蜓三维褶皱前翅的计算流体力学(computational fluiddynamics,CFD)模型和计算结构力学(computational structuralmechanics,CSD)模型,并通过模态分析验证了此模型有足够的精度.基于CFD方法和CFD/CSD双向流固耦合计算方法分别对蜻蜓滑翔飞行时刚性和柔性褶皱前翅的气动特性进行了数值模拟,结果表明,柔性褶皱前翅受气动载荷后,翅脉和翅膜产生形变,柔性前翅上下表面压力差相较于刚性前翅减小了,从而其升力和阻力也减小了,而在大攻角时,变形后的前缘脉诱导出比刚性前翅更强的前缘涡.因此在攻角小于10$^\circ$时刚性前翅的气动特性优于柔性前翅,继续增大攻角,柔性前翅的气动特性则优于刚性前翅.前翅受载后气动响应时间短,翅尖的变形最大,仅仅产生了垂直于翅膀所在平面方向上的变形,而没有发生扭转,翼根处受到应力最大,褶皱上凸部分承受蜻蜓滑翔时前翅的主要载荷.      

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.     

涡波一体乘波飞行器宽速域气动优化设计研究

涡波一体宽速域乘波飞行器通过在低速引入涡效应,显著改善了传统乘波体在低速状态下的升阻特性,具有在未来宽速域空天飞行器总体气动设计当中得到广泛应用的巨大潜力.但是,该设计方法的研究尚不完善,特别是在基准流场建立过程中忽略了三维效应、低速效应、黏性效应以及头部/前缘的钝化效应,因此其高低速气动特性均有优化设计的空间.针对此问题,本文结合高保真RANS求解器、自由变形参数化方法、鲁棒的结构网格变形方法、离散伴随方法以及序列二次规划算法,发展了基于离散伴随的宽速域飞行器气动优化设计方法.基于上述方法,针对涡波一体乘波飞行器开展了兼顾低速与高超声速气动性能的三维整机气动优化设计研究,获得了宽速域优化构型并对其进行了流动机理分析.结果表明,相较于初始构型,宽速域优化构型可以将飞行器高超声速状态下升阻特性略微提升的同时,显著增强低速状态飞行器背风面的旋涡效应,进而使飞行器低速状态的升力和升阻比均提升10%以上,改善了涡波一体宽速域乘波飞行器的高低速气动性能.