中等雷诺数下稳态运动单液滴滴外传质数值模拟.pdf

中等雷诺数下稳态运动单液滴滴外传质数值模犷 李天文毛在砂 陈家铺 (中国科学院化工冶金研究所，北京 100080) 关键词 单液滴传质 数值模拟 贴体正交坐标 引言 在萃取工程中，单液滴的传质构成整个传质过程的基础，几十年来人们对此进行了多 方面的研究，以期得到适用范围广的经验关联式或数学模型，但是遇到了实验和数学上的 重重困难。与此相成对照的是由于计算机技术的飞速进展和数值计算方法的不断创新，计 算流体力学作为一种崭新的 “实验手段”，在单液滴流体力学研究方面取得了长足的进展， 这为单液滴传质过程的数值模拟提供了必要的基础。本文把计算流体力学的研究方法引入 单液滴的传质过程研究，对连续相传质阻力控制控制的传质过程进行了数值模拟，展示了 单液滴传质研究的新思路。 1坐标变换 单液滴(分散相)在无限不互溶流体介质(连续相)中的运动属于自由界面问题 其求解 区域几何形状比较复杂，难以用一般的方法精确地施加边界条件。为此计算中采用了正交 贴体坐标变换，通过协变Laplace方程((1)将液滴内外的求解区域变换成计算平面((}.7)上几 何形状规整的1x1的正方形区域。变换时采用Ryskin和Leal的强限制法和弱限制法to 在对滴外区域进行坐标变换时，采用了毛在砂和陈家请(21提出的有界畸变函数(2)，将滴外 物理区域取为有限但足够大，可提高计算精度。图1为滴外区域坐标变换示意图。所用的 Laplace协变方程为 a}lfa卦J+an(fIAan))=。，aCf卦韶蔚。 1() 有界域畸变函数: ’国家自然科学基金资助项目(编号No和中石化总公司资助项目，并得到清华大学国家重点 化学工程联合实验室部分资助 ~」96一 | 勺 1 f(},71)二粤(1-O.Scos7rl7) (2) P 而其物理意义为 ) 3 ( fV,71)二督 ) ’ ( 、二A-;}T-I}Wz-lT，、=衙l0+0)} a()Externaldomain (b)Dropinterior (c)Unitsquarein(fn) Fig.lThecorrespondenceofexteriorandinteriorofadeformabledropwithcomputationaldomain 上面公式当中h,h。为尺度因子，刀为调节因子用来决定滴外物理平面的大小。 2数学模型及求解 数学模型考虑液滴(分散相)处在以恒定速度U流动的另 一无限不互溶流体介质中，并 且以下假定成立: (I)流体为粘性、不可压缩的牛顿型流体，作轴对称流动; (2)浓度对体系物性和流动的影响很小，因而可以对运动方程和浓度方程独立求解; (3)相界面处两侧的浓度处于热力学平衡状态; (4)体系纯净，无表面活性物质干扰。 2.1流动方程 流函数一涡函数法求解滴内外区域的Navier-Stokes方程轴对称的自由液滴运动的无 因次控制方程为其〔中下标1指滴外，下标2指滴内): 397 Re, 1Fua,, a 。‘.、a},a 。‘)l