플럭투에이션 기반 전유체 입자 시뮬레이션 FIMAT 동역학으로 보는 브라운 운동

초록

본 논문은 전체 영역을 변동 유체로 가정하고, 강체 입자 영역에 강체 운동 제약을 부여하는 FIMAT(F​luctuating Immersed Material) 동역학을 제시한다. 이 방법으로 장시간 확산 거동과 단시간 동역학을 모두 정확히 재현하며, 기존 분자동역학 대비 효율적인 브라운 입자 시뮬레이션이 가능함을 검증한다.

상세 요약

FIMAT 동역학은 변동 수리학(fluctuating hydrodynamics)의 기본 방정식에 열적 난류(stochastic stress)를 포함하고, 이를 전역 유체 영역에 적용한다는 점에서 기존 Immersed Boundary Method(IBM)와 차별화된다. 강체 입자 내부에서는 속도와 회전 운동을 강제적으로 일정하게 유지하도록 Lagrange multiplier 형태의 제약 방정식을 도입한다. 이 제약은 유체-입자 경계에서의 연속성 및 운동량 보존을 엄격히 만족시키며, 입자 표면에서 발생하는 난류 응력은 입자 전체에 균등하게 전달된다. 시간 통합은 두 단계로 나뉜다. 첫 번째는 저주파(장시간) 확산 현상을 포착하기 위해 무한소 시간 간격에서 평균화된 Stokes 방정식을 풀어 입자 위치와 회전각을 업데이트한다. 두 번째는 고주파(단시간) 진동을 재현하기 위해 직접적인 Navier‑Stokes‑Landau‑Lifshitz 방정식을 시간 중심 차분법으로 해결한다. 이중 시간 스케일링은 베르누이‑스미스(β‑S) 정리와 플럭투에이션‑디시피에이션 균형을 보장한다. 수치 실험에서는 구형 입자와 비구형(타원형, 다면체) 입자를 대상으로 확산계수, 전이 확산계수, 그리고 자기상관 함수 등을 측정했으며, 결과는 이론적 스토크스‑에인슈타인 관계와 기존 분자동역학 시뮬레이션과의 비교에서 오차가 5 % 이하임을 보였다. 또한, 격자 해상도와 시간 단계의 수렴성을 검증함으로써 제안된 방법이 수치적 안정성과 정확성을 동시에 확보함을 입증하였다. 이러한 접근은 입자 간 수소 결합이나 전기적 상호작용 등 복잡한 미세구조를 포함하지 않지만, 유체 매질의 열적 요동을 정확히 반영함으로써 마이크로·나노 입자의 브라운 운동을 효율적으로 모델링한다는 점에서 실용적 가치를 가진다.