Reciprocal swimming in granular media: the role of jamming and swimmer inertia

기술 요약: 입상 매질에서의 왕복 수영: 재밍(Jamming)과 수영체의 관성의 역할

문제 정의
입상 매질(Granular media)에서의 이동은 매질이 유체와 같은 흐름과 고체와 같은 재밍(jamming) 특성을 모두 나타낼 수 있기 때문에 독특한 과제를 안겨준다. 많은 생물학적 및 로봇 수영체는 이러한 환경을 항해하기 위해 복잡한 비왕복적 보행(예: 파동 운동)을 활용하지만, 단순한 왕복 운동(예: 조개 모양의 날개짓)은 저레이놀즈수 뉴턴 유체에서 순 변위를 생성하지 못할 것으로 이론적으로 예측된다("조개 정리(Scallop Theorem)"). 본 논문은 조개 모양의 수영체가 입상 매질에서 왕복적인 날개짓을 통해 어떻게 순 추진력을 얻는지 조사하며, 특히 **재밍 유도 이력 현상(jamming-induced hysteresis)**과 수영체의 관성의 역할을 분리하여 규명한다.

연구 방법론
저자들은 MercuryDPM 패키지를 이용한 이산 요소법(Discrete Element Method, DEM) 시뮬레이션을 사용하여 입상 저장소에 잠긴 조개 모양의 수영체를 모델링하였다.

• 수영체 기하학적 구조: 수영체는 간격 $S$를 가진 허브를 중심으로 회전하는 두 개의 정사각형 날개로 구성된다. 날개는 주기 $T$ 동안 닫힌 각도($\theta_c = 20^\circ$)와 열린 각도($\theta_o = 80^\circ$) 사이를 왕복하며 진동한다.
• 매질 특성: 입상 매질은 밀도 $1050 \text{ kg/m}^3$, 평균 직경 $1.0 \text{ mm}$인 약 $510,000$개의 입자로 구성된다. 입자 간 상호작용에는 헤르츠(Hertzian) 탄성 접촉, 점성 소산, 쿨롱 마찰이 포함된다.
• 매개변수 연구: 연구는 수영체의 질량($m_s$), 날개짓 주기($T$), 날개 간격($S$), 입자 마찰 계수($\mu$)를 체계적으로 변화시킨다.
• 검증: 시뮬레이션 결과는 수영체의 변위 및 입자 운동학에 관한 기존 실험 데이터(Xiao et al., 2025)와 비교 검증되었다(X-선 CT 이미징을 통해 확인).
• 분석 도구: 저자들은 정체 구역(stagnant zones)과 저항을 정량화하기 위해 입자 운동학(속도장, 변형률 속도), 구조적 지표(배위수 $Z$), 힘 전달(접촉력, 강한 접촉 네트워크)을 분석한다.

주요 기여 및 결과

1. 재밍 유도 추진 (준정적 영역)
준정적 영역(낮은 수영체 관성, 느린 날개짓)에서 순 추진력은 정체 구역(수영체와 함께 움직이는 재밍된 입자 영역)의 역사 의존적 형성에 의해 발생한다.

• 비대칭 정체 구역: 열리는 스트로크 동안에는 날개 끝 근처에 "추진 정체 구역(PSZ)"이 형성된다. 닫히는 스트로크 동안에는 날개 허브 근처에 "저항 정체 구역(RSZ)"이 형성된다.
• 힘의 이력 현상: 완전히 열린 상태와 완전히 닫힌 상태 사이에서 강한 입자 접촉 수($N_c$)가 크게 차이 난다. 구체적으로, 닫히는 스트로크는 허브 근처에서 더 높은 밀도의 강한 힘 사슬(높은 $N_c$)을 생성하여, 전방으로 열릴 때 발생하는 저항보다 후방 운동에 대한 더 큰 저항을 만들어낸다.
• 선형 상관관계: 주기당 순 변위($\Delta y$)는 강한 접촉 수의 차이($\Delta N_c = N_c(\theta_c) - N_c(\theta_o)$)와 선형 관계를 보인다. 이는 대칭성 파괴가 재밍된 접촉 네트워크의 이력적 진화에 의해 발생함을 확인시켜 준다.
• 마찰 의존성: 마찰이 없는 시뮬레이션($\mu=0$)에서는 정체 구역이 형성되지 않으며 순 이동이 사라지는데, 이는 입자 간 마찰이 이 메커니즘에 필수적임을 확인해 준다.

