基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

基于经纬仪的光束平差法定向方法改进

针对基于传统经纬仪工业测量系统使用光束平差法定向精度有待提高的问题，提出一种新的定向改进方法。该方法同时结合了精确互瞄法与光束平差法，简化了测量过程，提出新的约束，消减了光束平差法中定向阶段的基准点测量累积误差，实现了在非通视条件中完成多台经纬仪高精度定向的目标。试验数据表明：该方法可在10m范围内，在误差为±3″的测点交会精度下，维持0.1mm的定向精度。     

预旋系统稳态和非稳态计算对比研究

针对预旋系统由于旋转引起流场和温度场参数周期性瞬态变化的问题，分别采用滑移网格的瞬态法和固定转子相位的稳态法进行数值求解，并在稳态计算中通过改变转子相位来近似模拟非稳态问题的时空变化特性。通过稳态空间平均结果与非稳态时均结果的对比，以期为预旋系统非稳态问题的低成本求解提供方法依据。结果表明：稳态计算结果与非稳态计算结果相比，周期性波动频率一致，接受孔进口处，稳态计算的压力波动幅度小39%左右，温度波动幅度小15%左右；多个相位的稳态空间平均结果与非稳态时均结果相比，压力高0.2%，温度低0.1%；采用多个计算周期数与单个周期的非稳态计算结果差异微小。当采用喷嘴出口气流中心正对接受孔迎风面前缘的转子相位时，稳态计算结果与非稳态计算时均结果最接近。     

核热推进地面试验技术研究

核热推进具有比冲高、推力大等特点,在载人深空探测和星际货运任务上具有广阔应用前景。核热推进技术的研发需要进行大量地面试验。首先回顾了美国与俄罗斯的核热推进地面试验技术的发展,对地面试验技术进行分类总结。然后基于一种小型核火箭方案,研究了燃料元件非核试验、燃料元件辐照考验试验和带核整机地面试验等关键地面试验技术,并提出了初步试验方案。最后对我国核热推进地面试验的发展提出了一些建议。     

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

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

考虑非理想噪声干扰的直升机飞行动力学模型辨识方法

利用MATLAB/SIMULINK平台搭建了直升机仿真实验系统，并以UH-60直升机为对象实现了考虑非理想噪声干扰的仿真飞行实验。建立了非理想噪声数学模型并基于增广最小二乘法设计 了一套同时考虑模型参数和噪声参数的综合辨识算法，在此基础上，利用仿真实验数据实现了UH-60直升机纵向飞行动力学模型的辨识与验证。最后，通过与普通最小二乘法的对比验证了本文方法的优越性。      

高雷诺数湍流非预混火焰及NO生成的大涡模拟

采用两种亚网格湍流燃烧模型,即化学建表方法结合假定概率密度模型和稀疏拉格朗日过滤密度函数方法,对高雷诺数湍流非预混火焰Flame D进行数值研究,定量比较不同亚网格模型的差异,并对火焰特征、污染生物生成特性进行分析。结果表明,两类亚网格燃烧模型预测的温度及大组分分布相近,稀疏拉格朗日过滤密度函数方法可以更好地模拟CO质量分数分布。不同的假定概率分布均可合理描述湍流与火焰的相互作用,之间的差别主要体现在NO分布,Dirac函数远高估了NO生成,而Top-hat函数则略低估了NO生成,Beta函数表现最优。Flame D的高温区及NO质量分数均主要分布在当量混合线及富燃侧附近。受高温伴流的影响,NO质量分数与温度一直保持高度正相关,峰值主要集中在标量耗散率很小的区域。不同截面上,反应物中的O_2和生成物中的H_2O均与NO高度相关。     

热障涂层表面裂纹水浸超声检测

为有效检测热障涂层(TBCs)表面裂纹,进行了热障涂层表面裂纹水浸超声检测技术的应用研究。以常用的大气等离子喷涂(APS)热障涂层作为试验材料,采用水浸超声检测回波声程差和C扫描图像对不同深度、宽度和长度的热障涂层表面裂纹进行了表征。结果表明:热障涂层基体材料对热障涂层表面裂纹的检测效果无明显影响;通过优化检测参数和合理选用仪器探头,利用水浸超声方法可实现对热障涂层表面深度为0.1 mm、长度为2 mm裂纹的有效检测;水浸超声检测技术可以获得裂纹的长度、深度参数,是解决热障涂层表面裂纹检测问题的有效方法。      

转静盘腔流动瞬态响应特性研究

以航空发动机转静盘腔内非稳态流动为研究对象,基于用户自定义函数编程设置进口压力按斜坡函数变化时的非稳态边界条件,利用数值模拟方法,在滞后效应、流场演化和物性参数方面,研究了不同无量纲跃升幅值下转静盘腔流动瞬态响应机理。结果表明:无量纲跃升幅值由1.05增至1.2,盘腔内部的当地响应时间延长了81.1%,相对总温的超调量增加了157.0%,其峰值时间增加了31.3%,且进口压力波动的传播方向总是由进口指向出口;瞬态过程中,转静盘腔的容积效应和旋转效应共同作用使出口流量明显滞后进口流量,且其滞留质量流量与相对总温的超调现象密切相关;冲击射流、旋转贴壁射流和静子壁面侧的涡系构成了盘腔流场演化的主体,无量纲跃升幅值越大,流场演化过程越剧烈,响应时间越长;声速和热扩散率与盘腔响应速度密切相关。