R141b在圆形微通道内的沸腾换热实验研究

为提高换热强度、解决设备内部高热流密度散热问题,采用实验方法研究R141b在不同直径(D=0.5mm和1.0mm)水平圆形微通道内的沸腾换热特性,分析了热流密度(q=2.0kW/m~2～47.6kW/m~2)、质量干度(x=0～0.6)、质量流速(G=111.11kg/(m~2·s)～333.33kg/(m~2·s))的变化对平均传热系数h的影响,探究不同情况下影响沸腾换热的主导因素。实验研究表明:平均传热系数h随热流密度q的增加而减小,在不同范围内减小速率有明显差异;热流密度q=2kW/m~2～5kW/m~2时质量流速G对平均传热系数h影响较明显,热流密度较高时质量流速G对换热影响很小;在质量流速G=111.11kg/(m~2·s)～333.33kg/(m~2·s),质量干度x0.3时,平均传热系数h随质量干度x增加而明显下降,在设计微通道换热器时应尽量使R141b处于初始沸腾阶段以获得更好换热效果,并采取一定措施预防干度过高引起的换热恶化。     

扁管微小通道内R134a的沸腾换热试验

搭建适用于多种结构微小通道的沸腾换热试验系统,研究了制冷剂R134a在当量直径分别为0.63mm和0.72mm的多孔扁管微小通道内的沸腾换热特性。试验参数包括制冷剂质量流率为82~621kg/(m2·s),饱和压力为0.22~0.63MPa,干度为0~1;采用等热流密度方式加热,热流密度范围为9.7~64kW/m2。结果表明:R134a在扁管内沸腾换热中,当干度在0~0.6区间时,微小通道的传热系数明显高于常规通道,换热类型主要为核态沸腾,传热系数随热流密度和饱和压力的增大而增大,与质量流率关系不大;当干度大于0.6之后,传热系数随着干度的增大急剧减小,且在此干度区间,传热系数受热流密度和饱和压力影响较小,而受质量流率的影响相对较大。利用该结论和公开文献中R134a沸腾换热试验数据对Gungor-Winterton公式进行改进,改进后的公式对所有试验点的平均相对误差为-1.17%,平均绝对误差为19.24%,预测精度有了明显提高。      

R141b在矩形微尺度通道中的两相流传热特性

设计搭建水力直径分别为1mm和0.5mm的矩形微尺度通道实验台,研究了以R141b型制冷剂作为工质的两相流沸腾传热特性。实验取热流密度为1～16kW/m2、质量流速为111.1～333.3kg/(m2·s)和质量干度为0～1,分析了三者对平均传热系数的影响,探究影响换热的主导因素。结果表明:热流密度较高时,平均传热系数随热流密度增加而减小,流动换热主要受到沸腾传热的影响;当质量流速较大且热流密度较低时,平均传热系数随热流密度增加而有所增长;热流密度较低时,平均传热系数随质量流速变化明显,热流密度升高到一定值后,质量流速对平均传热系数的影响很小;当质量流速处于111.1～333.3kg/(m2·s)时,平均传热系数随质量干度的增加而减小。      

弯曲管道内湍流流动的数值模拟

为了对弯曲管道内的湍流流动特性进行有效的预测，用显示代数应力模式对－典型的弯曲管道内流动进行了数值模拟和分析，计算结果表明：（1）显式代数应力模式能对弯曲管道内流动的平均量（如速度、压力、摩阻等）进行较为准确的预测；（2）对模式中的ε方程做适当修正后，该模式能较好地模拟出流动曲率效尖对湍流特性的影响，是对雷诺应力模式的一个比较好的近似。此外，该模式在形式上与标准k－ε两方程模式类似，使用方便，具有较强的工程实用价值。     

微细管道内流体流动的黏性耗散分析

目前,还没有实验数据和理论分析可以得出在微细管道内流体的黏性对流体流动及换热的影响,而在微细管道(微米级)内流体高速流动中黏性耗散是不能够忽略的,但是没有合理的计算公式进行理论计算。因此,这里对粘性耗散在微细管道流中的影响进行了分析和理论验证。根据已有的分析给定了一个粘性耗散在微细管道流中造成明显影响的最小标准值。     

三维管道超声速气固两相流动的CFD/DSMC仿真

采用NND格式求解气相N-S方程，采用DSMC方法模拟固相颗粒运动，数值模拟了三维管道超声速气固两相流场，分析了固体颗粒对流场结构产生的影响。计算结果表明，将DSMC方法与CFD方法相结合是解决气固两相流动问题的一种有效途径。     

旋转楔形通道流动与换热特性的数值研究

