采用共轭梯度法的圆管内壁温度导热反问题的求解研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：083361 采用共轭梯度法的圆管内壁温度 导热反问题的求解 张有为 李辉 姜培学 热科学与动力工程教育部重点实验室，清华大学热能工程系，北京，100084； 联系电话：010E-mail：jiangpx@tsinghua.edu.cn 摘要：建立了采用共轭梯度法的非稳态二维导热反问题数学模型，并利用FORTRAN 语言编写了 程序进行计算求解，由同心圆管外壁温度随时间的变化得到管内壁温度随时间变化，。计算结果表明， 计算值与实验值相符较好，共轭梯度法能够准确地计算得到管内壁不同位置温度随时间的变化。 关键词：导热反问题、共轭梯度法、圆管内壁温度波动 0、引言 在工业生产、科学研究领域的许多场合都需要知道管内流动的管壁温度。对于管外 壁温度，有许多种实验测量方法；而对管内壁温度的实验测量非常困难。利用导热反问 题方法，可通过数值求解由管外壁温度得到管内壁温度。导热反问题（IHCP:Inverse Heat Conduction Problem ）是相对于导热正问题而言的。对于导热正问题，物体的热物性参 数和边界条件为已知，可以由此得出物体的温度分布。而当物体的热物性参数或部分边 界条件必须通过已知参数来确定时，就形成了导热反问题。 由于导热反问题的非适定性、计算量大等特点，使得导热反问题的求解比较困难。 国内外学者在导热反问题已进行了大量的研究，其中，针对一维问题、稳态问题的相关 研究工作较多，而对多维问题、非稳态问题的研究工作较少[1-6] 。本文所研究的导热反 问题是针对二维同心圆的管内流动，在已有研究工作的基础上，建立了采用共轭梯度法 的导热反问题数学模型，由圆管外壁不同位置温度的变化求解管内壁不同位置温度随时 间的变化。 1、数学物理模型 本文所研究的对象为二维同心圆管，圆管壁的二维非稳态导热方程为： 2 ∂T 1 ∂ ∂T 1 ∂ T ⎛ ⎞ α( ⎜r ⎟=+ 2 2 ) （1） ∂τ r ∂r ⎝ ∂r ⎠ r ∂φ 圆管外侧为第三类边界条件： ∂T λ h T T （2 ） − ( )w out ( =− sir ) ∂n 圆管内侧为第一类边界条件： * T T （3 ） 初始条件为： T x y ( T, ,0) 0 （4 ） 式中T * sir