对于涂覆液体板间蠕变拖曳层流中涡流的分析1.pdf

对于涂覆液体板间蠕变拖曳层流中涡流的 分析1) 2) 刘鑫 刘赵淼 王国斌 （北京工业大学机械工程及应用电子技术学院，100124 ） 摘要 本文采用复值函数理论和有限元方法，对波形板和平板之间的牛顿流体和非牛顿流体的蠕变拖曳层 流过程进行数值研究，分析雷诺数、剪切变稀系数和阻通比的变化对于流体在波形板上的薄膜流动过程中 涡流产生、发展的影响。研究结果表明：雷诺数在一定条件下也会产生涡流，且会影响涡流的作用区域， 雷诺数不是影响涡流的主要作用因素在于流场中存在着过大的速度梯度；对于非牛顿流体，剪切变稀系数 不会对流场的形态产生过多的影响，但随着剪切变稀系数的增大，发生流线分离的位置在纵向上逐渐升高， 涡流作用的区域增大，作用强度增强，使得流动更加不稳定；随着平均板距的增大，临界波形度的绝对值 也随之增大，阻通比则是随之减小的，临界阻通比随着波形度的增大而呈现减小的趋势。 关键词 雷诺数，剪切变稀系数，阻通比，涡流 引言 分离流和涡流运动是工程系统和自然界中薄膜流动较为常见的不稳定流动现象。流动分离发展 形成涡流，涡流是分离流动的极限，涡流的运动、发展反过来又会影响流场。涡流和流动分离将影 响涂覆过程中的物质混合，同时影响涂覆过程中的化学反应速率和其中涉及到的热量和物质传输速 [1-2] 率 。 因此，有必要对流动中涡流的生成和演化做深入的研究。本文采用有限元法，通过精确数值计 算研究二维情况下固定波形板和运动平板间的拖曳层流，通过计算结果所得到的流线图，在已有的 研究分析基础上，分析牛顿流体和非牛顿流体在拖曳涂覆流动过程中的流动分离和涡流现象。 1 物理模型 研究的物理模型如图 1 所示，模型顶部为平板，底板侧面轮廓曲线为三角形，底板保持固定， 平面板以恒定速度 U 移动。计算模型将底板侧面轮廓曲线取为三角形余弦函数曲线，其中振幅为A ， 波长为 λ。两板块之间的平均距离定义为H 。不考虑重力。x 轴放置在波形板侧面轮廓线的水平对称 轴线位置，z 轴置于波形板的波谷位置，与 x 轴垂直。波形板无限延伸，流动几何条件可以假设为 r ~ ~ ~ r ~ ~ ~ r （1） ( ) ( ) v u x z e w x z e , x + , z 1) 基金项目：国家自然科学基金“薄膜蠕变流动特性的动力学分析” 资助；北京市人事局留学回国 人员择优资助项目资助 2) E-mail: lxxiabaili@emails.bjut.edu.cn 图 1 物理模型 图2 计算模型 ~ ~ ~ ~ x x 2πz λ, z 2πu λu, U w w, U , 对模型参数进行无量纲化， 连续性方程不变，则 ∂u ∂w + 0 . （2 ） ∂x ∂z 进一步地，对压力参数也进行无量纲化，使得 ~ 。经过无量纲化处理和数学变换后， p p U λ 2πη 将整个系统几何形态参数简化到两个无量纲参数[2]，波形板的波形度：a = 2πA/λ，以及两板块之间 的平均距离：h = 2πH/λ。实际上， 波形度可以理解为基底波形板几何形态的纵横比，平均板距可视 为涂覆通道宽度