[논문 리뷰] Mathematical modelling of the cardiovascular system

Alfio Quarteroni|ArXiv.org|2003. 05. 01.

Elasticity and Material ModelingEngineering참고 문헌 15인용 수 18

한 줄 요약

이 논문은 유동-구조 상호작용 모델을 기반으로 한 임의 라그랑주형 유동(ARLE) 공식화된 라그랑주-뉴턴 방정식을 사용하여 대동맥 및 중대동맥에서의 혈류를 시뮬레이션하기 위한 수학적 프레임워크를 제시한다. 유체역학과 혈관벽 기계적 거동을 결합하기 위한 단일화된, 분할 반복 알고리즘을 도입하였으며, 3D, 1D, 그리고 집합 모델을 통합하는 다스케일 모델링 접근법을 제안하여 전체 순환계를 더 높은 정확도와 계산 효율성으로 시뮬레이션할 수 있도록 한다.

ABSTRACT

In this paper we will address the problem of developing mathematical models for the numerical simulation of the human circulatory system. In particular, we will focus our attention on the problem of haemodynamics in large human arteries.

연구 동기 및 목표

• 탄성 대동맥에서의 혈류를 정확하게 시뮬레이션하기 위한 강력한 수학적 및 수치적 프레임워크 개발.
• 동적 혈관벽 운동을 고려하여 이동하는 영역에서의 유동-구조 상호작용 문제를 정확한 영역 매핑을 위해 임의 라그랑주형(ALE) 공식화를 사용하여 해결.
• 세부적인 3D 모델과 단순화된 순환계 모델을 결합하여 유동 분리, 벽면 전단력, 수술 후 결과와 같은 임상적으로 관련 있는 혈역학 현상을 시뮬레이션할 수 있도록 하기.
• 3D 유동-구조 상호작용, 1D 초구형 시스템, 그리고 집합 모델을 통합하는 계층적 다스케일 모델링 전략을 수립하여 전체 순환계의 시뮬레이션을 가능하게 하기.
• 하향적 암시적 반복 분할 알고리즘을 통해 수치적 안정성과 정확성을 확보하고, 하향 반복 및 조화 확장에 의한 메쉬 변형을 통해 결합 시뮬레이션의 안정성을 높이기.

제안 방법

• 시간에 따라 변화하는 유체 영역과 이동 경계(특히 탄성 혈관벽)를 다룰 수 있도록 라그랑주-뉴턴 방정식을 ALE 프레임워크에 수식화.
• 유체의 속도와 압력, 혈관벽의 변위와 속도를 동시에 해결하는 단일화된 접근법을 사용하여 비압축성 라그랑주-뉴턴 방정식과 혈관벽의 구조 모델을 결합.
• 유체와 구조 간의 경계에서 운동학적 및 역학적 일致성을 확보하기 위해 하향 반복 기반의 분할 반복 알고리즘을 사용하며, 외삽된 벽의 운동과 속도를 이용해 ALE 매핑을 업데이트.
• 메쉬 변형의 품질을 유지하고 과도한 왜곡을 방지하기 위해 조화 확장을 사용하여 벽의 변위에서 ALE 매핑을 계산.
• 3D 유동-구조 상호작용 모델과 1D 축소 모델(초구형 보존 법칙), 그리고 집합 모델(전기 회로에 유사한 상미분 방정식계)을 통합하여 전체 순환계의 표현을 가능하게 함.
• 반복 간 변위 및 속도 차이의 L2 노름을 기반으로 수렴성 검사를 수행하며, 반복 과정의 안정성을 확보하기 위해 가속 계수를 적용함.

실험 결과

연구 질문

• RQ1 혈류의 박동성에 의해 변형되는 혈관벽을 고려할 때, 대동맥에서의 유동-구조 상호작용 문제를 어떻게 정확하게 모델링할 수 있는가?
• RQ2 이동 영역에서 라그랑주-뉴턴 방정식과 구조역학 모델을 결합하기 위한 가장 안정적이고 정확한 수치 전략은 무엇인가?
• RQ3 다스케일 모델링 기법을 어떻게 활용하여 지역 혈역학 정확도를 유지하면서 전체 순환계를 시뮬레이션할 수 있는가?
• RQ4 결합 알고리즘의 선택(예: 분할 대비 단일화)이 유동-구조 시뮬레이션의 안정성과 수렴성에 미치는 영향은 무엇인가?
• RQ5 경화증과 같은 병적 상태의 영향을 변화된 혈역학을 수치적으로 시뮬레이션하여 정량적으로 평가할 수 있는가?

주요 결과

• ALE 공식화는 특히 혈관벽이 최대 10%의 반경 방향 변형을 겪을 경우에도 정확하고 안정적인 시간 변화 영역에서의 혈류 시뮬레이션을 가능하게 한다.
• 가속 계수를 적용한 분할 반복 알고리즘은 수렴성과 안정성을 보장하며, 수치적 불안정성을 방지하기 위해 이산 수준에서 강한 결합이 필요함을 보여준다.
• 메쉬 변형에 조화 확장을 사용함으로써 메쉬 품질을 유지하고 장기간의 시뮬레이션 기간 동안 과도한 왜곡을 방지할 수 있다.
• 제안된 다스케일 프레임워크는 3D 유동-구조 상호작용 모델과 1D, 집합 모델을 성공적으로 결합하여 계산 비용을 줄인 전체 순환계 시뮬레이션을 가능하게 한다.
• 해부학적 기하학의 영향을 반영하여 혈류 패턴과 벽면 전단력에 영향을 미치는 관계를 포착함으로써 관상동맥 우회 수술의 혈역학적 특성을 임상적으로 관련 있는 방식으로 시뮬레이션할 수 있다.
• 하향 반복마다 유체, 구조, 메쉬 변형 문제를 반복적으로 해결해야 하므로 알고리즘이 계산적으로 비용이 많이 들지만, 복잡한 혈역학 시나리오에서 안정성을 확보하기 위해 필수적이다.