胀流型非牛顿幂律流体液滴铺展特性研究-苏州大学.docVIP

胀流型非牛顿幂律流体液滴铺展特性研究 鄂文汲 指导老师 魏琪 摘 要：为了研究胀流型非牛顿幂律流体液滴的铺展特性，我们利用广义Navier边界条件替代无滑移边界条件，对固体表面上薄液膜的一维简化模型采用润滑近似，构造了相似解法，采用数值方法计算了胀流型非牛顿幂律流体液滴的铺展速率、接触角和轮廓高度，得出了铺展半径、接触角和幂律指数的相关关系. 关键词：幂律流体；胀流型；液滴扩散 Abstract：In order to investigate the spreading of dilatant drops of power-law rheology , we compute the spreading rate of non-Newtonian drops using Navier-Stokes boundary condition instead of non-slip boundary condition and we construct a one-dimensional model on the substrate. To compute the spreading rate , apparent angles , height profiles and find the relationships between them , we use the thin film approximation , construct similarity solutions and use numerical methods. Key words : power-law fluids; dilatant fluid ; the spreading of drops 引言 液体在固体表面上的铺展，在诸如喷墨打印机、镀膜等许多领域中应用广泛，近年来在这个方面已经有大量的理论和实验研究。研究主要集中于牛顿流体。在非牛顿幂律流体方面，已有对完全浸润假塑性流体，部分浸润假塑性流体，完全浸润胀流型流体以及部分浸润胀流型流体的铺展规律的研究，在剪切变稀流体研究的基础上，本文对胀流型非牛顿幂律流体液滴的铺展特性做理论研究，给出了完全浸润胀流型液体铺展半径随时间演化的解析关系及动态接触角与接触线部分浸润移动速度关系的数值解，得出完全浸润胀流型流体液滴的铺展特性。 数学模型 根据幂率定律，假定有效粘度和切变率成幂指数关系：，m和λ的取值取决于流体的物性，当λ1时，流体称为剪切变稠流体，也称为胀流型流体。 牛顿流体有这样的问题——接触线佯谬，为解决这个问题要避免基板上的边界不滑移条件，引入滑移边界条件。 首先采用润滑近似建立固体表面上薄液膜的一维简化模型，构造一个对称的扩散解法，用来模拟液滴扩散。 在幂律流体中，切应力和切变率的关系为: (1) 这里，γ˙ij是切应变，将它定义为: (2) 是流体速度场，当λ = 1时,该流体即为牛顿流体。如果λ 1，代表剪切变稠流体。 采用润滑近似以得到运动方程。假设薄膜比它水平方向要薄的多，所以可以将运动近似看成水平方向的运动；同时忽略惯性效应，这样可以认为流体的流动仅取决于毛细作用力和粘性力间的平衡作用。 采用柱坐标，并假设液滴为圆形对称，在这个系统中，流体的速度有分量，且在基板处，和速度的r分量相比时，忽略它的z分量，即|u| ? |v|。假设应力主要是由于z方向上的水平速度的梯度决定的，假设分量和比其他任意方向的都大得多。 利用上面的假设，r分量的动量方程就简化为 ， （3） p 是压力，令z = h(r, t)为自由液面，然后将质量守恒方程用润滑近似写成 ， （4） U是z轴方向的水平速度的平均值。 假设液滴的尺寸比它的毛细长度要小得多，所以忽略重力，这样自由液面下的压力就可以近似写成 ， （5）是表面张力，是自由液面的平均曲率，是大气压力，采用润滑近似，将曲率近似写为 ， (6) 对式（3）积分，由于它满足自由液面接触角为零的应力条件，得到剪切应力和z