[論文レビュー] Voxelization based packing analysis for discrete element simulations of non-spherical particles

Venkata Rama Manoj Pola, Raghuram Karthik Desu+1 more|arXiv (Cornell University)|Oct 26, 2021

Textile materials and evaluations被引用数 5

ひとこと要約

本稿では、非球形粒子の局所的充填行動を分析するため、ボクセル化に基づく離散要素法（DEM）シミュレーション後の処理アルゴリズムを提案する。粒子集合体を3次元ボクセルに離散化することで、充填率の空間的変動を空間分解能をもって分析可能となり、壁の影響、積み重ね、空隙分布といった局所的効果が、スーパーキューブおよびマルチスフィア粒子において形状依存の充填不均一性について定量的な洞察を提供する。

ABSTRACT

A voxelization based post-processing algorithm is proposed to analyze the packing of non-spherical particle assemblies simulated using the Discrete Element Method. Voxelization of the particle data allows for isolating the geometric features of the granular assembly in various spatial sub-domains (2D surface or 3D region) and investigate the localized packing behaviour. Analyzing the local packing behaviour enables determining con-fined influences such as the wall-effect, stacking behaviour, local expansion/contraction and localized loading. The efficacy of the proposed technique to analyze practical granular assemblies is demonstrated through the packing analysis of different assemblies of superquadric cubes, superquadric ellipsoidals, and multi-spherical coffee beans.

研究の動機と目的

• DEM シミュレーションにおける非球形粒子の空間的分解能を有する充填解析の不足に対処すること。
• 粒子形状、壁に近い位置、幾何的制約による局所的充填変動の詳細な特徴付けを可能にすること。
• 複雑な顆粒集合体における充填構造を分析する汎用的で計算的に効率的な手法を開発すること。
• 粒子形状（例：立方体、楕円体、コーヒー豆型）が局所的充填密度および空隙分布に与える影響を定量化すること。
• ボクセルベースの空間分解を用いて、任意の空間的方向およびサブ領域における充填を分析するためのツールを提供すること。

提案手法

• ユーザーが指定した解像度の3次元グリッドにDEMでシミュレートされた粒子集合体をボクセル化し、離散的な空間セルにおける粒子幾何を表現する。
• 任意の法線ベクトルに沿ってボクセル化された粒子体積を2次元スライスに投影することで、任意の平面上の充填率を計算する。
• スーパーキューブおよびマルチスフィア粒子モデルを用いて、調整可能な幾何パラメータ（例：n₁, n₂, 球の重なり）を有する非球形形状を表現する。
• 空間フィルタリングおよび領域平均化を適用して、特定のサブドメイン（例：壁付近、コーナー、中心部）における充填行動を分離・分析する。
• ヒートマップおよび指定された方向に沿った1次元充填率プロファイルを生成し、空間的不均一性を可視化する。
• Pythonでアルゴリズムを実装し、再現性および再利用を目的として https://github.com/Manojpvr/VoxelDEM でオープンソース化した。

実験結果

リサーチクエスチョン

• RQ1 立方体、楕円体、コーヒー豆型などの粒子形状が、DEMでシミュレートされた集合体における局所的充填密度および空間的不均一性にどのように影響を与えるか？
• RQ2 壁に近い位置や幾何的閉塞が、非球形粒子系における局所的充填率にどの程度影響を及えるか？
• RQ3 ボクセル化に基づく後処理は、積み重ね、空隙形成、収縮／膨張といった局所的充填現象を効果的に捉え、可視化できるか？
• RQ4 任意の空間方向に沿って充填率はどのように変化し、異なる粒子形状ではどのようなパターンが現れるか？
• RQ5 粒子の曲率および不規則性は、秩序ある充填構造の促進または抑制に果たす役割は何か？

主な発見

• スーパーキューブ粒子の立方体形は、容器のコーナー付近で強い局所的密着を示し、ピーク充填率が約0.9に達する。
• 非球形粒子では壁付近の充填率が顕著に低く、コーヒー豆型および楕円体粒子では約0.1にまで低下する。これは平坦な壁面との接触性が悪いことが要因である。
• コーヒー豆型マルチスフィア粒子の平均充填率は0.6、楕円体スーパーキューブ粒子では0.62であり、全体的な密度は類似している。
• 楕円体およびコーヒー豆型粒子は中央部で平坦で均一な充填プロファイルを示すが、立方体粒子はより不均一な分布を示す。
• 楕円体およびコーヒー豆型粒子の曲率は、鋭い容器エッジとの幾何的不整合により、コーナー付近の充填が減少する原因となる。
• アルゴリズムは、おそらく自発的な整列により、楕円体粒子の周期的充填傾向を効果的に捉えている。一方、コーヒー豆型粒子は不規則性のため、より均一な空隙分布を示す。