相变材料热控系统内部接触热阻的辨识方法研究

在相变材料热控系统的吸热融化过程中,由于材料与容器壁面之间的空穴效应形成的接触热阻,将对熔化传热形成阻尼效应.采用反演热传导反问题的方法,将参数辨识中的灵敏度法和共轭梯度法应用到实验,建立了一种可用于材料熔化过程界面稳态和瞬态接触热阻的辨识方法.模拟辨识与初步应用表明:辨识方法计算精度高、稳定性好,为精细化分析相变材料热控系统中的熔化吸热过程奠定基础.     

融化过程自然对流影响作用的实验研究

采用实验的方法对相变材料融化过程中自然对流的影响进行了研究,获得了自然对流作用下的融化速率与平均Nusselt数的准则关系式.结果表明,与纯导热模型相比,自然对流作用加速了融化的进行,建立导热-对流耦合传热模型对于准确预测融化过程具有一定的必要性.     

进入火星大气的高温真实气体效应与气动加热研究

针对火星和地球大气分子热力学和化学行为的差异性,采用理论分析和数值模拟两种手段,研究探测器进入过程高超声速流动的分子振动激发、离解反应及热力学和化学非平衡等真实气体效应,获得不同气体模型条件下的高超声速气动加热规律,探究引起地火差异的根本原因。分析认为,探测器进入火星大气层的稀薄气体效应明显;激波层内发生CO 2气体为主的大规模离解,在极高温环境下O 2和CO也将离解;沿进入轨道的高超声速流动基本处于化学非平衡但热力学平衡状态;激波层内能量储存和分配模式因分子振动激发和化学反应而改变,分子振动激发会增强气动加热量,但均介于化学反应模型的完全非催化和完全催化壁结果之间;相同来流条件下CO 2介质高超声速气动加热强于空气介质,但真实的火星进入热载荷因大气稀薄而弱于地球再入环境。相关研究为我国未来火星探测器热防护系统设计提供技术支持。     

电加热过程的冰脊形成实验研究

利用结冰风洞设备和电加热装置,采用实验的方法研究了不同来流速度、环境温度和加热功率对冰脊形成的影响。研究表明：随着来流速度、加热功率的增加,冰防护区长度增加,冰脊往下游推后,而随着环境温度的降低,冰防护区长度减小,冰脊往上游靠近。冰脊的生长规律是从冰防护区外下游某位置开始冻结,逐步往冰防护区发展,从而为电热除冰系统加热模式的选取和传热特性的优化提供了实验依据。     

国家数值风洞(NNW)工程中的CFD基础科学问题研究进展

计算流体力学（CFD）基础理论方面的持续创新对数值模拟软件的功能拓展及推广应用具有十分重要的意义,国家数值风洞（NNW）工程在CFD相关的若干关键基础科学问题上开展了研究。简要综述了国内各参研团队近3年的研究进展,NNW工程在转捩与湍流模型及计算方法、多相多介质计算模型与方法、多物理场耦合计算模型与方法、高精度数值计算方法等方面形成了一系列阶段性的研究成果,建立了多种原创算法与模型。这些算法与模型中具有较高成熟度的将被集成至NNW工程的相关软件,并面向全国发布。     

高超声速飞行器界面多相催化数值研究进展

随着未来临近空间高超声速飞行器高速度、长航时新需求的提出,飞行器高温流动与热防护系统相互作用凸显,引发极端力学、热学条件下气固界面多相催化等高温界面效应。回顾了高超声速飞行器中界面多相催化理论建模和数值研究历程,重点综述了界面多相催化的给定速率系数模型、含微细观特征的唯象模型、基于微观理论模拟的跨尺度模型的研究进展。总结了作者团队在飞行器界面多相催化效应建模、机理和应用相关方面的研究结果。结合未来飞行器减重、增程、保形的设计需求,进一步提出了国内后续研究的重点方向,以期支撑热防护系统轻量化、低冗余设计。     

高超声速飞行器防热瓦结构的变厚度轻量化设计方法

针对非均匀气动加热条件下三维复杂防热瓦结构的轻量化需求,基于变厚度设计理念,发展了一种基于网格变形技术（ASD）和热固耦合分析相结合的防热瓦结构优化方法。优化结果表明,基于网格变形技术能够快速有效地解决优化过程中的网格自动更新问题,避免了网格重新划分的耗费及复杂结构网格重构的困难,并得到了光滑柔顺的厚度形状曲线;相比等厚设计,变厚度设计可以有效考虑载荷的非均匀效应,并极大地减轻结构重量;优化后可以更充分发挥防热瓦结构各层材料的承载能力。     

相变储能单元融化过程的传热强化

采用焓-多孔介质模型.对垂直热壁边界条件下相变储能单元的融化过程进行了数值模拟,研究了板式储能单元翻转对融化过程的传热强化作用及Rayleigh数、Stefan数和单元高宽比对翻转传热强化作用的影响.研究表明,在融化过程的适当阶段翻转储能单元,融化速率可提高约30%.Rayleigh数越大或Stefan数越小,传热强化作用越明显;单元高宽比的改变对翻转传热强化作用的影响相对较小.     

结冰风洞水滴直径标定方法研究

过冷水滴粒径大小是重要的结冰云雾参数,获知结冰风洞中的水滴直径,是得到定量结冰风洞实验结果的基础。对于结冰风洞内水滴直径单一或者分布比较集中的情况,提出了一种采用数值计算和结冰风洞实验相结合的手段标定水滴直径的方法。该方法首先采用拉格朗日法数值计算水滴运动轨迹,得到撞击极限随水滴直径变化的关系曲线,在此基础上,进行结冰风洞实验,测量实验得到的水滴撞击极限,通过在撞击极限与水滴直径关系曲线上进行插值,进而得到实验水滴直径大小。采用该方法对0.3m×0.2m结冰风洞内的水滴直径进行了标定,分别计算和测量了25m/s和35m/s两种速度条件下的水滴撞击极限,得到的水滴直径值相差不超过1μm,初步说明该方法的合理性。同时,对于结冰风洞内水滴粒径多尺寸分布的情况,还提出了相应的标定其容积平均直径MVD的方法,该方法在计算水滴收集率的基础上,通过测量驻点处的结冰厚度,实现对MVD的测量。采用本文提出的两种方法进行结冰风洞水滴粒径标定,只需要一般的长度测量工具即可进行,操作方便,成本低廉,克服了常规的水滴直径测量或标定需要专门设备的不足。     

考虑非平衡效应的过冷水滴凝固特性

非平衡凝固是过冷条件下水滴凝固过程的重要现象。本文针对飞机结冰过程过冷水滴的非平衡凝固效应,发展了改进的凝固特性预测模型及数值计算方法,并自行搭建了实验系统,开展了所建过冷水滴凝固模型与数值预测方法的实验验证。研究表明,所发展的改进模型可有效表征水滴过冷阶段的非平衡凝固效应,因而对冷水滴凝固速率的预测有较好的改进;当过冷度为0℃时,过冷模型退化为传统模型。基于所建方法,开展了过冷度及冷却条件对水滴凝固特性的影响分析,获得了不同条件下水滴凝固过程的温度分布及相界面变化特征。研究表明,过冷度越大或水滴尺度越小,凝固速率相对越高;在考虑非平衡凝固效应的条件下,过冷水滴凝固速率要高于不考虑非平衡凝固效应的工况。相关研究可为结冰热力学模型的改进,以及结冰特性的精细化预测提供参考。