平板加热数值模拟

针对平极加热工程问题,在PATRAN环境中,应用热流计算公式,通过PCL的二次开发,实现了全参数设计的平板热试验数值模拟,以此进行结构热强度虚拟试验。通过改变灯管与反射板的距离和灯管与试验件表面的距离计算出半板试验件的表面热流。采用全参数设计,可以进行优化计算,优化结果对物理试验装置的设计、加热器的功率选取提供理论参考。     

高热流密度燃烧室壁面阵列空气射流-燃油组合冷却结构研究（爆震专刊）

针对高热流密度燃烧室壁面热防护需求提出了一种空气阵列射流冲击和燃油冷却肋板的集成冷却方式,在射流平均雷诺数(Re_j)为10000至30000、燃油进口流速(v_f)为2.33m/s至5.23m/s的范围内,采用数值模拟方法对其传热特性进行了研究,并基于壁面加热侧当量对流换热系数的概念,分析了基准肋板以及燃油冷却肋板的传热增强作用。与无肋板靶面的阵列射流冲击相比,带肋板阵列射流冲击的面积平均当量对流换热系数是前者的1.6倍,压力损失系数相对提高了约25%;采用燃油冷却肋板,加热壁面综合传热能力进一步增强,在Re_j=10000时,采用燃油冷却肋板的面积平均当量对流换热系数是基准肋板的1.5倍以上,即使在Re_j=30000时,燃油冷却肋板的传热增强比也可以达到1.2;燃油冷却肋板的出口温度相对进口温度的提升在20K~50K范围内,其提升幅度随着射流雷诺数或燃油进口流速的增大而减小。     

方形通道内超临界碳氢燃料传热恶化数值研究

为了深入理解火箭发动机再生冷却过程中碳氢燃料的流动传热规律,采用RNG k-ε湍流模型结合增强壁面处理方法对非对称受热(上壁面外侧加热)方形冷却通道内超临界压力甲烷的对流换热进行了数值研究。重点考察了加热壁面内侧的传热恶化以及由传热恶化和固壁热传导共同作用引起的热流异常传递现象,拟合得到了传热恶化的临界热流密度和起始内壁温关系式、加热侧内壁面和侧壁面内侧平均热流密度的预测关系式。结果表明:当外壁热流密度和质量流速比值高于0.288 k J/kg时,近壁流体流动加速诱发了加热侧内壁面的传热恶化;同时,固体区域温度畸变导致加热侧内壁面热流密度减小,热流更多地向侧壁面内侧传导。运行压力越低,该现象越显著。     

机翼电加热防冰表面内外传热的耦合计算

分析了机翼防冰表面的传热传质现象,对电加热防冰表面内外传热耦合计算进行了研究.采用附面层积分法对机翼表面传热系数进行计算.根据传热传质比拟得到了防冰表面蒸发量.建立表面微元的能量与质量守恒差分方程组,利用高斯-赛德尔迭代法求解,将外部防冰载荷、蒙皮弦向导热和内部加热热流三者耦合,采用亚松弛系数稳定迭代过程,最后得到平衡表面温度、溢流水和溢流结冰的范围.计算结果和文献中试验以及其他数值仿真结果比较,吻合较好,表明了模型及计算方法的正确性.      

格林函数法求解含移动介质轴向导热的稳态温度场

用格林函数法求解了考虑移动介质轴向导热影响的稳态温度场,推导了第1类和第2类非齐次边界条件下求解稳态温度场的数学表达式,利用特征值和特征函数获得了格林函数解,证实格林函数法可以解决非齐次边界条件下的稳态传热问题.以平行板间和圆管内均匀速度移动介质的稳态传热为例,获得了在半无限长壁面和在有限长壁面有均匀热流密度两种情况下所对应的温度场的无量纲解析解,分析了移动介质轴向导热与Pe和x/H之间的关系.结果表明:在所分析的两种热流密度边界条件下,Pe越大无量纲温度越小;在半无限大区域有均匀热流密度条件下,无壁面热流密度区域温度变化越剧烈,有均匀壁面热流密度的区域无量纲温度随x/H的变化几乎不受Pe大小的影响;在有限区域有均匀热流密度条件下,与半无限大区域均匀加热的情况有显著差异,在全计算域内温度变化越剧烈.     

