A Lattice Boltzmann Method for Non-Newtonian Blood Flow in Coiled Intracranial Aneurysms

본 논문은 코일 덩어리를 불균질 다공성 매체로 취급하여 코일된 뇌동맥류 내의 비뉴턴성 혈류를 모델링하는 환자 맞춤형 격자 볼츠만 방법을 제시함으로써, 치료 후 혈역학을 평가하는 워크플로우를 검증합니다.

핵심

이 논문은 쉬운 언어와 일상적인 비유를 사용하여 설명합니다.

큰 그림: 누수하는 풍선 수리

머리 속에 있는 물 풍선에 약한 부분이 생겼다고 상상해 보세요. 의학 용어로 이를 뇌동맥류라고 합니다. 이것이 터지면 위험한 뇌졸중을 유발합니다.

이를 수리하기 위해 의사는 종종 코일링이라는 기술을 사용합니다. 그들은 풍선 안으로 아주 작은 와이어를 삽입하고 금속 코일 그물로 채웁니다. 이는 누수하는 양동이에 스펀지를 밀어 넣는 것과 같습니다. 목표는 풍선 내부의 급류하는 물 (혈액) 을 늦추어 약한 벽을 때리는 것을 멈추게 하고, 풍선이 치유되거나 흉터가 생기도록 기회를 주는 것입니다.

문제: 결과 예측의 어려움

어려운 점은 각 환자의 '누수하는 풍선' 모양이 다르고, 각 의사가 '스펀지'(코일) 를 약간 다른 방식으로 채운다는 것입니다. 수술 전에는 코일이 삽입된 후 물의 흐름이 어떻게 변할지 정확히 예측하기 매우 어렵습니다. 현재 의사는 치료가 성공할지 여부를 추측하기 위해 과거 경험에 크게 의존하고 있습니다.

해결책: 디지털 '비행 시뮬레이터'

이 논문의 저자들은 혈류에 대한 비행 시뮬레이터 역할을 하는 컴퓨터 프로그램을 개발했습니다. 추측 대신 코일이 삽입된 후 혈액이 어떻게 행동할지 정확히 보기 위해 가상 테스트를 실행할 수 있습니다.

그들이 이 시뮬레이터를 구축한 방법은 다음과 같습니다.

1. '다공성 스펀지' 트릭
보통 금속 와이어 (코일) 의 엉킨 더미를 통해 흐르는 물을 시뮬레이션하려면 컴퓨터가 모든 단일 와이어를 그려야 합니다. 이는 양동이에 있는 모든 모래 알갱이를 세어 보려는 것과 같습니다. 시간이 영원히 걸리고 슈퍼컴퓨터가 필요합니다.

저자들은 더 똑똑한 방법을 찾았습니다. 모든 와이어를 그리는 대신 코일 더미를 스펀지로 취급한 것입니다.

• 비유: 빽빽한 숲을 통해 물을 부은다고 상상해 보세요. 모든 나무 줄기를 돌아다니는 물의 흐름을 계산해 볼 수 있습니다 (매우 어려움). 또는 "이 전체 지역은 느리게 움직이는 진흙탕이다"라고 말할 수 있습니다 (훨씬 쉬움).
• 과학: 그들은 **부피 평균 나비에 - 스토크스 (VANSE)**라는 수학적 모델을 사용했습니다. 이는 코일로 채워진 영역을 스펀지가 얼마나 빽빽한지에 따라 물이 느려지는 '다공성 매체 (스펀지)'로 취급합니다.

2. '똑똑한 혈액' 요소
혈액은 물과 다릅니다. 혈액은 '끈적'하며 속도에 따라 흐름이 변합니다 (케첩이나 꿀과 같음). 컴퓨터 모델은 이러한 '비뉴턴' 행동을 고려하여 시뮬레이션이 단순한 물이 아닌 실제 혈액처럼 느껴지도록 합니다.

3. '빠른 엔진' (래티스 볼츠만 방법)
이러한 계산을 유용할 정도로 빠르게 실행하기 위해 그들은 **래티스 볼츠만 방법 (LBM)**이라는 특정 수학적 엔진을 사용했습니다.

• 비유: 이는 고속 비디오 게임 엔진과 같습니다. 다른 방법들이 전체 바다의 물리학을 한 번에 해결하려 한다면, LBM 은 바다를 작고 관리 가능한 타일로 나누고 입자가 그 사이에서 어떻게 튀는지 시뮬레이션합니다. 이를 통해 시뮬레이션은 현대 그래픽 카드 (GPU) 에서 놀라울 정도로 빠르게 실행될 수 있습니다.

연구에서 수행한 작업

