一种基于反问题算法的管道壁面热流密度测量技术

　　1　引　言

　　热流密度是热工计量及工业控制中的一个重要参数。由于热量是过程量,传统测量方法对过程量的计量本身存在较大的难度,精确地测量管道壁面的热流密度是一个比较困难的问题。在实验室和工业装置中,由于流体在管道壁面之间流体流动和传热的复杂性,热流密度的实测数据都很少。热流密度的传统测量方法是应用热传导的傅里叶定律q=-λ(ΔT/Δx),这种方法的缺陷是在实际测量时很难做到流过测量探头的热流与流过被测物体其它部分的热流相同,有时测量误差非常大,同时干扰因素多、干扰作用大,需要进行多项的修正和补偿,测量的准确度较低,因此传统测量方法无法满足许多工业测量与控制中对热流密度准确计量的要求。针对高温超导实验中管道壁面热流密度的精确测量问题,基于对流换热反问题的求解[1],本文提出了一种简明、通用的测量管道壁面稳态热流的新方法,并从数值模拟实验上分析了测量误差对热流密度估算的影响。

　　2　测量原理

　　考虑常物性流体(如机油)流过平行平板管道的层流受迫对流情况,流体速度边界层已充分发展,而热边界层正在发展中。其中管道的一个壁面受到随空间变化的热流,而另一壁面为绝热。流体进入管道的初始温度为T0。平行平板管道的几何形状如图1所示,L为平行平板管道间距,b为所测管道某段长度。

　　为了测量热流密度,在管道内壁附近y1处沿流线方向放置M个测温热电偶,从这M个温度测量值来计算未知的边界热流q(x)。其数学模型为

　　这是一个已知测量点的温度值求解未知的边界热流q(x)的二维对流换热反问题。反问题的求解在数学上可以转化为以下泛函变分的最优化控制问题,即

　　这里Tm(x,y1)是在测量点处计算出的温度值,Tc(x,y1)是在同一个测量点上实际测量得到的温度值。由于反问题对测量误差的敏感性,在对式(3)求最小化过程中使用了目前流行的共轭梯度法[2, 3],这种方法的原理是基于变分原理,把反问题转变成三个问题:正问题、灵敏度问题及伴随问题。关键过程在于确定搜索方向和步长,其迭代方法为

　　这里T为灵敏度函数。共轭梯度法的迭代计算过程为

　　(1)选择一个猜测的初值q,求解正问题并计算Tc(x,y1)。

　　(2)根据Tc(x,y1)和Tm(x,y1),判断是否满足收敛停止标准J(qk+1)<ε2,如满足则停止迭代,否则按下一步计算。

　　(3)计算共轭系数γk和下降方向Pk。

　　(4)计算灵敏度函数得到温度增量ΔT(Pk),从而计算步长βk。

　　(5)令k=k+1,计算新的壁面热流qk+1(x),计算新的Tc(x,y1),回到第2步。