三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

有限元传热分析中燃气辐射处理方法及其在飞机结构热分析中的应用

 研究了有限元传热分析中有吸收-发射性介质( 或称弱吸收性介质, 例如燃气) 参与辐射传热时的处理方法。通过对有、无吸收-发射性介质的封闭包壳的辐射能量交换方程的分析, 提出了用表面边界单元替代其辐射效应的准则; 研究了有吸收-发射性介质参与辐射时表面间辐射角系数计算方法及编制了相应的计算程序。最后, 对某型飞机机尾罩的温度场作了计算并对结果进行了分析。     

MTPS蜂窝夹芯结构传热性能及热应力分析

对于金属防热结构的蜂窝夹芯结构的稳态情况,基于材料的全灰体假设,同时考虑热传导和热辐射两种传热形式对温度场的耦合作用,利用热流量守恒建立了蜂窝芯层温度场的非线性积分方程。离散化后利用数值方法得到方程组的数值解。对于美国兰利研究中心的实验结果,与本文方法的对比计算结果基本吻合。进一步,利用计算结果讨论了给定面板温度边界情况下,下面板、柱体层的灰度、蜂窝结构长径比对夹芯温度场的影响。并根据温度场和近似的应力分析模型,用半解析结果讨论了稳态情况下蜂窝芯层上的热应力。     

激光焊接温度场数值模拟

建立焊接的温度场控制方程、小孔传热模型、求解方法、初始条件及边界条件的数学模型。在此基础上利用有限元系统软件ANSYS/Thermal模块,对激光焊的数学模型进行计算,并且获得了焊接过程的温度变化曲线及焊件整体瞬态温度场分布。     

多孔介质传热模型在多孔壁湍流中的适用性

为了考查不同的多孔介质传热模型在不同工况下的适用性,对带有高孔隙率多孔介质壁面槽道湍流及其传热进行了直接数值模拟研究。在多孔介质层外流体区域,通过有限差分方法求解不可压缩Navier-Stokes方程和温度对流扩散方程;在多孔介质层内,使用修正的Darcy-Brinkman-Forchheimer模型描述高孔隙率多孔介质阻力,以及分别采用局部热平衡（local thermal equilibrium, LTE）模型、局部非热平衡（local thermal non-equilibrium, LTNE）模型、理想金属（ideal metal foams, IMF）模型计算温度分布。通过对所得热场的统计特性的分析比较,探究了不同Biot数下水和空气两类流体介质的多孔介质传热模型的有效性。研究表明：LTE模型不足以准确预测金属泡沫多孔介质内传热问题,其等效导热系数因仅考虑孔隙率因素而低估了多孔介质层的传热能力;IMF模型在小比热容流体介质的算例中表现良好,可以代替LTNE描述多孔介质层内的传热,而在大比热容流体介质的算例中表现不佳,需要考虑比热容以及流固两相间的传热能力对预估的固体相温度分布...     

高温圆管内多孔介质耦合换热的辐射作用

通过数值模拟圆管内填充多孔介质的辐射对流耦合换热,研究高温固体骨架辐射效应对温度分布及换热的影响.基于局部非热平衡模型分别建立流体和固相能量方程,采用蒙特卡罗法求解固体骨架的辐射换热;对不考虑热辐射引起的温度场偏差和管壁发射率以及多孔结构参数的影响进行讨论.结果表明:多孔固体骨架的辐射效应对入口段温度场的影响明显,不考虑辐射将导致较大偏差,壁温为1500K时最大偏差为16%.管壁发射率对温度场影响较小,壁温为1500K时影响小于3%.孔隙率或孔径增大,壁面辐射热流密度比例增加,辐射效应体现明显.      

多层高温隔热结构的传热特性

建立了高温多层隔热结构传热计算模型,采用蒙特卡罗方法模拟每个半透明隔热材料层内的一维辐射传递,采用有限体积法对多层隔热结构辐射导热耦合换热能量方程进行求解,对多层隔热结构的瞬态传热及热响应进行了数值模拟分析,研究了反射屏个数等因素的影响.结果表明,反射屏对多层结构隔热性能的影响取决于隔热材料的辐射特性与导热性能,当隔热材料衰减系数小时,将屏布置于高温区可提高隔热性能,当衰减系数大时,反射屏布置在低温区进行蓄热则更加有利.     

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

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

轴对称喷管内外流场与结构温度场耦合计算

利用有限体积的对流项二次迎风插值格式和重整化群(RNG)k-ε湍流模型,二层增强型壁面函数,同时利用球形谐波法考虑热辐射的影响,以灰气体加权模型(WSGGM)确定气体介质的辐射性质,求解N-S方程、热传导方程、考虑吸收-发射性气体介质的辐射传输方程。采用流固耦合的流动与换热模型,流场与结构温度场互为边界条件交换数据,实现了流场解算与温度场解算的耦合数值分析。采用此计算模型对静止的轴对称收-扩喷管进行了数值模拟,计算得到的喷管壁温与试验数据吻合良好。在此基础上,进一步开展了巡航状态下轴对称收-扩喷管内外流场与结构温度场的耦合数值研究。     

