[논문 리뷰] A novel iterative penalty method to enforce boundary conditions in Finite Volume POD-Galerkin reduced order models for fluid dynamics problems

Sabrina Kelbij Star, Giovanni Stabile외 3명|arXiv (Cornell University)|2019. 12. 02.

Model Reduction and Neural Networks참고 문헌 60인용 수 11

한 줄 요약

이 논문은 유체역학을 위한 유한체적 기반 Proper Orthogonal Decomposition (POD)-Galerkin 감소 차원 모델(Rom)에서 시간에 따라 변화하는 경계 조건을 강제하기 위해 반복적 페널티 방법을 제안한다. 오차 허용 기반 반복 해법을 통해 페널티 인자를 자동으로 결정함으로써, 이 방법은 페널티 인자 수동 조정이 필요한 전통적 방법보다 더 빠르면서도 동일한 정확도를 달성한다. 리프팅 함수 방법과 동등한 정확도를 확보하면서도 리벳드 캐비티 유량의 경우 270–308×, Y-접합류 흐름의 경우 13–24×의 속도 향상을 제공한다.

ABSTRACT

A Finite-Volume based POD-Galerkin reduced order model is developed for fluid dynamics problems where the (time-dependent) boundary conditions are controlled using two different boundary control strategies: the lifting function method, whose aim is to obtain homogeneous basis functions for the reduced basis space and the penalty method where the boundary conditions are enforced in the reduced order model using a penalty factor. The penalty method is improved by using an iterative solver for the determination of the penalty factor rather than tuning the factor with a sensitivity analysis or numerical experimentation. The boundary control methods are compared and tested for two cases: the classical lid driven cavity benchmark problem and a Y-junction flow case with two inlet channels and one outlet channel. The results show that the boundaries of the reduced order model can be controlled with the boundary control methods and the same order of accuracy is achieved for the velocity and pressure fields. Finally, the reduced order models are 270-308 times faster than the full order models for the lid driven cavity test case and 13-24 times for the Y-junction test case.

연구 동기 및 목표

• 고정밀 유한체적(FV) 이산화로 유도된 감소 차원 모델(Rom)에서 시간에 따라 변화하는 경계 조건을 강제하는 데 도전하는 것.
• 감도 분석이나 수치 실험을 통해 수작업으로 페널티 인자를 조정하는 데 의존하는 전통적 페널티 방법의 한계를 극복하는 것.
• 정확도, 안정성 및 계산 효율성 측면에서 반복적 페널티 방법과 리프팅 함수 방법의 성능을 비교하는 것.
• 두 가지 벤치마크 유체역학 문제인 리벳드 캐비티와 다중 입구를 가진 Y-접합류 흐름에 대해 제안된 방법을 평가하는 것.
• 핵심 열역학 시스템 결합과 같은 산업 응용 분야에서 실시간 또는 근접 실시간 시뮬레이션 능력을 가능하게 하기 위해 Rom의 효율성과 경계 조건 제어 정밀도를 향상시키는 것.

제안 방법

• 압축성 없는 나비에-스토크스 방정식을 위한 유한체적 기반 POD-Galerkin Rom을 수립하며, 전체 순서 모델(FOM)을 적절한 정규직교 분해(POD)를 통해 구성된 감소 기저 공간에 투영한다.
• 페널티 항을 잔차에 추가하여 경계 조건을 약하게 강제하는 페널티 방법을 구현하며, 사용자가 정의한 오차 허용 범위를 목표로 하는 해법을 통해 페널티 인자를 반복적으로 결정한다.
• 수동 조정을 대체하기 위해 반복 알고리즘을 사용하여 페널티 인자를 자동으로 계산함으로써 감도 분석이나 시행착오 실험의 필요성을 줄인다.
• 두 가지 시험 케이스에 적용: 정체 및 비정상 상태인 리벳드 캐비티와 입구가 두 개인 유량이 하나인 Y-접합류 흐름으로, 모두 시간에 따라 변화하는 입구 경계 조건을 가짐.
• Rom의 안정성을 향상시키고 압력장을 재구성하기 위해 압력 푸아송 방정식을 통합함으로써 감소 차원 설정에서도 정확도와 안정성을 향상시킨다.
• 시간 이산화에 대해 Rom에서 2차 후진 차분법을 사용하여 FOM의 시간 정확도를 유지함으로써 일시적 시뮬레이션에서 일致성을 확보한다.

실험 결과

연구 질문

• RQ1 수동 조정이 아닌 반복적 페널티 방법을 사용하여 FV-POD-Galerkin Rom에서 시간에 따라 변화하는 경계 조건을 효과적으로 강제할 수 있는가?
• RQ2 반복적 페널티 방법이 유체역학 Rom의 속도 경계 조건을 강제할 때 리프팅 함수 방법과 비교해 정확도와 안정성 측면에서 어떻게 성능을 발휘하는가?
• RQ3 반복적 페널티 Rom이 기준 벤치마크 유동 문제에서 전체 순서 모델과 리프팅 함수 Rom에 비해 어떤 계산 효율성 향상을 달성하는가(속도 향상 비율)?
• RQ4 압력 푸아송 방정식은 Rom에서 재구성된 압력장의 정확도와 안정성에 얼마나 기여하는가?
• RQ5 다중 입구와 변동 가능한 입구 각을 가진 Y-접합류와 같은 복잡한 입구 구조에 적용했을 때, 이 방법의 견고성은 어느 정도인가?

주요 결과

• 반복적 페널티 방법은 리프팅 함수 방법과 동일한 정확도를 보이며, 속도장과 압력장 모두에서 신뢰할 수 있는 Rom 해를 생성한다.
• 반복 해법을 통해 각 시험 케이스당 약 1초 내에 페널티 인자를 자동으로 결정할 수 있으며, 수동 조정이나 감도 분석이 필요 없어진다.
• 리벳드 캐비티 시험 케이스에서 반복적 페널티 Rom은 전체 순서 모델 대비 308×의 속도 향상을 달성했고, 리프팅 함수 방법은 270×의 속도 향상을 기록했다.
• Y-접합류 시험 케이스에서는 반복적 페널티 Rom이 24×의 속도 향상을 기록했고, 리프팅 함수 방법은 13×의 속도 향상을 기록하여 복잡한 기하학적 구조에서 신규 방법의 성능 우월성을 보였다.
• 압력 푸아송 방정식 접근법은 Rom의 안정성을 효과적으로 향상시키고 압력장 재구성 정확도를 향상시켰지만, 압력의 상대 오차는 여전히 속도의 약 한 단계 높은 편이었다.
• 스냅샷 수집 시 사용된 입구 조건 범위 내에서 Rom을 테스트했을 때 방법은 안정적이고 정확했지만, 이 범위를 초월하면 성능 저하가 발생할 수 있었으며, 이는 대표성 있는 스퍼블릿 수집의 중요성을 강조한다.