A Modified Suspension-Balance Model for Deformable Particle Suspensions: Application to Blood Flows with Cell-Free Layer

이 논문은 미세혈관 채널 내 혈류를 효율적으로 시뮬레이션하기 위해 유체역학적 양력을 포함하는 수정된 현탁액-평형 모델을 제 제안하며, 세포 결핍층 형성, 헤마토크릿 및 속도 프로파일, 그리고 파라우스 및 파라우스-린드퀴스트 효과와 같은 주요 현상들을 성공적으로 포착한다.

핵심

바쁜 고속도로를 상상해 보세요. 자동차들은 사실 작고 말랑말랑한 빨간색 혈구(적혈구)이고, 도로는 미세한 혈관입니다. 일반적인 고속도로에서는 교통량이 고르게 퍼져 있을 것이라고 예상할 수 있습니다. 하지만 우리 몸속의 이 "자동차"들은 이상한 습성이 있습니다. 벽 근처에 있는 것을 싫어합니다. 대신 도로 한가운데로 모여드는 것을 선호하며, 길가 옆에는 텅 빈 깨끗한 차선을 남겨둡니다.

이 빈 차선을 **세포 결핍층(Cell-Free Layer, CFL)**이라고 부릅니다. 이는 혈액이 더 빠르게, 그리고 마찰을 줄이며 이동하는 데 도움을 주는 매우 중요한 특징입니다.

문제점: 기존 지도는 회전 구간을 놓치고 있었다

과학자들은 혈액이 어떻게 흐르는지 시뮬레이션하기 위해 수년 동안 컴퓨터 모델을 구축해 왔습니다. 그들은 이를 "현탁액 균형 모델(Suspension Balance Model, SBM)"이라고 부릅니다. 이 모델을 교통 시뮬레이션 소프트웨어라고 생각해 보세요.

기존 버전의 소프트웨어는 세포들이 서로 부딪히는 방식 덕분에 자동차들이 도로 중앙으로 이동한다는 점은 잘 예측했습니다. 하지만 왜 자동차들이 그토록 간절하게 벽을 떠나려 하는지는 설명하지 못했습니다. 즉, 세포 결핍층(CFL)이라는 깨끗한 구역을 만들어내지 못했습니다. 그것은 마치 자동차가 움직이고 있다는 것은 알지만, 자동차가 능동적으로 연석을 피하고 있다는 사실은 모르는 GPS와 같았습니다.

해결책: 새로운 "밀기" 버튼

Hugo Castillo-Sánchez와 Leonardo Liu가 이끄는 저자들은 이 소프트웨어를 수정하기로 했습니다. 그들은 적혈구가 말랑말랑하기(변형 가능하기) 때문에, 벽에 너무 가까워지면 특별한 종류의 보이지 않는 힘을 생성한다는 사실을 깨달았습니다.

그들은 이를 **양력(Lift Force)**이라고 부릅니다.

• 비유: 여러분이 수영장에서 벽 근처에서 수영하고 있다고 상상해 보세요. 여러분이 움직임에 따라, 물이 여러분을 벽으로부터 약간 밀어냅니다. 적혈구의 경우, 이 "밀어내는 힘"은 훨씬 강력합니다. 왜냐하면 그들은 말랑말랑하며 벽을 비집고 지나갈 때 모양이 변하기 때문입니다.
• 수정 사항: 연구팀은 이 "양력"을 컴퓨터 모델에 추가했습니다. 그들은 **수정된 현탁액 균형 모델(Modified Suspension Balance Model, MSBM)**을 만들었습니다. 이제 소프트웨어는 단순히 자동차를 관찰하는 것에 그치지 않고, 마치 물이 수영 선수를 밀어내는 것처럼 자동차를 벽으로부터 능동적으로 밀어냅니다.

시뮬레이션을 실행했을 때 어떤 일이 일어났나요?

컴퓨터에서 이 새로운 "양력" 기능을 켰을 때, 결과는 극적으로 변했습니다.

1. 빈 차선이 나타났습니다: 시뮬레이션은 실제 현상에서 보이는 것과 같은 벽 근처의 깨끗한 구역(CFL)을 성공적으로 만들어냈습니다.
2. 중앙의 교통 체증: 적혈구들이 중앙에 쌓이면서 밀도가 높은 핵(core)을 형성했습니다.
3. 흐름의 형태: 세포들이 가운데에 뭉쳐 있고 가장자리는 비어 있기 때문에, 혈액은 일반적인 강물처럼 매끄러운 곡선 형태로 흐르지 않았습니다. 대신, 평평한 윗부분을 가진 단단한 플러그(plug)나 피스톤처럼 흘렀습니다. 이는 실제 미세 혈관에서 일어나는 현상과 정확히 일치합니다.

새 모델의 테스트

연구팀은 단순히 추측만 한 것이 아니라, 실제 데이터 및 다른 복잡한 시뮬레이션들과 모델을 비교 테스트했습니다.

