[논문 리뷰] Voxelization based packing analysis for discrete element simulations of non-spherical particles

Venkata Rama Manoj Pola, Raghuram Karthik Desu외 1명|arXiv (Cornell University)|2021. 10. 26.

Textile materials and evaluations인용 수 5

한 줄 요약

이 논문은 비구형 입자를 위한 이산 요소법(DDEM) 시뮬레이션의 후처리 알고리즘으로, 볼륨 요소(voxelization) 기반의 방법을 제안한다. 입자 집합을 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 시뮬레이션에서 비구형 입자에 대한 공간 해상도가 확보된 밀도 분석의 부족을 보완하기 위해.
• 입자 형태, 벽과의 근접도, 기하학적 제약 등에 의해 유도되는 국소적 밀도 변화의 세부 특성 분석을 가능하게 하기 위해.
• 복잡한 입자 집합에서의 밀도 구조 분석을 위한 일반화 가능하고 계산 효율성이 높은 방법 개발을 위해.
• 입자 형태(예: 정육면체, 타원체, 커피 콩)가 국소적 밀도 및 공극 분포에 미치는 영향을 정량화하기 위해.
• 볼륨 기반 공간 분할을 이용해 임의의 공간 방향 및 부분 영역에서의 밀도 분석을 위한 도구 제공을 위해.

제안 방법

• 사용자 정의 해상도로 DEM 시뮬레이션 결과 입자 집합을 3차원 격자로 볼륨화하여, 이산 공간 셀 내에서 입자 기하구조를 표현하기 위해.
• 임의의 법선 벡터를 따라 3D 볼륨을 2D 슬라이스로 투영함으로써, 평면 상의 밀도 분율을 계산하기 위해.
• 초구형 및 다구형 입자 모델을 사용하여 기하학적 매개변수(예: n₁, n₂, 구의 겹침)를 조절할 수 있도록 비구형 형태를 표현하기 위해.
• 공간 필터링 및 지역 평균화를 적용하여 특정 하위 영역(예: 벽 근처, 모서리, 중심부)의 밀도 행동을 분리 및 분석하기 위해.
• 공간 비균형성을 시각화하기 위해 히트맵 및 특정 방향에 따른 1D 밀도 분율 프로파일 생성하기 위해.
• 알고리즘을 파이썬으로 구현하였으며, 재현성 및 재사용을 위해 https://github.com/Manojpvr/VoxelDEM 에 오픈소스로 배포하였다.

실험 결과

연구 질문

• RQ1 입자 형태(정육면체, 타원체, 커피 콩 등)가 DEM 시뮬레이션 결과 입자 집합의 국소적 밀도 및 공간 비균형성에 어떤 영향을 미치는가?
• RQ2 벽과의 근접도 및 기하학적 제약이 비구형 입자 시스템의 국소적 밀도 분율에 어느 정도의 영향을 미치는가?
• RQ3 볼륨화 기반 후처리가 적층, 공극 형성, 수축/팽창과 같은 국소적 밀도 현상을 효과적으로 포착하고 시각화할 수 있는가?
• RQ4 임의의 공간 방향을 따라 밀도 분율은 어떻게 변화하며, 다양한 입자 형태에서 어떤 패턴이 나타나는가?
• RQ5 입자 곡률과 비정형성은 정렬된 밀도 구조의 형성 여부에 어떤 역할을 하는가?

주요 결과

• 초구형 입자 중 정육면체 형태의 입자는 용기 모서리 부근에서 강한 국소적 밀도 증가를 보이며, 최대 밀도 분율이 약 0.9에 도달한다. 이는 적층 행동에 기인한다.
• 비구형 입자에서는 벽 부근의 평균 밀도 분율이 현저히 낮으며, 커피 콩 및 타원체 입자에서 약 0.1 수준이다. 이는 평탄한 벽면과의 접촉이 열악하기 때문이다.
• 커피 콩 형태의 다구형 입자 평균 밀도 분율은 0.6이며, 타원체 형태의 초구형 입자에서는 0.62로 나타나 전체 밀도는 유사하다.
• 타원체 및 커피 콩 입자는 중심 영역에서 더 평탄하고 균일한 밀도 프로파일을 보이며, 정육면체 입자에 비해 더 높은 비균형성 분포를 보인다.
• 타원체 및 커피 콩 입자의 곡률은 날카로운 용기 모서리와의 기하학적 불일치로 인해 모서리 부근의 밀도가 감소한다.
• 알고리즘은 타원체 입자에서 주기적인 밀도 경향을 성공적으로 포착하였으며, 이는 자발적 정렬에 기인할 가능성이 높다. 반면 커피 콩 입자는 비정형성으로 인해 더 균일한 공극 분포를 보인다.