단일 세포 해상도로 보는 원통형 및 협착 혈관 내 혈류 시뮬레이션

초록

본 연구는 혈액을 연성 입자 서스펜션으로 모델링하여, 100 µm 이하 직경의 혈관에서 혈액 세포를 단일 입자 수준으로 해석한다. 저혈소판(hematocrit < 30 %) 조건에서 혈관 직경에 따른 겉보기 점도 변화를 실험 데이터와 일치하도록 재현했으며, 특히 작은 혈관에서 관찰되는 Fahraeus‑Lindqvist 효과를 성공적으로 포착하였다.

상세 분석

이 논문은 혈액 흐름을 입자 기반으로 기술하는 새로운 모델을 제시한다. 기존의 연속체 모델은 혈구를 평균화된 점성 유체로 취급해 미세혈관에서의 비뉴턴 거동을 놓치기 쉽다. 반면, 전통적인 입자 시뮬레이션은 계산 비용이 급격히 증가해 실제 혈관 규모(수십~수백 µm)와 충분한 세포 수를 동시에 다루기 어렵다. 저자들은 이러한 딜레마를 해소하기 위해 ‘단순하면서도 효율적인’ 입자 모델을 설계했으며, 핵심 아이디어는 혈구를 변형 가능한 ‘소프트 입자’로 표현하고, 유체와의 상호작용을 격자볼츠만법(Lattice‑Boltzmann Method, LBM)으로 처리하는 것이다.

모델의 주요 구성 요소는 다음과 같다. 첫째, 혈구는 겉면 장력과 내부 압력 차이를 통해 변형을 제어하는 ‘스프링‑다이얼’ 메커니즘을 갖는다. 둘째, 혈구와 혈관벽 사이의 충돌은 반발력과 마찰을 포함한 포텐셜 함수로 구현되어, 세포가 혈관 내벽에 붙거나 슬립하는 현상을 자연스럽게 재현한다. 셋째, 유체 흐름은 3차원 D3Q19 격자볼츠만 스킴으로 계산되며, 입자와 유체 사이의 힘 교환은 ‘이중 분산’(immersed boundary) 기법을 통해 수행한다. 이러한 설계는 입자 수가 수천 개 수준에서도 시간당 수천 스텝을 실시간에 가깝게 수행할 수 있게 한다.

시뮬레이션 설정은 직경 20 µm, 50 µm, 100 µm인 원통형 혈관과, 중앙에 30 % 축소된 협착 구간을 가진 혈관을 포함한다. 혈소판 농도는 10 %~30 % 범위에서 변화를 주었으며, 입구에서 일정한 압력 구배를 가해 전형적인 파스칼 수치를 재현했다. 결과적으로, 혈관 직경이 감소함에 따라 겉보기 점도가 비선형적으로 감소하는 Fahraeus‑Lindqvist 효과가 관찰되었다. 특히 20 µm 이하의 초소관에서는 점도가 급격히 상승하는데, 이는 세포 간 충돌과 혈관벽 근접 효과가 지배적임을 시사한다. 협착 구간에서는 유속이 급증하고, 세포가 압축·정렬되는 현상이 뚜렷이 나타났으며, 이는 혈전 형성 위험을 예측하는 데 중요한 물리적 메커니즘으로 해석될 수 있다.

모델 검증을 위해 실험적 점도 데이터(Pries et al., 1992)와 비교했으며, 30 % 이하 혈소판 농도에서 평균 오차가 5 % 미만으로 매우 높은 일치도를 보였다. 또한, 기존의 전통적인 입자 기반 시뮬레이션과 비교했을 때 계산 시간은 10배 이상 단축되었으며, 메모리 사용량도 크게 감소했다. 한계점으로는 높은 혈소판 농도(>30 %)에서 세포 간 상호작용 모델이 단순화되어 있어 점도 과소평가가 발생할 수 있다는 점과, 혈관벽을 완전한 탄성체가 아닌 강체로 가정했기 때문에 실제 혈관의 비선형 탄성 효과를 완전히 포착하지 못한다는 점을 들 수 있다. 향후 연구에서는 혈관벽의 비선형 탄성 모델을 도입하고, 혈소판 활성화 및 응집 메커니즘을 추가함으로써 보다 정밀한 혈전 형성 시뮬레이션을 목표로 할 수 있다.