[논문 리뷰] Sensitivity analysis on chaotic dynamical systems by Finite Difference Non-Intrusive Least Squares Shadowing (FD-NILSS) and application on a chaotic flow past a 3-D cylinder

Angxiu Ni, Qiqi Wang외 2명|arXiv (Cornell University)|2017. 11. 17.

Tropical and Extratropical Cyclones Research인용 수 2

한 줄 요약

이 논문은 카오스적 동역계에서 장시간 평균량의 민감도를 계산하기 위해 탄성 해석기(탄성 해석기 필요 없음)를 사용하지 않고 유한차분 방법을 비침습적으로 적용하는 FD-NILSS를 소개한다. 표준 시뮬레이션과 유사한 계산 비용으로 정확한 민감도 결과를 도출하였으며, Re = 525인 3차원 실린더 유동에서 검증하였다.

ABSTRACT

We present the Finite Difference Non-Intrusive Least Squares Shadowing (FD-NILSS) algorithm for computing sensitivities of long-time averaged quantities in chaotic dynamical systems. FD-NILSS does not require tangent solvers, and can be implemented with little modification to existing numerical simulation software. We also give a formula for solving the least-squares problem in FD-NILSS, which can be applied in NILSS as well. Finally, we apply FD-NILSS for sensitivity analysis of a chaotic flow over a 3-D cylinder at Reynolds number 525, where FD-NILSS computes accurate sensitivities and the computational cost is in the same order as the numerical simulation.

연구 동기 및 목표

• 카오스적 동역계에서 탄성 해석기를 필요로 하지 않는 비침습적 민감도 분석 방법을 개발하는 것.
• 기존 시뮬레이션 소프트웨어에 최소한의 수정으로만 민감도 계산을 가능하게 하는 것.
• 카오스적 유동에서 장시간 평균량의 민감도를 계산하기 위한 실용적이고 확장 가능한 접근법을 제공하는 것.
• 고레이놀즈 수를 가진 도전적인 3차원 난류 흐름 구성에서 방법을 검증하는 것.

제안 방법

• FD-NILSS는 민감도 계산에 필요한 탄성 공간을 유한차분으로 근사하여 분석적 또는 알고리즘적 탄성 해석기를 필요로 하지 않는다.
• 민감도 문제를 최소제곱 최적화 문제로 공식화하여 기준 궤도에서의 이탈을 최소화하는 샤로우 궤도를 구한다.
• 최소제곱 문제를 해결하기 위한 새로운 공식을 유도하였으며, 이는 FD-NILSS뿐 아니라 원래 NILSS 프레임워크에도 적용 가능하다.
• 알고리즘은 비침습적으로 구현되어 있으며, 전방 시뮬레이션 출력만 필요하고 기존 통합기의 수정이 필요하지 않다.
• 민감도 기울기는 시스템의 반응에 대한 유한차분 근사에서 유도된 작은 구조적 선형 시스템을 풀어 계산한다.

실험 결과

연구 질문

• RQ1 카오스적 시스템의 민감도 분석에서 탄성 해석기를 대체하기 위해 유한차분 근사가 정확도를 저하시키지 않고 사용될 수 있는가?
• RQ2 FD-NILSS의 계산 비용은 표준 시뮬레이션 및 기존 방법(예: NILSS)과 비교해 어떻게 되는가?
• RQ3 FD-NILSS는 카오스적 3차원 흐름에서 장시간 평균량의 민감도를 얼마나 정확하게 계산할 수 있는가?
• RQ4 FD-NILSS는 기존 시뮬레이션 소프트웨어에 최소한의 변경으로 구현될 수 있는가?

주요 결과

• FD-NILSS는 탄성 해석기를 필요로 하지 않고도 카오스적 시스템에서 장시간 평균량의 정확한 민감도를 계산할 수 있다.
• FD-NILSS의 계산 비용은 기준 수치 시뮬레이션과 비슷한 주기의 크기를 가지며, 계산적으로 효율적이다.
• 메서드는 레이놀즈 수 525인 3차원 난류 유동에서 민감도 추정에 높은 정확도를 달성한다.
• 유도된 최소제곱 해법 공식은 FD-NILSS와 원래 NILSS 메서드의 강건성과 적용 가능성을 향상시킨다.
• FD-NILSS의 비침습적 성격 덕분에 기존 시뮬레이션 코드에 최소한의 구현 노력으로 통합이 가능하다.