[论文解读] Smooth particle methods without smoothing

Martin Campos Pinto|arXiv (Cornell University)|Dec 8, 2011

Fluid Dynamics Simulations and Interactions参考文献 19被引用 3

一句话总结

本文提出了一类新型粒子方法——线性变换粒子（LTP）与二次变换粒子（QTP）方法——通过后向特征流的多项式展开来变形粒子形状，实现无需重映射、扩展重叠或零矩条件的均匀收敛。关键结果表明，QTP粒子在极长模拟时间内无需重映射仍能保持高精度，且提出了一种动态重映射准则，可在最小化粒子间通信的同时减少数值耗散。

ABSTRACT

We present a new class of particle methods with deformable shapes that converge in the uniform norm without requiring remappings, extended overlapping or vanishing moments for the particles. The crux of the method is to use polynomial expansions of the backward characteristic flow to transport the numerical particles with improved accuracy: in the first order case the method consists of representing the transported density with linearly-transformed particles, the second order version computes quadratically-transformed particles, and so on. For programming purposes we provide explicit implementations of the resulting LTP and QTP schemes that only involve pointwise evaluations of the forward characteristic flow, and also come with local indicators for the accuracy of the corresponding transport scheme. Numerical tests using different transport problems demonstrate the accuracy of the proposed methods compared to standard particle schemes, and establish their robustness with respect to the remapping period. In particular, it is shown that QTP particles can be transported without remappings on very long periods of times, without hampering the accuracy of the numerical solutions. Finally, a dynamic criterion is proposed to automatically select the time steps where the particles should be remapped. The strategy does not require additional inter-particle communications, and it is validated by numerical experiments.

研究动机与目标

• 开发在L∞范数下实现均匀收敛的粒子方法，且无需重映射、扩展重叠或零矩条件。
• 通过后向特征流的多项式映射变形粒子形状，提升粒子格式在传输方程中的精度与鲁棒性。
• 提供与现有粒子代码兼容的实用实现，仅依赖于前向特征流的逐点求值。
• 设计一种动态重映射策略，无需额外粒子间通信即可选择最优重映射时机。
• 在多种传输问题中验证方法的精度与鲁棒性，尤其在无重映射的长时间模拟中。

提出的方法

• 该方法利用后向特征流的多项式展开定义粒子形状变换：LTP采用线性，QTP采用二次，一般r阶方法采用更高阶形式。
• 通过从后向流局部导数导出的系数，对粒子形状函数进行变形，从而在粒子传输中实现更高阶精度。
• LTP与QTP方法仅依赖于前向特征流的逐点求值实现，确保与标准粒子代码架构兼容。
• 基于局部误差指标（公式64与66）推导出动态重映射准则，可在无需全局通信的情况下估计传输误差。
• 重映射阈值由常数C_remapp控制，LTP取1，QTP取5，以在精度与计算成本之间取得平衡。
• 该方法避免平滑或卷积操作，从而保持粒子分辨率并最小化数值耗散。

实验结果

研究问题

• RQ1 粒子方法能否在不依赖重映射、扩展重叠或零矩条件的情况下，实现L∞范数下的均匀收敛？
• RQ2 通过后向流的多项式展开实现更高阶粒子形状变形，能否提升粒子传输中的精度与鲁棒性？
• RQ3 能否设计一种无需粒子间通信的动态重映射策略，以选择最优重映射时机？
• RQ4 LTP与QTP方法在精度与重映射频率方面，相较于标准重映射粒子方法表现如何？
• RQ5 QTP粒子能否在极长模拟时间内无需重映射保持高精度？

主要发现

• 在光滑性假设下，LTP与QTP方法在L∞范数下达到r阶收敛，且无需重映射、扩展重叠或零矩条件。
• QTP粒子可在无重映射条件下维持长达标准重映射粒子方法所需时间30倍的传输，同时保持高精度。
• 动态重映射策略在精度上可媲美最佳静态重映射策略，且重映射次数相近，同时无需额外粒子间通信。
• 在NLR测试案例中，动态策略表现略逊，原因在于一阶误差指标过于保守，提示采用二阶误差指标或可进一步提升性能。
• 该方法在多个二维传输问题中表现出鲁棒性，LTP的最优重映射周期比标准方法长10至50倍。
• 该实现仅需前向特征流的逐点求值，使得其与现有粒子代码集成简便，编程开销极小。