基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

柴油发动机热管理系统的高度特性

基于空中和高原环境的特点,研究了柴油发动机散热器传热性能的高度特性;进一步建立了发动机热管理系统的流动和传热模型,研究了热管理系统的高度特性;给出了基于最小功耗原则的热管理系统运行参数匹配方法.在固定系统热负荷条件下的计算结果表明:随着海拔高度的增加,在空中环境下,散热器的散热性能先上升后下降,热管理系统中冷却液和润滑油温度先下降后上升,热管理系统所需的最小功耗先减小后增加;在高原环境下,散热器的散热性能随高度增加大幅下降,热管理系统中冷却液和润滑油温度大幅上升,发动机出现过热问题,为了给发动机提供足够的散热能力,热管理系统所需要的最小功耗随着海拔高度的增加而大幅增加.      

离心叶轮全三维反问题方法设计误差稳定性研究

用流线曲率全三维反问题方法和引入滑移因子模型的准三维流线曲率反问题方法设计了叶片数不同的和增压比不同的小型高速离心叶轮.用商用计算流体力学(CFD)求解器估算了这些叶轮的性能.对比了这两种反问题方法的设计误差稳定性,考察了分流叶片对设计误差的影响.结果表明,流线曲率全三维反问题方法能够根据叶片负荷对叶片设计做出恰当调整.较准三维反问题方法,流线曲率全三维反问题方法所设计叶轮的增压比普遍偏高.随叶片负荷增加,较准三维反问题方法,流线曲率全三维反问题方法的设计误差变化较小,即其设计误差稳定性明显较好.分流叶片对准三维反问题方法设计结果的影响明显较小.      

细观周期性结构复合材料热固耦合双尺度渐进均匀化分析方法及有限单元法实现

针对周期性结构复合材料,根据稳态热固耦合平衡方程,将弹性矩阵和导热系数矩阵进行合并,将位移和温度作为空间场量同时进行小参数渐进展开,利用摄动理论和均匀化理论,推导并建立了周期性结构复合材料稳态热固耦合双尺度渐进均匀化分析方法。将复合材料力学分析从单一物理场推广到热固耦合场,为实现复合材料构件的热固耦合多尺度力学分析建立了基础。根据建立的偏微分方程式,利用变分原理推导得到其有限单元形式,用VC++6.0开发了热固耦合双尺度渐进均匀化分析程序。为验证该方法预测复合材料宏观等效性能的准确性,对纤维体积含量从20%到80%的单向纤维增强复合材料宏观等效参数进行了计算。当纤维体积含量低于65%时,得到的计算值与实验值吻合良好。更加真实地还原复合材料细观结构模型,计算精度将进一步提高。     

气相化学反应与声学振荡的耦合机理数值研究

通过数值仿真方法针对液体火箭发动机内的气相化学动力学与振荡声场的热声耦合机理进行了研究。采用任意拉格朗日算法解耦流动与化学源项间的刚性。采用的多步总包反应机理考虑了底层的准稳态组分脉动。通过入口流量边界的流量脉动向燃烧室中引入纵向声波,并建立了冷流声学相似场模型以分析热声耦合效应的强度。研究发现:在线性小振幅声场中,气相化学动力学控制的释热系统与声学振荡无明显耦合激励;在非线性有限幅值声场中,燃烧室压力与释热波动出现“突跃”并表现为陡峭前沿波,气相化学动力学控制的释热系统与声学振荡发生耦合激励,反应流较其冷流声学相似场的压力振荡振幅增强约300%。最后分析了耦合激励发生的可能原因,提出了气相化学动力学体系的“释热分岔”假说。     

基于OpenFOAM的真实流体模型库及均相求解器开发

为了正确认识高压环境下的流体非理想性特征,及其对于低温射流喷雾过程的影响作用,将多种不同真实流体模型写入到开源CFD程序OpenFOAM仿真平台中,编译形成具有广泛通用性的真实流体热物理模型库;并将高压下真实流体特征的影响引入到压力修正算法中;开发出适用于描述超临界环境下低温流体流动过程的均相求解器。对高压下流体热物理性质的非理想性,以及超临界环境下的低温液氮射流喷注过程进行研究。结果表明:高压下的真实流体效应主要体现在低温区域,本文所建立的模型库可以在任意热力学状态范围内正确预测流体热物性;基于此模型库所开发的均相求解器可正确预测超临界环境下低温射流的详细流动特征和结构,而理想气体模型则会大幅低估射流核心区密度分布,中心线上误差高达68.9%;综合考虑计算精度和效率,PR EoS及相关模型更适合于这类超临界流动问题的求解;低温射流跨临界喷注中局部区域所发生的虚假沸腾现象会导致射流温升速率降低、体积膨胀速率加快;射流核心区压缩因子明显1,在数值仿真时必须考虑高压环境下的真实流体效应。     

壁面催化效应对高超声速气动热影响研究

高超声速飞行器壁面催化效应会导致激波层中原子在壁面处复合释热,加剧周围气动热环境。针对高超声速流动壁面催化特性,选择不同飞行马赫数及高度条件,采用完全催化和非催化两种条件对球锥模型壁面热流率进行数值模拟计算,研究壁面催化效应对气动热的影响规律。结果表明,固定飞行高度时,壁面催化效应对气动热的影响随马赫数增加而加强,Ma=25条件下驻点处完全催化与非催化热流比值高达1.92;固定飞行马赫数时,在50km高空以上壁面催化效应对气动热的影响随高度增加而减弱;壁面催化效应不仅会影响壁面附近的流场特性及组分分布状态,而且对整个激波层都有一定的影响作用。     

可重复使用热防护材料应用与研究进展

可重复使用热防护系统是为高速重复使用飞行器而发展的关键性技术,涵盖了地球大气环境及非地球大气环境下的弹道式再入、高马赫数巡航等应用场景。根据现有高马赫数飞行器热防护现状,对高马赫数飞行器的主要热防护系统类型、特点和使用场景进行了简要介绍。在此基础上,结合国外里程碑式可重复使用飞行器（X-15、SR-71、航天飞机、X-33、X-37B、Spaceliner等）,梳理了可重复使用热防护材料的应用与研究进展,论述了代表性可重复使用热防护材料的发展、性能、研制进度、特点及应用前景。对国外在可重复使用热防护材料研制中的设计及发展思路,以及所存在的主要问题进行了总结归纳,为可重复使用热防护材料未来的发展提供了思路。