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

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

基于非结构嵌套网格的涵道螺旋桨数值模拟

基于非结构嵌套网格方法,通过求解可压缩Euler方程对涵道螺旋桨系统的非定常流场进行了数值模拟,并研究了桨尖间隙、螺旋桨位置、涵道尾扩张角等参数的影响规律.将计算区域分为相互重叠的旋转区域和静止区域,旋转区域包含螺旋桨和旋转的桨毂,静止区域包含涵道和部分静止的桨毂,采用嵌套网格算法真实地模拟了螺旋桨和涵道之间的相对旋转运动,得到了细致的螺旋桨流场.计算值和实验测量值的对比结果相一致.参数影响研究表明,减小桨尖间隙、在一定范围内后移螺旋桨或增加尾扩张角,均可以提高涵道螺旋桨的气动性能.     

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

基于伴随方程方法的非结构网格自适应技术及应用

对基于伴随方程的网格自适应技术的基本原理进行了研究,构造并实现了基于流场伴随方程的非结构网格自适应探测器,利用局部加密网格的办法,建立了基于Euler方程的无粘伴随自适应能力。对于粘性流动问题,提出首先应用无粘计算和伴随自适应技术从基准网格自动生成具有较好网格分辨率的非结构自适应网格,然后采用层推进法自动生成粘性计算网格的求解策略。采用RAE2822二维翼型和ONERA-M6三维机翼对所建立的无粘伴随自适应方法进行了验证,并应用所提出的方法和策略,对第二届国际涡流试验项目(VFE-2)中的65°尖前缘三角翼大攻角涡流场进行了数值模拟。对比计算结果表明,所建立的非结构网格伴随自适应方法是可行的和有效的。     

应用重迭网格技术求解复杂组合体无粘流场

本文数值求解了捆绑火箭这一复杂组合体的无粘干扰流场。该计算利用重迭网格(Overlapping Grid)技术求解了分离系数的Euler方程,获得了复杂组合体相互干扰的重要流场特性。通过与风洞试验的比较,计算结果与实验有较好的一致性。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

涵道螺旋桨与孤立螺旋桨气动特性的 数值模拟对比

基于非结构动态嵌套网格方法,对涵道螺旋桨与孤立螺旋桨的气动特性进行了非定常Euler方程数值模拟对比研究和分析.将计算区域分成旋转区与静止区两个子域,生成相互重叠的嵌套网格,有效地解决了螺旋桨与涵道以及静止部分桨毂之间的相对旋转问题.计算结果与实验结果相一致,验证了计算方法的正确性.涵道入口前缘形成的负压区是产生涵道附加拉力的原因所在.对比计算结果表明:与孤立螺旋桨相比,涵道的存在改善了螺旋桨桨尖区域的绕流特性,减小了桨尖损失;相同转速情况下,涵道螺旋桨产生更大的拉力,而所需功率略小,涵道螺旋桨系统具有更高的气动效率.      

全机复杂流场网格生成及N-S方程数值模拟

全机复杂流场网格生成对计算结果起着举足轻重的作用。本文首先通过合理的拓扑结构,分块对接的方法对复杂流场进行了网格划分,然后采用N-S方程有限体积法及分区对接法对翼—身—尾全机流场进行数值模拟。计算结果验证了网格生成的有效性和计算结果的合理性。     

螺旋桨非定常滑流数值模拟

运用非结构动态嵌套网格技术,通过求解三维非定常Euler方程,对螺旋桨飞机的螺旋桨滑流进行了真实数值模拟.动态嵌套网格方法通过建立运动区子域网格来模拟螺旋桨的旋转运动,静止的机身包含在静止区网格中,有效解决了螺旋桨滑流模拟中的相对运动问题.采用双时间方法求解控制方程,在每一个物理时间步上运用经典四步Runge-Kutt...     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

利用转移矩阵进行姿态控制

传统欧拉角控制方法在处理大姿态角机动时往往出现奇异和交联问题。针对这一问题对欧拉角的跳变现象进行了讨论,并从飞行器实际姿态与需求姿态之间的空间关系入手,给出了一种消除欧拉角跳变现象的有效方法－转移矩阵法,仿真计算表明,该方法在大型航天器的姿态控制中可以有效地消除欧拉角奇异、交联和跳跃现象,在欧拉角跳变的情况下仍然能够实现对飞行器的控制。     

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

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

基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

Unstructured Euler solutions for missile aerodynamics

The purpose of this paper is to validate an unstructured 3D Euler flow solver, USER3D, against the experimental data and to compare it with a classical, structured, Euler solver, FLU3M of ONERA, on two different missile geometries at a Mach number of 2 and at various angles of attack up to 20 ° . The first geometry is a conventional missile with an ogive nose, a cylindrical body and four straight tail fins. The second geometry is a unconventional missile with a lenticular body and no tail fins. The present results are found to be in good agreement with these available results and the discrepancies observed between them are explained in detail. The present unstructured Euler code has proved to be accurate, fast and reliable for the aerodynamic design of conventional and unconventional missiles.     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

双欧法与四元数法的应用比较

对解决欧拉方程奇异性的双欧拉和四元数法进行了对比研究，四元数法从理论上讲比较完美，但实际应用中存在较大的累积计算误差，从而影响计算精度；双欧法利用正、反欧拉方程间精华区倒挂关系进行分区交替运算，把精华区扩展到全域，不仅根除了奇异性，而且计算误差小。因此，对于解决欧拉方程奇异性来讲，双欧法要优于四元数法。     

旋翼表面非定常压力脉动计算的三维自由尾迹非定常面元法

利用最新发展的高分辨率三维自由尾迹非定常面元方法计算多叶片直升机旋翼的三维非定常气动力，将风洞实验数据与旋翼在悬停和出现严重桨涡干扰的低速前飞下滑状态下的计算结果进行了比较，在悬停时，还与欧拉方法的结果进行了比较，本方法与欧拉方法的结果以及风洞实验结果都吻合很好。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

Vega Prime实时视景仿真中实体姿态角研究与应用

实时视景仿真中准确描述实体的姿态十分重要,欧拉角就是一种描述实体姿态的有效方法,它广泛应用于各种系统,出于对各自问题解决的需要,不同系统对于欧拉角的描述不尽相同.实时视景仿真Vega Prime软件即采用欧拉角HPR来描述仿真实体的姿态,从欧拉角HPR的含义、如何将实体的姿态调整为给定姿态(h,p,r)和欧拉角HPR的求取等三个方面进行了深入研究,最后以实例的形式分析了欧拉角的求取过程和求取方法.     

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.