[논문 리뷰] A note on robust preconditioners for monolithic fluid-structure interaction systems of finite element equations

Ulrich Langer, Huidong Yang|arXiv (Cornell University)|2014. 12. 21.

Advanced Numerical Methods in Computational Mathematics참고 문헌 55인용 수 12

한 줄 요약

이 논문은 결합된 유한요소 행렬의 근사 $LDU$ 분해를 사용하여 단일체 유체-구조 상호작용(FSI) 시스템에 대한 강건한 조건부 행렬을 제안한다. $L$, $D$, $U$를 근사된 요소로 구성하고, 슈어 여집합을 근사하기 위해 특수화된 대수적 다중격자(AMG) 순환을 사용함으로써, 메esh, 유체 밀도, 시간 간격 크기 모두에 대해 강건성을 확보한다. 이로 인해 이전의 조건부 행렬 대비 GMRES 반복 횟수를 약 10배 감소시킨다.

ABSTRACT

In this note, we consider preconditioned Krylov subspace methods for discrete fluid-structure interaction problems with a nonlinear hyperelastic material model and covering a large range of flows, e.g, water, blood, and air with highly varying density. Based on the complete $LDU$ factorization of the coupled system matrix, the preconditioner is constructed in form of $\hat{L}\hat{D}\hat{U}$, where $\hat{L}$, $\hat{D}$ and $\hat{U}$ are proper approximations to $L$, $D$ and $U$, respectively. The inverse of the corresponding Schur complement is approximated by applying one cycle of a special class of algebraic multigrid methods to the perturbed fluid sub-problem, that is obtained by modifying corresponding entries in the original fluid matrix with an explicitly constructed approximation of the exact perturbation coming from the sparse matrix-matrix multiplications.

연구 동기 및 목표

• 기존의 단일체 FSI 시스템을 위한 조건부 행렬이 유체 밀도 변화에 대해 강건하지 못한 문제를 해결하기 위해.
• 블록 가우스-세이델 및 ILU 유형의 조건부 행렬에서 관찰되는 메쉬 의존성과 시간 간격 크기 민감성 문제를 극복하기 위해.
• 물과 공기와 같은 다양한 유체 성질에서도 효율성과 수렴 강건성을 유지하는 조건부 행렬을 개발하기 위해.
• 완전히 결합된 단일체 다중격자 솔버에 의존하지 않고도 강건성을 확보함으로써 구현 및 계산 오버헤드를 줄이기 위해.

제안 방법

• 조건부 행렬은 결합된 FSI 시스템 행렬에 대한 근사 $LDU$ 분해를 통해 구성되며, $L$, $D$, $U$를 각각 정확한 $L$, $D$, $U$ 요소의 근사로 사용한다.
• 슈어 여집합의 역행렬은 유체 하위행렬에 흐름을 가진 변형을 가한 후, 특수화된 대수적 다중격자(AMG) 방법을 한 번의 순환으로 적용하여 근사한다.
• 변형은 희소 행렬-행렬 곱셈을 통해 명시적으로 유도되며, 구조가 유체 시스템에 미치는 결합 효과를 반영한다.
• 완전한 $LDU$ 분해를 피하기 위해 저랭크 근사와 구조적 스무딩 기법을 사용함으로써 계산 효율성을 유지한다.
• 조건부 행렬은 GMRES 및 유연한 GMRES 솔버에 통합되었으며, 다양한 유체 밀도와 메쉬 정밀도에서 성능이 평가되었다.

실험 결과

연구 질문

• RQ1 근사 $LDU$ 분해 기반의 조건부 행렬이 다양한 유체 밀도에서 단일체 FSI 시스템에 대해 강건한 수렴 성능을 달성할 수 있는가?
• RQ2 제안된 조건부 행렬이 메쉬 정밀도와 시간 간격 크기 변화에 대해 강건성을 유지하는가?
• RQ3 반복 횟수와 확장성 측면에서 블록 가우스-세이델, SSOR, ILU(0) 유형과 비교해 새로운 조건부 행렬의 성능은 어떠한가?
• RQ4 흐름이 가미된 유체 행렬에 대해 단일 AMG 순환을 사용해 슈어 여집합의 역행렬을 효율적으로 근사할 수 있는가?

주요 결과

• 새로운 조건부 행렬은 $ {P}_{ILU}$ 및 $ {P}_{SSOR}$와 같은 이전의 조건부 행렬 대비 GMRES 반복 횟수를 약 10배 감소시키며, 특히 공기($ ho_f = 0.0011$ g/cm³)와 같은 저밀도 유동에서 뚜렷한 효과를 보인다.
• 유체 밀도 $ ho_f = 1.1$ g/cm³(물에 가까운 밀도)일 경우, 추가 질량 효과에 대해 강건성을 유지하며 낮은 반복 횟수를 유지한다.
• 시간 간격 크기 $ {Δ}t = 0.00125$ ms일 때, 세 가지 메쉬 수준에서 반복 횟수가 거의 일정하게 유지되어 메쉬 정밀도 변화에 대해 강건함을 보였다.
• 유체 밀도 $ { ho_f} \in \{1.1, 0.11, 0.011, 0.0011\}$ g/cm³의 광범위한 범위에서 안정적인 성능을 유지하였으며, 밀도가 감소함에 따라 반복 횟수는 안정되거나 감소하는 경향을 보였다.
• 유연한 GMRES는 표준 GMRES보다 略적으로 열악한 성능를 보였지만, 양쪽 모두 이전의 접근 방식보다 훨씬 효율적인 성능를 확보하였다.