考虑叶片径向和垂直于壁面方向导热的涡轮叶片对流冷却模型研究

为了提高涡轮叶片对流冷却模型预测精度,提出了一种在叶片固壁内同时考虑叶片径向和垂直于壁面方向(法向)导热的二维对流冷却模型。该模型在弦长方向划分多个元素,忽略元素内弦长方向叶片温度变化,在元素内的径向和法向建立二维导热方程作为叶片固壁温度场的控制方程,其边界条件包括叶表燃气绝热温度、燃气侧对流换热系数和叶片叶根、叶顶热流密度等。给出了该模型二维导热方程和边界条件的差分求解方法。以E~3涡轮高压导叶为例,将模型与CFD计算的叶片外壁面温度分布进行了对比。结果表明,该模型在给定冷气量下预测的叶片温度分布变化趋势与CFD相近,最大温度误差不超过6.5%,计算时间与CFD相比缩短了95%,能够快速、准确预测涡轮对流冷却叶片的冷气需求量。     

砂型材料导热系数的测试研究

本文在实际金属浇注法的基础上,运用参数优化技术设计了求解砂型材料导热系数的通用程序,测定了几种常用砂型材料的表观导热系数。实验结果表明:干砂型的表观导热系数是其密度、温度的函数;湿砂型水分凝聚区的表观导热系数是干砂区表观值的3～4倍。本研究为砂型材料导热系数的测试提供了一种新方法。     

金属蜂窝夹芯面板有效导热系数的数值计算

对具有漫灰体内壁的金属蜂窝夹芯面板,忽略蜂窝内腔的空气导热,综合考虑蜂窝结构的热传导和热辐射等热传递过程,利用有限元方法求解了周期性分布的蜂窝单胞稳态热传导控制方程,蜂窝内壁的边界条件是由净热量法得到的热辐射换热积分方程,由胞元的温度场分布数据及Fourier定律得到了蜂窝结构的有效导热系数。与现有文献相比,采用了较少近似的模型及较高精度的离散方法,计算结果表明计算模型及方法是可靠有效的。     

涂层材料瞬态传热实验与导热系数辨识

针对涂层-基体一体化的双层结构,为测试评估其中涂层材料的导热性能,提出基于瞬态平面热源法（transient plane source, TPS）的涂层材料导热系数反演辨识方法。根据Hot-Disk实验测试原理,建立基体-涂层-探头整体的二维非稳态传热模型;结合测量过程中的瞬时温升数据信息,采用粒子群优化算法反演辨识获得涂层材料的导热系数;并通过实验和数值模拟论证了上述方法的可靠性。结果表明：该测量方法能够有效获得涂层导热系数,测试反演的数值偏差小于4.0%。最后,实际测量和反演辨识获得了一种涂层材料常温至773 K的导热系数,随温度提高呈现增大趋势,数值范围为0.18～0.29 W/(m·K)。     

多孔介质有效导热系数的计算方法

高孔隙率多孔介质如泡沫陶瓷在新型多孔介质燃烧器技术中应用日益广泛,其重要传热特性参数--有效导热系数反映了两相流气、固相导热、对流和辐射的综合效应,对其研究尚非常缺乏.本文基于实验测定的温度分布,给出多孔介质有效导热系数的初始估值,用有限体积法求解二维控制方程,采用二维寻优搜索的办法,确定使测定点上测量与计算温度均方根误差为最小的径向与轴向有效导热系数,是一种逆计算方法.对球粒子颗粒床进行的有效性试验证明了方法的可行性.     

基于细观颗粒夹杂模型的固体推进剂导热系数预测

为准确预测固体推进剂的导热系数,采用了分子动力学方法对高填充比固体推进剂细观模型进行建模.利用有限元理论对固体推进剂细观模型稳态热传导进行计算,结合计算细观力学均匀化方法,计算了固体推进剂细观模型的均匀温度和热流密度.根据计算所得的平均温度和平均热流密度值求解稳态热传导方程,得到了两相和三相固体推进剂的宏观等效导热系数.其中,两相固体推进剂的仿真结果与试验结果的误差只有3.63%.结果表明:在预测固体推进剂导热系数时,采用固体推进剂的颗粒夹杂模型可以充分考虑粒径大小、颗粒随机分布的影响,更真实地反映固体推进剂的微结构特征,计算结果准确可靠.      

基于物理神经网络的导热问题无网格计算方法

建立了物理神经网络（PINNs）求解导热问题的通用框架，描述了三维非稳态问题、初始条件和三类边界条件及曲面边界的处理方法。使用PINNs求解了一个一维导热问题。求解结果与理论解的最大相对误差为0.001 7%，平均相对误差为0.001 1%。使用一个简化叶片的导热问题作为案例，将PINNs与传统有限元方法进行对比，探究了PINNs不同的网络架构和超参数对结果的影响。对于简化叶片的导热问题，有限元方法求解时间为11.7 s,PINNs平均求解时间为8.96 s，求解结果的最大误差为1.03%，平均误差为0.139%。微调实心叶片的内冷源强度，在训练收敛的PINNs基础上重新采样计算，新的计算收敛时间为1.41 s，证明了PINNs方法具有设计条件微调时的快速计算能力。     

