[논문 리뷰] High order finite element simulations for fluid dynamics validated by experimental data from the fda benchmark nozzle model

Vincent Chabannes, Christophe Prud'Homme외 2명|arXiv (Cornell University)|2017. 01. 09.

Cardiovascular Function and Risk Factors인용 수 3

한 줄 요약

이 논문은 개방 소스 Feel++ FEM 라이브러리를 사용하여 FDA 기준 노즐에서 혈류의 고차수 유한요소 시뮬레이션을 제시하며, 레이놀즈 수 Re = 500, 2000, 3500인 레이놀즈 수가 낮은, 전이 및 난류 영역에서 실험 데이터와의 검증을 수행한다. 축방향 속도와 압력 차이에 대해 실험 데이터와 뛰어난 일치를 보이며, 검증 지표 $E_z$와 $E_Q$는 특히 고차수 기하 및 속도-압력 이산화에서 높은 정확도와 강건성을 보여준다.

ABSTRACT

The objective of the present work is to construct a sound mathematical, numerical and computational framework relevant to blood flow simulations and to assess it through a careful validation against experimental data. We perform simulations of a benchmark proposed by the FDA for fluid flow in an idealized medical device, under different flow regimes. The results are evaluated using metrics proposed in the literature and the findings are in very good agreement with the validation experiment.

연구 동기 및 목표

• 환자 맞춤 혈류 시뮬레이션을 위한 신뢰성 있고 재현 가능하며 고정밀도 수학적 및 계산 프레임워크를 구축하기 위해.
• FDA의 기준 노즐 모델에서의 실험 데이터와 비교하여 개방 소스 Feel++ 라이브러리 내 고차수 유한요소법(FEM) 해석기의 검증을 위해.
• 메쉬 정밀도, 다항식 차수(P2, P3), 기하 구조 근사(G1, G2)가 복잡한 유동 영역에서 시뮬레이션 정확도에 미치는 영향을 평가하기 위해.
• 다양한 레이놀즈 수에서 핵심 검증 지표 $E_z$와 $E_Q$를 평가하여 실험 데이터와의 일치 정도를 정량화하기 위해.
• 의료 기기 개발에서 CFD의 신뢰성과 재현 가능성을 향상시키기 위해 기여하기 위해.

제안 방법

• 안정화된 갈루아킨 유한요소 방법을 사용하여 비압축성 비정상 나비에-스토크스 방정식을 보존 형태로 풀이한다.
• 속도, 압력 및 기하 구조 근사를 위해 일반화된 테일러-후드 유한요소 쌍 $P_{N+1}P_NG_{k_{\text{geo}}}$를 사용한다.
• PETSc 기반의 조절자(예: PCD)와 영역 분할을 사용하여 효율적인 병렬 해법을 위한 개방 소스 Feel++ 라이브러리를 활용한다.
• Gmsh를 사용하여 구조적이고 이방성인 테트라헤드론 메쉬를 생성하며, 목구장 및 벽 경계와 같은 중요한 영역에서 메쉬를 정밀화한다.
• 초기 조건으로 $\mathbf{u} = \mathbf{0}$을 설정하고, $\Gamma_{\text{inlet}}$에서 포아젤 유입 프로파일, $\Gamma_{\text{wall}}$에서 슬립 없는 경계 조건, $\Gamma_{\text{outlet}}$에서 자유 유출 조건을 적용한다.
• 순간적 행동을 해석하기 위해 시간 이산화에 $\Delta t = 10^{-3}$ s를 사용하여 Re = 500의 경우를, Re = 2000 및 3500의 경우 $\Delta t = 10^{-4}$ s를 사용한다.

실험 결과

연구 질문

• RQ1 고차수 유한요소 방법은 FDA 노즐 기준 시험에서 실험 데이터의 축방향 속도 및 압력 프로파일을 얼마나 정확하게 재현할 수 있는가?
• RQ2 메쉬 정밀도 및 다항식 차수(P2, P3)가 시뮬레이션 지표 $E_z$와 $E_Q$의 수렴성과 정확도에 미치는 영향은 무엇인가?
• RQ3 고차수 기하 구조 근사(G1 대비 G2)는 전이 및 난류 영역에서 실험 데이터와의 일치도에 어떤 영향을 미치는가?
• RQ4 왜 Re = 2000 시뮬레이션에서 수치적 제트 붕괴 지점이 실험 결과보다 더 하류에 위치하는가?
• RQ5 개방 소스 Feel++ FEM 프레임워크는 혈역학 시뮬레이션에서 실험 기준과 비교해 강건하고 정확하며 재현 가능한 결과를 달성할 수 있는가?

주요 결과

• Re = 500일 때, 모든 메쉬 정밀도 및 다항식 차수에서 정규화된 축방향 속도와 압력 차이가 실험 데이터와 뛰어난 일치를 보였다.
• 모든 시뮬레이션에서 검증 지표 $E_z$는 항상 낮게 유지되어(< 0.05) 계산된 축방향 속도와 실험 데이터 간 강력한 일치를 나타내었다.
• 질량 보존 오차 $E_Q$는 가장 굵은 메쉬(M0)를 제외한 모든 경우에서 2% 이하였으며, 이 경우 두 곳에서 약 10%에 이르렀다. 이는 메쉬 해상도에 민감함을 시사한다.
• $P_2P_1G_2$ 근사는 가장 굵은 메쉬에서 $E_Q < 2\%$를 달성했고, $P_3P_2G_1$는 M0에서 정확도 향상에 기여하지 않았다. 이는 굵은 메쉬에서 기하 구조 근사가 속도 차수보다 더 중요하다는 것을 시사한다.
• Re = 2000 시뮬레이션에서 수치적 제트 붕괴 지점은 실험 결과보다 더 하류에 위치해 있었으며, 이는 수치적 통합 정확도 또는 경계 조건 조정 문제일 수 있다.
• 본 연구는 고차수 FEM이 적절한 메쉬 및 기하 구조 근사를 통해 복잡한 혈역학 유동에서 매우 정확하고 재현 가능한 CFD 시뮬레이션을 가능하게 함을 확인하였다.