팀은 환자의 뇌동맥류에 대한 실제 CT 스캔을 가져와 3D 모델을 구축했습니다. 그런 다음 두 가지 유형의 시뮬레이션을 실행했습니다.

1. '초정밀' 버전: 코일의 모든 단일 와이어를 모델링했습니다 (숲의 모든 나무를 세는 것과 같음).
2. '스펀지' 버전: 새로운 '다공성 매체' 모델을 사용했습니다 (숲을 진흙탕으로 취급하는 것과 같음).

그들은 수술에서 의사가 실제로 사용하는 양인 15%, 20%, 25% 의 세 가지 다른 '꽉 찬 정도'(충전 밀도) 로 이를 테스트했습니다.

결과: 스펀지가 작동합니다!

결과는 고무적이었습니다.

• 정확도: '스펀지' 모델은 '초정밀' 모델과 거의 동일한 결과를 제공했습니다. '스펀지' 방법은 훨씬 더 빠르면서도 큰 그림을 잃지 않았습니다.
• 유량 감소: 코일을 더 많이 추가할수록 (스펀지를 더 밀도 있게 만들수록) 동맥류 내부의 혈류가 현저히 느려졌습니다.
• 안전성: '벽 전단 응력'(약한 벽을 문지르는 혈액의 힘) 이 극적으로 감소했습니다. 치료되지 않은 동맥류에서는 벽이 강하게 타격받고 있었지만, 코일을 사용하면 힘이 약 40% 감소하여 풍선이 터질 위험이 훨씬 낮아졌음을 시사합니다.

결론

이 논문은 뇌동맥류 치료를 시뮬레이션하는 새로운 더 빠른 방법을 제시합니다. 모든 와이어를 세는 대신 금속 코일을 '스펀지'로 취급함으로써 의사는 환자 맞춤형 시뮬레이션을 빠르게 실행할 수 있게 될 것입니다. 이 워크플로우는 특정 치료 계획이 혈류를 성공적으로 늦추고 환자를 보호할 수 있는지 더 잘 평가할 수 있게 하여, 순수한 추측에서 데이터 기반 결정으로 나아가게 합니다.

기술 요약

기술 요약: 뇌동맥류 코일 삽입 시 비뉴턴성 혈류에 대한 격자 볼츠만 방법

문제 제기
뇌동맥류는 출혈성 뇌졸중의 주요 원인이며, 종종 혈류 감소 및 벽 전단 응력 (WSS) 감소를 유도하기 위해 코일 색전술로 치료됩니다. 그러나 이러한 치료의 혈역학적 결과를 예측하는 것은 여전히 해결되지 않은 과제로 남아 있습니다. 정확한 예측을 위해서는 환자별 해부학적 구조와 매개변수, 특히 동맥류 내부의 불균질한 코일 분포를 표현하는 것이 필수적입니다. 전통적인 접근법은 다음과 같은 트레이드오프에 직면해 있습니다: 이산적인 코일 기하구조를 완전히 해석하는 시뮬레이션은 계산 비용이 매우 높고, 단순화된 모델들은 종종 코일 패킹의 이질성을 포착하지 못해 올바른 유동 예측에 실패합니다. 또한, 혈액은 비뉴턴성 유체이므로 시뮬레이션에 유변학적 복잡성이 추가됩니다.

방법론
저자들은 환자별 해부학적 구조와 부피 평균 접근법을 결합하여 코일이 삽입된 동맥류를 불균질한 다공성 매체로 모델링하는 워크플로우를 제안합니다. 핵심 방법론은 다음과 같습니다:

1. 지배 방정식: 이 모델은 비뉴턴성 유동에 대한 부피 평균 나비에 - 스토크스 방정식 (VANSE) 을 사용합니다. 유체 영역은 유체의 부피 분율을 나타내는 공극률 필드 ($\phi$) 로 특징지어집니다. 운동량 방정식에는 다공성 매체 (코일) 를 고려하기 위해 다르시 및 포르하이머 항력을 모델링하는 체적력 항 ($f(u)$) 이 포함됩니다.
2. 유변학: 혈액 점도는 전단 박리 (shear-thinning) 거동을 고려한 카루 - 야수다 (Carreau–Yasuda) 모델을 사용하여 모델링됩니다. 동점도는 전단 변형률 속도에 의존합니다.
3. 수치 기법: 방정식은 D3Q27 격자에서 문제 적응형 격자 볼츠만 방법 (LBM) 을 사용하여 해결됩니다.
   • 다공성 매체 처리: LBM 은 비뉴턴성 점도와 다공성 항력을 처리하기 위해 공간적으로 변화하는 공극률과 가변적인 완화율을 통합합니다.
   • 힘 구현: 속도에 이차적으로 의존하는 항력을 포함시키기 위해 구오 (Guo) 힘 스킴이 사용됩니다.
   • 속도 해: 등방성 투과도의 경우 속도는 명시적으로 해결됩니다. 이방성 투과도의 경우 고정점 반복법이 제안되지만, 저자들은 이것이 계산 시간을 크게 증가시킨다고 지적합니다.
