剪切流中乳液碰撞动力学行为的多维度解析与探究

剪切流中乳液碰撞动力学行为的多维度解析与探究一、绪论1.1研究背景与意义乳液作为一种将一种液相以液滴形式分散到另一种互不相溶的液相中形成的混合物，在能源、化工、生医、材料等前沿科技领域有着极为广泛的应用。在能源领域，乳液可用于提高石油采收率，通过调节乳液的性质，使其能够更好地驱替油藏中的原油，从而提高能源的开采效率；在化工领域，乳液常用于涂料、油墨、胶粘剂等产品的制备，例如水性涂料中的乳液能够使颜料均匀分散，提高涂料的稳定性和性能；在生医领域，乳液可作为药物载体，实现药物的靶向输送和控制释放，像一些抗癌药物通过乳液载体能够更精准地作用于肿瘤细胞，减少对正常细胞的伤害；在材料领域，乳液可用于制备纳米材料、多孔材料等，为材料科学的发展提供了新的途径。乳液微流控技术作为一种借助微流控过程中互不相溶流体间相互作用来实现微乳液液滴制备及精确操控的技术，在微混合、微反应、生物封装、功能材料制备、药物运输与释放等领域中极具发展前景，已成为国际上微流控技术领域发展的新热点。在微混合领域，乳液微流控技术能够实现微尺度下的高效混合，为化学反应提供更均匀的反应环境；在微反应领域，可通过精确控制乳液液滴的生成和反应条件，实现对反应过程的精细调控，提高反应的选择性和产率；在生物封装领域，能够将生物分子、细胞等封装在乳液液滴中，为生物研究和生物医学应用提供了有力的工具；在功能材料制备领域，可制备出具有特殊结构和性能的功能材料，满足不同领域的需求；在药物运输与释放领域，能够实现药物的精准运输和可控释放，提高药物的疗效和安全性。乳液微流控过程中涉及一系列复杂的界面现象及流体动力学问题，如乳液的形变、破碎及其间的相互作用（聚并、碰撞等）、多相离散与汇聚及其多界面演化等，这些问题直接影响到乳液微流控过程的操控精度与效率，进而决定了工程实际中乳液产品的外部形貌和内部微结构等重要品质指标。例如，在乳液液滴的制备过程中，液滴的尺寸分布和稳定性会受到乳液间相互作用的影响，如果液滴尺寸不均匀，可能会导致产品性能的不一致；在药物运输与释放过程中，乳液的形变和破碎可能会影响药物的释放速率和靶向性。因此，乳液体系的多相界面演化及流体动力学行为机理已成为乳液微流控技术基础研究的前沿热点。液滴碰撞是液滴群乳液体系中液滴间相互作用的基本形式，特别是在液滴密集度高的情况下，液滴碰撞行为直接影响乳液液滴在流场中的形态演变。在实际的乳液体系中，液滴之间的碰撞是不可避免的，而碰撞的结果可能导致液滴的聚并、破碎或者改变液滴的运动轨迹和形态。充分认识液滴碰撞动力学行为，不仅对多相流体动力学具有重要的学术价值，能够丰富和完善多相流体动力学的理论体系，加深对多相流体复杂流动现象的理解；而且对乳液液滴精密调控也具有重要的现实意义，有助于实现乳液液滴的精准制备和操控，提高乳液产品的质量和性能，满足不同领域对乳液的特殊需求，推动相关产业的发展。1.2国内外研究现状目前，国际上对单个液滴的流变行为开展了较为深入的研究。20世纪30年代，Taylor等人便对单个液滴在流场中变形与破裂行为开展了实验研究，为后续的研究奠定了基础。随后，Rallison、Youngren等学者对单一液滴的形变行为进行了数学建模和数值模拟，从理论层面进一步揭示了液滴形变的规律。Stone对单液滴在流场中流变行为的实验与理论工作进行了全面总结，综述了液滴界面动力学方面的研究进展，着重介绍了流体物性参数（如黏度、密度、界面张力等）和不同的流场类型（如剪切流场、拉伸流场、旋转流场等）等因素对液滴形变和液滴破碎的影响，为该领域的研究提供了系统的理论框架。最近，Salkin对液滴在流场中遇到障碍物的破碎行为进行了实验研究，探讨了障碍物对液滴破碎的影响机制；Ulloa实验研究了剪切流场的受限程度对液滴形变与破碎的影响，发现受限程度会改变液滴所受的作用力，从而影响其形变和破碎行为。刘向东等采用VOF方法对剪切流场中的单个双乳液液滴进行了液滴形变及破碎特性的研究，结果表明液滴的稳定形变随着毛细数的增大而加剧，并且在给定的雷诺数下存在一个临界毛细数，小于临界毛细数时，液滴只发生稳定变形，反之液滴破碎。相比于单个液滴在外流场中的流变行为，两个或者更多液滴在外流场的碰撞过程中存在液滴间的相互作用，使得其流变机理变得更加复杂。目前，已有少部分学者在此领域开展了相关的实验研究，针对液滴在剪切流场中碰撞过程的运动轨迹进行了可视化观测，但受到量测手段的限制，对液滴在剪切流场中相互碰撞过程的形变与碰撞的耦合机理以及压力变化特性尚未充分揭示。王程遥等基于VOF液/液相界面追踪方法，建立了不可压缩水/油单乳液液滴动力学模型并进行数值求解，模拟研究了在剪切流场作用下两个相同体积的球形液滴在碰撞过程中的变形、相互作用及压力演化，分析了液滴碰撞过程中毛细数对液滴间相互作用的影响，发现随着毛细数增大，界面挤压变形越明显，液滴变形系数也增大。在双乳液液滴动力学行为研究方面，其形变动力学行为和破裂动力学行为都受到了一定关注。Ha研究了双重乳液液滴在电场作用下的动力学特性，发现电场会影响液滴的形变和运动；刘向东等采用基于Cahn-Hilliard方程的相场方法捕捉液-液流体界面，并利用电流模型表征电场力作用，建立了电场和拉伸流场共同作用下双重乳液形变动力学模型，研究揭示了双重乳液液滴在剪切流场下的稳定特性。然而，对于双乳液液滴在复杂流场中，尤其是涉及碰撞过程时的动力学行为研究还相对较少，其内部多界面的相互作用机制以及内液滴对整个液滴碰撞行为的影响等方面仍有待深入探究。综合来看，当前对于乳液液滴动力学行为的研究取得了一定成果，但在液滴碰撞动力学行为方面，尤其是在复杂剪切流场中，不同类型乳液液滴（单乳液、双乳液等）碰撞过程的全面深入研究还存在不足。对于液滴碰撞过程中多物理场（如流场、电场、磁场等）耦合作用下的动力学行为研究还不够系统，缺乏对液滴碰撞过程中能量转化和传递机制的深入分析。此外，实验研究中对液滴碰撞过程的高精度测量和可视化技术仍有待进一步发展，数值模拟中对于复杂界面和多相流的精确模拟方法也需要不断完善，以更准确地揭示乳液液滴在剪切流中的碰撞动力学行为。1.3研究内容与方法本研究将综合运用实验观测、理论分析和数值模拟的方法，深入探究剪切流中乳液液滴碰撞过程的动力学行为。在理论分析方面，基于流体力学基本原理，结合液滴在剪切流场中的受力情况，建立适用于描述单乳液液滴和双乳液液滴在剪切流中碰撞过程的数学模型。对于单乳液液滴，考虑其在剪切流场中受到的黏性力、表面张力以及惯性力，通过建立动量方程、连续性方程，并引入合适的边界条件和无量纲参数（如雷诺数Re、毛细数Ca、黏度比A和密度比B等）来表征各相参数及各种作用力间的相对大小，从而构建起单乳液液滴碰撞的数学模型。对于双乳液液滴，除了考虑上述因素外，还需考虑内液滴与外液滴之间的相互作用，包括界面张力、压力差等对液滴整体行为的影响，建立更为复杂的数学模型来准确描述其在剪切流中的碰撞过程。同时，利用张量分析、摄动理论等数学工具对模型进行求解和分析，推导液滴在碰撞过程中的运动轨迹、变形程度、压力分布等关键物理量的理论表达式，为深入理解液滴碰撞动力学行为提供理论基础。在实验观测方面，设计并搭建一套高精度的剪切流场中单（双）乳液液滴碰撞的高速可视化实验平台。该平台主要包括剪切流场发生装置、乳液液滴注入系统、高速摄像及图像采集分析系统等部分。剪切流场发生装置采用双板平行剪切或旋转式剪切等方式，能够产生稳定且可控的剪切流场，通过调节两板的相对运动速度或旋转部件的转速来精确控制剪切率。乳液液滴注入系统能够精确控制单乳液液滴和双乳液液滴的生成和注入位置，确保液滴以特定的初始条件进入剪切流场。高速摄像及图像采集分析系统配备高分辨率、高帧率的高速摄像机，能够对液滴在剪切流场中的碰撞过程进行实时捕捉，帧率可达每秒数千帧甚至更高，以获取液滴碰撞瞬间及后续演化过程的详细图像信息。利用数字图像处理技术对采集到的图像进行分析，提取液滴的形状、尺寸、运动轨迹等关键信息，通过对不同实验条件下的大量实验数据进行统计和分析，总结液滴碰撞过程中的动力学行为规律，验证理论模型的准确性，并为数值模拟提供实验依据。在数值模拟方面，基于VOF（VolumeofFluid）相界面追踪方法，对建立的单（双）乳液液滴碰撞数学模型进行数值求解。利用计算流体力学软件（如ANSYSFluent、OpenFOAM等）对剪切流场进行数值模拟，通过合理划分计算网格，采用合适的数值算法（如有限体积法、有限差分法等）对控制方程进行离散求解，精确模拟液滴在剪切流场中的碰撞过程。在模拟过程中，详细分析液滴的变形、合并、破碎等现象，以及流场中的速度分布、压力分布等参数的变化情况。通过改变模拟参数（如剪切率、液滴尺寸、流体物性参数等），系统研究各因素对液滴碰撞动力学行为的影响规律，与理论分析和实验结果进行对比验证，进一步完善和优化数学模型。同时，利用数值模拟结果对液滴碰撞过程进行可视化展示，直观呈现液滴在剪切流中的复杂运动和相互作用过程，为深入理解液滴碰撞动力学行为提供直观的参考。具体研究内容如下：建立剪切流中乳液液滴碰撞的数学模型：根据流体力学基本方程，结合液滴在剪切流场中的受力分析，分别建立单乳液液滴和双乳液液滴碰撞的数学模型。考虑液滴间的相互作用、界面张力、黏性力、惯性力等因素，引入合适的边界条件和无量纲参数，如雷诺数（Re）、毛细数（Ca）、黏度比（A）和密度比（B）等，以准确描述液滴在剪切流中的碰撞过程。搭建剪切流场中乳液液滴碰撞的实验平台：设计并搭建包括剪切流场发生装置、乳液液滴注入系统、高速摄像及图像采集分析系统在内的实验平台。利用该平台进行单（双）乳液液滴在剪切流场中碰撞的实验研究，通过高速摄像记录液滴碰撞过程，运用数字图像处理技术提取液滴的形状、尺寸、运动轨迹等信息，分析液滴碰撞过程中的动力学行为。研究剪切流中乳液液滴碰撞的动力学行为：通过数值模拟和实验研究，深入分析单（双）乳液液滴在剪切流场中碰撞过程的动力学行为，包括液滴的变形、合并、破碎等现象，以及流场中的速度分布、压力分布等参数的变化情况。探讨不同因素（如剪切率、液滴尺寸、流体物性参数等）对液滴碰撞动力学行为的影响规律，揭示液滴碰撞过程中的内在物理机制。分析影响剪切流中乳液液滴碰撞的因素：系统研究剪切率、液滴尺寸、流体物性参数（如黏度、密度、界面张力等）、液滴初始位置和角度等因素对液滴碰撞动力学行为的影响。通过改变这些因素进行数值模拟和实验研究，分析各因素对液滴碰撞概率、碰撞后液滴的形态和运动状态等方面的影响，建立影响因素与液滴碰撞动力学行为之间的定量关系。二、剪切流中乳液碰撞动力学行为的理论基础2.1乳液与剪切流的基本概念乳液是一种多相体系，它是由一种液相以微小液滴的形式分散在另一种与之互不相溶的液相中所形成的。其中，被分散的液相称为分散相或内相，而连续的液相则称为连续相或外相。根据分散相和连续相的不同，乳液主要分为水包油（O/W）型和油包水（W/O）型两种类型。在水包油型乳液中，油相以液滴的形式分散在水相中，例如常见的牛奶，就是一种典型的水包油型乳液，其中油脂小液滴分散在水中；而在油包水型乳液中，水相则分散在油相中，如一些护肤霜，水相以微小液滴的形式均匀分散于油相基质中。乳液的稳定性是其重要特性之一，它受到多种因素的影响，包括界面张力、乳化剂的种类和浓度、液滴的大小和分布、以及连续相和分散相的性质等。界面张力越小，乳液越容易形成且稳定性越高；乳化剂能够降低界面张力，增加乳液的稳定性；较小且均匀的液滴尺寸分布也有助于提高乳液的稳定性。从微观结构来看，乳液中的液滴并非孤立存在，它们之间存在着复杂的相互作用。液滴表面通常带有电荷，这些电荷会在液滴周围形成双电层，双电层之间的静电斥力能够阻止液滴相互靠近并聚并，从而维持乳液的稳定性。此外，乳液中还可能存在其他添加剂，如增稠剂、防腐剂等，它们也会对乳液的微观结构和稳定性产生影响。增稠剂可以增加连续相的黏度，减小液滴的沉降速度，进一步提高乳液的稳定性；防腐剂则能防止微生物的生长繁殖，延长乳液的保质期。剪切流是一种流体流动形式，其特点是流体在流动过程中存在速度梯度。当流体在两个平行平板之间流动时，如果其中一个平板静止，另一个平板以恒定速度移动，那么平板间的流体就会形成剪切流。在这种情况下，靠近运动平板的流体速度较快，而靠近静止平板的流体速度较慢，从而在垂直于平板方向上产生速度梯度。剪切流的速度分布通常呈现出线性或非线性的变化规律，这取决于流体的性质和流动条件。对于牛顿流体，在简单的平行平板剪切流中，速度分布是线性的；而对于非牛顿流体，速度分布可能是非线性的，例如假塑性流体在剪切流中，随着剪切速率的增加，其黏度会降低，速度分布会变得更加复杂。常见的剪切流类型包括简单剪切流、Couette流和Poiseuille流等。简单剪切流是指流体在两个平行平板之间，由于平板的相对运动而产生的剪切流动，这种流动在许多微流控实验和理论研究中被广泛应用，是研究剪切流中乳液行为的基础模型之一；Couette流是一种特殊的简单剪切流，它通常由一个旋转的内筒和一个静止的外筒组成，流体在两筒之间的环形空间内流动，这种流动形式在研究流体的黏性、流动性以及液滴在旋转环境中的行为等方面具有重要应用；Poiseuille流则是指流体在圆形管道中，由于压力差而产生的剪切流动，在实际的管道输送、生物体内的血液流动等场景中较为常见。当剪切流作用于乳液体系时，会对乳液产生多方面的影响。剪切流会使乳液中的液滴受到剪切力的作用，从而发生形变和运动。在剪切力的作用下，液滴会被拉伸、扭曲，其形状会偏离球形，且液滴的运动轨迹也会发生改变，不再是简单的直线运动，而是会随着剪切流的方向和速度梯度进行复杂的曲线运动。剪切流还可能导致乳液液滴之间的碰撞频率增加。由于液滴在剪切流中的运动速度和方向各不相同，不同液滴之间更容易相互靠近并发生碰撞，这对乳液的稳定性和微观结构演化具有重要影响。如果液滴碰撞后发生聚并，乳液的液滴尺寸分布会发生变化，可能导致乳液的性质发生改变；而如果液滴碰撞后未发生聚并，液滴之间的相互作用也会影响乳液的流变性质。此外，剪切流还可能影响乳液中乳化剂的分布和吸附状态。乳化剂在液滴表面的吸附平衡可能会被打破，导致乳化剂重新分布，进而影响乳液的界面性质和稳定性。2.2碰撞动力学模型概述在研究乳液液滴在剪切流中的碰撞动力学行为时，多种撞击理论模型为我们提供了重要的分析基础，其中经典力学、弹性应力波传播、接触力学和塑性变形理论各自从不同角度描述了撞击过程。经典力学在碰撞研究中，主要应用基本力学定理来预测撞击后的速度，其核心是脉冲－动量定理。在两球碰撞的简单模型中，根据动量守恒定律m_1v_{1i}+m_2v_{2i}=m_1v_{1f}+m_2v_{2f}（其中m_1、m_2为两球质量，v_{1i}、v_{2i}为碰撞前速度，v_{1f}、v_{2f}为碰撞后速度），可以计算出碰撞后的速度状态。实际问题中，由于碰撞过程存在能量损失，这一损失通过恢复系数e来体现，恢复系数定义为碰撞后两物体分离速度与碰撞前接近速度的比值，即e=\frac{v_{2f}-v_{1f}}{v_{1i}-v_{2i}}。对于完全弹性碰撞，e=1，碰撞过程无能量损失；对于完全非弹性碰撞，e=0，两物体碰撞后粘在一起，能量损失最大；而一般的碰撞，0<e<1。经典力学方法具有简便和易于实现的特点，但它无法预报物体之间的接触力和物体的应力，对于乳液液滴碰撞中涉及的复杂界面力和内部应力分布难以准确描述。弹性应力波传播理论认为，撞击是通过以撞击点为起点，应力波在撞击物之间的传播来描述的。当两个乳液液滴发生碰撞时，碰撞瞬间产生的应力波会在液滴内部传播，总能量中的一部分会转化为振动能。在研究两杆的纵向碰撞时，应力波在杆中传播，其传播速度与杆的材料性质（如弹性模量、密度等）有关。这种理论考虑了碰撞过程中的能量转化和波动传播，弥补了经典力学在能量转化细节描述上的不足，然而它对于复杂的多液滴碰撞体系以及液滴内部复杂的流场和界面相互作用的处理仍存在一定困难。接触力学主要关注两个物体撞击产生的接触应力。对于球形接触面，Hertz理论常被用于获得撞击关系，从而计算撞击时间和最大变形。假设两个半径分别为R_1和R_2的球形液滴发生碰撞，接触力F与接触半径a之间的关系可通过Hertz理论表示为F=\frac{4}{3}E^*\sqrt{R^*}a^{\frac{3}{2}}（其中E^*为等效弹性模量，R^*为等效半径）。该理论在处理简单几何形状物体的碰撞时具有较高的准确性，但实际乳液液滴的形状并非完全规则的球形，且碰撞过程中液滴会发生形变，这使得Hertz理论的直接应用受到限制。为了更准确地描述碰撞过程，通常会假设材料有一个屈服点，当Hertz理论不适用时，可使用屈服区模型。此外，撞击力-变形关系常通过增加一个阻尼项来反映接触区域的能量耗散，从而将接触区作为一个弹簧-阻尼系统的模型，以更好地模拟实际碰撞过程中的能量损失和动态响应。塑性变形理论主要应用于塑性应变超过容许变形的高速撞击问题，如爆炸与侵彻等场景。在乳液液滴碰撞研究中，虽然通常不属于高速撞击范畴，但当液滴受到较大的外力作用，导致其内部发生较大的塑性变形时，也需要考虑塑性变形理论。Goldsmit提供了水动力学理论和塑性波传播理论两种方法来处理这类问题。水动力学理论中，假设物体密度发生变化，材料的状态方程与密度、温度的变化相关，同时利用能量、动量守恒等原理来分析碰撞过程；塑性波传播理论则侧重于研究塑性变形在材料中的传播规律，以及塑性变形对物体力学性能和碰撞行为的影响。然而，在乳液体系中，由于液滴的流动性和复杂的界面特性，将塑性变形理论直接应用于乳液液滴碰撞存在一定的复杂性，需要进一步结合乳液的特点进行修正和完善。恢复系数在碰撞动力学研究中具有重要意义，它不仅反映了碰撞过程的能量损失程度，还与多种因素密切相关。材料性质是影响恢复系数的关键因素之一，不同材料的弹性模量、屈服强度等力学性能不同，导致其在碰撞过程中的能量耗散方式和程度也不同，从而影响恢复系数的大小。对于硬度较高、弹性较好的材料制成的物体，碰撞时能量损失相对较小，恢复系数较大；而对于较软、易发生塑性变形的材料，碰撞时能量损失较大，恢复系数较小。碰撞速度也对恢复系数有显著影响，一般来说，随着碰撞速度的增加，物体在碰撞过程中发生的变形更加剧烈，能量损失增加，恢复系数会减小。碰撞角度同样会影响恢复系数，当两个物体以不同的角度碰撞时，接触点的受力情况和能量传递方式会发生变化，进而导致恢复系数的改变。在乳液液滴碰撞研究中，恢复系数可以帮助我们判断液滴碰撞后的运动状态和能量变化，通过实验测量和理论分析不同条件下乳液液滴碰撞的恢复系数，能够为建立准确的碰撞动力学模型提供重要依据，深入理解乳液液滴在剪切流中的碰撞行为。2.3相关物理参数及无量纲数在研究剪切流中乳液碰撞动力学行为时，雷诺数（Reynoldsnumber，Re）是一个重要的无量纲数，它在表征乳液碰撞过程中惯性力与黏性力的相对大小方面发挥着关键作用。雷诺数的定义为流体的惯性力与黏性力之比，其数学表达式为Re=\frac{\rhovL}{\mu}，其中\rho表示流体的密度，v是流体的特征速度，L为特征长度，\mu为流体的动力黏度。当雷诺数较小时，表明黏性力在流体运动中占据主导地位，流体的流动呈现出较为规则的层流状态，液滴在这种流场中的运动相对稳定，受到的黏性阻力较大，其运动轨迹和形变主要受黏性力的制约。例如，在一些微流控实验中，当雷诺数较低时，乳液液滴在剪切流中的运动较为平稳，液滴之间的相互作用主要通过黏性力传递，液滴的碰撞过程相对缓和，不易发生剧烈的变形和破碎。当雷诺数较大时，惯性力则成为主导因素，流体的流动转变为湍流状态，此时流场中的速度分布变得复杂，存在大量的漩涡和脉动，液滴在这种环境下会受到更强的惯性力作用，其运动变得更加不规则，碰撞时更容易发生剧烈的变形和破碎。在工业搅拌过程中，随着搅拌速度的增加，雷诺数增大，乳液液滴在强惯性力的作用下，碰撞频率增加，碰撞后的变形和破碎程度也更为显著，这对乳液的稳定性和微观结构产生重要影响。毛细数（Capillarynumber，Ca）也是一个关键的无量纲数，它用于衡量作用在液滴上的黏性力与表面张力的相对大小。毛细数的定义为Ca=\frac{\muv}{\sigma}，其中\sigma为界面张力。在乳液体系中，表面张力倾向于使液滴保持球形，以最小化液-液界面面积，从而降低体系的表面自由能；而黏性力则在剪切流的作用下，试图改变液滴的形状。当毛细数较小时，表面张力占据主导地位，液滴能够较好地保持球形，抵抗变形的能力较强，在碰撞过程中，液滴的变形程度较小，不易发生破碎。在低剪切速率的乳液体系中，液滴的毛细数较小，液滴之间的碰撞主要表现为弹性碰撞，液滴的形状基本保持不变。当毛细数较大时，黏性力超过表面张力，液滴更容易发生变形和破裂，在碰撞过程中，液滴会被显著拉伸、扭曲，甚至破碎成多个小液滴。在高剪切速率的情况下，乳液液滴受到较大的黏性力作用，毛细数增大，液滴在碰撞时容易发生破裂，导致乳液的粒径分布发生变化，进而影响乳液的性质和稳定性。黏度比（Viscosityratio，A）是分散相黏度\mu_d与连续相黏度\mu_c的比值，即A=\frac{\mu_d}{\mu_c}，它对乳液液滴在碰撞过程中的变形和运动行为有着重要影响。当黏度比较小时，意味着分散相的黏度相对连续相较低，液滴在剪切流中更容易受到连续相的作用而发生变形，在碰撞过程中，液滴的变形程度较大，且可能更容易与其他液滴发生聚并。在一些水包油型乳液中，如果油相（分散相）的黏度较低，在剪切流和碰撞作用下，油滴更容易被拉伸变形，与周围油滴聚并的可能性增加，从而影响乳液的稳定性。相反，当黏度比较大时，分散相黏度较高，液滴的刚性相对较大，抵抗变形的能力增强，在碰撞过程中，液滴的变形相对较小，运动轨迹也相对稳定，聚并的倾向相对较低。如果乳液中的分散相是高黏度的聚合物溶液，在碰撞时，液滴能够较好地保持自身形状，不易与其他液滴发生聚并，使得乳液的稳定性得以维持。密度比（Densityratio，B）是分散相密度\rho_d与连续相密度\rho_c的比值，即B=\frac{\rho_d}{\rho_c}，它在乳液碰撞过程中主要影响液滴的沉降和上浮行为，进而影响液滴的碰撞概率和碰撞后的运动状态。当密度比接近1时，分散相和连续相的密度相近，液滴在剪切流中的沉降或上浮速度较小，液滴在体系中分布较为均匀，碰撞的随机性相对较大，碰撞后的运动轨迹受密度因素的影响较小。在一些特殊的乳液体系中，通过调整分散相和连续相的成分，使密度比接近1，液滴在体系中能够长时间稳定悬浮，碰撞过程主要受剪切流和其他作用力的影响。当密度比与1相差较大时，液滴会由于密度差异而产生明显的沉降或上浮运动，这会改变液滴在体系中的分布情况，使得液滴在某些区域更容易发生碰撞，碰撞后的运动方向也会受到密度差导致的浮力或重力的影响。在油包水型乳液中，如果水相（分散相）的密度大于油相（连续相），水液滴会在重力作用下向下沉降，在沉降过程中与其他液滴发生碰撞的概率增加，且碰撞后液滴的运动方向可能会受到重力和剪切流的共同作用而发生改变。三、剪切流场中单乳液液滴碰撞的数值模拟3.1数学模型建立3.1.1物理模型描述考虑在双板平行剪切流场中，两个半径均为a的球形单乳液液滴在连续相中运动并发生碰撞的过程。连续相为牛顿流体，分散相（液滴）也视为牛顿流体，两液滴初始时刻在x-y平面内关于中心点对称分布，流动方向为x方向，速度梯度方向为y方向，x方向与y方向的初始质心距分别为\Deltax_0和\Deltay_0。两块平板在两个平行平面上分别以相同速度U向相反的方向滑动，从而产生剪切率G=\frac{2U}{H}（H为两板间距）的恒定剪切流场作用于液滴，促使液滴产生运动和变形。在该物理模型中，液滴与连续相之间存在明显的界面，界面处的张力对液滴的变形和运动起着关键作用。同时，液滴在剪切流场中受到黏性力和惯性力的作用，这些力的相互作用决定了液滴的碰撞动力学行为。例如，当剪切率G增大时，液滴受到的黏性力增大，液滴更容易发生变形和碰撞；而当液滴的初始位置发生变化时，液滴的碰撞轨迹和碰撞时刻也会相应改变。3.1.2控制方程连续性方程：在整个流场区域内，满足质量守恒定律，其连续性方程为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho为流体密度，\vec{v}为流体速度矢量，t为时间。该方程表明，在单位时间内，流入和流出某一控制体积的流体质量之差等于该控制体积内流体质量的变化率。在本文研究的不可压缩流体中，\frac{\partial\rho}{\partialt}=0，即流体密度不随时间变化。动量方程：考虑到流体的黏性、惯性以及界面张力的作用，动量方程为：\rho(\frac{\partial\vec{v}}{\partialt}+\vec{v}\cdot\nabla\vec{v})=-\nablap+\nabla\cdot(2\mu\vec{D})+\vec{F}_{vol}其中，p为压力，\mu为动力黏度，\vec{D}为应变率张量，\vec{F}_{vol}为界面张力源项。应变率张量\vec{D}的表达式为\vec{D}=\frac{1}{2}(\nabla\vec{v}+(\nabla\vec{v})^T)，它描述了流体微元的变形速率。界面张力源项\vec{F}_{vol}采用CSF（continuumsurfaceforce）方法计算，即\vec{F}_{vol}=\sigma\kappa\nabla\alpha，其中\sigma为界面张力系数，\kappa为界面曲率，\alpha为体积分数函数（对于连续相\alpha=1，对于分散相\alpha=0）。3.1.3边界条件周期性边界条件：在剪切流场的左右两侧边界采用周期性边界条件，即\vec{v}(x+L_x,y,z,t)=\vec{v}(x,y,z,t)，p(x+L_x,y,z,t)=p(x,y,z,t)，其中L_x为x方向的计算域长度。这意味着在x方向上，流场的物理量在经过一个周期长度L_x后保持不变，这种边界条件可以模拟无限长的剪切流场，避免边界效应的影响。壁面边界条件：对于上下平板壁面，采用无滑移边界条件，即\vec{v}(x,y=0,z,t)=\vec{v}(x,y=H,z,t)=0，表示流体在壁面处的速度为零，与壁面之间不存在相对滑动。同时，压力在壁面处满足法向应力连续条件，即p_{n}(x,y=0,z,t)=p_{n}(x,y=H,z,t)，其中p_{n}为法向压力。3.1.4数值求解方法采用有限体积差分法对控制方程进行离散求解。在求解过程中，压力场和速度场的耦合采用SIMPLE（Semi-ImplicitMethodforPressure-LinkedEquations）算法，该算法通过迭代求解压力修正方程，实现压力场和速度场的相互匹配，保证计算结果的准确性和稳定性。动量方程采用一阶迎风差分格式离散，这种格式能够较好地处理对流项，提高计算的稳定性。基于分段式界面计算（PLIC，Piecewise-LinearInterfaceCalculation）方法的Geo-Reconstruct格式对液/液界面进行重构，PLIC方法通过在每个网格单元内构建分段线性的界面来精确描述液滴的形状和位置，Geo-Reconstruct格式则在此基础上进一步优化了界面重构的精度和效率，能够准确地捕捉液滴在碰撞过程中的复杂界面变形。控制微分方程经离散成为线性代数方程组后，采用Gauss-Seidel方法进行迭代求解。在迭代计算中，为了保证计算的收敛性和稳定性，采用了欠松弛因子，其中压力项欠松弛因子为0.2，密度项为0.5，体积力源项为0.5，动量源项为0.2。通过调整这些欠松弛因子，可以使迭代过程更快地收敛到稳定解。时间步长取0.1ms，当每个迭代时间步的流体速度残差小于10^{-3}时，认为该时间步迭代计算收敛，即此时计算得到的速度场和压力场满足控制方程和边界条件的要求。3.1.5网格灵敏性分析为了确保所得的数值解均为网格独立的解，采用不同的网格尺寸进行网格独立性检测。分别设置网格尺寸为\Deltax=\Deltay=\Deltaz=0.05a、\Deltax=\Deltay=\Deltaz=0.1a和\Deltax=\Deltay=\Deltaz=0.2a（a为液滴半径）进行数值模拟，对比不同网格尺寸下液滴碰撞过程中的关键物理量，如液滴的变形系数、碰撞时间等。以液滴的变形系数D=(L-B)/(L+B)（L和B分别为发生变形后的椭圆形液滴的长轴与短轴）为例，当网格尺寸从\Deltax=\Deltay=\Deltaz=0.2a减小到\Deltax=\Deltay=\Deltaz=0.1a时，变形系数的计算结果变化较小，相对误差在5\%以内；当网格尺寸进一步减小到\Deltax=\Deltay=\Deltaz=0.05a时，变形系数的计算结果与\Deltax=\Deltay=\Deltaz=0.1a时相比，相对误差在2\%以内，表明此时网格尺寸对计算结果的影响已经很小，可以认为计算结果与网格无关。因此，在后续的数值模拟中，选择网格尺寸\Deltax=\Deltay=\Deltaz=0.1a，以在保证计算精度的同时，提高计算效率，减少计算资源的消耗。3.1.6模型验证为了验证所建立模型的准确性，对文献中Guido等开展的液滴碰撞实验进行了数值模拟对比。在相同的实验条件下，包括相同的剪切率、液滴初始位置、流体物性参数等，将数值模拟得到的液滴碰撞演化过程与实验结果进行比较。从液滴的运动轨迹来看，数值模拟结果与实验结果基本一致，液滴在剪切流场中的运动趋势、碰撞时刻以及碰撞后的分离方向等都能很好地吻合；从液滴的变形情况来看，数值模拟得到的液滴变形形状和变形程度与实验图像中的液滴变形也具有较高的相似性。通过对比验证，表明本文所建立的单乳液液滴碰撞的数值模型准确可靠，能够有效地模拟液滴在剪切流场中的碰撞动力学行为，为后续的研究提供了坚实的基础。3.1.7二维与三维模拟结果对比分别进行二维和三维的数值模拟，对比两者的结果，以分析维度对模拟结果的影响。在二维模拟中，假设液滴在x-y平面内运动，忽略z方向的变化；而在三维模拟中，考虑液滴在三维空间内的完整运动。在液滴的变形方面，二维模拟中液滴的变形主要集中在x-y平面内，变形形态相对较为简单；而三维模拟中，液滴在z方向也会发生一定的变形，变形形态更加复杂，呈现出三维的立体形状变化。例如，在二维模拟中，液滴在碰撞时可能只是在x-y平面内被压扁或拉长；而在三维模拟中，液滴除了在x-y平面变形外，还可能在z方向上出现一定的扭曲。在碰撞时间和碰撞力方面，二维模拟得到的碰撞时间和碰撞力与三维模拟结果也存在一定差异。由于三维模拟考虑了更多的自由度和相互作用，液滴之间的碰撞过程更加真实，碰撞时间可能会略有延长，碰撞力的分布也更加复杂。总体而言，三维模拟能够更全面地反映液滴在剪切流场中的碰撞动力学行为，但计算量相对较大；二维模拟虽然计算效率较高，但在描述液滴的复杂三维运动和相互作用时存在一定的局限性，在实际研究中需要根据具体问题和研究目的选择合适的模拟维度。3.2模拟结果与分析通过数值模拟，深入分析单乳液液滴在剪切流场中的运动行为。在剪切流场作用下，液滴的运动轨迹呈现出复杂的曲线形态，其运动方向不仅受到剪切流方向的影响，还与液滴之间的相互作用密切相关。初始时刻，液滴在剪切流的作用下开始发生变形，由球形逐渐变为椭圆形。随着时间的推移，液滴在流场拖曳力的作用下逐渐靠近，当液滴之间的距离足够小时，它们便会发生碰撞。在碰撞过程中，液滴受到强烈的相互挤推作用，产生剧烈的形变，在碰撞区形成明显的变形区域，液滴的形状发生显著改变，不再保持椭圆形，而是呈现出不规则的形状。碰撞后，液滴会相互绕过对方并分离，分离后的液滴形状逐渐恢复为椭圆状，但与碰撞前相比，其形状和运动轨迹仍存在一定的差异。以液滴在x方向与y方向的初始质心距分别为\Deltax_0/a=4，\Deltay_0/a=0.8为例，图1展示了液滴在剪切流场中经历碰撞过程的运动轨迹。从图中可以清晰地看到，随着剪切流场开始作用，液滴先变形为椭圆状（如t1时刻所示），随后在流场拖曳力的作用下，液滴相互靠近（t2-t3时刻），在t4时刻发生碰撞，此时液滴产生了剧烈形变，在碰撞作用处形成平滑界面。随后，两液滴以流场中心为对称轴旋转并相互绕过对方后分离（t5-t6时刻），当液滴完全分开后（t7时刻），液滴形状恢复为碰撞前的椭圆状，但液滴间的侧向质心间距\Deltay/a增大。[此处插入液滴在剪切流场中经历碰撞过程的运动轨迹图1][此处插入液滴在剪切流场中经历碰撞过程的运动轨迹图1]研究液滴间的相互作用发现，在碰撞过程中，液滴间的y向距离变化呈现出先增大后减小，随后稳定于一个值的规律，且最终的y向距离比初始的距离大。这是由于液滴在碰撞时，相互挤推作用使得液滴在y方向上产生了额外的速度分量，导致y向距离增大；而在碰撞后，随着液滴逐渐分离，这种额外的速度分量逐渐减小，y向距离又逐渐减小，最终稳定在一个新的值。为了更好地理解液滴碰撞过程中的流变行为，对碰撞过程中的压力分布进行分析。图2给出了与上述运动轨迹相对应时间点的碰撞过程压力分布时间序列图，图中的压力是以101325Pa为基准的相对压力值。从图中可以看出，随着施加外剪切流场，液滴开始变形，在表面张力的作用下，在变形液滴长轴方向的端部区域出现高压区（如t1时刻所示），从而在长轴端部的界面内外形成较高的压力梯度。在两液滴碰撞靠近过程中（t3-t4时刻），在碰撞区中心处产生一个高压区，这是由于液滴间的相互挤压，使得碰撞区域的流体受到压缩，压力升高。随着毛细数Ca的变化，界面挤压变形和液滴变形系数也会发生相应改变。当毛细数Ca从0.2增加到0.4时，界面挤压变形越明显，液滴变形系数D也从0.32增加到0.51，这表明毛细数越大，黏性力对液滴变形的影响越大，液滴越容易发生变形。[此处插入碰撞过程压力分布时间序列图图2][此处插入碰撞过程压力分布时间序列图图2]进一步探讨雷诺数、毛细数、液滴初始位置等因素对液滴碰撞动力学行为的影响。当雷诺数增大时，惯性力在液滴运动中所占的比重增加，液滴的运动速度加快，碰撞时的动能增大，导致液滴更容易发生破碎。在高雷诺数下，液滴在碰撞瞬间受到的惯性力可能超过表面张力和黏性力的合力，使得液滴无法保持完整而破碎成多个小液滴。毛细数对液滴变形和碰撞行为有着显著影响，随着毛细数的增大，黏性力逐渐超过表面张力，液滴的变形程度增大，碰撞时的接触面积和相互作用时间也会发生变化，从而影响液滴的碰撞结果。当毛细数较小时，表面张力主导液滴的形状，液滴在碰撞时变形较小；而当毛细数较大时，液滴在碰撞前就已经发生较大变形，碰撞时的相互作用更加复杂。液滴的初始位置对其碰撞动力学行为也至关重要。不同的初始质心距会导致液滴在剪切流场中的运动轨迹和碰撞时刻不同。当初始质心距较小时，液滴更容易在较短时间内发生碰撞，且碰撞时的相对速度和角度也会受到影响，进而影响碰撞后的运动状态。若两液滴初始质心距较小且处于同一流线附近，它们在剪切流作用下会很快相遇并发生碰撞，碰撞时的相对速度较小，可能导致液滴发生聚并；而当初始质心距较大时，液滴需要更长时间才能发生碰撞，碰撞时的相对速度可能较大，液滴更容易发生破碎或分离。为了更直观地展示各因素对液滴碰撞动力学行为的影响，绘制相图进行分析。在雷诺数-毛细数相图中（图3），可以清晰地看到不同区域对应着不同的液滴碰撞结果。在低雷诺数和低毛细数区域，液滴主要发生弹性碰撞，碰撞后液滴基本保持完整，变形较小；随着雷诺数和毛细数的增加，液滴逐渐进入变形和破碎区域，在该区域内，液滴碰撞时会发生明显的变形，甚至破碎成多个小液滴，且随着参数的进一步增大，液滴破碎的程度和概率也会增加。在液滴初始位置-碰撞结果相图中（图4），可以看出不同的初始位置会导致不同的碰撞概率和碰撞后液滴的运动状态。当初始质心距在一定范围内时，液滴碰撞概率较高，且碰撞后液滴的运动方向和速度会根据初始位置的不同而呈现出不同的分布规律。[此处插入雷诺数-毛细数相图图3][此处插入液滴初始位置-碰撞结果相图图4][此处插入雷诺数-毛细数相图图3][此处插入液滴初始位置-碰撞结果相图图4][此处插入液滴初始位置-碰撞结果相图图4]通过对相图的分析，可以总结出液滴碰撞动力学行为的特征。相图中的不同区域反映了各因素对液滴碰撞的综合影响，不同区域之间的边界可以看作是液滴碰撞行为发生转变的临界条件。在雷诺数-毛细数相图中，从弹性碰撞区域到变形和破碎区域的边界，对应着液滴从能够保持完整到开始发生破碎的临界状态，这个边界与雷诺数和毛细数的特定值相关，反映了惯性力、黏性力和表面张力之间的平衡关系发生改变时，液滴碰撞行为的变化。液滴初始位置-碰撞结果相图中，碰撞概率和碰撞后液滴运动状态的分布规律，为预测液滴在不同初始条件下的碰撞行为提供了直观的依据，有助于在实际应用中通过控制液滴的初始位置来调控乳液的性质和性能。四、剪切流场中单乳液液滴碰撞的可视化实验4.1实验系统与方法实验系统主要由剪切流场发生装置、乳液液滴注入系统、高速摄像及图像采集分析系统三大部分组成。剪切流场发生装置采用双板平行剪切结构，由两块尺寸为100mm\times50mm\times5mm的透明平行板组成，两板间距H可通过精密调节装置在1-5mm范围内精确调整。下板固定在实验平台上，上板由高精度电机驱动，可在水平方向以0.01-1m/s的速度匀速移动，从而产生稳定的剪切流场，剪切率G=\frac{U}{H}（U为上板移动速度）。为确保流场的稳定性和均匀性，在平行板的进出口处设置了整流装置，使流体能够平稳地进入和流出剪切区域。乳液液滴注入系统包括两个高精度注射泵（型号：YZB-100，流量范围：0.001-99.999ml/h）、两根内径为0.5mm的毛细管以及一个微流控芯片。两个注射泵分别用于输送连续相流体和分散相流体，通过精确控制注射泵的流量，能够实现对液滴生成频率和尺寸的精确控制。毛细管的一端与注射泵相连，另一端与微流控芯片的入口相接，微流控芯片采用聚二甲基硅氧烷（PDMS）材料制作，通过光刻和模塑工艺加工而成，其内部通道结构经过精心设计，能够在剪切流场中精确地生成单乳液液滴，并控制液滴的初始位置和运动方向。在微流控芯片的出口处，设置了一个液滴捕获装置，用于确保液滴能够准确地进入剪切流场的指定位置。高速摄像及图像采集分析系统配备一台高速摄像机（型号：Phantomv711，最高帧率：100000fps，分辨率：1280×800）以及一套图像采集分析软件（Image-ProPlus7.0）。高速摄像机安装在实验平台的正上方，通过调整镜头的焦距和角度，能够清晰地拍摄到剪切流场中液滴的碰撞过程。为了获得高质量的图像，在实验过程中采用了背光照明技术，通过在平行板下方放置一个高亮度的LED光源，使液滴在拍摄图像中呈现出清晰的轮廓。图像采集分析软件用于对拍摄到的图像进行处理和分析，能够自动识别液滴的位置、形状和尺寸，并跟踪液滴的运动轨迹。在图像分析过程中，采用了边缘检测算法和质心计算方法，以提高液滴参数测量的准确性。实验过程中，首先将连续相流体（如硅油，动力黏度\mu_c=0.1Pa\cdots，密度\rho_c=960kg/m^3）和分散相流体（如蓖麻油，动力黏度\mu_d=0.98Pa\cdots，密度\rho_d=961kg/m^3，界面张力系数\sigma=0.02N/m）分别装入对应的注射泵储液罐中。通过调节注射泵的流量，使连续相和分散相以一定的流速进入微流控芯片，在芯片内部的特定通道结构作用下，生成直径约为1mm的单乳液液滴。液滴生成后，在剪切流场发生装置产生的剪切流作用下，进入剪切区域。在实验操作中，通过调节上板的移动速度和两板间距，改变剪切率G，使其在10-100s^{-1}范围内变化；同时，通过调整注射泵的流量比，改变液滴的初始浓度和分布情况。在每个实验工况下，利用高速摄像机以5000fps的帧率对液滴在剪切流场中的碰撞过程进行连续拍摄，拍摄时间持续5s，以确保能够完整地记录液滴从进入剪切流场到发生碰撞及碰撞后一段时间内的运动和变形情况。数据采集主要通过图像采集分析软件完成，在拍摄结束后，将拍摄到的图像序列导入软件中，利用软件的自动识别和分析功能，提取每个液滴在不同时刻的位置坐标(x,y)、液滴的长轴长度L和短轴长度B。根据这些数据，计算液滴的变形系数D=\frac{L-B}{L+B}、液滴间的质心距\Deltar=\sqrt{(x_1-x_2)^2+(y_1-y_2)^2}以及液滴的运动速度v=\frac{\Deltar}{\Deltat}（\Deltat为相邻两帧图像的时间间隔）等关键参数。为了确保数据的准确性和可靠性，每个实验工况重复进行5次，对采集到的数据进行统计分析，取平均值作为最终结果，并计算数据的标准偏差，以评估数据的离散程度。4.2实验结果与讨论通过高速摄像记录的图像序列，能够清晰地观察到单乳液液滴在剪切流场中的碰撞过程。在剪切流场的作用下，液滴的运动轨迹呈现出复杂的曲线形状，这是由于液滴不仅受到剪切流的拖曳力，还受到周围流体的黏性力以及液滴间相互作用力的影响。初始时刻，液滴在剪切流的作用下开始逐渐变形，从球形变为椭圆形，这是因为剪切流在液滴表面产生了不均匀的应力分布，使得液滴沿流动方向被拉长。随着时间的推移，液滴在流场的作用下相互靠近，当液滴间的距离足够小时，便会发生碰撞。在碰撞瞬间，液滴会发生剧烈的形变，其形状变得极不规则。这是因为碰撞时液滴受到巨大的冲击力，导致液滴表面的应力急剧变化，使得液滴在碰撞点附近被压缩和扭曲。碰撞后，液滴可能会发生分离，也可能会发生聚并，这取决于液滴的表面性质、碰撞速度以及界面张力等多种因素。如果液滴表面的乳化剂能够有效地阻止液滴的聚并，且碰撞速度相对较低，液滴在碰撞后通常会相互分离；而当碰撞速度较大，或者乳化剂的作用较弱时，液滴可能会克服界面张力的阻碍，发生聚并形成更大的液滴。以实验中某一组工况为例，当剪切率G=50s^{-1}，液滴初始半径r=1mm时，通过对高速摄像图像的分析，得到了液滴在碰撞过程中的变形系数随时间的变化曲线，如图5所示。从图中可以看出，在碰撞前，液滴的变形系数随着时间逐渐增大，这是由于剪切流对液滴的持续作用导致液滴不断变形；在碰撞时刻，变形系数迅速增大，达到最大值，表明液滴在碰撞瞬间发生了剧烈的变形；碰撞后，变形系数逐渐减小，液滴的形状逐渐恢复，但最终仍会保持一定的变形，这是因为液滴在碰撞过程中消耗了部分能量，无法完全恢复到初始的球形状态。[此处插入液滴在碰撞过程中的变形系数随时间的变化曲线图图5][此处插入液滴在碰撞过程中的变形系数随时间的变化曲线图图5]进一步分析液滴的运动轨迹，通过图像分析软件对不同时刻液滴的质心位置进行追踪，得到了液滴在x-y平面内的运动轨迹图，如图6所示。从图中可以清晰地看到，液滴在剪切流场中沿着曲线运动，且两液滴的运动轨迹逐渐靠近，最终发生碰撞。在碰撞后，液滴的运动方向发生改变，这是由于碰撞过程中液滴间的相互作用力导致了动量的交换。同时，从运动轨迹中还可以观察到，液滴在y方向上的位移变化较为明显，这是因为剪切流在y方向上存在速度梯度，使得液滴在y方向上受到的作用力不均匀，从而产生了侧向位移。[此处插入液滴在x-y平面内的运动轨迹图图6][此处插入液滴在x-y平面内的运动轨迹图图6]为了研究不同因素对液滴碰撞行为的影响，对不同剪切率、液滴尺寸和流体物性参数等条件下的实验结果进行了对比分析。当剪切率增大时，液滴受到的剪切力增大，液滴的变形程度和碰撞速度都随之增加。在高剪切率下，液滴更容易发生破碎，这是因为较大的剪切力使得液滴表面的应力超过了液滴的承受能力，导致液滴破裂成多个小液滴。当剪切率从30s^{-1}增加到70s^{-1}时，液滴的破碎率从10\%增加到了30\%。液滴尺寸对碰撞行为也有显著影响，较大尺寸的液滴在碰撞时具有更大的惯性，碰撞后更难改变运动方向，且更容易发生聚并。这是因为大液滴的表面积与体积比较小，表面张力对其形状的维持作用相对较弱，在碰撞时更容易克服表面张力而发生聚并。实验结果表明，当液滴半径从0.8mm增大到1.2mm时，液滴的聚并率从20\%增加到了40\%。流体物性参数如黏度比和密度比也会影响液滴的碰撞行为。当黏度比增大时，分散相的黏度相对连续相增加，液滴的变形能力减弱，在碰撞过程中更倾向于保持自身形状，碰撞后的分离速度也会发生变化。如果连续相的黏度不变，分散相的黏度增大，液滴在碰撞时的变形程度会减小，碰撞后的分离速度可能会降低，这是因为较高的黏度使得液滴内部的流体流动受到限制，减少了液滴间的相互渗透和变形。密度比的变化会影响液滴在剪切流场中的沉降或上浮行为，进而影响液滴的碰撞概率和碰撞角度。当密度比偏离1较大时，液滴会在重力或浮力的作用下发生明显的垂直方向运动，改变了液滴在水平方向上的碰撞概率和碰撞角度，使得液滴的碰撞行为更加复杂。将实验结果与数值模拟结果进行对比，从液滴的运动轨迹来看，两者具有较好的一致性，数值模拟能够准确地预测液滴在剪切流场中的运动趋势和碰撞时刻。在液滴的变形系数方面，实验测量值与数值模拟结果在趋势上基本相同，但在具体数值上存在一定的偏差，这可能是由于实验过程中存在一些不可避免的误差，如液滴尺寸的测量误差、流体物性参数的不确定性以及实验装置的微小扰动等。数值模拟中采用的一些简化假设也可能导致与实际情况的差异。通过对比分析，进一步验证了数值模拟模型的可靠性，同时也为改进数值模拟方法提供了方向，例如在数值模拟中考虑更多的实际因素，提高模型的准确性，使其能够更精确地描述液滴在剪切流中的碰撞动力学行为。五、剪切流场中双乳液液滴碰撞的研究5.1双乳液液滴碰撞的数值模拟5.1.1数学模型建立物理模型描述：考虑在Couette流场中，两个相同的双乳液液滴发生碰撞的过程。双乳液液滴由外层连续相、中间层分散相和内层分散相组成，假设外层连续相和中间层分散相、内层分散相均为牛顿流体。内液滴半径为a_1，外液滴半径为a_2，初始时刻两双乳液液滴在x-y平面内关于中心点对称分布，流动方向为x方向，速度梯度方向为y方向，x方向与y方向的初始质心距分别为\Deltax_0和\Deltay_0。Couette流场由两个同心圆筒产生，内筒半径为R_1，以角速度\omega_1旋转，外筒半径为R_2，以角速度\omega_2旋转，两筒之间的流体形成剪切流场，作用于双乳液液滴，使其产生运动和变形。控制方程：在整个流场区域，质量守恒通过连续性方程来描述：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho为流体密度，\vec{v}为流体速度矢量，t为时间。对于不可压缩流体，\frac{\partial\rho}{\partialt}=0。动量方程考虑了流体的黏性、惯性以及界面张力的作用，表达式为：\rho(\frac{\partial\vec{v}}{\partialt}+\vec{v}\cdot\nabla\vec{v})=-\nablap+\nabla\cdot(2\mu\vec{D})+\vec{F}_{vol}其中，p为压力，\mu为动力黏度，\vec{D}为应变率张量，\vec{D}=\frac{1}{2}(\nabla\vec{v}+(\nabla\vec{v})^T)，\vec{F}_{vol}为界面张力源项。由于双乳液液滴存在两个界面，即外液滴与中间层的界面以及中间层与内液滴的界面，界面张力源项\vec{F}_{vol}采用CSF方法分别计算两个界面的贡献。对于外液滴与中间层的界面，\vec{F}_{vol1}=\sigma_1\kappa_1\nabla\alpha_1；对于中间层与内液滴的界面，\vec{F}_{vol2}=\sigma_2\kappa_2\nabla\alpha_2，其中\sigma_1和\sigma_2分别为两个界面的界面张力系数，\kappa_1和\kappa_2为相应界面的曲率，\alpha_1和\alpha_2为体积分数函数（对于外液滴与中间层界面，\alpha_1=1表示外液滴，\alpha_1=0表示中间层；对于中间层与内液滴界面，\alpha_2=1表示中间层，\alpha_2=0表示内液滴）。边界条件：在流场的进出口边界，采用周期性边界条件，即\vec{v}(x+L_x,y,z,t)=\vec{v}(x,y,z,t)，p(x+L_x,y,z,t)=p(x,y,z,t)，其中L_x为x方向的计算域长度，确保流场在进出口处的连续性，避免边界效应干扰液滴的碰撞过程。对于内外圆筒壁面，采用无滑移边界条件，即\vec{v}(r=R_1,\theta,z,t)=\vec{v}(r=R_2,\theta,z,t)=0，同时压力在壁面处满足法向应力连续条件，即p_{n}(r=R_1,\theta,z,t)=p_{n}(r=R_2,\theta,z,t)，保证流体在壁面处的速度和压力分布符合实际物理情况。数值求解方法：运用有限体积差分法对控制方程进行离散求解，通过将计算域划分为有限个控制体积，将连续的控制方程转化为离散的代数方程。压力场和速度场的耦合采用PISO（Pressure-ImplicitwithSplittingofOperators）算法，该算法相较于SIMPLE算法，在处理非稳态问题时具有更好的收敛性和计算效率，能够更准确地捕捉双乳液液滴在碰撞过程中的动态变化。动量方程采用二阶迎风差分格式离散，这种格式在处理对流项时具有更高的精度，能够更准确地描述流体的运动。基于PLIC方法的CompressiveInterfaceCapturingSchemeforArbitraryMeshes（CICSAM）格式对液/液界面进行重构，CICSAM格式在PLIC方法的基础上进行了改进，能够更好地处理复杂的界面形状和大变形问题，精确地捕捉双乳液液滴在碰撞过程中两个界面的演化。控制微分方程经离散成为线性代数方程组后，采用Bi-ConjugateGradientStabilized（BiCGSTAB）方法进行迭代求解。在迭代计算中，为了保证计算的收敛性和稳定性，采用了欠松弛因子，其中压力项欠松弛因子为0.3，密度项为0.6，体积力源项为0.6，动量源项为0.3。时间步长取0.05ms，当每个迭代时间步的流体速度残差小于10^{-4}时，认为该时间步迭代计算收敛，确保计算结果的准确性和可靠性。网格灵敏性分析：为确保数值解的网格独立性，采用不同的网格尺寸进行网格独立性检测。分别设置网格尺寸为\Deltax=\Deltay=\Deltaz=0.04a_2、\Deltax=\Deltay=\Deltaz=0.08a_2和\Deltax=\Deltay=\Deltaz=0.12a_2（a_2为外液滴半径）进行数值模拟，对比不同网格尺寸下双乳液液滴碰撞过程中的关键物理量，如液滴的变形系数、碰撞时间等。以液滴的变形系数D=(L-B)/(L+B)（L和B分别为发生变形后的椭圆形液滴的长轴与短轴）为例，当网格尺寸从\Deltax=\Deltay=\Deltaz=0.12a_2减小到\Deltax=\Deltay=\Deltaz=0.08a_2时，变形系数的计算结果变化较小，相对误差在4\%以内；当网格尺寸进一步减小到\Deltax=\Deltay=\Deltaz=0.04a_2时，变形系数的计算结果与\Deltax=\Deltay=\Deltaz=0.08a_2时相比，相对误差在1.5\%以内，表明此时网格尺寸对计算结果的影响已经很小，可以认为计算结果与网格无关。因此，在后续的数值模拟中，选择网格尺寸\Deltax=\Deltay=\Deltaz=0.08a_2，在保证计算精度的同时，提高计算效率，减少计算资源的消耗。模型验证：为验证所建立模型的准确性，对文献中相关的双乳液液滴实验进行数值模拟对比。在相同的实验条件下，包括相同的流场参数、液滴初始位置、流体物性参数等，将数值模拟得到的双乳液液滴碰撞演化过程与实验结果进行比较。从液滴的运动轨迹来看，数值模拟结果与实验结果基本一致，双乳液液滴在流场中的运动趋势、碰撞时刻以及碰撞后的分离方向等都能很好地吻合；从液滴的变形情况来看，数值模拟得到的双乳液液滴两个界面的变形形状和变形程度与实验图像中的液滴变形也具有较高的相似性。通过对比验证，表明本文所建立的双乳液液滴碰撞的数值模型准确可靠，能够有效地模拟双乳液液滴在剪切流场中的碰撞动力学行为，为后续的研究提供坚实的基础。5.1.2模拟结果与分析通过数值模拟，深入分析双乳液液滴在剪切流场中的运动行为。在剪切流场作用下，双乳液液滴的运动轨迹呈现出复杂的曲线形态，这是由于液滴不仅受到剪切流的拖曳力，还受到周围流体的黏性力以及液滴间相互作用力的影响。初始时刻，双乳液液滴在剪切流的作用下开始发生变形，外液滴和内液滴的形状都逐渐偏离球形。随着时间的推移，液滴在流场拖曳力的作用下逐渐靠近，当液滴之间的距离足够小时，它们便会发生碰撞。在碰撞过程中，双乳液液滴受到强烈的相互挤推作用，产生剧烈的形变。外液滴与中间层的界面以及中间层与内液滴的界面都会发生显著变形，在碰撞区形成复杂的变形区域。由于内液滴的存在，双乳液液滴的碰撞行为与单乳液液滴有所不同，内液滴会影响外液滴的变形和运动，使得碰撞过程更加复杂。碰撞后，双乳液液滴会相互绕过对方并分离，分离后的液滴形状逐渐恢复，但由于碰撞过程中的能量损失和界面变形，液滴的形状和运动轨迹仍会发生改变。以液滴在x方向与y方向的初始质心距分别为\Deltax_0/a_2=5，\Deltay_0/a_2=1为例，图7展示了双乳液液滴在剪切流场中经历碰撞过程的运动轨迹。从图中可以清晰地看到，随着剪切流场开始作用，双乳液液滴先变形，外液滴和内液滴的形状发生改变（如t1时刻所示），随后在流场拖曳力的作用下，液滴相互靠近（t2-t3时刻），在t4时刻发生碰撞，此时双乳液液滴产生了剧烈形变，两个界面都发生了明显的变形。随后，两液滴以流场中心为对称轴旋转并相互绕过对方后分离（t5-t6时刻），当液滴完全分开后（t7时刻），液滴形状逐渐恢复，但与碰撞前相比，液滴间的侧向质心间距\Deltay/a_2增大，且液滴的内部结构也发生了一定的变化。[此处插入双乳液液滴在剪切流场中经历碰撞过程的运动轨迹图7][此处插入双乳液液滴在剪切流场中经历碰撞过程的运动轨迹图7]研究双乳液液滴间的相互作用发现，在碰撞过程中，液滴间的y向距离变化呈现出先增大后减小，随后稳定于一个值的规律，且最终的y向距离比初始的距离大。这是因为液滴在碰撞时，相互挤推作用使得液滴在y方向上产生了额外的速度分量，导致y向距离增大；而在碰撞后，随着液滴逐渐分离，这种额外的速度分量逐渐减小，y向距离又逐渐减小，最终稳定在一个新的值。同时，由于内液滴的存在，液滴间的相互作用还会受到内液滴的位置和运动状态的影响。当内液滴位于外液滴的一侧时，会改变外液滴的受力分布，使得液滴间的相互作用更加复杂。为了更好地理解双乳液液滴碰撞过程中的流变行为，对碰撞过程中的压力分布进行分析。图8给出了与上述运动轨迹相对应时间点的碰撞过程压力分布时间序列图，图中的压力是以101325Pa为基准的相对压力值。从图中可以看出，随着施加外剪切流场，双乳液液滴开始变形，在表面张力的作用下，在外液滴长轴方向的端部区域和内液滴与中间层界面的端部区域出现高压区（如t1时刻所示），从而在这些区域的界面内外形成较高的压力梯度。在两液滴碰撞靠近过程中（t3-t4时刻），在碰撞区中心处产生一个高压区，这是由于液滴间的相互挤压，使得碰撞区域的流体受到压缩，压力升高。与单乳液液滴碰撞不同的是，双乳液液滴碰撞时，两个界面处的压力分布都发生了显著变化，内液滴与中间层界面处的压力变化也会影响外液滴的变形和运动。随着毛细数Ca的变化，界面挤压变形和液滴变形系数也会发生相应改变。当毛细数Ca从0.3增加到0.5时，界面挤压变形越明显，液滴变形系数D也从0.35增加到0.55，这表明毛细数越大，黏性力对双乳液液滴变形的影响越大，液滴越容易发生变形。[此处插入碰撞过程压力分布时间序列图图8][此处插入碰撞过程压力分布时间序列图图8]进一步探讨雷诺数、毛细数、液滴初始位置以及内液滴尺寸等因素对双乳液液滴碰撞动力学行为的影响。当雷诺数增大时，惯性力在液滴运动中所占的比重增加，液滴的运动速度加快，碰撞时的动能增大，导致液滴更容易发生破碎。在高雷诺数下，双乳液液滴受到的惯性力可能会使外液滴和内液滴的界面受到更大的拉伸和剪切作用，从而导致界面破裂，内液滴可能会从外液滴中释放出来，或者外液滴破碎成多个小液滴。毛细数对双乳液液滴变形和碰撞行为有着显著影响，随着毛细数的增大，黏性力逐渐超过表面张力，液滴的变形程度增大，碰撞时的接触面积和相互作用时间也会发生变化，从而影响液滴的碰撞结果。当毛细数较小时，表面张力主导液滴的形状，双乳液液滴在碰撞时变形较小；而当毛细数较大时，液滴在碰撞前就已经发生较大变形，碰撞时的相互作用更加复杂，内液滴与外液滴之间的相互作用也会更加明显，可能会导致内液滴在碰撞过程中发生位置偏移或变形。液滴的初始位置对其碰撞动力学行为也至关重要。不同的初始质心距会导致双乳液液滴在剪切流场中的运动轨迹和碰撞时刻不同。当初始质心距较小时，液滴更容易在较短时间内发生碰撞，且碰撞时的相对速度和角度也会受到影响，进而影响碰撞后的运动状态。若两双乳液液滴初始质心距较小且处于同一流线附近，它们在剪切流作用下会很快相遇并发生碰撞，碰撞时的相对速度较小，可能导致液滴发生聚并；而当初始质心距较大时，液滴需要更长时间才能发生碰撞，碰撞时的相对速度可能较大，液滴更容易发生破碎或分离。内液滴尺寸对双乳液液滴碰撞动力学行为也有重要影响。当内液滴尺寸增大时，内液滴的惯性增大，会对外液滴的变形和运动产生更大的阻碍作用。在碰撞过程中，较大尺寸的内液滴会使外液滴的变形更加困难，液滴间的相互作用也会发生改变。较大的内液滴可能会导致外液滴在碰撞时更容易发生破裂，因为内液滴的惯性会在碰撞瞬间对外液滴施加更大的作用力，使得外液滴的界面难以承受这种冲击力。内液滴尺寸的变化还会影响双乳液液滴的稳定性，当内液滴尺寸过大时，双乳液液滴可能会变得不稳定，更容易发生内液滴的泄漏或外液滴的破碎。为了更直观地展示各因素对双乳液液滴碰撞动力学行为的影响，绘制相图进行分析。在雷诺数-毛细数相图中（图9），可以清晰地看到不同区域对应着不同的双乳液液滴碰撞结果。在低雷诺数和低毛细数区域，双乳液液滴主要发生弹性碰撞，碰撞后液滴基本保持完整，内液滴和外液滴的界面变形较小；随着雷诺数和毛细数的增加，双乳液液滴逐渐进入变形和破碎区域，在该区域内，液滴碰撞时会发生明显的变形，外液滴和内液滴的界面都会发生显著变化，甚至可能导致内液滴的释放或外液滴的破碎，且随着参数的进一步增大，液滴破碎的程度和概率也会增加。在液滴初始位置-碰撞结果相图中（图10），可以看出不同的初始位置会导致不同的碰撞概率和碰撞后液滴的运动状态。当初始质心距在一定范围内时，双乳液液滴碰撞概率较高，且碰撞后液滴的运动方向和速度会根据初始位置的不同而呈现出不同的分布规律。同时，内液滴尺寸的变化也会在相图中体现出来，不同内液滴尺寸对应的碰撞结果区域会发生移动，反映出内液滴尺寸对双乳液液滴碰撞动力学行为的影响。[此处插入雷诺数-毛细数相图图9][此处插入液滴初始位置-碰撞结果相图图10][此处插入雷诺数-毛细数相图图9][此处插入液滴初始位置-碰撞结果相图图10][此处插入液滴初始位置-碰撞结果相图图10]通过对相图的分析，可以总结出双乳液液滴碰撞动力学行为的特征。相图中的不同区域反映了各因素对双乳液液滴碰撞的综合影响，不同区域之间的边界可以看作是液滴碰撞行为发生转变的临界条件。在雷诺数-毛细数相图中，从弹性碰撞区域到变形和破碎区域的边界，对应着双乳液液滴从能够保持完整到开始发生破碎5.2双乳液液滴碰撞的可视化实验双乳液液滴碰撞的可视化实验系统主要由微流控芯片、高速摄像系统、荧光显微镜以及数据采集与分析系统组成。微流控芯片采用聚二甲基硅氧烷（PDMS）材料制作，通过光刻和模塑工艺加工而成。芯片内部设计了特定的微通道结构，包括液滴生成通道和碰撞通道。液滴生成通道采用共流聚焦结构，通过控制连续相和分散相的流速，能够精确地生成大小均一、结构稳定的双乳液液滴。在共流聚焦结构中，连续相从两个对称的入口流入，将分散相夹在中间，在界面张力和剪切力的共同作用下，分散相被分割成微小的液滴，形成双乳液结构。碰撞通道则设计为T型结构，使生成的双乳液液滴能够在该通道中发生碰撞，通过调整T型通道的尺寸和液滴的注入速度，可以控制液滴的碰撞角度和速度。高速摄像系统选用高速摄像机（型号：Phantomv711，最高帧率：100000fps，分辨率：1280×800），其帧率可达10000fps，能够清晰地捕捉双乳液液滴碰撞的瞬间及后续的动态变化过程。为了获得高质量的图像，采用了高亮度的LED背光源，确保液滴在图像中具有清晰的轮廓。荧光显微镜则用于观察双乳液液滴的内部结构，通过对液滴内不同相添加荧光标记物，在荧光显微镜下可以清晰地分辨出内液滴、中间层和外液滴的形态和位置变化。数据采集与分析系统包括图像采集卡和专业的图像分析软件（Image-ProPlus7.0），能够实时采集高速摄像机拍摄的图像，并对图像进行处理和分析，提取液滴的尺寸、形状、运动轨迹等关键信息。实验过程中，连续相选用硅油，动力黏度\mu_c=0.1Pa\cdots，密度\rho_c=960kg/m^3；中间层分散相选用蓖麻油，动力黏度\mu_{d1}=0.98Pa\cdots，密度\rho_{d1}=961kg/m^3；内液滴选用水，动力黏度\mu_{d2}=0.001Pa\cdots，密度\rho_{d2}=1000kg/m^3。通过高精度注射泵（型号：YZB-100，流量范围：0.001-99.999ml/h）精确控制连续相、中间层分散相和内液滴的流速，从而实现对双乳液液滴尺寸和结构的精确调控。在液滴生成过程中，保持连续相流速不变，通过调整中间层分散相和内液滴的流速比，改变内液滴和中间层的相对尺寸，研究不同结构的双乳液液滴在碰撞过程中的行为差异。将生成的双乳液液滴引入碰撞通道，通过调整注射泵的流速，控制液滴的注入速度，从而改变液滴的碰撞速度。在碰撞通道中，液滴在剪切流场的作用下发生碰撞，利用高速摄像系统和荧光显微镜对碰撞过程进行实时观测和记录。在实验操作中，每个实验工况重复进行5次，以确保实验结果的可靠性和重复性。对采集到的图像数据进行处理和分析，利用图像分析软件自动识别液滴的轮廓和质心位置，计算液滴的变形系数、碰撞时间、碰撞角度等参数，并对这些参数进行统计分析，得到不同实验条件下双乳液液滴碰撞行为的规律。通过高速摄像记录的图像序列，清晰地观察到双乳液液滴在碰撞过程中的行为。在碰撞前，双乳液液滴在剪切流