液滴偏心碰撞的多体耗散粒子动力学研究.pdf

Modeling and Simulation 建模与仿真建模与仿真,2023,12(4),4213-4222 Published Online July 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.124384 文章引用文章引用:刘华建,曹达敏.液滴偏心碰撞的多体耗散粒子动力学研究J.建模与仿真,2023,12(4):4213-4222.DOI:10.12677/mos.2023.124384 液滴偏心碰撞的多体耗散粒子动力学研究液滴偏心碰撞的多体耗散粒子动力学研究 刘华建刘华建，曹达敏曹达敏 上海工程技术大学航空运输学院，上海 收稿日期：2023年4月13日；录用日期：2023年7月22日；发布日期：2023年7月28日 摘摘 要要 液滴碰撞细杆的现象在日常生活中随处可见液滴碰撞细杆的现象在日常生活中随处可见，本文采用了多体耗散粒子动力学本文采用了多体耗散粒子动力学(MDPD)模拟液滴以不同模拟液滴以不同角度偏心碰撞圆柱形杆角度偏心碰撞圆柱形杆，研究了不同碰撞速度，不同碰撞角度以及不同偏心率对于液滴薄片比例和临界研究了不同碰撞速度，不同碰撞角度以及不同偏心率对于液滴薄片比例和临界速度的影响。结果表明，分裂滴落的临界速度随着偏心率增大先减小而后增大，不同角度下，临界速度速度的影响。结果表明，分裂滴落的临界速度随着偏心率增大先减小而后增大，不同角度下，临界速度的改变趋势类似，接触细杆的临界速度随着偏心率增大而减小。的改变趋势类似，接触细杆的临界速度随着偏心率增大而减小。关键词关键词 多体耗散粒子动力学多体耗散粒子动力学，数值模拟数值模拟，液滴碰撞液滴碰撞 A Many-Body Dissipative Particle Dynamics Study of Droplets Eccentric Impacting on Fiber Huajian Liu,Damin Cao College of Air Transportation,Shanghai University of Engineering Science,Shanghai Received:Apr.13th,2023;accepted:Jul.22nd,2023;published:Jul.28th,2023 Abstract The phenomenon of droplet impacting on the fiber can be seen everywhere in our daily life.In this paper,we use Many-body Dissipative Particle Dynamics(MDPD)to simulate the droplet eccentric collision with a fiber at different angles.The results show that the critical velocity of splitting dropping decreases first and then increases with the increase of eccentricity.The changing trend of critical velocity is similar at different angles,and the critical velocity of contacting fiber de-creases with the increase of eccentricity.刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4214 建模与仿真 Keywords Many-Body Dissipative Particle Dynamics,Numerical Simulation,Droplet Impact Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 液滴碰撞圆柱体表面在自然界和日常生活中是很常见的现象，输电线路在寒冷天气情况中的结冰问题便是其中的一种，一直困扰着许多国家和地区，而防结冰是受关注的有效方法之一。液滴碰撞固体的实验一直被用于探究表面特性和液滴动力学，通过研究液滴碰撞圆柱体表面的铺展特征，对于输电线路防除冰的问题能够提供有效的应对方法和思路。已经有许多文章关于液滴碰撞固体进行了许多的实验和数值研究。在实验上，Banitabaei 1研究了液滴碰撞球形固体，发现在疏水性表面上，液滴扩散时会形成锥形薄片，随着韦伯数的增加，薄片拉伸长度和锥角会相应增加。Zhang 2发现研究了对液滴施加不同碰撞速度撞击不同直径比的超疏水圆柱体表面所产生的影响。Khurana 3研究了液滴撞击不同直径的亲水和超疏水圆柱表面上后的形态和动力学，提出了用于预测给定圆柱体表面上的润湿系数的理论模型。Liu 4研究了液滴碰撞亲水、疏水和超疏水圆柱形表面，发现圆柱体直径越小，液滴越容易发生破碎，反之会促进液滴的回缩和反弹。在仿真模拟方面，Wang 5通过 MDPD 方法，探究液滴碰撞圆柱体的冲击动力学特性，发现液滴碰撞疏水表面更容易引起接触线的回缩。Wang 6等人使用 MDPD 方法模拟了液滴碰撞不同倾斜角度的细杆，首先他们通过与实验比较来验证方法的可行性，然后研究了液滴偏心率、杆的倾斜角和润湿性对碰撞结果的影响，发现液滴临界速度、润湿长度、接触时间和液滴捕获率均随细杆倾角的增大而增大，随偏心距增大而减小。此外，他们还发现细杆的表面润湿性对接触时间的影响不大，但对临界速度的影响较大。液滴在圆柱体上的碰撞已经开展了广泛的研究，但大多数研究都集中于液滴以中心位置碰撞，对于液滴以一定角度的偏心碰撞涉及较少。在本文中，通过多体耗散粒子动力学的方法对液滴偏心碰撞细杆进行了模拟，系统了解碰撞角度和偏心率对于碰撞结果的影响。2.方法和模型设置方法和模型设置 2.1.多多体耗散粒子动力学方法体耗散粒子动力学方法 本文使用多体耗散粒子动力学7 8 9 10的方法进行仿真，是基于经典的耗散粒子动力学(DPD)改进而来的一种无网格数值方法，弥补了 DPD 方法不能解决液体和气体共存系统问题的不足。通过粒子间的相互作用来模拟系统，相互作用力ijF为三种不同力的总和，分别为保守力CijF、耗散力DijF和随机力RijF，方程如下：()()()CCijijijijdijijFAreBre=+(1)()()DDijijijijijFreve=(2)()RRijijij ijFre=(3)公式 1 中 A 和 B 分别为吸引力系数和排斥力系数，Open AccessOpen Access刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4215 建模与仿真 ()1ijCijcrrr=(4)()1ijdijdrrr=(5)是具有不同截止半径的权函数，ijr为粒子 i 和粒子 j 的之间的直线距离，ije为两个粒子之间的相对位置方向，密度权函数()ijr为i和j的总和，其公式为：()231512ijijddrrrr=(6)公式 2 中为耗散力系数，负号代表力的方向与粒子相对运动方向相反，()Dijr也是关于ijr的权函数，ijijvvv=，为粒子之间的相对速度。公式 3 中为正的随机力系数，()Rijr与()Cijr、()Dijr相同也是关于粒子相对距离的权函数。ij为随机变量，代表两粒子间产生的随机扰动。和满足关系式22Bk T=，其中Bk为 Boltzmann 常数，T 为系统的温度，这种关系耦合了耗散力和随机力。本文采用开源仿真软件 LAMMPS 实现，仿真流程图如图 1 所示。Figure 1.Simulation flowchat 图图 1.仿真流程图 2.2.模型建立模型建立 碰撞模型如图 2 所示，其中*02*eeDD=+，偏心长度 e 顺时针方向偏移为正，逆时针方向偏移为负，为液滴中心与细杆中心的角度，仿真中液滴粒子相关参数的设置在表 1 中列出。Table 1.Droplet particle related parameters 表表 1.液滴粒子相关参数 参数 符号 取值(MDPD 单位)吸引力参数(液体与液体)llA 80 吸引力参数(液体与固体)slA 49 刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4216 建模与仿真 Continued 排斥力参数 llslBB=25.0 时间步长 t 0.01 系统温度 BK T 1.0 随机力系数 4.5 耗散力系数 18.0 截断半径 rc 1.0 排斥力截断半径 rd 0.75 Figure 2.System model diagram 图图 2.系统模型示意图 3.结果与分析结果与分析 液滴以一定角度碰撞细杆时，由于施加的碰撞速度和重力的方向不一致，会导致碰撞状态有所不同。图 3 展示了液滴以不同角度碰撞细杆的形态变化，左侧为时刻t T，t 表示当前时刻，t 为整个碰撞过程的时间，施加的碰撞速度()5.7 m/sV=，此时液滴中心碰撞细杆会分裂为两团小液滴滴落，方便直观地看出液滴的形态变化。由图可知，液滴在碰撞过程的形态变化与中心碰撞相类似；与中心碰撞最直观的区别是，在拉伸断裂之前，两段薄片的扩散方向有所不同，两段薄片所包含的液滴体积相近似，且不会随着碰撞角度的增大而变化；0.84t T=时，当液滴以角度 =30、60碰撞细杆，两团小液滴在扩散的过程中汇聚成了一团液滴，可见最终状态会受到碰撞角度的影响。图 4 显示了不同时刻，液滴以 =0、60碰撞细杆的速度矢量分布图，以此了解影响碰撞状态的原因，右侧热力图通过速度magv与初始碰撞速度 V 的比值进行归一化。当液滴以 =0碰撞细杆，即中心碰撞时，两段薄片外部的速度较大，而靠近细杆的部分速度较小，此外左右两段薄片的速度分布也较为对称，图中0.4t T=，此时与细杆完全接触的液滴，受到细杆表面向上力的影响，液滴速度降低，而靠近细杆的两薄片则出现了向内收缩的趋势；接着在0.7t T=两段薄片继续向内收缩，从而液滴汇聚。当液滴以 =60碰撞细杆时，液滴整体的速度分布仍是靠近细杆的较小而远离细杆的较大，但此时液滴的速度分布不再对称，当0.4t T=，一段薄片由于受到重力的影响，速度略大于另一段，由速度矢量图可以看出，左右两段薄片的速度不再是镜像对称的关系，速度较大的部分向下滴落的趋势也较明显；当0.7t T=，偏移一侧的薄片向下的移动更明显，而另一段薄片仍以斜向下的方向扩散，接着两段薄片汇聚形成一团液滴。由此可以看出不同角度碰撞改变了分裂的两段薄片扩散方向，而导致了结果的不同。刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4217 建模与仿真 Figure 3.Morphological Changes of dropletat =0,30 and 60 图图 3.液滴以 =0，30和 60碰撞细杆的形态变化 Figure 4.Velocity vector distribution when droplet impact on fiber with(a)=0,(b)=60 at V=5.7(m/s)图图 4.V=5.7(m/s)时，液滴以(a)=0，(b)=60碰撞细杆时的速度矢量分布 液滴以不同角度碰撞的分裂滴落临界速度asV，如图 5 所示，进行了多次模拟，得到了不同角度下液滴分裂滴落临界速度的最大和最小误差。由图可以看出碰撞角度与临界速度的变化的规律如下：=10至 =40，随着角度的增加，临界速度随之增大，在 =45左右时，临界速度较大，表明此时扩散阶段的两段薄片比较容易在分裂后形成一团液滴；之后临界速度随 的增大而减小。刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4218 建模与仿真 Figure 5.Critical velocity of droplet splitting atdifferent angles 图图 5.不同角度碰撞下液滴分裂滴落的临界速度 碰撞的最终状态可分为：未接触，回弹，浸润滴落和分裂滴落，由于受到重力影响，当碰撞角度 大于 50时，若碰撞速度较小，液滴会垂直向下移动，此时则会出现液滴没有与细杆接触的情况，图 6汇总了不同碰撞速度和碰撞角度下液滴的最终状态，图中的线段是根据拟合得出的不同状态的临界曲线，其中：粉色曲线为未接触与液滴回弹状态的临界曲线，蓝色曲线为回弹和浸润滴落状态的临界曲线，绿色曲线为浸润滴落和分裂滴落状态的临界曲线，碰撞结果的变化规律可以归结为以下几种情况：当 V 和 都较小时，液滴在接触细杆后慢慢回弹并离开细杆表面；随着 V 和 增大，液滴接触并完全浸润细杆表面后滴落；继续增大 V，液滴接触细杆后沿着接触线的方向继续向外扩散，分裂成两团次级液滴并滴落，当碰撞速度 V=6(m/s)时，=40至 50的碰撞结果为浸润滴落，结合图 5 可知当 位于 45左右时，液滴需要较大的碰撞速度才能分裂滴落。Figure 6.The state of droplet at different angles and velocities 图图 6.不同碰撞角度和速度时，液滴的最终状态 刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4219 建模与仿真 图 7 显示了液滴以不同角度偏心碰撞细杆的形态变化。随着偏心率 e*的增大，两侧的薄片所含液滴比例不同，因此在随后的扩散阶段不呈对称的关系；当偏心率 e*较大时，在浸润阶段(如图中 =0，e*=0.7)，大部分液滴都聚集在偏移一侧，另一侧的液滴较少，所以在随后的扩散阶段不会形成薄片；比较图中 =30，e*=0.3 与 =60，e*=0.5 可以看出液滴与细杆接触阶段(0.4t T=)，两者形态类似，在随后的扩散阶段不同角度的偏心碰撞会引起两段薄片的形态有所不同。Figure 7.Morphological changes of droplet eccentrically impacting on fiber at different angles 图图 7.液滴以不同角度偏心碰撞细杆的形态变化 不同碰撞角度和偏心率左右两段薄片液滴比值定义为()0ee=，其中0为液滴初始所含颗粒数，e为偏移一侧薄片所含颗粒数。图 8 显示了在扩散阶段，图中的虚线表示 e*=10%和 e*=50%时，所对应的值，通过与 e*为负值时对比，随着 e*的增大，液滴比值 随之递减，大致上呈现如图中二次曲线的规律变化，变化如图中黑色曲线；此外当 e*较小时，碰撞角度对 的影响也较明显，随着 e*增加，不同碰撞角度下的 之间的差值减小。逆时针偏移时即 e*为负值，比顺时针偏移时较小，因为此时一侧薄片受到重力的影响更明显，在接触细杆之前偏移一侧包含更多液滴，进而导致在扩散阶段两段薄片更加不对称。与中心碰撞不同，由于偏移一侧薄片所含液滴较多，牵引另一侧的薄片，因此当偏心率较大时、速度较小时如 V=4.2(m/s)，两侧的薄片有短暂的汇聚，但两者的移动方向不同导致最终在随后的拉伸断裂阶段形成两团小液滴，如图 9(a)所示的速度矢量分布和热力图。当偏心率更大时，如图 9(b)所示，此时 更小，即偏移另一侧所含液滴更少，液滴接触杆后，液滴掠过细杆表面，不会形成拉伸的两段薄片，仍然以初始的液滴状态向着斜下方移动。当碰撞角度更大且 e*为负值时，若继续减小碰撞速度 V，则会出现液滴没有与细杆接触的情况。刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4220 建模与仿真 Figure 8.Droplet ratio of two lamellae under different impacting angles and eccentricity 图图 8.不同碰撞角度和偏心率下两段薄片所含液滴比例 Figure 9.Velocity vector distribution when droplet impact on fiber with(a)e*=30%,(b)e*=80%at V=4.2(m/s),=60 图图 9.V=4.2(m/s)，=60时，液滴以(a)e*=30%，(b)e*=80%碰撞细杆的速度矢量分布 图 10 显示了不同偏心率和碰撞角度下，液滴与细杆接触的临界速度ecV和分裂滴落的临界速度esV。如图 10(a)所示，当*30%e(或*30%e (或*70%e )，偏移一侧所含的液滴较少，因此在扩散阶段不会形成薄片，可见一定的 e*范围之间，能促进液滴分裂滴落。此外，对比 =0、60和 90，临界速度esV的变化规律类似，当 e*相同时esV之间的差值较为接近，因此不同的角度对于esV的影响较小。当 50时，才会出现刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4221 建模与仿真 液滴与细杆没有接触的情况，如图 10(b)所示，随着 e*的增加，ecV递减；同时，当碰撞角度 增加ecV也随之增大。(a)(b)Figure 10.Critical velocity Ves of droplet splitting and critical velocity Vec of droplet contact with fiber at different eccentric-ity and angles 图图 10.不同偏心率和碰撞角度下，液滴分裂滴落的临界速度 Ves(e*=70%70%)和接触细杆的临界速度 Vec 4.总结总结 本文采用 MDPD 方法模拟了液滴以一定角度偏心碰撞圆柱形细杆，单独考虑了碰撞角度对于结果影响，在重力的影响下，液滴以抛物线的轨迹向下移动，扩散阶段中分裂的两段薄片运动方向不同，造成最终的状态发生改变。在一定角度范围内，需要更大的初始速度，液滴才能分裂滴落；此外，当碰撞角度较大时，若碰撞速度较小，液滴不会接触细杆；接着综合考虑了碰撞角度和偏心率对于液滴两侧薄片比例、液滴分裂滴落临界esV和接触细杆临界速度ecV的影响。结果表明，分裂滴落时的两段薄片所含液滴会随着偏心率改变，而碰撞角度的影响则较小；分裂滴落的临界速度随着偏心率增大先减小而后增大，不同角度下，临界速度的改变趋势类似，接触细杆的临界速度随着偏心率增大而减小。参考文献参考文献 1 Banitabaei,S.A.and Amirfazli,A.(2017)Droplet Impact onto a Solid Sphere:Effect of Wettability and Impact Veloc-ity.Physics of Fluids,29,Article ID:062111.https:/doi.org/10.1063/1.4990088 2 Zhang,H.X.,Yi,X.,Du,Y.X.,et.al.(2019)Dynamic Behavior of Water Drops Impacting on Cylindrical Superhydro-phobic Surfaces.Physics of Fluids,31,Article ID:032104.https:/doi.org/10.1063/1.5083070 3 Khurana,G.,Sahoo,N.and Dhar,P.(2019)Post-Collision Hydrodynamics of Droplets on Cylindrical Bodies of Va-riant Convexity and Wettability.Physics of Fluids,31,Article ID:022008.https:/doi.org/10.1063/1.5064799 4 Liu,X.H.,Wang,K.M.,Fang,Y.Q.,et.al.(2019)Study of the Surface Wettability Effect on Dynamic Characteristics of Droplet Impacting a Tube with Different Curvature Ratios.Experimental Thermal and Fluid Science,115,Article ID:110060.https:/doi.org/10.1016/j.expthermflusci.2020.110060 5 Wang,Y.L.,Wang,Y.X.and Wang,S.R.(2020)Droplet Impact on Cylindrical Surfaces:Effects of Surface Wettabil-ity,Initial Impact Velocity and Cylinder Size.Journal of Colloid and Interface Science,578,207-217.https:/doi.org/10.1016/j.jcis.2020.06.004 6 Wang,C.Y.,Xiao,W.,Zhang,H.X.,et al.(2021)A Many-Body Dissipative Particle Dynamics Study of Eccentric Droplets Impacting Inclined Fiber.Physics of Fluids,33,Article ID:042001.https:/doi.org/10.1063/5.0044065 7 Warren,P.B.(2004)Vapor-Liquid Coexistence in Many-Body Dissipative Particle Dynamics.Physical Review E,68,Article ID:066702.https:/doi.org/10.1103/PhysRevE.68.066702 刘华建，曹达敏 DOI:10.12677/mos.2023.124384 4222 建模与仿真 8 Trofimov,S.Y.,Nies,E.L.F.and Michels,M.A.J.(2005)Constant-Pressure Simulations with Dissipative Particle Dy-namics.The Journal of Chemical Physics,123,Article ID:144102.https:/doi.org/10.1063/1.2052667 9 Liu,M.B.,Meakin,P.and Huang,H.(2006)Dissipative Particle Dynamics with Attractive and Repulsive Particle-Particle Interactions.Physics of Fluids,18,Article ID:017101.https:/doi.org/10.1063/1.2163366 10 Chen,S.,Liu,Y.,Khoo,B.C.,et.al.(2007)Mesoscopic Simulation of Binary Immiscible Fluids Flow in a Square Mi-crochannel with Hydrophobic Surfaces.Computer Modelingin Engineering and Sciences,19,181-196.