2. 관성 유도 추진 (동적 영역)
동적 영역(높은 수영체 질량, 빠른 날개짓)에서는 수영체의 관성과 **코스팅 시간(coasting times, $T_c$)**의 차이에 의해 구동되는 두 번째 메커니즘이 나타난다.

• 코스팅 시간 정의: $T_c$는 움직이는 수영체가 매질의 항력에 의해 정지하는 데 걸리는 특성 시간으로 정의된다. 이는 수영체 질량과 항력의 비율에 비례한다 ($T_c \propto m_s / P_s L T$).
• 감속 비대칭성: 날개짓 주기 $T$가 코스팅 시간 $T_c$에 근접할 때, 수영체는 날개가 방향을 바꿀 때 즉각적으로 멈추지 않는다. 대신 "감속(retardation)" 과정을 거친다.
• 방향 역전: 투영 면적이 더 큰 열리는 스트로크 동안의 저항 항력은 닫히는 스트로크 동안보다 크다. 결과적으로, 수영체는 열리는 스트로크 이후보다 닫히는 스트로크 이후에 더 긴 시간 동안 더 멀리 코스팅한다. 이러한 감속의 비대칭성은 후방 방향으로의 순 변위를 초래한다.
• 무마찰 확인: 이 관성 메커니즘은 마찰이 없는 매질에서도 지속되며, 이는 재밍 메커니즘과 구별되는 특징이다.

3. 통합 스케일링 법칙
본 논문은 이 두 메커니즘이 공존하며 서로 경쟁함을 보여준다.

• 통합 모델: 순 변위는 다음과 같은 통합 스케일링 법칙으로 설명될 수 있다:
  $$\Delta y/L = \alpha \Delta N_c - \beta (T_c/T)^n$$
  여기서 첫 번째 항은 재밍 기여(전방 운동 선호)를 나타내고, 두 번째 항은 관성 기여(후방 운동 선호)를 나타낸다.
• 영역 전이: 낮은 $T_c/T$에서는 재밍 메커니즘이 지배한다(전방 운동). $T_c/T$가 증가함에 따라 관성 메커니즘이 커지며, 결국 재밍 효과를 압도하여 후방 이동으로의 전환을 일으킨다.

의의
본 논문은 왕복 수영체가 이동할 수 있게 하는, 시간 가역성을 깨뜨리는 두 가지 뚜렷한 물리적 메커니즘을 밝혀냈다고 주장한다.

• 메커니즘 1 (재밍): 마찰이 있는 재밍된 매질 내에서 접촉 네트워크의 역사 의존적 재편성이 비대칭 저항을 생성하여 전방 추진을 가능하게 함을 확립한다. 이는 강한 접촉력의 차이로 정량화된다.
• 메커니즘 2 (관성): 수영체의 관성이 매질의 항력 응답 시간과 유사할 때, 코스팅 시간의 비대칭성을 기반으로 한 두 번째 대칭성 파괴 메커니즘을 식별하여 후방 추진을 가능하게 함을 밝힌다.
• 이론적 함의: 본 연구는 뉴턴 유체에서는 성립하지만 입상 매질의 복잡성(탄소성 및 관성)은 순 변위를 위한 다수의 경로를 제공한다는 점을 보여줌으로써, 입상 매질에서의 조개 정리(Scallop Theorem)에 대한 이해를 확장한다.
• 설계 통찰: 결과는 수영체의 기하학적 구조(예: 관성이 있는 중간 몸체 추가)와 질량을 조절함으로써 입상 환경에서 이동 방향과 효율을 최적화할 수 있음을 시사한다.