Efficient generation of large-scale non-equilibrium distributions of particles

技術的概要：大規模非平衡粒子分布の効率的生成

問題提起
周期的微構造を有する代表体積要素（RVE）は、複合材料において微視的不均質性と巨視的有効特性を結びつけるために不可欠である。従来のアルゴリズム（モンテカルロ法、ルバチェフスキー・スティリンガー事象駆動分子動力学など）は、平衡状態または準平衡状態の粒子配置の生成においては優れているが、糸巻き成形や射出成形などの現実の製造プロセスに特徴的な、不均一で異方的かつ履歴依存性の分布を捉えることにはしばしば失敗する。既存の手法は、実際の複合材料に見られる不規則性、欠陥、制御された無秩序性を抑制する傾向がある。さらに、高粒子数を伴う大規模微構造の生成は、衝突検出における二次スケーリングに起因する過大な計算コストに悩まされることが多い。ランダム性の調整可能な程度、制御された空間相関、およびクラスター化や凝集を含む現実的な非平衡特徴を備えた周期的 RVE を効率的に生成できるアルゴリズムの必要性が明確に存在する。

手法：膨張とランダム移動（SRM）アルゴリズム
本論文では、修正されたミラー・トルクエータ「膨張」技術と修正されたランダム逐次追加（RSA）方式を組み合わせたハイブリッド手法である、膨張とランダム移動（SRM）アルゴリズムを提案する。この手法は、単分散のインクルージョンを含む単位ボックス上で動作し、以下の反復ステップを経て進行する。

1. 初期化： 古典的な RSA を用いて、低体積分率において重なりのない初期配置を生成する。
2. 膨張： すべての粒子の半径を小さな増分（$c_w$）だけ同時に増加させる。
3. ランダム移動： 各粒子を小さなランダムベクトル（$c_m$）だけ変位させる。
4. 局所許容性チェック： 適応型セルベースの空間分割方式を用いて、隣接粒子までの距離を評価する。粒子が重なり合う場合、その粒子は以前の位置に戻され、新しいランダム変位が試行される。このループは、有効な位置が見つかるか、最大試行回数（$N_l$）に達するまで継続する。
5. シェイクステップ： 粒子が数回試行しても許容される位置を見つけられない場合、全球的な緩和（「シェイク」）が実行される。これは、すべての粒子をランダムに変位させ、特にクラスター化された配置における局所的なジャミングを克服するために重なりをテストするものである。
6. 終了： 目標体積分率が達成されるまで、このサイクルを繰り返す。

計算効率
手法の重要な構成要素は、セルベースの隣接探索の実装である。RVE を粒子直径よりわずかに大きいセルに分割することで、アルゴリズムは現在のセルおよびその直近の隣接セル（2 次元では $3 \times 3$ ブロック）内の粒子にのみ重なりチェックを制限する。これにより、衝突検出の計算複雑度が粒子数に対して二次から準線形（$O(N)$）に低減される。キャッシュミスを防ぐための粒子アクセスパターンの追加最適化により、生成時間は約 $t \propto N^{1.15}$ としてスケーリングし、2 次元および 3 次元の両方で最大 $10^7$ 個の粒子を用いたシミュレーションを可能にする。

主要な結果

• スケーラビリティ： 最適化された SRM アルゴリズムは、最大 $10^7$ 個の粒子を有する配置の生成に成功した。低から中程度の体積分率において、この手法は高い効率を維持する。
• 非平衡制御： 膨張率（$c_w$）と移動率（$c_m$）を調整することで、アルゴリズムは平衡に近い分布から強くクラスター化された配列まで、多様な微構造のスペクトルを生成できる。
  • 分子動力学（MD）との比較： 自然に平衡へと進化させる事象駆動分子動力学（MD）とは異なり、SRM は局所体積分率（LVF）分布が著しく広範で、より強いクラスター化のシグネチャを持つ配置を生成する。MD によって生成された「クラスター化」状態は、SRM によって生成されたものよりも平衡に近い。
  • 体積分率の限界： 2 次元では、この手法は古典的なランダム密充填（RCP）限界を超え（ディスクの場合、約 87% を達成）、3 次元では実用上の限界は球の場合、約 60% である。これを超えるとランダム移動は非効率的となり、事象駆動 MD が好まれる。
• 非球形インクルージョンへの拡張： このアルゴリズムは、薄い円形プレートレット（球状ディスク）に適用された。パラメータを調整することで、著者らは以下のものを生成した。
  • 「カードハウス（HoC）」ネットワーク： 高い膨張率と低い移動率によって形成される、高度に相互接続された構造。
  • 準ネマティックドメイン： 高い移動率で HoC 構造を緩和することによって形成される、準安定で部分的に秩序化した領域。
• 構造 - 特性への影響： プレートレット微構造のパーコレーション分析により、HoC 配置の臨界パーコレーション距離は、緩和された準ネマティック構造よりも約 1 桁低いことが明らかになった。これは、微視的配列が輸送特性に著しく影響することを示している。

意義と主張
本論文は、SRM アルゴリズムを複合材料の計算モデリングのための強力かつ柔軟なツールとして位置づけている。その主な意義は、実用的な製造において観察される統計的変動と非平衡特徴を反映する、現実的で多様かつ計算的にアクセス可能な粒子配置を生成する能力にある。

著者らは、SRM が純粋に幾何学的な手法であることを強調している。粒子の再配置は物理的な力や粒子間ポテンシャルによって支配されない。したがって、SRM は均質化研究、不確実性定量化、仮想試験のための統計的に代表的な初期微構造を生成するには理想的であるが、粒子系の真の物理的進化を再現したり、相転移を研究したりするものではない。この手法の価値は、制御された方法で微視的無秩序を変化させ、製造に起因する異方性を模倣することにより、不均質媒体におけるより正確な構造 - 特性関係分析を可能にする点にある。