基于离散数据点的变壁厚叶身参数化设计

在原始截面离散数据点的基础上,给出了面向工程实际的中弧线几何造型方法和基于离散数据点的内形曲线参数化方法;深入分析了UG(Unigraphic)系统参数化造型的实质;给出了可行的参数化方案;结合面向对象的思想和创成式的建模方式,实现了变壁厚叶身参数化设计.该模块基于UG/OpenAPI接口,采用即时记录数据点和标记回溯点的方式,在交互中可以实时动态修改,达到"所见即所得"的效果,并在实际应用中得到有效验证.      

广义多圈Lambert算法求解多脉冲最优交会问题

基于一种高效高精度的Battin多圈Lambert算法提出一种考虑轨道摄动的广义多圈Lambert算法.与现有算法相比,本算法虽然原理复杂但计算流程非常简单,效率极高,分别通过几次内外循环就可满足精度要求.广义多圈Lambert算法结合一种可行解迭代交会模型构成了一个通用的多圈多脉冲交会规划框架,应用两步法求解此多变量的复杂工程优化问题,首先利用高效率的进化全局优化算法以及解析轨道模型作全局搜索,然后利用序列二次规划算法以及简化高精度轨道计算模型作局部搜索,此方法可以保证高效高精度的求解多圈多脉冲交会问题.算例表明此方法特别适用于满足实际工程约束的交会规划问题.     

临近空间大气中波红外辐射特性的数值分析研究

介绍了临近空间大气环境特点及其非局域热平衡状态的辐射传输计算方法。利用战略高高度辐亮度代码SHARC（Strategic high-altitude radiance code）,计算了3—5 μm波段不同观测条件、光学现象及大气环境下的临近空间大气背景辐射。分析结果表明：临近空间中波红外辐射随观测天顶角的增大而增强,随太阳天顶角的增大而减弱;在38 km及75 km附近,临边背景辐射存在极大值;纬度越高,临边背景辐射的季节变化特征越明显。OH夜气辉和极光对背景辐射有重要增强作用。     

小半径弯管弯胀复合成形研究

为解决小半径弯管内侧起皱、外侧减薄的成形缺陷,提出了一种弯胀复合成形的方法。通过几何关系解析了初始弯曲半径和初始管坯直径与最终管径和最终弯曲半径之间的关系,确定了初始弯曲半径的取值范围。采用有限元与实验相结合的方式研究了不同初始弯曲半径对零件壁厚分布的影响,分析了胀形过程中零件弯曲处及直壁处的变形规律。最终发现:初始弯曲半径越大,零件的最终减薄越小,胀形过程中弯曲处内外侧发生均匀减薄,直壁处内侧减薄明显大于外侧减薄。     

针对旋翼流场模拟的自适应网格算法

针对基于格心格式求解器的旋翼流场模拟,提出了相应的自适应笛卡尔网格的数据存储结构及自适应算法。给出了相应的单元处理策略,简化了对自适应笛卡尔网格的处理;对于频繁的自适应加密过程中产生的大量重复点,采用了高效的交替数字树算法(Alternating digital tree,ADT)予以删除;对于自适应疏化过程中产生的大量无用点,提出了标记-删除-移动(Mark,delete,move,MDM)算法予以快速地删除,减少了不必要的计算资源消耗。对CaradonnaTung旋翼在不同悬停状态下进行了模拟验证,对比了压力分布系数与桨尖涡位置。之后对HELISHAPE 7A旋翼在前飞情况下进行了模拟验证,计算值与实验值吻合。此外,求解器对桨尖涡的捕捉效果得到了明显的提高,表明本文方法具有良好的有效性与鲁棒性。     

端部状态对斜置圆柱气动力分布的影响

