The power particle-in-cell method
The power particle-in-cell method
复制标题
DOI：
10.1145/3528223.3530066
复制
发表时间：
2022-07
期刊：
影响因子：
--
通讯作者：
Ziyin Qu;Minchen Li;F. D. Goes;Chenfanfu Jiang
中科院分区：
文献类型：
--
作者：
Ziyin Qu;Minchen Li;F. D. Goes;Chenfanfu Jiang
本文介绍了一种新的加权方案的粒子网格传输，产生混合拉格朗日/欧拉流体模拟均匀的粒子分布和精确的体积控制。在其核心，我们的方法通过计算体积约束的密度核来重新制定Power Particles的构建[de Goes et al. 2015]。我们采用这些优化的内核作为广义插值材料点方法（GIMP）中的粒子域，以便将Power Particles纳入Particle-In-Cell框架，因此命名为Power Particle-In-Cell方法。我们将体积约束密度核的构造作为正则化的最优运输问题，并描述了一种基于局部高斯卷积的迭代求解器，与[de Goes et al. 2015]相比，该迭代求解器可显著提高性能。我们还提出了新的扩展处理自由表面和固体的障碍，绕过细胞裁剪和鬼粒子的需要。我们证明了我们的传输权重的优点，通过改进混合方案的流体模拟，如流体隐式粒子（FLIP）方法和仿射粒子在细胞（APIC）方法与体积保持和鲁棒性变化的粒子每细胞的比例，同时保持低数值耗散，保存线性和角动量，并避免粒子重新播种或后处理松弛。
This paper introduces a new weighting scheme for particle-grid transfers that generates hybrid Lagrangian/Eulerian fluid simulations with uniform particle distributions and precise volume control. At its core, our approach reformulates the construction of Power Particles [de Goes et al. 2015] by computing volume-constrained density kernels. We employ these optimized kernels as particle domains within the Generalized Interpolation Material Point method (GIMP) in order to incorporate Power Particles into the Particle-In-Cell framework, hence the name the Power Particle-In-Cell method. We address the construction of volume-constrained density kernels as a regularized optimal transportation problem and describe an iterative solver based on localized Gaussian convolutions that leads to a significant performance speedup compared to [de Goes et al. 2015]. We also present novel extensions for handling free surfaces and solid obstacles that bypass the need for cell clipping and ghost particles. We demonstrate the advantages of our transfer weights by improving hybrid schemes for fluid simulation such as the Fluid Implicit Particle (FLIP) method and the Affine Particle-In-Cell (APIC) method with volume preservation and robustness to varying particle-per-cell ratio, while retaining low numerical dissipation, conserving linear and angular momenta, and avoiding particle reseeding or post-process relaxations.