竖直平板Blasius流层流边界层流动与传热耦合解

对考虑辐射传热、流体黏度随温度变化、有和无滑移边界条件下的可渗透竖直平板Blasius流层流边界层的无量纲速度场与温度场进行了深入研究.经相似变换将描述速度场与温度场耦合的偏微分方程组转换成非线性常微分方程组,用Runge-Kutta法对常微分方程组进行了数值求解.研究了无量纲参数对无量纲速度场及温度场的影响,着重分析了滑移边界条件下速度和温度随无量纲参数的变化规律.结果表明:吸入时边界层变薄,喷注时边界层加厚;对比于无滑移边界条件,滑移边界条件下速度、温度边界层变薄;随着变黏度参数a或喷注与吸入参数的增大,壁面摩擦因数、局部努塞尔数Nu增大,速度和温度边界层变薄;随着普朗特数Pr、热辐射参数R的增加,或毕渥数Bi,布林克曼数Br的减小, 温度边界层变薄.      

基于多孔介质局部非热平衡方法的刷式密封耦合传热特性

建立了刷式密封流动传热与摩擦接触的多物理场耦合分析的计算模型.推导了描述刷式密封在局部非热平衡条件下耦合传热过程的双能量方程,该方程既考虑了刷丝与气流的对流换热,同时也考虑了刷丝束内导热过程的各向异性.采用所建立的耦合计算模型数值预测了刷式密封在不同压比和转速下的温度分布.结果表明:刷式密封的摩擦热效应使刷丝与转子相接触的附近区域形成局部高温,刷丝的高温区域主要集中在围栏高度以下;泄漏气流与刷丝间强烈的对流换热使刷丝的温度分布和气流的温度分布非常相近;刷丝的最高温度随压比增大而升高,但围栏高度以下区域的平均温度随压比增大而降低;刷丝的最高温度和围栏高度以下区域的平均温度均随转速升高而升高;摩擦热效应引起的泄漏气流的温升效应使刷式密封泄漏量随转速升高而减小.      

超临界压力正十烷对流传热实验及计算研究

针对主动冷却超燃冲压发动机的实际工作条件,利用电加热管设备开展了超临界压力正十烷流动和传热特性实验,目的在于获得适用于发动机传热设计的燃料对流传热关联式。管子材料为不锈钢、内径1．5 mm、外径3．0mm、长度1300mm。采用K型热电偶测量管壁沿程外壁温。实验中正十烷压力约4．0～4．3MPa,温度约335～870K,流量分别为0．93、1．24和1．86g/s。正十烷的流动经历了层流、过渡和湍流3种流态。采用最小二乘法拟合实验数据,获得了正十烷在3种流态下的传热关联式。通过外壁温计算值与实验值的对比,验证了本文给出的传热关联式的适用性。     

高超声速多体分离过程气动加热计算技术

发展了无粘外流解与气动热工程方法相结合的用于全机外形复杂流动气动加热计算技术,并将其与动态多体分离过程流场解相结合,提出了一种高超声速飞行器多体分离过程气动加热特性计算技术,然后耦合结构传热,实现了飞行器热防护系统中结构材料温度分布特性的数值模拟。采用本文方法,针对假定高超声速多体分离布局及飞行条件开展了分离过程中的结构温度分布特性与热流密度分布特性计算,分析并给出了组合体分离过程中部件干扰对气动热的影响特性。结果表明本文发展的计算技术可为高超声速飞行器的防热方案设计、气动热特性定性分析及热防护系统设计等提供技术支持。     

超临界态碳氢燃料流固耦合传热及热裂解的计算方法研究

为深入研究主动再生冷却机理,采用固体热传导方程和流体控制方程,建立了超临界态碳氢燃料耦合传热及热裂解的数值方法。对两种实验状态下电加热管内RP-3航空煤油的耦合传热及热裂解过程进行了计算研究,其中RP-3的热裂解反应由包含1个总体反应和23个二次反应的Modified Kumar-Kunzru反应机理模拟。两种状态下计算得到的外壁温、燃料温度、RP-3裂解率以及裂解产物分布与测量值的最大误差分别小于10%和15%。计算与实验结果吻合较好,表明建立的数值方法是准确的,可以作为研究超临界态碳氢燃料耦合传热及热裂解特性的有效计算分析工具。     

轴承腔油滴含率及油滴相与空气能量传递分析

