基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

多面体网格在静气动弹性计算中的应用

介绍了多面体网格技术在静气动弹性问题中的应用.基于M6算例和径向基函数插值法,探究了多面体网格的CFD计算效率和变形能力,使用AGARD445.6机翼进行了静气动弹性数值计算.结果表明,多面体网格具有计算效率高、变形能力强、计算结果准确等优点,可适用于静气动弹性问题的数值计算.     

基于新型降阶模型的网格加筋圆柱壳频率分析

针对网格加筋圆柱壳频率分析计算量大的问题,提出一种基于多项式及梁单元形函数的模型降阶方法,即运用多项式及梁单元形函数将复杂结构有限元的节点位移转化为主节点位移,实现模型降阶。并通过控制多项式阶数及梁单元个数来调整降阶模型的分析精度,而降阶模型的计算效率较精细模型有显著提高。以网格加筋圆柱壳为例,对比本文提出的模型降阶方法与3-D实体模型及其他模型降阶方法的频率分析结果,结果表明,提出的模型降阶方法不仅能捕捉到结构的整体模态,还能反映出结构部分局部模态,适用范围广,能够为工程结构设计提供简单有效的计算模型。     

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

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

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

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

桨间距对某型共轴对转螺旋桨气动性能的影响

为掌握共轴对转螺旋桨的桨间气动干扰规律，降低桨间气动干扰强度，提升共轴对转螺旋桨的气动性能，基于非定常雷诺平均Navier-Stocks方程耦合湍流模型的计算方法，并使用了滑移网格技术，研究了4种不同桨间距的6×6构型对转螺旋桨的桨间气动干扰对其气动性能的影响，桨间距选择了4种不同方案。研究结果表明：在4种不同的轴向桨间距中，当桨间距为0.25倍螺旋桨直径时，共轴对转螺旋桨的平均推进效率最高，并且气动干扰导致的效率脉动幅度较小；随着桨间距的增大，前、后排桨受到的气动干扰强度都会减小，相比于后桨，前桨因气动干扰造成的脉动对桨间距更加敏感。可见共轴对转螺旋桨的桨间距会对螺旋桨的气动干扰，及气动性能产生较为明显的影响，在设计共轴对转螺旋桨时选择合适的桨间距，有利于提高螺旋桨气动性能。     

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

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

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

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

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

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。