牛顿流体液滴成形过程与断裂原因,本文是关于过程方面学术论文怎么写和牛顿和成形和原因研究方面本科毕业论文范文.
过程论文参考文献：
摘 要:利用高速摄影机对实验中液滴成形过程进行拍摄,获取一系列成形过程的瞬时图片.运用Fluent软件中VOF(Volumeoffluid)方法,对低流量情况下毛细管末端液滴成形过程进行数值模拟,采用几何重构方法(Geometricreconstruction)对液滴成形的气液界面进行捕捉,并与高速摄影获得的成形实验结果比较,结果表明数值模拟结果与实验一致,表明数值模拟的可行性.在此基础上,计算出液滴整个成形过程中所受合力的变化,分析了成形过程液滴内部的压力,结果表明:液滴成形过程可分为成形、颈缩、断裂三个阶段;断裂阶段重力和表面张力将达到平衡,平衡之后重力占主导地位,微小的体积增量使得颈缩处急剧变细,液滴断裂;成形和颈缩阶段沿液滴中心线压力逐渐增大,断裂阶段颈缩处压力最大,达到３６０Pa.
关键词:液滴成形;几何重构;压力;合力;断裂
中图分类号:O３５９＋ ．１ 文献标志码:A 文章编号:１６７３Ｇ３８５１(２０１８)０５Ｇ０３１９Ｇ０６
０ 引 言
液滴成形一直是国内外学者的研究焦点,其研究成果广泛应用于日常生活和工业生产中,如喷墨印刷[１Ｇ２]、增材制造[３]等.为了充分认识液滴成形及其断裂的过程,国内外学者进行了大量研究.Zhang等[４]研究了影响初始液滴和卫星液滴的体积以及颈缩线长度的参数,包括流体粘度、流速、电场等,结果表明:粘度和流速增加时形成的颈缩线长度相应增加,生成的初始液滴和卫星液滴的体积相应增加;随着电场强度的增加,初级液滴的体积减小,而颈缩线的长度和直径以及卫星液滴的体积相应增加.Shi等[５]通过实验研究发现,液体粘度对液滴形状和长度有着显著的影响,液滴在高粘度下更容易脱落成单个液滴.王定标等[６]使用VOF模型对液滴形成过程进行了数值模拟,结果表明:不同阶段液滴内部速度由不同条件决定,颈缩阶段由入口速度决定,轴线压力线性增长;断裂阶段由重力决定,轴线压力的变化与速度变化情况一致,最后得到液滴完整长度与管内径的比值和韦伯数(We)之间的关系式.Ambraneswaran等[７]研究发现:毛细管力占主导地位时,１D算法相比较２D算法的相对误差不超过几个百分点;惯性力或粘性力大于表面张力时,１D算法的相对误差最高可达１５％;Oh(奥内佐格数)较小时,１D 算法由于无法精确地预测流体的界面而导致相对误差较大.Eggers等[８]使用泰勒展开来简化二维模型并获得一维有限元模型.一维有限元法具有以下两个优点:一是节省计算时间;二是在初始液滴破裂后可以不间断地计算卫星液滴.Schulkes[９]计算了液滴体积随着时间的推移稳定增加时喷嘴悬挂液滴末端的演变情况,研究发现:当液滴下降变得不稳定和增长速度变慢时,连接两端的颈缩线很快出现.当增长速度变大时,液滴及时分离,且能够计算断裂后的液滴体积.Xu等[１０]研究正弦函数流量下无量纲数(We,Oh)对液滴动力特性的影响,给出了不同流向不同We 数时液滴的断裂情况.贺丽萍等[１１]运用了VOF(Volumeoffluid)方法和几何重构方法成功模拟了液滴成形的过程.研究发现:液滴内部的涡环将经历生成Ｇ湮灭、合并Ｇ分开的过程,得到其与液滴形状的关系,得到了压力与液滴不同部位的关系.
液滴成形的影响参数(粘度、流速等)和形成特征(初始体积、卫星体积、流线长度等)一直是众多学者的研究对象,而液滴成形过程中所受力的作用和其内部流场、压力的研究较少.本文以毛细管末端液滴为研究对象,通过研究低流速下液滴内部流场、压力及作用在液滴上的力,从理论上分析液滴的成形过程和解释液滴断裂原因,为液滴在科学研究和工业领域上的应用提供更多的参考.
１ 实验和数值研究方法
１．１ 实验方法
实验研究了流体从垂直毛细管喷嘴末端以恒定流速流入空气中成形液滴的过程.在初始液滴脱离之前拍摄一系列成形过程的图像,实验台的示意如图１所示.实验装置由不锈钢毛细管(外直径３mm,内直径２mm)、蠕动泵(BT１００Ｇ１F)、高速摄影机等组成.实验中管子为不锈钢毛细管,管子末端加工成平面,并使其平直且内部和边角无毛刺.所选择的管子长度超过５cm,具有较大值的长度与直径之比,以保证所述管子出口上游部分有充分发展的分布均匀的流速.实验设备最重要的就是高速摄影系统,包括高速摄影机、光源、视频系统以及相关硬件,使用具有最高分辨率为１９２０×１０８０ 像素的高级CMOS传感器的高速摄影机(PCO．dimax HD)拍摄液滴成形序列图像,使用２０００fps进行拍摄,以便在液滴脱落之前能拍摄到颈部破裂的图像.背光由５００W 灯泡提供.此外,该系统还可以在极短的时间内一边显示一边存储液滴外形.由蠕动泵提供连续且稳定流量,该蠕动泵的流量范围为０．２μL/min至５００mL/min,流量控制操作流程如图２所示.
在该实验中流量调为１mL/min.毛细管插在流体中,整个装置除了泵以外都放在隔震台上,质量分数为５０％的甘油溶液的性质[４]见表１所示.为避免实验中存在气泡,要彻底清理毛细管,消除影响液体流速的不确定因素.液滴成形在低流速情况下进行,在拍摄前使液体先流动１h,然后进行拍摄,这样能得到可靠的液滴成形过程照片.
１．２ 控制方程
本文借助Fluent商业软件中的VOF模块开展数值模拟研究.整个计算域中不存在热的交换,设置主相(Primaryphase)为气相,第二相(Secondaryphase)为液相,并且假设气相和液相均为黏性不可压的.流场由NierＧStokes方程控制:
其中:ρ 表示密度;v 表示流场中速度的矢量;p 表示压力;μ 表示动力粘度;σ 表示表面周线单位长度上的表面张力值,即为表面张力系数;C 表示VOF方法定义的体积分数.本文选用的是连续表面张力模型[１２](Thecontinuumsurfaceforcemodel)加入到动量方程来处理表面张力,κ为界面曲率,表现形式如下:
其中:C＝０表示单元网格中都是气相;C＝１表示单元网格中都是液相;０＜C＜１表示单元网格中存在气相和液相的交界面.每个单元网格内的粘度和密度由式(６)和式(７)得到:
１．３ 边界条件与计算方法
数值模拟的几何结构与实验一致,因为所研究的对象为三维轴对称模型,故将其简化为二维对称模型,计算域和网格如图３所示.边界１为速度进口(VelocityＧinlet),U ＝０．００５３m/s,相当于１mL/min流量.边界２、３、４为壁面(wall),其中壁面２和３的接触角设为６０°,壁面４的接触角设为１７９°,５、６、７设为压力出口(PressureＧoutlet),８ 设为对称轴(Axis),９为气相和液相的界面.边界６和７的交点设为参考压力点,压力值设为１０１３２５Pa.
对重力作用下毛细管末端液滴成形过程的控制方程进行离散求解,用PISO(Pressureimplicitwithsplittingofoperators)算法对速度压力耦合进行迭代计算,PISO使用了预测—修正—再修正的处理过程,可加快单个迭代步的收敛速度.压力场的求解选择的是体积力分数(BodyＧforceＧweighted)计算方法.选择乘方格式(Powerlaw)求解动量.对于气液界面的处理,采用的是VOF中基于PLIC(PiecewiseＧlinearinterfacecalculation)的几何重构法.在ANSYluent中这个方法是最准确的,适用于一般的非结构化网格,Youngs[１３]对几何重构方法进行了详细描述.
１．４ 液滴受力分析
奥内佐格数(Oh＝μ/(ρRσ)１/２,R 为管外半径)为度量黏性力和表面张力效应之比的无量纲数,韦伯数(We＝ρUD/μ,D 为管外直径)为度量惯性力和表面张力效应之比的无量纲数.由于低流速下奥内佐格数和韦伯数量级都在１０－２,惯性力和粘性力相比表面张力非常小,故可以忽略,液滴主要受重力和表面张力作用,
２ 结果与讨论
２．１ 不同成形阶段分析
通过对实验和模拟图片的二次处理得到了液滴气液界面和长度的对比图,气液界面如图４所示.由图４可知,液滴成形的每个阶段实验和模拟在外观上非常相似.图５比较了实验和模拟液滴随时间变化的长度,实验和模拟的液滴成形长度很接近,即实验与模拟有很好的一致性,验证了数值模拟的可行性.
图４显示了液滴成形到断裂的全过程,成形阶段液体逐渐流入,悬挂在管口处的液滴体积逐渐增加,如图４(a)—(b)成形阶段所示.当液滴的体积到达临界体积时,液滴运动开始不稳定,即出现颈缩现象,如图４(b)—(g)颈缩阶段所示,液滴的大部分在喷嘴口处悬挂,颈缩处直径逐渐变小.断裂阶段液滴体积微小的增量将使液滴不稳定,体积微小增量将使颈缩直径迅速变小致其断裂,如图４(g)—(i)所示.
２．２ 受力分析
图６是气液界面和接触角变化示意,接触角如图６所示(虚线和虚水平线之间),图７为绘制的液滴湿润直径和接触角的变化情况.图７中０~０．４６s是一个短暂且变化剧烈的阶段,仅仅持续约０．４６s,湿润直径和接触角变化明显.液体从毛细管顶部流出形成弯液面,悬挂在毛细管壁.接触角急剧下降并达到临界值,湿润直径逐渐增加,考虑到毛细管壁厚度有限,湿润直径达到毛细管的外缘(图６中的线１的界面)时不再增加.
图７中０．４６~２．５０s这个阶段在时间上占据了大部分成形过程.液滴的形状从弯液面转化为类似葫芦的形状,如图６所示的线２和线３.接触角的值基本保持不变,湿润直径也仍然保持着毛细管外半径的值.时间２．５００~２．７４８s也是一个短暂且变化剧烈的阶段,持续大约０．２４８s.接触角从８５°左右快速增加到近１２０°,湿润直径也是快速减小.如图６中的线４变化到线５,然后液滴脱落成一个单独的液滴.
重力和表面张力的合力如图８所示.颈缩现象出现之前重力和表面张力的大小由式(８)和式(９)计算得出,颈缩现象出现之后,重力和表面张力由式(１０)和式(１１)计算得出.由图８中可知,重力和表面张力的合力在时间为２．７２４s时为０,即重力和表面张力大小相等.当时间超过２．７２４s后,重力大于表面张力,使得液滴不断拉长,颈缩部位越来越细直至液滴断裂,液滴脱落.
２．３ 液滴内部压力分析
图９为成形液滴中的压力场和流线的变化.每个时刻流线在中心线的右侧显示,压力场及其等值线的值显示在其左侧.从图９(a)—(e)明显可以看出,与喷嘴出口附近相比,液滴顶端处的压力更大,各个时间段从入口处到液滴顶端压力逐渐增加,即呈现静水压的状态.从图９(f)可以看出,在断裂阶段颈缩处急剧变细,压力显著增加,压力达到３６０Pa,明显高于液滴顶端和管口附近的压力.此时颈缩处以下近球部分在压差的作用下有向下断裂的趋势,而颈缩处以上部分有向上运动的趋势,致使颈缩处越来越细.最后当颈缩处的直径趋近于零时,颈缩处液滴的径向截面上拉应力趋向于无穷大,液滴断裂.图９(a)—(e)中随着液滴的长度增加其内部的流线形成的流环也逐渐加大,随着颈缩直径越来越小,如图９(f)所示,颈缩部分下方流入的液体体积变少,流环变小,直至颈缩处断裂,下方液滴内各点流线到达一致自由下落.
３ 结 论
本文通过高速摄影机对液滴成形过程进行拍摄,得到其演变的一系列图片,并对其进行了数值模拟,从受力情况和液滴内部压力解释了液滴断裂的原因,得到如下结论:
a)与实验结果比较验证了VOF方法模拟液滴成形和断裂的可行性.
b)根据液滴形态变化将成形过程分为三个阶段:成形阶段,颈缩阶段,断裂阶段.在断裂阶段t＝２．７２４s时重力和表面张力达到平衡,之后重力占主导地位,使得液滴不断拉长,颈缩处越来越细直至液滴断裂,液滴脱落.
c)液滴即将断裂时,颈缩处压力最大,变化急剧,其值达到３６０Pa,颈缩处上部分液体和下部分液体内存在压力差,使得颈缩处直径越来越细,颈缩处液滴的径向截面上拉应力趋向于无穷大,液滴从颈缩处断裂.
括而言之:本文是一篇关于经典过程专业范文可作为牛顿和成形和原因研究方面的大学硕士与本科毕业论文过程论文开题报告范文和职称论文论文写作参考文献.
张启云:宁流千滴汗不坏一粒粮
余锦洲 李 莉日前,安徽省铜陵市粮食局北斗山粮库主任张启云荣获全国粮食工作先进个人称号 安徽省仅有5 名,这是铜陵市粮食局建局以来个人获得的最高荣誉 张启云1993 年7 月毕业于安徽省蚌埠粮食学校仓.
小儿血液灌流中出现血流量不足的原因分析和护理
摘要目的分析小儿血液灌流中血流量不足症状出现的原因,并探讨护理对策 方法将我院2010年10月2015年10月收治的小儿血液灌流患儿16例作为研究对象,回顾性分析影响患儿血液灌流的因素,并针对其原因采.
硅胶泪道引流管泪道环路插管治疗泪小管断裂效果观察
摘要目的针对硅胶泪道引流管泪道环路插管治疗泪小管断裂效果展开讨论,为日后的临床治疗提供参考与指导 方法将我院收治的泪小管断裂患者作为本次实验研究对象,所有患者的入选年份,均集中在2015年5月2017.
滴答秋雨等
滴答秋雨刘福宇夜临年幼的女儿哭闹着久久不愿睡去她害怕闭上眼睛整个世界都变黑了滴答,滴答窗外传来落雨的声音外婆说你不要哭了你把天上住的老爷爷都惹伤心了我说宝贝你听,滴答,滴答…&helli.