[论文解读] Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part II. Numerical Results

Anthony J. C. Ladd|arXiv (Cornell University)|Jul 1, 1993

Lattice Boltzmann Simulation Studies被引用 9

一句话总结

本文提出了一种基于离散化玻尔兹曼方程方法的颗粒悬浮液数值模拟，该方法的计算成本与颗粒数量呈线性关系，从而实现了对胶体系统中流体动力学相互作用的高效建模。该方法准确模拟了无布朗运动和有布朗运动的蠕动流，以及有限雷诺数流，实现了最多1024个颗粒的稳定动力学，且传输系数和短时动力学均经过验证。

ABSTRACT

A new and very general technique for simulating solid-fluid suspensions has been described in a previous paper (Part I); the most important feature of the new method is that the computational cost scales with the number of particles. In this paper (Part II), extensive numerical tests of the method are described; for creeping flows, both with and without Brownian motion, and at finite Reynolds numbers. Hydrodynamic interactions, transport coefficients, and the short-time dynamics of random dispersions of up to 1024 colloidal particles have been simulated.

研究动机与目标

• 验证一种新型数值方法在颗粒数量上具有线性计算量的颗粒悬浮液模拟能力。
• 研究在无热涨落和有热涨落条件下，胶体悬浮液在蠕动流条件下的流体动力学相互作用。
• 将该方法扩展至有限雷诺数流，并评估其在捕捉复杂悬浮液动力学方面的准确性。
• 计算传输系数并分析最多1024个颗粒的无序分散体系中的短时动力学。
• 为软物质和流体动力学研究提供一种计算高效且通用的固体-流体悬浮液模拟框架。

提出的方法

• 采用格子玻尔兹曼方法，通过离散速度空间来模拟悬浮体系中的流体-颗粒相互作用。
• 实施一种边界条件方案，通过在规则格点上使用反弹规则，在颗粒表面强制实现无滑移条件。
• 通过在格子玻尔兹曼方程中引入随机力项来模拟布朗运动。
• 在D3Q19格点上求解离散玻尔兹曼方程，通过碰撞和流体步骤模拟流体流动和颗粒运动。
• 采用粒子解析方法，对单个颗粒进行追踪，并通过流场计算其流体动力学相互作用。
• 采用时间推进算法，以离散时间步长推进系统演化，确保质量与动量守恒。

实验结果

研究问题

• RQ1 所提出的方法能否准确模拟低雷诺数下胶体悬浮液中的流体动力学相互作用？
• RQ2 该方法在模拟悬浮液中颗粒动力学的布朗运动效应方面表现如何？
• RQ3 在有限雷诺数条件下，该方法在多大程度上保持稳定与准确？
• RQ4 模拟预测的传输系数（如自扩散系数、黏度）是多少，与理论预期相比如何？
• RQ5 计算成本如何随颗粒数量变化，该方法能否高效处理大规模系统（如1024个颗粒）？

主要发现

• 该方法实现了计算成本与颗粒数量的线性增长，支持最多1024个胶体颗粒的模拟。
• 在蠕动流中，流体动力学相互作用被准确捕捉，预测的自扩散系数与理论预测一致。
• 引入随机力项成功再现了布朗运动效应，验证了热涨落模型的正确性。
• 有限雷诺数条件下的模拟显示出稳定且物理解释一致的流场和颗粒轨迹。
• 无序悬浮液的短时动力学，包括颗粒扩散和相关函数，与理论模型良好吻合。
• 该方法在从低到中等雷诺数的各类流态下均表现出强健性与准确性。