水平均质表面上液滴形状的数值模拟.pdfVIP

中固T程热物理学会 传热传质学 学术会议论文 编号：063268 水平均质表面上液滴形状的数值模拟 朱恂范勇廖强木 蕈庆人学丁程热物理研究所重庆400044 Tel：023—651 02474 Emaii：lqzx@cqu．edu．cn 摘要本文对水、I，．均质表血．I：静止液滴的形状进行J，数值模拟。对给定体积的液滴，由假定的初始自由 表面…发，采用迭代算法．使计算过程中系统的总能量／1i断减小，最终扶得静止液滴的彤状。分别模 拟．r接触角为75。．90。和120。时水’I，均顾表面j：的液滴形状。 ‘ 关键词水’|，均质表面液滴形状能量最小化‘ 1引言 在科学研究和【：业实践的许多领域中，都涉及剑了关丁静I卜液滴形状的问题。同体 表面上重力和表面张力作用卜．的液滴，其自由液面的结构可以用Young—Laplace方程来 描述。…住轴对称情况卜。可以简化为～个常微分方程，Paddayl2‘，Brien¨’等人已经给 山了满意的结果。针对二维情况，国外已有学者虑川有限元的数值方法进行模拟，Broun 等Hl住一个具有径向润湿区域的斜板上模拟了与之接触的液滴的平衡形状。’ 2数学模型 图l：均质表面上的静J}：液滴 考虑如图1所示的I卉卜·液一气系统，系统总能姑U的表达式如卜．： U2 “) 上式中，Q、Q2表示液滴和空气·’i据的体积：P、pp分别为液体和空气的密度：r、 ‘，、k和盯、矿。，、仃黔分别表示气液、川液、州气界面和各界面的表面张力系数： g为重力加速度。 根据能鼙最小原理，当液滴住均质表面上处丁．平衡状态时，系统的总能苗取得最小 值。即【，一min时，可以获得液滴的形状。 式(1-)应满足以卜_条{；，l：： (1)气液不相溶：万k施=-,3L艘P (2) 1205 (2)液滴体积守恒：y(Q)=L撒=Vo=const (3) (3)I州体表面积不变：万}．皿，=一万}订0+ (4) 二相接触线上的表面张力满足Young力程： 。 (5) 、 cosej=b罨s—os11}o 式中，0为Young接触角。 对(1)式取微分，将(2)～(4)式代入得： (6) 记 —Osl--—Ogs)d‘， (7) 仃 盯：可(Q，r，L，)：L土里三≠至2弘dQ+}dr+t 可以看出，6U：面。当U一_驭最小值时，U也取得最小值。丁．是问题转化为求可 的最小值，即U—}min。 记 ～bsi一仃gsl 6：鱼二丝超’Z2。 由式(5)可知Z=一cosO。 在图l所示的系统中，当同液界面r确定时，Q和L，便唯一的确定-卜．来，于是(7) 式中