脉动热管的非线性混沌分析.docxVIP

脉动热管的非线性混沌分析 0 脉动热管模型 动脉冷热压缩性（pm）也被称为学术板运动管，由科学家akuchi制成。 脉动热管虽结构简单, 但内部气液两相流动和能量传递复杂, 通过模型的简化难以描述管内状态, 管内汽泡的生成, 运动, 破裂以及气液两相流动这些复杂因素影响, 导致流动行为异常复杂。宋延熙、徐进良等 牛蔚然 1 实验设备 图1为脉动热管模型和附属实验装置, 系统主要由闭式脉动热管测试件, 温度采集系统, 电加热和冷却系统组成。具体参数详情同文献 2 分析与讨论的结果 2.1 混合特性描述 结果表明, 热管传热性能主要取决于蒸发段温度, 冷凝段波动较小, 为避免热电偶的测量误差, 本文分析只取蒸发段温度测点时间序列。 2.1.1 大尺度变化与加热负荷之间的关系 实验数据初步分析发现, 随热负荷的增加, 波动呈现越来越强的趋势。如图2所示, 以丙酮和去离子水为工质, 当加热负荷较低时, 温度波动具有一定随机性 (如图2-a-a、2-b-b所示) , 看不出有周期的成分, 但随加热负荷的增加, 壁温时间序列出现类周期性波动 (如图2-a-e、2-b-e所示) , 这种波动并非严格的周期运动, 而是同时出现一定的不规则运动。波动与管内大汽泡生成运动有关。每个大尺度的波动之间存在着小尺度的波动, 这种高频小尺度波动与小汽泡运动有关。随负荷增加, 温度时间序列波动越剧烈, 在大剧烈波动同时又看似有周期成分。故而可以称脉动热管壁温波动是一种介于周期与随机之间的运动。 2.1.2 放空频率的测试 图3为三种工况热负荷维20 W的功率谱, 其中图3-c为热管管内不充装任何工质的空管, 所示为一条平滑的曲线, 没有尖峰, 表明在空管时壁温只有热负荷一种影响因素, 也验证不存在振荡频率。从如3-a和图3-b可以发现, 无论是去离子水还是丙酮都存在尖峰和连续谱, 不同工况存在不同频谱走势, 在频谱上观测管内流动存在不少随机成分, 但依然与随机信号的谱图不一样。且功率幅值在低频段以指数次幂急剧降低, 有明显混沌特征。 2.1.3 时间延迟对系统可测性的影响 图4所示为几种工况在55%充液率时, 不同热负荷、不同工质下的自相关函数, 自相关函数表示信号一个时刻与另一个时刻的依赖关系, 本文采用Matlab程序计算。图4-a可见, 对于同种工质下, 自相关函数随时间延迟的增加而降低, 表明管内具有混沌特性, 即随时间的增加, 系统越来越不可测。随热负荷的增加, 自相关系数随时间延迟的降低速度越来越大, 可预测性减少。图4-b在20 W热负荷下发现, 空管自相关系数呈现线性关系, 丙酮和去离子水都随时间延迟而下降, 丙酮的下降速度大于去离子水, 表明丙酮的可预测性更低这与后文分析以去离子水为工质的关联维数小结果一致。 2.1.4 壁温信号分析 图5为两种工况热端壁温的lg (R/S) -lgτ关系图, 曲线斜率即为Hurst指数。由图可知曲线的线性关系较好, 表明壁温信号有很好的自相似分析特征。图中两种工况都在τ=103.5左右曲线斜率开始突变, 在此之后均出现几个大的波浪缺口。根据Mandelbrot的理论, 当线性关系出现几个缺口就对应几个周期成分, 故而在功率谱上很难观测到尖峰。 2.2 壁温热负荷的系显关系 脉动热管温度序列存在混沌特性以及分形特性, 表明温度波动既有混沌又有分形的特征, 关联维数是定量刻画分形特征的参数, 是分形最基本的特征。对于混沌吸引子而言, 表征描述动力学系统最低的基本变量个数。因此, 分析关联惟数与传热的关系显的比较重要。计算方法采用G-P程序算法。 如图6-a所示, 以丙酮为工质 (Q=40 W) , 取热端温度计算得到的关联积分与尺度的关系曲线。计算发现当嵌入维m增加时曲线的线性区域的斜率也相应增大, 这是混沌时间序列区别其他时间序列的特征。结合吸引子可知, 在热负荷40 W时, 混沌形态已充分展开, 故不存在多尺度性, 只有在流态混杂的时候才会出现多尺度性, 计算得到关联维数在5.75左右, 说明在这一工况影响管内流动的因素达到比较稳定的状态, 关联维数反映壁温信号波动的影响因素。如图6-b所示, 给出了不同热负荷下的关联维数随嵌入维数增加的收敛图。无论在何种热负荷下关联维数都随嵌入维先增加后维持稳定的趋势。当dm≥14时, 关联维数的收敛于稳定的值即为关联维数达到饱和值。 2.3 多尺度分析对壁温时间序列的相关性 2.3.1 小波温度波动信号的分解 对壁温信号进行小波分解十尺度, 小波函数有多重选择方式, 本文对各种小波函数的误差分析如图7所示: 小波分析有多种小波函数, Daubechies (db) 小波由于有良好的紧支性, 正则性, 使得重构过程比较光滑, 在离散信号分析中得到广泛应用。但不同紧支