• 시간 여행: 그들은 "빈 차선"이 시간이 지남에 따라 어떻게 형성되는지 관찰했습니다. 처음에는 세포들이 곳곳에 퍼져 있었으나, "양력"이 세포들을 벽으로부터 서서히 밀어내면서 차선이 깨끗해졌습니다. 이는 고속 카메라 실험에서 보이는 속도 및 행동 양상과 일치했습니다.
• "파레우스 효과(Fåhræus Effect)": 이것은 단순한 관찰을 일컫는 멋진 용어입니다. 혈액은 아주 작은 관을 통과할 때 예상보다 더 빠르게 흐르며, 중앙의 세포 농도는 출구에서의 농도와 다릅니다. 그들의 새로운 모델은 이를 완벽하게 예측했습니다.
• "파레우스-린드퀴스트 효과(Fåhræus-Lindqvist Effect)": 이것은 혈액이 매우 작은 관을 통해 흐를 때 더 "묽게"(덜 끈적하게) 변한다는 관찰 결과입니다. 그들의 모델 역시 이를 포착하여, 벽 근처의 빈 차선이 마찰을 줄여 혈액 흐름을 더 쉽게 만든다는 것을 보여주었습니다.

핵심 요약

이 논문은 자신들의 컴퓨터 모델에 "양력"이라는 단순한 "밀기"를 추가함으로써, 이제 혈액이 미세한 혈관에서 어떻게 행동하는지 정확하게 시뮬레이션할 수 있다고 주장합니다.

• 수행하는 기능: 세포 결핍층의 형성, 플러그 형태의 흐름, 그리고 우리 몸에서 혈액 흐름을 효율적으로 만드는 유명한 "파레우스" 효과들을 포착합니다.
• 아직 하지 못하는 것(현재): 저자들은 관이 매우 큰 경우(40 마이크로미터 이상) 모델이 세포를 약간 너무 과하게 밀어낸다는 점을 인정했습니다. 그들은 이것이 세포들이 빽빽하게 모여 있을 때 서로를 벽으로부터 "차폐(shielding)"하는 방식을 모델이 아직 충분히 반영하지 못했기 때문이라고 추측합니다. 그들은 향후 연구를 통해 이를 수정할 계획입니다.

요약하자면, 그들은 적혈구가 단순히 수동적인 승객이 아니라, 고속도로를 깨끗하게 유지하기 위해 스스로 벽으로부터 밀어내는 능동적인 수영 선수라는 점을 이해하는, 더 나은 혈류 디지털 트윈을 구축했습니다.

기술 요약

기술 요약: 변형 가능한 입자 현탁액을 위한 수정된 현탁액 균형 모델 (Modified Suspension Balance Model)

문제 정의
제한된 기하학적 구조 내에서 적혈구(RBC) 및 기타 변형 가능한 입자 현탁액의 흐름을 모델링하는 것은 본질적인 비선형성, 다체 상호작용, 그리고 비뉴턴 유체 역학적 거동으로 인해 상당한 어려움을 수반한다. 미세 순환에서 나타나는 핵심적인 현상은 세포 결핍층(Cell-Free Layer, CFL)의 형성인데, 이는 벽 근처에서 적혈구가 결핍된 영역을 의미한다. 이 층은 낮은 겉보기 점도를 유지하고 혈소판의 가장자리 배치(margination)를 촉진하는 데 있어 생리학적으로 필수적이다. 세포 규모의 직접적인 계산 모델은 CFL 형성을 포착할 수 있지만, 혈관 네트워크나 장기 규모의 시뮬레이션에는 계산 비용이 너무 많이 든다. 반대로, 기존의 연속체 모델인 전통적인 현탁액 균형 모델(Suspension Balance Model, SBM)은 입자를 벽으로부터 밀어내는 변형 유도 유체역학적 양력(lift force)을 충분히 고려하지 못하기 때문에 CFL을 포착하는 데 실패하는 경우가 많다.

방법론
저자들은 벽-적혈구 유체역학적 상호작용을 설명하기 위해 특정 양력 항을 추가하여 기존 SBM 프레임워크를 보완한 수정된 현탁액 균형 모델(MSBM)을 제안한다.

