[논문 리뷰] Stability of the kinematically coupled $\beta$- scheme for fluid-structure interaction problems in hemodynamics

Sunčica Čanić, Boris Muha외 1명|arXiv (Cornell University)|2012. 05. 31.

Advanced Numerical Methods in Computational Mathematics참고 문헌 22인용 수 11

한 줄 요약

이 논문은 혈액역학에서 유체-구조 상호작용에 대해 고정된 β-스키마의 조건부 안정성을 분석적으로 증명한다. 이는 고전적인 딜리클레-뉴턴 스킴과 달리, 로빈 유형의 경계 조건을 통해 구조 관성의 암묵적 통합을 통해 부가 질량 효과를 제거함으로써 이루어진다. 스킴은 모든 β ∈ [0, 1] 및 모든 물리적으로 타당한 매개변수에서 안정하며, 기준 비선형 FSI 문제에 대한 수치 시뮬레이션을 통해 검증되었다.

ABSTRACT

It is well-known that classical Dirichlet-Neumann loosely coupled partitioned schemes for fluid-structure interaction (FSI) problems are unconditionally unstable for certain combinations of physical and geometric parameters that are relevant in hemodynamics. It was shown in \cite{causin2005added} on a simple test problem, that these instabilities are associated with the so called ``added-mass effect''. By considering the same test problem as in \cite{causin2005added}, the present work shows that a novel, partitioned, loosely coupled scheme, recently introduced in \cite{MarSun}, called the kinematically coupled $\beta$-scheme, does not suffer from the added mass effect for any $\beta \in [0,1]$, and is unconditionally stable for all the parameters in the problem. Numerical results showing unconditional stability are presented for a full, nonlinearly coupled benchmark FSI problem, first considered in \cite{formaggia2001coupling}.

연구 동기 및 목표

• 부가 질량 효과로 인해 고전적인 딜리클레-뉴턴 느슨하게 결합된 스킴에서 발생하는 불안정성 문제를 해결한다.
• 문헌 [11]에서 제안한 운동학적으로 결합된 β-스키마가 이러한 불안정성을 극복하는지 조사한다.
• 모든 β ∈ [0, 1]에 대해 스킴의 조건부 안정성에 대한 이론적 및 수치적 증거를 확립한다.
• 고전적인 스킴이 실패하는 실질적인 혈액역학 매개변수에서 스킴의 강건성을 입증한다.
• 혈액 흐름의 비선형이고 완전히 결합된 FSI 문제에 적용 가능한 안정성 분석 프레임워크를 제공한다.

제안 방법

• 부가 질량 효과를 분리하여 분석하기 위해 [18]에서 제시한 단순화된 시험 문제를 채택한다.
• 매개변수 β를 통해 유체와 구조 하위문제 간에 유체 압력을 분배함으로써 운동학적으로 결합된 β-스키마를 수립한다.
• 유체 하위문제에 구조 관성을 암묵적으로 포함하는 로빈 유형의 경계 조건을 도입하여 명시적 딜리클레 결합을 피한다.
• 단순화된 문제에 대해 본 네만 안정성 분석을 수행하여 부가 질량 연산자가 스킴을 불안정하게 만들지 않음을 보여준다.
• 점성 유체와 점탄성 구조를 갖는 비선형 FSI 기준 문제에 스킴을 구현한다.
• 시간에 암묵적이고 분할된 반복 구조를 사용하여 매 시간 단계에서 유체 및 구조 하위문제를 순차적으로 해결한다.

실험 결과

연구 질문

• RQ1 부가 질량 효과가 존재하는 상황에서 운동학적으로 결합된 β-스키마는 모든 β ∈ [0, 1]에 대해 조건부 안정성을 유지하는가?
• RQ2 로빈 경계 조건을 통한 암묵적 구조 관성 통합이 명시적 딜리클레-뉴턴 결합보다 불안정성을 방지하는 방식은 무엇인가?
• RQ3 저밀도 구조 및 높은 유체 관성가치를 갖는 혈액역학에서 흔히 나타나는 극한 매개변수 영역에서도 스킴이 안정성을 유지하는가?
• RQ4 비선형 기준 문제에서 스킴의 수치 성능을 고전적인 딜리클레-뉴턴 스킴 및 단일 스킴과 비교하면 어떻게 되는가?
• RQ5 기준 문제에 대해 스킴의 시간 수렴 특성은 어떠한가? 최적 순서에 도달하는가?

주요 결과

• 운동학적으로 결합된 β-스키마는 모든 β ∈ [0, 1]에 대해 조건부 안정하며, 이는 유체 및 구조 매개변수에 관계없이 성립한다. 이는 고전적인 딜리클레-뉴턴 스킴이 불안정하게 만드는 매개변수 조건에서도 마찬가지다.
• 스킴은 로빈 경계 조건을 통해 암묵적으로 구조 관성을 통합함으로써 부가 질량 효과를 피한다. 이는 명시적 딜리클레 데이터를 통한 직접적 통합을 방지하기 때문이다.
• [31]에서 제시한 기준 FSI 문제에 대한 수치 시뮬레이션은 고전적인 스킴이 실패하는 임계 영역의 매개변수에서도 조건부 안정성을 확인한다.
• 스킴은 시간에 대해 2차 수렴을 달성하며, Δt가 감소함에 따라 압력, 속도 및 변위의 L2 오차 순서는 약 1.0~1.2 수준을 유지한다.
• 운동학적으로 결합된 β-스키마는 바디아, 콰이나, 쿼아티وني의 단일 스킴 결과와 매우 유사한 결과를 생성하여 정확도를 검증한다.
• 그림 11의 블록도 비교는 운동학적으로 결합된 스킴이 딜리클레-뉴턴 스킴에 존재하는 불안정성 유발 피드백 루프를 피하고 있음을 시각적으로 나타낸다.