基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

共轴刚性旋翼悬停状态地面效应气动特性

为了研究地面效应下共轴刚性旋翼的气动特性,建立了一套基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,采用运动嵌套网格模拟双旋翼的反转运动。地面采用无滑移边界条件,并对旋翼和地面附近的网格进行加密,以更好地捕捉旋翼的流场细节和尾迹特征。计算结果与Lynx尾桨试验结果进行对比,验证了所建立方法的有效性。对地面效应下共轴刚性旋翼的气动性能和流场进行分析,结果发现：相对于单独的上下旋翼而言,共轴旋翼地面效应下的拉力增益更大,这是由于上下旋翼桨叶表面的压强干扰受地面高压的影响而减弱;地面的干扰主要影响双旋翼尾迹的径向位置,对其轴向位置影响不大,上下旋翼尾迹在地面附近相互融合、分裂,形成复杂的桨尖涡尾迹;双旋翼在地效下的尾迹径向扩张半径比单旋翼大,这是由于双旋翼的径向射流速度更大;随着旋翼距地面高度的增加,双旋翼间的气动干扰强度逐渐恢复,因此下旋翼拉力增益的下降速度比上旋翼更大;共轴旋翼桨尖涡相对卷起高度和扩张半径均随离地高度增加而减小。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

尾迹对压气机转子性能影响的非定常数值模拟

采用三维粘性非定常数值模拟方法,研究了上游尾迹对轴流压气机转子性能及其尖部非定常流动的影响.结果表明,在一定情况下,上游静子尾迹与转子内部流动的非定常相互作用,有可能改善近失速点的气动性能,如转子压比和效率升高,工作范围增大.其原因主要为:上游静子尾迹使转子尖区一次泄漏涡强度减弱,减少了二次泄漏涡强度或抑制了二次泄漏涡的产生,最终导致尖区损失减少;此外,尾迹使尖区激波位置后移,改变了尖部弦向的负荷分布,最终导致压气机稳定工作范围增大.     

非定常预处理方法在倾转旋翼飞行器悬停状态气动干扰模拟中的应用

针对倾转旋翼飞行器的多种速度尺度并存的复杂流动问题,发展了能够统一求解同时包含不可压缩流动和可压缩流动的非定常预处理方法和一种考虑了流场物理信息传播方向的预处理远场边界条件。结合适用于倾转旋翼飞行器的复杂外形及大位移问题的非结构混合网格技术和重叠插值技术,建立了倾转旋翼飞行器的数值模拟方法。在此基础上,首先通过悬停算例验证了数值方法的有效性,确定了选择预处理截断参数的原则。然后探究了倾转旋翼飞行器直升机模式悬停状态的尾迹涡结构演化过程,研究了“喷泉效应”和“地面效应”对倾转旋翼气动性能的影响。计算结果表明,旋翼尾迹涡的耗散在不同方位角存在较大差异。“地面效应”抑制了旋翼尾迹的发展,提高了旋翼整个方位角区域的升力系数,而不仅限于机翼上方方位角区域。 “地面效应”作用强于“喷泉效应”作用,使旋翼相对于孤立旋翼拉力系数略微增大了5.26%,但同时整个倾转旋翼飞行器的总的向上的力只有孤立旋翼拉力的82.46%。结论表明预处理方法和预处理远场边界条件能够较好地模拟倾转旋翼飞行器悬停状态的干扰流场,兼顾了数值计算的耗散性和稳定性。气动干扰对倾转旋翼飞行器整体气动性能有重大影响。      

浸入式贴体网格边界方法的对转桨扇非定常流动数值模拟

应用所提出的浸入式贴体网格边界方法和开发的CFD程序,进行了对转桨扇跨声速全周叶片流场的非定常数值模拟.介绍了程序流程、网格生成策略和求解器原理.分析了前、后桨叶尖涡、尾缘涡的非定常干扰对叶片表面流场的影响.统计了对转桨扇的推力、功率和推进效率的非定常波动,并同NUMECA模拟的定常结果进行了对比,结果发现两种模拟得到...     

悬停状态倾转旋翼机非定常气动干扰研究

悬停状态旋翼/机翼机身干扰流场的高精度数值模拟对准确预估倾转旋翼飞行器气动性能具有十分重要的意义,开展该干扰流场的高精度数值模拟一直是直升机空气动力学领域的研究热点和难点之一。本文基于运动嵌套网格技术,建立一套适用于倾转旋翼非定常流场的CFD方法。采用雷诺平均Navier-Stokes(RANS)方程作为主控方程,湍流模型选用Spalart-Allmaras模型,时间推进上采用高效的隐式LU-SGS格式。在此基础上,首先开展某新型倾转旋翼无人机的旋翼/机翼机身非定常干扰流场数值模拟研究,得到了与文献相符的“喷泉效应”现象,计算结果显示干扰作用使总拉力损失16.6%。然后着重研究了不同襟副翼预置角对降低悬停状态下旋翼/机翼机身气动干扰作用的影响,研究结果显示45°襟副翼预置角效果最佳,使全机总拉力损失从原来的16.6%降到13.8%。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

某型倾转旋翼机的旋翼桨叶气动优化设计

