基于传热特性与撞击流动进行液滴与壁面撞击分析.docx

基于传热特性与撞击流动进行液滴与壁面撞击分析喷雾冷却作为一种高效换热技术受到国内外学者的广泛关注，喷雾过程中大量液滴撞击壁面，传热包括液滴与壁面、液滴与表面液膜，是属于复杂的多相流传热问题目前关于喷雾冷却的研究主要以传热计算分析为主，仅限于从宏观方面去分析传热影响因素及可能的原因，而对于传热机理的研究还未有统一的定论；由于喷雾液滴直接与壁面接触，因此涉及微观领域的喷雾液滴撞击壁面传热分析作为一种较为先进和有效的技术手段，将对喷雾冷却机理的研究起到关键性的促进作用[1,2]该研究也引起了国内外学者的广泛关注，并从机理分析、实验验证及数值模拟等方面开展了液滴撞击表面的动力学及热力学特性研究以帮助理解喷雾冷却强化换热的原因[3,4]梁刚涛等人[5]搭建单个液滴撞击壁面实验台，利用高速摄像机记录液滴形态变化，然后与VOF方法仿真得到的结果进行对比验证，结果表明该方法得到了结果很好的与实验结论相匹配，说明液滴撞击研究用VOF方法是可行的廖斌等人[6]也使用VOF方法建立单个液滴撞击壁面三维模型，分析撞击过程中液滴运行轨迹，结果发现该方法可以直观的将液滴在空中的运行轨迹表示出来，并且当液滴撞击表面发生变形的过程也可以观察到。

周龙玉[7]研究的重点在于分析液滴撞击过程接触角、液滴初始直径以及不同撞击表面带来张力系数对撞击过程的影响，他先使用高速摄像机拍摄其中一组工况以验证VOF方法的准确性，然后就参与该仿真方法改变不同的参数变量得到影响液滴撞击壁面的因素刘红等人[8]采用数值仿真手段分析液滴撞击初期及扩散阶段特性参数对表面形成的液膜产生的影响，重点考虑了液滴初始动能、表面张力系数以及粘度系数等隋涛等人[9]采用VOF数值仿真方法评估当改变压力环境和液滴初始速度时对撞击壁面造成的影响蒋昌波和邓斌等人[10]为了精确的描述单个液滴在下落及撞击过程中流动情况，使用基于Navier-Stokes方程的VOF方法成功捕捉到了液滴撞击变形的整个过程通过调研国内外学者目前的工作发现，关于液滴撞击壁面的研究以仿真为主，辅以实验并已经取得了一些成果，但是涉及动力学特性研究较多，而结合传热的分析还需要继续推进，尤其结合喷雾冷却等实际应用背景下的研究；因此本文在此基础上采用LevelSet-VOF相界面追踪的方法，搭建液滴撞击固体壁面的数值模型，重点分析液滴尺寸、初始动能及表面液膜厚度对液滴流动及传热的影响1、液滴撞击壁面模型及模型验证1.1控制方程及传热模型考虑目前的VOF模型存在追踪液滴运动轨迹不准确的弊端，因此本文在建立仿真模型过程中加入LevelSet方法以液滴边界相界面为追踪目标，提高了追踪精度，液滴撞击过程是液相与气相相结合的过程，因此涉及两相流的问题，其控制方程为：式（1)-（2）中，u为速度；φ为LevelSet函数；p为压力强度；表面张力系数；为表面张力项；为黏性应力张量；为曲率；为Heaviside函数。

Heaviside函数表达式如下：式（3）中，h为网格大小，黏性应力张量，应变率张量S表达式为：在模型中将液滴撞击整个过程分成不同的区域，比如在空中运动及撞击固壁灯，而不同区域r和m可由Heaviside函数计算：液滴在撞击固壁过程中除了形态发生改变存在着热量的传递，接下来将对液滴撞击壁面过程传热模型进行描述液滴接触高温壁面时蒸发吸热，液滴表面有一层饱和空气边界层，水蒸气与液滴之间的传热传质取决于计算区域内主体蒸汽与边界层内饱和蒸汽间的温差，以及蒸汽分压力的大小，其中，边界层的蒸汽分压力取决于液滴表面温度[11,12]在液滴撞击壁面过程中单位时间内液滴热量等于壁面热传导以及液滴自身蒸发吸热的差值如式（7）所示：式（7）中，md为液滴质量，Td为液滴温度，cp,d为液滴的比定压热容，As为液滴与壁面接触面积，K为液滴与加热壁面间表面传热系数，Ts为加热壁面温度，r为液滴汽化潜热，qm,d为液滴蒸发速率液滴受热蒸发时蒸汽离开液滴表面的质量通量如式（8）所示：式（8）中，Md为液滴的摩尔质量，hD为表面传质系数，ρv0为液滴表面蒸汽密度，ρv为环境蒸汽密度根据理想气体蒸汽压力与浓度的关系如式（9）所示：式（9）中，Pv0为液滴表面蒸汽压力，Pv为环境蒸汽压力，Tg为环境蒸汽温度，R为水蒸汽摩尔气体常数。