• 이론적 프레임워크: 이 모델은 현탁액 벌크 및 입자 상에 대한 표준 질량 및 운동량 보존 방정식을 유지한다. 핵심적인 수정 사항은 입자 상의 운동량 균형에 있다. 전통적인 SBM은 입자 상 응력의 발산($\nabla \cdot \langle \Sigma \rangle_p$)과 저해된 침강 항력(hindered settling drag)에 의해 구동되는 입자 이동을 설명한다. MSBM은 변형 가능한 소포체 및 적혈구의 실험적 관찰로부터 유도된 추가적인 체적 평균 양력 항 $\phi L_\perp$를 도입한다.
• 양력 공식화: 양력은 국소 전단율($\dot{\gamma}$), 벽으로부터의 거리($h$), 그리고 환원 부피 매개변수($\nu$)의 함수로 모델링된다. 이 항은 $L_\perp \propto \frac{\dot{\gamma}}{(h+h_0)^\beta}$ 형태를 취하며, 여기서 $\beta$는 입자 의존성(특히 적혈구)을 반영하도록 조정된 멱법칙 계수이고, $h_0$는 벽에서의 특이점 발생을 방지하기 위한 기하학적 인자이다.
• 수치적 구현: MSBM은 OpenFOAM에 구현되었으며, 기존의 SbmFoam 솔버를 확장하여 SbmBloodFoam을 생성하였다. 이 솔버는 유한체적법(Finite-Volume-Method, FVM) 접근 방식을 사용한다. 연속체 모델을 검증하기 위해, 저자들은 격자 볼츠만법(Lattice Boltzmann method, LB)과 조립된 형태의 변형 가능한 캡슐 모델 기반의 실험적으로 검증된 다중 규모 혈류 솔버를 활용하였다.

주요 결과
본 연구는 미세 혈관 채널 및 관 내의 혈류를 시뮬레이션하여 과도 상태 및 정상 상태 거동을 포착하는 MSBM의 능력을 입증한다:

1. CFL 형성 및 속도 프로파일: de conventional SBM이 포물선형 속도 프로파일과 벽에서의 비제로(non-zero) 농도를 예측하는 것과 달리, MSBM은 CFL 형성과 그에 따른 혈액 특유의 "플러그 흐름(plug-flow)" 유사 속도 프로파일을 성공적으로 포착한다. 모델은 양력 매개변수($\beta$)가 증가할수록 CFL 두께가 증가함을 보여준다.
2. 매개변수 민감도:
   • 벌크 헤마토크릿 ($\phi_b$): $\phi_b$가 증가함에 따라 CFL 두께가 감소하며, 이는 문헌 내용과 일치한다.
   • 최대 충전 농도 ($\phi_m$): 적혈구 응집 또는 충전 경향을 나타내는 더 높은 $\phi_m$ 값은 더 두꺼운 CFL과 더 높은 중앙 헤마토크릿 피크를 유발한다.
3. 시간적 진화: 모델은 CFL의 과도적 발달 과정을 포착한다. 초기에는 헤마토크릿이 균일하지만, 시간이 지남에 따라 적혈구는 중심부로 이동하여 벽 근처에 세포 결핍층을 형성하고 농축된 핵(core)을 형성하며, 결국 정상 상태에 도달한다.
4. 이론 및 실험 검증:
   • 단순 전단 흐름: 모델은 단일 입자 이동 실험에서 관찰된 스케일링 법칙과 일치하게, 짧은 시간 척도에서 CFL 두께의 세제곱($\delta_{CFL}^3$)이 무차원 시간($t\dot{\gamma}$)에 대해 선형적으로 스케일링됨을 재현한다.
   • 채널 흐름: MSBM의 속도 및 CFL 두께 예측은 다양한 채널 높이와 벌크 헤마토크릿에 대해 Losserand 및 Salame의 실험 데이터와 잘 일치한다.
   • 관형 흐름: 모델은 파러 효과(Fåhræus Effect) (배출 헤마토크릿이 관 내부 헤마토크릿보다 높은 현상)와 파러-린드퀴스트 효과(Fåhræus-Lindqvist Effect) (혈관 직경이 감소함에 따라 겉보기 점도가 감소하는 현 현상)를 성공적으로 포착한다. 예측된 겉보기 점도와 헤마토크릿 비율은 특히 직경 40 µm 미만의 혈관에 대해 경험적 상관관계(예: Pries 등) 및 실험 데이터와 좋은 일치도를 보인다.

의의 및 주장
본 논문은 제한된 환경 하에서 농축된 변형 가능한 입자 현탁액의 미세구조적 불균질성과 비뉴턴 유체 거동을 효율적으로 포착하는 새로운 연속체 계산 프레임워크를 확립한다. 저자들은 변형 유도 양력을 통합함으로써 MSBM이 계산 비용이 많이 드는 이산 입자 시뮬레이션과 정확도가 낮은 연속체 모델 사이의 간극을 메운다고 주장한다.

본 연구는 대규모 생의학 및 산업적 응용을 위한 효율적인 단계로서 제시된다. 저자들은 모델이 작은 혈관에서는 잘 작동하지만, 큰 관(>40 µm)에서는 CFL과 파러-린드퀴스트 효과를 약간 과대 예측한다는 점을 겸허히 언급하였다. 저자들은 이를 현재의 양력 공식에 입자 밀집 지역에서의 국소적 입자-벽 윤활 상호작용에 의한 양력 감소를 반영하는 "입자 차폐 효과(particle-screening effect)"가 결여되어 있기 때문이라고 설명하였다. 이러한 한계는 현재 모델의 실패가 아닌 향후 연구 과제로 식별되었다.