为了模拟涡轮叶片尾缘冷却通道的流动换热特性,建立了带侧向出流的旋转楔形冷却通道的数值模型,采用k-ωSST两方程模型研究了流动与旋转对于楔形通道内换热分布与不同出口冷气出流比例的影响。在静止情况下,由于沿程侧向的冷气出流,换热沿程降低。而旋转压缩了通道顶部内侧区域的流动滞止涡,导致内侧区域的换热显著增强。静止情况下沿程冷气出流占比不断提高,但雷诺数的影响有限。随着旋转数的提高,通道冷气出流比例受到旋转离心力的主导作用,前半段冷气出流比例下降,后半段出流比例上升,在50%径向高度位置为平均值。在大旋转数下(Ro0.4),通道前半段冷气出流比例逐渐下降为0。通道安装角度的变化对冷气出流比例的影响有限。     

基于气体射流的气液两相流动减阻特性

为了明确射流表面减阻特性,基于气液两相流动原理,建立具有射流功能的气液两相流动减阻计算模型,采用数值计算方法,通过RNG k-ε湍流模型对射流表面不同射流方向角、以及不同气体射流速度情况下的流场进行数值模拟及实验验证,分析气体射流对壁面黏性阻力、减阻率及流场的影响,研究气体射流的气液两相的减阻特性。结果表明:气体射流表面能够减小壁面黏性阻力,具有较好的减阻效果;四种射流方向角下,射流表面减阻率均在气体射流速度为5m/s时达到极大值,并且当射流方向角在xoz平面为-30°时,减阻率最大,为7.18%;数值模拟的减阻率效果略高于实验结果;气体射流的气液两相流动对壁面边界层进行了有效控制,减小了壁面的黏性剪应力和雷诺应力,导致壁面的黏性阻力减小,使得射流表面具有减阻性能。     

摆盘流动与换热性能数值仿真

针对受燃烧室出口高温燃气冲刷的摆盘装置的冷却水回路,建立二维、三维物理仿真模型,模拟了冷却水在内部流道的 流动与换热过程,考察了冷却水进口压力、冷却水流量对流阻和换热性能的影响.结果表明:①摆盘冷却水进口压力由2.3×105Pa提高到8.3×105Pa,摆盘壁温变化微小,结构1壁温升高5K,结构2壁温升高7K;②冷却水进口速度由0.5m/s提高到5.3m/s,结构1壁温降低约120K,结构2壁温降低约100K,冷却效果明显;③结构2通过缩小流道的流通面积,能在更小的冷却水流量的工况下得到更好的换热效果.      

水中弹药引信外置涡轮高速旋转两相流动数值仿真

为预测空化条件下的涡轮转速,获得涡轮叶片附近空化情况、压力系数、阻力系数等数据,采用SchnerrSauer空化模型,进行了水中弹药引信涡轮高速旋转两相流动的三维流场数值仿真研究。计算得到不同水深(10m、30m及50m)和弹速(5~100m/s)条件下的流场数据。结果表明:引信涡轮表面在弹速30m/s左右时空化初生,且空化尺度随着弹速的增加而增大;弹速100m/s时,头部粘性阻力占总阻力的18.7%;涡轮转速与水深、静压、弹速相关,高速时涡轮转速随水深减小明显,与弹速基本成正比线性关系,非线性度约为0.998。     

不同流动配置下涡轮叶片交错肋冷却流动传热特性数值模拟

为掌握交错肋冷却结构应用在涡轮叶片不同区域的流动传热性能,针对一种交错肋冷却结构在三种不同流动配置中在等质量流量和子通道雷诺数工况下进行了数值计算研究,三种流动配置包含了径向流动配置（RFC）, 横向流动配置（CFC）和转折流动配置(TFC)。通过比较本研究得到的数值模拟结果与公开文献中的实验数据,定性定量地验证了本次数值计算的有效性。在等冷却质量流量下,RFC配置拥有最高的平均努塞尔数和压力损失,而CFC和TFC配置的平均传热性能相似且明显降低,但压力损失大大减少。在相同的子通道雷诺数下,三种流动配置下的交错肋通道展现出相似的传热强化性能,但TFC配置的压力损失最小。在研究范围内,在RFC配置中肋表面的平均换热比主表面的平均换热约高出16.3%,而在CFC配置和TFC配置中该值则分别高出38.2%以及 30.6%。不同的流动配置会引发子通道内不同的流动特性,包括流动转折和子通道间的交互作用。     

Numerical Study of Flow and Heat Transfer Characteristics of Impingement/Effusion Cooling

Three-dimensional numerical simulation is carried out to investigate the flow and heat transfer characteristics of impingement/effusion cooling systems. The impingement/effusion holes are arranged on two parallel perforated plates respectively in a staggered manner. Every effusion hole has an inclined angle of 30° with respect to the surface. The two parallel plates are spaced three times the diameter of the effusion hole. The ratio of center-to-center spacing of adjacent holes to the diameter of the effusi...     

壁面导热对微型向心涡轮性能及流动影响