4. 기하구조 및 시뮬레이션 설정:
   • 기하구조: 익명화된 CT 스캔에서 유래된 환자별 동맥류 기하구조가 사용되었습니다.
   • 코일 표현: 코일은 15%, 20%, 25% 의 패킹 밀도를 가진 이산 탄성 막대 모델로 시뮬레이션되었습니다. 결과적으로 생성된 이산 코일 구조는 윈도우 컨볼루션을 통해 연속적인 공극률 필드로 변환되었습니다.
   • 경계 조건: 벽면에는 미끄럼 없는 조건 (half-way bounce-back) 이 적용되었습니다. 유입부는 1D-0D 동맥 네트워크 모델에서 유래된 시간 분해 속도 프로파일에 의해 구동되었으며, 유출부는 외삽법을 사용했습니다.
   • 구현: 시뮬레이션은 병렬 계산을 위해 GPU 가속 (Nvidia RTX A6000) 을 활용하는 XLB 프레임워크를 사용하여 구현되었습니다.

주요 기여

• 하이브리드 모델링 접근법: 이 논문은 코일 삽입된 동맥류 내 혈류를 시뮬레이션하기 위해 VANSE 를 위한 LBM 과 비뉴턴성 유변학을 통합하며, 개별 와이어를 해석하는 대신 코일 덩어리를 다공성 매체로 취급합니다.
• 평균화의 검증: 이 연구는 완전히 해석된 코일 시뮬레이션과 부피 평균 다공성 매체 접근법 간의 직접적인 비교를 제공합니다.
• 환자별 워크플로우: 저자들은 의료 영상 (CT 스캔) 에서 혈역학적 평가에 이르기까지, 시뮬레이션된 코일 배치에서 공극률 필드를 생성하는 완전한 파이프라인을 시연합니다.

결과

• 유동 감소: 시뮬레이션은 코일 삽입이 모혈관 내 유동을 유지하면서 동맥류낭 내 유동 속도를 현저히 감소시킨다는 것을 확인시켜 줍니다. 감소는 더 높은 패킹 밀도 (최대 25%) 에서 더 두드러집니다.
• 모델 유효성: 부피 평균 (다공성 매체) 시뮬레이션은 완전히 해석된 시뮬레이션과 좋은 일치를 보입니다. 평균화 접근법은 국소적인 미세 규모의 효과를 잃지만, 전체적인 유동 거동과 거시적 양은 일관됩니다.
• 벽 전단 응력 (WSS): 치료는 WSS 크기를 극적으로 감소시킵니다. 치료받지 않은 동맥류에서 WSS 는 약 2 Pa 에 달했으나, 코일 삽입된 동맥류에서는 1.5 Pa 미만의 크기를 보여 파열 위험이 감소했음을 나타냅니다. 세 가지 패킹 밀도 (15%, 20%, 25%) 간의 WSS 차이는 미미한 것으로 밝혀졌으며, 이는 패킹 밀도 단독보다는 코일 배치 분포가 더 큰 영향을 미칠 수 있음을 시사합니다.
• 정량적 지표: 동맥류 영역 전체에 걸쳐 부피 평균 모델은 치료받지 않은 경우 대비 유동 크기 약 50% 감소와 WSS 약 40% 감소를 예측했습니다.

의의 및 주장
저자들은 이 워크플로우가 잠재적인 코일 치료로 인한 혈류 변화에 대한 환자별 평가를 가능하게 한다고 주장합니다. 주요 의의는 완전히 해석된 시뮬레이션에 대한 계산적으로 효율적인 대안을 제공하면서도 전역 유동 지표의 정확성을 희생하지 않는 부피 평균 다공성 매체 접근법의 유효성에 있습니다.

이 논문은 이방성 투과도를 고려하는 것이 이론적으로 정확도를 높일 수 있지만, 정확한 코일 배치에 대한 불확실성과 신속한 의사결정의 필요성을 고려할 때 현재 등방성 가정이 임상 실무에 충분할 것이라고 겸손하게 결론 내립니다. 저자들은 향후 작업으로 상세한 매개변수 보정, 불확실성 정량화, 그리고 이방성 투과도 확장에 대한 분석이 필요하다고 지적합니다.