[论文解读] Particle-based simulations of steady-state mass transport at high Péclet numbers

Thomas Müller, Paolo Arosio+6 more|arXiv (Cornell University)|Oct 17, 2015

Groundwater flow and contamination studies被引用 3

一句话总结

本文提出一种基于粒子的模拟方法，用于在高佩克莱特数条件下模拟稳态质量传输，通过根据通量分布而非随机浓度初始化粒子。通过将这些通量加权粒子沿蒙特卡罗轨迹传播，并检测其在预设区域的通过情况，该方法能准确预测下游浓度分布，与使用荧光胶体的实验微流控数据高度一致，并通过25 nm的粒子半径验证了该方法的有效性。

ABSTRACT

Conventional approaches for simulating steady-state distributions of particles under diffusive and advective transport at high Péclet numbers involve solving the diffusion and advection equations in at least two dimensions. Here, we present an alternative computational strategy by combining a particle-based rather than a field-based approach with the initialisation of particles in proportion to their flux. This method allows accurate prediction of the steady state and is applicable even at high Péclet numbers where traditional particle-based Monte-Carlo methods starting from randomly initialised particle distributions fail. We demonstrate that generating a flux of particles according to a predetermined density and velocity distribution at a single fixed time and initial location allows for accurate simulation of mass transport under flow. Specifically, upon initialisation in proportion to their flux, these particles are propagated individually and detected by summing up their Monte-Carlo trajectories in predefined detection regions. We demonstrate quantitative agreement of the predicted concentration profiles with the results of experiments performed with fluorescent particles in microfluidic channels under continuous flow. This approach is computationally advantageous and readily allows non-trivial initial distributions to be considered. In particular, this method is highly suitable for simulating advective and diffusive transport in microfluidic devices.

研究动机与目标

• 解决传统基于粒子的蒙特卡罗方法在高佩克莱特数下因初始化不佳而导致的计算效率低下和失效问题。
• 开发一种无需求解福克-普朗克方程的模拟策略，而是直接从基于通量的初始条件出发模拟粒子轨迹。
• 实现在复杂几何结构中，特别是具有泊肃叶流的微流体器件中的稳态浓度分布的准确预测。
• 通过微通道中荧光胶体的实验验证该方法，并将模拟结果与实测数据进行对比。
• 通过使用实测的初始概率分布作为输入，最大限度减少拟合参数，证明该方法的鲁棒性。

提出的方法

• 粒子在 t=0 且 x=0 处根据通量加权分布 Φ(x=0,y,z) = C(x=0,y,z)·v(y,z)/v̄ 初始化，其中 C 为初始浓度，v 为局部流体速度，v̄ 为平均速度。
• 每个粒子的轨迹通过朗之万方程模拟，结合了热运动的随机性与基于流场的确定性输运。
• 通过统计粒子在下游预设检测区域中停留的时间步数，重建稳态浓度分布。
• 该方法通过仅从单个初始时间步模拟所有粒子，避免了持续的粒子注入，从而降低了计算成本。
• 该方法考虑了微通道中的抛物线速度分布，并消除了靠近通道壁面粒子运动缓慢所导致的检测偏差。
• 模拟基于第一性原理物理，仅将粒子半径作为自由参数，且已与实验测得的初始分布进行验证。

实验结果

研究问题

• RQ1 基于粒子的模拟方法是否能准确预测高佩克莱特数下传统蒙特卡罗方法失效的稳态质量传输？
• RQ2 与随机浓度初始化相比，根据粒子通量分布初始化是否能提高模拟对输运-扩散过程的准确性与效率？
• RQ3 在不进行拟合的情况下，实验测得的初始概率分布能在多大程度上用作输入以预测下游浓度分布？
• RQ4 模拟的粒子分布与微流体通道中实验的荧光显微镜数据匹配程度如何？
• RQ5 模拟结果与实测的荧光胶体束扩散展宽之间在定量上的一致性如何？

主要发现

• 该方法在微流体通道中实现了模拟结果与实验测得浓度分布之间的出色定量一致性，当使用已知的25 nm粒子半径时匹配最佳。
• 在下游20 mm处，25 nm粒子模拟的归一化平方误差为380，表明其与实验噪声水平的差异约为20倍以内。
• 使用10 nm或40 nm粒子半径的模拟得到的归一化误差达到数千至数万，证实了当使用正确半径时，该方法具有高度敏感性与准确性。
• 该方法成功捕捉了荧光胶体束在80 mm传播距离内的扩散展宽，与实验图像在多个下游位置高度一致。
• 该方法计算高效且可扩展，因为计算粒子通量的成本与系统几何形状和边界条件的依赖性极弱。
• 验证结果表明，该模拟并非对数据的拟合，而是仅基于实测初始分布和粒子尺寸的第一性原理预测。