为了探讨传热对于微型向心涡轮的性能和流动的影响规律,根据微型向心涡轮的特点,针对涡轮的轮毂侧壁面传热现象,建立了涡轮的气动性能随传热量变化的一维模型,并利用数值模拟对模型进行了验证,分析了不同传热量下微型向心涡轮的流动变化。结果表明所发展的模型有较好的精确度,可以用以评估传热对微型向心涡轮性能的影响;较大的传热会改变转子进口相对气流角并影响转子整个通道内的流动;涡轮对外传热使得涡轮输出功降低,涡轮效率明显降低。      

压电风扇激励非定常流动和换热特性数值研究

 利用动网格技术对压电谐振风扇产生的非定常流场进行了数值模拟,以期进一步揭示压电风扇的流场特征和换热特性。研究表明:压电风扇上下两个区域均出现涡流,涡向相反(压电风扇上方为逆时针,下方为顺时针),涡对的尺度、位置和扰动范围随时间呈周期性变化规律;时均速度并非随着与压电风扇自由端距离的增大呈现单调衰减的趋势,而是在距离压电风扇自由端距离为一倍振幅的截面上出现峰值速度最大的速度分布型,该位置正是涡环达到最大型面的瞬间涡核所处位置。涡串发展、运动过程中与周围流体发生干涉融合形成的射流起到了强化换热的效果,换热效果最好的地方不是出现在平衡位置,而是涡对破碎、流体紊流度最强的地方。     

超临界压力下正十烷流动传热的数值模拟

为了深入理解主动再生冷却过程中碳氢燃料的超临界传热特性,基于SIMPLE算法建立了数值模拟方法,考虑了碳氢燃料物性随温度的剧烈变化,并利用电加热管实验结果验证了计算方法。针对超临界压力下细管道内正十烷的流动传热现象进行了系统的数值计算研究,考察了计算网格无关性和超临界流动传热过程中的压力效应。结果表明:网格选择与正十烷的状态有关;在超临界压力下,较低的正十烷压力引起临界温度附近的努赛尔数减小,导致传热效率下降;目前常用的传热经验公式在正十烷临界区域附近与数值计算结果差别较大。     

叶片扩压器内两相流场的数值模拟

湍流模型采用 kk-εk-θp 模型 ,为改善计算的稳定性和收敛性 ,根据 Rhie,Chow提出的动量插值概念 ,以协变物理速度分量作为求解变量 ,在压力修正方程中考虑了密度及浓度修正的影响 ,完成了叶片扩压器内部的两相流场的数值计算 ,计算结果与实际情况较为吻合 ,表明了所建立的计算方法的可行性。结果表明 :叶片扩压器吸力面上的颗粒拟温度高于压力面上的颗拟温度 ,而且靠近壁面处的大于叶道中的 ,可采用颗粒拟温度来标志磨损部分。这为改进扩压器的设计和性能提供了依据     

有预旋进气转静盘腔中的流动和换热特性数值研究

运用RNG湍流模型对有静盘外缘预旋进气和轴向中心进气、转盘外缘轴向出气和外围屏径向出气的转-静盘腔中流动和传热过程进行了数值研究。在盘腔结构中,静止内隔片和旋转内隔片将盘腔分为内转-静盘腔室和外预旋腔室。由于盘腔流动结构具有周向周期性,取盘腔的1/30作为计算域。研究发现:旋转雷诺数Re_θ、预旋喷嘴进气无量纲流动速度C_w,_p、静盘中心轴向进气的无量纲流动速度C_w,_d和预旋比β_p都不同程度地影响盘腔内流动和换热。预旋和隔片的协同作用可以很好地起到阻止热燃气入侵转-静盘腔和提供低静温冷却气体流入涡轮叶片冷却通道。      

高转速转静盘腔流动换热相似性数值研究

为了得到高转速转静盘腔流动换热实验中的相似准则,在定几何、变物性、可压缩和有耗散的情况下,对转静盘腔中的控制方程进行无量纲化,得到了除无量纲位置外的8个无量纲准则数.并以高转速转静盘腔模型为基础,采用数值方法对影响发动机工况(高温)和实验工况(常温)换热相似性的主要无量纲准则数进行了研究.研究表明:流动雷诺数、旋转雷诺...     

歼击机座舱空气流动和传热的数值模拟与实验

根据歼击机座舱内空气流动和传热的特点,建立了座舱内壁面传热边界条件计算模型和空气分配系统供气边界条件计算模型,在此基础上建立了歼击机座舱流场和温度场数值仿真平台,并利用地面模拟实验验证了仿真平台的有效性。计算结果与实验结果的对比表明：该平台能较真实地反映座舱壁面传热的不均匀性,以及空气分配系统供气孔口流出气流流量和温度的不均匀性与非对称性;速度场和温度场计算结果与实验结果的误差分别约为15%和6%,这证明了该平台具有较高的模拟精度和工程应用价值。