三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好     

自适应直角切割网格民机增升装置绕流数值模拟

提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明所述方法的正确性.     

基于直接力控制的阵风响应及阵风减缓研究

考察某无人机模型的阵风气动响应特性和基于直接力操纵的阵风减缓效果.在求解非定常Euler方程时引入“网格速度”方法模拟阵风边界条件,通过动态嵌套网格技术实现舵面运动.首先对NACA0006翼型迎角阶跃型阵风的气动力响应进行了验证计算,计算结果与理论结果、文献计算结果吻合良好.在此基础上对某无人机模型在迎角阶跃型、One-minus-cosine型阵风作用下的气动力响应过程进行了数值模拟,并分析评估了不同舵面运动方式对飞机气动性能的影响.最终研究比较了不同舵面运动方式下阵风气动响应的减缓效果.结果表明:通过合理设计舵面运动,可以有效达到抑制阵风引起的非定常气动干扰的目的.     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

鸭式导弹舵面后缘后掠角对滚转影响数值研究

鸭式布局导弹的滚转特性是抑制其性能提升的重要约束.采用数值模拟方法对鸭式导弹滚转特性进行研究,分析舵面后缘后掠角对其影响.以三维可压缩雷诺平均N-S方程为主控方程来模拟流场特性,湍流模型采用SST的κ-ω二方程,基于气动特性和流场特性分析研究了固定尾翼鸭式气动布局导弹的横滚特性.仿真结果表明,鸭式布局导弹的舵面后缘后掠...     

CFD在概念-初步设计阶段三维气动外形优化中的应用

在飞行器概念-初步设计阶段,建立基于CFD的气动优化链对于提高优化计算的效率具有较好的工程应用价值。使用德国宇航院开发的CPACS数据格式给出飞行器平面形状,结合NURBS翼型参数化方法对飞行器几何外形进行参数化;自动生成计算网格并求解Euler方程数值模拟流场以评估参数化气动外形的气动特性,进而构建响应面模型;使用SQP梯度算法搜索响应面模型以获取满足约束的最优解。以Onera M6机翼为例,对该优化链进行验证。结果表明:在满足约束的条件下,基于CFD的气动优化链能够成功地进行气动外形优化。     

螺旋桨非定常滑流的高效数值模拟研究

基于自主研发的大型"亚跨超CFD软件平台"(TRIP3.0),开发了针对螺旋桨飞机的旋转机械非定常求解模块。控制方程采用曲线坐标系下的非定常Euler方程,综合运用动态拼接结构网格技术和双时间步推进方法,重点模拟了双发螺旋桨飞机在螺旋桨转动时的非定常滑流效应。为了保证计算模拟的高效性,引入多重网格技术和大规模并行计算技术。通过模拟某单独螺旋桨验证算例,验证了该方法的可行性。模拟了不同迎角和螺旋桨转速下全机的气动特性,分析了滑流对流场的影响情况,结果表明:螺旋桨滑流对全机流场和气动特性影响明显,并使得升力系数增大。     

带控制舵弹体气动特性分析

本文用数值模拟方法研究了带控制舵弹体的气动特性。该问题外形复杂，在控制舵和弹身之间会出现明显的缝隙。为了解决此类外形的气动力计算问题，我们发展了结构网格和非结构网格的混合网格技术。这样一方面充分发挥了贴体结构网格质量好、效率高的优点，另一方面充分发挥了非结构网格几何灵活性好、易于模拟复杂外形的优点。数值计算与风洞实验相当吻合，说明本方法具有良好的精度。最后对两种外形进行大量计算并进行气动特性分析，给出了静稳定特性和配平特性的研究结果。     

三维复杂无粘流场的自适应八叉树结构直角网格算法模拟

运用自适应八叉树结构直角网格对三维复杂流场进行了划分 ,并采用点面判别法进行外法向量的判断 ,使与之关联的与物面相交网格的体积计算更为简便准确。在此基础上 ,采用以中心差分为基础的 Jameson的有限体积法进行 Euler方程求解 ,在通量计算过程中 ,采用以面为基础的算法 ,减少了计算工作量。对多个外挂、翼身组合体及全机等复杂外形绕流流场进行了模拟 ,均取得了较好的效果     

外吹式襟翼动力增升绕流数值分析

基于Euler方程来数值模拟运输机吹气襟翼的流动机理,通过优化发动机及襟翼系统的安装,以此评估其动力增升效果.计算结果分析表明本文方法对动力增升概念设计有很好的适用性.     

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

考虑喷流效应的机载导弹发射及气动干扰数值模拟

研究比较了机载导弹发射时是否考虑导弹尾喷流情况下导弹运动特性及机-弹气动干扰影响.通过动态结构嵌套网格技术实现了导弹运动的模拟,结合刚体六自由度运动方程求解三维非定常Euler方程得到流场信息.导弹推力由发动机燃烧室喷流计算得到;当不考虑喷流效应时,推力通过直接在导弹尾部给定非定常作用力实现.应用此方法模拟了基于类"全...     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

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

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

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。     

1

Ｄｉｒｅｃｔｏｐｅｒａｔｉｎｇｃｏｓｔｓｆｏｒｔｒａｎｓｐｏｒｔａｉｒｃｒａｆｔｈａｖｅｉｎｃｒｅａｓｅｄｒａｐｉｄｌｙｓｉｎｃｅｔｈｅｏｉｌｃｒｉｓｉｓｏｆｔｈｅ１９７０ｓ．Ｔｈｅｐｏｔｅｎｔｉａｌｆｏｒｌａｒｇｅｆｕｅｌｓａｖｉｎｇｓａｎｄｔｈｅ...     

Numerical studies of static aeroelastic effects on grid fin aerodynamic performances

The grid fin is an unconventional control surface used on missiles and rockets. Although aerodynamics of grid fin has been studied by many researchers, few considers the aeroelastic effects.In this paper, the static aeroelastic simulations are performed by the coupled viscous computational fluid dynamics with structural flexibility method in transonic and supersonic regimes. The developed coupling strategy including fluid–structure interpolation and volume mesh motion schemes is based on radial basis functions. Results are presented for a vertical and a horizontal grid fin mounted on a body. Horizontal fin results show that the deformed fin is swept backward and the axial force is increased. The deformations also induce the movement of center of pressure, causing the reduction and reversal in hinge moment for the transonic flow and the supersonic flow,respectively. For the vertical fin, the local effective incidences are increased due to the deformations so that the deformed normal force is greater than the original one. At high angles of attack, both the deformed and original normal forces experience a sudden reduction due to the interference of leeward separated vortices on the fin. Additionally, the increment in axial force is shown to correlate strongly with the increment in the square of normal force.     

非结构动网格在三维可动边界问题中的应用研究

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(Arb itrary Lagrangian-Eu lerian ALE)描述的Eu ler方程,模拟了做谐和振动的刚性无后掠机翼及柔性M6机翼的跨音速绕流,计算结果与实验结果吻合良好。