[논문 리뷰] Scalable smoothing strategies for a geometric multigrid method for the immersed boundary equations

Amneet Pal Singh Bhalla, Matthew G. Knepley외 3명|arXiv (Cornell University)|2016. 12. 07.

Lattice Boltzmann Simulation Studies참고 문헌 34인용 수 5

한 줄 요약

이 논문은 스토크스 및 레이놀즈 수가 낮거나 중간 수준인 유동 조건에서 잠긴 경계(IB) 방정식에 적용된 기하학적 멀티그리드 솔버를 위한 확장 가능하고 병렬 처리 가능한 스무딩 전략을 개발한다. 추가적인 반복자(Vanka) 스무딩기와 슈어 여부 기반 스무딩기를 도입하여 저-중간 레이놀즈 수 및 높은 재료 강성 조건에서 효과적인 성능을 입증하였으며, 복잡한 유체-구조 상호작용 문제에서의 확장성도 확보하였다.

ABSTRACT

The immersed boundary (IB) method is a widely used approach to simulating fluid-structure interaction (FSI). Although explicit versions of the IB method can suffer from severe time step size restrictions, these methods remain popular because of their simplicity and generality. In prior work (Guy et al., Adv Comput Math, 2015), some of us developed a geometric multigrid preconditioner for a stable semi-implicit IB method under Stokes flow conditions; however, this solver methodology used a Vanka-type smoother that presented limited opportunities for parallelization. This work extends this Stokes-IB solver methodology by developing smoothing techniques that are suitable for parallel implementation. Specifically, we demonstrate that an additive version of the Vanka smoother can yield an effective multigrid preconditioner for the Stokes-IB equations, and we introduce an efficient Schur complement-based smoother that is also shown to be effective for the Stokes-IB equations. We investigate the performance of these solvers for a broad range of material stiffnesses, both for Stokes flows and flows at nonzero Reynolds numbers, and for thick and thin structural models. We show here that linear solver performance degrades with increasing Reynolds number and material stiffness, especially for thin interface cases. Nonetheless, the proposed approaches promise to yield effective solution algorithms, especially at lower Reynolds numbers and at modest-to-high elastic stiffnesses.

연구 동기 및 목표

• 잠긴 경계(IB) 방법을 위한 기하학적 멀티그리드 솔버에서 기존의 Vanka 유형 스무딩기의 병렬 확장성에 한계가 있음을 해결한다.
• 스토크스 및 레이놀즈 수가 중간 수준인 조건에서 반응형 IB 수식의 병렬 고성능 계산 환경에 적합한 새로운 스무딩 전략을 개발한다.
• 다양한 재료 강성과 레이놀즈 수를 가진 유체-구조 상호작용 문제의 솔버 효율성과 확장성을 향상시킨다.
• 반응형 시간 이산화에 기반한 IB 방정식에서 발생하는 선형 시스템을 안정적이고 확장 가능한 방식으로 해결할 수 있도록 한다.
• 박막 및 두꺼운 구조 모델을 포함한 다양한 물리적 영역에서 신규 스무딩기의 성능을 조사한다.

제안 방법

• 스토크스-IB 시스템에서 수렴성은 유지하면서도 병렬 처리 성능을 향상시키기 위해 Vanka 스무딩기의 추가적 변형을 채택한다.
• 선형 시스템의 블록 구조를 활용하여 수렴 속도를 향상시키기 위해 슈어 여부 기반 스무딩기를 도입한다.
• 국소적으로 정렬된 카르테시안 격자에서 표준 갈레르킨 조정 및 보간 연산자를 사용한 기하학적 멀티그리드를 구현한다.
• 유럽계 및 라그랑주계 방정식의 공간 이산화에 대해 2차 유한 차분법을 사용한다.
• 정확하고 안정적인 라그랑주-유럽계 결합을 위해 페스킨의 4점 정규화 델타 함수를 적용한다.
• 카르레프 하위공간 솔버를 위한 기하학적 멀티그리드 조건부 해법에 스무딩기를 통합하여, 사다리꼴 시스템의 확장 가능한 해법을 가능하게 한다.

실험 결과

연구 질문

• RQ1 추가적 Vanka 스무딩기는 기하학적 멀티그리드 프레임워크에서 잠긴 경계 방정식에 효과적으로 적용될 수 있으며, 효율적인 병렬 처리를 가능하게 할 수 있는가?
• RQ2 전통적인 스무딩기와 비교해 슈어 여부 기반 스무딩기는 스토크스-IB 시스템에서 더 뛰어난 수렴 성능을 보일 수 있는가?
• RQ3 잠긴 경계 수식에서 레이놀즈 수와 재료 강성이 증가함에 따라 솔버 성능은 어떻게 저하되는가?
• RQ4 다양한 구조 두께와 물리적 매개변수에서 제안된 스무딩기의 확장성과 안정성은 어떠한가?
• RQ5 신규 스무딩기는 저-레이놀즈 수 유동뿐만 아니라 중간-고레이놀즈 수 유동에서도 유체-구조 상호작용 시뮬레이션에서 효율성을 유지할 수 있는가?

주요 결과

• 추가적 Vanka 스무딩기는 기존의 점 기반 Vanka 스무딩기의 대안으로 효과적이며, 분산 메모리 아키텍처에서 더 나은 확장성을 제공한다.
• 슈어 여부 기반 스무딩기는 특히 고강성 재료와 복잡한 인터페이스 기하학을 가진 문제에서 뛰어난 수렴 성능을 보였다.
• 선형 솔버 성능은 레이놀즈 수와 재료 강성이 증가함에 따라 저하되며, 특히 박막 인터페이스의 경우 두드러지게 나타나지만, 신규 스무딩기는 저-중간 레이놀즈 수 범위에서는 여전히 효과적이다.
• 제안된 솔버는 스토크스 및 비-스토크스 유동을 포함한 광범위한 물리적 매개변수에서 확장 가능한 성능을 달성한다.
• 스무딩기는 두꺼운 및 얇은 구조 모델 모두에서 안정성과 효율성을 유지하며, 다양한 탄성 강성 조건에서도 수렴 속도가 안정된 상태를 유지한다.
• 성능 프로파일링 결과, RAS 및 RMS 스무딩기는 잘 최적화되어 있으며 고성능 컴퓨팅 플랫폼에서 대규모 시뮬레이션에 적합하다.