用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

凹腔/支板结构亚燃冲压燃烧室性能

为了避免基于凹腔火焰稳定器的亚燃冲压燃烧室壁面喷注时燃料与主流空气掺混非均匀性问题和提高燃烧室的性能,提出在亚燃冲压燃烧室中使用支板喷注代替壁面喷注的方案,数值模拟了凹腔/支板结构亚燃冲压燃烧室中燃料分布及流场结构,并分析了支板结构对燃料空气混合及燃烧室性能的影响。研究表明:支板虽然使燃烧室出口的总压恢复系数相对于壁面喷注方式下的降低了63%,但能使燃料均匀分布于整个流道内,增强了燃料与空气掺混,使燃烧室出口的混合效率和燃烧效率分别提高了21.4%和20.5%。燃烧效率的提高弥补了采用支板导致的燃烧室内气流的额外总压损失所带来的机械能损失,使得支板喷注时燃烧室出口的比冲提高了39.6%。因此,在亚燃冲压燃烧室中设置凹腔/支板结构,有利于提高燃烧室整体性能。      

2198和5A90铝锂合金脉冲阳极氧化膜制备及耐蚀性

在18%(质量分数,下同)H_2SO_4+5%C_2H_2O_4水溶液中,采用脉冲(PC)电流对2198和5A90两种铝锂合金进行阳极氧化处理。用扫描电子显微镜(SEM)观察铝锂合金阳极氧化膜表面和截面形貌;用能谱仪(EDS)对其成分进行面扫描和线扫描;用动电位极曲线检测氧化膜在3.5%NaCl水溶液中的耐蚀性。结果表明:2198和5A90铝锂合金阳极氧化膜主要由Al的氧化物组成;2198合金氧化膜表面存在细小颗粒,厚度约为150μm;5A90合金氧化膜表面存在微孔,为后续封孔处理提供结构条件,厚度约为180μm;用脉冲方法在两种铝锂合金表面生成的较厚阳极氧化膜具有较高的耐蚀性。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

基于混合粒子群算法的上升段交会弹道快速优化设计

基于梯度搜索的高效性和粒子群搜索的随机性,提出了一种混合粒子群算法,并应用该算法研究了运载火箭上升段交会弹道快速优化设计问题.以运载火箭与目标飞行器在交会时刻的距离最小为目标函数,设计了运载火箭飞行程序,建立了运载火箭上升段交会弹道优化模型,同时分别采用混合粒子群算法、遗传算法和粒子群算法进行求解.仿真结果表明:基于本文算法对运载火箭上升段交会弹道进行优化设计,平均交会位置误差为4.137m,较遗传算法减少了17.940m,平均优化耗时488.922s,较粒子群算法缩短了2342.125s.混合粒子群算法搜索速度较快,收敛精度较高,可用于运载火箭上升段交会弹道的快速优化设计.      

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

超临界压力RP-3在竖直细圆管内混合对流研究

研究了超临界压力下碳氢燃料航空煤油RP-3在竖直细圆管内混合对流,分析了浮升力及热物性对碳氢燃料在垂直管中对流换热的影响。实验中控制热流密度从200~500 kW/m2变化,进口压力变化范围为3~5 MPa,进口雷诺数从5 000~10 500范围内变化。研究表明:在向上流动情况中进口段存在较为明显的入口效应,换热出现恶化现象,而在向下流动中未出现;对于向上和向下流动,由于热物性的综合影响,换热系数沿流动方向增大;在较低进口雷诺数(Re=5 700)时,对于向下流动,随着浮升力影响的增大,浮升力改变了流体径向速度分布,出现了换热强化;在较高进口雷诺数(Re=10 500)时,浮升力对换热的影响依然显著;判别式Bo*数小于5.6×10-7未能预测浮升力对碳氢燃料换热影响。      

改进的多处理器混合关键性系统可调度性分析

针对混合关键性系统的多重认证需求,研究多核处理器平台中全局调度算法fixed-priority and Earliest Deadline First by Virtual Deadline（fpEDF-VD）的可调度性分析问题。fpEDF-VD结合处理器利用率和虚拟截止期两个方面来计算任务优先级,系统可调度性取决于是否存在可行的虚拟截止期调整参数。考虑到现有可调度分析方法仅测试有限数量的调整参数候选值,不能有效地判定系统可调度性,故提出了一种改进的判定方法。该方法基于传统（非混合关键）任务调度算法fpEDF的可调度利用率约束条件,利用函数图像分析研究不同关键性级别的系统可调度性需求,并在此基础上给出有效虚拟截止期调整参数的确切范围。通过实例分析及与现有判定方法的比较,验证了该方法的正确性和高效性。与理论分析一致,基于随机生成任务集的仿真实验结果表明改进后的方法具有更优越的可调度性能,能显著地提高任务集的可调度接受率。