根据质量守恒定律有：式（10）中，ρd为液滴密度1.2计算区域及边界条件设定撞击固体壁面的液滴尺寸模型如图1所示计算区域为长方形，区域大小为3mm×6mm，使用四边形均匀网格液滴撞击固体壁面区域顶部设定为速度入口，两边设定为压力出口，底部设定为固体作为壁面处理；开始时刻液滴刚刚与壁面接触，液滴圆形底部刚好与固壁相切；初始直径为D0，初始速度为U0，环境压力0.1MPa，考虑重力g的影响，环境气体为空气，壁面为光滑铜表面且温度恒定考虑到液滴撞壁后的运动过程较复杂，对模型作出如下假设和简化：（1）初始液滴为球形；（2）忽略液滴与环境气体的热辐射效应；（3）忽略环境气体与液滴间的剪切力；（4）气液两相均为不可压缩流体且环境气体物性参数为恒定值图1初始时刻液滴撞击壁面模型网格数量及精细程度会影响仿真结果的准确性，但是对用于计算的机器性能要求也较高，因此需要兼顾结果准确性和计算设备的性能确定网格的梳理，因此需要对网格进行无关性验证；分别采用4万5千、8万、12万5千的结构性网格，液滴初始速度为0.1m/s，液滴直径0.5mm，分别计算四个工况下液滴撞击铺展因子变化规律，结果如图2所示从图2可以发现，液滴撞击过程也是液滴铺展的过程，因此随着时间推移，铺展因子均是增大的，横向对比发现，当网格数密度增大，铺展因子也随着增大，当网格数量达到12万5千时，铺展因子增加趋势开始平缓，此时网格数量对液滴铺展结果影响较小。

验证了液滴撞击壁面结果与网格数之间无关联性图2网格无关性验证综合考虑计算精度及计算效率本文模拟采用12万5千的网格数量1.3VOF模型验证Liang等[13]在Leidenfrost温区（壁面温度Tw>TDL）对水滴撞击壁面的运动特性及传热行为进行了实验研究为验证复合LevelSet-VOF的可行性，选取其中一组实验数据（工况参数为：壁面温度384K，液滴初始直径2mm，液滴初始温度20℃，We=2）和模拟结果进行对比图3示出液滴撞击光滑壁面的过程，对比实验结果发现，模拟液滴撞壁后的运动行为与实验结果较吻合，这说明本文的数值模拟方法应用于液滴撞击过程形态变化是可行的图3实验结果与模拟结果对比2、单液滴撞击壁面运动特性研究本文以水滴为对象，研究其在不同初速度、初直径撞击不同壁面温度的运动特性，采用的相关参数见表1表1VOF仿真计算相关参数图4初始速度0.1m/s、5m/s，初始直径0.2mm、0.5mm撞击高温（373K)壁面的（a)液滴铺展形态及（b)铺展因子随时间的变化曲线图4分别给出了初始速度0.1m/s、5m/s，初始直径0.2mm、0.5mm撞击高温（373K>TDL）壁面的液滴运动形态及铺展因子随时间的变化曲线。

从图4（a）发现，直径为0.2mm的液滴以速度0.1m/s与壁面接触时，液滴底部由于先接触过热的壁面，因此先蒸发，液态变成气态，生成空穴，同时，液滴底部沿着壁面水平方向向四周运动铺展，形成薄液膜在373K的壁面上快速蒸发，在向上的蒸汽压力及液滴与壁面之间表面张力共同作用下驱使液滴有回缩的趋势，液滴悬浮在壁面上，最终导致液滴反弹脱离壁面；直径为0.5mm的液滴以速度0.1m/s与壁面接触时，与直径0.2mm的液滴相比均受到向上的蒸汽压力及表面张力作用使得液滴先铺展后回缩，在回缩开始阶段液滴在壁面铺展成中间液膜薄四围厚的圆盘状，且随着时间推移中间越来越薄，受表面张力作用导致中间断裂，形成破碎的更小的小液滴，余下的液滴继续回缩直至形成稳定的液滴[14]；当液滴初始速度继续增加至5m/s时，液滴在0.8ms内已经铺展破碎完全，这是因为此时液滴动能较大，惯性力起主导作用，液滴撞壁后与壁面局部区域持续接触，此时又因为液滴持续蒸发导致大量二次液滴溅射，并且由于液滴铺展时受到周围气体的阻滞作用，边缘的小液滴会向上运动，液滴的铺展行为演变为铺展破碎由于喷雾冷却研究中有一部分是分析液膜的厚度，该特性与液滴撞击壁面的铺展因子有极大关系，因此引起众学者的关注。