气动加热下光学窗口传热特性的数值模拟

针对气动加热条件下光学窗口辐射导热耦合传热开展研究,以了解光学窗口的热辐射发射规律及热防护性能。采用有限体积法、蒙特卡罗法、谱带模型建立了具有镜反射界面的光学窗口辐射传输及辐射导热耦合传热计算模型,对具有两层SiO2玻璃的光学窗口传热特性进行了模拟。结果表明处于高温区玻璃向外发射热辐射,低温区玻璃吸收热辐射;除了加热面附近外的玻璃区域,辐射热流密度远大于导热热流密度,在2.9～4.2μm谱带内的辐射热流密度大于0.4～2.9μm谱带内的;在0.4～2.9μm谱带内,可见光波长内的出射辐射热流密度最小,中红外波长范围内出射辐射热流密度最大。在光学窗口总厚度一定的情况下,降低玻璃层厚度和玻璃夹层气体热导率能够有效降低光学窗口非加热面稳态温度。     

燃油冷却面板传热特性试验与计算分析研究

采用热电偶测温、气壁红外测温及燃油样品裂解度测量等多种手段,在DJ-21电弧加热器上进行了燃油冷却面板传热特性试验。进行了共计19次燃油冷却面板传热特性试验,试验高状态对应平均热流为1.6MW/m2,低状态对应平均热流为1.1MW/m2;用于冷却的燃油质量流率为1.84～5.8g/s。为了反映冷却面板热流密度分布,以喷管三维流动计算结果作为输入条件,将计算得到的热流密度与试验测量的冷壁热流密度比较,用以确定流场计算方案、流场切取方案和热流密度计算方案。发展了冷却面板稳态准三维热分析程序,将等效热流对应的冷壁对流换热系数和燃气总温作为高温燃气侧的边界条件。使用热分析程序完成了相应的计算。通过试验与计算数据对比研究,表明热分析计算的可信性。试验验证了冷却面板的设计与加工是可行的。     

狭窄空间内沸腾的几个特征

将石英或瓷质板用电热膜技术加工成加热板。将两块加热板对称地放入饱合氟利昂－１１３（Ｒ－１１３）中形成竖直窄通道。对通道内沸腾的始沸点、流型和传热规律等基本特性进行实验研究。结果表明：始沸热流随加热板通道间距减小而降低;沸腾流型与加热过程、通道间距等因素有关;存在最佳传热通道间距;最后提出小通道内始沸热流、传热规律的半经验关系式。     

高热流密度燃烧室壁面阵列空气射流-燃油组合冷却结构研究

针对高热流密度燃烧室壁面热防护需求，提出了一种空气阵列射流冲击和燃油冷却肋板的集成冷却方式，在射流平均雷诺数Re_j为1×10⁴~3×10⁴，燃油进口流速v_f为2.33～5.23m/s内，采用数值模拟方法对其传热特性进行了研究，并基于壁面加热侧当量对流换热系数的概念，分析了基准肋板以及燃油冷却肋板的传热增强作用。与无肋板靶面的阵列射流冲击相比，带肋板阵列射流冲击的面积平均当量对流换热系数是前者的1.6倍，压力损失系数相对提高了约25%；采用燃油冷却肋板，加热壁面综合传热能力进一步增强，在Re_j=1×10⁴时，采用燃油冷却肋板的面积平均当量对流换热系数是基准肋板的1.5倍以上，即使在Re_j=3×10⁴时，燃油冷却肋板的传热增强比也可以达到1.2；燃油冷却肋板的出口温度相对进口温度的提升在20~50K内，其提升幅度随着射流雷诺数或燃油进口流速的增大而减小。     

