基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

倾转旋翼机巡航状态旋翼滑流影响

倾转旋翼机由于需要兼顾垂直起降和高速平飞2种典型工况下的动力需求，采用大直径旋翼作为推进装置会使机翼大部分处于旋翼滑流区内，这与常规螺旋桨飞机存在较大差异。为评估不同数值计算方法并研究旋翼滑流对倾转旋翼机气动特性的影响，针对选取两叶旋翼的某倾转旋翼机方案，利用激励盘模型、多参考系（MRF）模型、滑移网格模型分别进行了巡航状态下旋翼滑流对全机气动特性影响的数值模拟研究。结果表明:相对于无滑流状态，滑流定常影响使全机阻力增大，最大升阻比降低了7.5%，尾翼产生的升力增大，纵向静稳定度增加了17.1%，全机低头力矩增大；当迎角较小时，滑流虽然改变了机翼表面的升力分布，但是全机升力变化不大；滑流非定常影响会使全机气动特性产生周期性波动，升力系数波动幅度为9.0%，阻力系数波动幅度为10.8%，并且随着迎角的增大，波动幅度也越大。      

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

螺旋桨滑流对发动机进气道气动性能的影响

针对螺旋桨/短舱/进气道一体化构型,采用多重参考系(MRF)法和滑移网格非定常法进行了数值模拟研究,实现了螺旋桨/短舱/进气道一体化内外流耦合流场仿真,研究了不同工况下螺旋桨滑流对进气道气动性能的影响,并与无桨构型的进气道性能进行了对比,最后通过数值分析解释了螺旋桨桨盘后总压上升的原因。结果表明:多重参考系法和滑移网格非定常法对进气道出口的时均参数的计算结果吻合较好,多重参考系法可用于螺旋桨及进气道的气动性能参数计算;螺旋桨做功使进气道总压恢复系数增高,螺旋桨转速越大,总压恢复系数的增量就越大,且螺旋桨滑流使进气道出口畸变指数明显增大;通过数值分析推导的桨盘前后总压关系公式与仿真结果的误差均在3%以内。      

eLoran系统ASF网格应用算法研究

ASF网格是提高增强型罗兰（Enhanced Long Range Navigation，eLoran）系统精度的重要方法。根据所在网格四个顶点的附加二次时延（Added Secondary Factor，ASF）值通过网格应用算法得到待测试点的ASF值，一定程度上比公式计算的ASF值更准确。在计算过程中，用户四个顶点ASF值的测量误差会传递给用户，同时内插算法本身也会引入误差，此时内插算法的选择尤为重要。本文通过仿真、分析比较几种常用的内插算法，得出结论：反距离插值算法引入的误差最大，双线应插值算法引入的误差最小，而且误差呈现一定层次感，网格从内到外误差逐渐减小。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

三角网格面自动铺丝定角度路径规划算法

研究三角网格曲面上的路径规划问题,以局部测地线为补充,提出了改进的定角度初始路径规划算法。该算法避免了传统参数化求解过程的无解和多解情况,尤其适用于曲率变化较大的复杂曲面和拼接曲面。采用等距覆盖算法来获得均匀的满铺路径,并分别提出基于路径样条线末端延伸和芯模曲面整体延拓的两种边界处理思路,以解决偏移路径无法到达曲面边界的问题。以某机翼模型为算例对上述算法进行验证,利用计算机辅助三维交互应用(Computer aided three-dimensional interactive application,CATIA)的逆向工程模块实现了完整的算法可视化,建立了三维建模软件与路径规划算法的系统关系。验证表明本文提出的算法可靠、系统、适应性强,具有工程应用价值。     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。