[논문 리뷰] MC-DEM: a novel simulation scheme for modeling dense granular media

Nicolas Brodu, J.A. Dijksman외 1명|arXiv (Cornell University)|2014. 10. 23.

Granular flow and fluidized bedsEngineering참고 문헌 16인용 수 19

한 줄 요약

이 논문은 다중 접촉 간의 상호작용을 선형 탄성 기반 변형을 통해 명시적으로 모델링하는 새로운 이산요소법 MC-DEM를 소개한다. 이는 밀도가 높은 입자계에서 힘 예측을 향상시키기 위한 것이다. 이 방법은 3차원 실험 데이터의 힘 분포 및 스케일링 법칙과 뛰어난 일치를 보이며, 유한요소법(FEM)을 사용하지 않고도 전통적인 DEM보다 뛰어난 성능을 발휘한다.

ABSTRACT

This article presents a new force model for performing quantitative simulations of dense granular materials. Interactions between multiple contacts (MC) on the same grain are explicitly taken into account. Our readily applicable method retains all the advantages of discrete element method (DEM) simulations and does not require the use of costly finite element methods. The new model closely reproduces our recent experimental measurements, including contact force distributions in full 3D, at all compression levels up to the experimental maximum limit of 13\%. Comparisons with traditional non-deformable spheres approach are provided, as well as with alternative models for interactions between multiple contacts. The success of our model compared to these alternatives demonstrates that interactions between multiple contacts on each grain must be included for dense granular packings.

연구 동기 및 목표

• 밀도가 높은 입자계에서 전통적인 DEM 시뮬레이션과 실험 데이터 간의 낮은 정량적 일치도를 해결하기 위해.
• 표준 DEM에서 간과되는 한 입자 위의 다중 접촉 간 상호 영향을 모델링하기 위해.
• DEM의 단순성은 유지하면서도 밀도가 높고, 거의 정지 상태에 가까운 시스템에 대해 정확도를 향상시킬 수 있는 계산 효율성이 높은 방법을 개발하기 위해.
• 접촉력 및 거시적 반응에 대한 전체 3차원 실험 측정치와 모델을 검증하기 위해.
• 입자 체계의 변형과 접촉 간 힘 결합이 밀도가 높은 입자 매체의 정확한 시뮬레이션에 필수적임을 입증하기 위해.

제안 방법

• 구형 입자에 작용하는 점력에 대한 해석적 선형 탄성 해를 사용하여 국소 입자 변형을 계산한다.
• 다른 접촉 k에서 발생한 힘이 동일한 입자 위의 접촉 c에서 유도하는 법선 변위 δk→c 를 보정 항으로 도입한다.
• 효과적 오버랩은 δc + Σk δk→c 로 갱신되며, 이는 후속 허츠 힘 법칙 F ∝ (δc + Σk δk→c)^3/2 에 사용된다.
• 교착 상태에 영향을 주는 경계 조건 효과를 보완하기 위해 경험적 비율 계수 γ 를 사용하여 실험 데이터에 대한 적합도를 향상시킨다.
• 모든 변형은 장거리 상호작용이나 유한요소법(FEM) 계산 없이 입자 수준에서 계산된다.
• 형상-변형 관계가 알려진 다른 입자계 시스템, 예를 들어 에멀전 또는 거품 등에도 일반화 가능하다.

실험 결과

연구 질문

• RQ1 한 입자 위의 다중 접촉 상호작용을 고려한 수정된 DEM 모델이 밀도가 높은 3차원 입자계에서 실험적 힘 분포를 더 잘 재현할 수 있는가?
• RQ2 탄성 변형을 통한 접촉 간 힘 결합을 포함할 경우, 밀도가 높은 입자 매체의 거시적·미시적 반응에 어떤 영향을 미치는가?
• RQ3 표준 DEM에서의 독립적 접촉 가정이 임계 정지 상태 근처의 밀도가 높은 입자계의 진정된 기계적 거동을 얼마나 잘못 반영하는가?
• RQ4 비용이 많이 드는 FEM 기반 시뮬레이션을 대체할 수 있는 단순하고 해석적인 다중 접촉 효과 보정이 고정밀도를 유지할 수 있는가?
• RQ5 입자 변형과 힘 전파가 밀도 및 변형률에 따라 접촉력의 스케일링에 어떤 영향을 미치는가?

주요 결과

• MC-DEM는 3DXP 연구에서 확보한 3차원 실험 데이터, 즉 전체 3차원 힘 분포 및 거시적 힘-변형률 반응과 뛰어난 일치를 보였다.
• 모델은 실험 스케일링 법칙 F ∝ f Z φ b 를 높은 정밀도로 재현하였으며, 표준 DEM보다 뚜렷한 편차를 보이지 않았다.
• MC-DEM 시뮬레이션은 해제 단계에서 약간의 곡률을 보이며, 표준 DEM보다 실험 데이터에 더 가까운 경향을 보였지만, 미미한 히스테리시스는 여전히 존재했다.
• MC-DEM에서 최대 상단 힘 F 는 실험적 Rtop 값과 작은 Fmin 오프셋 내에서 일치하였으며, 강력한 정량적 일치를 나타냈다.
• 모델은 입자 변형으로 인한 새로운 접촉 형성 현상을 성공적으로 재현하였으며, 이는 표준 DEM에서 누락된 핵심 메커니즘이었다.
• 경험적 비율 계수 γ 는 경계 조건 효과를 효과적으로 보상하여, 향후 모델에서 복잡한 구속 조건을 입자 수준 통계로 모델링할 수 있음을 시사한다.