비등방성 입자와 액체계면을 위한 라티스 볼츠만 시뮬레이션

초록

본 논문은 라티스 볼츠만(LB) 방법과 분자동역학(MD)을 결합해 다중성분 유체와 비구형(타원형) 입자의 상호작용을 시뮬레이션한다. IBM Blue Gene/P 슈퍼컴퓨터에서 수십만 코어를 활용한 강력한 스케일링을 구현하고, 단일 입자의 계면 흡착 및 Pickering 에멀전·bijel 형성 과정을 검증한다.

상세 분석

이 연구는 라티스 볼츠만(LB) 방법을 기반으로 Shan‑Chen 다중성분 모델을 사용해 두 개 이상의 유체를 동시에 다루면서, 입자와 유체 사이의 경계조건을 ‘moving bounce‑back’ 방식으로 구현한다. 입자는 격자 위에 이산화된 강체로 취급되며, 입자‑유체 상호작용을 보존하기 위해 충격량 보정(force C)과 추가 힘 F(t)을 도입한다. 구형 입자에 적용되던 Hertz 강체 포텐셜을 일반화해 타원형 입자(주축 Rₚ, 부축 Rₒ)의 방향성 (ô)과 거리 (rᵢⱼ)를 포함하는 식(7‑8)으로 확장했으며, 근접 입자 간 윤활(lubrication) 보정을 식(9)로 추가해 격자 해상도가 충분치 않을 때 발생하는 인공적인 마찰을 보정한다.

접촉각 Θ는 입자 색상 Δρ (밀도 차)으로 조절되며, 실험적으로 Θ와 Δρ 사이에 선형 관계 Θ = 243.2·Δρ + 90°가 관찰된다. 이는 입자 표면 친수·소수성을 자유롭게 설정할 수 있게 해, 실제 클레이 입자와 같은 비등방성 입자의 계면 친화성을 정밀히 모델링한다.

병렬 구현 측면에서는 3‑차원 도메인 분할을 이용해 LB와 MD 연산을 동일한 MPI 프로세스에 매핑했으며, 네트워크 토폴로지를 고려해 MPI_Cart_create와 수동 도메인 매핑을 적용해 통신 병목을 최소화했다. 최적화 전후의 스케일링 결과는 1024 × 2048 격자와 4 × 10⁸ 입자를 262 144 코어에서 8.97 × 10⁹ LUPS(Lattice Updates per Second)를 달성, 코어당 3.42 × 10⁴ LUPS, 전체 효율 93 %에 도달했음을 보여준다. 특히 Shan‑Chen 힘 계산이 전체 시간의 41.8 %를 차지하고, 입자‑유체 결합이 29.6 %를 차지한다는 점에서 향후 최적화 여지가 명확히 드러난다.

물리적 결과로는, 종축 비율 m = 2인 타원형 입자가 접촉각 Θ = 90°인 경우, 초기 기울기 φ = 0°에서 시작해 인터페이스에 접근하면서 z(중심‑계면 거리)와 φ가 동시에 변하는 복합적인 흡착 궤적을 보였다. 안정된 평형은 주축이 인터페이스에 수직(φ = π/2)인 상태이며, 메타안정 상태는 φ = 0°에 존재한다. 또한 입자 농도와 접촉각을 조절함으로써 Pickering 에멀전에서 bijel로의 전이 현상을 재현했으며, 높은 종축비(>2)의 경우 적은 입자 수만으로도 인터페이스를 안정화시키는 ‘팁‑투‑팁’ 배열이 자연스럽게 형성되는 것을 확인했다. 이러한 결과는 실험적 관찰과 일치하며, 비등방성 입자의 형태가 계면 변형 및 캡illary 상호작용에 미치는 영향을 정량적으로 설명한다.

전반적으로, 본 논문은 대규모 병렬 LB‑MD 프레임워크를 통해 비구형 입자의 복잡한 계면 동역학을 정밀히 시뮬레이션할 수 있음을 입증하고, 고성능 컴퓨팅 환경에서 멀티스케일 콜로이드 시스템을 연구하는 데 필요한 알고리즘적·구조적 설계 원칙을 제시한다.