数值分析环形池热毛细对流的不稳定性.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093300 数值分析环形池热毛细对流的不稳定性1 石万元，李友荣，彭岚 (重庆大学动力工程学院，重庆400044) Tel ：023； Email ： shiwy@ 摘要: 通过线性稳定性分析，确定了环形浅池内硅熔体和Pr数介于6.7至57.9之间的几种硅油热毛细对 流失稳的临界Marangoni数Ma 、临界波数m 、临界相速ω ，获得了不同Pr数下的热流体波，特别是在 c c c 高Pr数下的热流体波。结果发现：当Pr≤15.9时，随着Pr数的增加，发生热流体波的临界Ma 数迅速上 升。Pr=23.9时是一个转折点，它所对应的Ma 小于Pr=15.9时的Ma 。当Pr≥23.9后，临界Ma数又随Pr c c 数的增而缓慢增加。临界波数m 和临界相速ω 随Pr 数的增加而降低。Pr≥23.9后的热流体波表现出与 c c 中、低Pr数流体不同的特征。 关键词: 热毛细对流, 线性稳定性分析, 热流体波 0. 引言 热毛细对流的一个显著特点是当界面温差超过一定的阈值后会失去稳定性，产生各种 各样的耗散结构，如热流体波、滚胞等。Smith 和 Davis[1]通过线性稳定性分析预测了 无限大浅液层热流体波的特点。Azami 等[2]在硅熔体(普朗特数 Pr=0.011)表面观察到这 样的热流体波。Schwabe 等[3-5]在地面和空间微重力环境测量了硅油(Pr=6.8)热毛细对流 失稳的临界条件，观察到热流体波和滚胞等失稳结构。Garnier 等[6,7]观察到环形池内硅 油(Pr=10)螺旋状的热流体波。Li 等[8,9]通过数值模拟得到硅熔体和硅油表面热流体波， 验证了 Schwabe 等人的实验结果。通过实验和数值模拟，我们[10,11]在深度为 1mm 的硅 油浅液层表面也获得热流体波。 [12-14] 线性稳定性分析是研究热毛细对流不稳定性特征的重要手段，Hoyas 等 通过线 性稳定性分析确定了发生热流体波的临界条件，计算得到了热流体波，但他假设流体 的Pr 数为无限大。最近，Gelfgat 等[15-18]也通过线性稳定性分析确定了热毛细对流失稳 的临界条件，但他没有计算热流体波。我们[19] 以前用线性稳定性分析预测了不同深度 硅油液层热毛细对流失稳的临界条件。 Smith 和Davis[1] 的线性稳定性分析认为热流体波的特征与流体的Pr 数有关，但我 们注意到，已有的实验或数值分析往往局限于较低 Pr 数流体，如硅熔体或Pr 数在 10 左右的硅油，少有分析高Pr 数流体热流体波的报道。最近，Teitel 等人[20]测量了Pr 数 从6.8 到 102 之间的硅油热毛细对流失稳的临界条件，但他们没有观察这些流体的热流 体波。本文中，我们通过线性稳定性分析确定了Pr 数较小的硅熔体(Pr=0.011) 和Pr 数 介于 6.7 至 57.9 之间的硅油在较大Pr 数范围内热毛细对流失稳的临界条件，获得了它 们的热流体波，特别是在高Pr 数下的热流体波。 1. 数学物理模型和数值方法 我们以如图 1 所示的物理模型为分析对象，它是一内径R =20mm 、外径R =40mm、深 i o 度d=1mm 的环形池，池内流体为硅熔体(Pr=0.011)或Pr 数从 6.7 到 57.9 的几种硅油， 它们的物性参数见表 1。环形池外壁维持高温 T ，内壁低温 T ，液层上表面为自由