包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

大迎角涡流的同步迭代算法

 基于多重线涡模型(MLVM),本文建立了一种计算机时少、收敛性强的适用于大迎角涡流势流数值模拟的迭代算法。采用本文算法,首次得到了迎角直到60°的旋转体涡流的收敛解。计算出的空气动力数据与实验结果吻合。     

翼型跨音速小扰动势流-边界层计算

翼型绕流的外部无粘流,采用跨音速小扰动速势方程和改进小扰动速势方程对三种差分格式进行对比计算:(1)Murman-Cole非守恒式;(2)守恒式;(3)Engquist-Osher式。在相同的初值条件下,它们的计算稳定性、收敛性和收敛解基本相同,但激波位置(1)在前,(2)在后,(3)在中间。对人工时间阻尼项εφ_(xt)的作用进行了线化分析,所得结论与数值实验结果相符。二级精度格式计算表明:采用这种方案能在几乎不增加计算机时的前提下,大大提高差分计算精度。本文最后采用边界层积分方程法,计算翼型边界层,并与势流进行迭代计算。在弱激波条件下,迭代收敛解与实验值接近。     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

基于流水避石原理的无人机三维航路规划方法

借鉴自然界流水避石现象,提出一种基于流体计算的无人机(UAV)三维(3D)航路规划方法.首先介绍了球心位于坐标原点时,球形障碍三维绕流问题的解析解.之后采用旋转平移矩阵与流线数据叠加方法生成了任意位置多障碍同时存在的三维流线.为验证解析解的有效性同时给出该方法基于数值模拟的计算过程,对适合无人机三维航路规划的流体模型和数值求解方法进行了分析,并给出了通过数值模拟求解航路的方法.最后,根据无人机机动约束对流线进行处理得到可飞航路,将航路长度、纵向和横侧向机动次数作为子目标函数对航路进行综合评价.仿真结果表明:解析法航路规划中,圆球障碍的地形建模简单计算量小,航路集中在由起点至终点的航路带间;数值法航路规划适合障碍分布复杂的地形,航路分布于规划空间中.这两种方法的航路平滑,能够满足无人机飞行约束,航路具有绕流意义的最优性,可以避免势场法的局部极小问题,并且可以提供多条备选航路.     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。     

一种边界非协调数值方法

研究了一种边界非协调的数值方法。该方法的提出是为了能够方便的求解复杂的运动边界绕流问题。空间上采用Galerkin有限体积空间离散格式。数值离散中,网格点分为内部计算点、IB点和外点三种类型。其中临近物面的IB点在嵌入物面边界条件时起到了作用,其流动变量值则是通过流场内解的线性插值获得。时间推进上采用结合人工压缩项的双时间步推进格式。静止和振荡圆柱不可压绕流的数值实验初步验证了本文方法的可靠性和准确性。     

开放式气动力数值模拟系统研究

气动力数值模拟系统是ＣＦＤ流场解算技术、网格生成技术、数据可视化技术和网络技术相结合的产物。描述了一种运行于计算机网格环境下多任务、多用户、交互式的数值模拟软件环境－ 气动力数值模拟系统,该环境能计算模拟亚、跨、超音速绕流的气动力,也能方便地集成模拟不同复杂气动布局飞机流场的解算程序。介绍了其网格生成、数据可视化、系统接口等核心部分的结构和功能。     

Parametric analyses on dynamic stall control of rotor airfoil via synthetic jet

The effects of synthetic jet control on unsteady dynamic stall over rotor airfoil are investigated numerically. A moving-embedded grid method and an Unsteady Reynolds Averaged Navier-Stokes(URANS) solver coupled with k-x Shear Stress Transport(SST) turbulence model are established for predicting the complex flowfields of oscillatory airfoil under jet control. Additionally, a velocity boundary condition modeled by sinusoidal function has been developed to fulfill the perturbation effect of periodic jet. The validity of present CFD method is evaluated by comparisons of the calculated results of baseline dynamic stall case for rotor airfoil and jet control case for VR-7 B airfoil with experimental data. Then, parametric analyses are conducted emphatically for an OA212 rotor airfoil to investigate the effects of jet control parameters(jet location, dimensionless frequency, momentum coefficient, jet angle, jet type and dual-jet) on dynamic stall characteristics of rotor airfoil. It is demonstrated by the calculated results that efficiency of jet control could be improved with specific momentum coefficient and jet angle when the jet is located near separation point of rotor airfoil. Furthermore, the dual-jet could improve control efficiency more obviously on dynamic stall of rotor airfoil with respect to the unique jet, and the influence laws of dual-jet's angles and momentum coefficients on control effects are similar to those of the unique jet. Finally,unsteady aerodynamic characteristics of rotor via synthetic jet which is located on the upper surface of rotor blade in forward flight are calculated, and as a result, the aerodynamic characteristics of rotor are improved compared with the baseline. The results indicate that synthetic jet has the capability in improving aerodynamic characteristics of rotor.     