基于有限元法的平纹编织C/SiC复合材料等效导热系数预测方法

提出了一种采用有限元法计算平纹编织C/SiC复合材料等效导热系数的方法.首先研究了材料的细观结构,根据材料显微照片建立了带基体碳纤维束复合材料的单胞模型,采用基体中孔隙分布随机生成的单胞模型,计算了孔隙率对基体等效导热系数的影响,通过施加3组边界条件计算出带基体碳纤维束和复合材料的等效导热系数.最后应用提出的方法计算分析了碳纤维体积分数和孔隙率对复合材料等效导热系数的影响规律.结果显示:复合材料等效导热系数随碳纤维体积分数增大而线性下降,碳纤维体积分数从54%增加到78%的过程中,复合材料y轴方向的等效导热系数下降了12.8%,x与z轴方向的等效导热系数同时下降了8.6%;复合材料等效导热系数随孔隙率增大呈加速下降趋势,孔隙率从0增加到30%的过程中,材料的x与z轴方向的等效导热系数下降了22.91%,y轴方向的等效导热系数下降了34.66%.     

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

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

预测平纹编织C/SiC复合材料等效导热系数的三维纤维随机模型

提出了一种能够刻画平纹编织C/SiC复合材料编织结构的三维纤维随机模型,探讨了微观结构影响宏观等效导热系数的机理。采用蒙特卡罗法模拟纤维位置的随机分布,用有限元法计算等效导热系数,通过数理统计方法分析在不同纤维体积分数下等效导热系数的随机分布规律;并且比较了相同纤维体积分数下,纤维随机模型平均等效导热系数与纤维周期排列模型等效导热系数的差别。结果表明:在不同体积分数下,纤维位置随机等效导热系数都呈高斯分布;纤维随机模型中纤维会出现聚集现象,形成局部“导热热障”,所以纤维随机模型平均等效导热系数小于纤维周期排列模型,这种差异在中等纤维体积分数下最明显。      

界面层及其建模方式对编织结构CMC材料导热系数影响研究

针对编织结构陶瓷基复合材料（Ceramic matrix composite,CMC）中基体、纤维束和界面层等组分的不同传热特征,以及考虑到界面层结构极薄的尺寸特征,探究了界面层及其在细观结构代表单元中的引入方式对编织结构CMC材料内部热量传输特征和各向异性导热系数的影响。研究中对比分析了不考虑界面层、含隐式界面层和含显式界面层等三种代表性体积单元模型的温度场、热流密度场及各向异性等效导热系数,获取了界面层导热系数对CMC材料整体导热系数的影响规律。研究结果表明：编织结构CMC材料内部温度场存在明显的不均匀性,不同模型计算获取的热流密度场具有明显区别。同时基于三个模型获取的各向异性导热系数也具有较大差异,显示界面层方法预估精度较高。此外,随着界面层导热系数增加,CMC材料整体各向异性导热系数明显增加,其对水平经纱Y方向上等效导热系数的影响最大。     

求半透明材料导热系数的复合换热反问题

在激光脉冲法测不透明材料导热系数的基础上,本文提出了一种已知半透明材料辐射特性,求其导热系数的复合换热反问题方法。      

基于物面弱边界条件的并行全隐式流场求解方法

提出了一种基于物面弱边界条件的Navier-Stokes方程并行全隐式求解方法。首先,实现了非结构网格格点格式有限体积法的粘性物面弱边界条件。以此为基础,推导了相应的粘性Jacobian矩阵。相比传统的物面强边界条件,弱边界条件的应用简化了该推导过程。接着,为了解决Spalart-Allmaras(S-A)模型的收敛和保正性问题,提出了基于物面弱边界条件的S-A模型无条件保正性收敛方法,通过构造特别的隐式矩阵同时保证S-A模型的收敛和湍流工作变量的正值。流场方程全隐式时间推进和无条件保正性收敛方法所得到的线性方程组采用多色高斯-塞德尔迭代法进行求解,并由MPI实现并行化。通过若干二维和三维算例对本文所提出的方法进行了测试,其结果表明:1)随着罚系数的增大,基于弱边界条件的结果趋向于基于强边界条件的结果;2)对于力的计算,弱边界条件自动满足通量一致关系而无需特别处理。     

基于流-固耦合传热的热气防冰系统干空气飞行蒙皮温度场计算研究

