관상동맥 무밸브 심장관의 동적 흡입 펌핑과 혈구 흐름

초록

본 연구는 배아 단계에서 밸브가 없는 원통형 심장을 모델링하여, 동적 흡입 펌핑(DSP)과 연동 수축(peristalsis)이 혈류에 미치는 영향을 수치적으로 비교한다. IBM 기반 적응 격자(IBAMR)를 이용해 혈액 세포(혈구)의 부피분율(헤마토크리트)과 와머스리 수(Wo) 변화를 고려했으며, Wo ≤ 10에서는 혈구가 흐름에 큰 영향을 주지 않지만, 전 범위의 Wo에서 연동 수축은 일관된 순환을 유도한다는 결과를 얻었다.

상세 분석

이 논문은 배아 심장의 무밸브 구조가 어떻게 혈액을 순환시키는지를 물리·수학적으로 규명하려는 시도이다. 저자들은 먼저 고전적인 Liebau 펌프(동적 흡입 펌핑, DSP)의 기본 메커니즘을 재조명한다. 비대칭적인 수축 부위가 유연한 관에 가해질 때, 파동이 강체 경계에서 반사되어 비대칭적인 압력 구배를 만든다. 이 압력 구배가 순환 흐름을 유도한다는 것이 DSP의 핵심이다. 그러나 기존 모델들은 대부분 점성 무시, 1차원 흐름, 작은 변형 등 과도한 단순화를 전제로 했으며, 혈구와 같은 큰 입자의 존재를 배제했다.

저자들은 이러한 한계를 극복하기 위해 Immersed Boundary Method(IBM)를 기반으로 한 IBAMR 프레임워크를 채택했다. IBM은 유연한 구조와 점성 유체를 동시에 시뮬레이션할 수 있는 강력한 도구이며, 적응 격자 기법을 통해 혈구와 관벽 근처에서 높은 해상도를 유지하면서도 전체 계산 비용을 절감한다. 특히, 혈구를 탄성 스프링 네트워크로 모델링함으로써 실제 혈구의 변형과 상호작용을 정량화했다.

시뮬레이션 설정은 두 가지 주요 시나리오로 나뉜다. 첫 번째는 DSP 모델로, 관의 하부 평탄 부분만을 유연하게 하고, 그 안에 10% 길이의 비대칭적인 근육 스프링을 부착해 수축을 구현한다. 스프링의 휴식 길이는 시간에 따라 사인파 형태로 변조되어, 실제 근육 수축 주기를 모사한다. 두 번째는 연동 수축 모델로, 관 전체에 파동을 전파시키는 방식으로 압축-이완을 순차적으로 진행한다. 두 모델 모두 동일한 폐쇄형 ‘레이스트랙’ 형태의 관을 사용했으며, 혈구 부피분율을 0%부터 80%까지 10% 간격으로 변화시켰다.

핵심 변수인 와머스리 수(Wo)는 관 직경, 유체 점성, 진동 주파수의 조합으로 정의되며, 저자들은 Wo ≤ 10(저주파, 점성 지배)과 Wo > 10(관성 지배) 두 영역을 집중적으로 탐색한다. 결과는 흥미롭다. 낮은 Wo 영역에서는 혈구가 존재해도 평균 유량에 유의미한 차이가 없었다. 이는 점성력이 지배적인 상황에서 혈구가 흐름을 방해하거나 촉진하는 효과가 상쇄된 것으로 해석될 수 있다. 반면, 연동 수축은 모든 Wo와 모든 헤마토크리트에서 지속적인 순환을 생성했으며, 특히 높은 Wo에서 유량이 크게 증가했다. 이는 연동 수축이 파동 전파와 관강성의 조합을 통해 보다 효율적인 압력 구배를 만들기 때문이다.

또한, 저자들은 파동 반사와 비대칭 수축 위치가 DSP 효율에 미치는 민감도를 분석했다. 강체 경계와 유연 부위의 비율, 수축 부위의 길이(Ls)와 위치(LA) 등을 변형했을 때, 유량이 급격히 변하거나 흐름이 역전되는 현상이 관찰되었다. 이는 DSP가 설계 파라미터에 매우 민감함을 시사한다.

전반적으로 이 연구는 두 가지 주요 기여를 제공한다. 첫째, 혈구와 같은 큰 입자를 포함한 2D IBM 시뮬레이션을 최초로 구현함으로써 배아 심장 내 혈류 역학을 보다 현실적으로 재현했다. 둘째, 다양한 Wo와 헤마토크리트 조건에서 DSP와 연동 수축의 효율을 정량적으로 비교함으로써, 실제 배아 심장에서 어느 메커니즘이 우세할지에 대한 근거를 제시했다. 결과는 연동 수축이 더 안정적이고 효율적인 펌핑 메커니즘임을 뒷받침하지만, 특정 저주파 환경에서는 DSP도 충분히 순환을 유지할 수 있음을 보여준다. 이러한 통찰은 향후 미세유체 펌프 설계, 조직공학 심장 모형 개발, 그리고 배아 심장 발달 연구에 중요한 지침이 될 것이다.