同时反演材料热传导系数和比热的算法.pdfVIP

第 28卷 第 5期 计 算 物 理 Vo1．28．No．5 2011年 9月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Sep．，2011 文章编号 ：1001—246X(2011)05-0719-06 同时反演材料热传导系数和比热的算法 周 宇， 钱炜祺， 何开锋， 桂业伟 (空气动力学国家重点实验室，中国空气动力研究与发展 中心，四川 绵阳 621000) 摘 要 ：讨论 同时反演随温度变化 的热传导系数和 比热的算法．首先将材料 的热传导系数和比热表示成随位置和 时间变化的函数 ，利用伴随方程法获得 目标函数对 k(，)，c (，f) 的梯度 ；然后将材料的热传导系数和 比 热直接按温度 区间分段离散 ，建立 目标函数对k ()，c ()梯度与 目标函数对 k(，t)，C (，t) 的梯度的 关系；随后利用这种关系进行反演计算．算例表明，这样建立的将热传导系数和比热表示为温度的分段函数进 行反演的方法是可靠有效的，并且具有 良好 的抗噪性能． 关键词：热传导系数 ；比热；热传导逆问题；反演 ；伴随方程法 中图分类号 ：0482．2 文献标识码 ：A 0 引言 碳碳复合材料是再人飞行器防热系统使用的主要材料，在飞行器再入的高温环境下，这类材料的热传导 系数和比热会随温度变化 ，直接影响飞行器防热系统和飞行器舱 内的温度分布．因此，准确确定这类材料 的热传导系数和比热随温度变化的规律对于再入飞行器防热系统的设计具有重要意义．热传导逆问题研究 在航天领域 ，在核物理、冶金等工业研究领域 中有着广泛的应用背景 ，取得 了不少研究成果 ，如对表面热 流反演 、对 内部热源的识别 、黑体辐射反问题的研究 等．热传导系数和 比热 的确定实际上是一类典型 的热传导逆问题，现有的大多数研究都是固定热传导系数反演 比热，或是固定 比热后反演热传导系数，而在 工程实际中通常这二者都是未知的，因此 ，开展同时反演材料热传导系数和比热的算法研究更具工程实用价 值．文 [6—8]开展了反演材料热传导系数和比热的研究工作 ，但只考虑了材料热传导系数和比热为常数或 是温度 的已知函数形式的情况．对于假设材料参数为温度 已知函数的情况，首先需要对材料性质有足够的了 解才能给出较为准确的函数形式 ，以多项式函数为例 ，若选取 函数阶数过高容易引入虚假的震荡 ，而阶数过 低又难 以描述材料参数随温度 的变化 ；并且不同的函数形式要求重新推导灵敏度方程 ，缺乏通用性．文 [9] 采用的方法是将材料热传导系数和比热表示为随位置和时间变化的函数 ，再通过材料 内部温度场 的测量结 果用共轭梯度法来同时进行反演，最后整理出热物性参数与温度之间的关系．但是，这样得出的解有可能不 唯一 ，文 [9]对这一点没有进行深入论述．本文对这一问题进行深入讨论，给出如何同时反演热传导系数和 比热随温度变化函数的算法和算例． 1 反演材料热传导系数和比热的数学模型 要确定材料的热物性参数 ，需要开展相应的实验 ，其示意 图如图 1，受测材料试样 的左右两端受常值均 匀热流q。和q，其余两端近似为绝热边界条件 ，则此问题可以简化为一维传热问题．温度传感器置于测点 (m=1，…，Ⅳ)处 ，测量这些点的温度历程 ( ，)．将热传导系数和 比热表示为随位置和时间变化的函 数 ，则这一传热问题可简化为如下无量纲形式 ： c(x,tOT= ((x,tOT)， (1) 边界条件为 收稿 日期 ：2010—09—13；修回13期 ：2010—12—23 基金项 目：国家 自然科学基金 11072258)及空气动力学国家重点实验室基金(SKLA2009A0103)资助项 目 作者简介 ：周宇(1984一)，男，安徽怀远，实习研究员，主要从事计算流体力学，数理逆 问题研究，E-mail：yuzhou84@gmail．corn 720 计 算 物 理