非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

三角翼大迎角风洞试验支架干扰数值模拟研究

现代战争要求战斗机能够在大迎角（AOA）状态下进行过失速飞行,对飞机大迎角绕流流场的研究主要的方法有风洞试验和数值模拟。在大迎角风洞试验中,常用的是尾支撑方法,支架的存在会对模型的试验结果产生一定的影响,本文通过数值模拟来对这个影响进行研究。以开源计算流体力学软件OpenFOAM 2.3为平台,采用PIMPLE算法求解Navier-Stokes（N-S）方程,PIMPLE算法是SIMPLE（Semi-Implicit Method for Pressure-linked Equations）算法和PISO（Pressure Implicit with Splitting of Operator）算法的结合体;采用基于有限体积的空间离散方法和空间二阶精度的线性插值方法,时间离散采用后向差分方法,湍流模型采用SA-DDES（Spalart-Allmaras-Delayed Detached Eddy Simulation）模型。为了验证方法的可靠性,首先对0°、10°、30°、50°、70°以及90°迎角下的有支架三角翼绕流流场进行计算,并将计算结果与试验结果进行对比,两者吻合较好。在此基础上,数值模拟了无支架的三角翼绕流流场,对比有/无支架情况下数值模拟结果,得到支架对三角翼绕流流场、背风面压强分布和气动力的影响。计算结果表明：大迎角情况下,有支架与无支架时相比,支架的存在会影响三角翼附近的流场（但是不会改变涡系等流动结构）、改变翼表面压强分布,从而导致三角翼的法向力系数和俯仰力矩系数发生明显变化。     

轴流压气机机匣变形对多排转子流场特性的影响

在实际轴流压气机加工、装配、使用过程中,机匣的圆度控制、叶片的高度控制、转子和机匣之间的同轴度控制等都普遍存在误差,意味着叶尖间隙的不均匀性必然地存在于每台压气机中。随着对轴流压气机流动研究的不断深入,叶尖周向非均匀间隙成为了需要考虑的问题。针对由机匣变形造成的非均匀叶尖间隙,研究了其对多排转子流场特性的影响。首先介绍了表征机匣变形程度的新参数,进而引进了非轴对称压气机模型的建模方法。针对3种不同的间隙周向布局,采用定常和非定常方法进行了数值模拟,对比了均匀与非均匀间隙下多排转子的气动性能,分析了机匣变形对气动损失分布和传播的影响,并研究了非均匀间隙下非定常压力和气动力的脉动特征。结果表明：非均匀间隙会降低转子性能,转子叶排叶尖区域的流场在非均匀间隙下呈现明显的周向非对称性,在上游叶排“尾迹”的干扰下,后排转子叶尖流场的周向非对称性要大于前排转子。非均匀间隙布局与各叶片气动力分布呈现明显的对应关系,小间隙区域叶片的气动力高,大间隙区域叶片气动力相对较低,从而增加了非均匀间隙下转子叶片的气动力幅值。     

连续旋转爆轰发动机冷流场的混合特性研究

连续旋转爆轰发动机（CRDE）中燃料和氧化剂的快速掺混是实现爆轰波成功起爆和稳定传播的重要前提,然而目前国际上关于这方面的研究还相对较少。本文采用大涡模拟（LES）方法,对非预混CRDE中燃料和氧化剂的混合过程及其主要机理开展深入研究。研究结果表明,非预混CRDE流场中存在欠膨胀特征、大尺度涡结构,以及回流区等复杂的流动现象,其中由于Kelvin-Helmholtz（K-H）不稳定性产生的大尺度湍流涡结构是促进氢/氧混合的主要机制。此外,本文还考察了氧气喷注位置对非预混CRDE的流场结构和混合特征的影响,发现氧气喷注位置会影响射流剪切层形态、涡尺度,以及回流区分布等,进而影响氢气和氧气射流的混合过程和混合程度。与其他进气位置相比,氧气在靠近内壁面喷注时更有利于氢/氧的快速掺混。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

考虑控制系统故障的重复使用运载器轨迹优化

传统的重复使用运载器的制导系统与控制系统的设计是分开进行的。当运载器出现故障时,需要进行轨迹和控制的重构,故障下气动特性和配平能力的变化使得重构后的制导系统与控制系统难以良好匹配。考虑到控制系统对轨迹设计和优化的影响,通过将三自由度轨迹优化与多舵面控制分配结合,提出一种考虑控制舵面气动力和故障的重复使用运载器轨迹优化方法。仿真结果表明,多舵面控制分配算法可以将运载器故障下的配平能力转化为轨迹优化的约束条件,从而避免了传统三自由度轨迹优化方法结果无法被控制系统跟踪的可能性,并为制导与控制一体化设计提供参考依据。     

基于三维非结构混合网格的离散伴随优化方法

伴随优化方法在优化过程中不受设计变量个数的限制,有巨大的工程应用前景。基于非结构混合网格技术,采用雷诺平均Navier-Stokes方程、离散伴随优化方法以及自由形面变形（FFD）技术,建立了飞行器气动外形优化设计系统。离散伴随方程采用公式推导的方法直接获取,并运用LU-SGS迭代求解。通过与差分结果进行比较的方式,对通量的雅可比矩阵和目标函数敏感导数开展了确认,并对文中黏性通量项、限制器以及湍流黏性系数的近似处理方法做了分析对比,验证了近似处理方法的合理性和可行性。采用建立的优化系统,完成了ONERA M6机翼在跨声速条件下的减阻优化,并比较了有无容积约束下优化结果的差别。模拟验证结果表明,建立的飞行器气动外形优化设计系统具有较高的可靠性和有效性,可以用于三维飞行器外形的减阻优化。     

在轨服务中非合作目标空间圆检测方法研究

为自主完成空间服务任务,需要满足各种功能需求、突破多种关键技术。针对在轨服务中具有圆特征的非合作目标空间圆（星箭对接环、发动机喷管等）,首先分析了非合作目标空间圆检测技术在在轨服务系统中的主要应用;然后提出了一种非合作目标空间圆的检测方法,通过Canny算子检测边缘,并用Freeman链码法对边缘进行提取分类,再利用RED算法进行非合作目标空间圆检测;最后给出仿真结果,在保证检测精度的前提下,较传统的RED算法明显降低了算法的耗时量。     

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

基于GPM的高超声速飞行器再入轨迹优化

为保证高超声速飞行器具有最大再入距离,应用高斯伪谱方法求取了存在路径约束及终端约束的飞行器最优再入轨迹。首先选取迎角及倾斜角作为最优控制量,并分析再入轨迹的最优控制模型;然后应用高斯伪谱方法将最优控制问题离散成非线性规划问题,并采用SNOPT软件包对非线性问题进行求解,得到最优轨迹。通过仿真证明了该方法求得的最优轨迹对初始值不敏感,且具有良好的实时性。