다중입도 전단증가 현상의 강성 전이와 임계 현상

초록

본 연구는 2차원 LF‑DEM 시뮬레이션을 이용해 비브라운 입자들의 다중입도 전단증가 현상을 조사한다. 펩플 게임 알고리즘으로 강성 클러스터를 분석한 결과, 포장률 ϕ가 임계값 ϕ_c(≈ shear‑jamming 전이 ϕ_J^μ)보다 클 때 강성 클러스터가 급격히 성장한다는 전이 현상을 발견하였다. 전이 지표 f_rig와 그 민감도 χ_rig는 2D 퍼콜레이션과 동일한 임계 지수 β, γ, ν를 보이며, 다중입도 시스템과 통계적으로 동등한 이분산 시스템 사이에 차이가 없음을 확인하였다.

상세 분석

이 논문은 고밀도 비브라운 입자들의 전단증가(Shear Thickening) 현상을 미세구조 수준에서 규명하고자, 2차원 Lubrication‑Flow Discrete Element Method(LF‑DEM)와 펩플 게임(pebble game) 알고리즘을 결합한 새로운 분석 프레임워크를 제시한다. 먼저, 입자 크기 분포를 다중입도(polydisperse)와 통계적으로 동등한 이분산(bidisperse) 두 경우로 설정하고, 비정상적인 고응력(σ≫σ₀) 조건에서 시뮬레이션을 수행한다. LF‑DEM은 저레놀즈, 저스트okes 수 영역에서 입자 간의 윤활(lubrication) 상호작용, 마찰 접촉(contact) 및 짧은 거리 전기반발(repulsion)을 포함해 실제 실험과 일치하는 전단‑점도 곡선을 재현한다.

핵심은 전단에 의해 형성되는 마찰 접촉 네트워크가 강성 클러스터(rigid cluster)로 전이되는 과정을 정량화한 것이다. 펩플 게임은 2D 평면에서 각 입자를 2개의 ‘펩’(degree of freedom)으로 모델링하고, 접촉을 ‘힌지’(constraint)로 취급해 전체 시스템의 자유도 감소를 계산한다. 이때, 전체 자유도가 충분히 제한되어 고정된 구조가 형성되면 해당 입자 집합을 ‘강성 클러스터’라 정의한다.

연구 결과, 포장률 ϕ가 특정 임계값 ϕ_c(ϕ_c<ϕ_J^μ)보다 작을 때는 강성 클러스터가 미미하거나 국소적으로만 존재한다. ϕ가 ϕ_c를 초과하면 f_rig≈(ϕ−ϕ_c)^β 형태로 급격히 증가하고, 그 변동성 χ_rig≈|ϕ−ϕ_c|^{−γ} 로 발산한다. 여기서 얻어진 β≈0.14, γ≈2.4는 2D 퍼콜레이션 이론의 알려진 값과 일치한다. 또한, 시스템 크기 L을 변화시켜 유한크기 스케일링을 수행한 결과, 상관길이 ξ∝|ϕ−ϕ_c|^{−ν} 로, ν≈1.33이라는 퍼콜레이션 차원을 재현하였다.

특히 주목할 점은 다중입도 시스템과 통계적으로 동등한 이분산 시스템 사이에 f_rig와 χ_rig의 전이 곡선이 거의 동일하다는 것이다. 이는 입자 크기 분포가 복잡해지더라도 마찰 접촉 네트워크의 토폴로지는 전체 부피분율과 평균 입자 크기만으로 결정된다는 암시를 제공한다. 또한, ϕ_c와 ϕ_J^μ이 입자 강성(k_n) 및 폴리디스퍼시티 지수와 비단조(non‑monotonic) 관계를 보이며, 강성 파라미터가 클수록 전이점이 낮은 ϕ에서 나타난다. 이는 입자 변형 가능성이 강성 클러스터 형성에 중요한 역할을 함을 시사한다.

결론적으로, 전단‑유도 마찰 네트워크의 강성 전이는 전통적인 전단‑점도 전이(CST→DST)와 shear‑jamming을 연결하는 비평형 임계 현상이며, 2D 퍼콜레이션과 동일한 보편적 임계 거동을 보인다. 이러한 통계적 동등성은 복잡한 다중입도 현상을 단순화된 이분산 모델로도 충분히 설명할 수 있음을 보여준다.