低速预处理方法在悬停旋翼模拟中的应用

在可压缩方程的求解当中引入时间导数低速预处理技术,构造了能够同时处理低速和高速流场的统一数值方法,采用该方法对Caradonna-Tung旋翼的亚声速悬停状态进行了模拟,其中方程采用在旋转坐标系下求解,网格采用重叠网格进行模拟。数值结果表明,在加入低速预处理方法以后,流场的收敛速度得到了明显提高,更为适合对下洗速度场等细节的模拟;另外,在加入低速预处理以后,方程的各个特征值都处在同一量级,避开了因特征值差异带来的数值耗散过大的问题,提高了数值计算的精度。

基于非定常低速预处理和DES的三角翼数值模拟

发展出一套基于低速预处理技术采用SA-DES和SST-DES模型,对三维低速非定常流动进行数值模拟的方法。采用非结构混合网格有限体积方法求解非定常流场,时间离散采用基于LU-SGS隐式格式的双时间步长方法。对非定常低速预处理方法进行了推导和分析,并应用该方法对大攻角三角翼的非定常流场进行了数值模拟,对不同站位压强分布与试验值和参考文献值进行了对比,吻合较好。

基于低速预处理的风力机翼型外形优化设计

针对经典的S809翼型,耦合基于低速预处理的流场求解方法和序列二次规划方法,开展针对翼型升阻比的翼型气动外形优化设计研究。优化结果显示优化翼型具有较大的翼型前缘半径和较平坦的上表面。数值计算结果表明,优化翼型在设计点1的状态下升阻比提高43.3%,在设计点2的状态下升阻比提高48.9%。进一步数值验证表明,优化翼型在雷诺数为5.0×10^(5)状态下的最大升力系数从S809翼型的1.140增大到1.297,在雷诺数为1.0×10^(6)状态下的最大升力系数从1.236增大到1.418。在优化翼型的基础上,开展翼型气动外形人工修型研究,数值模拟表明修型翼型能更好地消除气流分离,从而进一步增大翼型升力系数、减小翼型阻力系数。

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出。本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证。

旋翼/机身气动干扰的数值模拟

通过求解粘性非定常可压缩N-S方程,发展了旋翼/机身气动干扰的数值模拟方法,主要技术措施包括:(1)通过ADT数据结构,建立了高效的多块对接网格动态重叠技术;(2)在可压缩方程求解中引入低速预处理方法,克服低速流场求解时的刚性问题;(3)加入了多重网格方法,提高收敛效率;(4)湍流模型方程求解需要反复计算壁面最近距离,采用Wigton的优化算法很大程度上提高了计算速度。利用建立的计算方法对佐治亚大学(GIT)旋翼/机身干扰模型进行了数值模拟,分别比较了机身表面的平均压力分布、旋翼的下洗流场以及非定常压力脉动,对干扰现象及其流动机理进行了分析。

基于混合笛卡尔网格的风力机问题数值模拟

基于混合笛卡尔网格方法,对S809翼型和Phase VI叶片的绕流问题进行数值模拟研究。计算网格物面附近生成贴体结构网格,其余部分使用笛卡尔网格进行填充,两者之间通过查找"贡献单元"的方法来进行流场信息的传递。计算方法采用低速预处理方法消除可压缩方程在处理低速问题时的刚性,同时结合旋转坐标系方法考虑风力机旋转效应,进行风力机低速问题数值模拟研究。另外,在Phase VI叶片的计算中,通过笛卡尔网格的局部加密来捕捉尾涡。通过与实验数据及他人数值计算的结果进行对比分析,结果表明该文开发的这套基于混合笛卡尔网格的低速求解器能较为准确模拟风力机叶片轴流状态的流场,并能准确捕捉叶尖涡和叶根涡在风力机下游的发展。