本文也将从铺展因子变化规律上进行相关的研究分析从图4（b）可以看出，当液滴撞击壁面后将出现三种形态铺展行为（反弹、回缩破碎、铺展破碎），对应的其铺展因子的变化规律也各有不同：当撞击壁面的液滴直径和初始动能均较小时，壁面的表面张力能起主导作用，因此液滴铺展因子随时间先增加后减小，对应液滴撞击先后经历铺展、回缩、反弹阶段；发生回缩破碎的液滴其铺展因子在1ms时增至最大，而后增加趋势减缓对应着液滴先铺展至最大之后开始回缩；当液滴动能足够大时，液滴在短时间内铺展至最大，液滴发生铺展破碎，铺展因子随时间的增幅增大3、液滴撞击壁面传热特性研究在喷雾冷却过程中最关心的是传热能力的问题，因此本文针对液滴撞击壁面进行传热分析，以帮助理解喷雾冷却传热强化的机理，选取液滴撞击壁面过程特定参数如液滴直径、初始速度及壁面液膜厚度进行分析3.1液滴直径与速度对传热的影响设定液滴初始直径分别为0.2mm、0.5mm、0.7mm，其他计算条件为液滴速度0.5m/s，壁面液膜厚度0.02mm，壁面热流密度300W/cm2，采用VOF方法进行仿真计算图5不同直径液滴以32.11m/s速度撞入液膜时表面传热系数变化趋势图5为不同直径液滴撞击到壁面时传热系数的变化规律；从图中可以发现，随着撞击时间推移传热系数均不断增大，为了得到液滴直径对传热的影响，横向对比三个直径下传热系数，选取撞击时间10ms为节点，直径0.2mm的液滴以0.1m/s速度撞击液膜，表面传热系数为5.63W/cm2·℃；当液滴直径增大到0.4mm，表面传热系数达到5.97W/cm2·℃；当液滴直径继续增大到0.7mm，表面传热系数达到6.17W/cm2·℃。

上述现象表明液滴撞击壁面时如果初始直径较大将会促进传热；原因在于相同撞击时间，液滴直径大，其与壁面接触的面积就大，铺展更加均匀，溅射的水花较大，从而带走的热量增加图6初始直径0.7mm液滴以不同速度撞击液膜时换热系数变化趋势图6为初始速度分别为0.1m/s、0.5m/s和5m/s的液滴撞击到壁面上时传热系数变化规律；从图中看出液膜表面的传热系数随液滴速度的增大而增大产生上述现象的原因在于液滴在撞击壁面过程中除了受到壁面表面张力还受初始速度带来的动能影响，当液滴速度大表明液滴初始速度大，当撞击壁面时有利于脱离表面张力对液滴的束缚，扰动加强，更有利于促进热量迅速被液滴带走，并且液滴初始动能大有利于液滴在短时间内就迅速在壁面铺展，增大了换热面积；这些均可以提升传热能力；该结论与文献[12]中的实验得到的结论基本一致3.2液膜厚度对传热的影响本文对单液滴撞击厚度分别为0.02mm、0.05mm和0.08mm的液膜过程进行了VOF数值仿真，设定液滴初始直径统一为0.2mm、初始速度为0.1m/s图7液滴撞击不同厚度液膜时表面传热系数随时间变化趋势图7为在不同液膜厚度下液滴撞击时表面传热系数随时间的变化情况，从图中看出在液膜厚度为0.02mm、0.05mm、0.08mm情况下，整个壁面的表面传热系数均随时间的增大单调递增加。

并且当接触时间超过0.4ms后表面传热系数递增趋势越大；这是由于两个方面原因所造成的：初始时刻，液滴撞击液膜厚度薄的壁面时，由于液膜自身对液滴的束缚力较小，因此液滴撞击壁面带来的扰动较大，此外液膜厚度薄，液滴在表面铺展越。