以机翼热气防冰系统为研究对象,建立了包含热气防冰系统防冰腔内外流场对流换热和固体结构导热的三维稳态流-固耦合传热物理模型,对整个计算区域生成混合网格,边界条件为第三类边界条件,采用计算流体力学方法以 FLUENT 软件为工具,对干空气飞行状态下流-固耦合传热模型进行了求解,获得防冰腔蒙皮内外表面对流换热系数分布和温度场结果,并对计算结果进行了分析。结果表明：防冰腔铝合金蒙皮沿展向和厚度方向导热显著,温度分布较均匀,防冰引气温度为200℃时,防冰腔蒙皮内外表面上最高温度为101℃,最低温度为21℃,3 mm厚的蒙皮同一点处内外表面最大温差仅为4℃,防冰腔排气口处气体的平均温度为63℃。热气防冰系统蒙皮温度场计算方法和计算结果,能够为热气防冰系统干空气飞行试验设计和测试中温度传感器的选型与布置提供依据。     

多孔陶瓷等效导热系数的近似计算

提出了一个近似计算多孔陶瓷等效导效系数的公式。该公式综合考虑了孔隙大小，孔隙度以及气相辐射与导热的影响，可用于在任意温度下具有任意孔隙度的多孔陶瓷的等效导热系数的计算，用文献数据检验表明，该公式在高温区有较高的计算精度。     

基于高温瞬态热响应的石英窗口导热系数反演

提出一种利用高温加热下瞬态热响应信息来同时反演石英窗口在不同温度下导热系数的方法。结合实验测试过程,首先建立石英窗口-石墨板-石英窗口三明治结构高温非灰体辐射-导热耦合传热模型。再通过实际测试温度响应结果和材料光谱辐射特性参数,构建石英窗口导热系数的遗传算法反演辨识模型。采用1 100 K加热条件下实验数据,反演辨识获得了常温至1 100 K石英窗口导热系数介于1.35~2.58 $\mathrm{W}/(\mathrm{m}?\mathrm{K})$,并拟合出高温导热系数关联式。结果表明:该材料辐射对导热系数影响随温度升高而增大,在1 100 K下辐射影响占比16.12%。最终通过变工况的高温传热过程测试,验证了方法的可靠性及导热系数数据的准确性。      

柔性导热纳米复合材料研究进展

柔性导热纳米复合材料在微电子、航空航天、智能驱动等领域具有广泛的应用。综述了以零维纳米粒子、一维纳米纤维和二维纳米片层为填料的柔性导热纳米复合材料的最新研究进展。概述了柔性导热纳米复合材料的研究中存在的主要问题及解决方法。纳米填料的分散性差是阻碍柔性导热纳米复合材料发展的重要问题,可采用物理包覆法、偶联剂法、接枝改性等方法对纳米填料进行表面修饰,提高纳米填料的分散性及与基体的结合性,从而提高材料的导热性能;另外完善导热机制的理论研究,建立导热微观模型,可为柔性导热纳米复合材料的实际应用提供重要依据。开发新型纳米填料,特别是混合填料体系是改善柔性导热纳米复合材料综合性能的重要途径。     

导热-辐射耦合传热的多尺度分析和数值模型

复合材料高温传热特性的准确预测对飞行器热防护结构的设计有重要意义,相关模型也是国家数值风洞工程中多相多介质计算模型的重要组成部分。针对周期性结构复合材料的高温传热问题,利用多尺度渐进分析方法,对包含导热方程和辐射传输方程的导热-辐射耦合传热模型开展了研究。建立了表征单元尺度模型及宏观平均导热-辐射耦合传热模型,获得了材料宏观等效导热系数与表征单元模型间的关系,发现宏观等效辐射吸收和散射等系数可通过表征单元内的体积平均求取。结合有限容积方法与格子Boltzmann方法,建立了复合材料导热-辐射耦合传热多尺度数值模型。采用二维常物性材料传热过程的模拟验证了多尺度模型的有效性,结果表明所建立的多尺度模型能够对温度场给出准确高效的计算结果。该方法有助于为复合材料传热特性的预测提供数值手段。     

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

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

基于Gama-Theta转捩模型的气冷涡轮气热耦合研究

基于非结构化混合网格,采用有限体积法求解了N-S方程,对流项格式为AUSM+,梯度的求解方法为加权最小二乘法,求解器采用LUSGS隐式求解方法,并进行了预处理和网格重排序。与导热方程耦合,搭建了气热耦合平台,采用直接耦合方法,交界面互相传递温度,并且采用能保证通量守恒的面积加权类的插值方式实现数据传递,将γ-Reθt两方程转捩模型应用在流场程序中。通过MARKII叶片4311和5411工况的实验结果进行了验证,结果表明采用转捩模型对压力分布的影响不大,与实验值吻合较好,通过间歇因子的分布可以看出,γ-Reθt模型成功地预测了从层流到湍流的转捩过程,在层流区采用γ-Reθt模型计算的涡粘系数与真实流动情况更加吻合,由于涡粘系数对传热的影响很大,从而有效地提高了边界层的层流区和转捩区的传热计算的精度,得到的温度和换热系数与实验值吻合更好,但是该模型是和SST模型耦合在一起的,由于SST模型的局限性,所以对激波和边界层干扰区域的模拟产生了7%左右的温度误差。