基于共轭梯度法的三维弯管内壁温度的反演研究.pdfVIP

中国工程热物理年会 传热传质学 学术会议论文 编号：093080 基于共轭梯度法的三维弯管 内壁温度的反演 1 1 2 刘波 ，卢涛 ，姜培学 （1.北京化工大学机电工程学院，北京 100029） （2.清华大学热科学与技术教育部重点实验室，北京 100084） （TelEmail：likesurge@） 摘要：用外壁测量温度与计算温度构建反问题目标函数，以有限单元法作为导热正问题求解方法，基 于共轭梯度算法建立了三维弯管非稳态的导热反问题模型，利用弯管外壁的测量温度反演出内壁的温 度，并与内壁温度实验值对比。对比结果显示，反演值与实验值符合较好，利用该模型可以精确地反 演出弯管内壁面上各点的温度波动值。 关键词：共轭梯度法，弯管导热，内壁温度反演 0 前 言 导热反问题在核能工程、航空航天、冶金及化工等领域有广泛的实际应用背景和价 值。导热反问题（Inverse Heat conduction Problem-IHCP ）是指利用研究对象边界的测量 温度，通过一定的反演算法，确定其它边界上的温度、热流密度、边界形状或源项强度、 热物性等未知参量的一类问题。由于导热反问题的不适定性及非线性，使得导热反问题 的求解远比导热正问题复杂。 国内外学者利用共轭梯度法、最大熵法、正则化算法及神经网络法等不同的反演算 法，对不同的导热问题进行了分析。Chen 等[1] [2] 和Lin 等 利用共轭梯度法分别求解了肋 片和平行板通道的导热反问题，并探讨了测量误差对反演结果的影响。Huang 和Wang[3] ， [4] [5] Huang 和Chen 等人应用共轭梯度法求解了三维的热传导反问题。薛奇天等 ，Xue 和 [6] [7] Yan ，杨海天和胡国俊 等应用共轭梯度法求解了多宗量的热传导反问题，得到了较满 意的结果。范春利等[8]应用共轭梯度法求解热传导反问题，识别试件内壁的缺陷。针对 在求解存在急剧骤变热流密度或温度的导热问题时，共轭梯度法等传统算法很难获得较 好的反演结果，Kim 和Lee[9] [10] ，Kim 等 在共轭梯度法的基础上应用最大熵法进行求解， 得到了较为精确的反演值。由于共轭梯度法在迭代过程中具有一定的抗不适定性[11] ，目 前仍然是导热反问题求解的最常用方法之一。 弯管是工业中用于流体输送最为普遍的装置，如化工装置、核动力装置的管道系统。 当管内壁由于热冷流体掺混、气液共存时，流体温度波动激烈，容易导致管壁产生热疲 劳，对管道损害非常大，严重时可能造成事故，因此人们尤其关心管内壁的温度波动情 况。由于高温或受测温条件限制，想要直接测量管内壁面温度波动数据比较困难，而管 外壁面的温度测量则较容易，因此可利用外壁面的测量温度，通过求解导热反问题 （IHCP ）来预测管内壁面的温度波动。本文使用有限单元法建立了适于进行敏感度分析 的导热正问题模型，应用共轭梯度法求解导热反问题，反演出弯管内壁面上各点的温度 随时间波动情况。 1 正问题 物理模型如图 1 所示，本文使用有限单元法求解导热问题，离散单元为六面体八节 点等参单元。控制方程为 2 2 2 ∂T ∂ T ∂ T ∂ T ρc λ( 2 + 2 + 2 ) (x ,y , z