类乘波体防热壁板气动加热-温度场耦合特性分析

针对类乘波体防热壁板,采用松耦合算法推进气动热和传热迭代计算,研究气动加热-温度场的耦合特性.结果表明:在壁板达到辐射平衡前,冷壁热流和辐射平衡热流与真实气动加热的误差分别达+55.1%和-15.4%以上,必须将气动加热和温度场进行耦合分析;当地绝热壁面温度不随时间变化,表面传热系数是耦合效应的关键参数;采用平均表面传热系数进行瞬态气动加热-温度场耦合计算只进行2次气动加热计算,壁面温度预测误差在2.5%以内,可有效提高气动加热-温度场耦合计算的效率.      

气动加热下飞行器油箱温度场耦合分析

气动加热与结构热传递耦合问题在航空航天领域非常重要.分析了飞行器持续受热过程中的耦合性,指出了耦合性在分析持续气动加热问题中的必要性.针对存在的耦合性,提出了一种处理方法,即在“飞机结构三维温度场分析”程序中增加气动热流工程算法,同时采用对流传热来模拟油箱内燃油与油箱内壁面传热.对流换热系数以及出口油温通过与油箱壁面结构温度耦合计算来确定,即在每个时间点上,将当前结构内壁面温度咒代入对流换热系数α以及油温L的计算公式,求得该时间点的α和t,以此实现气动/温度场/流场三个物理场的实时耦合.计算结果表明,耦合计算结果与试验结果吻合得很好.     

侧壁开孔管式点火器点火热流与固体药柱表面间传热

分析了固体火箭发动机侧壁开孔长管式点火器点火热流与固体药柱表面间传热机理，根据轴对称边界层原理，借鉴平面间阵列圆孔射流传热模型，在考虑管式点火器与药柱表面同轴圆面间射流传热过程的基础上推导得到传热系数表达式，仿真研究表明模型合理，能对试验现象进行合理解释。     

防冰部件表面流动换热与温度计算分析

在考虑防冰部件湿表面水膜流动与传热相互影响的基础上,提出了一种基于湿表面传热与流动耦合作用的防冰部件表面温度预测方法.应用此方法,计算了一种全蒸发防冰翼型工作时的表面流动换热和温度分布.分析了湿表面上各项热流的作用效果.结果表明:部件表面换热系数对壁面温度的计算最为敏感,影响很大:水膜流动能吸收并传递加热热量;全蒸发防...     

一种用于高分辨扇扫成像的MIMO阵列

为了提高扇扫成像系统的方位分辨率,设计了一种具有高方位分辨率的多输入多输出(MIMO)阵列。首先,建立了同时适用于窄带和宽带信号的MIMO阵列与虚拟阵列之间的等效模型。在分析匹配滤波输出成分的基础上,证明了当发射信号具有良好相关特性时,MIMO阵列的匹配滤波输出可等效为单输入多输出(SIMO)阵列上经过脉冲压缩后的回波。接着,根据这一等效关系,获得了与MIMO阵列等效的虚拟阵列中虚拟发射阵元与虚拟接收阵元坐标解析解。最后,在约束阵列物理尺寸为固定值的前提下,利用该解析解设计出一种由2个发射阵元和一条接收均匀线阵列(ULA)组成的MIMO阵列。理论推导和仿真结果表明:所设计的MIMO阵列可以在恒定的物理尺寸和虚拟阵列为均匀直线阵的前提下,获得优于其他MIMO阵列和SIMO阵列的方位分辨率。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

钝体外形气动加热与结构传热一体化数值模拟

