Pore-scale distribution and transport of active particles in a two-dimensional lattice

기술 요약: 2차원 격자 내 활성 입자의 기공 규모 분포 및 수송

문제 정의
박테리아 및 합성 활성 콜로이드와 같은 운동성 미세 유영체(microswimmers)는 고체 장애물과 불균일한 전단 흐름이 존재하는 복잡한 다공성 환경에 자주 서식한다. 단순한 채널 구조와 같은 기하학적 구조에서는 표면 축적(surface accumulation)과 전단에 대한 반응(rheotaxis)이 잘 알려져 있으나, 구조화된 다공성 매질에서의 수송 메커니즘은 여전히 덜 이해되어 있다. 구체적으로, 입자의 활성도(activity), 매질의 공극률(porosity), 다분산성(polydispersity), 그리고 배경 흐름의 세기와 같은 지배 매개변수들이 질서 있는 격자 내 활성 현탁액의 시간적 진화와 정상 상태 통계에 어떻게 집합적으로 영향을 미치는지에 대한 이해에 공백이 존재한다. 또한, 이러한 복잡한 전단 환경에서 상류(upstream) 및 하류(downstream) 유영 체제 사이의 전이 과정 중 극성장(polarization field)의 위상 결함(topological defects)이 수행하는 역할에 대한 체계적인 조사가 필요하다.

방법론
저자들은 2차원 정사각형 기둥 격자 내에 존재하는 희박한 활성 브라운 입자(Active Brownian Particles, ABPs)의 이중 모델링 접근법을 사용하여 연구를 수행한다.