在有限长度圆柱气动力测试风洞试验中,端部状态的变化对气动力有很大影响。在雷诺数为6.43×104时,通过刚性模型风洞测压试验对3种端部状态下的斜置圆柱气动力特性进行研究,基于平均风压系数沿斜置圆柱的轴向和周向分布规律,分析端部状态对斜置圆柱气动力测试结果的影响。结果表明:斜置圆柱表面风压沿轴向可分为近上游端部区、中间区和近下游端部区;由于马蹄涡的存在,在近上游端部区的背风面,风压沿轴向有一个突变,端部状态对马蹄涡强度有很大影响;在中间区,背风面风压沿轴向出现交替变化,可能是由于在近上游端部区形成的旋涡受到轴向流作用向中间区流动导致的,端部状态对交替风压变化的幅度和阻力系数有很大影响;当下游端部封闭时,近下游端部区的迎风面和背风面风压都会出现突变,且迎风面的突变幅度更大,下游端板对近下游端部区迎风面的风压分布影响更大。     

高焓电弧风洞试验热化学非平衡流场数值模拟

针对高焓电弧风洞内部流动的热化学非平衡效应及气体组分和振动能量冻结效应导致的试验数据外推困难问题,基于高焓风洞喷管/试验段/试验模型一体化数值模拟的思路,通过数值求解三维热化学非平衡Navier-Stokes方程,开展了FD-15高焓电弧风洞典型运行状态下流场的数值模拟,与典型试验状态的气动热数据进行了对比验证,研究了试验数据外推飞行条件的方法及有效性问题,分析了提高驻室总压对试验数据外推的影响。研究表明:（1）风洞试验段来流离解度高,热化学非平衡效应及其冻结现象严重;（2）热流校核试验测量数据位于一体化数值模拟的完全催化热流和非催化热流之间,分布合理,验证了计算方法和程序的正确性;（3）试验模型安放位置对模型表面压力和热流存在影响,模型与喷管出口的距离越大,模型表面压力和热流越低;（4）当驻室总压较低时,通过双尺度模拟准则（模拟飞行条件总焓和双尺度参数ρ_∞L）外推热流失效,使用部分模拟准则（模拟飞行条件总焓和驻点压力）外推热流也会出现较大差异,在非催化条件下这一现象更加明显;（5）当驻室总压较高时,使用双尺度模拟准则或部分模拟准则外推飞行条件,产生的热流差异明显减小。     

氮气辐射强度的激波管测量与验证

超高速飞行器的热防护设计必须考虑激波层内高温气体发射与吸收的辐射能量,需要有效的辐射加热评估手段。相应飞行条件下的光谱辐射强度地面实验测量是验证数值模型和方法、理解高焓流动的重要手段。基于燃烧驱动激波管,发展辐射强度标定技术,针对富氮气环境,开展高温气体光谱辐射强度的高分辨定量化测试,掌握辐射特征,为数值验证提供基础数据。实验获得了激波速度5.70和6.20km/s条件下的气体光谱辐射强度精细结构,数据表明激波波后的非平衡过程对辐射强度存在很大影响。通过求解耦合化学反应动力学模型的Navier-Stokes方程和辐射特性模型,得到对应实验条件下的流场参数和辐射强度,计算结果和实验数据符合很好,验证了数值模拟方法。     

一种非介入式高超声速边界层不稳定波的测量方法

地面风洞实验是开展高超声速边界层转捩研究的主要手段之一,但是目前可用于高超声速边界层三维空间测量的实验技术仍极为缺乏,且已有测量技术的动态响应频率普遍较低。基于光的折射和干涉原理,搭建了一套非介入式聚焦激光差分干涉仪测量系统（Focused Laser Differential Interferometry,FLDI）,可有效获取三维流场空间点的密度变化。在马赫数为8的常规高超声速风洞中,使用FLDI开展了来流雷诺数107/m、7°半锥角尖锥标模边界层的不稳定波测量实验。结果显示FLDI成功捕获到频率在327 kHz的第二模态不稳定波及其谐波（645 kHz）。通过与PCB测试结果进行对比,FLDI的高信噪比、高解析频率（本文实验有效解析频率1.5 MHz）、高空间分辨率（沿流向小于1 mm）等优点得以体现。鉴于FLDI的高时空分辨率等优良特性,其可用于高超声速边界层不稳定波行为以及感受性等问题的研究,为深入认识高超声速边界层转捩机制以及感受性问题提供了有效手段。     

基于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格式具有更低的数值耗散。