동맥경화 위험을 가늠하는 카루티 분기점 정적 유동 CFD 분석

초록

본 연구는 실제 환자의 카루티 동맥 분기 기하학을 기반으로 정적 흐름 조건에서 전산유체역학(CFD) 시뮬레이션을 수행하였다. 저벽전단응력(WSS), 흐름 정체 및 분리 영역을 확인했으며, 특히 분기 하류에서 저 WSS가 관찰되어 동맥경화 발생 위험 부위와의 연관성을 제시한다.

상세 분석

이 논문은 최근 25년간 카루티 동맥이 왜 동맥경화에 취약한지를 규명하려는 흐름역학적 접근을 시도한다. 우선 ‘in vivo’ 영상으로부터 양측 경동맥 협착 환자의 실제 기하학을 추출하고, 이를 기반으로 3차원 메쉬를 생성하였다. 메쉬 품질에 대한 언급이 부족하지만, 경계층 해상도를 확보하기 위해 최소 0.1 mm 이하의 셀 크기를 적용했을 것으로 추정된다. 혈액은 뉴턴 유체로 가정하고, 점도 0.0035 Pa·s, 밀도 1060 kg/m³를 사용하였다. 유입 경계조건은 평균 평균 유속(≈0.3 m/s)으로 설정한 정적 흐름이며, 출구는 정압 조건을 부여하였다. 이러한 정적 가정은 실제 맥박 흐름의 시간적 변동성을 무시함으로써 저 WSS 영역을 과소 혹은 과대 평가할 위험이 있다.

수치해석은 Navier‑Stokes 방정식의 비압축성 형태를 사용했으며, 수렴 기준은 연속 방정식 잔차 10⁻⁶ 이하로 설정하였다. 결과는 저 WSS가 분기 하부와 내측벽에서 집중되는 것을 보여준다. 저 WSS는 혈관 내피세포의 기능 저하와 염증 반응을 촉진해 플라크 형성을 유도한다는 기존 연구와 일치한다. 또한, 흐름 정체와 분리 영역이 동일 부위에 겹쳐 나타나는 점은 혈액 내 입자(LDL 등)의 장기 체류를 의미하며, 이는 동맥경화 진행 메커니즘을 뒷받침한다.

하지만 논문은 몇 가지 한계점을 명확히 제시하지 않는다. 첫째, 혈관벽을 강체로 가정했기 때문에 혈관 순응성에 따른 WSS 변동을 반영하지 못한다. 둘째, 단일 환자 사례에 국한되어 통계적 일반화가 어렵다. 셋째, 경계조건이 실제 심박동에 비해 단순화되어 있어, 급격한 흐름 가속·감속 단계에서 발생할 수 있는 고 WSS 피크를 놓쳤을 가능성이 있다. 향후 연구에서는 펄스 흐름, 유동‑구조 상호작용(FSI), 비뉴턴 혈액 모델, 그리고 다수 환자 코호트를 활용한 통계적 분석이 필요하다.

전반적으로, 이 연구는 정적 CFD를 통해 카루티 분기점에서 저 WSS와 흐름 정체가 동맥경화 위험 부위와 일치함을 시각적으로 확인시켜 주었으며, 임상적 위험 평가와 맞춤형 스텐트 설계에 활용될 수 있는 기초 데이터를 제공한다.