1. 미시적 모델링: 입자는 속도 $v$와 회전 확산 $D_R$을 갖는 가늘고 긴 점 크기의 비상호작용 개체로 모델링된다. 이들의 역학은 병진 및 회전 노이즈를 포함하는 과감쇠 랑주뱅 방정식(overdamped Langevin equations)에 의해 지배된다.
2. 연속체 모델링: 확률 분포 함수를 이용한 포커-플랑크 방정식(Fokker-Planck equation)을 사용하여 거시적 기술을 도출한다. 이는 밀도($n$) 및 극성($m$) 장의 진화를 풀기 위해 2-모멘트 근사(네마틱 질서 $Q=0$ 무시)를 사용하여 폐쇄된다.
3. 수치 시뮬레이션:
   • 단일 기둥: 기준이 되는 축적 및 극성 거동을 확립하기 위해, 무한 도메인 내 고립된 기둥에 대한 연속체 운동론적 방정식을 수치적으로 해결한다(Method of Lines).
   • 주기적 격자: 다수의 입자를 포함하는 이중 주기적 정사각형 격자에서 브라운 역학(Brownian Dynamics, BD) 시뮬레이션을 수행한다. 이 접근법은 연속체 폐쇄가 깨질 수 있는 복잡한 기하학적 구조와 다수의 기둥을 분석할 수 있게 한다.
   • 흐름 조건: 배경 압력 구동 흐름을 부과하여, 기본 해법 방법(Method of Fundamental Solutions)으로 계산된 2차원 스토크스 흐름장을 생성한다. 입자의 역학은 흐름 이류(advection)와 전단 내 방향성을 위한 제프리 방정식(Jeffery's equation)을 포함하도록 수정된다.

주요 무차원 매개변수는 스윔 페클레 수($Pe_s = v/aD_R$), 흐름 페클레 수($Pe_f = u_\infty/aD_R$), 그리고 매질 공극률($\epsilon$)이다.

주요 기여 및 결과

• 정적 매질에서의 표면 축적: 흐름이 없는 경우, 자기 추진력은 입자들이 기둥 표면에 축적되도록 유도하여, 증가된 수 밀도와 방사형 극성(표면을 향해 정렬하는 입자들)을 가진 경계층을 형성한다. 시스템은 초기 균일 상태에서 지수적 성장을 거쳐 정상 상태로 진화한다.
• 과도적 진동: 구속된 격자($\epsilon < 1$)에서 정상 상태로의 접근은 표면 수 밀도의 과도하고 감쇠된 진동에 의해 특징지어진다. 이는 기둥 사이의 탄도 운동이 생성하는 결핍 파동(depletion waves)이 전파되고 상호작용하기 때문이다. 이러한 진동의 주파수는 도메인 길이에 반비례하여 스케일링된다.
• 공극률 및 활성도의 영향:
  • 낮은 공극률(높은 구속)은 이동 거리가 짧아지기 때문에 평형 시간을 단축시키지만, 강한 활성도 조건에서 정상 상태 표면 축적 비율($\eta = \bar{n}_A/\bar{n}_B$)을 감소시킨다.
  • 높은 활성도($Pe_s$)는 정상 상태 축적을 증가시키고 구배를 날카롭게 만든다.
  • 다분산 격자(두 가지 크기의 기둥)에서, 작은 기둥은 큰 기둥으로부터의 결핍 효과로 인해 정상 상태 밀도는 낮지만 초기 축적 속도는 더 빠르다. 다분산 매질은 동일한 공극률을 가진 단분산 매질에 비해 유체 벌크로부터 유영체를 필터링하는 데 더 효과적인 것으로 나타났다.
• 흐름 유도 역학 및 위상 결함:
  • 레오탁시스(Rheotaxis): 배경 흐름이 존재할 때, 입자들은 기둥 표면 근처에서 양의 레오탁시스(상류 유영)를 보이며, 기둥의 후류(wake)에 우선적으로 축적되는 경향을 보인다.
  • 유영 전이: 시스템은 낮은 흐름 세기($Pe_f$)에서의 전역적 상류 극성에서 높은 $Pe_f$에서의 상류 및 하류 유영 영역의 공존으로 전이된다.
  • 위상 결함: 이 전이는 극성장의 순수 운동학적 위상 결함(kinematic topological defects)의 핵 생성에 의해 매개된다. 중간 정도의 $Pe_f$에서, 나선형 결함(+1 전하)과 쌍곡선 결함(-1 전하) 쌍이 출현한다. 이 결함들은 격자 프레임 내에서 정지해 있으며, 유체 역학적 불안정성에 의해 구동되는 활성 극성 유체의 결함과는 구별되는 2차원 전단 흐름의 특성에서 기인한다.
  • 스윔 힘(Swim Force): 입자의 불균일한 분포는 기둥에 순 스율 힘을 발생시킨다. 이 힘은 중간 정도의 흐름 세기에서 정점에 도달한 후, 이류가 표면을 고갈시킴에 따라 감소한다.

의의
본 논문은 복잡한 기하학적 구조에서 활성 수송을 지배하는 물리적 메커니즘을 격리하여 연구할 수 있는 통제된 프레임워크를 제공한다. 연구는 활성 콜로이드의 수송이 자기 추진과 복잡한 전단 프로파일의 상호작용에 의해 구동됨으로써 수동적 수송과 근본적으로 다르다는 것을 보여준다. 주요 발견은 운동학적 위상 결함이 상류 및 하류 유영 체제 사이의 전이를 촉진하는 메커니즘이라는 것이며, 이는 격자의 2차원 전단 환경에 특유한 현상이다. 구조화된 격자는 구조화된 미세 유체 환경에서의 수송을 이해하기 위한 모델 시스템 역할을 하며, 구속, 활성도, 흐름 세기가 활성 물질의 공간적 분포와 집합적 행동을 어떻게 결정하는지에 대한 통찰을 제공한다. 본 연구는 희박하고 비상호작용하는 현탁액에서도 복잡한 전단 흐름이 결함 핵 생성 및 과도적 진동을 포함한 풍부한 역학적 거동을 유도할 수 있음을 강조하며, 이는 생물학적 식생(colonization) 및 다공성 환경에서의 합성 수송을 이해하는 데 중요하다.