论文部分内容阅读
在微流体动力学领域中,介观尺度下包含了一系列流体分子的微泡,在诸多领域,如药物的制备与提纯,香水的制造,化工制剂的分离过滤,药物输运,声孔效应,超声空化清洁,空化初生,梢涡空化等,有着极大的研究和应用价值,其中很多的动力学问题还未解决并且亟待解决。 耗散粒子动力学(DPD: Dissipative Particle Dynamics)作为一种自下而上基于粒子法的粗粒化介观数值方法,已经成功应用于各种介观系统的研究。通过使用单个粒子来表示多个分子或分子团的粗粒化过程,从而抛弃了与我们所关心的性质无关的更小空间和时间尺度的细节,消除了它们所产生的不必要的计算资源消耗。DPD粒子之间的相互作用力以成对的形式出现,由保守力、耗散力和随机力三种组成。耗散力和随机力通过刻意地选择参数和权重函数来满足涨落耗散定理,构成了DPD方法中内置的恒温器。多体耗散粒子动力学(MDPD: Many-body Dissipative Particle Dynamics)方法作为标准DPD方法扩展版本,通过修正保守力的形式以添加了与密度相关的成分,使得DPD流体能够模拟出气-液两相共存状态。 本文通过引入分子动力学(MD: Moleculars Dynamics)方法中的Berendsen恒压器,发现原始的Berendsen恒压器在DPD方法和MDPD方法中能够实现快速的压力变化控制,对非平衡动力学过程的研究有着独具一格的优势。对于单相系统能够成为很好的压力控制器,但是作用于多相系统则会对两相界面造成“非物理”的干扰甚至破坏。针对该问题作者提出局部的Berendsen恒压器,很好地弥补了原始Berendsen恒压器的该弱点,使Berendsen恒压器适用于多相系统,可以实现准确的两相界面的模拟。两相界面上不会因恒压器进行压力控制而形成“非物理”过渡,并且局部的Berendsen恒压器仍然保持原始Berendsen恒压器的非平衡动力学过程的特性优势。 为了实现微气泡的(M)DPD方法模拟,通过分析不包含任何粒子的真空泡模型在模拟微气泡动力学问题存在的不足,无法捕捉较为复杂的动力学行为,如微气泡的振荡和溃灭等,论证了包含粒子的气体泡模型对于复杂微气泡动力学行为模拟的必要性。通过对比在包含粒子的(M)DPD气体泡模型保守力中排斥力的影响参数对液气密度比的影响,发现采用DPD气体泡模型通过增加保守力斥力截断半径更容易达到微气泡模拟的高液气密度比的要求,本文中得到的最大液气密度比约为11。 通过对于微泡在平板泊肃叶流中迁移运动的研究,模拟的结果表明较大的粘度比导致近壁的平衡位置,同时较大的毛细数(Ca)使微泡迁移到通道中心线。对相同Ca数的微泡迁移,通道有一个临界宽度,如果通道宽度小于该临界值,则微泡的平衡位置只能定位在通道的中心线处。