地面对带后缘分离的翼型绕流影响的计算

本文提出了用在地面上分布奇点来模拟地面的方法。计算结果与精确解和镜像法结果很吻合。     

基于LBM-LES方法的多段翼型流场数值模拟

运用Lattice Boltzmann(LB)方法与Large Eddy Simulation(LES)结合的方法数值模拟研究了二维翼型绕流流场,其中,粒子速度模型采用D2Q9模型,松弛时间中对应的湍流粘性系数由Smagorinsky模型计算,网格生成方法采用无网格技术。对NACA0012翼型进行数值模拟,验证其有效性;为进一步研究其模拟复杂物体绕流的能力,以三段翼型30P30N为研究对象进行数值研究,计算出了给定攻角条件下的气动特性,并与实验值进行比较。计算结果表明,LBM-LES方法的计算结果能够与实验值一致,是一种可行、有效且较新的数值模拟方法。     

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

基于CST参数化方法的翼型快速设计

为弥补传统飞机翼型设计周期长、代价高的缺点,将CST翼型参数化方法,与机器学习中的高斯过程回归方法相结合,通过对已有的翼型数据的学习,实现对未知翼型气动性能或者外形数据的快速准确预测。选取一组NACA四位族翼型,获得其CST参数描述数据,并分别计算其在一定条件下的升力系数、阻力系数和压力分布数据。利用这些数据对高斯过程回归模型进行训练,实现了翼型的快速正设计以及反设计系统。并将实验结果与采用NACA翼型参数表示法得到的预测结果进行了对比。实验结果表明,基于CST参数化方法的翼型快速设计准确度高、速度快,具有很大的应用价值。     

基于CFD方法的机翼阵风响应研究

通过引入＂网格速度＂方法模拟阵风条件,求解非定常欧拉方程实现了不同展弦比0012平直机翼阵风响应的数值模拟。首先采用该方法对NACA0006翼型迎角阶跃型阵风的气动力响应进行了计算,计算结果与文献结果、理论值吻合良好。进一步对展弦比分别为5、10的0012平直机翼在迎角阶跃型、One-Minus-Cosine型阵风作用下的气动力响应过程进行了模拟分析。研究结果表明,当展弦比增大时,阵风作用下机翼的升力系数响应会增加,机翼翼根部位气动响应幅值大于翼尖部位响应特性。     

翼型结冰过程数值模拟验算与分析

应用FENSAP-ICE结冰计算软件,对NACA0012翼型进行了流动特性、水滴撞击特性以及冰型生成过程的计算;同时,对结冰前后的翼型进行气动力特性计算对比分析,其中包括升力特性对比、阻力特性对比、流场细节分析以及压力系数分布对比。计算结果表明：翼型前缘结冰后,导致翼型前缘气流提前分离,最大升力系数、失速攻角大幅减小,...     

结冰对翼型流场影响的研究

本文以NACA23012翼型为研究对象,采用CFD软件模拟了结冰对翼型流场的影响。结果表明：结冰对翼型的气动性能有着明显的影响,尤其是在大攻角时,升力系数显著下降,而阻力系数则明显增加;对于光滑翼型和结冰翼型．计算值与试验值在‘定攻角范围内吻合良好：马赫数对两种翼型的升力系数和阻力系数随攻角的变化影响均比较小,雷诺数也一样。     

海上湿气对翼型气动特性的影响初探

研究高湿度空气对风力机翼型气动性能影响对海上风力机翼型优化、提高发电效率具有重要意义。首先采用组分输运模型数值计算湿空气条件下翼型升阻力特性,其次在组分输运模型基础上,考虑湿空气中水蒸气的凝结,建立用户自定义函数定义凝结体积源项,求解输运方程。对湿空气流过翼型表面气动性能的初步研究发现：升阻力系数随湿度变化而改变,不考虑凝结时,升阻力系数随湿度增大而减小,当发生凝结时,升力系数增加。翼型表面前缘点附近凝结流量最大,翼型表面凝结流量随湿度增大而增大,温度越高,凝结流量越大。     

基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。