非交错网格上不可压流动中的近似投影方法

对于压力增量形式的投影方法,在非交错网格上,仅仅采用近似投影过程,通常不足以抑制压力场的数值振荡.本文通过选取合适的投影变量,构造了一种可以在非交错网格上避免压力振荡的新方法.该方法在不增加计算量的前提下,让前一时刻的压力参与投影过程,有效地割断了前一时刻中可能存在的压力数值振荡向下一时刻传递的途径,避免了数值振荡的积累,提高了数值方法的稳定性.此外,我们在投影变量中考虑了"几何网格雷诺数"的影响,拓宽了本方法的适用范围.数值试验验证了本文提出的近似投影方法是可行的.     

旋翼翼型定常-非定常特性综合优化设计新方法

旋翼翼型对直升机旋翼及全机气动特性有至关重要的影响。结合直升机飞行特性及旋翼工作气动环境,本文提出了一种新的旋翼翼型设计理念,即定常优化设计同非定常设计相结合的方法来开展直升机旋翼专用翼型的优化设计。首先,建立旋翼翼型非定常气动特性的高精度CFD求解方法,以获得旋翼翼型在相应状态下的升力、阻力和力矩等气动参数。其次,针对旋翼翼型优化状态的特点,对于翼型静态优化及动态失速优化分别采用遗传算法和序列二次规划算法。在上述方法建立基础上,首先针对定常状态下的原始翼型(SC1095)进行优化设计,获得一个满足设计要求的静态优化翼型,进一步着重对该翼型在非定常状态下进行优化设计,成功地得到一种新翼型,并具有非常规外观。结果表明,在设计状态下,新翼型在保持良好定常气动特性的同时,明显减弱了动态失速状态下的分离涡,从而显著改善了翼型的动态失速特性。     

基于代理模型的飞翼多目标气动优化设计

针对飞翼布局飞行器,采用雷诺平均N-S方程（ RANS）计算流场,使用基于代理模型的多目标优化方法进行了同时考虑起飞性能和巡航性能的多点多目标气动优化设计。在设计过程中,将飞翼的平面形状、剖面形状及扭转角同时作为设计变量（共58个设计变量）,将提高起飞时的升力系数和提高巡航升阻比为设计目标,以起飞状态和巡航状态的力矩系数作为气动约束,并以飞翼平面面积不减和剖面厚度不减作为几何约束。通过采用基于Kriging模型的多目标优化方法,以较小的计算花费得到了较好的Perato前沿。取Pareto前沿中一个最优解与基准外形的性能进行了对比,结果显示,优化外形的性能较基准外形的气动性能得到全面大幅提高且所有约束得到严格满足。     

盒式机翼布局气动特性研究

采用基于S-A一方程湍流模型的雷诺平均法求解N-S方程,研究了盒式机翼布局前后翼几何参数对气动特性的影响规律,为盒式布局飞机气动设计提供参考。结果显示:随着前后翼垂直距离的增大,盒翼升力和升阻比均增大;增大前后翼的纵向距离,盒翼升力增加;盒翼气动性能对两翼纵向距离的敏感程度不及对两翼的垂直距离。     

超声速进气道与弹体一体化外形研究

通过数值求解N S方程和风洞测力试验,研究了超声速进气道的布局、剖面形状及几何参数对全弹气动力和进气道内流的影响。结果表明:细长旋成体弹身背风面和侧边不宜安放进气道;矩形进气道能够产生涡升力,升阻比明显高于半圆形进气道,而且内部流动接近二维流动,流场畸变及流动损失情况均比半圆形进气道要好;矩形进气道横截面积增加,对升力系数影响不大,但阻力系数增加明显;其高宽比增加,升力系数增大,阻力系数减小。     

基于多点约束的大展弦比机翼静气动弹性计算

基于雷诺平均N-S方程和多块结构网格技术计算大展弦比机翼的气动力.机翼采用梁模型结构,利用有限元方法计算结构变形.建立基于多点约束(MPC,Multi-Point Constrain)的气动、结构数据的双向传递方法.采用MPC数据插值方法快捷地实现自主研发的计算流体力学(CFD,Computational Fluid Dynamics)程序和NASTRAN计算结构力学(CSD,Computational Structural Dynamics)软件的耦合计算,发展出可以考虑结构非线性的静气动弹性CFD/CSD耦合计算方法,该方法比只考虑单方向变形的柔度法更精确.开展了某大展弦比机翼的静气动弹性数值计算分析,结果显示只考虑机翼纵向变形而忽略展向变形对计算结果会有一定影响,而该算例几何非线性影响较小.     

二股流喷射控制推力矢量数值计算研究

利用求解雷诺平均的N-S方程的数值模拟方法,对一种新型的推力矢量喷管———基于次流喷射控制的二维收缩-扩张喷管(2DCD)———的推力矢量性能进行了研究。根据数值模拟结果分析了外流马赫数、喷管压强比和次流总压与主流总压之比对矢量偏角的影响。     

后向台阶横向喷流干扰流场的数值研究

 用有限差分法求解层流全N-S方程,对超音速后向台阶绕流与台阶后横向喷流干扰流场进行了数值模拟。不带喷流的超音速后向台阶绕流数值模拟结果与实验及其他计算结果吻合较好。最后给出横向喷流压力比p_j/p_∞=3.033～7.033和喷口宽度w=0.1h～0.2h(h为台阶高度)时的数值模拟结果。得到了合理的流场结构;并可看到喷流压力比和喷口宽度对流场特性有显著影响,在计算参数范围内,底部压力与系数Q(=wp_j/hp_∞)几乎成线性关系。     

吸气式内外流一体化飞行器动导数数值模拟

动导数是飞行器轨道及姿态控制系统设计时的重要参数,对飞行器开环系统受到扰动时振荡的敛散特性起重要作用。在标模验证的基础上开展吸气式内外流一体化飞行器WR-A动态特性分析,采用N-S方程模拟强迫简谐振动的非定常流场,获得了飞行器的直接阻尼导数、加速度导数和旋转导数,并对WR-A外形大振幅振动时的进气道性能参数进行了分析。研究表明,在某些情况下反映流动时滞效应的加速度导数占直接阻尼导数的比例最高可以达到40%以上,而在某些情况与直接阻尼导数符号相反。此外,在动导数研究基础上,采用参数化的动导数概念初步开展了非定常气动参数建模研究。针对WR-A大振幅强迫振动的研究显示,采用所建立的模型预测的气动参数与完全非定常计算吻合较好。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。