[論文レビュー] Smooth particle methods without smoothing

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

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

ひとこと要約

本稿では、後退的特徴線流れの多項式展開を用いて粒子形状を変形させる、線形および2次変換粒子（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は2次、一般にはr次スキーム。
• 粒子形状関数は、後退的流れの局所的微分から得られる係数を用いて変形され、粒子輸送における高次精度を実現する。
• LTPおよびQTPスキームは、前方特徴線流れの点関数評価のみに依存して実装されており、標準的な粒子コードアーキテクチャと互換性を持つ。
• 動的リマップ基準は、局所誤差指標（式64および66）に基づいて導出され、グローバル通信を要せず輸送誤差を推定する。
• リマップしきい値は定数C_remappによって制御され、LTPでは1、QTPでは5として設定されており、精度と計算コストのバランスを取る。
• 平滑化や畳み込みを避けることで、粒子分解能を保持し、数値的拡散を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 リマップ、拡張重なり、モーメントの消滅を要件としない粒子法が、L∞ノルムにおいて一様収束を達成できるか？
• RQ2 後退的流れの多項式展開による高次粒子形状変形が、粒子輸送における精度と頑健性を向上させられるか？
• RQ3 粒子間通信を伴わず、最適なリマップタイミングを選択できる動的リマップ戦略を設計できるか？
• RQ4 LTPおよびQTPスキームの性能は、標準的なリマップ付き粒子法と比較して、精度とリマップ頻度の面でどのように異なるか？
• RQ5 QTP粒子は、リマップなしで非常に長いシミュレーション時間にわたり高い精度を維持できるか？

主な発見

• LTPおよびQTPスキームは、滑らかさの仮定の下で、リマップ、拡張重なり、モーメントの消滅を要件とせず、r次のL∞収束を達成する。
• QTP粒子は、標準的なリマップ付き粒子法が要求する時間の約30倍まで、リマップなしで輸送可能であり、高い精度を維持する。
• 動的リマップ戦略は、最良の静的リマップ戦略と同等の精度を達成し、ほぼ同じ回数のリマップで実現可能であり、追加の粒子間通信を伴わない。
• NLRテストケースでは、誤差指標が過度に楽観的である一次誤差指標のため、動的戦略の性能がわずかに劣る。これは、二次誤差指標の導入で性能向上が期待できる。
• 本手法は複数の2次元輸送問題において頑健であり、LTPでは標準手法と比較して最適なリマップ周期が10〜50倍長くなる。
• 実装は前方特徴線流れの点関数評価のみを必要とし、既存の粒子コードへの統合は、最小限のプログラミングオーバヘッドで容易に実現可能。