管道内壁侵蚀形状识别的无网格法研讨.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083349 管道内壁侵蚀形状识别的无网格法研究 1 2 谭建宇 ，刘林华 ，杨建国 （哈尔滨工业大学(威海) 汽车工程学院，山东 威海 264209 ） (Tel.: 0631-2650375，Email: tanjy@hit.edu.cn) 摘要：本文使用最小二乘配点无网格法结合共轭梯度法对管道内壁面侵蚀形状进行了反演识别。最小 二乘配点无网格法使用节点和辅助点离散求解区域，采用移动最小二乘近似构造试函数，进行导热正 问题求解；反演过程使用共轭梯度法计算目标函数，Akima 三次样条函数被用来将离散的未知内边界 点连成光滑曲线。经算例验证，使用最小二乘配点无网格法结合共轭梯度法反演具有较高精度。 关键词：无网格法 共轭梯度法 几何反演 导热 0 引言 火力发电厂、大型化工厂液体输送管道在使用过程中，由于长期受到高温、高压以 及化学侵蚀作用，管道内壁容易发生破损，影响管道的正常运行。因此，对该类管道的 内壁面侵蚀状况的进行及时、准确的检测，对火力发电厂、大型化工厂安全与经济运行 都有着重要意义。 由于输送管道表面为曲面，输送距离较长，位置变化不定，这给管道检测带来了一 定的难度，现在检测多采用红外热像仪测量管道外壁面温度，而后采用数值方法反演内 [1] [2] 壁面形状。国内外学者针对该类问题，使用边界元法 、有限差分法、有限单元法 等 数值方法结合最优化反演方法对其进行了研究。传统数值方法需要在反演过程中对变化 的求解区域进行反复的网格划分，而边界元法较为复杂，需要首先求出其基本解。 为避免以上问题，本文使用最小二乘配点无网格法[3]结合共轭梯度法对管道侵蚀形 状进行了反演识别。最小二乘配点无网格法使用节点和辅助点离散求解区域，采用移动 最小二乘近似构造试函数，并进行导热正问题的求解；反演过程使用共轭梯度法计算目 标函数，Akima三次样条函数[4]被用来将变化的待识别边界节点连成光滑曲线。文中使 用两个典型算例对最小二乘配点无网格法结合共轭梯度法反演进行了验证。 1 数值方法 1.1 能量方程 由于管道内液体流量、温度基本稳定，管道材料确定，可针对某一剖面进行研究， 因此该问题可视为二维、稳态、常物性、无内热源导热问题，其控制方程为： 2 2 ∂ T ∂ T 2 + 2 0 在Ω上 (1) ∂x ∂y 其边界条件为： 资助项目：国家自然科学基金(No T T , (x , y ) ∈Γ (2) w 1 −λ∂T ∂n qw , (x, y ) ∈Γ2 (3) −λ∂T ∂n h(T =−T ) , (x , y )