波形隔板穿流孔孔径对通道流动和换热的影响

采用数值模拟的方法,对波形隔板结构复合通道的换热和流阻特性进行研究.设计了6种不同穿流孔孔径的隔板,研究孔径对通道中流动和换热的影响.将大小相间孔隔板结构的数值结果与已有实验结果进行对比,趋势符合较好.数值研究结果表明:在均孔隔板结构中,当孔径小时,通道换热好,同时流阻大;当孔径增大时,通道换热变差,同时流阻降低.当隔板上穿流孔总穿流面积一定时,大小相间孔的隔板结构与均孔隔板结构相比,换热略好些.     

离子发动机羽流二维轴对称数值模型与验证

为了解离子发动机羽流特别是交换电荷离子(CEX,Charge-Exchange)的分布和流动特性,建立了离子发动机羽流的物理模型,采用粒子网格(PCI,Particle in Cell)方法对2种型号的离子发动机羽流场进行数值模拟计算,与其地面实验数据进行对比分析.结果表明:在CEX离子密度大小及分布、电势的大小及最大电势梯度的位置、CEX离子流动角方面,模拟结果同实验结果符合得相当好.在电势结构方面,由于舱壁电势的影响,模拟结果同实验结果相比有一些差别.羽流模型和计算结果为相关羽流实验和数值模拟研究提供参考.      

层板内冷通道辐射换热影响

分析了层板换热过程,使用了"去除复杂表面"和"无限大平板" 等假设,给出了该假设下层板内腔辐射换热计算公式.又以一特定几何结构的层板模型为例,使用了3维流体力学计算程序求解了流-固耦合情况下的层板内部换热过程,研究了其在4种航空发动机典型工作状态下,内冷通道中的对流和辐射换热情况,得到了辐射与对流换热强度之比θ随冷气入口Re数和燃气加热功率的变化曲线,并对该曲线进行了指数拟合.建立了一套快速简捷计算层板内腔辐射状况的方法,为层板内腔换热研究是否应该忽略辐射,以及对层板冷却效果的修正,提供了一定的判断和计算依据.      

旋转光滑U形通道内流动和换热的数值模拟

数值模拟了旋转状态下涡轮叶片U形内冷通道湍流流场和温度场的分布,分析了流阻和换热的变化规律.结果表明,旋转状态下哥氏力、离心力和浮升力的共同作用使得流场发生了复杂的变化.旋转强化了换热,减小了流阻.但旋转使得换热在各个面换热能力分布不均,增加了温度梯度.      

开槽交错肋通道换热和流阻特性

对一种特殊结构的涡轮叶片内冷通道——开槽交错肋通道进行了换热和流阻的测量实验.全部的实验都在静止的状态下完成,其中换热实验采用了水蒸气凝结换热方案.实验的 Re数范围从5 000~45 000,得到了大量关于开槽交错肋涡轮叶片内冷通道的实验数据,并将结果与不开槽交错肋通道作出对比,得出一定槽宽的开槽交错肋通道相比于相同尺寸的不开槽交错肋通道,换热效果有了明显的提高,而流阻系数相应有所减小,所以整体换热效果优于不开槽通道.此外,对于不同槽宽通道的各项特性也进行了实验对比.得出 4 mm开槽交错肋通道的换热效果最好,接下来依次是2 mm、6 mm和8mm的通道.      

基于润滑理论的二维积冰数值模拟

将基于润滑理论的水膜-冰层模型推广到在任意二维截面上建立的正交曲线坐标系,得到描述二维型线表面水膜流动和积冰的偏微分控制方程.采用一种隐式-显式有限差分格式建立控制方程组的代数形式,给出求解非稳态控制方程组的数值方法.为验证模型和数值解法的有效性,在典型的航空积冰和传输线积冰环境下对翼型和传输线表面的积冰算例进行数值模拟.将水膜-冰层模型的冰形计算结果与传统Messinger模型的模拟结果以及冰风洞试验结果进行对比,低温明冰条件下采用当前方法计算的翼型表面冰形接近于Messinger模型的模拟冰形曲线;相对高温条件下的计算结果比传统Messinger模型更为精确.低速积冰环境下Messinger模型难以模拟的传输线表面积冰同样可以采用水膜-冰层模型进行有效 预测.     

轴对称羽流流场数值模拟和实验验证

