竖壁液膜流壁面热流率对流动影响的实验的研究.pdfVIP

中国工程热物理学会 传热传质学术会议 竖壁液膜流壁面热流率对流动影响的实验研究+ 王补宣张金涛彭晓峰 清华大学热能工程系，北京100084 ‘ 摘要：本文用相空间重构的方法，对竖壁薄液膜流动的液膜厚度时间序列进行了重构， 并甩分维数对奇怪吸引予的动力特征进行了描述，由此描述了壁暖热流率对竖壁薄液 膜流动特性的影响。 1．引言： 竖壁液膜流动的动力特性受多种因素的影响，已有大量的文献[1-6]从实验和数值计 算等不同方面探讨了质量流率(雷诺数)对流动特征的影响。理论及实验研究的结果表 明，竖壁液膜流动是本征地不稳定的，在较小的质量流率时为规则的波状流。当质量 流率大于某个临界值时，流动演变成湍流状。我们通过理论分析邛】表明薄液膜内熟状 态不平衡也是导致液膜流动状态演化的重要因素。在此基础上，通过实验并结合分维 数描述方法，本文定量地描述了壁面热流率的不同对薄液膜流动状态的影响。 2．实验设计 管，其表面光洁度与外径误差商上。液膜由一个0．5mm间隙的环形分配器形成。实验 段上方有高位水箱，由溢流管保持液位恒定。 实验采用HP34970A数据采集／开关单元及34902A16通道多路转换器。该系统具 有6。的采样精度，采样速率达600读数／秒。在液膜厚度的测量中，采样间隔是15ms。 3．实验结果及流态的描述 3．1实验结果 本文报道了两种质量流率及不同热流率下液膜厚度的测量结果。图2(a)和(b)分 ． 度振幅与液膜厚度平均振幅之比)随时间t的变化曲线。从图2(a)及(b)中可以看踏 液膜波动是混乱，而且存在弧立大波及很深的波谷。 3．2时问序列的相空间重构 Takens指出，离散的时间序列可以通过的时间延迟的办法重构系统的相空间， 国家自然科学基金资助课题 H—13 从而得到该系统的定量描述。采用时闻延迟构造相空问的方法是] (1) X。(t)={u(t)，u(th)，…，u(t竹(k‘1)) 其中u(t)是时闯序列，x。(t)是构造的相空间，t是延迟时间，k是系统的基层 维数(beddilnensiOn)。如果时间序列不含噪声，则其是时问延迟可以是任意值。对 于含有噪声的序列，应该通过选取延迟时间，使得噪声的影响降至最小。 本文采用文献[8]的方法选择延迟时间，选取的大致步骤是，假设一个延迟时问， 构造相空间，从构造的相空间中任取两列，计算此两个序列信息熵。使信息熵达到第 一个晟小值的延迟时间就是相空间重构所需要的。例3(a)和(b)是质量流率O．017kg／s 和0．033kg／s的冷态液膜的重构相空间的相圈。 3．3系统的分维数 相空间重构所得的系统的特征可用系统的维数定量地表征出来∞。“。如果该系统 是混沌的．则其维数不是整数。根据文献[12]，为了从时间序列中计算出分维数，在进 行相空间重构时，基层维数k必须满足条件： k_2m+l (2) 其中m为混沌吸引子的维数。 从相空间计算系统的分维数已有多种方法，本文采纳最常用的关联维数的计算方 法以及文献[13]所给出的分维数计算方法。关联维数的核心内容如文献[i0，11]所指 出的，一个由时间序列构成的相空阔的关联维数可表示为 驴脚警 ㈤ 其中： c(班嘞而‰‰¨对]IX，--X]怿) 在关联维数的计算中，N应该尽可能地丈。受计算机能力的限制，本文的实验数据处 理采用N=5000和k=36。用PII400PC计算一个：}Ic态的关联维数需用15小时。 目4(a)给出了质量流率的液膜流动，其关鞋维数与壁面加热率的关系。 由于关联维数的计算太耗时，故本文也采用了另一种分维数计算方法。该方法由 文献[13]提出，被文献[14]所采用，其核心足：