基于Level set方法的双层流体热毛细对流的数值研讨.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093038 基于Level set方法的双层流体 热毛细对流的数值研究 周小明，黄护林，蔡亮，曹鑫徫 (南京航空航天大学高新技术研究院，南京，210016) (Tel：02584896183，Email:hlhuang@) 摘要：基于Level set方法建立了双层流体热毛细对流的数学模型，通过变密度二阶投影法求解控制方 程，C-N隐式技术用于扩散项更新，三阶龙格库塔技术用于对流项的更新，采用连续表面张力模型（CSF） 模拟Marangoni效应。三维数值模拟了微重力环境下双层流体系统中交界面变形与流场、温度场相互耦 合的热毛细对流，结果显示，在Marangoni效应的作用下，交界面在热端凸起，在冷端凹陷；随着 Marangoni数增大，双层流体交界面的变形率随之增大，对流强度也随之增大；交界面与壁面的接触条 件对交界面变形和流体的流型结构具有很大影响。 关键词：Level set 方法；热毛细对流；双层流体；界面变形，数值模拟 0引言 热毛细对流是由于流体自由表面温度分布不均匀而引起表面张力梯度驱动的流体 [1] 流动，这一对流现象存在于很多应用领域, 如晶体生长和镀膜等 。常重力环境下，热 毛细对流通常被自然对流所覆盖，而在微重力环境下，热毛细对流成为流体热对流的主 要形式，因此开展微重力环境下热毛细对流流动特性的研究对空间晶体生长具有重要意 义。然而，由于空间微重力实验非常昂贵，通过数值模拟的方法对热毛细对流进行研究 显得尤为重要。目前，在双层流体系统热毛细对流的研究中，很多研究者为了简化对自 [2-4] [5] 由表面的处理，都假设双层流体的交界面不发生变形 。Koster 的微重力实验结果表 明，微重力环境下多层流体热毛细对流的数值研究必须考虑交界面存在的变形和交界面 与壁面的动力学接触条件。虽然一些研究者已在双层流体热毛细对流的研究中考虑了交 [6-7] 界面的变形 ，但他们的研究中忽略了上层流体粘度，并仅限于两维研究。因此，本 文基于Level set方法建立了双层流体热毛细对流的数学模型，三维数值研究了微重力环 境下双层流体系统中交界面变形与流场、温度场相互耦合的热毛细对流。 1物理数学模型 物理模型如图 1 所示，其为一个长(x)为40mm，宽(y )为20mm ，高(z)为20mm 的长 方形腔体，特征长度L=40mm 。初始状态为静止的双层非混合流体系统，上下层分别为 等厚的流体，其中上层为Silicone oil 封闭液，下层为InB 熔液，左右壁面分别维持恒定 i 的温度为T 和T (T T ) 。对微重力条件下，物理模型的简化如下： c h h c 1）流体为不可压缩的牛顿型粘性流体，流体物性不发生变化； 基金项目：国家自然科学基金资助（No） 2 ）双层流体交界面受到表面张力梯度（即Marangoni 效应）的作用，流动仍是层流； 3 ）流体表面张力是温度的线性函数； 4 ）研究模型左右壁面为恒温条件，其它壁面为绝热条件。 图1 物理模型 基于以上假设，无量纲控制方程如下： u 0 (1) f u f 0