为了改善迄今轴承腔油滴运动分析未包含油滴与空气热量传递以及油滴相/空气介质均相处理中油滴含率确定方式不够准确的不足,借助单个油滴运动分析和油滴尺寸分布提出了轴承腔油滴含率及油滴相与空气能量传递分析方法。首先在考虑单个油滴与空气对流换热条件下,将油滴能量方程嵌入运动方程,并同步离散求解油滴运动方程和能量方程,实现了单个油滴运动速度和温度的联立求解,通过与试验结果对比证明考虑温度效应有助于提升油滴速度计算的准确性;然后确定了轴承腔油滴尺寸分布,并离散油滴尺寸分布范围以及轴承腔流场空间,进行油滴运动分析确定了轴承腔不同径向位置处油滴的体积和质量含率及油滴相与空气的动能和热能传递量。构建了圆盘腔油滴含率试验台,开展了油滴体积含率的试验测量,理论计算值与试验值较为吻合,验证了轴承腔油滴含率分析方法的有效性。     

具有针肋的狭窄空间冲击冷却实验和数值计算

对具有全高度针肋扰流的狭窄空间冲击冷却进行了实验和数值计算,并与平板靶板冲击冷却传热性能进行了对比分析.射流冲击雷诺数范围为15000~30000.实验采用瞬态液晶热像技术获得了冲击靶板上详细的传热分布,并通过数值计算获得了冲击冷却系统中的流场和传热特征.实验研究表明:狭窄空间冲击冷却中的针肋靶板端壁上的平均传热性能比平板靶板提高约7.0%,压力损失提高约17.9%,并且针肋改善了靶板端壁上传热均匀性.另一方面,数值计算分析表明近壁面射流以及空间中的上洗涡流与针肋表面发生强烈相互作用,并且针肋显著地增加了换热面积,因此具有针肋扰流的冲击冷却系统具有显著增强的总体传热性能,比平板冲击冷却提高约27.0%.      

先进旋涡燃烧室多场协同分析

为了对先进旋涡燃烧室(AVC)的流动和传热性能进行综合评价,对不同来流速度、来流温度、壁面温度下AVC燃烧性能进行了数值模拟,对速度场、温度场、压力场进行了多场协同分析。结果表明:随着燃烧室来流速度的增大,速度与速度梯度的协同角α、速度与温度梯度的协同角β,速度与压力梯度的协同角θ均减小,温度梯度与速度梯度的协同角γ及压力梯度与速度梯度的协同角φ增大。随着来流温度的提高,α,θ增大,β,γ,φ减小。随着壁面温度的提高,α,θ减小,β,γ,φ增大。对于AVC湍流流场的传热强化问题,增大来流速度和来流温度,降低壁面温度会强化流动换热;增大来流速度和壁面温度,降低来流温度可以减少传热功耗;增大来流速度和壁面温度,降低来流温度能够提高强化传热的综合性能。     

非正交界面下的插值耦合传热方法

针对复杂热防护结构,开发了一种求解非正交性流/固界面的耦合传热程序。流体和固体区域采用同一积分、守恒型的RANS(Reynolds-Averaged Navier-Stock)方程,通过有限体积法进行离散求解。为了保证非正交界面上的温度和热流密度连续,提出了一种结合网格周边信息计算界面热流密度的插值方法。利用该插值耦合方法模拟了双层复合平板和喷管的热传导。数值结果表明:该插值方法在空间上具有2阶精度;喷管壁面上的对流换热系数沿轴向先增大后减小,在喉部上游达到最大值,当喷管入口压力增加3.38倍时,对流换热系数的最大值相应增加了3.13倍;喉衬与壳体界面上温度的计算值和试验结果存在一定差异,这是由于本文数值计算未考虑接触热阻引起的。     

液氧甲烷发动机台阶型冷却通道的耦合传热特性

为了研究液氧甲烷发动机再生冷却通道中跨临界甲烷的流动和传热特性,以及冷却通道较大幅度的突扩突缩对冷却效果的影响,采用整场直接耦合的方法对推力室三维耦合传热进行了数值模拟,考虑了燃气的非平衡流动.通过计算得到了推力室三维温度场和流场.计算结果表明：由于喉部截面附近存在较强的二次流,燃气侧壁面温度的最大值出现在喉部上游.由于突扩突缩处存在较强的旋涡运动,冷却剂的湍流强度增强,冷却剂侧表面传热系数显著提高,燃气侧壁面温度出现局部极小值,同时也产生了较大的局部损失.由于铜内衬热阻比镍外套热阻小得多,从燃烧室进入的大部分热量在冷却通道底面和侧面被冷却剂吸收.冷却通道底面的温度和热流密度沿程变化比顶面更剧烈.