功能梯度材料瞬态热传导数值模拟.pdfVIP

中国工程热物理学会 传质传热学 123053 学术会议论文 编号： 功能梯度材料瞬态热传导数值模拟 龚京风，明平剑，宣领宽，张文平 哈尔滨工程大学动力与能源工程学院，哈尔滨 150001 Tel: E-mail: gongjingfeng@126.com 摘 要 本文介绍一种格点型交错非结构有限体积法(FVM)，用于求解三维不规则区域功能梯度材料(FGM) 瞬态热传导问题。该方法基于FVM 离散瞬态热传导方程，温度存储在网格节点上，材料属性及温度梯度 存储在网格单元中心上。数值模拟FGM 立方体、FGM 转子内的瞬态热传导，并分别与解析解、数值解进 行比较，验证本文方法的正确性。数值结果表明本文方法可以有效的计算不规则FGM 的瞬态热传导问题。 关键词 功能梯度材料；瞬态热传导；有限体积法；交错网格 0 引言 功能梯度材料(FGM)是为满足某些工程领域对材料性能的特殊需求而发展起来的一种 新型材料。其概念由日本学者新野正之、平井敏雄1984 年提出。FGM 的设计思想是，通过 控制材料性质在空间的连续变化，使结构元件不同位置具有相应的性能，从而能够在苛刻的 环境下工作。 由不同材料构成的FGM 可作为新型耐热材料，用于存在高温和大温度梯度的环境中， 如航空航天器的机头尖端、发动机燃烧室内壁。研究FGM 中的三维瞬态热传导问题具有重 要意义。近年来，学者们提出并发展了不同的数值方法用于研究热传导问题，如无网格方法 (MLPG)[1-3]、边界元方法(BEM)[4] 、有限差分方法(FDM)[5] 、有限体积方法(FVM)[6] 。其中FVM 能像FEM 一样适用于非结构网格和复杂边界，且处理效率与FDM 相似，远远高于FEM 。 特别是控制容积有限元法(CVFEM)，它兼有FVM 能保证满足守恒性的特点及有限元法对不 规则区域适应性强的优势[7] 。CVFEM 由Baliga 和Patankar[8] 1980 年提出，用于研究对流扩 散问题。随后得到迅速发展，并被用于不同领域[6,9-10]，CVFEM 是数值研究的有效手段。 本文基于CVFEM 发展一种显格式的格点型交错非结构FVM ，用于求解三维FGM 瞬态 热传导问题。热传导方程采用显格式推进求解。前处理和后处理是基于内部程序代码 GTEA(General Transport Equation Analyzer ，即通用输运方程求解器) 中的CGNS 系统接口[11] 实现。该方法适用于材料属性随空间任意变化的FGM 热传导问题。通过计算FGM 立方体、 FGM 转子的瞬态热传导，验证本文方法的正确性。 1 数值方法 1.1 控制方程 本文考虑各向同性FGM 中瞬态热传导问题，FGM 的材料属性是空间变化的，其微分 控制方程为 ∂T (x, y, z, t) ∇⋅ k x, y, z ∇T x, y, z,t +Q x, y, z,t ρ x, y, z c x, y, z ( ( ) ( )) ( ) ( ) ( ) ∂t \* MERGEFORMAT (1) 式中T(x , y , z, t)为温度，k(x , y , z)为导热系数，Q(x , y , z , t)为外部热源( 以下公式推导中忽略) ， ρ(x , y , z)为密度，c(x , y , z)为比热。初始温度场可表示为 T T \* MERGEFORMAT (2)