두 단계 유체역학 모델로 본 활성 콜로이드 운동

초록

본 논문은 구형 마이크로스위머가 포함된 얇은 활성 콜로이드 층의 집단 운동을 재현하기 위해, 입자와 연속 유체를 각각 라그랑지안·오일러식으로 다루는 두 단계 유체역학 모델을 제시한다. 입자의 무작위 구동, 탄성 충돌, 그리고 유체와의 스톡스 마찰을 핵심 메커니즘으로 포함시키고, 각 메커니즘의 특성 시간 스케일을 정확히 해석한다. 모델은 브라운 운동, 집단 운동, 그리고 그 전이 상태를 재현하며, 미세스위머 속도가 임계값을 초과하면 유체와의 수소역학적 상호작용이 지배적이 되고, 속도가 더 증가하면 충돌 효과가 다시 브라운‑유사 상태로 전환되는 현상을 보여준다.

상세 분석

이 연구는 활성 콜로이드 시스템을 두 개의 상(phase)으로 분리해 각각의 동역학을 직접 해석한다는 점에서 기존의 단일‑상 모델과 차별화된다. 입자는 라그랑지안 입자(Lagrangian particle)로 취급되어, (1) 일정한 크기의 활성 힘 F_active 을 받으며, 회전 확산과 백색 잡음에 의해 방향이 지속적으로 변한다. (2) 입자 간 충돌은 완전 탄성으로 가정하고, 충돌 시 순간적인 속도 교환을 통해 에너지를 보존한다. (3) 입자는 주변 유체와 스톡스 마찰 F_drag = u_p / τ_St 으로 연결되며, τ_St = m_p / (3πμd_p) 로 정의된 스톡스 시간에 의해 입자와 유체 사이의 운동량 교환 속도가 결정된다.

유체는 연속체 방정식(비압축성 Navier‑Stokes)으로 기술되며, 입자에 의해 발생하는 체적 힘 f_b = (1/ρV)∑ m_p du_p/dt 가 소스 항으로 추가된다. 저마흐수( low‑Mach ) 흐름을 가정하고, 속도·압력 결합을 피하기 위해 포텐셜 H = |u|²/2 − p/ρ 를 도입해 모멘텀 방정식을 ∂u/∂t + F + ∇H = 0 형태로 변형한다. 이후 포아송 방정식 ∇²H = −∇·F 를 멀티그리드 방법으로 푸는 예측‑보정(predictor‑corrector) 알고리즘을 적용한다.

시간 스케일 해석이 핵심이다. 입자 움직임과 유체 반응 사이의 시간 비율을 τ = Δt_f/Δt_p 로 정의하고, τ = 1이면 입자와 유체가 동일한 시간 스케일로 상호작용한다. τ > 1이면 입자는 유체에 에너지를 전달하기 전에 여러 번 방향을 바꿀 수 있어, 충돌과 수소역학적 상호작용을 보다 정확히 포착한다. 저자들은 τ = 1, 10, 100에 대해 테스트 문제(원형 도메인 내에서 입자가 원형 궤도를 따라 움직이는 경우)를 수행했으며, τ가 클수록 입자가 목표 궤적을 유지하고, 유체에 뚜렷한 와류가 형성되는 것을 확인했다. 이는 수치적 시간 해상도가 물리적 현상 재현에 결정적임을 보여준다.

주요 결과는 매개변수 u_p (입자 속도)와 τ, 그리고 부피분율 n_p 에 대한 상도(phase diagram)이다. 세 가지 동역학 모드가 식별된다. (1) 브라운‑유사 모드: 입자는 거의 정지한 유체 속에서 무작위하게 움직이며, 집단 흐름이 거의 없으며, τ가 작을 때 나타난다. (2) 전이 모드: 초기에는 집단 흐름이 형성되지만, 시간이 지나면서 다시 브라운‑유사 상태로 전환된다. (3) 집단 모드: 입자와 유체가 강하게 결합해 지속적인 와류와 유향 흐름을 만든다. 특히 u_p 가 약 0.5 mm/s 이상일 때 수소역학적 상호작용이 급격히 강화되어 집단 모드가 나타나며, u_p 가 ~1.3 mm/s 이상을 초과하면 충돌 빈도가 증가해 다시 브라운‑유사 상태로 되돌아간다. 부피분율 n_p 가 0.2에서 0.6으로 증가하면 전이 구간이 넓어지고, 집단 모드가 유지되는 u_p 범위가 축소된다.

이 모델은 실험적 마이크로스위머(예: 박테리아, 마그네틱 나노입자 기반 스위머)에서 보고된 비선형 전이 현상을 정량적으로 재현한다는 점에서 의미가 크다. 또한, 두 단계 접근법은 입자‑유체 상호작용을 별도 소스 항 없이 자연스럽게 포함시켜, 향후 복잡한 경계 조건(예: 파이프, 혈관)이나 외부 구동(전기·자기장)과 결합한 시뮬레이션에 확장 가능성을 제공한다.