针对在轨航天器羽流效应对航天器整体寿命的影响,在"DFH"系列卫星姿控发动机喷管羽流实验研究的基础上,建立喷管的轴对称模型,选取矩形计算区域,设定一定的边界条件和计算条件,采用轴对称单组分直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法对相应实验状态的羽流流场进行数值模拟研究.通过对比数值模拟结果和实验测量结果,表明DSMC数值模拟羽流流场的计算结果符合实际流动规律,实验验证了所编制的轴对称DSMC数值模拟程序的正确性和有效性.      

旋转状态下涡轮叶片前缘的流动与换热

用数值模拟的方法对旋转状态下涡轮叶片前缘冷却结构进行了数值研究,该结构由进气腔、叶片尾缘块和前缘块构成,对此结构不同的旋转速度情况进行了计算,根据计算结果分析了旋转对涡轮叶片前缘流动与换热的影响.计算结果表明,旋转状态下带气膜出流的冲击流动中,前尾缘冲击面的换热随着转速的增加而减小,且尾缘冲击面的换热比前缘冲击面的换热要好;同时前尾缘冲击面换热的差别随着转速的增加将越来越小.     

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

复合材料层板低速倾斜冲击损伤分析

纤维增强复合材料层板对低速冲击事件敏感，冲击产生的损伤会导致材料结构承载性能及使用寿命大幅下降。基于此，提出了一种基于连续介质损伤力学的有限元模型，研究了复合材料层板低速倾斜冲击力学行为。采用Hashin准则结合渐进退化模型预测层内损伤起始和演化；采用界面单元结合双线性Traction-Separation本构关系模拟层间分层；编写用户材料VUMAT子程序，实现基于ABAQUS/Explicit软件平台的数值求解。数值计算结果与现有正冲击下实验数据吻合较好，验证了模型的有效性。探讨了冲击角度、冲击能量对复合材料层板倾斜冲击力学性能的影响，分析倾斜冲击下层板损伤模式及失效机理，为复合材料结构倾斜冲击问题数值分析提供参考。      

直升机旋翼防/除冰电加热控制律仿真

电热防除/冰系统的控制涉及到电加热与外流场的传热耦合,计算较为复杂,可利用的数据资料较为稀少。为探索电热防/除冰系统工作时与外流场的耦合传热规律,建立了二维电热除冰的数学模型,该模型在Messinger模型和改进的焓法模型基础上耦合了外表面与环境的复杂换热以及融冰和重新结冰过程的相变换热;采用控制容积法对控制微分方程进行离散后,使用TDMA（Tri-Diagonal Matrix Algorithm）和ADI（Alternating Direction Implicit）方法对离散得到的线性方程组进行求解,进而得到了除冰表面温度分布,同时揭示了电热防/除冰系统的耦合传热规律;分析了不同结冰条件下,加热时间控制律和加热热流密度对除冰表面温度的影响。计算发现合理设计加热热流密度大小及分布和加热时间控制律,可实现电热除冰系统能源的高效利用,进而确保飞行安全。     

超临界压力RP-3壁面结焦对流阻的影响

实验研究了长时间加热条件下航空煤油RP-3在微细不锈钢管内流动过程中结焦对流动及换热的影响规律。实验中系统压力保持为5 MPa,燃油质量流量为3 g/s。在燃油溶解氧达到饱和的条件下,实验段进出口油温分别为130℃和450℃。实验从开始到终止持续36 h。实验结果表明,随着时间的增长,管内的结焦量不断增加。由于壁面结焦现象对管内流动和换热产生严重的影响,管内各段换热在实验前期迅速恶化逐渐趋于稳定,管内流阻随着实验时间的增加持续增长。管内流动阻力随着时间的增长呈现出"快速增长→平稳增长→急速增长"的过程。另外,基于实验结果,提出了一种影响系数作为判断结焦对换热器单管影响的工程模型。      

临近空间环境下封闭方腔内耦合换热特性

以临近空间浮空器载荷舱为应用背景,对复杂热边界条件下含热源的三维封闭方腔内自然对流、表面辐射和导热的耦合问题进行了数值模拟。综合考虑对流换热、长波辐射、太阳辐射等因素的影响,建立了临近空间热环境模型。通过Fluent软件用户自定义函数（UDF）引入外部非定常的辐射-对流耦合热边界条件,对腔内换热特性的昼夜变化进行研究,并分析了腔壁厚度、发射率和导热系数对其的影响。数值结果表明,腔内平均温度昼夜变化很小,约为12.9 K,但温度场分布随太阳方位变化而变化;腔内对流换热较弱,同一时刻最大温差约为71.3 K;腔壁热阻和发射率增加会削弱自然对流的强度。      