为了提高气动加热与热防护结构传热多物理场数值模拟的稳态计算效率与计算精度,发展了一种基于有限体积法的气动加热与结构传热一体化数值计算方法。该方法将高速流场与结构温度场统一到同一物理场,基于统一的控制方程组,采用基于LU-SGS隐式时间迭代和自适应时间步长的有限体积方法进行求解,避开了传统气动加热与结构传热耦合求解方法在时间域内的所需繁琐数据交替迭代策略。对二维/三维钝体进行一体化数值计算分析,计算结果表明:二维钝体非稳态下,得到2s时圆管驻点温度最高达到390.2K,驻点热流密度和结构温度与参考文献和实验值吻合较好,证明了方法的可靠性和可行性。同时分析了三维钝体应用算例的流-固-热稳态计算特征,计算得到稳态时钝头体结构外壁表面最高温度达到535.6K,表明一体化计算方法可用于长航时飞行条件下的气动加热-结构传热多物理场耦合计算分析,为高速飞行器热防护结构设计与选材提供一定的理论与技术支持。     

超临界压力下航空煤油流动与传热特性试验

本文针对超燃冲压发动机再生冷却系统运行条件,实验研究了大庆RP-3煤油在超临界压力下的流动和传热特性,目的在于研究煤油在各种工况下的流动参数变化以及对流传热规律。煤油通过二级煤油加热/输运系统加热,试验的煤油压力约2.6M～5.0MPa,油温约300～800K。相应的壁面热流密度为10～300kW/m2。通过油温与壁温的同步测量,结合非定常传热分析,获得了超临界压力下、亚/超临界温度范围内煤油的流动和传热物理参量的变化曲线。     

再入飞行器鼻锥逆向喷流对流场及气动热的影响

 使用计算流体力学（CFD）方法研究逆向喷流热防护系统对降低再入飞行器鼻锥物面热流的效果,获得了流场参数,回流再附点位置,物面压力分布以及热流分布。分析了逆向喷流对降低物面热流的物理机理,喷流通过与来流相互作用形成马赫盘,将来流导流到四周,不与物面直接作用形成气动加热,同时喷流回流形成低温区,降低物面与接触气体的温差,进而降低了物面热流。随着总压比率增大,这种效果越明显,气动加热越轻。为更合理分析喷流强度对流场及传热量的影响,将总压比率和流量相结合,提出了新的参数R PA。分析该参数的应用效果,结果发现不同的流量与总压比率组合成相同的参数R PA,可以实现相同的激波位置、再附点位置、表面热流峰值位置和总传热量。这说明该参数可用于表征喷流强度,用以分析喷流对流场及传热量的影响。     

气气喷注器推力室传热特性研究

设计了单喷嘴气气喷注器容热式推力室,进行了0.92～6.1 MPa范围内7个燃烧室压力工况,共17次热试车;采用壁面测温方法获得推力室沿轴向方向的内壁热流分布,得到各室压工况下燃烧室内壁热流分布曲线,不同工况的热流曲线显示出相似的分布;并从中分离出对流传热热流,得到对流传热热流与室压的关系。为拓展范围,采用多组分湍流N-S方程描述推力室内燃烧流动,采用6组分9步反应模型来描述氢氧反应,反应速率由Arrhenius公式计算,进行了5～20 MPa更高室压范围内的燃烧内流场的数值模拟,并耦合计算了各工况燃气与室壁之间的传热,获得了与试验相同规律的结果。     

石英灯阵模拟非均匀气动加热的功率优化

针对结构热试验非均匀热环境模拟的精细化需求,以石英灯阵热流分布理论解析方法为基础,建立了基于线性分析模型的石英灯阵功率快速优化方法,对石英灯阵模拟非均匀气动加热的功率优化问题进行了分析研究。结果表明：线性分析方法可以快速有效地优化石英灯阵功率,对于比较平缓的热流分布可以获得较好的优化结果,对于大梯度热流分布也能获得与遗传算法较为一致的热流优化结果和规律;模拟非均匀气动加热的功率优化结果对热流输入比较敏感,需要通过优化来选择合适的测点范围和测点数量得到相对稳定的系数矩阵。相关研究结果对模拟非均匀气动加热的石英灯阵功率优化设计具有重要参考意义。