입자 크기와 공극률을 지정한 확률적 입자 포장 기법

초록

본 논문은 입도곡선을 경험적 누적분포함수로 활용해 입자 크기를 생성하고, 목표 공극률을 만족하도록 입자를 무작위로 배치하는 방법을 제시한다. 입도곡선은 감마, 로그정규, 와이블, 하이퍼볼릭 분포의 혼합으로 모델링하고, 최적 적합 파라미터를 이용해 이론적 분포에서 샘플을 추출한다. 추출된 원형 입자는 Simple Sequential Inhibition(SSI) 포인트 프로세스로 겹치지 않게 제한된 영역에 배치된다. 결과는 이산요소법(DEM) 등 입자 기반 시뮬레이션의 현실적인 초기 조건을 제공한다.

상세 요약

이 연구는 입자 기반 시뮬레이션, 특히 Discrete Element Method(DEM)와 같은 수치 해석에 필요한 입자 구성을 효율적으로 생성하는 새로운 프레임워크를 제시한다. 핵심 아이디어는 현장 실험이나 실험실 측정으로 얻은 입도곡선을 경험적 누적분포함수(empirical CDF)로 변환하고, 이를 기반으로 입자 크기의 확률분포를 추정하는 것이다. 저자는 입도곡선을 단순히 히스토그램으로 처리하는 대신, 감마, 로그정규, 와이블, 하이퍼볼릭 등 네 가지 연속 확률분포를 후보로 삼아 최소제곱법과 최대우도법을 이용해 파라미터를 추정한다. 각 분포는 서로 다른 꼬리 특성을 가지고 있어, 실제 토양이나 점토와 같이 비대칭적이고 장측면이 긴 입도분포를 보다 정밀하게 모델링할 수 있다.

파라미터 추정 후, 추정된 이론분포를 이용해 입자 반지름을 샘플링한다. 여기서 중요한 점은 샘플링 과정이 균등분포의 혼합으로 구현된다는 점이다. 즉, 전체 구간을 여러 구간으로 나누고 각 구간마다 균등분포를 적용해 복합적인 입도분포를 재현한다. 이 방식은 복잡한 분포를 직접 샘플링하는 데 따르는 계산 비용을 크게 낮추면서도, 원본 입도곡선과의 차이를 최소화한다.

다음 단계는 목표 공극률(porosity)을 만족하도록 입자 수를 결정하는 것이다. 공극률은 전체 영역 면적 대비 입자들이 차지하는 면적의 비율로 정의되며, 목표값에 도달하기 위해서는 적절한 입자 수와 크기 조합이 필요하다. 저자는 전체 영역 면적 A와 목표 공극률 φ를 이용해 필요한 입자 총 면적 Σπr_i^2 = (1−φ)A 를 만족하도록 입자 수 N을 조정한다. 이때 입자 크기의 변동성을 고려해 평균 입자 면적과 분산을 반영한 통계적 접근을 사용한다.

입자 배치는 Simple Sequential Inhibition(SSI) 포인트 프로세스를 적용한다. SSI는 새로운 입자를 무작위 위치에 제안하고, 기존 입자와 겹치지 않을 경우에만 수용하는 절차적 알고리즘이다. 이 과정은 입자 간 겹침을 방지하면서도 포장 밀도를 높이는 장점을 가진다. 저자는 SSI의 수렴성을 검증하기 위해 배치 시도 횟수와 성공률을 모니터링하고, 필요 시 재시도 횟수를 제한하여 계산 시간을 제어한다. 또한, 경계 조건을 고려해 입자가 영역 외부로 넘어가지 않도록 반사 또는 재시도 전략을 적용한다.

실험에서는 다양한 모래와 점토 시료를 대상으로 입도곡선 데이터를 수집하고, 각각에 대해 네 가지 분포 중 최적 모델을 선택하였다. 감마분포는 비교적 균일한 입도분포를, 로그정규와 와이블은 비대칭적인 꼬리를, 하이퍼볼릭은 극단적인 대형 입자 비중을 잘 포착한다. 각 시료에 대해 추정된 파라미터를 바탕으로 10,000개 이상의 입자를 생성하고, 목표 공극률 0.35~0.45 범위 내에서 SSI를 통해 포장하였다. 결과적으로 생성된 입자 집합은 시각적으로 실제 토양 구조와 유사했으며, 입자 간 거리와 접촉 네트워크 통계량도 실험 데이터와 일치하는 수준을 보였다.

이 논문의 주요 기여는 (1) 입도곡선을 확률적 모델링하여 이론분포와 경험분포를 연결한 방법, (2) 목표 공극률을 만족하도록 입자 수를 통계적으로 추정하는 절차, (3) SSI를 활용한 효율적인 비중첩 입자 배치 알고리즘이다. 이러한 접근은 기존에 수작업으로 수행되던 입자 초기화 과정을 자동화하고, 다양한 입도특성을 가진 재료에 대해 일반화 가능한 프레임워크를 제공한다. 향후 연구에서는 3차원 구형 입자 확장, 비구형 입자 형태 모델링, 그리고 동적 포장(예: 압축·진동) 과정을 통합하는 방향으로 발전시킬 수 있다.