무질서 매질에서 미세수영자들의 지속적 무작위 운동과 트래핑

초록

이 연구는 2차원 무작위 다공성 매질 속에서 스퀴머 모델을 이용해 미세수영자들의 전수체 유체역학을 풀어낸다. 활동성, 유체 상호작용, 그리고 장애물의 무질서가 결합하면 결정론적이면서도 확산적인 이동이 나타나며, 푸셔와 풀러에 따라 서로 다른 정적·동적 트래핑 메커니즘이 발생한다. 짧은 거리의 반발 포텐셜을 조절하면 트래핑 확률과 지속 시간이 크게 변하며, 이는 근거리 상호작용의 민감성을 보여준다.

상세 분석

본 논문은 저레옹즈 수( Reynolds number ≪1 ) 조건에서 2차원 평면에 배치된 원형 장애물 군집 속을 이동하는 스퀴머(squirmer)를 전산 유한요소법(FEM)으로 정밀 시뮬레이션하였다. 스퀴머는 반경 a인 원판 형태이며, 표면 전단속도 u_S = B₁(1+β p·n)(nn−I)·p 로 정의된다. 여기서 B₁은 첫 번째 스퀴밍 모드로 자유 공간에서 평균 속도 U = B₁/2 를 결정하고, β는 스트레스릿 강도 파라미터로 β<0은 푸셔, β>0은 풀러, β=0은 중성 스퀴머를 의미한다. 장애물은 반경 R=4a인 원형 디스크 N개가 무작위로 배치된 형태이며, 전체 부피분율 ϕ는 0.15(희박)와 0.45(조밀) 두 값을 사용한다. 스퀴머와 장애물 사이에는 최소 거리 h_cut = δ·a 를 보장하는 짧은 거리 반발 포텐셜이 적용되며, δ는 1/20과 1/4 두 경우를 탐색한다. 이러한 설정은 효과적인 장애물 반경 R_eff = R + a + h_cut 를 만들고, ϕ=0.45일 때 실제 유효 부피분율 ϕ_eff≈0.69 로 거의 퍼콜레이션 임계점에 도달한다는 점이 핵심이다.

시뮬레이션은 주기적 경계조건을 적용해 무한히 큰 매질을 모사하고, 각 파라미터 조합에 대해 128개의 궤적을 4개의 매질 샘플링으로 평균하였다. 스퀴머의 위치와 방향은 힘·토크 자유 조건(∮_S n·σ dS =0, ∮_S r×n·σ dS =0) 하에 Stokes 방정식(μ∇²u = ∇p, ∇·u=0)을 풀어 얻는다.

주요 결과는 크게 두 가지 동역학적 상태로 구분된다. 첫째, “탐색 상태”에서는 스퀴머가 장애물 사이를 자유롭게 이동하며 평균제곱변위(MSD) ⟨|Δr(t)|²⟩가 초기에는 U²t²(영구적 직진)에서 t≈τ = (R+a)/U 이후 선형(t) 혹은 초과선형(t^α, α>1)으로 전이한다. 특히, 희박 매질(ϕ=0.15)에서 풀러(β>0)는 슬라이딩 스캔이 적어 방향 재조정이 적고, MSD가 t²에 가까운 초확산을 보인다. 반면 푸셔(β<0)는 장애물 표면에 접선 정렬되는 경향이 있어 약간의 서브확산을 나타낸다. 두 번째, “트래핑 상태”는 정적 트래핑(장애물 코너에 고정)과 동적 트래핑(두 개 혹은 그 이상의 장애물 사이에서 준주기적 궤도)으로 나뉜다. 동적 트래핑은 궤도가 완전히 닫히지 않아 매 순간 약간씩 위치가 이동하며, 충분히 긴 시간 후 탈출(‘hopping’)이 가능하다. 트래핑 확률 P(β)는 β의 절대값이 클수록 증가하지만, δ=1/4와 같이 반발 범위가 넓어지면 푸셔의 트래핑이 억제되고 풀러의 트래핑이 강화되는 비대칭성이 역전된다. 이는 근거리 반발력이 유동장과의 상호작용을 차폐하거나 강화함으로써 스트레스릿 효과를 조절한다는 것을 의미한다.

생존 확률 S(t) 분석에서는 트래핑이 일어나지 않은 비율이 시간에 따라 지수적 감소를 보이며, 평균 트래핑 시간 ⟨T⟩는 β와 ϕ, δ에 따라 크게 변한다. 특히, 강한 푸셔(β=−8)는 조밀 매질에서 거의 즉시 정적 트래핑에 빠지는 반면, 중성(β=0)과 약한 풀러(β=2)는 장시간 탐색 후에야 트래핑에 도달한다. 동적 트래핑 궤도는 반경 ρ≈a~2a 수준으로 측정되며, 이는 스퀴머가 장애물 사이에서 제한된 영역 내에서 지속적으로 회전한다는 것을 보여준다.

결론적으로, 이 연구는 난류가 없는 결정론적 유체-구조 상호작용만으로도 비평형 확산과 트래핑 현상이 자연스럽게 발생함을 입증한다. 푸셔·풀러 비대칭성, 장애물 밀도, 그리고 근거리 반발 포텐셜은 모두 미세수영자들의 장거리 전파와 지역화에 결정적인 역할을 하며, 이러한 메커니즘은 자연계 미생물의 서식지 적응 및 인공 마이크로로봇 설계에 직접적인 시사점을 제공한다.