[논문 리뷰] Numerical study of 3D vesicles under flow: discovery of new peculiar behaviors

Thierry Biben, Alexander Farutin외 1명|ArXiv.org|2009. 12. 23.

Blood properties and coagulation인용 수 4

한 줄 요약

이 연구는 간단한 비틀림 유동 조건에서의 3차원 수치 해석을 통해, 유동 속도가 증가함에 따라 진동-흡착 모드 영역이 크게 넓어지고, 뒤집힘(TB) 및 진동-흡착(VB) 진동 주기가 유동 속도에 거의 영향을 받지 않는 두 가지 새로운 행동을 규명하였다. 경계 적분 시뮬레이션과 분석 모델링을 통해 저자들은 큰 변형에서 고차 조화성분이 중요한 역할을 함을 입증하였으며, 이는 강한 유동 조건에서 고전적 Keller-Skalak 이론의 타당성을 도전한다.

ABSTRACT

The study of vesicles under flow, a model system for red blood cells (RBCs), is an essential step in understanding various intricate dynamics exhibited by RBCs in vivo and in vitro. Quantitative 3D analyses of vesicles under flow are presented. The regions of parameters to produce tumbling (TB), tank-treating (TT), vacillating-breathing (VB) {and even Kayaking (K)} modes are determined. New qualitative features are found: (i) a significant widening of the VB mode region in parameter space upon increasing shear rate $\dotγ$ and (ii) a striking robustness of period of TB and VB with $\dotγ$. Analytical support is also provided. These findings shed new light on the dynamics of blood flow.

연구 동기 및 목표

• 적혈구(RBCs)의 모델로 삼는 간단한 비틀림 유동 조건에서 3차원 적혈구 동역학을 체계적으로 연구하기 위해.
• 뒤집힘(TB), 탱크트레딩(TT), 진동-흡착(VB), 카이킹(K) 모드의 매개변수 영역을 식별하고 특성화하기 위해.
• 특히 진동-흡착(VB) 모드의 유동 속도 의존적 넓이 증가와 TB 및 VB 진동 주기가 유동 속도에 영향을 받지 않는 현상을 포함해 이전에 관측되지 않은 적혈구 동역학 행동을 밝혀내기 위해.
• 강한 유동 조건이나 큰 변형에서 Keller-Skalak(KS) 이론의 한계를 평가하기 위해.
• 관측된 수치 결과를 설명하는 분석 이론을 개발하고 검증하기 위해, 특히 형태 동역학에서 고차 조화성분의 역할을 설명하기 위해.

제안 방법

• 스토크스 방정식을 해결하기 위해 경계 적분(BI) 방법을 사용한 수치 시뮬레이션을 통해 간단한 비틀림 유동에서의 적혈구를 분석하였다.
• 막의 탄성 특성을 모델링하기 위해 헬프리치 굽힘력과 굽힘 모odulus κ, 면적 불가약성 제약 조건을 통합하였다.
• 유동 조건에서 적혈구 형태 변화를 특성화하기 위해 무차원 변형도 매개변수 $ C_{\bar{k}} $ 를 사용하였다.
• 특히 적혈구 주축의 각속도 $ \dot{\theta} $ 에서의 진동 특성에 대해 시뮬레이션 결과를 Keller-Skalak(KS) 이론과 비교하였다.
• 큰 변형에서 고차 조화성분(예: $ \cos(4\theta) $)의 발생을 설명하기 위해 펌진 이론을 적용하였다.
• 장기 시뮬레이션을 통해 일시적 동역학과 안정된 상태 모드로의 회복을 평가하였으며, 특히 TT-VB, TB-K 등의 모드 전이선 근처에서의 동역학을 분석하였다.

실험 결과

연구 질문

• RQ1 유동 속도 $ \dot{\gamma} $ 가 증가함에 따라 매개변수 공간에서 진동-흡착(VB) 모드 영역의 너비는 어떻게 변화하는가?
• RQ2 이전 시뮬레이션과는 달리, 뒤집힘(TB) 및 진동-흡착(VB) 모드의 진동 주기가 왜 유동 속도에 거의 영향을 받지 않는가?
• RQ3 강한 유동 또는 큰 적혈구 변형에서 Keller-Skalak 이론의 붕괴 정도는 어느 정도인가?
• RQ4 고차 조화성분(예: $ \cos(4\theta) $)은 크게 변형된 적혈구의 동역학에서 어떤 역할을 하는가?
• RQ5 모드 전이선 근처에서 적혈구가 안정된 상태 동역학으로 회복하는 데 걸리는 시간은 얼마이며, 이 회복 시간 척도는 무엇인가?

주요 결과

• 유동 속도 $ \dot{\gamma} $ 가 증가함에 따라 매개변수 공간에서 진동-흡착(VB) 모드가 발생하는 영역이 유의미하게 넓어지며, 이는 이전 이론이나 시뮬레이션으로는 기록되지 않은 행동이다.
• 뒤집힘(TB) 및 진동-흡착(VB) 모드의 진동 주기는 유동 속도 변화 범위가 10배에 이르는 동안 거의 일정하게 유지되며, 이는 $ \dot{\gamma} $ 에 독립적인 강건한 동역학적 행동을 나타낸다.
• Keller-Skalak(KS) 이론은 변형이 작거나 점성비 $ \lambda $ 가 클 경우에만 TB 및 TT 동역학을 정확히 기술하며, 높은 $ \dot{\gamma} $ 또는 중간 수준의 $ \lambda $ 에서는 심각한 이격 현상이 발생한다.
• 고차 조화성분, 특히 $ \cos(4\theta) $ 는 큰 변형도($ C_{\kappa} \geq 8 $)에서 각운동의 동역학에서 두드러지게 나타나며, 이는 TB-VB 또는 TB-K 분기 전이와 같은 복잡한 전이의 시작을 시사한다.
• TT-VB 전이선 근처에서 안정된 상태 동역학으로의 회복 시간은 $ 30 \eta_{in}R^3 / \kappa $ 수준이며, 약 125회의 진동에 해당하여 느린 일시적 영역임을 나타낸다.
• 본 논문에서 개발한 분석 이론은 수치적 관측 결과를 성공적으로 설명하며, 특히 고차 항의 발생과 진동 주기의 강건성에 대한 설명이 뚜렷하다.