空间遥感器用环路热管瞬态数值模拟与在轨验证

为满足空间遥感器环路热管（LHP）在轨应用需求，建立了高分九号卫星电荷耦合器件（CCD）用LHP瞬态数值模型。模型采用了节点-网络法和流动与传热关系式耦合的方法，考虑了蒸发器与储液器之间的传热传质过程。通过仿真与在轨数据的对比，发现LHP内部组件温度偏差为0.2~0.4℃，冷凝器测点温度偏差为0.5~2.0℃；预热器通过干度的变化调节了冷凝器外热流和热源工作模式的影响；热源的工作模式对蒸发器向储液器漏热、回路流阻及两相段长度均有影响。所提模型可用于不同轨道外热流及热源工作模式下，研究LHP内部各参数的变化规律，预测LHP系统的瞬态工作特性，并指导后续产品的设计与研发。      

涡轮叶片尾缘复合通道的流动与换热

实验研究了涡轮叶片尾缘带纵向隔板内冷通道的流动与换热特性.采用相变加热的方法为实验提供了等壁温边界条件,模型实验件采用了扰流柱和微小通道结构,并且分别与直隔板和波形隔板进行组合实验,以期望得到综合换热效果最优的组合.为了获得通道的局部换热信息,实验件的一侧外壁面被细分成了10个集水区域分别进行研究.结果表明:当Re在10 000~60 000之间时,综合换热效果最差的是直隔板加光通道的组合结构,综合换热效果最佳的是波形隔板加微小通道的组合结构.     

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

轴向通流旋转盘腔流动换热稳态实验与数值模拟

通过开展稳态实验及数值模拟探究了轴向通流旋转盘腔的流动结构与换热特性。通过改变流量系数、旋转雷诺数等参数,探究不同工况下旋转盘两侧及盘罩内侧壁面温度和努塞尔数的径向分布规律。结果表明:在轮缘加热的状态下,旋转盘两侧温度径向分布均呈凹函数形态,且旋转盘迎风面换热强度普遍高于背风面;后轴颈盘罩向两端旋转盘导热,其壁面温度径向分布呈"中间高、两侧低"的状态;随着轴向流量系数的增大,盘腔内部气体对流加剧,径向臂及涡对结构更加明显,旋转盘及轴颈表面换热效果增强;旋转盘腔内的流动换热特性受强迫对流和类Rayleigh-Benard对流2种机理的共同影响。      

贮箱内低温推进剂汽化过程的CFD数值仿真

为研究贮箱内低温推进剂相变对推进剂温度和贮箱压力的影响,对贮箱内的传热传质过程进行了仿真.仿真涉及的物理过程包括贮箱与外界环境的换热、推进剂的自然对流、推进剂与贮箱内壁面的换热以及低温推进剂的相变过程等.根据热力学平衡原理建立了低温推进剂相变模型,使用CFD(Computational Fluid Dynamic)方法对处于地面常压停放状态的液氢贮箱进行了450 s的仿真.研究表明随着贮箱壁面传热过程的稳定,推进剂的温度分布、流动状态以及相变情况会趋于稳定;通过仿真获得了推进剂单位时间的汽化量;影响相变的主要因素是贮箱壁面漏热以及推进剂自身的对流运动.     

热管内部流动与传热问题的模化分析

针对毛细芯圆管热管内部的耦合流动与传热问题进行模化分析,建立了二维轴对称模型,并进行了数值模拟。与前人一些关于高温钠热管的试验结果比较表明,数值计算结果与试验结果偏差不大,而且它们的变化规律也相符得很好。     

矩形微槽乙醇水溶液传热特性数值模拟

由于影响微槽传热特性的因素较多,实验研究结果相互矛盾,使得实验研究和数值仿真相结合成为进一步深入认识微槽传热特性的有效手段.在实验研究的基础上,采用有限体积法对矩形微槽中体积浓度30％的乙醇水溶液传热特性进行了三维数值模拟. 矩形微槽采用两种不同对称结构和热边界条件,进口采用流动完全发展边界条件.对不同对称结构和物性条件(常物性和变物性)的计算结果进行了比较,并与实验结果进行了对比.研究表明,对所研究微槽结构和工质,采用流动完全发展进口条件,以及随温度变化的热物性参数,数值计算结果与实验结果基本吻合.