倾转旋翼机的旋翼气动外形设计需要对其直升机模式和固定翼模式下的不同要求进行综合考虑。通过基于Kriging模型的多目标遗传算法建立一套适用于倾转旋翼桨叶气动外形优化设计方法,采用拉丁超立方抽样试验设计方法得到样本点,并建立Kriging模型替代费时的流动数值模拟。以最大化旋翼地面悬停拉力和高空巡航效率为目标,以地面悬停功率不增和巡航拉力不减为约束条件,进行倾转旋翼桨叶平面形状优化设计,并经过非定常数值模拟和风洞试验验证。结果表明:数值模拟结果和风洞试验结果吻合良好,优化结果满足设计指标。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

An investigation of tip vortices unsteady interaction for Fokker 29 propeller with swirl recovery vane

A numerical study was performed to explore the unsteady interaction between the upstream propeller and the downstream swirl recovery vane (SRV) by transient simulations. Much larger fluctuations of thrust coefficient were observed on the vane, which indicates that the varia-tions of the total efficiency depend mainly on the working performance of the stator. The harmonic loads of the decomposed unsteady blade-surface pressures show that the stator experiences about ten times higher of unsteadiness compared with the rotor. Notable changes appear at the vane lead-ing edge due to the potential disturbance as well as the sweeping effects from the wake of the upstream propeller, whereas more significant unsteadiness occurs at the stator tip region as a result of the interaction between the rotor/stator tip vortices. The visualization of vortex structures addresses that the rotor tip vortex has a dominant effect on the stator tip vortex since the latter one starts right at the impingement location on the vane top in this configuration. Furthermore, a longer and a shorter SRV were investigated based on the original case to explore different inter-action patterns for the rotor/stator tip vortices. Weaker effects have been observed as expected.     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

基于运动嵌套网格的前飞旋翼绕流N-S方程数值计算

 通过求解 Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间上采用中心平均的有限体积法进行离散,时间方向应用含隐式子迭代的双时间法推进求解。为了验证程序的正确性,数值计算了一有升力悬停流场,旋翼桨叶表面压力分布的计算值与实验值吻合很好。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

A high-order compact finite-difference scheme for large-eddy simulation of active flow control

The purpose of this article is to summarize a computational approach, which developed and matured over an extended period of time, and has been shown to be useful for performing large-eddy simulation (LES) of flows with active control. Because of the nature of active flow control, simulation of this class of problems typically cannot be carried out accurately by methods less sophisticated than LES. Active control flowfields are highly unsteady, and can be characterized by small-scale fluid structures which are produced by the control process, but may also be inherent in the original uncontrolled situation. The numerical scheme is predicated upon an implicit time-marching algorithm, and utilizes a high-order compact finite-difference approximation to represent spatial derivatives. Robustness of the scheme is maintained by employing a low-pass Pade-type nondispersive spatial filter, which also accounts for the fine-scale turbulent dissipation that otherwise is traditionally provided by an explicitly added subgrid-scale (SGS) stress model. Geometrically complex applications are accommodated by an overset grid technique, where spatial accuracy is preserved through use of high-order interpolation. Utility of the method is illustrated by specific computational examples, including suppression of acoustic resonance in supersonic cavity flow, leading-edge vortex control of a delta wing, efficiency enhancement of a transitional highly loaded low-pressure turbine blade, and separation control of a wall-mounted hump model. Control techniques represented in these examples are comprised of both steady and pulsed mass injection or removal, as well as plasma-based actuation. For each case, features of the flowfield are elucidated and the solutions are compared to the baseline situation where no control was enforced. Where available, comparisons are also made with experimental data.     

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

直升机前飞旋翼桨尖涡定常质量射流控制计算

针对直升机前飞状态开展了旋翼非定常桨尖涡的模拟,采用定常桨尖空气质量射流(TAMI)控制方式对旋翼桨尖涡进行了控制。建立了一个适用于前飞旋翼桨尖涡高精度捕捉和质量射流控制模拟的数值方法,在该方法中,采用有限体积法进行空间离散,应用5阶Roe-WENO(weighted essentially non-oscillatory)格式进行流场重构及控制面对流通量计算;采用双时间方法进行时间推进,伪时间步上应用隐式LU-SGS(lower upper symmetric Gauss-Seidal)格式;引入射流边界条件对质量射流进行模拟;采用运动嵌套网格方法对前飞旋翼桨尖的挥舞运动进行模拟,并对桨叶网格和背景网格进行针对性加密。基于所建立的方法对前飞状态旋翼非定常桨尖涡及其质量射流控制进行了模拟,计算结果表明:前飞状态下旋翼桨尖涡存在较大的前后差异,桨盘前侧的桨尖涡涡核强度远小于桨盘后侧;桨盘前侧旋翼桨尖涡的涡核强度很难由定常质量射流控制来减弱,而桨盘后侧的旋翼桨尖涡则可以通过定常质量射流得到有效控制。     

跨声速轴流压气机近失速状态的间隙泄漏流流动特性

为了研究间隙泄漏流以及泄漏涡自身的非定常性对轴流压气机的旋转失速的影响.对跨声速轴流压气机NASA转子37进行全三维定常、非定常的数值模拟,对比了最大效率和近失速工况的实验和数值模拟结果,定常计算所获得的总性能与试验结果符合良好.对于非定常计算,详细分析了NASA转子37近失速工况下流场结果,揭示了NASA转子37在近失速工况点,间隙泄漏流存在较明显的非定常性,这种非定常性表现为间隙泄漏流激波干涉引起间隙泄漏涡的周期性破碎.