采用混合场方法计算二维介质-导体复合散射体雷达散射截面

在二维介质导体复合散射体计算中,对导体部分应用积分方程时,不论电场积分方程还是磁场积分方程,去求解导体的散射问题时,在内谐振条件下将出现解答的不唯一性或不稳定性,产生谐振伪解,给计算工作带来困难.为避免谐振伪解,对导体部分采用混合场积分方程方法结合介质部分体积分方程和矩量法,对TM入射平面波情形建立了该模型求解方程,并编程求解出了该模型雷达散射截面(RCS).     

基于STL文件和有限元法的在轨航天器关键部件热计算

针对STL(stereo lithography)文件在传递复杂几何实体模型信息方面具有精度高的特点,提出了一种基于STL文件计算复杂结构角系数和外热流的方法.根据单元和节点的拓扑关系识别六面体网格边界单元并进行辐射换热计算.详细阐述了基于STL文件和有限元法进行复杂结构的角系数和外热流的计算方法.研究了导热-辐射耦合计算的有限元方法.最后利用有限元法计算航天器关键部件在轨的三维瞬态温度场.      

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

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

高温下加热/冷却管内参与性介质的辐射与对流耦合换热研究

本文通过数值模拟方法研究了高温下加热/冷却管内层流充分发展流动的参与性介质热起始段的辐射与对流耦合换热。采用离散坐标法求解圆柱坐标系下的非线性积分微分辐射传递方程,采用有限差分法求解辐射与对流的耦合换热能量方程。通过计算,分析了不同壁温及入口温度下,介质导热系数、壁面发射率、介质散射反照率及散射相函数等因素对局部努谢尔数、壁面热流密度及有效加热/冷却长度的影响,并同时采用努谢尔数和壁面热流密度评价了耦合换热中辐射对对流换热的影响。      

复合喷管热结构耦合计算的一种策略

基于流固耦合完成复合喷管流动换热、传热过程温度场的数值模拟,编写了计算流体动力学(CFD)与有限单元法(FEM)的接口程序,快速实现温度载荷和机械载荷数据的传递,提出了复合喷管热结构耦合计算一种策略.利用直接约束法考虑了复合喷管部件界面问高温碳化带来的边界非线性因素.计算结果表明耦合方法得到的复合喷管温度场是合理的,高...     

一维和二维电热除冰相变传热特性的参数影响分析

基于飞机电热除冰过程冰层融化的特征,开展了一维和二维电热除冰相变传热特性数值计算研究和参数影响分析,重点考虑了加热模式、冷却时间、加热功率和加热单元间隔等参数对冰层相变传热的影响。采用了基于焓-多孔介质方法的热焓模型,将计算区域看作是包括多层材料和冰、水及其混合区的多孔介质,采用了结构化网格拓扑结构对计算区域进行划分,采用了有限体积方法对控制方程组进行离散,采用线性插值的方法获得混合区的物性参数,耦合能量方程和液态水体积分数公式,迭代求解了计算域的温度分布,获得了不同材料间界面温度的变化,重点分析了冰-保护层界面的温度变化。基于二维电热除冰模型的冰-保护层界面温度不均匀特征,提出了耦合考虑冰-保护层界面热点和冷点温度的冰脱落温度判断准则。研究表明:高功率的周期性加热模式要优于低功率的连续性加热模式,采用合理的冷却时间和加热功率,可获得更低的能量消耗和更好的除冰效果。合理设置加热单元间隔可以提高周期性的除冰效率,但也会形成冷点和热点,造成冰-保护层界面温度分布的不均匀。冰-保护层界面的冷点类似于锚点,即使此刻热点的冰已经融化,整个冰层也无法脱落和剥离。因此,冰脱落的判断要耦合考虑界面热点和冷点的温度特征。      

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

Parallel algorithm and its convergence of spatial domain decomposition of discrete ordinates method for solving radiation heat transfer problem

To improve the computational efficiency and hold calculation accuracy at the same time,we study the parallel computation for radiation heat transfer. In this paper, the discrete ordinates method(DOM) and the spatial domain decomposition parallelization(DDP) are combined by message passing interface(MPI) language. The DDP–DOM computation of the radiation heat transfer within the rectangular furnace is described. When the result of DDP–DOM along one-dimensional direction is compared with that along multi-dimensional directions, it is found that the result of the latter one has higher precision without considering the medium scattering. Meanwhile, an in-depth study of the convergence of DDP–DOM for radiation heat transfer is made. Analyzing the cause of the weak convergence, we relate the total number of iteration steps when the convergence is obtained to the number of sub-domains. When we decompose the spatial domain along one-,two- and three-dimensional directions, different linear relationships between the number of total iteration steps and the number of sub-domains will be possessed separately, then several equations are developed to show the relationships. Using the equations, some phenomena in DDP–DOM can be made clear easily. At the same time, the correctness of the equations is verified.     

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

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